Experimental Crossing Confirms Reproductive Isolation between Cryptic Species within Eulimnogammarus verrucosus (Crustacea: Amphipoda) from Lake Baikal

Ancient lakes are known speciation hotspots. One of the most speciose groups in the ancient Lake Baikal are gammaroid amphipods (Crustacea: Amphipoda: Gammaroidea). There are over 350 morphological species and subspecies of amphipods in Baikal, but the extent of cryptic variation is still unclear. One of the most common species in the littoral zone of the lake, Eulimnogammarus verrucosus (Gerstfeldt, 1858), was recently found to comprise at least three (pseudo)cryptic species based on molecular data. Here, we further explored these species by analyzing their mitogenome-based phylogeny, genome sizes with flow cytometry, and their reproductive compatibility. We found divergent times of millions of years and different genome sizes in the three species (6.1, 6.9 and 8 pg), further confirming their genetic separation. Experimental crossing of the western and southern species, which are morphologically indistinguishable and have adjacent ranges, showed their separation with a post-zygotic reproductive barrier, as hybrid embryos stopped developing roughly at the onset of gastrulation. Thus, the previously applied barcoding approach effectively indicated the separate biological species within E. verrucosus. These results provide new data for investigating genome evolution and highlight the need for precise tracking of the sample origin in any studies in this morphospecies.

Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato

Structural variants (SVs) underlie important crop improvement and domestication traits. However, resolving the extent, diversity, and quantitative impact of SVs has been challenging. We used long-read nanopore sequencing to capture 238,490 SVs in 100 diverse tomato lines. This panSV genome, along with 14 new reference assemblies, revealed large-scale intermixing of diverse genotypes, as well as thousands of SVs intersecting genes and cis-regulatory regions. Hundreds of SV-gene pairs exhibit subtle and significant expression changes, which could broadly influence quantitative trait variation. By combining quantitative genetics with genome editing, we show how multiple SVs that changed gene dosage and expression levels modified fruit flavor, size, and production. In the last example, higher order epistasis among four SVs affecting three related transcription factors allowed introduction of an important harvesting trait in modern tomato. Our findings highlight the underexplored role of SVs in genotype-to-phenotype relationships and their widespread importance and utility in crop improvement.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

A Very Oil Yellow1 Modifier of the Oil Yellow1-N1989 Allele Uncovers a Cryptic Phenotypic Impact of Cis-regulatory Variation in Maize

Forward genetics determines the function of genes underlying trait variation by identifying the change in DNA responsible for changes in phenotype. Detecting phenotypically-relevant variation outside protein coding sequences and distinguishing this from neutral variants is not trivial; partly because the mechanisms by which DNA polymorphisms in the intergenic regions affect gene regulation are poorly understood. Here we utilized a dominant genetic reporter to investigate the effect of cis and trans-acting regulatory variation. We performed a forward genetic screen for natural variation that suppressed or enhanced the semi-dominant mutant allele Oy1-N1989, encoding the magnesium chelatase subunit I of maize. This mutant permits rapid phenotyping of leaf color as a reporter for chlorophyll accumulation, and mapping of natural variation in maize affecting chlorophyll metabolism. We identified a single modifier locus segregating between B73 and Mo17 that was linked to the reporter gene itself, which we call very oil yellow1 (vey1). Based on the variation in OY1 transcript abundance and genome-wide association data, vey1 is predicted to consist of multiple cis-acting regulatory sequence polymorphisms encoded at the wild-type oy1 alleles. The vey1 locus appears to be a common polymorphism in the maize germplasm that alters the expression level of a key gene in chlorophyll biosynthesis. These vey1 alleles have no discernable impact on leaf chlorophyll in the absence of the Oy1-N1989 reporter. Thus, the use of a mutant as a reporter for magnesium chelatase activity resulted in the detection of expression-level polymorphisms not readily visible in the laboratory.

The devil is in the details: genetic variation in introduced populations and its contributions to invasion

