Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

An alternative 3' splice site of PeuHKT1;3 improves the response to salt stress through enhancing affinity to K+ in Populus

Alternative splicing (AS) serves as a crucial post-transcriptional regulator in plants that contributes to the resistance to salt stress. However, the underlying mechanism is largely unknown. In this research, we identified an important AS transcript in Populus euphratica, PeuHKT1:3a, generated by alternative 3' splice site splicing mode that resulted in the removal of 252 bases at the 5' end of the first exon in PeuHKT1:3. Protein sequence comparison showed that the site of AS occurred in PeuHKT1:3 is located at a crucial Ser residue within the first pore-loop domain, which leads to inefficient K+ transport in HKT I-type transporters. Expressing PeuHKT1;3a in an axt3 mutant yeast strain can effectively compensate for the lack of intracellular K+, whereas the expression of PeuHKT1;3 cannot yield the effect. Furthermore, in transgenic Arabidopsis and poplar plants, it was observed that lines expressing PeuHKT1;3a exhibited greater salt tolerance compared to those expressing the PeuHKT1;3 strain. Analysis of ion content and flux demonstrated that the transgenic PeuHKT1;3a line exhibited significantly higher K+ content compared to the PeuHKT1;3 line, while there was no significant difference in Na+ content. In conclusion, our findings revealed that AS can give rise to novel variants of HKT I-type proteins in P. euphratica with modified K+ selectivity to keep a higher K+/Na+ ratio to enhanced salt tolerance.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Comparative genomic analysis of the Growth-Regulating Factors-Interacting Factors (GIFs) in six Salicaceae species and functional analysis of PeGIF3 reveals their regulatory role in Populus heteromorphic leaves

BACKGROUND

The growth-regulating factor-interacting factor (GIF) gene family plays a vital role in regulating plant growth and development, particularly in controlling leaf, seed, and root meristem homeostasis. However, the regulatory mechanism of heteromorphic leaves by GIF genes in Populus euphratica as an important adaptative trait of heteromorphic leaves in response to desert environment remains unknown.

RESULTS

This study aimed to identify and characterize the GIF genes in P. euphratica and other five Salicaceae species to investigate their role in regulating heteromorphic leaf development. A total of 27 GIF genes were identified and characterized across six Salicaceae species (P. euphratica, Populus pruinose, Populus deltoides, Populus trichocarpa, Salix sinopurpurea, and Salix suchowensis) at the genome-wide level. Comparative genomic analysis among these species suggested that the expansion of GIFs may be derived from the specific Salicaceae whole-genome duplication event after their divergence from Arabidopsis thaliana. Furthermore, the expression data of PeGIFs in heteromorphic leaves, combined with functional information on GIF genes in Arabidopsis, indicated the role of PeGIFs in regulating the leaf development of P. euphratica, especially PeGIFs containing several cis-acting elements associated with plant growth and development. By heterologous expression of the PeGIF3 gene in wild-type plants (Col-0) and atgif1 mutant of A. thaliana, a significant difference in leaf expansion along the medial-lateral axis, and an increased number of leaf cells, were observed between the overexpressed plants and the wild type.

CONCLUSION

PeGIF3 enhances leaf cell proliferation, thereby resulting in the expansion of the central-lateral region of the leaf. The findings not only provide global insights into the evolutionary features of Salicaceae GIFs but also reveal the regulatory mechanism of PeGIF3 in heteromorphic leaves of P. euphratica.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

The Euphrates Poplar Responses to Abiotic Stress and Its Unique Traits in Dry Regions of China (Xinjiang and Inner Mongolia): What Should We Know?

