The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

Genomic Divergence and the Evolution of Ecotypes in Bottlenose Dolphins (Genus Tursiops)

Climatic changes have caused major environmental restructuring throughout the world's oceans. Marine organisms have responded to novel conditions through various biological systems, including genomic adaptation. Growing accessibility of next-generation DNA sequencing methods to study nonmodel species has recently allowed genomic changes underlying environmental adaptations to be investigated. This study used double-digest restriction-site associated DNA (ddRAD) sequence data to investigate the genomic basis of ecotype formation across currently recognized species and subspecies of bottlenose dolphins (genus Tursiops) in the Southern Hemisphere. Subspecies-level genomic divergence was confirmed between the offshore common bottlenose dolphin (T. truncatus truncatus) and the inshore Lahille's bottlenose dolphin (T. t. gephyreus) from the southwestern Atlantic Ocean (SWAO). Similarly, subspecies-level divergence is suggested between inshore (eastern Australia) Indo-Pacific bottlenose dolphin (T. aduncus) and the proposed Burrunan dolphin (T. australis) from southern Australia. Inshore bottlenose dolphin lineages generally had lower genomic diversity than offshore lineages, a pattern particularly evident for T. t. gephyreus, which showed exceptionally low diversity. Genomic regions associated with cardiovascular, musculoskeletal, and energy production systems appear to have undergone repeated adaptive evolution in inshore lineages across the Southern Hemisphere. We hypothesize that comparable selective pressures in the inshore environment drove similar adaptive responses in each lineage, supporting parallel evolution of inshore bottlenose dolphins. With climate change altering marine ecosystems worldwide, it is crucial to gain an understanding of the adaptive capacity of local species and populations. Our study provides insights into key adaptive pathways that may be important for the long-term survival of cetaceans and other organisms in a changing marine environment.

Spatially and temporally varying selection influence species boundaries in two sympatric Mimulus

Spatially and temporally varying selection can maintain genetic variation within and between populations, but it is less well known how these forces influence divergence between closely related species. We identify the interaction of temporal and spatial variation in selection and their role in either reinforcing or eroding divergence between two closely related Mimulus species. Using repeated reciprocal transplant experiments with advanced generation hybrids, we compare the strength of selection on quantitative traits involved in adaptation and reproductive isolation in Mimulus guttatus and Mimulus laciniatus between two years with dramatically different water availability. We found strong divergent habitat-mediated selection on traits in the direction of species differences during a drought in 2013, suggesting that spatially varying selection maintains species divergence. However, a relaxation in divergent selection on most traits in an unusually wet year (2019), including flowering time, which is involved in pre-zygotic isolation, suggests that temporal variation in selection may weaken species differences. Therefore, we find evidence that temporally and spatially varying selection may have opposing roles in mediating species boundaries. Given our changing climate, future growing seasons are expected to be more similar to the dry year, suggesting that in this system climate change may actually increase species divergence.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

Large scale DNA barcoding of the subfamily Culterinae (Cypriniformes: Xenocyprididae) in East Asia unveils geographic scale effect, taxonomic warnings and cryptic diversity

Geographical scale might be expected to impact significantly the efficiency of DNA barcoding as spatially comprehensive sampling provides opportunities to uncover intricate relationships among closely related species and to detect cryptic diversity for widespread taxa. Here, we present a DNA barcoding study on a Xencyprididae subfamily (Culterinae) involving the production of 998 newly generated DNA barcodes from East Asian drainages (80 localities). Together with 513 barcodes mined from BOLD and GenBank, a reference library consisting of 1511 DNA barcodes (116 localities) for 42 species was assembled, accounting for 66% of known Culterinae species. Intraspecific genetic distances are positively correlated to geographical scale, while a negative correlation is detected between interspecific genetic distances and geographical scale. The present study demonstrates that geographical scale influences the efficiency of DNA barcoding by narrowing the width of the barcoding gap. DNA‐based species delimitation analyses delimited 44 molecular operational taxonomic units (MOTUs). Rampant cryptic diversity is detected within eight species with multiple MOTUs, whereas 25 species present mismatch between morphological and molecular delimitations. A total of 18 species are lumped into nine MOTUs due to low interspecific divergence and/or mixed lineages. Several MOTU divergences are hypothesized to relate to known biogeographical barriers and geological events during the Pliocene and Pleistocene. This study provides new insights into the taxonomy and phylogeography of the subfamily Culterinae.

Divergent selection on flowering phenology but not on floral morphology between two closely related orchids

