Temporal Dynamics of Avian Populations during Pleistocene Revealed by Whole-Genome Sequences

Genetically based demographic reconstructions require careful consideration of generation time

Bakker et al. use Robinson et al.'s reconstruction of three species of vulture to illustrate how incorrect generation time estimates can yield inaccurate results, underscoring the importance of generation time specification for genetically based reconstructions, especially for comparisons and species of conservation concern.

River network rearrangements promote speciation in lowland Amazonian birds

Large Amazonian rivers impede dispersal for many species, but lowland river networks frequently rearrange, thereby altering the location and effectiveness of river barriers through time. These rearrangements may promote biotic diversification by facilitating episodic allopatry and secondary contact among populations. We sequenced genome-wide markers to evaluate the histories of divergence and introgression in six Amazonian avian species complexes. We first tested the assumption that rivers are barriers for these taxa and found that even relatively small rivers facilitate divergence. We then tested whether species diverged with gene flow and recovered reticulate histories for all species, including one potential case of hybrid speciation. Our results support the hypothesis that river rearrangements promote speciation and reveal that many rainforest taxa are micro-endemic, unrecognized, and thus threatened with imminent extinction. We propose that Amazonian hyper-diversity originates partly from fine-scale barrier displacement processes-including river dynamics-which allow small populations to differentiate and disperse into secondary contact.

Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater

Background

The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required.

Results

We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene.

Conclusions

The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.

Searching for genetic evidence of demographic decline in an arctic seabird: beware of overlapping generations

Genetic data are useful for detecting sudden population declines in species that are difficult to study in the field. Yet this indirect approach has its own drawbacks, including population structure, mutation patterns, and generation overlap. The ivory gull (Pagophila eburnea), a long-lived Arctic seabird, is currently suffering from rapid alteration of its primary habitat (i.e., sea ice), and dramatic climatic events affecting reproduction and recruitment. However, ivory gulls live in remote areas, and it is difficult to assess the population trend of the species across its distribution. Here we present complementary microsatellite- and SNP-based genetic analyses to test a recent bottleneck genetic signal in ivory gulls over a large portion of their distribution. With attention to the potential effects of population structure, mutation patterns, and sample size, we found no significant signatures of population decline worldwide. At a finer scale, we found a significant bottleneck signal at one location in Canada. These results were compared with predictions from simulations showing how generation time and generation overlap can delay and reduce the bottleneck microsatellite heterozygosity excess signal. The consistency of the results obtained with independent methods strongly indicates that the species shows no genetic evidence of an overall decline in population size. However, drawing conclusions related to the species' population trends will require a better understanding of the effect of age structure in long-lived species. In addition, estimates of the effective global population size of ivory gulls were surprisingly low (~1000 ind.), suggesting that the evolutionary potential of the species is not assured.

Varied diversification patterns and distinct demographic trajectories in Ethiopian montane forest bird (Aves: Passeriformes) populations separated by the Great Rift Valley

Taxon-specific characteristics and extrinsic climatic and geological forces may both shape population differentiation and speciation. In geographically and taxonomically focused investigations, differentiation may occur synchronously as species respond to the same external conditions. Conversely, when evolution is investigated in taxa with largely varying traits, population differentiation and speciation is complex and shaped by interactions of Earth's template and species-specific traits. As such, it is important to characterize evolutionary histories broadly across the tree of life, especially in geographic regions that are exceptionally diverse and under pressures from human activities such as in biodiversity hotspots. Here, using whole-genome sequencing data, we characterize genomic variation in populations of six Ethiopian Highlands forest bird species separated by a lowland biogeographic barrier, the Great Rift Valley (GRV). In all six species, populations on either side of the GRV exhibited significant but varying levels of genetic differentiation. Species' dispersal ability was negatively correlated with levels of population differentiation. Isolation with migration models indicated varied patterns of population differentiation and connectivity among populations of the focal species. We found that demographic histories-estimated for each individual-varied by both species and population but were consistent between individuals of the same species and sampling region. We found that genomic diversity varied by half an order of magnitude across species, and that this variation could largely be explained by the harmonic mean of effective population size over the past 200,000 years. Overall, we found that even in highly dispersive species like birds, the GRV acts as a substantial biogeographic barrier.

Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans)

Lotus (Nelumbo Adans.), a relict plant, is the testimony of long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remains unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus, respectively. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10-9 base/generation in lotus shoot stem cells. The past history of lotus revealed that the ancestors of lotus in eastern and southern Asia could be traced back ~20 million years ago (Mya) and experienced twice significant bottlenecks and population splits. We further identified the selected genes among three lotus groups in different habitats, suggesting that 453 genes between tropical and temperate group and 410 genes between two subgroups from Northeastern China and the Yangtze River - Yellow River Basin might play important roles in natural selection in lotus's adaptation and resilience. Our findings not only improve an understanding of the lotus evolutionary history and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for the relict plants.

The formation of avian montane diversity across barriers and along elevational gradients

Tropical mountains harbor exceptional concentrations of Earth's biodiversity. In topographically complex landscapes, montane species typically inhabit multiple mountainous regions, but are absent in intervening lowland environments. Here we report a comparative analysis of genome-wide DNA polymorphism data for population pairs from eighteen Indo-Pacific bird species from the Moluccan islands of Buru and Seram and from across the island of New Guinea. We test how barrier strength and relative elevational distribution predict population differentiation, rates of historical gene flow, and changes in effective population sizes through time. We find population differentiation to be consistently and positively correlated with barrier strength and a species' altitudinal floor. Additionally, we find that Pleistocene climate oscillations have had a dramatic influence on the demographics of all species but were most pronounced in regions of smaller geographic area. Surprisingly, even the most divergent taxon pairs at the highest elevations experience gene flow across barriers, implying that dispersal between montane regions is important for the formation of montane assemblages.

Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus)

Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.