At the moment, drought, salinity, and low-temperature stress are ubiquitous environmental issues. In arid regions including Xinjiang and Inner Mongolia and other areas worldwide, the area of tree plantations appears to be rising, triggering tree growth. Water is a vital resource in the agricultural systems of countries impacted by aridity and salinity. Worldwide efforts to reduce quantitative yield losses on Populus euphratica by adapting tree plant production to unfavorable environmental conditions have been made in response to the responsiveness of the increasing control of water stress. Although there has been much advancement in identifying the genes that resist abiotic stresses, little is known about how plants such as P. euphratica deal with numerous abiotic stresses. P. euphratica is a varied riparian plant that can tolerate drought, salinity, low temperatures, and climate change, and has a variety of water stress adaptability abilities. To conduct this review, we gathered all available information throughout the Web of Science, the Chinese National Knowledge Infrastructure, and the National Center for Biotechnology Information on the impact of abiotic stress on the molecular mechanism and evolution of gene families at the transcription level. The data demonstrated that P. euphratica might gradually adapt its stomatal aperture, photosynthesis, antioxidant activities, xylem architecture, and hydraulic conductivity to endure extreme drought and salt stress. Our analyses will give readers an understanding of how to manage a gene family in desert trees and the influence of abiotic stresses on the productivity of tree plants. They will also give readers the knowledge necessary to improve biotechnology-based tree plant stress tolerance for sustaining yield and quality trees in China's arid regions.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

Genome-wide analysis of R2R3-MYB transcription factors reveals their differential responses to drought stress and ABA treatment in desert poplar (Populus euphratica)

The R2R3-MYB transcription factors are widely involved in the regulation of plant growth, biotic and abiotic stress responses. Meanwhile, seed germination, which is stimulated by internal and external environments, is a critical stage in the plant life cycle. However, the identification, characterization, and expression profiling of the Populus euphratica R2R3-MYB family in drought response during seed germination have been unknown. Our study attempted to identify and characterize the R2R3-MYB genes in P. euphratica (PeR2R3-MYBs) and explore how R2R3-MYBs trigger the drought and abscisic acid (ABA) response mechanism in its seedlings. Based on the analysis of comparative genomics, 174 PeR2R3-MYBs were identified and expanded driven by whole genome duplication or segment duplication events. The analysis of Ka/Ks ratios showed that, in contrast to most PeR2R3-MYBs, the other PeR2R3-MYBs were subjected to positive selection in P. euphratica. Further, the expression data of PeR2R3-MYBs under drought stress and ABA treatment, together with available functional data for Arabidopsis thaliana MYB genes, supported the hypothesis that PeR2R3-MYBs involved in response to drought are dependent or independent on ABA signaling pathway during seed germination, especially PeR2R3-MYBs with MYB binding sites (MBS) cis-element and/or tandem duplication. This study is the first report on the genome-wide analysis of PeR2R3-MYBs, as well as the other two Salicaceae species. The duplication events and differential expressions of PeR2R3-MYBs play important roles in enhancing the adaptation to drought desert environment. Our results provide a reference for prospective functional studies of R2R3-MYBs of poplars and lay the foundation for new breeding strategies to improve the drought tolerance of P. euphratica.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

Repeated turnovers keep sex chromosomes young in willows

Background

Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively.

Results

Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the “deleterious mutation load” and “sexually antagonistic selection” theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration.

Conclusions

Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

Populus × canescens root suberization in reaction to osmotic and salt stress is limited to the developing younger root tip region

Populus is a valuable and fast-growing tree species commonly cultivated for economic and scientific purposes. But most of the poplar species are sensitive to drought and salt stress. Thus, we compared the physiological effects of osmotic stress (PEG8000) and salt treatment (NaCl) on poplar roots to identify potential strategies for future breeding or genetic engineering approaches. We investigated root anatomy using epifluorescence microscopy, changes in root suberin composition and amount using gas chromatography, transcriptional reprogramming using RNA sequencing, and modifications of root transport physiology using a pressure chamber. Poplar roots reacted to the imposed stress conditions, especially in the developing younger root tip region, with remarkable differences between both types of stress. Overall, the increase in suberin content was surprisingly small, but the expression of key suberin biosynthesis genes was strongly induced. Significant reductions of the radial water transport in roots were only observed for the osmotic and not the hydrostatic hydraulic conductivity. Our data indicate that the genetic enhancement of root suberization processes in poplar might be a promising target to convey increased tolerance, especially against toxic sodium chloride.

