Parallel Speciation of Wild Rice Associated with Habitat Shifts

Genomics of hybrid parallel origin in Aquilegia ecalcarata

BACKGROUND

The parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear.

RESULTS

In this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata.

CONCLUSIONS

Our results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

Multiple domestications of Asian rice

The origin of domesticated Asian rice (Oryza sativa L.) has been controversial for more than half a century. The debates have focused on two leading hypotheses: a single domestication event in China or multiple domestication events in geographically separate areas. These two hypotheses differ in their predicted history of genes/alleles selected during domestication. Here we amassed a dataset of 1,578 resequenced genomes, including an expanded sample of wild rice from throughout its geographic range. We identified 993 selected genes that generated phylogenetic trees on which japonica and indica formed a monophyletic group, suggesting that the domestication alleles of these genes originated only once in either japonica or indica. Importantly, the domestication alleles of most selected genes (~80%) stemmed from wild rice in China, but the domestication alleles of a substantial minority of selected genes (~20%) originated from wild rice in South and Southeast Asia, demonstrating separate domestication events of Asian rice.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

Mitochondrial composition of and diffusion limiting factors of three social wasp genera Polistes, Ropalidia, and parapolybia (Hymenoptera: Vespidae)

Background

Social wasps Polistes, Ropalidia, and Parapolybia, belonging to the subfamily Polistinae, have obviously different distribution patterns, yet the factors leading to this difference remain unknown.

Results

The 17 newly sequenced mitogenomes of Polistes, Ropalidia, and Parapolybia contain 37 genes, and there are obvious differences among the compositions of the three genera. The monophyly of the genus Polistes and a monophyletic Ropalidiini: (Ropalidia + Parapolybia) are concordant with previous morphological analysis of the subfamily Polistinae. Our inferred divergence time demonstrates Polistes (at around 69 Ma) was diverged earlier than Ropalidia and Parapolybia (at around 61 Ma). The rearrangement of both trnY and trnL1 are shared by all the Polistinae. In addition, the unique rearrangement of TDRL derived at 69 Ma is detected in Polistes, and Ropalidia contains a Reversal which may derive at 61 Ma. Hereafter, the possibility is elaborated that Polistes originated in Aisa and then dispersed from Africa to South America, and Polistes and Ropalidia spread from Southeast Asia to Australia. At last, continental drift and Quaternary Ice Ages are inferred to be two main limiting factors in the current distributions of the three genera.

Conclusions

Obvious differences occur in the mitochondrial composition of Polistes, Ropalidia, and Parapolybia. According to the reconstructed time-calibrated framework, it is inquired that the continental drifts and the climate are mainly diffusion limiting factors of the three genera.

Supplementary information

The online version contains supplementary material available at 10.1186/s12862-022-02017-6.

Molecular signatures of parallel adaptive divergence causing reproductive isolation and speciation across two genera

Parallel evolution of reproductive isolation (PERI) provides strong evidence for natural selection playing a fundamental role in the origin of species. However, PERI has been rarely demonstrated for well established species drawn from different genera. In particular, parallel molecular signatures for the same genes in response to similar habitat divergence in such different lineages is lacking. Here, based on whole-genome sequencing data, we first explore the speciation process in two sister species of Carpinus (Betulaceae) in response to divergence for temperature and soil-iron concentration in habitats they occupy in northern and southwestern China, respectively. We then determine whether parallel molecular mutations occur during speciation in this pair of species and also in another sister-species pair of the related genus, Ostryopsis, which occupy similarly divergent habitats in China. We show that gene flow occurred during the origin of both pairs of sister species since approximately 9.8 or approximately 2 million years ago, implying strong natural selection during divergence. Also, in both species pairs we detected concurrent positive selection in a gene (LHY) for flowering time and in two paralogous genes (FRO4 and FRO7) of a gene family known to be important for iron tolerance. These changes were in addition to changes in other major genes related to these two traits. The different alleles of these particular candidate genes possessed by the sister species of Carpinus were functionally tested and indicated likely to alter flowering time and iron tolerance as previously demonstrated in the pair of Ostryopsis sister species. Allelic changes in these genes may have effectively resulted in high levels of prezygotic reproductive isolation to evolve between sister species of each pair. Our results show that PERI can occur in different genera at different timescales and involve similar signatures of molecular evolution at genes or paralogues of the same gene family, causing reproductive isolation as a consequence of adaptation to similarly divergent habitats.

•

PERI provides strong evidence for natural selection playing a fundamental role in the origin of species

•

PERI is rarely demonstrated for well-established species drawn from different genera

•

We detected PERI across two genera (Carpinus and Ostryopsis) in the family Betulaceae

•

PERI can occur in different genera at different timescales and involve molecular signatures at similar pathways

Evidence for two types of Aquilegia ecalcarata and its implications for adaptation to new environments

Spurs have played an important role in the radiation of the genus Aquilegia, but little is known about how the spurless state arose in A. ecalcarata. Here we aim to characterize the genetic divergence within A. ecalcarata and gain insights into the origin of this species. A total of 19 populations from A. ecalcarata and 23 populations from three of its closest relatives (Aquilegia kansuensis, Aquilegia rockii and Aquilegia yabeana) were sampled in this study. We sequenced fifteen nuclear gene fragments across the genome and three chloroplast loci to conduct phylogenetic, PCoA and STRUCTURE analyses. Our analyses indicate that A. ecalcarata may not be monophyletic and can be divided into two distinct lineages (A. ecalcarata I and A. ecalcarata II). A. ecalcarata I is genetically close to A. kansuensis, whereas A. ecalcarata II is close to A. rockii. Isolation-with-migration analysis suggested that historical gene flow was low between A. ecalcarata I and A. rockii, as well as between A. ecalcarata II and A. kansuensis. The two distinct lineages of A. ecalcarata show significant divergence in 13 floral traits and also have distinct distributions. In addition, both A. ecalcarata I and II are adapted to a stony environment that differs from that of their closest relatives, indicating a habitat shift may have driven new adaptations. Our findings enrich the understanding of how floral evolution contributes to species diversification.