Rates of Protein Evolution across the Marsupial Phylogeny: Heterogeneity and Link to Life-History Traits

Despite the importance of effective population size (Ne) in evolutionary and conservation biology, it remains unclear what factors have an impact on this quantity. The Nearly Neutral Theory of Molecular Evolution predicts a faster accumulation of deleterious mutations (and thus a higher dN/dS ratio) in populations with small Ne; thus, measuring dN/dS ratios in different groups/species can provide insight into their Ne. Here, we used an exome data set of 1,550 loci from 45 species of marsupials representing 18 of the 22 extant families, to estimate dN/dS ratios across the different branches and families of the marsupial phylogeny. We found a considerable heterogeneity in dN/dS ratios among families and species, which suggests significant differences in their Ne. Furthermore, our multivariate analyses of several life-history traits showed that dN/dS ratios (and thus Ne) are affected by body weight, body length, and weaning age.

Invasion genomics uncover contrasting scenarios of genetic diversity in a widespread marine invader

Invasion rates have increased in the past 100 y irrespective of international conventions. What characterizes a successful invasion event? And how does genetic diversity translate into invasion success? Employing a whole-genome perspective using one of the most successful marine invasive species world-wide as a model, we resolve temporal invasion dynamics during independent invasion events in Eurasia. We reveal complex regionally independent invasion histories including cases of recurrent translocations, time-limited translocations, and stepping-stone range expansions with severe bottlenecks within the same species. Irrespective of these different invasion dynamics, which lead to contrasting patterns of genetic diversity, all nonindigenous populations are similarly successful. This illustrates that genetic diversity, per se, is not necessarily the driving force behind invasion success. Other factors such as propagule pressure and repeated introductions are an important contribution to facilitate successful invasions. This calls into question the dominant paradigm of the genetic paradox of invasions, i.e., the successful establishment of nonindigenous populations with low levels of genetic diversity.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Genomic consequences of colonisation, migration and genetic drift in barn owl insular populations of the eastern Mediterranean

The study of insular populations was key in the development of evolutionary theory. The successful colonisation of an island depends on the geographic context, and specific characteristics of the organism and the island, but also on stochastic processes. As a result, apparently identical islands may harbour populations with contrasting histories. Here, we use whole genome sequences of 65 barn owls to investigate the patterns of inbreeding and genetic diversity of insular populations in the eastern Mediterranean Sea. We focus on Crete and Cyprus, islands with similar size, climate and distance to mainland, that provide natural replicates for a comparative analysis of the impacts of microevolutionary processes on isolated populations. We show that barn owl populations from each island have a separate origin, Crete being genetically more similar to other Greek islands and mainland Greece, and Cyprus more similar to the Levant. Further, our data show that their respective demographic histories following colonisation were also distinct. On the one hand, Crete harbours a small population and maintains very low levels of gene flow with neighbouring populations. This has resulted in low genetic diversity, strong genetic drift, increased relatedness in the population and remote inbreeding. Cyprus, on the other hand, appears to maintain enough gene flow with the mainland to avoid such an outcome. Our study provides a comparative population genomic analysis of the effects of neutral processes on a classical island‐mainland model system. It provides empirical evidence for the role of stochastic processes in determining the fate of diverging isolated populations.

Domestication obscures genomic estimates of population history

Domesticated species are valuable models to examine phenotypic evolution, and knowledge on domestication history is critical for understanding the trajectories of evolutionary changes. Sequentially Markov Coalescent models are often used to infer domestication history. However, domestication practices may obscure the signal left by population history, affecting demographic inference. Here we assembled the genomes of a recently domesticated species-the society finch-and its parent species-the white-rumped munia-to examine its domestication history. We applied genomic analyses to two society finch breeds and white-rumped munias to test whether domestication of the former resulted from inbreeding or hybridization. The society finch showed longer and more runs of homozygosity and lower genomic heterozygosity than the white-rumped munia, supporting an inbreeding origin in the former. Blocks of white-rumped munia and other ancestry in society finch genomes showed similar genetic distance between the two taxa, inconsistent with the hybridization origin hypothesis. We then applied two Sequentially Markov Coalescent models-psmc and smc++-to infer the demographic histories of both. Surprisingly, the two models did not reveal a recent population bottleneck, but instead the psmc model showed a specious, dramatic population increase in the society finch. Subsequently, we used simulated genomes based on an array of demographic scenarios to demonstrate that recent inbreeding, not hybridization, caused the distorted psmc population trajectory. Such analyses could have misled our understanding of the domestication process. Our findings stress caution when interpreting the histories of recently domesticated species inferred by psmc, arguing that these histories require multiple analyses to validate.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Life-history traits and habitat availability shape genomic diversity in birds: implications for conservation

