Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Comparative Mitogenome of Phylogenetic Relationships and Divergence Time Analysis within Potamanthidae (Insecta: Ephemeroptera)

Potamanthidae belongs to the superfamily Ephemeroidea but has no complete mt genome released in the NCBI (except for two unchecked and one partial mt genome). Since the sister clade to Potamanthidae has always been controversial, we sequenced seven mt genomes of Potamanthidae (two species from Rhoenanthus and five species from Potamanthus) in order to rebuild the phylogenetic relationships of Potamanthidae in this study. The divergence time of Potamanthidae was also investigated by utilizing five fossil calibration points because of the indeterminate origin time. In addition, because Rhoenanthus coreanus and Potamanthus luteus are always in low-temperature environments, we aimed to explore whether these two species were under positive selection at the mt genome level. Amongst the 13 PCGs, CGA was used as the start codon in COX1, whereas other genes conformed to initiating with an ATN start codon. From this analysis, UUA (L), AUU (I), and UUU (F) had the highest usage. Furthermore, the DHU arm was absent in the secondary structure of S1 in all species. By combining the 13 PCGs and 2 rRNAs, we reconstructed the phylogenetic relationship of Potamanthidae within Ephemeroptera. The monophyly of Potamanthidae and the monophyly of Rhoenanthus and Potamanthus were supported in the results. The phylogenetic relationship of Potamanthidae + (Ephemeridae + Polymitarcyidae) was also recovered with a high prior probability. The divergence times of Potamanthidae were traced to be 90.44 Mya (95% HPD, 62.80-121.74 Mya), and the divergence times of Rhoenanthus and Potamanthus originated at approximately 64.77 Mya (95% HPD, 43.82-88.68 Mya), thus belonging to the late Pliocene Epoch or early Miocene Epoch. In addition, the data indicated that R. coreanus was under negative selection and that ATP8 and ND2 in Potamanthidae had a high evolutionary rate.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

A delayed and unsynchronized ovary development as revealed by transcriptome of brain and pituitary of Coilia nasus

Coilia nasus is an anadromous fish that has been successfully domesticated in the last decade due to its high economic value. The fish exhibits a delayed ovary development during the reproductive season, despite breeding and selection for five to six offspring. The molecular mechanism of the delayed ovary development is still unknown, so the obstacles have not been removed in the large-scale breeding program. This study aims to investigate the key genes regulating ovarian development by comparing the transcriptomes of ovarian-stage IV and stage II brain/pituitary of Coilia nasus. Ovarian stages were validated by histological sections. A total of 75,097,641 and 66,735,592 high-quality reads were obtained from brain and pituitary transcriptomes, respectively, and alternatively spliced transcripts associated with gonadal development were detected. Compared to ovarian Ⅱ- brain, 515 differentially expressed genes (DEGs) were upregulated and 535 DEGs were downregulated in ovarian Ⅳ- brain, whereas 470 DEGs were upregulated and 483 DEGs were downregulated in ovarian Ⅳ- pituitary compared to ovarian Ⅱ- pituitary. DEGs involved in hormone synthesis and secretion and in the GnRH signaling pathway were screened. Weighted gene co-expression network analysis identified gene co-expression modules that were positively correlated with ovarian phenotypic traits. The hub genes Smad4 and TRPC4 in the modules were co-expressed with DEGs including Kiss1 receptor and JUNB, suggesting that ovarian development is controlled by a hypothalamic-pituitary-gonadal axis. Our results have provided new insights that advance our understanding of the molecular mechanism of C. nasus reproductive functions and will be useful for future breeding.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Exploring invertebrate indicators of ecosystem health by focusing on the flow transitional zones in a large, shallow eutrophic lake

The river-lake transitional zone provides a unique environment for the biological community and can reduce pollution inputs in lake ecosystems from their catchments. To explore environmental conditions with high purification potential in Lake Taihu and indicator species, we examined the river-to-lake changes in water and sediment quality and benthic invertebrate communities in the transitional zone of four regions. The spatial variations in the environment and invertebrate community observed in this study followed the previously reported patterns in Taihu; the northern and western regions were characterized by higher nutrient concentrations in water, higher heavy metal concentrations in sediment, and higher total invertebrate density and biomass dominated by pollution-tolerant oligochaetes and chironomids. Although nutrient concentrations were low and transparency was high in the eastern region, the taxon richness was the lowest there, which disagreed with the previous findings and might be due to a poor cover of macrophytes in this study. The river-to-lake change was large in the southern region for water quality and the invertebrate community. Water circulation induced by strong wind-wave actions in the lake sites of the southern region is assumed to have promoted photosynthetic and nutrient uptake activities and favored invertebrates that require well-aerated conditions such as polychaetes and burrowing crustaceans. Invertebrates usually adapted to brackish and saline environments are suggested to be indicators of a well-circulated environment with active biogeochemical processes and a less eutrophic state in Taihu, and wind-wave actions are key to maintaining such a community and natural purifying processes.

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

