Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi)

Genomic structural variation is associated with hypoxia adaptation in high-altitude zokors

Zokors, an Asiatic group of subterranean rodents, originated in lowlands and colonized high-elevational zones following the uplift of the Qinghai-Tibet plateau about 3.6 million years ago. Zokors live at high elevation in subterranean burrows and experience hypobaric hypoxia, including both hypoxia (low oxygen concentration) and hypercapnia (elevated partial pressure of CO2). Here we report a genomic analysis of six zokor species (genus Eospalax) with different elevational ranges to identify structural variants (deletions and inversions) that may have contributed to high-elevation adaptation. Based on an assembly of a chromosome-level genome of the high-elevation species, Eospalax baileyi, we identified 18 large inversions that distinguished this species from congeners native to lower elevations. Small-scale structural variants in the introns of EGLN1, HIF1A, HSF1 and SFTPD of E. baileyi were associated with the upregulated expression of those genes. A rearrangement on chromosome 1 was associated with altered chromatin accessibility, leading to modified gene expression profiles of key genes involved in the physiological response to hypoxia. Multigene families that underwent copy-number expansions in E. baileyi were enriched for autophagy, HIF1 signalling and immune response. E. baileyi show a significantly larger lung mass than those of other Eospalax species. These findings highlight the key role of structural variants underlying hypoxia adaptation of high-elevation species in Eospalax.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

RegCloser: a robust regression approach to closing genome gaps

BACKGROUND

Closing gaps in draft genomes leads to more complete and continuous genome assemblies. The ubiquitous genomic repeats are challenges to the existing gap-closing methods, based on either the k-mer representation by the de Bruijn graph or the overlap-layout-consensus paradigm. Besides, chimeric reads will cause erroneous k-mers in the former and false overlaps of reads in the latter.

RESULTS

We propose a novel local assembly approach to gap closing, called RegCloser. It represents read coordinates and their overlaps respectively by parameters and observations in a linear regression model. The optimal overlap is searched only in the restricted range consistent with insert sizes. Under this linear regression framework, the local DNA assembly becomes a robust parameter estimation problem. We solved the problem by a customized robust regression procedure that resists the influence of false overlaps by optimizing a convex global Huber loss function. The global optimum is obtained by iteratively solving the sparse system of linear equations. On both simulated and real datasets, RegCloser outperformed other popular methods in accurately resolving the copy number of tandem repeats, and achieved superior completeness and contiguity. Applying RegCloser to a plateau zokor draft genome that had been improved by long reads further increased contig N50 to 3-fold long. We also tested the robust regression approach on layout generation of long reads.

CONCLUSIONS

RegCloser is a competitive gap-closing tool. The software is available at https://github.com/csh3/RegCloser . The robust regression approach has a prospect to be incorporated into the layout module of long read assemblers.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

Altitudinal variation in organ mass from three mountain systems: The case of mesquite lizard Sceloporus grammicus

High altitude environments provide a fertile ground for investigating the benefits of phenotypic adjustments at several levels of biological organization. Low oxygen partial pressure and low environmental temperature are the main limiting factors that promote phenotypic variation in different organs, such as the lung and heart. Although high-altitude environments act like natural laboratories, most morphological studies conducted to date lack replication. Here, we evaluated organ mass variation in nine populations of Sceloporus grammicus, throughout three altitudinal gradients (mountains) from the Trans-Mexican volcanic belt. A total of 84 individuals from three different altitudes at three different mountains were collected. Then, we used generalized linear models to analyze the pattern of variation in internal organs mass as a function of altitude and temperature. We observed a striking pattern of altitudinal variation in the size of cardiorespiratory organs: while heart mass increased with altitude and decreased with temperature, the lung showed a significant statistical interaction between mountain transect and temperature. Overall, our results support the hypothesis that cardiorespiratory organs should be bigger in populations occurring at higher altitudes. Moreover, the study of different mountain systems allowed us to observe some differences in one mountain in relation to the other two.