Advanced Breeding for Biotic Stress Resistance in Poplar

Poplar is one of the most important forest trees because of its high economic value. Thanks to the fast-growing rate, easy vegetative propagation and transformation, and availability of genomic resources, poplar has been considered the model species for forest genetics, genomics, and breeding. Being a field-growing tree, poplar is exposed to environmental threats, including biotic stresses that are becoming more intense and diffused because of global warming. Current poplar farming is mainly based on monocultures of a few elite clones and the expensive and long-term conventional breeding programmes of perennial tree species cannot face current climate-change challenges. Consequently, new tools and methods are necessary to reduce the limits of traditional breeding related to the long generation time and to discover new sources of resistance. Recent advances in genomics, marker-assisted selection, genomic prediction, and genome editing offer powerful tools to efficiently exploit the Populus genetic diversity and allow enabling molecular breeding to support accurate early selection, increasing the efficiency, and reducing the time and costs of poplar breeding, that, in turn, will improve our capacity to face or prevent the emergence of new diseases or pests.

PtoNF-YC9-SRMT-PtoRD26 module regulates the high saline tolerance of a triploid poplar

BACKGROUND

Sensing and responding to stresses determine the tolerance of plants to adverse environments. The triploid Chinese white poplar is widely cultivated in North China because of its adaptation to a wide range of habitats including highly saline ones. However, its triploid genome complicates any detailed investigation of the molecular mechanisms underlying its adaptations.

RESULTS

We report a haplotype-resolved genome of this triploid poplar and characterize, using reverse genetics and biochemical approaches, a MYB gene, SALT RESPONSIVE MYB TRANSCRIPTION FACTOR (SRMT), which combines NUCLEAR FACTOR Y SUBUNIT C 9 (PtoNF-YC9) and RESPONSIVE TO DESICCATION 26 (PtoRD26), to regulate an ABA-dependent salt-stress response signaling. We reveal that the salt-inducible PtoRD26 is dependent on ABA signaling. We demonstrate that ABA or salt drives PtoNF-YC9 shuttling into the nucleus where it interacts with SRMT, resulting in the rapid expression of PtoRD26 which in turn directly regulates SRMT. This positive feedback loop of SRMT-PtoRD26 can rapidly amplify salt-stress signaling. Interference with either component of this regulatory module reduces the salt tolerance of this triploid poplar.

CONCLUSION

Our findings reveal a novel ABA-dependent salt-responsive mechanism, which is mediated by the PtoNF-YC9-SRMT-PtoRD26 module that confers salt tolerance to this triploid poplar. These genes may therefore also serve as potential and important modification targets in breeding programs.

An eco-evo-devo genetic network model of stress response

The capacity of plants to resist abiotic stresses is of great importance to agricultural, ecological and environmental sustainability, but little is known about its genetic underpinnings. Existing genetic tools can identify individual genetic variants mediating biochemical, physiological, and cellular defenses, but fail to chart an overall genetic atlas behind stress resistance. We view stress response as an eco-evo-devo process by which plants adaptively respond to stress through complex interactions of developmental canalization, phenotypic plasticity, and phenotypic integration. As such, we define and quantify stress response as the developmental change of adaptive traits from stress-free to stress-exposed environments. We integrate composite functional mapping and evolutionary game theory to reconstruct omnigenic, information-flow interaction networks for stress response. Using desert-adapted Euphrates poplar as an example, we infer salt resistance-related genome-wide interactome networks and trace the roadmap of how each SNP acts and interacts with any other possible SNPs to mediate salt resistance. We characterize the previously unknown regulatory mechanisms driving trait variation; i.e. the significance of a SNP may be due to the promotion of positive regulators, whereas the insignificance of a SNP may result from the inhibition of negative regulators. The regulator-regulatee interactions detected are not only experimentally validated by two complementary experiments, but also biologically interpreted by their encoded protein-protein interactions. Our eco-evo-devo model of genetic interactome networks provides an approach to interrogate the genetic architecture of stress response and informs precise gene editing for improving plants' capacity to live in stress environments.