Selective effects of a short transient environmental fluctuation on a natural population

Natural populations experience continuous and often transient changes of environmental conditions. These in turn may result in fluctuating selection pressures leading to variable demographic and evolutionary population responses. Rapid adaptation as short-term response to a sudden environmental change has in several cases been attributed to polygenic traits, but the underlying genomic dynamics and architecture are poorly understood. In this study, we took advantage of a natural experiment in an insect population of the non-biting midge Chironomus riparius by monitoring genome-wide allele frequencies before and after a cold snap event. Whole genome pooled sequencing of time series samples revealed 10 selected haplotypes carrying ancient polymorphisms, partially with signatures of balancing selection. By constantly cold exposing genetically variable individuals in the laboratory, we could demonstrate with whole genome resequencing (i) that among the survivors, the same alleles rose in frequency as in the wild, and (ii) that the identified variants additively predicted fitness (survival time) of its bearers. Finally, by simultaneously sequencing the genome and the transcriptome of cold exposed individuals we could tentatively link some of the selected SNPs to the cis- and trans-regulation of genes and pathways known to be involved in cold response of insects, such as cytochrome P450 and fatty acid metabolism. Altogether, our results shed light on the strength and speed of selection in natural populations and the genomic architecture of its underlying polygenic trait. Population genomic time series data thus appear as promising tool for measuring the selective tracking of fluctuating selection in natural populations.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes

Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Reduced olfactory acuity in recently flightless insects suggests rapid regressive evolution

Background

Insects have exceptionally fast smelling capabilities, and some can track the temporal structure of odour plumes at rates above 100 Hz. It has been hypothesized that this fast smelling capability is an adaptation for flying. We test this hypothesis by comparing the olfactory acuity of sympatric flighted versus flightless lineages within a wing-polymorphic stonefly species.

Results

Our analyses of olfactory receptor neuron responses reveal that recently-evolved flightless lineages have reduced olfactory acuity. By comparing flighted versus flightless ecotypes with similar genetic backgrounds, we eliminate other confounding factors that might have affected the evolution of their olfactory reception mechanisms. Our detection of different patterns of reduced olfactory response strength and speed in independently wing-reduced lineages suggests parallel evolution of reduced olfactory acuity.

Conclusions

These reductions in olfactory acuity echo the rapid reduction of wings themselves, and represent an olfactory parallel to the convergent phenotypic shifts seen under selective gradients in other sensory systems (e.g. parallel loss of vision in cave fauna). Our study provides evidence for the hypothesis that flight poses a selective pressure on the speed and strength of olfactory receptor neuron responses and emphasizes the energetic costs of rapid olfaction.

Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles

Relative contributions of pre-existing vs de novo genomic variation to adaptation are poorly understood, especially in polyploid organisms. We assess this in high resolution using autotetraploid Arabidopsis arenosa, which repeatedly adapted to toxic serpentine soils that exhibit skewed elemental profiles. Leveraging a fivefold replicated serpentine invasion, we assess selection on SNPs and structural variants (TEs) in 78 resequenced individuals and discover significant parallelism in candidate genes involved in ion homeostasis. We further model parallel selection and infer repeated sweeps on a shared pool of variants in nearly all these loci, supporting theoretical expectations. A single striking exception is represented by TWO PORE CHANNEL 1, which exhibits convergent evolution from independent de novo mutations at an identical, otherwise conserved site at the calcium channel selectivity gate. Taken together, this suggests that polyploid populations can rapidly adapt to environmental extremes, calling on both pre-existing variation and novel polymorphisms.

Relative contributions of pre-existing versus de novo genomic variation to adaptation remain unclear. Here, the authors address this problem by examining the adaptation of autotetraploid Arabidopsis arenosa to serpentine soils and find that both types of variations contribute to rapid adaptation.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.