From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation.

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Evolutionary Implications of the RNA N6-Methyladenosine Methylome in Plants.

Epigenetic modifications play important roles in genome evolution and innovation. However, most analyses have focused on the evolutionary role of DNA modifications, and little is understood about the influence of posttranscriptional RNA modifications on genome evolution. To explore the evolutionary significance of RNA modifications, we generated transcriptome-wide profiles of N6-methyladenosine (m6A), the most prevalent internal modification of mRNA, for 13 representative plant species spanning over half a billion years of evolution. These data reveal the evolutionary conservation and divergence of m6A methylomes in plants, uncover the preference of m6A modifications on ancient orthologous genes, and demonstrate less m6A divergence between orthologous gene pairs with earlier evolutionary origins. Further investigation revealed that the evolutionary divergence of m6A modifications is related to sequence variation between homologs from whole-genome duplication and gene family expansion from local-genome duplication. Unexpectedly, a significant negative correlation was found between the retention ratio of m6A modifications and the number of family members. Moreover, the divergence of m6A modifications is accompanied by variation in the expression level and translation efficiency of duplicated genes from whole- and local-genome duplication. Our work reveals new insights into evolutionary patterns of m6A methylomes in plant species and their implications, and provides a resource of plant m6A profiles for further studies of m6A regulation and function in an evolutionary context.

Evolution of a Human-Specific Tandem Repeat Associated with ALS.

Tandem repeats are proposed to contribute to human-specific traits, and more than 40 tandem repeat expansions are known to cause neurological disease. Here, we characterize a human-specific 69 bp variable number tandem repeat (VNTR) in the last intron of WDR7, which exhibits striking variability in both copy number and nucleotide composition, as revealed by long-read sequencing. In addition, greater repeat copy number is significantly enriched in three independent cohorts of individuals with sporadic amyotrophic lateral sclerosis (ALS). Each unit of the repeat forms a stem-loop structure with the potential to produce microRNAs, and the repeat RNA can aggregate when expressed in cells. We leveraged its remarkable sequence variability to align the repeat in 288 samples and uncover its mechanism of expansion. We found that the repeat expands in the 3'-5' direction, in groups of repeat units divisible by two. The expansion patterns we observed were consistent with duplication events, and a replication error called template switching. We also observed that the VNTR is expanded in both Denisovan and Neanderthal genomes but is fixed at one copy or fewer in non-human primates. Evaluating the repeat in 1000 Genomes Project samples reveals that some repeat segments are solely present or absent in certain geographic populations. The large size of the repeat unit in this VNTR, along with our multiplexed sequencing strategy, provides an unprecedented opportunity to study mechanisms of repeat expansion, and a framework for evaluating the roles of VNTRs in human evolution and disease.

Tek (Tie2) is not required for cardiovascular development in zebrafish.

Angiopoietin/TIE signalling plays a major role in blood and lymphatic vessel development. In mouse, Tek (previously known as Tie2) mutants die prenatally due to a severely underdeveloped cardiovascular system. In contrast, in zebrafish, previous studies have reported that although embryos injected with tek morpholinos (MOs) exhibit severe vascular defects, tek mutants display no obvious vascular malformations. To further investigate the function of zebrafish Tek, we generated a panel of loss-of-function tek mutants, including RNA-less alleles, an allele lacking the MO-binding site, an in-frame deletion allele and a premature termination codon-containing allele. Our data show that all these mutants survive to adulthood with no obvious cardiovascular defects. MO injections into tek mutants lacking the MO-binding site or the entire tek locus cause similar vascular defects to those observed in MO-injected +/+ siblings, indicating off-target effects of the MOs. Surprisingly, comprehensive phylogenetic profiling and synteny analyses reveal that Tek was lost in the largest teleost clade, suggesting a lineage-specific shift in the function of TEK during vertebrate evolution. Altogether, these data show that Tek is dispensable for zebrafish development, and probably dispensable in most teleost species.

Assortative mating and mate-choice contributes to the maintenance of a developmental dimorphism in Streblospio benedicti.

