Genetic Constraints Research Articles

Quantifying the Mutational Robustness of Protein-Coding Genes.

We use large-scale mutagenesis data and computer simulations to quantify the mutational robustness of protein-coding genes by taking into account constraints arising from protein function and the genetic code. Analyses of the distribution of amino acid substitutions from 18 mutagenesis studies revealed an average of 45% of neutral variants; while mutagenesis data of 12 proteins artificially designed under no other constraints but stability, reach an average of 60%. Simulations using a lattice protein model allow us to contrast these estimates to the expected mutational robustness of protein families by generating unbiased samples of foldable sequences, which we find to have 30% of neutral variants. In agreement with mutagenesis data of designed proteins, the model shows that maximally robust protein families might access up to twice the amount of neutral variants observed in the unbiased samples (i.e. 60%). A biophysical model of protein-ligand binding suggests that constraints associated to molecular function have only a moderate impact on robustness of approximately 5 to 10% of neutral variants; and that the direction of this effect depends on the relation between functional performance and thermodynamic stability. Although the genetic code constraints the access of a gene's nucleotide sequence to only 30% of the full distribution of amino acid mutations, it provides an extra 15 to 20% of neutral variants to the estimations above, such that the expected, observed, and maximal robustness of protein-coding genes are approximately 50, 65, and 75%, respectively. We discuss our results in the light of three main hypothesis put forward to explain the existence of mutationally robust genes.

Homage to Felsenstein 1981, or why are there so few/many species?

If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build-up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between "two-allele" models, which are subject to this effect, and "one-allele" models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.

Quantitative trait locus mapping reveals an independent genetic basis for joint divergence in leaf function, life-history, and floral traits between scarlet monkeyflower (Mimulus cardinalis) populations.

Across taxa, vegetative and floral traits that vary along a fast-slow life-history axis are often correlated with leaf functional traits arrayed along the leaf economics spectrum, suggesting a constrained set of adaptive trait combinations. Such broad-scale convergence may arise from genetic constraints imposed by pleiotropy (or tight linkage) within species, or from natural selection alone. Understanding the genetic basis of trait syndromes and their components is key to distinguishing these alternatives and predicting evolution in novel environments. We used a line-cross approach and quantitative trait locus (QTL) mapping to characterize the genetic basis of twenty leaf functional/physiological, life history, and floral traits in hybrids between annualized and perennial populations of scarlet monkeyflower (Mimulus cardinalis). We mapped both single and multi-trait QTLs for life history, leaf function and reproductive traits, but found no evidence of genetic co-ordination across categories. A major QTL for three leaf functional traits (thickness, photosynthetic rate, and stomatal resistance) suggests that a simple shift in leaf anatomy may be key to adaptation to seasonally dry habitats. Our results suggest that the co-ordination of resource-acquisitive leaf physiological traits with a fast life-history and more selfing mating system results from environmental selection rather than functional or genetic constraint. Independent assortment of distinct trait modules, as well as a simple genetic basis to leaf physiological traits associated with drought escape, may facilitate adaptation to changing climates.

Evolutionary History of Plant Metabolism.

Tremendous chemical diversity is the hallmark of plants and is supported by highly complex biochemical machinery. Plant metabolic enzymes originated and were transferred from eukaryotic and prokaryotic ancestors and further diversified by the unprecedented rates of gene duplication and functionalization experienced in land plants. Unlike microbes, which have frequent horizontal gene transfer events and multiple inputs of energy and organic carbon, land plants predominantly rely on organic carbon generated from CO2 and have experienced very few, if any, gene transfers during their recent evolutionary history. As such, plant metabolic networks have evolved in a stepwise manner and on existing networks under various evolutionary constraints. This review aims to take a broader view of plant metabolic evolution and lay a framework to further explore evolutionary mechanisms of the complex metabolic network. Understanding the underlying metabolic and genetic constraints is also an empirical prerequisite for rational engineering and redesigning of plant metabolic pathways.

Introduced populations of ragweed show as much evolutionary potential as native populations.