Identification and Comparative Analysis of Conserved and Species-Specific microRNAs in Four Populus Sections

The conservation and diversity of microRNA (miRNA) families provide insights into the evolution of miRNA genes. However, there are few studies to explore the miRNA genes at the genus level in plants. Here, we identified 1194 miRNA loci in four Populus species P. deltoides, P. euphratica, P. tremula, and P. trichocarpa from Aigeiros, Turanga, Populus, and Tacamahaca sections, respectively, by combining de novo and homolog-based approaches. Our results indicated that a similar number of miRNA loci exist in each species (296–301 miRNA loci). Among the identified 143 miRNA families, 68 families are shared by the studied four species, and 31 families are species-specific, which might be related to local adaptation. Additionally, multiple miRNA-related single nucleotide polymorphisms (SNPs) were found, indicating that polymorphisms in pre-miRNA hairpins were likely to affect miRNA biogenesis. This study expanded the breadth and depth of miRNA annotations and provided valuable resources for further exploring the diversity and function of poplar miRNAs.

The assembled and annotated genome of the masked palm civet (Paguma larvata)

Background

The masked palm civet (Paguma larvata) acts as an intermediate host of severe acute respiratory syndrome coronavirus (SARS-CoV), which caused SARS, and transfered this virus from bats to humans. Additionally, P. larvata has the potential to carry a variety of zoonotic viruses that may threaten human health. However, genome resources for P. larvata have not been reported to date.

Findings

A chromosome-level genome assembly of P. larvata was generated using PacBio sequencing, Illumina sequencing, and Hi-C technology. The genome assembly was 2.44 Gb in size, of which 95.32% could be grouped into 22 pseudochromosomes, with contig N50 and scaffold N50 values of 12.97 Mb and 111.81 Mb, respectively. A total of 21,582 protein-coding genes were predicted, and 95.20% of the predicted genes were functionally annotated. Phylogenetic analysis of 19 animal species confirmed the close genetic relationship between P. larvata and species belonging to the Felidae family. Gene family clustering revealed 119 unique, 243 significantly expanded, and 58 significantly contracted genes in the P. larvata genome. We identified 971 positively selected genes in P. larvata, and one known human viral receptor gene PDGFRA is positively selected in P. larvata, which is required for human cytomegalovirus infection.

Conclusions

This high-quality genome assembly provides a valuable genomic resource for exploring virus–host interactions. It will also provide a reliable reference for studying the genetic bases of the morphologic characteristics, adaptive evolution, and evolutionary history of this species.

Identification of genes involved in drought tolerance in seedlings of the desert grass, Psammochloa villosa (Poaceae), based on full-length isoform sequencing and de novo assembly from short reads

Psammochloa villosa is a perennial herbaceous plant that is dominant within arid regions of the Inner Mongolian Plateau and the Qinghai-Tibet Plateau in China, where it is an endemic species and exhibits strong drought tolerance and wind resistance. To study drought tolerance in P. villosa and determine its molecular basis, we simulated high and moderate drought stress in a controlled environment and then analyzed transcriptome sequences by combining long-read sequences from a representative, wild-grown individual with short reads from the treatment groups. We obtained 184,076 high-quality isoforms as a reference and 168,650 genes (91.6%), which we were able to annotate according to public databases. Ultimately, we obtained 119,005 unigenes representing the transcriptome of P. villosa under drought stress and, among these, we identified 3089 differentially expressed genes and 1484 transcription factors. Physiologically, P. villosa that was exposed to high and moderate drought stress had reduced germination rates and shorter buds but generated more chlorophyll, which is atypical under drought stress and possibly reflects an adaptation of these plants to their arid environment. We inferred that significantly upregulated genes were annotated as 'Chlorophyll a-b binding protein' and 'Light-harvesting chlorophyll-protein' among drought and control groups. Broadly, our analyses revealed that drought stress triggered many genome-level responses, especially related to mitigation of radical oxygen species (ROS), which increase in concentration under drought stress. In particular, in the high drought stress group compared with the control, GO enrichment analysis revealed a significant enrichment of upregulated genes (n = 10) involved in mitigation of oxidative stress. Similarly, using KEGG we found significant enrichment of genes in the phenylpropanoid biosynthesis pathway (11 genes), which yields phenols that scavenge ROS. We also inferred that many genes involved in metabolism of arginine and proline, which may serve as both scavengers of ROS and osmoprotectants that interact with stress response genes based on our protein-protein interaction network analysis. We verified the relative expression levels of eight genes associated with mitigation of ROS, DNA repair, and transmembrane transporter activity using qRT-PCR, and the results were consistent with our inferences from transcriptomes. This study provides insights into the genomic and physiological basis of drought tolerance in P. villosa and represents a resource for development of the species as a forage crop via molecular breeding within arid lands.

