Computational dissection of genetic variation modulating the response of multiple photosynthetic phenotypes to the light environment.

BACKGROUND: The expression of biological traits is modulated by genetics as well as the environment, and the level of influence exerted by the latter may vary across characteristics. Photosynthetic traits in plants are complex quantitative traits that are regulated by both endogenous genetic factors and external environmental factors such as light intensity and CO2 concentration. The specific processes impacted occur dynamically and continuously as the growth of plants changes. Although studies have been conducted to explore the genetic regulatory mechanisms of individual photosynthetic traits or to evaluate the effects of certain environmental variables on photosynthetic traits, the systematic impact of environmental variables on the dynamic process of integrated plant growth and development has not been fully elucidated. RESULTS: In this paper, we proposed a research framework to investigate the genetic mechanism of high-dimensional complex photosynthetic traits in response to the light environment at the genome level. We established a set of high-dimensional equations incorporating environmental regulators to integrate functional mapping and dynamic screening of gene‒environment complex systems to elucidate the process and pattern of intrinsic genetic regulatory mechanisms of three types of photosynthetic phenotypes of Populus simonii that varied with light intensity. Furthermore, a network structure was established to elucidate the crosstalk among significant QTLs that regulate photosynthetic phenotypic systems. Additionally, the detection of key QTLs governing the response of multiple phenotypes to the light environment, coupled with the intrinsic differences in genotype expression, provides valuable insights into the regulatory mechanisms that drive the transition of photosynthetic activity and photoprotection in the face of varying light intensity gradients. CONCLUSIONS: This paper offers a comprehensive approach to unraveling the genetic architecture of multidimensional variations in photosynthetic phenotypes, considering the combined impact of integrated environmental factors from multiple perspectives.

Enhanced genome-wide association reveals the role of YABBY11-NGATHA-LIKE1 in leaf serration development of Populus.

Leaf margins are complex plant morphological features that contribute to leaf shape diversity, which affects plant structure, yield, and adaptation. Although several leaf margin regulators have been identified to date, the genetic basis of their natural variation has not been fully elucidated. In this study, we profiled two distinct leaf morphology types (serrated and smooth) using the persistent homology mathematical framework (PHMF) in two poplar species (Populus tomentosa and Populus simonii, respectively). A combined genome-wide association study (GWAS) and expression quantitative trait nucleotide (eQTN) mapping were applied to create a leaf morphology control module using data from P. tomentosa and P. simonii populations. Natural variation in leaf margins was associated with YABBY11 (YAB11) transcript abundance in poplar. In P. tomentosa, PtoYAB11 carries a premature stop codon (PtoYAB11PSC), resulting in the loss of its positive regulation of NGATHA-LIKE1 (PtoNGAL-1) and RIBULOSE BISPHOSPHATE CARBOXYLASE LARGE SUBUNIT (PtoRBCL). Overexpression of PtoYAB11PSC promoted serrated leaf margins, enlarged leaves, enhanced photosynthesis, and increased biomass. Overexpression of PsiYAB11 in P. tomentosa promoted smooth leaf margins, higher stomatal density, and greater light damage repair ability. In poplar, YAB11-NGAL1 is sensitive to environmental conditions, acts as a positive regulator of leaf margin serration, and may also link environmental signaling to leaf morphological plasticity.

Identification of key genes involved in lignin and flavonoid accumulation during Tilia tuan seed maturation.

