Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation

Major specialized natural products from the endangered plant Heptacodium miconioides, potential medicinal uses and insights into its longstanding unresolved systematic classification

A comprehensive phytochemical investigation of the flower buds and leaves/twigs of Heptacodium miconioides, a cultivated ornamental plant native to China and categorized as 'vulnerable', has led to the isolation of 45 structurally diverse compounds, which comprise 18 phenylpropanoids (1-4, 7-20), 11 pentacyclic triterpenoids (5, 6, 21-29), eight secoiridoid glycosides (30-37), three quinic acid derivatives (38-40), and a few miscellaneous components (41-45). Among them, (+)-α-intermedianol (1), (+)-holophyllol A (2), and (-)-pseudolarkaemin A (3) represent previously unreported enantiomeric lignans, while (+)-7'(R)-hydroxymatairesinol (4) is an undescribed naturally occurring lignan. Heptacoacids A (5) and B (6) are undescribed 24-nor-urs-28-oic acid derivatives. Their chemical structures were determined by 2D-NMR, supplemented by evidence from specific rotations and circular dichroism spectra. Given the uncertainty surrounding the systematic position of Heptacodium, integrative taxonomy (ITA), a method utilized to define contentious species, is applied. Chemotaxonomy, a vital aspect of ITA, becomes significant. By employing hierarchical clustering analysis (HCA) and syntenic pattern analysis methods, a taxonomic examination based on the major specialized natural products from the flower buds of H. miconioides and two other Caprifoliaceae plants (i.e., Lonicera japonica and Abelia × grandiflora) could offer enhanced understanding of the systematic placement of Heptacodium. Additionally, compounds 39 and 40 displayed remarkable inhibitory activities against ATP-citrate lyase (ACL), with IC50 values of 0.11 and 1.10 μM, respectively. In summary, the discovery of medical properties and refining systematic classification can establish a sturdy groundwork for conservation efforts aimed at mitigating species diversity loss while addressing human diseases.

Chromosome-level and haplotype-resolved genome assembly of Dracaena cambodiana (Asparagaceae)

Dracaena cambodiana Pierre ex Gagn. (Asparagaceae) is the source plant of Dragon's blood and has high ornamental values in gardening. Currently, this species is classified as the second-class state-protected species in the National Key Protected Wild Plants (NKPWP) of China. However, limited genomic data has hindered a more comprehensive scientific understanding of the processes involved in the production of Dragon's blood and the related conservation genomics research. In this study, we assembled a haplotype-resolved genome of D. cambodiana. The haploid genomes, haplotype A and haplotype B, are 1,015.22 Mb and 1,003.13 Mb in size, respectively. The completeness of haplotype A and haplotype B genomes was 98.60% and 98.20%, respectively, using the "embryophyta_10" dataset. Haplotype A and haplotype B genomes contained 27,361 and 27,066 protein-coding genes, respectively, with nearly all being functionally annotated. These findings provide new insights into the genomic characteristics of D. cambodiana and will offer additional genomic resources for studying the biosynthesis mechanism of Dragon's blood and the horticultural application of Dragon trees.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Transcriptome data analysis provides insights into the conservation of Michelia lacei, a plant species with extremely small populations distributed in Yunnan province, China

BACKGROUND

Michelia lacei W.W.Smith (Magnoliaceae), was classified as a Plant Species with Extremely Small Populations (PSESP) by the Yunnan Provincial Government in both action plans of 2012 and 2021. This evergreen tree is known for its high ornamental and scientific value, but it faces significant threats due to its extremely small population size and narrow geographical distribution. The study aims to understand the genetic structure, diversity, and demographic history of this species to inform its conservation strategies.

RESULTS

The analysis of transcriptome data from 64 individuals across seven populations of M. lacei identified three distinct genetic clusters and generated 104,616 single-nucleotide polymorphisms (SNPs). The KM ex-situ population, originating from Longling County, exhibited unique genetic features, suggesting limited gene flow. The genetic diversity was substantial, with significant differences between populations, particularly between the KM lineage and the OTHER lineage. Demographic history inferred from the data indicated population experienced three significant population declines during glaciations, followed by periods of recovery. We estimated the effective population size (Ne) of the KM and OTHER lineages 1,000 years ago were 85,851 and 416,622, respectively. Gene flow analysis suggested past gene flow between populations, but the KM ex-situ population showed no recent gene flow. A total of 805 outlier SNPs, associated with four environmental factors, suggest potential local adaptation and showcase the species' adaptive potential. Particularly, the BZ displayed 515 adaptive loci, highlighting its strong potential for adaptation within this group.