Assortative mating, where individuals non-randomly mate with respect to phenotype or genotype, can occur when preferences between potential mates have evolved. When such mate preferences occur in a population it can drive evolutionary and phenotypic divergence. But the extent to which assortative mating, mate preference, and development are evolutionarily linked remains unclear. Here we use Streblospio benedicti, a marine annelid with a rare developmental dimorphism, to investigate if mate-choice could contribute to developmental evolution. For S. benedicti two types of ecologically and phenotypically similar adults persist in natural populations, but they give rise to distinctly different offspring with alternative life-histories. This dimorphism persists despite the absence of post-zygotic reproductive barriers, where crosses between the developmental types can produce phenotypically intermediate offspring. How this life-history strategy evolved remains unknown, but assortative mating is a typical first step in evolutionary divergence. Here we investigate if female mate-choice is occurring in this species. We find that mate preferences could be contributing to the maintenance of alternative developmental and life-history strategies.

Soft selection reduces loss of heterozygosity in asexual reproduction.

The adaptive value of sexual reproduction is still debated in evolutionary theory. It has been proposed that the advantage of sexual reproduction over asexual reproduction is to promote genetic diversity, to prevent the accumulation of harmful mutations or to preserve heterozygosity. Since these hypothetical advantages depend on the type of asexual reproduction, understanding how selection affects the taxonomic distribution of each type could help us discriminate between existing hypotheses. Here, I argue that soft selection, competition among embryos or offspring in selection arenas prior to the hard selection of the adult phase, reduces loss of heterozygosity in certain types of asexual reproduction. Since loss of heterozygosity leads to the unmasking of recessive deleterious mutations in the progeny of asexual individuals, soft selection facilitates the evolution of these types of asexual reproduction. Using a population genetics model, I calculate how loss of heterozygosity affects fitness for different types of apomixis and automixis, and I show that soft selection significantly reduces loss of heterozygosity, hence increases fitness, in apomixis with suppression of the first meiotic division and in automixis with central fusion, the most common types of asexual reproduction. Therefore, if sexual reproduction evolved to preserve heterozygosity, soft selection should be associated with these types of asexual reproduction. I discuss the evidence for this prediction and how this and other observations on the distribution of different types of asexual reproduction in nature is consistent with the heterozygosity hypothesis.

Understanding the contemporary high obesity rate from an evolutionary genetic perspective.

The topic of obesity is gaining increasing popularity globally. From an evolutionary genetic perspective, it is believed that the main cause of the high obesity rate is the mismatch between environment and genes after people have shifted toward a modern high-calorie diet. However, it has been debated for over 60 years about how obesity-related genes become prevalent all over the world. Here, we review the three most influential hypotheses or viewpoints, i.e., the thrifty gene hypothesis, the drifty gene hypothesis, and the maladaptation viewpoint. In particular, genome-wide association studies in the recent 10 years have provided rich findings and evidence to be considered for a better understanding of the evolutionary genetic mechanisms of obesity. We anticipate this brief review to direct further studies and inspire the future application of precision medicine in obesity treatment.

Genomic Instability Evolutionary Footprints on Human Health: Driving Forces or Side Effects?

In this work, we propose a comprehensive perspective on genomic instability comprising not only the accumulation of mutations but also telomeric shortening, epigenetic alterations and other mechanisms that could contribute to genomic information conservation or corruption. First, we present mechanisms playing a role in genomic instability across the kingdoms of life. Then, we explore the impact of genomic instability on the human being across its evolutionary history and on present-day human health, with a particular focus on aging and complex disorders. Finally, we discuss the role of non-coding RNAs, highlighting future approaches for a better living and an expanded healthy lifespan.

Geonomics: Forward-Time, Spatially Explicit, and Arbitrarily Complex Landscape Genomic Simulations.

