Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica

Comprehensive Genomic Analysis of Trihelix Transcription Factor Genes and Their Expression Underlying Abiotic Stress in Euphrates Poplar (Populus euphratica)

Trihelix transcription factors (TTFs) are light-sensitive proteins characterized by a triple-helix structure that play a crucial role in regulating plant growth and development, especially in response to abiotic stressors, such as drought and salinity. This intriguing family of proteins has been the focus of extensive functional studies across various plant species. Despite their recognized significance, the trihelix family in Populus euphratica has not been thoroughly explored, warranting more attention. This study identifies 35 full-length trihelix genes in Populus euphratica, which are grouped into five categories (GT-1, GT-γ, GT-2, SIP1, and SH4) based on their conserved motifs and structural similarities, and these genes are unevenly distributed across 19 linkage groups on the chromosomes. A syntenic analysis was conducted in P. euphratica, comparing it to various other species. The promoters of P. euphratica contain numerous stress-responsive cis-elements, indicating the potential for these trihelix genes to respond to abiotic stress. RT-qPCR analysis discovered significant induction of the trihelix gene family in response to drought and salt stress, with 21 PeuTTF genes exhibiting distinct expression levels under drought conditions and five PeuTTF genes responsive to salt stress. Notably, heightened expression of PeuTTF6, PeuTTF9, and PeuTTF20 was observed in both roots and leaves during drought stress, suggesting that TTF expression is connected to the plant's response to such conditions. Additionally, significant increases in expression were noted for PeuTTF2, PeuTTF31, and PeuTTF32, which may be convoluted in the response to salt stress. These discoveries highlight the role that PeuTTF genes play in improving drought tolerance in P. euphratica plants. We offer new perspectives on the evolutionary trends and variants of PeuTTF genes in P. euphratica, and we establish the groundwork for understanding the functional properties of PeuTTF genes under salt-stressed and drought-stressed conditions. This study provides opportunities for the advancement of desert poplar agriculture and may have wider ramifications for tree plant breeding techniques targeted at improving tree performance and durability, particularly in dry areas.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

Phylogeographic analysis reveals multiple origins of the desert shrub Reaumuria songarica in northern Xinjiang, involving homoploid and tetraploid hybrids

Hybrid speciation plays an important role in species diversification. The establishment of reproductive isolation is crucial for hybrid speciation, and the identification of diverse types of hybrids, particularly homoploid hybrid species, contributes to a comprehensive understanding of this process. Reaumuria songarica is a constructive shrub widespread in arid Central Asia. Previous studies have inferred that the R. songarica populations in the Gurbantunggut Desert (GuD) originated from homoploid hybridizations between its eastern and western lineages and may have evolved into an incipient species. To further elucidate the genetic composition of different hybrid populations and to determine the species boundary of this hybrid lineage, we investigated the overall phylogeographic structure of R. songarica based on variation patterns of five cpDNA and one nrITS sequences across 32 populations. Phylogenetic analyses demonstrated that within the GuD lineage, the Wuerhe population evolved directly from ancestral lineages, whereas the others originated from hybridizations between the eastern and western lineages. PCoA and genetic barrier analysis supported the subdivision of the GuD lineage into the southern (GuD-S) and northern (GuD-N) groups. Populations in the GuD-S group had a consistent genetic composition and the same ancestral female parent, indicating that they belonged to a homoploid hybrid lineage. However, the GuD-N group experienced genetic admixture of the eastern and western lineages on nrITS and cpDNA, with some populations inferred to be allopolyploid based on ploidy data. Based on cpDNA haplotypes, BEAST analyses showed that the GuD-S and GuD-N groups originated after 0.5 Ma. Our results suggest that multiple expansions and contractions of GuD, driven by Quaternary climatic oscillations and the Kunlun-Yellow River tectonic movement, are important causes of the complex origins of R. songarica populations in northern Xinjiang. This study highlights the complex origins of the Junggar Basin flora and the underappreciated role of hybridization in increasing its species diversity.

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.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

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

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

Differences in ecophysiological responses of Populus euphratica females and males exposed to salinity and alkali stress

Soil salinity is usually accompanied by alkalization in northwest China, and they both negatively impact plant growth and result in severe ecological problems. Some studies have reported tree responses to salinity or alkali stress alone, however, the interactive salinity and alkali effects are still unclear, especially in dioecious trees. In this study, we measured growth, morphology, leaf stomata, gas exchange, carbon isotope composition (δ¹³C), total soluble sugar and starch contents, Na⁺ accumulation and allocation, oxidative stress, and antioxidants of female and male Populus euphratica seedlings in response to salinity, alkali and their interaction. Our study showed no significant sexual differences in studied traits under control conditions. In addition, P. euphratica females showed greater inhibitory and negative effects, such as bigger decreases in growth and gas exchange, lower stomatal density and water use efficiency (as described by δ¹³C), and lower levels of soluble sugars and antioxidant enzyme activities compared with males under salinity, alkali and interactive stress conditions. Furthermore, P. euphratica males had a greater ability of ion exclusion and Na ⁺ transport restriction. For example, males allocated more Na⁺ to stems and roots than females, whereas females had higher Na⁺ contents in leaves under stress conditions. In conclusion, our results indicated that P. euphratica males have superior resistance and they perform better than females under salinity, alkali and their interactive stress conditions.

PeCLH2 Gene Positively Regulate Salt Tolerance in Transgenic Populus alba × Populus glandulosa