KEY MESSAGE: Transcriptomics and phenotypic data analysis identified 24 transcription factors (TFs) that play key roles in regulating the competitive accumulation of lignin and flavonoids. Tilia tuan Szyszyl. (T. tuan) is a timber tree species with important ecological and commercial value. However, its highly lignified pericarp results in a low seed germination rate and a long dormancy period. In addition, it is unknown whether there is an interaction between the biosynthesis of flavonoids and lignin as products of the phenylpropanoid pathway during seed development. To explore the molecular regulatory mechanism of lignin and flavonoid biosynthesis, T. tuan seeds were harvested at five stages (30, 60, 90, 120, and 150 days after pollination) for lignin and flavonoid analyses. The results showed that lignin accumulated rapidly in the early and middle stages (S1, S3, and S4), and rapid accumulation of flavonoids during the early and late stages (S1 and S5). High-throughput RNA sequencing analysis of developing seeds identified 50,553 transcripts, including 223 phenylpropanoid biosynthetic pathway genes involved in lignin accumulation grouped into 3 clusters, and 106 flavonoid biosynthetic pathway genes (FBPGs) grouped into 2 clusters. Subsequent WGCNA and time-ordered gene co-expression network (TO-GCN) analysis revealed that 24 TFs (e.g., TtARF2 and TtWRKY15) were involved in flavonoids and lignin biosynthesis regulation. The transcriptome data were validated by qRT-PCR to analyze the expression profiles of key enzyme-coding genes. This study revealed that there existed a competitive relationship between flavonoid and lignin biosynthesis pathway during the development of T. tuan seeds, that provide a foundation for the further exploration of molecular mechanisms underlying lignin and flavonoid accumulation in T. tuan seeds.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa.

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Dynamic physiological and transcriptome changes reveal a potential relationship between the circadian clock and salt stress response in Ulmus pumila.

Despite the important role the circadian clock plays in numerous critical physiological responses in plants, such as hypocotyl elongation, leaf movement, stomatal opening, flowering, and stress responses, there have been no investigations into the effect of the circadian clock on physiological and transcriptional networks under salt stress. Ulmus pumila L. has been reported to tolerate 100-150 mM NaCl treatment. We measured the diurnal variation in photosynthesis and chlorophyll fluorescence parameters and performed a time-course transcriptome analysis of 2-years-old U. pumila seedlings under salt treatment to dissect the physiological regulation and potential relationship between the circadian network and the salt stress response. Seedlings in 150 mM NaCl treatment exhibited salt-induced physiological enhancement compared to the control group. A total of 7009 differentially expressed unigenes (DEGs) were identified under salt stress, of which 16 DEGs were identified as circadian rhythm-related DEGs (crDEGs). Further analysis of dynamic expression changes revealed that DEGs involved in four crucial pathways-photosynthesis, thiamine metabolism, abscisic acid synthesis and metabolism, and the hormone-MAPK signal crosstalk pathway-are closely related to the circadian clock. Finally, we constructed a co-expression network between the circadian clock and these four crucial pathways. Our results help shed light on the molecular link between the circadian network and salt stress tolerance in U. pumila.

The Genetic Basis of Phosphorus Utilization Efficiency in Plants Provide New Insight into Woody Perennial Plants Improvement.

Soil nutrient restrictions are the main environmental conditions limiting plant growth, development, yield, and quality. Phosphorus (P), an essential macronutrient, is one of the most significant factors that vastly restrains the growth and development of plants. Although the total P is rich in soil, its bio-available concentration is still unable to meet the requirements of plants. To maintain P homeostasis, plants have developed lots of intricate responsive and acclimatory mechanisms at different levels, which contribute to administering the acquisition of inorganic phosphate (Pi), translocation, remobilization, and recycling of Pi. In recent years, significant advances have been made in the exploration of the utilization of P in annual plants, while the research progress in woody perennial plants is still vague. In the meanwhile, compared to annual plants, relevant reviews about P utilization in woody perennial plants are scarce. Therefore, based on the importance of P in the growth and development of plants, we briefly reviewed the latest advances on the genetic and molecular mechanisms of plants to uphold P homeostasis, P sensing, and signaling, ion transporting and metabolic regulation, and proposed the possible sustainable management strategies to fasten the P cycle in modern agriculture and new directions for future studies.

Integration of genome wide association studies and co-expression networks reveal roles of PtoWRKY 42-PtoUGT76C1-1 in trans-zeatin metabolism and cytokinin sensitivity in poplar.