The influence of genetic variation on invasion success has captivated researchers since the start of the field of invasion genetics 50 years ago. We review the history of work on this question and conclude that genetic variation-as surveyed with molecular markers-appears to shape invasion rarely. Instead, there is a significant disconnect between marker assays and ecologically relevant genetic variation in introductions. We argue that the potential for adaptation to facilitate invasion will be shaped by the details of genotypes affecting phenotypes, and we highlight three areas in which we see opportunities to make powerful new insights. (i) The genetic architecture of adaptive variation. Traits shaped by large-effect alleles may be strongly impacted by founder events yet more likely to respond to selection when genetic drift is strong. Large-effect loci may be especially relevant for traits involved in biotic interactions. (ii) Cryptic genetic variation exposed during invasion. Introductions have strong potential to uncover masked variation due to alterations in genetic and ecological environments. (iii) Genetic interactions during admixture of multiple source populations. As divergence among sources increases, positive followed by increasingly negative effects of admixture should be expected. Although generally hypothesized to be beneficial during invasion, admixture is most often reported among sources of intermediate divergence, supporting the possibility that incompatibilities among divergent source populations might be limiting their introgression. Finally, we note that these details of invasion genetics can be coupled with comparative demographic analyses to link genetic changes to the evolution of invasiveness itself.

How much of the variation in the mutation rate along the human genome can be explained?

It has been claimed recently that it may be possible to predict the rate of de novo mutation of each site in the human genome with a high degree of accuracy [Michaelson et al. (2012), Cell 151: 1431−1442]. We show that this claim is unwarranted. By considering the correlation between the rate of de novo mutation and the predictions from the model of Michaelson et al., we show there could be substantial unexplained variance in the mutation rate. We investigate whether the model of Michaelson et al. captures variation at the single nucleotide level that is not due to simple context. We show that the model captures a substantial fraction of this variation at CpG dinucleotides but fails to explain much of the variation at non-CpG sites.

Laboratory evolution of adenylyl cyclase independent learning in Drosophila and missing heritability

Gene interactions are acknowledged to be a likely source of missing heritability in large-scale genetic studies of complex neurological phenotypes. However, involvement of rare variants, de novo mutations, genetic lesions that are not easily detected with commonly used methods and epigenetic factors also are possible explanations. We used a laboratory evolution study to investigate the modulatory effects of background genetic variation on the phenotypic effect size of a null mutation with known impact on olfactory learning. To accomplish this, we first established a population that contained variation at just 23 loci and used selection to evolve suppression of the learning defect seen with null mutations in the rutabaga adenylyl cyclase. We thus biased the system to favor relatively simplified outcomes by choosing a Mendelian trait and by restricting the genetic variation segregating in the population. This experimental design also assures that the causal effects are among the known 23 segregating loci. We observe a robust response to selection that requires the presence of the 23 variants. Analyses of the underlying genotypes showed that interactions between more than two loci are likely to be involved in explaining the selection response, with implications for the missing heritability problem.

Stressful environments induce novel phenotypic variation: hierarchical reaction norms for sperm performance of a pervasive invader

Genetic variation for phenotypic plasticity is ubiquitous and important. However, the scale of such variation including the relative variability present in reaction norms among different hierarchies of biological organization (e.g., individuals, populations, and closely related species) is unknown. Complicating interpretation is a trade-off in environmental scale. As plasticity can only be inferred over the range of environments tested, experiments focusing on fine tuned responses to normal or benign conditions may miss cryptic phenotypic variation expressed under novel or stressful environments. Here, we sought to discern the presence and shape of plasticity in the performance of brown trout sperm as a function of optimal to extremely stressful river pH, and demarcate if the reaction norm varies among genotypes. Our overarching goal was to determine if deteriorating environmental quality increases expressed variation among individuals. A more applied aim was to ascertain whether maintaining sperm performance over a wide pH range could help explain how brown trout are able to invade diverse river systems when transplanted outside of their native range. Individuals differed in their reaction norms of phenotypic expression of an important trait in response to environmental change. Cryptic variation was revealed under stressful conditions, evidenced through increasing among-individual variability. Importantly, data on population averages masked this variability in plasticity. In addition, canalized reaction norms in sperm swimming velocities of many individuals over a very large range in water chemistry may help explain why brown trout are able to colonize a wide variety of habitats.