Closely related species often differ in traits that influence reproductive success, suggesting that divergent selection on such traits contribute to the maintenance of species boundaries. Gymnadenia conopsea ss. and Gymnadenia densiflora are two closely related, perennial orchid species that differ in (a) floral traits important for pollination, including flowering phenology, floral display, and spur length, and (b) dominant pollinators. If plant-pollinator interactions contribute to the maintenance of trait differences between these two taxa, we expect current divergent selection on flowering phenology and floral morphology between the two species. We quantified phenotypic selection via female fitness in one year on flowering start, three floral display traits (plant height, number of flowers, and corolla size) and spur length, in six populations of G. conopsea s.s. and in four populations of G. densiflora. There was indication of divergent selection on flowering start in the expected direction, with selection for earlier flowering in two populations of the early-flowering G. conopsea s.s. and for later flowering in one population of the late-flowering G. densiflora. No divergent selection on floral morphology was detected, and there was no significant stabilizing selection on any trait in the two species. The results suggest ongoing adaptive differentiation of flowering phenology, strengthening this premating reproductive barrier between the two species. Synthesis: This study is among the first to test whether divergent selection on floral traits contribute to the maintenance of species differences between closely related plants. Phenological isolation confers a substantial potential for reproductive isolation, and divergent selection on flowering time can thus greatly influence reproductive isolation and adaptive differentiation.

Analysis of Genome-Wide Differentiation between Native and Introduced Populations of the Cupped Oysters Crassostrea gigas and Crassostrea angulata

The Pacific cupped oyster is genetically subdivided into two sister taxa, Crassostrea gigas and Crassostrea angulata, which are in contact in the north-western Pacific. The nature and origin of their genetic and taxonomic differentiation remains controversial due the lack of known reproductive barriers and the high degree of morphologic similarity. In particular, whether the presence of ecological and/or intrinsic isolating mechanisms contributes to species divergence is unknown. The recent co-introduction of both taxa into Europe offers a unique opportunity to test how genetic differentiation is maintained under new environmental and demographic conditions. We generated a pseudochromosome assembly of the Pacific oyster genome using a combination of BAC-end sequencing and scaffold anchoring to a new high-density linkage map. We characterized genome-wide differentiation between C. angulata and C. gigas in both their native and introduced ranges, and showed that gene flow between species has been facilitated by their recent co-introductions in Europe. Nevertheless, patterns of genomic divergence between species remain highly similar in Asia and Europe, suggesting that the environmental transition caused by the co-introduction of the two species did not affect the genomic architecture of their partial reproductive isolation. Increased genetic differentiation was preferentially found in regions of low recombination. Using historical demographic inference, we show that the heterogeneity of differentiation across the genome is well explained by a scenario whereby recent gene flow has eroded past differentiation at different rates across the genome after a period of geographical isolation. Our results thus support the view that low-recombining regions help in maintaining intrinsic genetic differences between the two species.

Phylogeography of two closely related species of Allium endemic to East Asia: Population evolution in response to climate oscillations

This study investigated the effects of climate oscillations on the evolution of two closely related Allium species, A. neriniflorum and A. tubiflorum. We sequenced three cp DNA (cpDNA) fragments (rps16, rpl32-trnL, and trnD-trnT, together approximately 2,500 bp in length) of two closely related Allium species, with samples from 367 individuals in 47 populations distributed across the total range of these species. The interspecific and intraspecific divergence times of the two species were in the Quaternary glaciation. The population divergence was high for the cpDNA variation, suggesting a significant phylogeographic structure (NST = 0.844, GST = 0.798, p < 0.05). Remarkable ecological differentiation was also revealed by Niche models and statistical analyses. Our results suggest the speciation event of the two species was triggered by violent climatic changes during the Quaternary glaciation.

A comprehensive genomic history of extinct and living elephants

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

Introgression between divergent corn borer species in a region of sympatry: Implications on the evolution and adaptation of pest arthropods

The Asian corn borer, Ostrinia furnacalis, and European corn borer, O. nubilalis (Lepidoptera: Crambidae), cause damage to cultivated maize in spatially distinct geographies and have evolved divergent hydrocarbons as the basis of sexual communication. The Yili area of Xinjiang Uyghur Autonomous Region in China represents the only known region where O. furnacalis has invaded a native O. nubilalis range, and these two corn borer species have made secondary contact. Genetic differentiation was estimated between Ostrinia larvae collected from maize plants at 11 locations in Xinjiang and genotyped using high-throughput SNP and microsatellite markers. Maternal lineages were assessed by direct sequencing of mitochondrial cytochrome c oxidase subunit I and II haplotypes, and a high degree of genotypic diversity was demonstrated between lineages based on SNP genotypes. Furthermore, historical introgression was predicted among SNP genotypes only at sympatric locations in the Yili area, whereas in Xinjiang populations only O. furnacalis haplotypes were detected and no analogous introgressed genotypes were predicted. Our detection of putative hybrids and historical evidence of introgression defines Yili area as a hybrid zone between the species in normal ecological interactions and furthermore, might indicate that adaptive traits could spread even between seemingly divergent species through horizontal transmission. Results of this study indicate there may be a continuum in the degree of reproductive isolation between Ostrinia species and that the elegance of distinct and complete speciation based on modifications to the pheromone communication might need to be reconsidered.