Highly Replicated Evolution of Parapatric Ecotypes

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio

The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.

Phylogeographic analysis of Saxifraga fortunei complex (Saxifragaceae) reveals multiple origins of morphological and ecological variations in the Japanese Archipelago

Phenotypic polymorphism within a species is a notable phenomenon in evolutionary biology to understand the process of adaptive speciation and other historical events. The Saxifraga fortunei complex is a widespread herb found in East Asia. It includes several ecotypic taxa corresponding to their habitat environments. The distribution of the various ecotypes in a limited area of the Japanese Archipelago makes the species a suitable model to investigate the impact of population demographic history and natural selection on lineage diversification. Here, Sanger-based sequencing was used to estimate the divergence timeframe between populations of the Eurasian continent and Japan. Genome-wide SNPs obtained by ddRAD sequencing were used to investigate the phylogeographic origins of ecotypic taxa. The phylogenetic analyses revealed the divergence of the Japanese population from the continental population in the late Miocene. Two distinct regional clades of North and South Japan were identified; phenotypic diversification was evident only in the southern clade. The South Japan clades displayed a historical distribution expansion from north to south. The phenotypic variations appeared to have generated during the expansion. The ecotypic boundaries were incongruent with the genetic grouping. We propose that morphological and ecological specialization in Japanese populations was repeatedly generated by local natural selection.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara)

Background

Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown.

Results

Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding–noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding–noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding–noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare.

Conclusion

We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07754-4.

The landscape of gene-CDS-haplotype diversity in rice: Properties, population organization, footprints of domestication and breeding, and implications for genetic improvement

Polymorphisms within gene coding regions represent the most important part of the overall genetic diversity of rice. We characterized the gene-coding sequence-haplotype (gcHap) diversity of 45 963 rice genes in 3010 rice accessions. With an average of 226 ± 390 gcHaps per gene in rice populations, rice genes could be classified into three main categories: 12 865 conserved genes, 10 254 subspecific differentiating genes, and 22 844 remaining genes. We found that 39 218 rice genes carry >255 179 major gcHaps of potential functional importance. Most (87.5%) of the detected gcHaps were specific to subspecies or populations. The inferred proto-ancestors of local landrace populations reconstructed from conserved predominant (ancient) gcHaps correlated strongly with wild rice accessions from the same geographic regions, supporting a multiorigin (domestication) model of Oryza sativa. Past breeding efforts generally increased the gcHap diversity of modern varieties and caused significant frequency shifts in predominant gcHaps of 14 266 genes due to independent selection in the two subspecies. Low frequencies of "favorable" gcHaps at most known genes related to rice yield in modern varieties suggest huge potential for rice improvement by mining and pyramiding of favorable gcHaps. The gcHap data were demonstrated to have greater power than SNPs for the detection of causal genes that affect complex traits. The rice gcHap diversity dataset generated in this study would facilitate rice basic research and improvement in the future.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.

Divergence in flowering time is a major component contributing to reproductive isolation between two wild rice species (Oryza rufipogon and O. nivara)

It is of critical importance for our understanding of speciation process to determine the forms of reproductive isolation and their relative importance in species divergence. Oryza nivara and O. rufipogon are direct ancestors of Asian cultivated rice and a progenitor-daughter species pair. Investigating the reproductive isolation between them provides insights into plant speciation and helps understanding of the rice domestication. Here, we quantitatively measured the major components of reproductive isolation between the two species based on common garden and crossing experiments for three pairs of sympatric populations in Nepal, Cambodia and Laos. We revealed significant differences in the flowering times between species pairs, with O. nivara flowering much earlier than O. rufipogon. A very weak reduction in seed set but no reduction in F1 viability and fertility were detected for the crosses between species relative to those within species. Moreover, we detected asymmetrical compatibility between species and found that emasculation significantly decreased pollination success in O. nivara but not in O. rufipogon. Our study demonstrates that the divergence between O. nivara and O. rufipogon is maintained almost entirely by the difference in flowering times and suggests that differential flowering times contribute to both habitat preferences and reproductive isolation between species.

A chromosome-level genome assembly of the wild rice Oryza rufipogon facilitates tracing the origins of Asian cultivated rice

Oryza rufipogon Griff. is a wild progenitor of the Asian cultivated rice Oryza sativa. To better understand the genomic diversity of the wild rice, high-quality reference genomes of O. rufipogon populations are needed, which also facilitate utilization of the wild genetic resources in rice breeding. In this study, we generated a chromosome-level genome assembly of O. rufipogon using a combination of short-read sequencing, single-molecule sequencing, BioNano and Hi-C platforms. The genome sequence (399.8 Mb) was assembled into 46 scaffolds on the 12 chromosomes, with contig N50 and scaffold N50 of 13.2 Mb and 20.3 Mb, respectively. The genome contains 36,520 protein-coding genes, and 49.37% of the genome consists of repetitive elements. The genome has strong synteny with those of the O. sativa subspecies indica and japonica, but containing some large structural variations. Evolutionary analysis unveiled the polyphyletic origins of O. sativa, in which the japonica and indica genome formations involved different divergent O. rufipogon (including O. nivara) lineages, accompanied by introgression of genomic regions between japonica and indica. This high-quality reference genome provides insight on the genome evolution of the wild rice and the origins of the O. sativa subspecies, and valuable information for basic research and rice breeding.