More than 25% of species assessed by the International Union for Conservation of Nature (IUCN) are threatened with extinction. Understanding how environmental and biological processes have shaped genomic diversity may inform management practices. Using 68 extant avian species, we parsed the effects of habitat availability and life-history traits on genomic diversity over time to provide a baseline for conservation efforts. We used published whole-genome sequence data to estimate overall genomic diversity as indicated by historical long-term effective population sizes (Ne) and current genomic variability (H), then used environmental niche modelling to estimate Pleistocene habitat dynamics for each species. We found that Ne and H were positively correlated with habitat availability and related to key life-history traits (body mass and diet), suggesting the latter contribute to the overall genomic variation. We found that H decreased with increasing species extinction risk, suggesting that H may serve as a leading indicator of demographic trends related to formal IUCN conservation status in birds. Our analyses illustrate that genome-wide summary statistics estimated from sequence data reflect meaningful ecological attributes relevant to species conservation.

Displaced clines in an avian hybrid zone (Thamnophilidae: Rhegmatorhina) within an Amazonian interfluve

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi.

Genomic Evidence for Rare Hybridization and Large Demographic Changes in the Evolutionary Histories of Four North American Dove Species

Introductions and invasions provide opportunities for interaction and hybridization between colonists and closely related native species. We investigate this phenomenon using the mitochondrial DNA COI and 81,416 base-pairs of overlapping nuclear variation to examine the evolutionary histories and signatures of hybridization among introduced feral Rock Pigeon and Eurasian Collared-Dove and native White-winged and Mourning doves in southwestern North America. First, we report all four species to be highly divergent across loci (overall pair-wise species Φ_ST range = 0.17-0.70) and provide little evidence for gene flow at evolutionary timescales. Despite this, evidence from multiple population genetics analyses supports the presence of six putative contemporary late-stage hybrids among the 182 sampled individuals. These putative hybrids contain various ancestry combinations, but all involve the most populous species, the Mourning Dove. Next, we use a novel method to reconstruct demographic changes through time using partial genome sequence data. We identify recent, species-specific fluctuations in population size that are likely associated with changing environments since the Miocene and suggest that these fluctuations have influenced the genetic diversity of each dove species in ways that may impact their future persistence. Finally, we discuss the importance of using multiple marker types when attempting to infer complex evolutionary histories and propose important considerations when analyzing populations that were recently established or of domestic origins.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum

Pleistocene glacial cycles influenced the diversification of high‐latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger‐scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20–30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre‐LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high‐latitude species.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F_ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F_ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F_ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species.

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Genome of the Southern Giant Petrel Assembled Using Third-Generation DNA Sequencing and Linked Reads Reveals Evolutionary Traits of Southern Avian

The southern giant petrel Macronectes giganteus, a large seabird of the southern oceans, is one of only two members of the genus Macronectes and is the largest species in the order Procellariiformes. Although these two families account for the vast majority of the avian fauna inhabiting the Antarctic and sub-Antarctic regions, studies on the status of some populations and the associated genetic data are currently extremely limited. In this study, we assembled the genome of M. giganteus by integrating Pacific Biosciences single-molecule real-time sequencing and the Chromium system developed by 10x Genomics. The final M. giganteus genome assembly was 1.248 Gb in size with a scaffold N50 length of 27.4 Mb and a longest scaffold length of 120.4 Mb. The M. giganteus genome contains 14,993 predicted protein-coding genes and has 11.06% repeat sequences. Estimated historical effective population size analysis indicated that the southern giant petrel underwent a severe reduction in effective population size during a period coinciding with the early Pleistocene. The availability of this newly sequenced genome will facilitate more effective genetic monitoring of threatened species. Furthermore, the genome will provide a valuable resource for gene functional studies and further comparative genomic studies on the life history and ecological traits of specific avian species.

Temporal dynamics of teleost populations during the Pleistocene: a report from publicly available genome data

Background

Global climate oscillation, as a selection dynamic, is an ecologically important element resulting in global biodiversity. During the glacial geological periods, most organisms suffered detrimental selection pressures (such as food shortage and habitat loss) and went through population declines. However, during the mild interglacial periods, many species re-flourished. These temporal dynamics of effective population sizes (N_e) provide essential information for understanding and predicting evolutionary outcomes during historical and ongoing global climate changes.

Results

Using high-quality genome assemblies and corresponding sequencing data, we applied the Pairwise Sequentially Markovian Coalescent (PSMC) method to quantify N_e changes of twelve representative teleost species from approximately 10 million years ago (mya) to 10 thousand years ago (kya). These results revealed multiple rounds of population contraction and expansion in most of the examined teleost species during the Neogene and the Quaternary periods. We observed that 83% (10/12) of the examined teleosts had experienced a drastic decline in N_e before the last glacial period (LGP, 110–12 kya), slightly earlier than the reported pattern of N_e changes in 38 avian species. In comparison with the peaks, almost all of the examined teleosts maintained long-term lower N_e values during the last few million years. This is consistent with increasingly dramatic glaciation during this period.

Conclusion

In summary, these findings provide a more comprehensive understanding of the historical N_e changes in teleosts. Results presented here could lead to the development of appropriate strategies to protect species in light of ongoing global climate changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07816-7.

Phylogeny and Comparative Analysis for the Plastid Genomes of Five Tulipa (Liliaceae)