Understanding the drivers of spatial patterns of genomic diversity has emerged as a major goal of evolutionary genetics. The flexibility of forward-time simulation makes it especially valuable for these efforts, allowing for the simulation of arbitrarily complex scenarios in a way that mimics how real populations evolve. Here, we present Geonomics, a Python package for performing complex, spatially explicit, landscape genomic simulations with full spatial pedigrees that dramatically reduces user workload yet remains customizable and extensible because it is embedded within a popular, general-purpose language. We show that Geonomics results are consistent with expectations for a variety of validation tests based on classic models in population genetics and then demonstrate its utility and flexibility with a trio of more complex simulation scenarios that feature polygenic selection, selection on multiple traits, simulation on complex landscapes, and nonstationary environmental change. We then discuss runtime, which is primarily sensitive to landscape raster size, memory usage, which is primarily sensitive to maximum population size and recombination rate, and other caveats related to the model's methods for approximating recombination and movement. Taken together, our tests and demonstrations show that Geonomics provides an efficient and robust platform for population genomic simulations that capture complex spatial and evolutionary dynamics.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster.

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Course-based undergraduate research experiences (CURES) as a pathway to diversify science.

There is widespread recognition of the need to increase research opportunities in biomedical science for undergraduate students from underrepresented backgrounds. Here, we describe the implementation of team-based science combined with intensive mentoring to conduct a large-scale project examining the evolution of behavior. This system can be widely applied in other areas of STEM to promote research-intensive opportunities in STEM fields and to promote diversity in science.

Background selection under evolving recombination rates.

Background selection (BGS), the effect that purifying selection exerts on sites linked to deleterious alleles, is expected to be ubiquitous across eukaryotic genomes. The effects of BGS reflect the interplay of the rates and fitness effects of deleterious mutations with recombination. A fundamental assumption of BGS models is that recombination rates are invariant over time. However, in some lineages, recombination rates evolve rapidly, violating this central assumption. Here, we investigate how recombination rate evolution affects genetic variation under BGS. We show that recombination rate evolution modifies the effects of BGS in a manner similar to a localized change in the effective population size, potentially leading to underestimation or overestimation of the genome-wide effects of selection. Furthermore, we find evidence that recombination rate evolution in the ancestors of modern house mice may have impacted inferences of the genome-wide effects of selection in that species.

The road not taken: Evolution of tetrodotoxin resistance in the Sierra garter snake (Thamnophis couchii) by a path less travelled.

The repeated evolution of tetrodotoxin (TTX) resistance provides a model for testing hypotheses about the mechanisms of convergent evolution. This poison is broadly employed as a potent antipredator defence, blocking voltage-gated sodium channels (Nav ) in muscles and nerves, paralysing and sometimes killing predators. Resistance in taxa bearing this neurotoxin and a few predators appears to come from convergent replacements in specific Nav residues that interact with TTX. This stereotyped genetic response suggests molecular and phenotypic evolution may be constrained and predictable. Here, we investigate the extent of mechanistic convergence in garter snakes (Thamnophis) that prey on TTX-bearing newts (Taricha) by examining the physiological and genetic basis of TTX resistance in the Sierra garter snake (Th. couchii). We characterize variation in this predatory adaptation across populations at several biological scales: whole-animal TTX resistance; skeletal muscle resistance; functional genetic variation in three Nav encoding loci; and levels of gene expression for one of these loci. We found Th. couchii possess extensive geographical variation in resistance at the whole-animal and skeletal muscle levels. As in other Thamnophis, resistance at both levels is highly correlated, suggesting convergence across the biological levels linking organism to organ. However, Th. couchii shows no functional variation in Nav loci among populations or difference in candidate gene expression. Local variation in TTX resistance in Th. couchii cannot be explained by the same relationship between genotype and phenotype seen in other taxa. Thus, historical contingencies may lead different species of Thamnophis down alternative routes to local adaptation.

Daisy-chain gene drives for the alteration of local populations.