Molecular mechanisms of adaptive evolution in wild animals and plants

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

Gut microbiota facilitates adaptation of the plateau zokor (Myospalax baileyi) to the plateau living environment

Gut microbiota not only helps the hosts to perform many key physiological functions such as food digestion, energy harvesting and immune regulation, but also influences host ecology and facilitates adaptation of the host to extreme environments. Plateau zokors epitomize successful physiological adaptation to their living environment in the face of the harsh environment characterized by low temperature, low pressure and hypoxia in the Tibetan plateau region and high concentrations of CO₂ in their burrows. Therefore, here we used a metagenomic sequencing approach to explore how gut microbiota contributed to the adaptive evolution of the plateau zokor on the Qinghai-Tibet Plateau. Our metagenomic results show that the gut microbiota of plateau zokors on the Tibetan plateau is not only enriched in a large number of species related to energy metabolism and production of short-chain fatty acids (SCFAs), but also significantly enriched the KO terms that involve carbohydrate uptake pathways, which well address energy uptake in plateau zokors while also reducing inflammatory responses due to low pressure, hypoxia and high CO₂ concentrations. There was also a significant enrichment of tripeptidyl-peptidase II (TPPII) associated with antigen processing, apoptosis, DNA damage repair and cell division, which may facilitate the immune response and tissue damage repair in plateau zokors under extreme conditions. These results suggest that these gut microbiota and their metabolites together contribute to the physiological adaptation of plateau zokors, providing new insights into the contribution of the microbiome to the evolution of mammalian adaptation.

ChIP-seq profiling of H3K4me3 and H3K27me3 in an invasive insect, Bactrocera dorsalis

Introduction: While it has been suggested that histone modifications can facilitate animal responses to rapidly changing environments, few studies have profiled whole-genome histone modification patterns in invasive species, leaving the regulatory landscape of histone modifications in invasive species unclear. Methods: Here, we screen genome-wide patterns of two important histone modifications, trimethylated Histone H3 Lysine 4 (H3K4me3) and trimethylated Histone H3 Lysine 27 (H3K27me3), in adult thorax muscles of a notorious invasive pest, the Oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), using Chromatin Immunoprecipitation with high-throughput sequencing (ChIP-seq). Results: We identified promoters featured by the occupancy of H3K4me3, H3K27me3 or bivalent histone modifications that were respectively annotated with unique genes key to muscle development and structure maintenance. In addition, we found H3K27me3 occupied the entire body of genes, where the average enrichment was almost constant. Transcriptomic analysis indicated that H3K4me3 is associated with active gene transcription, and H3K27me3 is mostly associated with transcriptional repression. Importantly, we identified genes and putative motifs modified by distinct histone modification patterns that may possibly regulate flight activity. Discussion: These findings provide the first evidence of histone modification signature in B. dorsalis, and will be useful for future studies of epigenetic signature in other invasive insect species.

Vitamin D3 Metabolic Enzymes in Plateau Zokor (Myospalax baileyi) and Plateau Pika (Ochotona curzoniae): Expression and Response to Hypoxia

Vitamin D3 (D3) is produced endogenously from 7-dehydrocholesterol by irradiation and is an important secosteroid for the absorption of calcium and phosphate. Lithocholic acid (LCA) increases intestinal paracellular calcium absorption in a vitamin D receptor-dependent manner in vitamin D-deficient rats. The plateau zokor (Myospalax baileyi), a strictly subterranean species, and plateau pika are endemic to the Qinghai-Tibet Plateau. To verify whether the zokors were deficient in D3 and reveal the effects of hypoxia on D3 metabolism in the zokors and pikas, we measured the levels of 25(OH)D3, calcium, and LCA, and quantified the expression levels of D3 metabolism-related genes. The results showed an undetectable serum level of 25(OH)D3 and a significantly higher concentration of LCA in the serum of plateau zokor, but its calcium concentration was within the normal range compared with that of plateau pika and Sprague-Dawley rats. With increasing altitude, the serum 25(OH)D3 levels in plateau pika decreased significantly, and the mRNA and protein levels of CYP2R1 (in the liver) and CYP27B1 (in the kidney) in plateau pika decreased significantly. Our results indicate that plateau zokors were deficient in D3 and abundant in LCA, which might be a substitution of D3 in the zokor. Furthermore, hypoxia suppresses the metabolism of D3 by down-regulating the expression of CYP2R1 and CYP27B1 in plateau pika.