Intraspecific evolution of the intercellular signaling network underlying a robust developmental system

Many biological systems produce an invariant output when faced with stochastic or environmental variation. This robustness of system output to variation affecting the underlying process may allow for "cryptic" genetic evolution within the system without change in output. We studied variation of cell fate patterning of Caenorhabditis elegans vulva precursors, a developmental system that relies on a simple intercellular signaling network and yields an invariant output of cell fates and lineages among C. elegans wild isolates. We first investigated the system's genetic variation in C. elegans by means of genetic tools and cell ablation to break down its buffering mechanisms. We uncovered distinct architectures of quantitative variation along the Ras signaling cascade, including compensatory variation, and differences in cell sensitivity to induction along the anteroposterior axis. In the unperturbed system, we further found variation between isolates in spatio-temporal dynamics of Ras pathway activity, which can explain the phenotypic differences revealed upon perturbation. Finally, the variation mostly affects the signaling pathways in a tissue-specific manner. We thus demonstrate and characterize microevolution of a developmental signaling network. In addition, our results suggest that the vulva genetic screens would have yielded a different mutation spectrum, especially for Wnt pathway mutations, had they been performed in another C. elegans genetic background.

HSP90 affects the expression of genetic variation and developmental stability in quantitative traits

Modulation of the activity of the molecular chaperone HSP90 has been extensively discussed as a means to alter phenotype in many traits and organisms. Such changes can be due to the exposure of cryptic genetic variation, which in some instances may also be accomplished by mild environmental alteration. Should such polymorphisms be widespread, natural selection may be more effective at producing phenotypic change in suboptimal environments. However, the frequency and identity of buffered polymorphisms in natural populations are unknown. Here, we employ quantitative genetic dissection of an Arabidopsis thaliana developmental response, hypocotyl elongation in the dark, to detail the underpinnings of genetic variation responsive to HSP90 modulation. We demonstrate that HSP90-dependent alleles occur in continuously distributed, environmentally responsive traits and are amenable to quantitative genetic mapping techniques. Furthermore, such alleles are frequent in natural populations and can have significant effects on natural phenotypic variation. We also find that HSP90 modulation has both general and allele-specific effects on developmental stability; that is, developmental stability is a phenotypic trait that can be affected by natural variation. However, effects of revealed variation on trait means outweigh effects of decreased developmental stability, and the HSP90-dependent trait alterations could be acted on by natural selection. Thus, HSP90 may centrally influence canalization, assimilation, and the rapid evolutionary alteration of phenotype through the concealment and exposure of cryptic genetic variation.

HSP90-buffered genetic variation is common in Arabidopsis thaliana

HSP90 is a protein chaperone particularly important in the maturation of a diverse set of proteins that regulate key steps in a multitude of biological processes. Alterations in HSP90 function produce altered phenotypes at low penetrance in natural populations. Previous work has shown that at least some of these phenotypes are due to genetic variation that remains phenotypically cryptic until it is revealed by the impairment of HSP90 function. Exposure of such "buffered" genetic polymorphisms can also be accomplished by environmental stress, linking the appearance of new phenotypes to defects in protein homeostasis. Should such polymorphisms be widespread, natural selection may be more effective at producing phenotypic change in suboptimal environments. In evaluating this hypothesis, a key unknown factor is the frequency with which HSP90-buffered polymorphisms occur in natural populations. Here, we present Arabidopsis thaliana populations suitable for genetic mapping that have constitutively reduced HSP90 levels. We employ quantitative genetic techniques to examine the HSP90-dependent polymorphisms affecting a host of plastic plant life-history traits. Our results demonstrate that HSP90-dependent natural variation is present at high frequencies in A. thaliana, with an expectation that at least one HSP90-dependent polymorphism will affect nearly every quantitative trait in progeny of two different wild lines. Hence, HSP90 is likely to occupy a central position in the translation of genotypic variation into phenotypic differences.