CONCLUSIONS

The comprehensive genomic analysis of M. lacei provides valuable insights into its genetic background and highlights the urgent need for conservation efforts. The study underscores the importance of ex-situ conservation methods, such as seed collection and vegetative propagation, to safeguard genetic diversity and promote population restoration. The preservation of populations like MC and BZ is crucial for maintaining the species' genetic diversity. In-situ conservation measures, including the establishment of in-situ conservation sites and community engagement, are essential to enhance protection awareness and ensure the long-term survival of this threatened plant species.

Regional environmental differences significantly affect the genetic structure and genetic differentiation of Carpinus tientaiensis Cheng, an endemic and extremely endangered species from China

Differences in topography and environment greatly affect the genetic structure and genetic differentiation of species, and endemic or endangered species with limited geographic ranges seem to be more sensitive to changes in climate and other environmental factors. The complex topography of eastern China is likely to affect genetic differentiation of plants there. Carpinus tientaiensis Cheng is a native and endangered plants from China, and exploring its genetic diversity has profound significance for protection and the collection of germplasm resources. Based on AFLP markers, this study found that C. tientaiensis has low genetic diversity, which mainly came from within populations, while Shangshantou and Tiantai Mountain populations have relatively high genetic diversity. The Nei genetic distance was closely related to geographical distance, and temperature and precipitation notablely affected the genetic variation and genetic differentiation of C. tientaiensis. Based on cpDNA, this study indicated that C. tientaiensis exhibits a moderate level of genetic diversity, and which mainly came from among populations, while Tiantai Mountain population have the highest genetic diversity. It demonstrated that there was genetic differentiation between populations, which can be divided into two independent geographical groups, but there was no significant phylogeographic structure between them. The MaxEnt model showed that climate change significantly affects its distribution, and the suitable distribution areas in Zhejiang were primarily divided into two regions, eastern Zhejiang and southern Zhejiang, and there was niche differentiation in its suitable distribution areas. Therefore, this study speculated that the climate and the terrain of mountains and hills in East China jointly shape the genetic structure of C. tientaiensis, which gived rise to an obvious north-south differentiation trend of these species, and the populations located in the hilly areas of eastern Zhejiang and the mountainous areas of southern Zhejiang formed two genetic branches respectively.

Transcriptome data reveals the conservation genetics of Cypripedium forrestii, a plant species with extremely small populations endemic to Yunnan, China

The Cypripedium forrestii is an orchid species with extremely small populations (PSESP) in Yunnan, China. C. forrestii is range-restricted and less-studied than many orchid species, and it is exposed to various threats to its survival. We investigated its potential habitats and collected 52 samples from eight locations, as well as two outgroup species for reference. We developed genetic markers (SNPs) for C. forrestii based on transcriptome sequencing (RNA-seq) data, and analyzed the genetic diversity, population structure, gene flow and demographic history of C. forrestii in detail. C. forrestii is a taxonomically independent species to protect. We found that the genetic diversity of C. forrestii was very low (1.7e-4) compared with other endangered species. We identified three genetic clusters, and several populations with distinct genetic backgrounds. Most genetic diversity was found within sampling sites (87.87%) and genetic clusters (91.39%). Gene flow has been greatly limited over the most recent generations, probably due to geographical distance, historical climate change and habitat fragmentation. We also detected a severe bottleneck event brought about by the recent population constraints. These factors, together with its reproductive characteristics, contribute to the population fragmentation and low genetic diversity of C. forrestii. Based on our findings, we suggest an integrative conservation strategy to protect and recover the genetic diversity of C. forrestii and a further comprehensive study of its ecological traits in the future.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

GAEP: a comprehensive genome assembly evaluating pipeline