Invasive species are a global economic and ecological problem. They also offer an opportunity to understand evolutionary processes in a colonizing context. The impacts of evolutionary factors, such as genetic variation, on the invasion process are increasingly appreciated, but there remain gaps in the empirical literature. The adaptive potential of populations can be quantified using genetic variance–covariance matrices (G), which encapsulate the heritable genetic variance in a population. Here, we use a multivariate Bayesian approach to assess the adaptive potential of invasive populations of ragweed (Ambrosia artemisiifolia), a serious allergen and agricultural weed. We compared several aspects of genetic architecture and the structure of G matrices between three native and three introduced populations, based on phenotypic data collected in a field common garden experiment. We found moderate differences in the quantitative genetic architecture among populations, but we did not find that introduced populations suffer from a limited adaptive potential or increased genetic constraint compared with native populations. Ragweed has an annual life history, is an obligate outcrosser, and produces very large numbers of seeds and pollen grains. These characteristics, combined with the significant additive genetic variance documented here, suggest ragweed will be able to respond quickly to selection pressures in both its native and introduced ranges.

The genetic architecture of sexual dimorphism in the moss Ceratodon purpureus.

A central problem in evolutionary biology is to identify the forces that maintain genetic variation for fitness in natural populations. Sexual antagonism, in which selection favours different variants in males and females, can slow the transit of a polymorphism through a population or can actively maintain fitness variation. The amount of sexually antagonistic variation to be expected depends in part on the genetic architecture of sexual dimorphism, about which we know relatively little. Here, we used a multivariate quantitative genetic approach to examine the genetic architecture of sexual dimorphism in a scent-based fertilization syndrome of the moss Ceratodon purpureus. We found sexual dimorphism in numerous traits, consistent with a history of sexually antagonistic selection. The cross-sex genetic correlations (rmf) were generally heterogeneous with many values indistinguishable from zero, which typically suggests that genetic constraints do not limit the response to sexually antagonistic selection. However, we detected no differentiation between the female- and male-specific trait (co)variance matrices (Gf and Gm, respectively), meaning the evolution of sexual dimorphism may be constrained. The cross-sex cross-trait covariance matrix B contained both symmetric and asymmetric elements, indicating that the response to sexually antagonistic or sexually concordant selection, and the constraint to sexual dimorphism, are highly dependent on the traits experiencing selection. The patterns of genetic variances and covariances among these fitness components is consistent with partly sex-specific genetic architectures having evolved in order to partially resolve multivariate genetic constraints (i.e. sexual conflict), enabling the sexes to evolve towards their sex-specific multivariate trait optima.

Distinct genetic basis for root responses to lipo-chitooligosaccharide signal molecules from different microbial origins.

Lipo-chitooligosaccharides (LCOs) were originally found as symbiotic signals called Nod Factors (Nod-LCOs) controlling the nodulation of legumes by rhizobia. More recently, LCOs were also found in symbiotic fungi and, more surprisingly, very widely in the kingdom Fungi, including in saprophytic and pathogenic fungi. The LCO-V(C18:1, fucosylated/methyl fucosylated), hereafter called Fung-LCOs, are the LCO structures most commonly found in fungi. This raises the question of how legume plants such as Medicago truncatula can discriminate between Nod-LCOs and Fung-LCOs. To address this question, we performed a genome-wide association study on 173 natural accessions of M. truncatula, using a root branching phenotype and a newly developed local score approach. Both Nod-LCOs and Fung-LCOs stimulated root branching in most accessions, but the root responses to these two types of LCO molecules were not correlated. In addition, the heritability of the root response was higher for Nod-LCOs than for Fung-LCOs. We identified 123 loci for Nod-LCO and 71 for Fung-LCO responses, of which only one was common. This suggests that Nod-LCOs and Fung-LCOs both control root branching but use different molecular mechanisms. The tighter genetic constraint of the root response to Fung-LCOs possibly reflects the ancestral origin of the biological activity of these molecules.

Distinct genetic pathways define pre-malignant versus compensatory clonal hematopoiesis in Shwachman-Diamond syndrome