Increased cuticular wax deposition does not change residual foliar transpiration

The effect of contrasting environmental growth conditions (in vitro tissue culture, ex vitro acclimatisation, climate chamber, greenhouse and outdoor) on leaf development, cuticular wax composition, and foliar transpiration of detached leaves of the Populus × canescens clone 84 K were investigated. Our results show that total amounts of cuticular wax increased more than 10-fold when cultivated in different growth conditions, whereas qualitative wax composition did not change. With exception of plants directly taken from tissue culture showing rapid dehydration, rates of water loss (residual foliar transpiration) of intact but detached leaves were constant and independent from growth conditions and thus independent from increasing wax amounts. Since cuticular transpiration measured with isolated astomatous P. × canescens cuticles was identical to residual foliar transpiration rates of detached leaves, our results confirm that cuticular transpiration of P. × canescens leaves can be predicted with high accuracy from residual transpiration of detached leaves after stomatal closure. Our results convincingly show that more than 10-fold increased wax amounts in P. × canescens cuticles do not lead to decreased rates of residual (cuticular) transpiration.

Transcriptome Analysis of Populus euphratica under Salt Treatment and PeERF1 Gene Enhances Salt Tolerance in Transgenic Populus alba × Populus glandulosa

Populus euphratica is mainly distributed in desert environments with dry and hot climate in summer and cold in winter. Compared with other poplars, P. euphratica is more resistant to salt stress. It is critical to investigate the transcriptome and molecular basis of salt tolerance in order to uncover stress-related genes. In this study, salt-tolerant treatment of P. euphratica resulted in an increase in osmo-regulatory substances and recovery of antioxidant enzymes. To improve the mining efficiency of candidate genes, the analysis combining both the transcriptome WGCNA and the former GWAS results was selected, and a range of key regulatory factors with salt resistance were found. The PeERF1 gene was highly connected in the turquoise modules with significant differences in salt stress traits, and the expression levels were significantly different in each treatment. For further functional verification of PeERF1, we obtained stable overexpression and dominant suppression transgenic lines by transforming into Populus alba × Populusglandulosa. The growth and physiological characteristics of the PeERF1 overexpressed plants were better than that of the wild type under salt stress. Transcriptome analysis of leaves of transgenic lines and WT revealed that highly enriched GO terms in DEGs were associated with stress responses, including abiotic stimuli responses, chemical responses, and oxidative stress responses. The result is helpful for in-depth analysis of the salt tolerance mechanism of poplar. This work provides important genes for poplar breeding with salt tolerance.

Genome-Wide Comprehensive Analysis of the GASA Gene Family in Populus