With the rapid development of sequencing technologies, especially the maturity of third-generation sequencing technologies, there has been a significant increase in the number and quality of published genome assemblies. The emergence of these high-quality genomes has raised higher requirements for genome evaluation. Although numerous computational methods have been developed to evaluate assembly quality from various perspectives, the selective use of these evaluation methods can be arbitrary and inconvenient for fairly comparing the assembly quality. To address this issue, we have developed the Genome Assembly Evaluating Pipeline (GAEP), which provides a comprehensive assessment pipeline for evaluating genome quality from multiple perspectives, including continuity, completeness, and correctness. Additionally, GAEP includes new functions for detecting misassemblies and evaluating the assembly redundancy, which performs well in our testing. GAEP is publicly available at https://github.com/zy-optimistic/GAEP under the GPL3.0 License. With GAEP, users can quickly obtain accurate and reliable evaluation results, facilitating the comparison and selection of high-quality genome assemblies.

The genetic structure and demographic history revealed by whole-genome resequencing provide insights into conservation of critically endangered Artocarpus nanchuanensis

Introduction

Whole-genome resequencing technology covers almost all nucleotide variations in the genome, which makes it possible to carry out conservation genomics research on endangered species at the whole-genome level.

Methods

In this study, based on the whole-genome resequencing data of 101 critically endangered Artocarpus nanchuanensis individuals, we evaluated the genetic diversity and population structure, inferred the demographic history and genetic load, predicted the potential distributions in the past, present and future, and classified conservation units to propose targeted suggestions for the conservation of this critically endangered species.

Results

Whole-genome resequencing for A. nanchuanensis generated approximately 2 Tb of data. Based on abundant mutation sites (25,312,571 single nucleotide polymorphisms sites), we revealed that the average genetic diversity (nucleotide diversity, π) of different populations of A. nanchuanensis was relatively low compared with other trees that have been studied. And we also revealed that the NHZ and QJT populations harboured unique genetic backgrounds and were significantly separated from the other five populations. In addition, positive genetic selective signals, significantly enriched in biological processes related to terpene synthesis, were identified in the NHZ population. The analysis of demographic history of A. nanchuanensis revealed the existence of three genetic bottleneck events. Moreover, abundant genetic loads (48.56% protein-coding genes) were identified in Artocarpus nanchuanensis, especially in genes related to early development and immune function of plants. The predication analysis of suitable habitat areas indicated that the past suitable habitat areas shifted from the north to the south due to global temperature decline. However, in the future, the actual distribution area of A. nanchuanensis will still maintain high suitability.

Discussion

Based on total analyses, we divided the populations of A. nanchuanensis into four conservation units and proposed a number of practical management suggestions for each conservation unit. Overall, our study provides meaningful guidance for the protection of A. nanchuanensis and important insight into conservation genomics research.

A haplotype-resolved genome assembly of Rhododendron vialii based on PacBio HiFi reads and Hi-C data

Rhododendron vialii (subgen. Azaleastrum) is an evergreen shrub with high ornamental value. This species has been listed as a plant species with extremely small populations (PSESP) for urgent protection by China's Yunnan provincial government in 2021, due to anthropogenic habitat fragmentation. However, limited genomic resources hinder scientifically understanding of genetic threats that the species is currently facing. In this study, we assembled a high-quality haplotype-resolved genome of R. vialii based on PacBio HiFi long reads and Hi-C reads. The assembly contains two haploid genomes with sizes 532.73 Mb and 521.98 Mb, with contig N50 length of 35.67 Mb and 34.70 Mb, respectively. About 99.92% of the assembled sequences could be anchored to 26 pseudochromosomes, and 14 gapless assembled chromosomes were included in this assembly. Additionally, 60,926 protein-coding genes were identified, of which 93.82% were functionally annotated. This is the first reported genome of R. vialii, and hopefully it will lay the foundations for further research into the conservation genomics and horticultural domestication of this ornamentally important species.

Ploidy variation in Rhododendron subsection Maddenia and its implications for conservation