To understand the mechanisms that mediate germline genetic leukemia predisposition, we studied the inherited ribosomopathy Shwachman-Diamond syndrome (SDS), a bone marrow failure disorder with high risk of myeloid malignancies at an early age. To define the mechanistic basis of clonal hematopoiesis in SDS, we investigate somatic mutations acquired by patients with SDS followed longitudinally. Here we report that multiple independent somatic hematopoietic clones arise early in life, most commonly harboring heterozygous mutations in EIF6 or TP53. We show that germline SBDS deficiency establishes a fitness constraint that drives selection of somatic clones via two distinct mechanisms with different clinical consequences. EIF6 inactivation mediates a compensatory pathway with limited leukemic potential by ameliorating the underlying SDS ribosome defect and enhancing clone fitness. TP53 mutations define a maladaptive pathway with enhanced leukemic potential by inactivating tumor suppressor checkpoints without correcting the ribosome defect. Subsequent development of leukemia was associated with acquisition of biallelic TP53 alterations. These results mechanistically link leukemia predisposition to germline genetic constraints on cellular fitness, and provide a rational framework for clinical surveillance strategies.

Titin Circular RNAs Create a Back-Splice Motif Essential for SRSF10 Splicing.

TTN (Titin), the largest protein in humans, forms the molecular spring that spans half of the sarcomere to provide passive elasticity to the cardiomyocyte. Mutations that disrupt the TTN transcript are the most frequent cause of hereditary heart failure. We showed before that TTN produces a class of circular RNAs (circRNAs) that depend on RBM20 to be formed. In this study, we show that the back-splice junction formed by this class of circRNAs creates a unique motif that binds SRSF10 to enable it to regulate splicing. Furthermore, we show that one of these circRNAs (cTTN1) distorts both localization of and splicing by RBM20. We calculated genetic constraint of the identified motif in 125 748 exomes collected from the gnomAD database. Furthermore, we focused on the highest expressed RBM20-dependent circRNA in the human heart, which we named cTTN1. We used shRNAs directed to the back-splice junction to induce selective loss of cTTN1 in human induced pluripotent stem cell-derived cardiomyocytes. Human genetics suggests reduced genetic tolerance of the generated motif, indicating that mutations in this motif might lead to disease. RNA immunoprecipitation confirmed binding of circRNAs with this motif to SRSF10. Selective loss of cTTN1 in human induced pluripotent stem cell-derived cardiomyocytes induced structural abnormalities, apoptosis, and reduced contractile force in engineered heart tissue. In line with its SRSF10 binding, loss of cTTN1 caused abnormal splicing of important cardiomyocyte SRSF10 targets such as MEF2A and CASQ2. Strikingly, loss of cTTN1 also caused abnormal splicing of TTN itself. Mechanistically, we show that loss of cTTN1 distorts both localization of and splicing by RBM20. We demonstrate that circRNAs formed from the TTN transcript are essential for normal splicing of key muscle genes by enabling splice regulators RBM20 and SRSF10. This shows that the TTN transcript also has regulatory roles, besides its well-known signaling and structural function. In addition, we demonstrate that the specific sequence created by the back-splice junction of these circRNAs has important functions. This highlights the existence of functionally important sequences that cannot be recognized as such in the human genome but provides an as-yet unrecognized source for functional sequence variation.

Mapping the adaptive landscape of a major agricultural pathogen reveals evolutionary constraints across heterogeneous environments

The adaptive potential of pathogens in novel or heterogeneous environments underpins the risk of disease epidemics. Antagonistic pleiotropy or differential resource allocation among life-history traits can constrain pathogen adaptation. However, we lack understanding of how the genetic architecture of individual traits can generate trade-offs. Here, we report a large-scale study based on 145 global strains of the fungal wheat pathogen Zymoseptoria tritici from four continents. We measured 50 life-history traits, including virulence and reproduction on 12 different wheat hosts and growth responses to several abiotic stressors. To elucidate the genetic basis of adaptation, we used genome-wide association mapping coupled with genetic correlation analyses. We show that most traits are governed by polygenic architectures and are highly heritable suggesting that adaptation proceeds mainly through allele frequency shifts at many loci. We identified negative genetic correlations among traits related to host colonization and survival in stressful environments. Such genetic constraints indicate that pleiotropic effects could limit the pathogen’s ability to cause host damage. In contrast, adaptation to abiotic stress factors was likely facilitated by synergistic pleiotropy. Our study illustrates how comprehensive mapping of life-history trait architectures across diverse environments allows to predict evolutionary trajectories of pathogens confronted with environmental perturbations.

Virulence during Newcastle Disease Viruses Cross Species Adaptation.