Gibberellic acid-stimulated Arabidopsis (GASA) proteins, as cysteine-rich peptides (CRPs), play roles in development and reproduction and biotic and abiotic stresses. Although the GASA gene family has been identified in plants, the knowledge about GASAs in Populus euphratica, the woody model plant for studying abiotic stress, remains limited. Here, we referenced the well-sequenced Populus trichocarpa genome, and identified the GASAs in the whole genome of P. euphratica and P. trichocarpa. 21 candidate genes in P. trichocarpa and 19 candidate genes in P. euphratica were identified and categorized into three subfamilies by phylogenetic analysis. Most GASAs with signal peptides were located extracellularly. The GASA genes in Populus have experienced multiple gene duplication events, especially in the subfamily A. The evolution of the subfamily A, with the largest number of members, can be attributed to whole-genome duplication (WGD) and tandem duplication (TD). Collinearity analysis showed that WGD genes played a leading role in the evolution of GASA genes subfamily B. The expression patterns of P. trichocarpa and P. euphratica were investigated using the PlantGenIE database and the real-time quantitative PCR (qRT-PCR), respectively. GASA genes in P. trichocarpa and P. euphratica were mainly expressed in young tissues and organs, and almost rarely expressed in mature leaves. GASA genes in P. euphratica leaves were also widely involved in hormone responses and drought stress responses. GUS activity assay showed that PeuGASA15 was widely present in various organs of the plant, especially in vascular bundles, and was induced by auxin and inhibited by mannitol dramatically. In summary, this present study provides a theoretical foundation for further research on the function of GASA genes in P. euphratica.

Chromosome-length genome assemblies of six legume species provide insights into genome organization, evolution, and agronomic traits for crop improvement

INTRODUCTION

Legume crops are an important source of protein and oil for human health and in fixing atmospheric N₂ for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement.

OBJECTIVES

Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution.

METHODS

We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies.

RESULTS

We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of "QTL-hotspot" identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, 'Legumepedia'.

CONCLUSION

This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Overexpression of antisense phosphatase 2C affords cold resistance in hybrid Populus davidiana × Populus bolleana

BACKGROUND

Overexpression of the abiotic and biotic stress-resistance genes of the plant signaling pathway is well known for its significant role in the regulation of plant growth and enhancement of the productivity of agricultural land under changing climatic conditions.

OBJECTIVES

This research aimed to clone Populus davidiana × Populus bolleana PP2C (PdPP2C) gene and analyze its structure and function, and downregulate PdPP2C by overexpression of its antisense PdPP2C (AS-PdPP2C) gene for enhancing cold resistance in transgenic lines of hybrid P. davidiana × P. bolleana.

METHODS

PdPP2C was cloned and transformed to identify its function, and its antisense was overexpressed via downregulation to increase the cold resistance in transgenic lines of hybrid P. davidiana × P. bolleana.

RESULTS

Antisense inhibition of protein phosphatase 2C accelerates the cold acclimation of Poplar (P. davidiana × P. bolleana) gene in terms of antifreeze.

CONCLUSION

PdPP2C was expressed in the roots, stems, and leaves of P. davidiana × P. bolleana, and the expression was higher in the leaves. The expression of PdPP2C was also significantly downregulated at low-temperature (0 °C and 4 °C) stress. The relative conductivity and malondialdehyde content of non-transgenic lines were higher than those of AS-PdPP2C lines after 2 days of cold treatment at - 1 °C. The leaves of the transgenic lines were not wilted and showed no chlorosis compared with those of the non-transgenic lines. The AS-PdPP2C transgenic lines also showed higher freezing resistance than the non-transgenic lines. AS-PdPP2C participated in the regulation of freezing resistance.

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Genome Assembly of Salicaceae Populus deltoides (Eastern Cottonwood) I-69 Based on Nanopore Sequencing and Hi-C Technologies

Populus deltoides has important ecological and economic values, widely used in poplar breeding programs due to its superior characteristics such as rapid growth and resistance to disease. Although the genome sequence of P. deltoides WV94 is available, the assembly is fragmented. Here, we reported an improved chromosome-level assembly of the P. deltoides cultivar I-69 by combining Nanopore sequencing and chromosome conformation capture (Hi-C) technologies. The assembly was 429.3 Mb in size and contained 657 contigs with a contig N50 length of 2.62 Mb. Hi-C scaffolding of the contigs generated 19 chromosome-level sequences, which covered 97.4% (418 Mb) of the total assembly size. Moreover, repetitive sequences annotation showed that 39.28% of the P. deltoides genome was composed of interspersed elements, including retroelements (23.66%), DNA transposons (6.83%), and unclassified elements (8.79%). We also identified a total of 44 362 protein-coding genes in the current P. deltoides assembly. Compared with the previous genome assembly of P. deltoides WV94, the current assembly had some significantly improved qualities: the contig N50 increased 3.5-fold and the proportion of gaps decreased from 3.2% to 0.08%. This high-quality, well-annotated genome assembly provides a reliable genomic resource for identifying genome variants among individuals, mining candidate genes that control growth and wood quality traits, and facilitating further application of genomics-assisted breeding in populations related to P. deltoides.