Cytokinins are important for in vitro shoot regeneration in plants. Cytokinin N-glucosides are produced via an irreversible glycosylation pathway, which regulates the endogenous cytokinin content. Although cytokinin N-glucoside pathways have been uncovered in higher plants, no regulator has been identified to date. We performed a metabolome genome-wide association study (mGWAS), weighted gene co-expression network analysis (WGCNA), and expression quantitative trait nucleotide (eQTN) mappings to build a core triple genetic network (mGWAS-gene expression-phenotype) for the trans-zeatin N-glucoside (ZNG) metabolite using data from 435 unrelated Populus tomentosa individuals. Variation of the ZNG level in poplar is attributed to the differential transcription of PtoWRKY42, a member of WRKY multigene family group IIb. Functional analysis revealed that PtoWRKY42 negatively regulated ZNG accumulation by binding directly to the W-box of the UDP-glycosyltransferase 76C 1-1 (PtoUGT761-1) promoter. Also, PtoWRKY42 was strongly induced by leaf senescence, 6-BA, wounding, and salt stress, resulting in a reduced ZNG level. We identified PtoWRKY42, a negative regulator of cytokinin N-glucosides, which contributes to the natural variation in ZNG level and mediates ZNG accumulation by directly modulating the key glycosyltransferase gene PtoUGT76C1-1.

Miniature inverted repeat transposable elements cis-regulate circular RNA expression and promote ethylene biosynthesis, reducing heat tolerance in Populus tomentosa.

Transposable elements (TEs) and their reverse complementary sequence pairs (RCPs) are enriched around loci that produce circular RNAs (circRNAs) in plants. However, the function of these TE-RCP pairs in modulating circRNA expression remains elusive. Here, we identified 4609 circRNAs in poplar (Populus tomentosa) and showed that miniature inverted repeat transposable elements (MITEs)-RCPs were enriched in circRNA flanking regions. Moreover, we used expression quantitative trait nucleotide (eQTN) mapping to decipher the cis-regulatory role of MITEs. eQTN results showed that 14 single-nucleotide polymorphisms (SNPs) were significantly associated with Circ_0000408 and Circ_0003418 levels and the lead associated SNPs were located in MITE-RCP regions, indicating that MITE-RCP sequence variations affect exon circularization. Overexpression and knockdown analysis showed that Circ_0003418 positively modulated its parental gene, which encodes the RING-type E3 ligase XBAT32, and specifically increased the expression of the PtoXBAT32.5 transcript variant, which lacks the E3 ubiquitin ligase domain. Under heat stress, PtoXBAT32.5 expression was induced with up-regulation of Circ_0003418, resulting in increased production of ethylene and peroxidation of membrane lipids. Our findings thus reveal the cis-regulatory mechanism by which a MITE-RCP pair affects circRNA abundance in poplar and indicate that Circ_0003418 is a negative regulator of poplar heat tolerance via the ubiquitin-mediated protein modification pathway.

Transcriptome analysis and association mapping reveal the genetic regulatory network response to cadmium stress in Populus tomentosa.

Long non-coding RNAs (lncRNAs) play essential roles in plant abiotic stress responses, but the response of lncRNA-mediated genetic networks to cadmium (Cd) treatment remain elusive in trees, the promising candidates for phytoremediation of Cd contamination. We identified 172 Cd-responsive lncRNAs and 295 differentially expressed target genes in the leaves of Cd-treated Populus tomentosa. Functional annotation revealed that these lncRNAs were involved in various processes, including photosynthesis, hormone regulation, and phenylalanine metabolism. Association studies identified 78 significant associations, representing 14 Cd-responsive lncRNAs and 28 target genes for photosynthetic and leaf physiological traits. Epistasis uncovered 83 pairwise interactions among these traits, revealing Cd-responsive lncRNA-mediated genetic networks for photosynthesis and leaf physiology in P. tomentosa. We focused on the roles of two Cd-responsive lncRNA-gene pairs, MSTRG.22608.1-PtoMYB73 and MSTRG.5634.1-PtoMYB27, in Cd tolerance of Populus, and detected insertions/deletions within lncRNAs as polymorphisms driving target gene expression. Genotype analysis of lncRNAs and heterologous overexpression of PtoMYB73 and PtoMYB27 in Arabidopsis indicated their effects on enhancing Cd tolerance, photosynthetic rate, and leaf growth, and the potential interaction mechanisms of PtoMYB73 with abiotic stresses. Our study identifies the genetic basis for the response of Populus to Cd treatment, facilitating genetic improvement of Cd tolerance in trees.