Effects of Mountain Uplift and Climatic Oscillations on Phylogeography and Species Divergence in Four Endangered Notopterygium Herbs

Mountain uplift and climatic fluctuations are important driving forces that have affected the geographic distribution and population dynamics history of organisms. However, it is unclear how geological and climatic events might have affected the phylogeographic history and species divergence in high-alpine herbal plants. In this study, we analyzed the population demographic history and species differentiation of four endangered Notopterygium herbs on the high-altitude Qinghai-Tibetan Plateau (QTP) and adjacent areas. We combined phylogeographic analysis with species distribution modeling to detect the genetic variations in four Notopterygium species (N. incisum, N. franchetii, N. oviforme, and N. forrestii). In total, 559 individuals from 74 populations of the four species were analyzed based on three maternally inherited chloroplast fragments (matK, rbcL, and trnS-trnG) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-five chloroplast DNA (cpDNA) and 48 ITS haplotypes were identified in the four species. All of the cpDNA and ITS haplotypes were species-specific, except N. franchetii and N. oviforme shared one cpDNA haplotype, H32. Phylogenetic analysis suggested that all four species formed a monophyletic clade with high bootstrap support, where N. franchetii and N. oviforme were sisters. In addition, each Notopterygium species generated an individual clade that corresponded to their respective species in the ITS tree. Population dynamics analyses and species distribution modeling showed that the two widely distributed herbs N. incisum and N. franchetii exhibited obvious demographic expansions during the Pleistocene ice ages. Molecular dating suggested that the divergence of the four Notopterygium species occurred approximately between 3.6 and 1.2 Mya, and it was significantly associated with recent extensive uplifts of the QTP. Our results support the hypothesis that mountain uplift and Quaternary climatic oscillations profoundly shaped the population genetic divergence and demographic dynamics of Notopterygium species. The findings of this and previous studies provide important insights into the effects of QTP uplifts and climatic changes on phylogeography and species differentiation in high altitude mountainous areas. Our results may also facilitate the conservation of endangered herbaceous medicinal plants in the genus Notopterygium.

Inferring species divergence times using pairwise sequential Markovian coalescent modelling and low-coverage genomic data

Understanding when species diverged aids in identifying the drivers of speciation, but the end of gene flow between populations can be difficult to ascertain from genetic data. We explore the use of pairwise sequential Markovian coalescent (PSMC) modelling to infer the timing of divergence between species and populations. PSMC plots generated using artificial hybrid genomes show rapid increases in effective population size at the time when the two parent lineages diverge, and this approach has been used previously to infer divergence between human lineages. We show that, even without high coverage or phased input data, PSMC can detect the end of significant gene flow between populations by comparing the PSMC output from artificial hybrids to the output of simulations with known demographic histories. We then apply PSMC to detect divergence times among lineages within two real datasets: great apes and bears within the genus Ursus Our results confirm most previously proposed divergence times for these lineages, and suggest that gene flow between recently diverged lineages may have been common among bears and great apes, including up to one million years of continued gene flow between chimpanzees and bonobos after the formation of the Congo River.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Molecular phylogeny, biogeography and ecological niche modelling of Cardiocrinum (Liliaceae): insights into the evolutionary history of endemic genera distributed across the Sino-Japanese floristic region

BACKGROUND AND AIMS

The patterns of evolutionary assembly in the Sino-Japanese floristic region (SJFR) remain largely unknown due to a lack of integrative multidimensional studies throughout the region. To address this issue, we elucidated the evolutionary history of Cardiocrinum (Liliaceae), a genus containing four taxa distributed across the SJFR.

METHODS

Fifty-four populations were sampled throughout the geographical range of Cardiocrinum to assess genetic structure, analyse phylogenetic relationships and reconstruct ancestral area based on six chloroplast DNA (cpDNA) fragments and three low copy nuclear genes (LCNG). Ecological niche modelling was used to examine the potential range shifts of Cardiocrinum in response to climatic change.