Species of Tulipa (Liliaceae) are of great horticultural importance and are distributed across Europe, North Africa, and Asia. The Tien Shan Mountain is one of the primary diversity centres of Tulipa, but the molecular studies of Tulipa species from this location are lacking. In our study, we assembled four Tulipa plastid genomes from the Tien Shan Mountains, T. altaica, T. iliensis, T. patens, and T. thianschanica, combined with the plastid genome of T. sylvestris to compare against other Liliaceae plastid genomes. We focussed on the species diversity and evolution of their plastid genomes. The five Tulipa plastid genomes proved highly similar in overall size (151,691-152,088 bp), structure, gene order, and content. With comparative analysis, we chose 7 mononucleotide SSRs from the Tulipa species that could be used in further population studies. Phylogenetic analyses based on 24 plastid genomes robustly supported the monophyly of Tulipa and the sister relationship between Tulipa and Amana, Erythronium. T. iliensis, T. thianschanica, and T. altaica were clustered together, and T. patens was clustered with T. sylvestris, with our results clearly demonstrating the relationships between these five Tulipa species. Our results provide a more comprehensive understanding of the phylogenomics and comparative genomics of Tulipa.

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Genetic load has potential in large populations but is realized in small inbred populations

Populations with higher genetic diversity and larger effective sizes have greater evolutionary capacity (i.e., adaptive potential) to respond to ecological stressors. We are interested in how the variation captured in protein-coding genes fluctuates relative to overall genomic diversity and whether smaller populations suffer greater costs due to their genetic load of deleterious mutations compared with larger populations. We analyzed individual whole-genome sequences (N = 74) from three different populations of Montezuma quail (Cyrtonyx montezumae), a small ground-dwelling bird that is sustainably harvested in some portions of its range but is of conservation concern elsewhere. Our historical demographic results indicate that Montezuma quail populations in the United States exhibit low levels of genomic diversity due in large part to long-term declines in effective population sizes over nearly a million years. The smaller and more isolated Texas population is significantly more inbred than the large Arizona and the intermediate-sized New Mexico populations we surveyed. The Texas gene pool has a significantly smaller proportion of strongly deleterious variants segregating in the population compared with the larger Arizona gene pool. Our results demonstrate that even in small populations, highly deleterious mutations are effectively purged and/or lost due to drift. However, we find that in small populations the realized genetic load is elevated because of inbreeding coupled with a higher frequency of slightly deleterious mutations that are manifested in homozygotes. Overall, our study illustrates how population genomics can be used to proactively assess both neutral and functional aspects of contemporary genetic diversity in a conservation framework while simultaneously considering deeper demographic histories.

Whole-genome resequencing provides insights into the population structure and domestication signatures of ducks in eastern China

Background

Duck is an ancient domesticated animal with high economic value, used for its meat, eggs, and feathers. However, the origin of indigenous Chinese ducks remains elusive. To address this question, we performed whole-genome resequencing to first explore the genetic relationship among variants of these domestic ducks with their potential wild ancestors in eastern China, as well as understand how the their genomes were shaped by different natural and artificial selective pressures.

Results

Here, we report the resequencing of 60 ducks from Chinese spot-billed ducks (Anas zonorhyncha), mallards (Anas platyrhnchos), Fenghua ducks, Shaoxing ducks, Shanma ducks and Cherry Valley Pekin ducks of eastern China (ten from each population) at an average effective sequencing depth of ~ 6× per individual. The results of population and demographic analysis revealed a deep phylogenetic split between wild (Chinese spot-billed ducks and mallards) and domestic ducks. By applying selective sweep analysis, we identified that several candidate genes, important pathways and GO categories associated with artificial selection were functionally related to cellular adhesion, type 2 diabetes, lipid metabolism, the cell cycle, liver cell proliferation, and muscle functioning in domestic ducks.

Conclusion

Genetic structure analysis showed a close genetic relationship of Chinese spot-billed ducks and mallards, which supported that Chinese spot-billed ducks contributed to the breeding of domestic ducks. During the long history of artificial selection, domestic ducks have developed a complex biological adaptation to captivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07710-2.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.

Genome-wide diversity in the California condor tracks its prehistoric abundance and decline