Chromosomal-Level Reference Genome of the Neotropical Tree Jacaranda mimosifolia D. Don

Jacaranda mimosifolia D. Don is a deciduous tree widely cultivated in the tropics and subtropics of the world. It is famous for its beautiful blue flowers and pinnate compound leaves. In addition, this tree has great potential in environmental monitoring, soil quality improvement, and medicinal applications. However, a genome resource for J. mimosifolia has not been reported to date. In this study, we constructed a chromosome-level genome assembly of J. mimosifolia using PacBio sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼707.32 Mb in size, 688.76 Mb (97.36%) of which could be grouped into 18 pseudochromosomes, with contig and scaffold N50 values of 16.77 and 39.98 Mb, respectively. A total of 30,507 protein-coding genes were predicted, 95.17% of which could be functionally annotated. Phylogenetic analysis among 12 plant species confirmed the close genetic relationship between J. mimosifolia and Handroanthus impetiginosus. Gene family clustering revealed 481 unique, 103 significantly expanded, and 16 significantly contracted gene families in the J. mimosifolia genome. This chromosome-level genome assembly of J. mimosifolia will provide a valuable genomic resource for elucidating the genetic bases of the morphological characteristics, adaption evolution, and active compounds biosynthesis of J. mimosifolia.

Modeling genome-wide by environment interactions through omnigenic interactome networks

How genes interact with the environment to shape phenotypic variation and evolution is a fundamental question intriguing to biologists from various fields. Existing linear models built on single genes are inadequate to reveal the complexity of genotype-environment (G-E) interactions. Here, we develop a conceptual model for mechanistically dissecting G-E interplay by integrating previously disconnected theories and methods. Under this integration, evolutionary game theory, developmental modularity theory, and a variable selection method allow us to reconstruct environment-induced, maximally informative, sparse, and casual multilayer genetic networks. We design and conduct two mapping experiments by using a desert-adapted tree species to validate the biological application of the model proposed. The model identifies previously uncharacterized molecular mechanisms that mediate trees' response to saline stress. Our model provides a tool to comprehend the genetic architecture of trait variation and evolution and trace the information flow of each gene toward phenotypes within omnigenic networks.

Improved de novo chromosome-level genome assembly of the vulnerable walnut tree Juglans mandshurica reveals gene family evolution and possible genome basis of resistance to lesion nematode

Manchurian walnut (Juglans mandshurica Maxim.) is a synonym of J. cathayensis, a diploid, vulnerable, temperate deciduous tree valued for its wood and nut. It is also valued as a rootstock for Juglans regia because of its reported tolerance of lesion nematode. Reference genomes are available for several Juglans species, our goal was to produce a de novo, chromosome-level assembly of the J. mandshurica genome. Here, we reported an improved assembly of J. mandshurica with a contig N50 size of 6.49 Mb and a scaffold N50 size of 36.1 Mb. The total genome size was 548 Mb encoding 29,032 protein coding genes which were annotated. The collinearity analysis showed that J. mandshurica and J. regia originated from a common ancestor, with both species undergoing two WGD events. A genomic comparison showed that J. mandshurica was missing 1657 genes found in J. regia, and J. mandshurica includes 2827 genes not found in of the J. regia genome. The J. mandshurica contained 1440 unique paralogues that were highly enriched for flavonoid biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction. Four gene families related to disease resistance notable contraction (rapidly evolving; LEA, WAK, PPR, and PR) in J. mandshurica compared to eight species. JmaPR10 and JmaPR8 contained three orthologous gene pairs with J. regia that were highly expressed in root bark. JmaPR10 is a strong candidate gene for lesion nematodes resistance in J. mandshurica. The J. mandshurica genome should be a useful resource for study of the evolution, breeding, and genetic variation in walnuts (Juglans).

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Bioinformatic analysis of chromatin organization and biased expression of duplicated genes between two poplars with a common whole-genome duplication

