Domestication obscures genomic estimates of population history

Origin of the Laurentian Great Lakes fish fauna through upward adaptive radiation cascade prior to the Last Glacial Maximum

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex ("ciscoes") radiated rapidly prior to the Last Glacial Maximum (80-90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr (Salvelinus namaycush), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation.

Evolution of Chromosomal Inversions across an Avian Radiation

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

The genomes of Hercules beetles reveal putative adaptive loci and distinct demographic histories in pristine North American forests

Beetles, despite their remarkable biodiversity and a long history of research, remain lacking in reference genomes annotated with structural variations in loci of adaptive significance. We sequenced and assembled high-quality chromosome-level genomes of four Hercules beetles which exhibit divergence in male horn size and shape and body colouration. The four Hercules beetle genomes were assembled to 11 pseudo-chromosomes, where the three genomes assembled using Nanopore data (Dynastes grantii, D. hyllus and D. tityus) were mapped to the genome assembled using PacBio + Hi-C data (D. maya). We demonstrated a striking similarity in genome structure among the four species. This conservative genome structure may be attributed to our use of the D. maya assembly as the reference; however, it is worth noting that such a conservative genome structure is a recurring phenomenon among scarab beetles. We further identified homologues of nine and three candidate-gene families that may be associated with the evolution of horn structure and body colouration respectively. Structural variations in Scr and Ebony2 were detected and discussed for their putative impacts on generating morphological diversity in beetles. We also reconstructed the demographic histories of the four Hercules beetles using heterozygosity information from the diploid genomes. We found that the demographic histories of the beetles closely recapitulated historical changes in suitable forest habitats driven by climate shifts.

The Genetic Selection of HSPD1 and HSPE1 Reduce Inflammation of Liver and Spleen While Restraining the Growth and Development of Skeletal Muscle in Wuzhishan Pigs

Wuzhishan (WZS) pigs, which are minipigs native to Hainan Province in China, are characterized by strong resistance to extreme hot temperatures and humidity. The relationship between their immune response and growth still needs to be clarified. In this study, we used whole genome sequencing (WGS) to detect variations within 37 WZS pigs, 32 Large White (LW) pigs, and 22 Xiangxi black (XXB) pigs, and ~2.49 GB of SNPs were obtained. These data were combined with those of two other pig breeds, and it was found that most of the genes detected (354) were located within the distinct genetic regions between WZS pigs and LW pigs. The network that was constructed using these genes represented a center including 12 hub genes, five of which had structural variations (SVs) within their regulatory regions. Furthermore, RNA-seq and RT-qPCR data for 12 genes were primarily consistent in liver, spleen, and LDM tissues. Notably, the expression of HSPs (HSPD1 and HSPE1) was higher while that of most genes involved in the JAK3-STAT pathway were lower in liver tissue of WZS pigs, compared with LW pigs. This likely not only reduced inflammation-related immune response but also impaired their growth. Our findings demonstrated the role of HSPs in the connection between inflammation and growth rate, while also providing the fundamental genetic selection of the adaptability of WZS pigs.

Chromosomal-level reference genome of a wild North American mallard (Anas platyrhynchos)

The mallard (Anas platyrhynchos) is one of the most common, economically, and socially important birds around the world. Mallards were not only an important food source for early humans but eventually becoming intimately linked with people as they were domesticated over the last 2,000 years. To date, mallard genomes are largely reconstructed from samples of domestic or unknown genetic heritage. Here, we report the first high-quality genome assembly and annotation of a genetically vetted wild mallard from North America (NAwild_v1.0). The genome was assembled using a combination of shotgun libraries, proximity ligation Chicago, and Dovetail Hi-C libraries. The final assembly is ∼1.04 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome-level scaffolds, and with a N50/L50 of 79.1 Mb/4 scaffolds. We used a combination of gene prediction and similarity approaches to annotate a total of 23,584 functional genes, of which 19,242 were associated to GO terms. The genome assembly and the set of annotated genes yielded a 95.4% completeness score when compared with the BUSCO aves_odb10 dataset. Next, we aligned 3 previously published mallard genomes to ours, and demonstrate how runs of homozygosity and nucleotide diversity are substantially higher and lower, respectively, to ours and how these artificially changed genomes resulted in profoundly different and biased demographic histories. Our wild mallard assembly not only provides a valuable resource to shed light onto genome evolution, speciation, and other adaptive processes, but also helping with identifying functional genes that have been significantly altered during the domestication process.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

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

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