Whole genome assembly of the armored loricariid catfish Ancistrus triradiatus highlights herbivory signatures

The catfish Ancistrus triradiatus belongs to the species-rich family Loricariidae. Loricariids display remarkable traits such as herbivory, a benthic lifestyle, the absence of scales but the presence of dermal bony plates. They are exported as ornamental fish worldwide, with escaped fishes becoming a threat locally. Although genetic and phylogenetic studies are continuously increasing and developmental genetic investigations are underway, no genome assembly has been formally proposed for Loricariidae yet. We report a high-quality genome assembly of Ancistrus triradiatus using long and short reads, and a newly assembled transcriptome. The genome assembly is composed of 9530 scaffolds, including 85.6% of ray-finned fish BUSCOs, and 26,885 predicted protein-coding genes. The genomic GC content is higher than in other catfishes, reflecting the higher metabolism associated with herbivory. The examination of the SCPP gene family indicates that the genes presumably triggering scale loss when absent, are present in the scaleless A. triradiatus, questioning their explanatory role. The analysis of the opsin gene repertoire revealed that gene losses associated to the nocturnal lifestyle of catfishes were not entirely found in A. triradiatus, as the UV-sensitive opsin 5 is present. Finally, most gene family expansions were related to immunity except the gamma crystallin gene family which controls pupil shape and sub-aquatic vision. Thus, the genome of A. triradiatus reveals that fish herbivory may be related to the photic zone habitat, conditions metabolism, photoreception and visual functions. This genome is the first for the catfish suborder Loricarioidei and will serve as backbone for future genetic, developmental and conservation studies.

Expression Plasticity of Transposable Elements is Highly Associated with Organismal Re-adaptation to Ancestral Environments

Understanding the roles of phenotypic plasticity in adaptive evolution has gained recognition for decades. Studies involving multiple taxa have shown that gene expression plasticity serves as “long-term memory” to facilitate re-adaptations to ancestral environments. Nevertheless, the general pattern and the underlying genetic basis of expression plasticity remain unclear. The transposable elements (TEs) play crucial roles in gene expression regulation and are widely distributed within the genome. Given this, we re-analyzed the transcriptomic data of chicken (Gallus gallus) generated from a reciprocal transplant experiment to examine whether expression shifts of TEs are involved in the re-adaptation process. Similar to the protein-coding genes, the plastic changes of TEs overwhelmingly exceed the genetic changes in the re-adaptation process. Further, the associated TEs co-expressed with diverse genes to perform a regulatory activity. Thus, our study supports the general function of phenotypic plasticity in adaptive evolution, and suggests a regulatory functions of TEs in this process.

Genomic insights into zokors' phylogeny and speciation in China

SignificancePhylogenies are the basis of many ecological and evolutionary studies. However, zokor phylogeny and speciation patterns are heavily debated. This study disentangled the phylogeny and speciation of zokors genomically. Six species of the Eospalax were separated into high-altitude E. baileyi and E. smithi and the rest four low-altitude species by recent Qinghai-Tibet Plateau uplift 3.6 million y ago. E. rothschildi and E. smithi speciated south of the Qinling-Huaihe Line, where refuges were supplied during glaciation. Introgression and incomplete lineage sorting led to the complex phylogeny of zokors. Genomic islands were formed due to ancient polymorphisms and divergence hitchhiking. This study concluded that climatic, geological, and tectonic events shaped the phylogeny and speciation of zokors in China.