Due to their small population sizes, threatened and endangered species frequently suffer from a lack of genetic diversity, potentially leading to inbreeding depression and reduced adaptability.¹ During the latter half of the twentieth century, North America's largest soaring bird,² the California condor (Gymnogyps californianus; Critically Endangered³), briefly went extinct in the wild. Though condors once ranged throughout North America, by 1982 only 22 individuals remained. Following decades of captive breeding and release efforts, there are now >300 free-flying wild condors and ∼200 in captivity. The condor's recent near-extinction from lead poisoning, poaching, and loss of habitat is well documented,⁴ but much about its history remains obscure. To fill this gap and aid future management of the species, we produced a high-quality chromosome-length genome assembly for the California condor and analyzed its genome-wide diversity. For comparison, we also examined the genomes of two close relatives: the Andean condor (Vultur gryphus; Vulnerable³) and the turkey vulture (Cathartes aura; Least Concern³). The genomes of all three species show evidence of historic population declines. Interestingly, the California condor genome retains a high degree of variation, which our analyses reveal is a legacy of its historically high abundance. Correlations between genome-wide diversity and recombination rate further suggest a history of purifying selection against linked deleterious alleles, boding well for future restoration. We show how both long-term evolutionary forces and recent inbreeding have shaped the genome of the California condor, and provide crucial genomic resources to enable future research and conservation.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kākā (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kākā in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

Complete chloroplast genomes of Chinese wild-growing Vitis species: molecular structures and comparative and adaptive radiation analysis

As a basalmost family of Vitaceae, Chinese wild Vitis species offer key insights into the demographic history of grapes. In this study, we obtained 10 complete chloroplast (cp) genomes from Chinese wild-growing Vitis species based on our whole genome re-sequencing data. These chloroplast genomes ranged from 160,838 to 232,020 bp in size and exhibited typical quadripartite structures. Comparative analyses revealed that inverted repeat (IR) regions are especially abundant and contribute to cp genome arrangements. Phylogenetic analysis of the whole Vitis cp genomes supported three clearly partitioned main origins, in keeping with their geographic distributions, among which East Asian species from China were found to be sister species with Eurasian Vitis species but exhibited significant divergence from the North American group. Two well-supported subgroups were observed within the Chinese wild-growing Vitis species. Among these species, Vitis piasezkii and Vitis betulifolia were closely related species, exhibiting a support rate of 100%. The molecular clock-based divergence time suggested that the earliest split subspecies was Vitis pseudoreticulata, which further indicated that the origin and initial gene pool are located in southern China (the habitat of V. pseudoreticulata is located in the region). Coincidentally, the divergence time was during the Pleistocene period (2.6-0.1 Ma). Due to glacial/interglacial temperature fluctuations, cold-adapted subspecies, e.g., Vitis amurensis, could re-colonize new habitats. Our results may help to elucidate the adaptive radiation of Chinese wild Vitis species in different environments.

Population genomic, climatic and anthropogenic evidence suggest the role of human forces in endangerment of green peafowl (Pavo muticus)

Both anthropogenic impacts and historical climate change could contribute to population decline and species extinction, but their relative importance is still unclear. Emerging approaches based on genomic, climatic and anthropogenic data provide a promising analytical framework to address this question. This study applied such an integrative approach to examine potential drivers for the endangerment of the green peafowl (Pavo muticus). Several demographic reconstructions based on population genomes congruently retrieved a drastic population declination since the mid-Holocene. Furthermore, a comparison between historical and modern genomes suggested genetic diversity decrease during the last 50 years. However, climate-based ecological niche models predicted stationary general range during these periods and imply the little impact of climate change. Further analyses suggested that human disturbance intensities were negatively correlated with the green peafowl's effective population sizes and significantly associated with its survival status (extirpation or persistence). Archaeological and historical records corroborate the critical role of humans, leaving the footprint of low genomic diversity and high inbreeding in the survival populations. This study sheds light on the potential deep-time effects of human disturbance on species endangerment and offers a multi-evidential approach in examining underlying forces for population declines.

Vulnerability to Fishing and Life History Traits Correlate with the Load of Deleterious Mutations in Teleosts

Understanding why some species accumulate more deleterious substitutions than others is an important question relevant in evolutionary biology and conservation sciences. Previous studies conducted in terrestrial taxa suggest that life history traits correlate with the efficiency of purifying selection and accumulation of deleterious mutations. Using a large genome data set of 76 species of teleostean fishes, we show that species with life history traits associated with vulnerability to fishing have an increased rate of deleterious mutation accumulation (measured via dN/dS, i.e., nonsynonymous over synonymous substitution rate). Our results, focusing on a large clade of aquatic species, generalize previous patterns found so far in few clades of terrestrial vertebrates. These results also show that vulnerable species to fishing inherently accumulate more deleterious substitutions than nonthreatened ones, which illustrates the potential links among population genetics, ecology, and fishing policies to prevent species extinction.

Dynamic evolution of transposable elements, demographic history, and gene content of paleognathous birds

Palaeognathae includes ratite and tinamou species that are important for understanding early avian evolution. Here, we analyzed the whole-genome sequences of 15 paleognathous species to infer their demographic histories, which are presently unknown. We found that most species showed a reduction of population size since the beginning of the last glacial period, except for those species distributed in Australasia and in the far south of South America. Different degrees of contraction and expansion of transposable elements (TE) have shaped the paleognathous genome architecture, with a higher transposon removal rate in tinamous than in ratites. One repeat family, AviRTE, likely underwent horizontal transfer from tropical parasites to the ancestor of little and undulated tinamous about 30 million years ago. Our analysis of gene families identified rapid turnover of immune and reproduction-related genes but found no evidence of gene family changes underlying the convergent evolution of flightlessness among ratites. We also found that mitochondrial genes have experienced a faster evolutionary rate in tinamous than in ratites, with the former also showing more degenerated W chromosomes. This result can be explained by the Hill-Robertson interference affecting genetically linked W chromosomes and mitochondria. Overall, we reconstructed the evolutionary history of the Palaeognathae populations, genes, and TEs. Our findings of co-evolution between mitochondria and W chromosomes highlight the key difference in genome evolution between species with ZW sex chromosomes and those with XY sex chromosomes.

Revisiting the evolutionary history of domestic and wild ducks based on genomic analyses

Although domestic ducks have been important poultry species throughout human history, their origin remains enigmatic, with mallards and/or Chinese spot-billed ducks being proposed as the direct wild ancestor(s) of domestic ducks. Here, we analyzed 118 whole genomes from mallard, Chinese spot-billed, and domestic ducks to reconstruct their evolutionary history. We found pervasive introgression patterns among these duck populations. Furthermore, we showed that domestic ducks separated from mallard and Chinese spot-billed ducks nearly 38 thousand years ago (kya) and 54 kya, respectively, which is considerably outside the time period of presumed duck domestication. Thus, our results suggest that domestic ducks may have originated from another wild duck population that is currently undefined or unsampled, rather than from present-day mallard and/or Chinese spot-billed ducks, as previously thought. Overall, this study provides new insight into the complex evolution of ducks.

Comparative chloroplast genomes and phylogenetic analysis of Aquilegia

PREMISE

Aquilegia is an ideal taxon for studying the evolution of adaptive radiation. Current phylogenies of Aquilegia based on different molecular markers are inconsistent, and therefore a clear and accurate phylogeny remains uncertain. Analyzing the chloroplast genome, with its simple structure and low recombination rate, may help solve this problem.

METHODS

Next-generation sequencing data were generated or downloaded for Aquilegia species, enabling their chloroplast genomes to be assembled. The assemblies were used to estimate the genome characteristics and infer the phylogeny of Aquilegia.

RESULTS

In this study, chloroplast genome sequences were assembled for Aquilegia species distributed across Asia, North America, and Europe. Three of the genes analyzed (petG, rpl36, and atpB) were shown to be under positive selection and may be related to adaptation. The phylogenetic tree of Aquilegia showed that its member species formed two clades with high support, North American and European species, with the Asian species being paraphyletic; A. parviflora and A. amurensis clustered with the North American species, while the remaining Asian species were found in the European clade. In addition, A. oxysepala var. kansuensis should be considered as a separate species rather than a variety.

DISCUSSION

The complete chloroplast genomes of these Aquilegia species provide new insights into the reconstruction of the phylogeny of related species and contribute to the further study of this genus.

Persistent panmixia despite extreme habitat loss and population decline in the threatened tricolored blackbird (Agelaius tricolor)

Habitat loss and alteration has driven many species into decline, often to the point of requiring protection and intervention to avert extinction. Genomic data provide the opportunity to inform conservation and recovery efforts with details about vital evolutionary processes with a resolution far beyond that of traditional genetic approaches. The tricolored blackbird (Agelaius tricolor) has suffered severe losses during the previous century largely due to anthropogenic impacts on their habitat. Using a dataset composed of a whole genome paired with reduced representation libraries (RAD-Seq) from samples collected across the species' range, we find evidence for panmixia using multiple methods, including PCA (no geographic clustering), admixture analyses (ADMIXTURE and TESS conclude K = 1), and comparisons of genetic differentiation (average FST = 0.029). Demographic modeling approaches recovered an ancient decline that had a strong impact on genetic diversity but did not detect any effect from the known recent decline. We also did not detect any evidence for selection, and hence adaptive variation, at any site, either geographic or genomic. These results indicate that species continues to have high vagility across its range despite population decline and habitat loss and should be managed as a single unit.

Whole-genome assembly and resequencing reveal genomic imprint and key genes of rapid domestication in narrow-leafed lupin

Shifting from a livestock-based protein diet to a plant-based protein diet has been proposed as an essential requirement to maintain global food sustainability, which requires the increased production of protein-rich crops for direct human consumption. Meanwhile, the lack of sufficient genetic diversity in crop varieties is an increasing concern for sustainable food supplies. Countering this concern requires a clear understanding of the domestication process and dynamics. Narrow-leafed lupin (Lupinus angustifolius L.) has experienced rapid domestication and has become a new legume crop over the past century, with the potential to provide protein-rich seeds. Here, using long-read whole-genome sequencing, we assembled the third-generation reference genome for the narrow-leafed lupin cultivar Tanjil, comprising 20 chromosomes with a total genome size of 615.8 Mb and contig N50 = 5.65 Mb. We characterized the original mutation and putative biological pathway resulting in low seed alkaloid level that initiated the recent domestication of narrow-leafed lupin. We identified a 1133-bp insertion in the cis-regulatory region of a putative gene that may be associated with reduced pod shattering (lentus). A comparative analysis of genomic diversity in cultivars and wild types identified an apparent domestication bottleneck, as precisely predicted by the original model of the bottleneck effect on genetic variability in populations. Our results identify the key domestication genetic loci and provide direct genomic evidence for a domestication bottleneck, and open up the possibility of knowledge-driven de novo domestication of wild plants as an avenue to broaden crop plant diversity to enhance food security and sustainable low-carbon emission agriculture.

The demography of extinction in eastern North American birds

Species are being lost at an unprecedented rate during the Anthropocene. Progress has been made in clarifying how species traits influence their propensity to go extinct, but the role historical demography plays in species loss or persistence is unclear. In eastern North America, five charismatic landbirds went extinct last century, and the causes of their extinctions have been heavily debated. Although these extinctions are most often attributed to post-colonial human activity, other factors such as declining ancestral populations prior to European colonization could have made these species particularly susceptible. We used population genomic data from these extinct birds and compared them with those from four codistributed extant species. We found extinct species harboured lower genetic diversity and effective population sizes than extant species, but both extinct and non-extinct birds had similar demographic histories of population expansion. These demographic patterns are consistent with population size changes associated with glacial-interglacial cycles. The lack of support for overall population declines during the Pleistocene corroborates the view that, although species that went extinct may have been vulnerable due to low diversity or small population size, their disappearance was driven by human activities in the Anthropocene.

De novo assembly of a chromosome-scale reference genome for the northern flicker Colaptes auratus

The northern flicker, Colaptes auratus, is a widely distributed North American woodpecker and a long-standing focal species for the study of ecology, behavior, phenotypic differentiation, and hybridization. We present here a highly contiguous de novo genome assembly of C. auratus, the first such assembly for the species and the first published chromosome-level assembly for woodpeckers (Picidae). The assembly was generated using a combination of short-read Chromium 10× and long-read PacBio sequencing, and further scaffolded with chromatin conformation capture (Hi-C) reads. The resulting genome assembly is 1.378 Gb in size, with a scaffold N50 of 11  and a scaffold L50 of 43.948 Mb. This assembly contains 87.4-91.7% of genes present across four sets of universal single-copy orthologs found in tetrapods and birds. We annotated the assembly both for genes and repetitive content, identifying 18,745 genes and a prevalence of ∼28.0% repetitive elements. Lastly, we used fourfold degenerate sites from neutrally evolving genes to estimate a mutation rate for C. auratus, which we estimated to be 4.007 × 10-9 substitutions/site/year, about 1.5× times faster than an earlier mutation rate estimate of the family. The highly contiguous assembly and annotations we report will serve as a resource for future studies on the genomics of C. auratus and comparative evolution of woodpeckers.

Fine-scale genetic structure in the critically endangered red-fronted macaw in the absence of geographic and ecological barriers

Behavioural and socio-cultural traits are recognized in the restriction of gene flow in species with high cognitive capacity and complex societies. This isolation by social barriers has been generally overlooked in threatened species by assuming disrupted gene flow due to population fragmentation and decline. We examine the genetic structure and ecology of the global population of the Critically Endangered red-fronted macaw (Ara rubrogenys), an endemic species to the inter-Andean valleys of Bolivia. We found a fine-scale genetic structuring in four genetic clusters. Genetic diversity was higher in wild compared to captive-bred macaws, but similar to that of captive wild-caught macaws. We found no clear evidence of severe genetic erosion in the population in recent decades, but it was patent in historic times, overlapping with drastic human habitat transformation and macaw persecution over millennia. We found no evidence of geographical and ecological barriers, owing to the high dispersal ability, nesting and foraging habits between genetic clusters. The lack of genetic intermixing despite long-distance foraging and seasonal movements suggests recruitment in natal colonies and other social factors reinforcing philopatry-related genetic structure. Conservation efforts should be specifically focussed on major threats in each genetic cluster as independent conservation units, and also considered in ex-situ management.

Contrasting Trends of Population Size Change for Two Eurasian Owlet Species—Athene brama and Glaucidium radiatum From South Asia Over the Late Quaternary

Climatic oscillations over the Quaternary have had a lasting impact on species’ distribution, evolutionary history, and genetic composition. Many species show dramatic population size changes coinciding with the last glacial period. However, the extent and direction of change vary across biogeographic regions, species-habitat associations, and species traits. Here we use genomic data to assess population size changes over the late Quaternary using the Pairwise Sequential Markovian Coalescent (PSMC) approach in two Eurasian Owlet species—the Spotted Owlet, Athene brama, and the Jungle Owlet, Glaucidium radiatum. While Spotted Owlets are typically associated with open habitats, Jungle Owlets are found in deciduous forests and scrublands. We find that the effective population size for the Spotted Owlet increased after the Interglacial period till the Last Glacial Maxima and subsequently declined toward the Mid-Holocene. On the other hand, effective population size estimates for the Jungle Owlet increased gradually throughout this period. These observations are in line with climatic niche model-based predictions for range size change for both species from a previous study and suggest that habitat associations at the local scale are important in determining responses to past climatic and vegetational changes. The Spotted Owlet result also aligns well with the expectation of open habitat expansion during the arid Glacial Maxima, whereas for the Jungle Owlet the contrasting expectation does not hold. Therefore, assessing the impacts of glacial history on population trajectories of multiple species with different habitat associations is necessary to understand the impacts of past climate on South Asian taxa.

Are population isolations and declines a threat to island endemic water striders? A lesson from demographic and niche modelling of Metrocoris esakii (Hemiptera: Gerridae)

Genetic stochasticity and bottlenecking in the course of Pleistocene glaciations have been identified as threatening the survival of local endemics. However, the mechanisms by which local endemic species balance the influences of these two events remain poorly understood. Here, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) data set, mined mitochondrial sequences and constructed ecological niche models for the island endemic water strider Metrocoris esakii (Hemiptera: Gerridae). We found that M. esakii comprised three divergent lineages (i.e., north, central and south) isolated by geographical barriers and generally experienced population declines with the constriction of suitable areas during the Last Glacial Maximum (LGM). Further demographic model testing and stairway plots revealed a history of recent gene flow among the neighbouring lineages and rapid recovery at the end of the LGM, indicating that M. esakii at least had the potential for an adaptive response to population fragmentation and bottlenecking. The northern lineage did not show genetic bottlenecking during the LGM, which was probably due to its large effective population size (N_e ) from migration, which improved its adaptive potential. Relative to the ddRAD-seq data set, the demographic results based on mitochondrial sequences were less conclusive, showing weak differentiation and oversimplified demographic trajectories for the three genetic lineages. Overall, this study provides some degree of optimism for the survival of island endemic water striders from a demographic perspective, but further evaluation of their extinction risk under the impacts of human activities is required.