KEY RESULTS

The molecular data showed high genetic similarity in the cpDNA (98·37 %) and LCNG (94·53 %) sequences. The biogeographical analyses revealed that the ancestor of Cardiocrinum diversified during the late Miocene (approx. 7·32 Mya) in Central China. The ancestor of the C. giganteum lineage dispersed westward to the Himalayas and south-west China with the split between C. giganteum and C. giganteum var. yunnanense occurring around 4·11 Mya consistent with the period of orogeny of the Hengduan Mountains. Some populations of the C. cathayanum lineage dispersed eastward to south Japan via the land bridge approx. 4·97 Mya, providing opportunities for allopatric speciation of C. cordatum The predicted suitable habitats of Cardiocrinum have become smaller and more fragmented since the Last Glacial Maximum.

CONCLUSIONS

Our study provides evidence of a biogeographical pattern of dispersal from Central China to the Himalayas in the west and Japan in the east for genera distributed across the SJFR, and highlights that the orogeny of the Hengduan Mountains and fluctuations of the sea level of the East China Sea played important roles in promoting species divergence.

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai-Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.

Plant regeneration from seeds responds to phylogenetic relatedness and local adaptation in Mediterranean Romulea (Iridaceae) species

Seed germination is the most important transitional event between early stages in the life cycle of spermatophytes and understanding it is crucial to understand plant adaptation and evolution. However, so far seed germination of phylogenetically closely related species has been poorly investigated. To test the hypothises that phylogenetically related plant species have similar seed ecophysiological traits thereby reflecting certain habitat conditions as a result of local adaptation, we studied seed dormancy and germination in seven Mediterranean species in the genus Romulea (Iridaceae). Both the across-species model and the model accounting for shared evolutionary history showed that cool temperatures (≤ 15°C) were the main factor that promoted seed germination. The absence of embryo growth before radicle emergence is consistent with a prompt germination response at cool temperatures. The range of temperature conditions for germination became wider after a period of warm stratification, denoting a weak primary dormancy. Altogether these results indicate that the studied species exhibit a Mediterranean germination syndrome, but with species-specific germination requirements clustered in a way that follows the phylogenetic relatedness among those species. In addition, species with heavier seeds from humid habitats showed a wider range of conditions for germination at dispersal time than species from dry habitats possessing lighter seeds. We conclude that while phylogenetically related species showed very similar germination requirements, there are subtle ecologically meaningful differences, confirming the onset of adaptation to local ecological factors mediated by species relatedness.

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai-Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai-Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD-psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 - 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.

Complete mitochondrial genome of Chinese bamboo rat, Rhizomys sinensis and species divergence comparison

The Chinese bamboo rat Rhizomys sinensis belongs to family Spalacidae, and is distributed in China, Myanmar, and Vietnam. In this study, the entire mitochondrial genome of R. sinensis was firstly determined. The genome is 16,564 bases in length, containing 13 protein coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and a putative control region. The composition and arrangement of its genes are identical to most other mammals. The whole base composition of the mitogenome is A 31.9%, G 12.0%, T 31.2% and C 24.9%, with an A+T rich pattern. The species divergence between R. sinensis and R. pruinosus is 0.113, in accordance with a "10× rule", as the intraspecific barcode variation averaged 1%. The mitogenome data of R. sumatrensis is required to better understand the phylogenetic relationships within Rhizomys.

Polymorphism and Divergence in Two Willow Species, Salix viminalis L. and Salix schwerinii E. Wolf

We investigated species divergence, present and past gene flow, levels of nucleotide polymorphism, and linkage disequilibrium in two willows from the plant genus Salix. Salix belongs together with Populus to the Salicaceae family; however, most population genetic studies of Salicaceae have been performed in Populus, the model genus in forest biology. Here we present a study on two closely related willow species Salix viminalis and S. schwerinii, in which we have resequenced 33 and 32 nuclear gene segments representing parts of 18 nuclear loci in 24 individuals for each species. We used coalescent simulations and estimated the split time to around 600,000 years ago and found that there is currently limited gene flow between the species. Mean intronic nucleotide diversity across gene segments was slightly higher in S. schwerinii (π_i = 0.00849) than in S. viminalis (π_i = 0.00655). Compared with other angiosperm trees, the two willows harbor intermediate levels of silent polymorphisms. The decay of linkage disequilibrium was slower in S. viminalis compared with S. schwerinii, and we speculate that this is due to different demographic histories as S. viminalis has been partly domesticated in Europe.