The nonrandom three-dimensional organization of chromatin plays an important role in the regulation of gene expression. However, it remains unclear whether this organization is conserved and whether it is involved in regulating gene expression during speciation after whole-genome duplication (WGD) in plants. In this study, high-resolution interaction maps were generated using high-throughput chromatin conformation capture (Hi-C) techniques for two poplar species, Populus euphratica and Populus alba var. pyramidalis, which diverged ~14 Mya after a common WGD. We examined the similarities and differences in the hierarchical chromatin organization between the two species, including A/B compartment regions and topologically associating domains (TADs), as well as in their DNA methylation and gene expression patterns. We found that chromatin status was strongly associated with epigenetic modifications and gene transcriptional activity, yet the conservation of hierarchical chromatin organization across the two species was low. The divergence of gene expression between WGD-derived paralogs was associated with the strength of chromatin interactions, and colocalized paralogs exhibited strong similarities in epigenetic modifications and expression levels. Thus, the spatial localization of duplicated genes is highly correlated with biased expression during the diploidization process. This study provides new insights into the evolution of chromatin organization and transcriptional regulation during the speciation process of poplars after WGD.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

A chromosome-level genome assembly for the tertiary relict plant Tetracentron sinense oliv. (trochodendraceae)

Tetracentron sinense and Trochodendron aralioides are two Tertiary relict species of large trees in the family Trochodendraceae with narrow distributions on the mainland and islands of eastern Asia. They belong to the order Trochodendrales, which is one of the four early-diverged eudicot lineages. These two relict species provide a good system in which to examine genomic changes that occurred as they survived during repeated climatic oscillations in the Quaternary. We sequenced the genome of Te. sinense and compared it with that of Tr. aralioides. We found that Te. sinense has a smaller genome size (986.3 Mb) than that of Tr. aralioides (1610 Mb). Repetitive elements made the major contribution to the contrasting genome sizes in the two species, with most bursts of repeats occurring within the past four million years when the climate oscillated greatly. These species share two rounds of whole-genome duplications. The mainland species Te. sinense had a larger effective population size than the island species Tr. aralioides after the largest glaciation during the Quaternary climatic oscillation. However, soon after this recovery stage, the effective population sizes of both species continued to decrease, although the current effective population size of Te. sinense is still larger than that of Tr. aralioides. We recovered three distinctly diverged clades through resequencing the genomes of 50 individuals across the distributional range of Te. sinense in China. Our results provide an important genomic resource with which to examine early trait evolution in the core eudicots and assist efforts to conserve this relict tree species.

The Genome Sequence of Alpine Megacarpaea delavayi Identifies Species-Specific Whole-Genome Duplication

Megacarpaea delavayi (Brassicaceae), a plant found the high mountains of southwest China at high altitudes (3000–4800 m), is used as a vegetable or medicine. Here, we report a draft genome for this species. The assembly genome of M. delavayi is 883 Mb, and 61.59% of the genome is composed of repeat sequences. Annotation of the genome identified a total of 41,114 protein-coding genes. We found that M. delavayi experienced an independent whole-genome duplication (WGD), paralleling those independent WGDs in Iberis, Biscutella, and Anastatica in the early Miocene. Phylogenetic analyses based on the single-copy genes confirmed the position of the genus Megacarpaea within the expanded lineage II of the family and resolved its basal divergence to a subclade consisting of Anastatica, Iberis, and Biscutella. Species-specific and fast-evolving genes in M. delavayi are mainly involved in “DNA repair” and “response to UV-B radiation.” These genetic changes may together help this species survive in high-altitude environments. The reference genome reported here provides a valuable resource for studying adaptation of this and other alpine plants to the high-altitude habitats.