Comparative phylogeographic and demographic analyses reveal a congruent pattern of sister relationships between bird populations of the northern and south-central Atlantic Forest

The Pernambuco Center of Endemism (PCE) is the northernmost strip of the Atlantic Forest (AF). Biogeographic affinities among avifaunas in the PCE, the southern-central Atlantic Forest (SCAF), and Amazonia (AM) have not been studied comprehensively, and current patterns of genetic diversity in the PCE remain unclear. The interplay between species' ecological attributes and historical processes, such as Pleistocene climate fluctuations or the appearance of rivers, may have affected population genetic structures in the PCE. Moreover, the role of past connections between the PCE and AM and the elevational distribution of species in assembling the PCE avifauna remain untested. Here, we investigated the biogeographic history of seven taxa endemic to the PCE within a comparative phylogeographic framework based on a mean of 3,618 independent single nucleotide polymorphisms (SNPs) extracted from flanking regions of ultraconserved elements (UCEs) and one mitochondrial gene. We found that PCE populations were more closely related to SCAF populations than they were to those in AM, regardless of their elevational range, with divergence times placed during the Mid-Pleistocene. These splits were consistent with a pattern of allopatric divergence with gene flow until the upper Pleistocene and no signal of rapid changes in population sizes. Our results support the existence of a Pleistocene refugium driving current genetic diversity in the PCE, thereby rejecting the role of the São Francisco River as a primary barrier for population divergence. Additionally, we found that connections with Amazonia also played a significant role in assembling the PCE avifauna through subsequent migration events.