Gene Coexpression Network Analysis Indicates that Hub Genes Related to Photosynthesis and Starch Synthesis Modulate Salt Stress Tolerance in Ulmus pumila.

Ulmus pumila L. is an excellent afforestation and biofuel tree that produces high-quality wood, rich in starch. In addition, U. pumila is highly adaptable to adverse environmental conditions, which is conducive to its utilization for vegetating saline soils. However, little is known about the physiological responses and transcriptional regulatory network of U. pumila under salt stress. In this study, we exposed five main cultivars in saline-alkali land (Upu2, 5, 8, 11, and 12) to NaCl stress. Of the five cultivars assessed, Upu11 exhibited the highest salt resistance. Growth and biomass accumulation in Upu11 were promoted under low salt concentrations (<150 mM). However, after 3 months of continuous treatment with 150 mM NaCl, growth was inhibited, and photosynthesis declined. A transcriptome analysis conducted after 3 months of treatment detected 7009 differentially expressed unigenes (DEGs). The gene annotation indicated that these DEGs were mainly related to photosynthesis and carbon metabolism. Furthermore, PHOTOSYNTHETIC ELECTRON TRANSFERH (UpPETH), an important electron transporter in the photosynthetic electron transport chain, and UpWAXY, a key gene controlling amylose synthesis in the starch synthesis pathway, were identified as hub genes in the gene coexpression network. We identified 25 and 62 unigenes that may interact with PETH and WAXY, respectively. Overexpression of UpPETH and UpWAXY significantly increased the survival rates, net photosynthetic rates, biomass, and starch content of transgenic Arabidopsis plants under salt stress. Our findings clarify the physiological and transcriptional regulators that promote or inhibit growth under environmental stress. The identification of salt-responsive hub genes directly responsible for photosynthesis and starch synthesis or metabolism will provide targets for future genetic improvements.

High-Temperature-Responsive Poplar lncRNAs Modulate Target Gene Expression via RNA Interference and Act as RNA Scaffolds to Enhance Heat Tolerance.

High-temperature stress is a threat to plant development and survival. Long noncoding RNAs (lncRNAs) participate in plant stress responses, but their functions in the complex stress response network remain unknown. Poplar contributes to terrestrial ecological stability. In this study, we identified 204 high-temperature-responsive lncRNAs in an abiotic stress-tolerant poplar (Populus simonii) species using strand-specific RNA sequencing (ssRNA-seq). Mimicking overexpressed and repressed candidate lncRNAs in poplar was used to illuminate their regulation pattern on targets using nano sheet mediation. These lncRNAs were predicted to target 185 genes, of which 100 were cis genes and 119 were trans genes. Gene Ontology enrichment analysis showed that anatomical structure morphogenesis and response to stress and signaling were significantly enriched. Among heat-responsive LncRNAs, TCONS_00202587 binds to upstream sequences via its secondary structure and interferes with target gene transcription. TCONS_00260893 enhances calcium influx in response to high-temperature treatment by interfering with a specific variant/isoform of the target gene. Heterogeneous expression of these two lncRNA targets promoted photosynthetic protection and recovery, inhibited membrane peroxidation, and suppressed DNA damage in Arabidopsis under heat stress. These results showed that lncRNAs can regulate their target genes by acting as potential RNA scaffolds or through the RNA interference pathway.

Changes in DNA Methylation in Response to 6-Benzylaminopurine Affect Allele-Specific Gene Expression in Populus Tomentosa.

