A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

First whole-genome sequence and assembly of the Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps), a critically endangered species, using Oxford Nanopore Technologies

The Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps) is currently considered one of the most endangered primates in the world and is classified as critically endangered [International union for conservation of nature (IUCN)]. It faces multiple threats, the most significant one being habitat loss due to deforestation in western Ecuador. Genomic tools are keys for the management of endangered species, but this requires a reference genome, which until now was unavailable for A. f. fusciceps. The present study reports the first whole-genome sequence and assembly of A. f. fusciceps generated using Oxford Nanopore long reads. DNA was extracted from a subadult male, and libraries were prepared for sequencing following the Ligation Sequencing Kit SQK-LSK112 workflow. Sequencing was performed using a MinION Mk1C sequencer. The sequencing reads were processed to generate a genome assembly. Two different assemblers were used to obtain draft genomes using raw reads, of which the Flye assembly was found to be superior. The final assembly has a total length of 2.63 Gb and contains 3,861 contigs, with an N50 of 7,560,531 bp. The assembly was analyzed for annotation completeness based on primate ortholog prediction using a high-resolution database, and was found to be 84.3% complete, with a low number of duplicated genes indicating a precise assembly. The annotation of the assembly predicted 31,417 protein-coding genes, comparable with other mammal assemblies. A reference genome for this critically endangered species will allow researchers to gain insight into the genetics of its populations and thus aid conservation and management efforts of this vulnerable species.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae)

KEY MESSAGE

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

A high-quality chromosome-level Eutrema salsugineum genome, an extremophile plant model

BACKGROUND

Eutrema salsugineum (2n = 14), a halophyte in the family Brassicaceae, is an attractive model to study abiotic stress tolerance in plants. Two versions of E. salsugineum genomes that previously reported were based on relatively short reads; thus, the repetitive regions were difficult to characterize.

RESULTS

We report the sequencing and assembly of the E. salsugineum (Shandong accession) genome using long-read sequencing and chromosome conformation capture data. We generated Oxford Nanopore long reads at high depth (> 60X) of genome coverage with additional short reads for error correction. The new assembly has a total size of 295.5 Mb with 52.8% repetitive sequences, and the karyotype of E. salsugineum is consistent with the ancestral translocation Proto-Calepineae Karyotype structure in both order and orientation. Compared with previous assemblies, this assembly has higher contiguity, especially in the centromere region. Based on this new assembly, we predicted 25,399 protein-coding genes and identified the positively selected genes associated with salt and drought stress responses.

CONCLUSION

The new genome assembly will provide a valuable resource for future genomic studies and facilitate comparative genomic analysis with other plants.

Genomic convergence underlying high-altitude adaptation in alpine plants

Evolutionary convergence is one of the most striking examples of adaptation driven by natural selection. However, genomic evidence for convergent adaptation to extreme environments remains scarce. Here, we assembled reference genomes of two alpine plants, Saussurea obvallata (Asteraceae) and Rheum alexandrae (Polygonaceae), with 37,938 and 61,463 annotated protein-coding genes. By integrating an additional five alpine genomes, we elucidated genomic convergence underlying high-altitude adaptation in alpine plants. Our results detected convergent contractions of disease-resistance genes in alpine genomes, which might be an energy-saving strategy for surviving in hostile environments with only a few pathogens present. We identified signatures of positive selection on a set of genes involved in reproduction and respiration (e.g., MMD1, NBS1, and HPR), and revealed signatures of molecular convergence on genes involved in self-incompatibility, cell wall modification, DNA repair and stress resistance, which may underlie adaptation to extreme cold, high ultraviolet radiation and hypoxia environments. Incorporating transcriptomic data, we further demonstrated that genes associated with cuticular wax and flavonoid biosynthetic pathways exhibit higher expression levels in leafy bracts, shedding light on the genetic mechanisms of the adaptive "greenhouse" morphology. Our integrative data provide novel insights into convergent evolution at a high-taxonomic level, aiding in a deep understanding of genetic adaptation to complex environments.

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae)

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

The chromosome-scale genome of Kobresia myosuroides sheds light on karyotype evolution and recent diversification of a dominant herb group on the Qinghai-Tibet Plateau

Kobresia species are common in meadows on the Qinghai-Tibet Plateau. They are important food resources for local livestock, and serve a critical foundation for ecosystem integration. Genetic resources of Kobresia species are scarce. Here, we generated a chromosome-level genome assembly for K. myosuroides (Cyperaceae), using PacBio long-reads, Illumina short-reads, and Hi-C technology. The final assembly had a total size of 399.9 Mb with a contig N50 value of 11.9 Mb. The Hi-C result supported a 29 pseudomolecules model which was in consistent with cytological results. A total of 185.5 Mb (44.89% of the genome) transposable elements were detected, and 26,748 protein-coding genes were predicted. Comparative analysis revealed that Kobresia plants have experienced recent diversification events during the late Miocene to Pliocene. Karyotypes analysis indicated that the fission and fusion of chromosomes have been a major driver of speciation, which complied with the lack of whole-genome duplication (WGD) in K. myosuroides genome. Generally, this high-quality reference genome provides insights into the evolution of alpine sedges, and may be helpful to endemic forage improvement and alpine ecosystem preservation.