Polyploidy, which is common in plants, can confound taxon recognition and hence conservation assessments. In the taxonomically complex genus Rhododendron, 25 % of the over 1,300 taxa are considered under threat and 27 % Near Threatened or Data Deficient, with their taxonomy needing to be resolved urgently. Although ploidy levels of Rhododendron taxa range from diploid (2x) to dodecaploid (12x) according to previous reports, the extent of polyploidy across the genus has not been examined. We first summarized the taxonomic distribution of polyploids in the genus based on the literature. Then as a case study, we estimated ploidy levels of 47 taxa in subsection Maddenia (subgenus Rhododendron, section Rhododendron) using flow cytometry, together with verification of meiotic chromosome counts for representative taxa. The summary of reported ploidy in Rhododendron indicates that polyploidy is most common in subgenera Pentanthera and Rhododendron. In subsection Maddenia, all examined taxa are diploids except for the R. maddenii complex that shows a high ploidy variation (2-8x, 12x). We investigated ploidy level of 12 taxa in subsection Maddenia for the first time, and estimated genome sizes of two Rhododendron species. Knowledge of ploidy levels will inform phylogenetic analysis of unresolved species complexes. Overall, our study of subsection Maddenia provides a model for examining multiple issues including taxonomic complexity, ploidy variation and geographic distribution in relation to biodiversity conservation.

IRAPs in Combination with Highly Informative ISSRs Confer Effective Potentials for Genetic Diversity and Fidelity Assessment in Rhododendron

The species belonging to the Rhododendron genus are well-known for their colorful corolla. Molecular marker systems have the potential to elucidate genetic diversity as well as to assess genetic fidelity in rhododendrons. In the present study, the reverse transcription domains of long terminal repeat retrotransposons were cloned from rhododendrons and used to develop an inter-retrotransposon amplified polymorphism (IRAP) marker system. Subsequently, 198 polymorphic loci were generated from the IRAP and inter-simple sequence repeat (ISSR) markers, of which 119 were derived from the IRAP markers. It was shown that in rhododendrons, IRAP markers were superior to the ISSRs in some polymorphic parameters, such as the average number of polymorphic loci (14.88 versus 13.17). The combination of the IRAP and ISSR systems was more discriminative for detecting 46 rhododendron accessions than each of the systems on their own. Furthermore, IRAP markers demonstrated more efficiency in genetic fidelity detection of in-vitro-grown R. bailiense Y.P.Ma, C.Q.Zhang and D.F.Chamb, an endangered species recently recorded in Guizhzhou Province, China. The available evidence revealed the distinct properties of IRAP and ISSR markers in the rhododendron-associated applications, and highlighted the availability of highly informative ISSR and IRAP markers in the evaluation of genetic diversity and genetic fidelity of rhododendrons, which may facilitate preservation and genetic breeding of rhododendron plants.

Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly

BACKGROUND

The genus Rhododendron (Ericaceae), a species-rich and widely distributed genus of woody plants, is distinguished for the beautiful and diverse flowers. Rhododendron delavayi Franch. and Rhododendron irroratum Franch., are highly attractive species widely distributed in south-west China and abundant new varieties have been selected from their genetic resources.

METHODS

We constructed chromosome-scale genome assemblies for Rhododendron delavayi and Rhododendron irroratum. Phylogenetic and whole-genome duplication analyses were performed to elucidate the evolutionary history of Rhododendron. Further, different types of gene duplications were identified and their contributions to gene family expansion were investigated. Finally, comprehensive characterization and evolutionary analysis of R2R3-MYB and NBS-encoding genes were conducted to explore their evolutionary patterns.

RESULTS

The phylogenetic analysis classified Rhododendron species into two sister clades, 'rhododendrons' and 'azaleas'. Whole-genome duplication (WGD) analysis unveiled only one WGD event that occurred in Rhododendron after the ancestral γ triplication. Gene duplication and gene family expansion analyses suggested that the younger tandem and proximal duplications contributed greatly to the expansion of gene families involved in secondary metabolite biosynthesis and stress response. The candidate R2R3-MYB genes likely regulating anthocyanin biosynthesis and stress tolerance in Rhododendron will facilitate the breeding for ornamental use. NBS-encoding genes had undergone significant expansion and experienced species-specific gain and loss events in Rhododendron plants.

CONCLUSIONS

The reference genomes presented here will provide important genetic resources for molecular breeding and genetic improvement of plants in this economically important Rhododendron genus.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Creating novel ornamentals via new strategies in the era of genome editing