Cytokinins play important roles in the growth and development of plants. Physiological and photosynthetic characteristics are common indicators to measure the growth and development in plants. However, few reports have described the molecular mechanisms of physiological and photosynthetic changes in response to cytokinin, particularly in woody plants. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression in response to the external environment. In this study, we examined genome-wide DNA methylation variation and transcriptional variation in poplar (Populus tomentosa) after short-term treatment with the synthetic cytokinin 6-benzylaminopurine (6-BA). We identified 460 significantly differentially methylated regions (DMRs) in response to 6-BA treatment. Transcriptome analysis showed that 339 protein-coding genes, 262 long non-coding RNAs (lncRNAs), and 15,793 24-nt small interfering RNAs (siRNAs) were differentially expressed under 6-BA treatment. Among these, 79% were differentially expressed between alleles in P. tomentosa, and 102,819 allele-specific expression (ASE) loci in 19,200 genes were detected showing differences in ASE levels after 6-BA treatment. Combined DNA methylation and gene expression analysis demonstrated that DNA methylation plays an important role in regulating allele-specific gene expression. To further investigate the relationship between these 6-BA-responsive genes and phenotypic variation, we performed SNP analysis of 460 6-BA-responsive DMRs via re-sequencing using a natural population of P. tomentosa, and we identified 206 SNPs that were significantly associated with growth and wood properties. Association analysis indicated that 53% of loci with allele-specific expression had primarily dominant effects on poplar traits. Our comprehensive analyses of P. tomentosa DNA methylation and the regulation of allele-specific gene expression suggest that DNA methylation is an important regulator of imbalanced expression between allelic loci.

Indole-3-acetic acid has long-term effects on long non-coding RNA gene methylation and growth in Populus tomentosa.

DNA methylation and long non-coding RNAs (lncRNAs) regulate plant growth and development, but their relationship and effect on responses to the auxin phytohormone indole-3-acetic acid (IAA) remain largely unknown, particularly in woody plants such as poplar (Populus tomentosa). Following treatment of 1-year-old clonal plants with 100 µM IAA, key poplar lncRNA genes showed changes in methylation, but whole-genome methylation levels showed no significant change. Moreover, 100 µM IAA inhibited growth of the 1-year-old poplar clones, possibly through the suppression of photosynthesis. This inhibition had a long-term effect, persisting at 1 month after removal of the exogenous IAA. Transcriptome analysis identified two candidate lncRNA genes that show changes in expression following IAA treatment, TCONS_00003480 and TCONS_00004832. TCONS_00003480 contains the same microRNA target sites of ptc-miR6464 as the 4-coumarate: CoA ligase 2 transcript, which encodes a lignin biosynthesis enzyme. And TCONS_00004832 shares the same target sites of ptc-miR6437a with the Photosystem II reaction center protein D and Cytochrome C Oxidase 17 transcripts, which are related to photosynthesis. The two lncRNAs as the mimics to corresponding target genes of miRNAs to prevent them from degrading. Examination of lncRNA gene expression and methylation revealed a negative relationship (r = - 0.29, P < 0.05); moreover, hypermethylation of the two candidate lncRNA genes remained 1 month after IAA treatment, suggesting that changes in methylation might be involved in the long-term effects of plant hormones. Therefore, our study reveals a long-term effect of IAA on the growth of P. tomentosa, possibly via methylation-mediated epigenetic changes in lncRNA gene expression and the interaction with corresponding miRNAs, leading to regulation of genes related to photosynthesis and growth.

Osmotic stress-responsive promoter upstream transcripts (PROMPTs) act as carriers of MYB transcription factors to induce the expression of target genes in Populus simonii.