The hypothesis that host adaptation in virulent Newcastle disease viruses (NDV) has been accompanied by virulence modulation is reviewed here. Historical records, experimental data, and phylogenetic analyses from available GenBank sequences suggest that currently circulating NDVs emerged in the 1920–1940′s from low virulence viruses by mutation at the fusion protein cleavage site. These viruses later gave rise to multiple virulent genotypes by modulating virulence in opposite directions. Phylogenetic and pathotyping studies demonstrate that older virulent NDVs further evolved into chicken-adapted genotypes by increasing virulence (velogenic-viscerotropic pathotypes with intracerebral pathogenicity indexes [ICPIs] of 1.6 to 2), or into cormorant-adapted NDVs by moderating virulence (velogenic–neurotropic pathotypes with ICPIs of 1.4 to 1.6), or into pigeon-adapted viruses by further attenuating virulence (mesogenic pathotypes with ICPIs of 0.9 to 1.4). Pathogenesis and transmission experiments on adult chickens demonstrate that chicken-adapted velogenic-viscerotropic viruses are more capable of causing disease than older velogenic-neurotropic viruses. Currently circulating velogenic–viscerotropic viruses are also more capable of replicating and of being transmitted in naïve chickens than viruses from cormorants and pigeons. These evolutionary virulence changes are consistent with theories that predict that virulence may evolve in many directions in order to achieve maximum fitness, as determined by genetic and ecologic constraints.

Genetic variation underlies the plastic response to shade of snapdragon plants (Antirrhinum majus L.)

ABSTRACT A classic example of phenotypic plasticity in plants is the set of traits that change in response to shade. There is widespread evidence that plants in low light conditions often avoid shade by growing taller or by increasing their photosynthetic efficiency, i.e. the shade avoidance syndrome. Whether this plasticity might evolve in response to natural selection depends upon the presence of its standing genetic variation in wild populations. There is limited evidence for heritable standing variation in the plastic response of plants to shade. In this study, we used an experimental common garden approach to investigate this plastic response in snapdragon plants (Antirrhinum majus L.) originating from four natural populations from the Mediterranean region. Our results showed that individual plants reacted strongly to the presence of shade by growing longer shoots, longer internodes, and increasing their specific leaf area in these four populations. Our results also revealed genetic variation for the plastic response within these populations, as well as little genetic constraints to its evolution. Our findings imply that natural populations of A. majus harbour standing genetic variation for phenotypic plasticity in response to shade, providing them the potential to evolve in response to selection.

How does the strength of selection influence genetic correlations?

Genetic correlations between traits can strongly impact evolutionary responses to selection, and may thus impose constraints on adaptation. Theoretical and empirical work has made it clear that without strong linkage and with random mating, genetic correlations at evolutionary equilibrium result from an interplay of correlated pleiotropic effects of mutations, and correlational selection favoring combinations of trait values. However, it is not entirely clear how change in the overall strength of stabilizing selection across traits (breadth of the fitness peak, given its shape) influences this compromise between mutation and selection effects on genetic correlation. Here, we show that the answer to this question crucially depends on the intensity of genetic drift. In large, effectively infinite populations, genetic correlations are unaffected by the strength of selection, regardless of whether the genetic architecture involves common small‐effect mutations (Gaussian regime), or rare large‐effect mutations (House‐of‐Cards regime). In contrast in finite populations, the strength of selection does affect genetic correlations, by shifting the balance from drift‐dominated to selection‐dominated evolutionary dynamics. The transition between these domains depends on mutation parameters to some extent, but with a similar dependence of genetic correlation on the strength of selection. Our results are particularly relevant for understanding how senescence shapes patterns of genetic correlations across ages, and genetic constraints on adaptation during colonization of novel habitats.

Birdsong Learning and Culture: Analogies with Human Spoken Language

Unlike many species, song learning birds and humans have independently evolved the ability to communicate via learned vocalizations. Both birdsong and spoken language are culturally transmitted across generations, within species-specific constraints that leave room for considerable variation. We review the commonalities and differences between vocal learning bird species and humans, across behavioral, developmental, neuroanatomical, physiological, and genetic levels. We propose that cultural transmission of vocal repertoires is a natural consequence of the evolution of vocal learning and that at least some species-specific universals, as well as species differences in cultural transmission, are due to differences in vocal learning phenotypes, which are shaped by genetic constraints. We suggest that it is the balance between these constraints and features of the social environment that allows cultural learning to propagate. We describe new opportunities for exploring meaningful comparisons of birdsong and human vocal culture.