Salt is an important environmental stress factor, which seriously affects the growth, development and distribution of plants. Chlorophyllase plays an important role in stress response. Nevertheless, little is known about the physiological and molecular mechanism of chlorophyll (Chlase, CLH) genes in plants. We cloned PeCLH2 from Populus euphratica and found that PeCLH2 was differentially expressed in different tissues, especially in the leaves of P. euphratica. To further study the role of PeCLH2 in salt tolerance, PeCLH2 overexpression and RNA interference transgenic lines were established in Populus alba × Populus glandulosa and used for salt stress treatment and physiologic indexes studies. Overexpressing lines significantly improved tolerance to salt treatment and reduced reactive oxygen species production. RNA interference lines showed the opposite. Transcriptome analysis was performed on leaves of control and transgenic lines under normal growth conditions and salt stress to predict genes regulated during salt stress. This provides a basis for elucidating the molecular regulation mechanism of PeCLH2 in response to salt stress and improving the tolerance of poplar under salt stress.

Transcription factor PagLBD21 functions as a repressor of secondary xylem development in Populus

During secondary growth in trees, xylem cells differentiated from cambium cells begin to synthesize secondary cell walls that are primarily composed of cellulose, hemicellulose and lignin, and are deposited between primary cell walls and plasma membranes, leading to wood formation. Identification of regulatory genes functioning during this developmental process is valuable for increasing wood production. In this study, we functionally characterized an LBD (LATERAL ORGAN BOUNDARIES DOMAIN) transcription factor PagLBD21 as a repressor of secondary xylem development in Populus. Compared to wild type plants, transgenic plants overexpressing PagLBD21 (PagLBD21OE) exhibited decreased xylem widths in cross-sections. Consistent with the functional analysis, RNA sequencing (RNA-seq) analysis revealed that genes functioning in xylem development and secondary cell wall biosynthesis pathways were significantly down-regulated in PagLBD21OE plants. We also performed DNA affinity purification followed by sequencing (DAP-seq) to identify genome-wide target genes of PagLBD21. Furthermore, we compared the RNA-seq and DAP-seq datasets of PagLBD21 and PagLBD3, and the results showed that there was a significant overlap between their target genes, suggesting these two LBD transcription factors are functionally redundant during secondary growth.

Diversification and historical demography of Haloxylon ammodendron in relation to Pleistocene climatic oscillations in northwestern China

The influence of aridification and climatic oscillations on the genetic diversity and evolutionary processes of organisms during the Quaternary in northwestern China is examined using Haloxylon ammodendron. Based on the variation of two cpDNA regions (trnS-trnG and trnV) and one nDNA sequence (ITS1-ITS4) in 420 individuals from 36 populations, the spatial genetic structure and demographic history of H. ammodendron in arid China is examined. Median-joining network and Bayesian inference trees enabled the identification of three diverged lineages within H. ammodendron from 24 different haplotypes and 16 ribotypes, distributed across western (Xinjiang), eastern (Gansu and Inner Mongolia) and southern (Qinghai) regions. AMOVA analysis demonstrated that more than 80% of observed genetic variation related to lineage split was based on cpDNA and nDNA variation. Allopatric divergence among the three groups was mainly triggered by geographical isolation due to Xingxingxia rock and uplift of the Qilian Mountains during the Quaternary. Local adaptive differentiation among western, eastern and southern groups occurred due to gene flow obstruction resulting from arid landscape fragmentation accompanied by local environmental heterogeneity of different geographical populations. The southern margin of the Junggar Basin and the Tengger Desert possibly served as two independent glacial refugia for H. ammodendron. The distribution of genetic variation, coupled with SDMs and LCP results, indicated that H. ammodendron probably moved northward along the Junggar Basin and westward along Tengger Desert at the end of the last glacial maximum; postglacial re-colonization was probably westward and southward along the Hexi Corridor.

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.

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.

Phylogeography Reveals Geographic and Environmental Factors Driving Genetic Differentiation of Populus sect. Turanga in Northwest China

Populus sect. Turanga (hereafter referred to as "Populus"), including Populus euphratica and Populus pruinosa, are the predominant tree species in desert riparian forests in northwestern China. These trees play key roles in maintaining ecosystem balance, curbing desertification, and protecting biodiversity. However, the distribution area of Populus forests has been severely diminished and degraded in recent years due to increased habitat destruction and human activity. Understanding the genetic diversity among Populus individuals and populations is essential for designing conservation strategies, but comprehensive studies of their genetic diversity in northwest China are lacking. Here, we assessed the population structures and genetic diversity of 1,620 samples from 85 natural populations of Populus (59 P. euphratica and 26 P. pruinosa populations) covering all of northwestern China using 120 single nucleotide polymorphism (SNP) markers. Analysis of population structure revealed significant differentiation between these two sister species and indicated that strong geographical distribution patterns, a geographical barrier, and environmental heterogeneity shaped the extant genetic patterns of Populus. Both P. euphratica and P. pruinosa populations in southern Xinjiang had higher genetic diversity than populations in other clades, perhaps contributing to local geographic structure and strong gene flow. Analysis of molecular variance (AMOVA) identified 15% variance among and 85% variance within subpopulations. Mantel tests suggested that the genetic variation among P. euphratica and P. pruinosa populations could be explained by both geographical and environmental distance. The genetic diversity of P. euphratica showed a significant negative correlation with latitude and longitude and a positive correlation with various environmental factors, such as precipitation of warmest quarter and driest month, temperature seasonality, and annual mean temperature. These findings provide insights into how the genetic differentiation of endangered Populus species was driven by geographical and environmental factors, which should be helpful for designing strategies to protect these genetic resources in the future.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.