Phased telomere-to-telomere reference genome and pangenome reveal an expansion of resistance genes during apple domestication.

The cultivated apple (Malus domestica Borkh.) is a cross-pollinated perennial fruit tree of great economic importance. Earlier versions of apple reference genomes were unphased, fragmented, and lacked comprehensive insights into the apple's highly heterozygous genome, which impeded advances in genetic studies and breeding programs. In this study, we assembled a haplotype-resolved telomere-to-telomere (T2T) reference genome for the diploid apple cultivar Golden Delicious. Subsequently, we constructed a pangenome based on 12 assemblies from wild and cultivated species to investigate the dynamic changes of functional genes. Our results revealed the gene gain and loss events during apple domestication. Compared with cultivated species, more gene families in wild species were significantly enriched in oxidative phosphorylation, pentose metabolic process, responses to salt, and abscisic acid biosynthesis process. Our analyses also demonstrated a higher prevalence of different types of resistance gene analogs (RGAs) in cultivars than their wild relatives, partially attributed to segmental and tandem duplication events in certain RGAs classes. Structural variations, mainly deletions and insertions, have affected the presence and absence of TIR-NB-ARC-LRR, NB-ARC-LRR, and CC-NB-ARC-LRR genes. Additionally, hybridization/introgression from wild species has also contributed to the expansion of resistance genes in domesticated apples. Our haplotype-resolved T2T genome and pangenome provide important resources for genetic studies of apples, emphasizing the need to study the evolutionary mechanisms of resistance genes in apple breeding.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape.

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Genetic diversity and gene expression diversity shape the adaptive pattern of the aquatic plant Batrachium bungei along an altitudinal gradient on the Qinghai-Tibet plateau.

It is an intriguing issue of evolutionary biology how genetic diversity and gene expression diversity shape the adaptive patterns. Comparative transcriptomic studies of wild populations in extreme environments provide critical insights into the relative contribution of genetic and expressive components. In this study, we analyzed the genetic diversity and gene expression diversity of 20 populations of the aquatic plant Batrachium bungei along elevations ranging from 2690 to 4896 m on the Qinghai-Tibet plateau (QTP). Based on single nucleotide polymorphisms (SNPs) and gene expression data from 100 individuals of B. bungei, we found that variation in genetic sequence was more sensitive to detect weak differentiation than gene expression. Using 292,613 high-quality SNPs, we documented a significant phylogeographical structure, a low within-population genetic diversity, and a high inter-population genetic differentiation in B. bungei populations. Analysis of relationship between geographic distance, genetic distance, and gene expression similarity showed that geographic isolation shaped gene flow patterns but not gene expression patterns. We observed a negative relationship between genetic diversity and gene expression diversity within and among B. bungei populations, and we demonstrated that as environmental conditions worsen with increasing altitude, genetic diversity played an increased role in maintaining the stability of populations, while the corresponding role of gene expression diversity decreased. These results suggested that genetic diversity and gene expression diversity might act as a complementary mechanism contributing to the long-term survival of B. bungei in extreme environments.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences.

BACKGROUND: Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates. RESULTS: The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP. CONCLUSION: The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Population transcriptomics uncover the relative roles of positive selection and differential expression in Batrachium bungei adaptation to the Qinghai-Tibetan plateau.

KEY MESSAGE: Positive selection genes are related to metabolism, while differentially expressed genes are related to photosynthesis, suggesting that genetic adaptation and expression regulation may play independent roles in different gene classes. Genome-wide investigation of the molecular mechanisms for high-altitude adaptation is an intriguing topic in evolutionary biology. The Qinghai-Tibet Plateau (QTP) with its extremely variable environments is an ideal site for studying high-altitude adaptation. Here, we used transcriptome data of 100 individuals from 20 populations collected from various altitudes on the QTP to investigate the adaptive mechanisms of the aquatic plant Batrachium bungei at both the genetic and transcriptional level. To explore genes and biological pathways that may contribute to QTP adaptation, we employed a two-step approach, in which we identified positively selected genes and differentially expressed genes using the landscape genomic and differential expression approaches. The positive selection analysis showed that genes involved in metabolic regulation played a crucial role in B. bungei adaptation to the extreme environments of the QTP, especially intense ultraviolet radiation. Altitude-based differential expression analysis suggested that B. bungei could increase the rate of energy dissipation or reduce the efficiency of light energy absorption by down regulating the expression of photosynthesis-related genes to adapt to the strong ultraviolet radiation. Weighted gene co-expression network analysis identified ribosomal genes as hubs of altitude adaptation in B. bungei. Only a small part of genes (about 10%) overlapped between positively selected genes and differentially expressed genes in B. bungei, suggesting that genetic adaptation and gene expression regulation might play relatively independent roles in different categories of functional genes. Taken together, this study enriches our understanding of the high-altitude adaptation mechanism of B. bungei on the QTP.