Ornamental breeding has traditionally focused on improving novelty, yield, quality, and resistance to biotic or abiotic stress. However, achieving these goals has often required laborious crossbreeding, while precise breeding techniques have been underutilized. Fortunately, recent advancements in plant genome sequencing and editing technology have opened up exciting new frontiers for revolutionizing ornamental breeding. In this review, we provide an overview of the current state of ornamental transgenic breeding and propose four promising breeding strategies that have already proven successful in crop breeding and could be adapted for ornamental breeding with the help of genome editing. These strategies include recombination manipulation, haploid inducer creation, clonal seed production, and reverse breeding. We also discuss in detail the research progress, application status, and feasibility of each of these tactics.

Does the critically endangered Rhododendron amesiae deserve top priority for conservation?

1.If a threatened plant has the problem of inaccurate species delimitation, its conservation programs that have previously been implemented might be debated.2.We made a comprehensive comparison of the critically endangered R. amesiae and its close relative R. concinnum, employing both morphological and population genomic data (ddRAD-seq).3.We suggest that the critically endangered R. amesiae can be merged into R. concinnum. Hence, the threatened status of R. amesiae is needed to be reevaluated.

Endophytic Pestalotiopsis species associated with Rhododendron in Cangshan Mountain, Yunnan Province, China

Rhododendron is an essential ornamental plant that is abundant in Yunnan province. In Cangshan Mountain, Yunnan, China, 61 species of Rhododendron have been reported. Endophytic fungi are internal plant tissue inhabitants that do not harm the host. It has emerged as an exciting research topic as they have the potential to provide numerous secondary metabolites. This study is focused on taxonomic novelties and new host records of endophytic fungi associated with Rhododendron plants collected from Cangshan Mountain in Yunnan Province, China. Pestalotiopsis fungi are associated with a vast array of plant species worldwide. In this study, fresh leaves of Rhododendron cyanocarpum, Rhododendron decorum, and Rhododendron delavayi were collected from Cangshan Mountain, Yunnan Province, China. Endophytic Pestalotiopsis fungi associated with Rhododendron were characterized based on phylogenetic analyses of combined ITS, TEF1-α, and TUB genes along with morphological characteristics. Six new species (Pestalotiopsis appendiculata, Pestalotiopsis cangshanensis, Pestalotiopsis daliensis, Pestalotiopsis fusoidea, Pestalotiopsis rosarioides, and Pestalotiopsis suae) and a new host record (Pestalotiopsis trachicarpicola) are described. Detailed descriptions and color photo plates of these species are provided. It is the first time that the endophytic fungi of Rhododendron plants in Cangshan Mountain have been studied.

Combined de novo transcriptomic and physiological analyses reveal RyALS3-mediated aluminum tolerance in Rhododendron yunnanense Franch

Rhododendron (Ericaceae) not only has ornamental value, but also has great medicinal and edible values. Many Rhododendron species are native to acid soils where aluminum (Al) toxicity limits plant productivity and species distribution. However, it remains unknown how Rhododendron adapts to acid soils. Here, we investigated the physiological and molecular mechanisms of Al tolerance in Rhododendron yunnanense Franch. We found that the shoots of R. yunnanense Franch did not accumulate Al after exposure of seedlings to 50 μM Al for 7 days but predominantly accumulated in roots, suggesting that root Al immobilization contributes to its high Al tolerance. Whole-genome de novo transcriptome analysis was carried out for R. yunnanense Franch root apex in response to 6 h of 50 μM Al stress. A total of 443,639 unigenes were identified, among which 1,354 and 3,413 were up- and down-regulated, respectively, by 6 h of 50 μM Al treatment. Both Gene Ontology (GO) enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that genes involved in "ribosome" and "cytoskeleton" are overrepresented. Additionally, we identified Al-tolerance homologous genes including a tonoplast-localized ABC transporter RyALS3; 1. Overexpression of RyALS3; 1 in tobacco plants confers transgenic plants higher Al tolerance. However, root Al content was not different between wild-type plants and transgenic plants, suggesting that RyALS3; 1 is responsible for Al compartmentalization within vacuoles. Taken together, integrative transcriptome, physiological, and molecular analyses revealed that high Al tolerance in R. yunnanense Franch is associated with ALS3; 1-mediated Al immobilization in roots.

Genomic Data Reveals Population Genetic and Demographic History of Magnolia fistulosa (Magnoliaceae), a Plant Species With Extremely Small Populations in Yunnan Province, China