Complex RNA transcription and processing produces a diverse range catalog of long noncoding RNAs (lncRNAs), important biological regulators that have been implicated in osmotic stress responses in plants. Promoter upstream transcript (PROMPT) lncRNAs share some regulatory elements with the promoters of their neighbouring protein-coding genes. However, their function remains unknown. Here, using strand-specific RNA sequencing, we identified 209 differentially regulated osmotic-responsive PROMPTs in poplar (Populus simonii). PROMPTs are transcribed bidirectionally and are more stable than other lncRNAs. Co-expression analysis of PROMPTs and protein-coding genes divided the regulatory network into five independent subnetworks including 27 network modules. Significantly enriched PROMPTs in the network were selected to validate their regulatory roles. We used delaminated layered double hydroxide lactate nanosheets (LDH-lactate-NS) to transport synthetic nucleic acids into live tissues to mimic overexpression and interference of a specific PROMPT. The altered expression of PROMPT_1281 induced the expression of its cis and trans targets, and this interaction was governed by its secondary structure rather than just its primary sequence. Based on this example, we proposed a model that a concentration gradient of PROMPT_1281 is established, which increases the probability of its interaction with targets near its transcription site that shares common motifs. Our results firstly demonstrated that PROMPT_1281 act as carriers of MYB transcription factors to induce the expression of target genes under osmotic stress. In sum, our study identified and validated a set of poplar PROMPTs that likely have regulatory functions in osmotic responses.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus.

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

Adaptive evolution and functional innovation of Populus-specific recently evolved microRNAs.

Lineage-specific microRNAs (miRNAs) undergo rapid turnover during evolution; however, their origin and functional importance have remained controversial. Here, we examine the origin, evolution, and potential roles in local adaptation of Populus-specific miRNAs, which originated after the recent salicoid-specific, whole-genome duplication. RNA sequencing was used to generate extensive, comparable miRNA and gene expression data for six tissues. A natural population of Populus trichocarpa and closely related species were used to study the divergence rates, evolution, and adaptive variation of miRNAs. MiRNAs that originated in 5' untranslated regions had higher expression levels and their expression showed high correlation with their host genes. Compared with conserved miRNAs, a significantly higher proportion of Populus-specific miRNAs appear to target genes that were duplicated in salicoids. Examination of single nucleotide polymorphisms in Populus-specific miRNA precursors showed high amounts of population differentiation. We also characterized the newly emerged MIR6445 family, which could trigger the production of phased small interfering RNAs from NAC mRNAs, which encode a transcription factor with primary roles in a variety of plant developmental processes. Together, these observations provide evolutionary insights into the birth and potential roles of Populus-specific miRNAs in genome maintenance, local adaptation, and functional innovation.

Stable methylation of a non-coding RNA gene regulates gene expression in response to abiotic stress in Populus simonii.

DNA methylation plays important roles in responses to environmental stimuli. However, in perennial plants, the roles of DNA methylation in stress-specific adaptions to different abiotic stresses remain unclear. Here, we present a systematic, comparative analysis of the methylome and gene expression in poplar under cold, osmotic, heat, and salt stress conditions from 3h to 24h. Comparison of the stress responses revealed different patterns of cytosine methylation in response to the four abiotic stresses. We isolated and sequenced 1376 stress-specific differentially methylated regions (SDMRs); annotation revealed that these SDMRs represent 1123 genes encoding proteins, 16 miRNA genes, and 17 long non-coding RNA (lncRNA) genes. The SDMR162 region, consisting of Psi-MIR396e and PsiLNCRNA00268512, is regulated by epigenetic pathways and we speculate that PsiLNCRNA00268512 regulates miR396e levels by acting as a target mimic. The ratios of methylated cytosine declined to ~35.1% after 1 month of recovery from abiotic stress and to ~15.3% after 6 months. Among methylated miRNA genes, only expression of the methylation-regulated gene MIRNA6445a showed long-term stability. Our data provide a strong basis for future work and improve our understanding of the effect of epigenetic regulation of non-coding RNA expression, which will enable in-depth functional analysis.

Variation in genomic methylation in natural populations of Populus simonii is associated with leaf shape and photosynthetic traits.