Comparative chloroplast genome and transcriptome analysis on the ancient genus Isoetes from China.

Isoetes is a famous living fossil that plays a significant role in the evolutionary studies of the plant kingdom. To explore the adaptive evolution of the ancient genus Isoetes from China, we focused on Isoetes yunguiensis (Q.F. Wang and W.C. Taylor), I. shangrilaensis (X. Li, Y.Q. Huang, X.K. Dai & X. Liu), I. taiwanensis (DeVol), I. sinensis (T.C. Palmer), I. hypsophila_GHC (Handel-Mazzetti), and I. hypsophila_HZS in this study. We sequenced, assembled, and annotated six individuals' chloroplast genomes and transcriptomes, and performed a series of analyses to investigate their chloroplast genome structures, RNA editing events, and adaptive evolution. The six chloroplast genomes of Isoetes exhibited a typical quadripartite structure with conserved genome sequence and structure. Comparative analyses of Isoetes species demonstrated that the gene organization, genome size, and GC contents of the chloroplast genome are highly conserved across the genus. Besides, our positive selection analyses suggested that one positively selected gene was statistically supported in Isoetes chloroplast genomes using the likelihood ratio test (LRT) based on branch-site models. Moreover, we detected positive selection signals using transcriptome data, suggesting that nuclear-encoded genes involved in the adaption of Isoetes species to the extreme environment of the Qinghai-Tibetan Plateau (QTP). In addition, we identified 291-579 RNA editing sites in the chloroplast genomes of six Isoetes based on transcriptome data, well above the average of angiosperms. RNA editing in protein-coding transcripts results from amino acid changes to increase their hydrophobicity and conservation in Isoetes, which may help proteins form functional three-dimensional structure. Overall, the results of this study provide comprehensive transcriptome and chloroplast genome resources and contribute to a better understanding of adaptive evolutionary and molecular biology in Isoetes.

Comparative Analysis of Chloroplast Genomes within Saxifraga (Saxifragaceae) Takes Insights into Their Genomic Evolution and Adaption to the High-Elevation Environment.

Saxifraga species are widely distributed in alpine and arctic regions in the Northern hemisphere. Highly morphological diversity within this genus brings great difficulties for species identification, and their typical highland living properties make it interesting how they adapt to the extreme environment. Here, we newly generated the chloroplast (cp) genomes of two Saxifraga species and compared them with another five Saxifraga cp genomes to understand the characteristics of cp genomes and their potential roles in highland adaptation. The genome size, structure, gene content, GC content, and codon usage pattern were found to be highly similar. Cp genomes ranged from 146,549 bp to 151,066 bp in length, most of which comprised 130 predicted genes. Yet, due to the expansion of IR regions, the second copy of rps19 in Saxifraga stolonifera was uniquely kept. Through sequence divergence analysis, we identified seven hypervariable regions and detected some signatures of regularity associated with genetic distance. We also identified 52 to 89 SSRs and some long repeats among seven Saxifraga species. Both ML and BI phylogenetic analyses confirmed that seven Saxifraga species formed a monophyletic clade in the Saxifragaceae family, and their intragenus relationship was also well supported. Additionally, the ndhI and ycf1 genes were considered under positive selection in species inhabiting relatively high altitudes. Given the conditions of intense light and low CO2 concentration in the highland, the products of these two genes might participate in the adaptation to the extreme environment.

Biased gene expression reveals the contribution of subgenome to altitude adaptation in allopolyploid Isoetes sinensis.

Allopolyploids are believed to inherit the genetic characteristics of its progenitors and exhibit stronger adaptability and vigor. The allotetraploid Isoetes sinensis was formed by the natural hybridization and polyploidization of two diploid progenitors, Isoetes taiwanensis and Isoetes yunguiensis, and was believed to have the potential to adapt to plateau environments. To explore the expression pattern of homoeologous genes and their contributions to altitude adaptation, we transplanted natural allotetraploid I. sinensis (TnTnYnYn) along the altitude gradient for a long-term, and harvested them in summer and winter, respectively. One year after transplanting, it still lived well, even in the extreme environment of the Qinghai-Tibet Plateau. Then, we performed high-throughput RNA sequencing to measure their gene expression level. A total of 7801 homoeologous genes were expressed, among which 5786 were identified as shared expression in different altitudes and seasons. We further found that altitude variations could change the subgenome bias trend of I. sinensis, but season could not. Moreover, the functions of uniquely expressed genes indicated that temperature might be an important restrictive factor during the adaptation process. Through the analysis of DEGs and uniquely expressed genes, we found that Y subgenome provided more contributions to high altitude adaptation than T subgenome. These adaptive traits to high altitude may be inherited from its plateau progenitor I. yunguiensis. Through weighted gene co-expression network analysis, pentatricopeptide repeats gene family and glycerophospholipid metabolism pathway were considered to play important roles in high-altitude adaptation. Totally, this study will enrich our understanding of allopolyploid in environmental adaptation.