Avian phenotypic convergence is subject to low genetic constraints based on genomic evidence

BackgroundPhenotypic convergence between distinct species provides an opportunity to examine the predictability of genetic evolution. Unrelated species sharing genetic underpinnings for phenotypic convergence suggests strong genetic constraints, and thus high predictability of evolution. However, there is no clear big picture of the genomic constraints on convergent evolution. Genome-based phylogenies have confirmed many cases of phenotypic convergence in birds, making them a good system for examining genetic constraints in phenotypic convergence. In this study, we used hierarchical genomic approaches to estimate genetic constraints in three convergent avian traits: nocturnality, raptorial behavior and foot-propelled diving.ResultsPhylogeny-based hypothesis tests and positive selection tests were applied to compare 16 avian genomes, representing 14 orders, and identify genes with strong convergence signals. We found 43 adaptively convergent genes (ACGs) associated with the three phenotypic convergence cases and assessed genetic constraints in all three cases, from (amino acid) site mutations to genetic pathways. We found that the avian orders shared few site mutations in the ACGs that contributed to the convergent phenotypes, and that these ACGs were not enriched in any genetic pathways. In addition, different pairs of orders with convergent foot-propelled diving or raptorial behaviors shared few ACGs. We also found that closely related orders that shared foot-propelled diving behavior did not share more ACGs than did distinct orders, suggesting that convergence among these orders could not be explained by their initial genomic backgrounds.ConclusionsOur analyses of three avian convergence events suggest low constraints for phenotypic convergence across multiple genetic levels, implying that genetic evolution is unpredictable at the phylogenetic level of avian order. Ours is one of first studies to apply hierarchical genomic examination to multiple avian convergent cases to assess the genetic constraints in life history trait evolution.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Multi-biomarker approach and IBR index to evaluate the effects of different contaminants on the ecotoxicological status of Apis mellifera

The honeybee, Apis mellifera L. (Hymenoptera: Apidae), a keystone pollinator of wild plant species and agricultural crops, is disappearing globally due to parasites and diseases, habitat loss, genetic constraints, beekeeper management issues and to the widespread use of pesticides. Besides insecticides, widely studied in this species, honeybees are also exposed to herbicides and fungicides and heavy metals whose lethal and sublethal effects need to be investigated. In this context, our study aimed to evaluate the effects of fungicides and of heavy metals on honeybees and to develop and apply a multi-biomarker approach that include an Integrated Biological Index (IBRv2) to assess the toxicological status of this species. Biomarkers of neurotoxicity (AChE and CaE), metabolic alteration (ALP, and GST) and immune system (LYS, granulocytes) were measured, following honeybees’ exposure to cadmium or to a crop fungicide, using the genotoxic compound EMS as positive control. A biomarker of genotoxicity (NA assay) was developed and applied for the first time in honeybees. At the doses tested, all the contaminants showed sublethal toxicity to the bees, highlighting in particular genotoxic effects. The data collected were analyzed by an IBRv2 index, which integrated the seven biomarkers used in this study. IBRv2 index increased with increasing cadmium or fungicide concentrations. The IBRv2 represents a simple tool for a general description of honeybees ecotoxicological health status. Results highlight the need for more in-depth investigations on the effects of fungicides on non-target organisms, such as honeybees, using sensitive methods for the determination of sublethal effects. This study contributes to the development of a multi-biomarker approach to be used for a more accurate ecotoxicological environmental monitoring of these animals.

A Bayesian assessment of the conservation status of humpback whales (Megaptera novaeangliae) in the western South Atlantic Ocean