DNA methylation, one of the best-studied types of chromatin modification, suppresses the expression of transposable elements, pseudogenes, repetitive sequences, and individual genes. However, the extent and variation of genome-wide DNA methylation in natural populations of plants remain relatively unknown. To investigate variation in DNA methylation and whether this variation associates with important plant traits, including leaf shape and photosynthesis, 20 413 DNA methylation sites were examined in a poplar association population (505 individuals) using methylation-sensitive amplification polymorphism (MSAP) technology. Calculation of epi-population structure and kinships assigned individuals into subsets (K=3), revealing that the natural population of P. simonii consists of three subpopulations. Population epigenetic distance and geographic distance showed a significant correlation (r=0.4688, P<0.001), suggesting that environmental factors may affect epigenetics. Single-marker approaches were also used to identify significant marker-trait associations, which found 1087 high-confidence DNA methylation markers associated with different phenotypic traits explaining ~5-15% of the phenotypic variance. Among these loci, 147 differentially methylated fragments were obtained by sequencing, representing 130 candidate genes. Expression analysis of six candidate genes indicated that these genes might play important roles in leaf development and regulation of photosynthesis. This study provides association analysis to study the effects of DNA methylation on plant development and these data indicate that epigenetics bridges environmental and genetic factors in affecting plant growth and development.

Population genomic analysis of gibberellin-responsive long non-coding RNAs in Populus.

Long non-coding RNAs (lncRNAs) participate in a wide range of biological processes, but lncRNAs in plants remain largely unknown; in particular, we lack a systematic identification of plant lncRNAs involved in hormone responses. Moreover, allelic variation in lncRNAs remains poorly characterized at a large scale. Here, we conducted high-throughput RNA-sequencing of leaves from control and gibberellin (GA)-treated Populus tomentosa and identified 7655 reliably expressed lncRNAs. Among the 7655 lncRNAs, the levels of 410 lncRNAs changed in response to GA. Seven GA-responsive lncRNAs were predicted to be putative targets of 18 miRNAs, and one GA-responsive lncRNA (TCONS_00264314) was predicted to be a target mimic of ptc-miR6459b. Computational analysis predicted 939 potential cis-regulated target genes and 965 potential trans-regulated target genes for GA-responsive lncRNAs. Functional annotation of these potential target genes showed that they participate in many different biological processes, including auxin signal transduction and synthesis of cellulose and pectin, indicating that GA-responsive lncRNAs may influence growth and wood properties. Finally, single nucleotide polymorphism (SNP)-based association analysis showed that 112 SNPs from 52 GA-responsive lncRNAs and 1014 SNPs from 296 potential target genes were significantly associated with growth and wood properties. Epistasis analysis also provided evidence for interactions between lncRNAs and their potential target genes. Our study provides a comprehensive view of P. tomentosa lncRNAs and offers insights into the potential functions and regulatory interactions of GA-responsive lncRNAs, thus forming the foundation for future functional analysis of GA-responsive lncRNAs in P. tomentosa.

Methylation of microRNA genes regulates gene expression in bisexual flower development in andromonoecious poplar.

Previous studies showed sex-specific DNA methylation and expression of candidate genes in bisexual flowers of andromonoecious poplar, but the regulatory relationship between methylation and microRNAs (miRNAs) remains unclear. To investigate whether the methylation of miRNA genes regulates gene expression in bisexual flower development, the methylome, microRNA, and transcriptome were examined in female and male flowers of andromonoecious poplar. 27 636 methylated coding genes and 113 methylated miRNA genes were identified. In the coding genes, 64.5% of the methylated reads mapped to the gene body region; by contrast, 60.7% of methylated reads in miRNA genes mainly mapped in the 5' and 3' flanking regions. CHH methylation showed the highest methylation levels and CHG showed the lowest methylation levels. Correlation analysis showed a significant, negative, strand-specific correlation of methylation and miRNA gene expression (r=0.79, P <0.05). The methylated miRNA genes included eight long miRNAs (lmiRNAs) of 24 nucleotides and 11 miRNAs related to flower development. miRNA172b might play an important role in the regulation of bisexual flower development-related gene expression in andromonoecious poplar, via modification of methylation. Gynomonoecious, female, and male poplars were used to validate the methylation patterns of the miRNA172b gene, implying that hyper-methylation in andromonoecious and gynomonoecious poplar might function as an important regulator in bisexual flower development. Our data provide a useful resource for the study of flower development in poplar and improve our understanding of the effect of epigenetic regulation on genes other than protein-coding genes.