Elucidating the genetic background of threatened species is fundamental to their management and conservation, and investigating the demographic history of these species is helpful in the determination of the threats facing them. The woody species of the genus Magnolia (Magnoliaceae) have high economic, scientific and ecological values. Although nearly half of all Magnolia species have been evaluated as threatened, to date there has been no population genetic study employing Next Generation Sequencing (NGS) technology in this genus. In the present study, we investigate the conservation genomics of Magnolia fistulosa, a threatened species endemic to the limestone area along the Sino-Vietnamese border, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. To increase the reliability of our statistical inferences, we employed two approaches, Stacks and ipyrad, for SNP calling. A total of 15,272 and 18,960, respectively, putatively neutral SNPs were generated by Stacks and ipyrad. Relatively high genetic diversity and large population divergence were detected in M. fistulosa. Although higher absolute values were calculated using the ipyrad data set, the two data sets showed the same trends in genetic diversity (π, H e), population differentiation (F ST) and inbreeding coefficients (F IS). A change in the effective population size of M. fistulosa within the last 1 Ma was detected, including a population decline about 0.5-0.8 Ma ago, a bottleneck event about 0.2-0.3 Ma ago, population fluctuations during the last glacial stage, and the recovery of effective population size after the last glacial maximum. Our findings not only lay the foundation for the future conservation of this species, but also provide new insights into the evolutionary history of the genus Magnolia in southeastern Yunnan, China.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Demographic history and identification of threats revealed by population genomic analysis provide insights into conservation for an endangered maple

Recent advancements in whole genome sequencing techniques capable of covering nearly all the nucleotide variations of a genome would make it possible to set up a conservation framework for threatened plants at the genomic level. Here we applied a whole genome resequencing approach to obtain genome-wide data from 105 individuals sampled from the 10 currently known extant populations of Acer yangbiense, an endangered species with fragmented habitats and restricted distribution in Yunnan, China. To inform meaningful conservation action, we investigated what factors might have contributed to the formation of its extremely small population sizes and what threats it currently suffers at a genomic level. Our results revealed that A. yangbiense has low genetic diversity and comprises different numbers of genetic groups based on neutral (seven) and selected loci (13), with frequent gene flow between populations. Repeated bottleneck events, particularly the most recent one occurring within ~10,000 years before present, which decreased its effective population size (N_e ) < 200, and severe habitat fragmentation resulting from anthropogenic activities as well as a biased gender ratio of mature individuals in its natural habitat, might have together contributed to the currently fragmented and endangered status of A. yangbiense. The species has suffered from inbreeding and deleterious mutation load, both of which varied among populations but had similar patterns; that is, populations with higher FROH (frequency of runs of homozygosity) always carried a larger number of deleterious mutations in the homozygous state than in populations with lower FROH. In addition, based on our genetic differentiation results, and the distribution patterns of homozygous deleterious mutations in individuals, we recommend certain conservation actions regarding the genetic rescue of A. yangbiense. Overall, our study provides meaningful insights into the conservation genetics and a framework for the further conservation for the endangered A. yangbiense.

The chromosome-scale genome assembly, annotation and evolution of Rhododendron henanense subsp. lingbaoense

Rhododendron henanense subsp. lingbaoense (hereafter referred to as R. henanense) is an endemic species naturally distributed in the Henan province, China, with high horticultural, ornamental and medicinal value. Herein, we report a de novo genome assembly for R. henanense using a combination of PacBio long read and Illumina short read sequencing technologies. In total, we assembled 634.07 Mb with a contig N50 of 2.5 Mb, representing ~96.93% of the estimated genome size. By applying Hi-C data, 13 pseudochromosomes of R. henanense genome were assembled, covering ~98.21% of the genome assembly. The genome was composed of ~65.76% repetitive sequences and 31,098 protein-coding genes, 88.77% of which could be functionally annotated. Rhododendron henanense displayed a high level of synteny with other Rhododendron species from the Hymenanthes subgenus. Our data also suggests that R. henanense genes related to stress responses have undergone expansion, which may underly the unique abiotic and biotic stress resistance of the species. This alpine Rhododendron chromosome-scale genome assembly provides fundamental molecular resources for germplasm conservation, breeding efforts, evolutionary studies, and elucidating the unique biological characteristics of R. henanense.