The population of humpback whales (Megaptera novaeangliae) wintering off the eastern coast of South America is referred to by the International Whaling Commission as ‘Breeding Stock A’ (BSA). This population was heavily exploited in 20th century modern commercial whaling operations. After more than 30 years of protection, its present status remains unknown. A deterministic sex and age-aggregated population dynamics model was used to estimate the pre-exploitation population size (K), the maximum net recruitment rate (rmax), the maximum depletion level (Nmin/K), and other quantities of interest of BSA. Input data included modern whaling catch series, absolute estimates of abundance, observed growth rates and indices of relative abundance. A Bayesian statistical method was used to calculate probability distributions for the model parameters. Prior distributions were set on rmax – an uninformative (Uniform [0, 0.106]) and an informative (Normal [0.067, 0.042]) – and on the population size in 2005 – N2005 (Uniform [500, 22,000]). A total of 10,000 samples were used to compute the joint posterior distribution of the model parameters using the Sampling-Importance-Resampling algorithm. Sensitivity of model outputs to the priors on rmax, a genetic constraint, data inclusion and catch allocation scenarios was investigated. Medians of the posterior probability distributions of quantities of interest for the base case scenario were: rmax = 0.069 (95% probability intervals [PI] = 0.013–0.104), K = 24,558 (95% PI = 22,791–31,118), Nmin/K = 2% (PI = 0.31%–12.5%), N2006/ K = 27.4% (PI = 18.3%–39.5%), N2020/K = 61.8% (PI = 23.8%–88.6%), and N2040/K = 97.3% (PI = 31.6%–99.9%). Despite apparent recovery in the past three decades, the western South Atlantic humpback whale population is still low relative to its pre-exploitation size and requires continued conservation efforts.

Genomic differentiation in an endemic Philippine genus (Aves: Sarcophanops) owing to geographical isolation on recently disassociated islands.

Phylogeographical studies of Philippine vertebrates have demonstrated that genetic variation is broadly partitioned by Pleistocene island aggregation. Contemporary island discontinuity is expected to influence genetic differentiation but remains relatively undocumented, perhaps because the current episode of island isolation started in relatively recent times. We investigated inter- and intra-island population structure in a Philippine endemic bird genus (Sarcophanops) to determine whether genetic differentiation has evolved during the recent period of isolation. We sequenced thousands of genome-wide restriction site associated DNA (RAD) markers from throughout the Mindanao group to assess fine-scale genetic structure across islands. Specifically, we investigated patterns of gene flow and connectivity within and between taxonomic and geographical bounds. A previous assessment of mitochondrial DNA detected deep structure between Sarcophanops samarensis and a sister species, Sarcophanops steerii, but was insufficient to detect differentiation within either species. Analysis of RAD markers, however, revealed structure within S. samarensis between the islands of Samar/Leyte and Bohol. This genetic differentiation probably demonstrates an effect of recent geographical isolation (after the Last Glacial Maximum) on the genetic structure of Philippine avifauna. We suggest that the general lack of evidence for differentiation between recently isolated populations is a failure to detect subtle population structure owing to past genetic sampling constraints, rather than the absence of such structure.

Temporal, spatial, and genetic constraints contribute to the patterning and penetrance of murine neurofibromatosis-1 optic glioma.

Brain tumors are the most common solid tumors of childhood, but little is understood about the factors that influence their development. Pediatric low-grade gliomas in particular display unique temporal and spatial localization associated with different genetic mutations (eg, BRAF genomic alterations, mutations in the neurofibromatosis type 1 [NF1] gene) for reasons that remain unclear. NF1 low-grade gliomas typically arise in the optic pathway of young children as optic pathway gliomas (OPGs), likely from a cell of origin that resides within the third ventricular zone (TVZ). However, the factors that contribute to their distinct temporal patterning and penetrance have not been adequately explored. TVZ neuroglial progenitor cells (NPCs) were analyzed over the course of mouse brain development. Progenitors isolated by fluorescence-activated cell sorting (FACS) were assessed for functional and molecular differences. The impact of different germline Nf1 mutations on TVZ NPC properties was analyzed using genetically engineered mice. We identify 3 individual factors that could each contribute to Nf1 optic glioma temporal patterning and penetrance. First, there are 3 functionally and molecularly distinct populations of mouse TVZ NPCs, one of which ("M" cells) exhibits the highest clonogenic incidence, proliferation, and abundance during embryogenesis. Second, TVZ NPC proliferation dramatically decreases after birth. Third, germline Nf1 mutations differentially increase TVZ NPC proliferation during embryogenesis. The unique temporal patterning and penetrance of Nf1 optic glioma reflects the combined effects of TVZ NPC population composition, time-dependent changes in progenitor proliferation, and the differential impact of the germline Nf1 mutation on TVZ NPC expansion.