De Novo Assembly of the Northern Cardinal (Cardinalis cardinalis) Genome Reveals Candidate Regulatory Regions for Sexually Dichromatic Red Plumage Coloration

Northern cardinals (Cardinalis cardinalis) are common, mid-sized passerines widely distributed in North America. As an iconic species with strong sexual dichromatism, it has been the focus of extensive ecological and evolutionary research, yet genomic studies investigating the evolution of genotype-phenotype association of plumage coloration and dichromatism are lacking. Here we present a new, highly-contiguous assembly for C. cardinalis We generated a 1.1 Gb assembly comprised of 4,762 scaffolds, with a scaffold N50 of 3.6 Mb, a contig N50 of 114.4 kb and a longest scaffold of 19.7 Mb. We identified 93.5% complete and single-copy orthologs from an Aves dataset using BUSCO, demonstrating high completeness of the genome assembly. We annotated the genomic region comprising the CYP2J19 gene, which plays a pivotal role in the red coloration in birds. Comparative analyses demonstrated non-exonic regions unique to the CYP2J19 gene in passerines and a long insertion upstream of the gene in C. cardinalis Transcription factor binding motifs discovered in the unique insertion region in C. cardinalis suggest potential androgen-regulated mechanisms underlying sexual dichromatism. Pairwise Sequential Markovian Coalescent (PSMC) analysis of the genome reveals fluctuations in historic effective population size between 100,000-250,000 in the last 2 millions years, with declines concordant with the beginning of the Pleistocene epoch and Last Glacial Period. This draft genome of C. cardinalis provides an important resource for future studies of ecological, evolutionary, and functional genomics in cardinals and other birds.

Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi

Background

The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai–Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history.

Results

Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche.

Conclusions

The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.

Genome-wide analyses reveal drivers of penguin diversification

Penguins are the only extant family of flightless diving birds. They currently comprise at least 18 species, distributed from polar to tropical environments in the Southern Hemisphere. The history of their diversification and adaptation to these diverse environments remains controversial. We used 22 new genomes from 18 penguin species to reconstruct the order, timing, and location of their diversification, to track changes in their thermal niches through time, and to test for associated adaptation across the genome. Our results indicate that the penguin crown-group originated during the Miocene in New Zealand and Australia, not in Antarctica as previously thought, and that Aptenodytes is the sister group to all other extant penguin species. We show that lineage diversification in penguins was largely driven by changing climatic conditions and by the opening of the Drake Passage and associated intensification of the Antarctic Circumpolar Current (ACC). Penguin species have introgressed throughout much of their evolutionary history, following the direction of the ACC, which might have promoted dispersal and admixture. Changes in thermal niches were accompanied by adaptations in genes that govern thermoregulation and oxygen metabolism. Estimates of ancestral effective population sizes (N _e ) confirm that penguins are sensitive to climate shifts, as represented by three different demographic trajectories in deeper time, the most common (in 11 of 18 penguin species) being an increased N _e between 40 and 70 kya, followed by a precipitous decline during the Last Glacial Maximum. The latter effect is most likely a consequence of the overall decline in marine productivity following the last glaciation.

Whole-genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole-genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long-term population dynamics of Branchiostoma. The inferred ancestral N e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10-9) or third glaciation (mutation rate = 9 × 10-9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long-term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.

Paleoclimatic evolution as the main driver of current genomic diversity in the widespread and polymorphic Neotropical songbird Arremon taciturnus

Several factors have been proposed as drivers of species diversification in the Neotropics, including environmental heterogeneity, the development of drainage systems and historical changes in forest distribution due to climatic oscillations. Here, we investigate which drivers contributed to the evolutionary history and current patterns of diversity of a polymorphic songbird (Arremon taciturnus) that is widely distributed in Amazonian and Atlantic forests as well as in Cerrado gallery and seasonally-dry forests. We use genomic, phenotypic and habitat heterogeneity data coupled with climatic niche modelling. Results suggest the evolutionary history of the species is mainly related to paleoclimatic changes, although changes in the strength of the Amazon river as a barrier to dispersal, current habitat heterogeneity and geographic distance were also relevant. We propose an ancestral distribution in the Guyana Shield, and recent colonization of areas south of the Amazon river at ~380 to 166 kya, and expansion of the distribution to southern Amazonia, Cerrado and the Atlantic Forest. Since then, populations south of the Amazon River have been subjected to cycles of isolation and possibly secondary contact due to climatic changes that affected habitat heterogeneity and population connectivity. Most Amazonian rivers are not associated with long lasting isolation of populations, but some might act as secondary barriers, susceptible to crossing under specific climatic conditions. Morphological variation, while stable in some parts of the distribution, is not a reliable indicator of genetic structure or phylogenetic relationships.