If they are able to spread in wild populations, CRISPR-based gene-drive elements would provide new ways to address ecological problems by altering the traits of wild organisms, but the potential for uncontrolled spread tremendously complicates ethical development and use. Here, we detail a self-exhausting form of CRISPR-based drive system comprising genetic elements arranged in a daisy chain such that each drives the next. "Daisy-drive" systems can locally duplicate any effect achievable by using an equivalent self-propagating drive system, but their capacity to spread is limited by the successive loss of nondriving elements from one end of the chain. Releasing daisy-drive organisms constituting a small fraction of the local wild population can drive a useful genetic element nearly to local fixation for a wide range of fitness parameters without self-propagating spread. We additionally report numerous highly active guide RNA sequences sharing minimal homology that may enable evolutionarily stable daisy drive as well as self-propagating CRISPR-based gene drive. Especially when combined with threshold dependence, daisy drives could simplify decision-making and promote ethical use by enabling local communities to decide whether, when, and how to alter local ecosystems.

Fitness effects but no temperature-mediated balancing selection at the polymorphic Adh gene of Drosophila melanogaster.

Polymorphism in the alcohol dehydrogenase (ADH) protein of Drosophila melanogaster, like genetic variation in many other enzymes, has long been hypothesized to be maintained by a selective trade-off between thermostability and enzyme activity. Two major Adh variants, named Fast and Slow, are distributed along latitudinal clines on several continents. The balancing selection trade-off hypothesis posits that Fast is favored at high latitudes because it metabolizes alcohol faster, whereas Slow is favored at low latitudes because it is more stable at high temperatures. Here we use biochemical and physiological assays of precisely engineered genetic variants to directly test this hypothesis. As predicted, the Fast protein has higher catalytic activity than Slow, and both the Fast protein and regulatory variants linked to it confer greater ethanol tolerance on transgenic animals. But we found no evidence of a temperature-mediated trade-off: The Fast protein is not less stable or active at high temperatures, and Fast alleles increase ethanol tolerance and survivorship at all temperatures tested. Further, analysis of a population genomic dataset reveals no signature of balancing selection in the Adh gene. These results provide strong evidence against balancing selection driven by a stability/activity trade-off in Adh, and they justify caution about this hypothesis for other enzymes except those for which it has been directly tested. Our findings tentatively suggest that environment-specific selection for the Fast allele, coupled with demographic history, may have produced the observed pattern of Adh variation.

Comparing signals of natural selection between three Indigenous North American populations.

While many studies have highlighted human adaptations to diverse environments worldwide, genomic studies of natural selection in Indigenous populations in the Americas have been absent from this literature until very recently. Since humans first entered the Americas some 20,000 years ago, they have settled in many new environments across the continent. This diversity of environments has placed variable selective pressures on the populations living in each region, but the effects of these pressures have not been extensively studied to date. To help fill this gap, we collected genome-wide data from three Indigenous North American populations from different geographic regions of the continent (Alaska, southeastern United States, and central Mexico). We identified signals of natural selection in each population and compared signals across populations to explore the differences in selective pressures among the three regions sampled. We find evidence of adaptation to cold and high-latitude environments in Alaska, while in the southeastern United States and central Mexico, pathogenic environments seem to have created important selective pressures. This study lays the foundation for additional functional and phenotypic work on possible adaptations to varied environments during the history of population diversification in the Americas.

A hapless mathematical contribution to biology : Chromosome inversions in Drosophila, 1937-1941.

This is the story, told in the light of a new analysis of historical data, of a mathematical biology problem that was explored in the 1930s in Thomas Morgan's laboratory at the California Institute of Technology. It is one of the early developments of evolutionary genetics and quantitative phylogeny, and deals with the identification and counting of chromosomal inversions in Drosophila species from comparisons of genetic maps. A re-analysis of the data produced in the 1930s using current mathematics and computational technologies reveals how a team of biologists, with the help of a renowned mathematician and against their first intuition, came to an erroneous conclusion regarding the presence of phylogenetic signals in gene arrangements. This example illustrates two different aspects of a same piece: (1) the appearance of a mathematical in biology problem solved with the development of a combinatorial algorithm, which was unusual at the time, and (2) the role of errors in scientific activity. Also underlying is the possible influence of computational complexity in understanding the directions of research in biology.

Fine-Scale Haplotype Structure Reveals Strong Signatures of Positive Selection in a Recombining Bacterial Pathogen.

Identifying genetic variation in bacteria that has been shaped by ecological differences remains an important challenge. For recombining bacteria, the sign and strength of linkage provide a unique lens into ongoing selection. We show that derived alleles <300 bp apart in Neisseria gonorrhoeae exhibit more coupling linkage than repulsion linkage, a pattern that cannot be explained by limited recombination or neutrality as these couplings are significantly stronger for nonsynonymous alleles than synonymous alleles. This general pattern is driven by a small fraction of highly diverse genes, many of which exhibit evidence of interspecies horizontal gene transfer and an excess of intermediate frequency alleles. Extensive simulations show that two distinct forms of positive selection can create these patterns of genetic variation: directional selection on horizontally transferred alleles or balancing selection that maintains distinct haplotypes in the presence of recombination. Our results establish a framework for identifying patterns of selection in fine-scale haplotype structure that indicate specific ecological processes in species that recombine with distantly related lineages or possess coexisting adaptive haplotypes.

DNA Methylation Profiles Suggest Intergenerational Transfer of Maternal Effects.

The view of maternal effects (nongenetic maternal environmental influence on offspring phenotype) has changed from one of distracting complications in evolutionary genetics to an important evolutionary mechanism for improving offspring fitness. Recent studies have shown that maternal effects act as an adaptive mechanism to prepare offspring for stressful environments. Although research into the magnitude of maternal effects is abundant, the molecular mechanisms of maternal influences on offspring phenotypic variation are not fully understood. Despite recent work identifying DNA methylation as a potential mechanism of nongenetic inheritance, currently proposed links between DNA methylation and parental effects are indirect and primarily involve genomic imprinting. We combined a factorial breeding design and gene-targeted sequencing methods to assess inheritance of methylation during early life stages at 14 genes involved in growth, development, metabolism, stress response, and immune function of Chinook salmon (Oncorhynchus tshawytscha). We found little evidence for additive or nonadditive genetic effects acting on methylation levels during early development; however, we detected significant maternal effects. Consistent with conventional maternal effect data, maternal effects on methylation declined through development and were replaced with nonadditive effects when offspring began exogenous feeding. We mapped methylation at individual CpG sites across the selected candidate genes to test for variation in site-specific methylation profiles and found significant maternal effects at selected CpG sites that also declined with development stage. While intergenerational inheritance of methylated DNA is controversial, we show that CpG-specific methylation may function as an underlying molecular mechanism for maternal effects, with important implications for offspring fitness.

Inference of Chromosome-Length Haplotypes Using Genomic Data of Three or a Few More Single Gametes.

Compared with genomic data of individual markers, haplotype data provide higher resolution for DNA variants, advancing our knowledge in genetics and evolution. Although many computational and experimental phasing methods have been developed for analyzing diploid genomes, it remains challenging to reconstruct chromosome-scale haplotypes at low cost, which constrains the utility of this valuable genetic resource. Gamete cells, the natural packaging of haploid complements, are ideal materials for phasing entire chromosomes because the majority of the haplotypic allele combinations has been preserved. Therefore, compared with the current diploid-based phasing methods, using haploid genomic data of single gametes may substantially reduce the complexity in inferring the donor's chromosomal haplotypes. In this study, we developed the first easy-to-use R package, Hapi, for inferring chromosome-length haplotypes of individual diploid genomes with only a few gametes. Hapi outperformed other phasing methods when analyzing both simulated and real single gamete cell sequencing data sets. The results also suggested that chromosome-scale haplotypes may be inferred by using as few as three gametes, which has pushed the boundary to its possible limit. The single gamete cell sequencing technology allied with the cost-effective Hapi method will make large-scale haplotype-based genetic studies feasible and affordable, promoting the use of haplotype data in a wide range of research.