Deqiang Zhang

The sucrose synthase gene family in Populus: structure, expression, and evolution

Sucrose synthase is a key enzyme in sucrose metabolism in plant cells, and it is involved in the synthesis of cell wall cellulose. Although the sucrose synthase gene (SUS) family in the model plants Arabidopsis thaliana has been characterized, little is known about this gene family in trees. This study reports the identification of two novel SUS genes in the economically important poplar tree. These genes were expressed predominantly in mature xylem. Using molecular cloning and bioinformatics analysis of the Populus genome, we demonstrated that SUS is a multigene family with seven members that each exhibit distinct but partially overlapping expression patterns. Of particular interest, three SUS genes were preferentially expressed in the stem xylem, suggesting that poplar SUSs are involved in the formation of the secondary cell wall. Gene structural and phylogenetic analyses revealed that the Populus SUS family is composed of four main subgroups that arose before the separation of monocots and dicots. Phylogenetic analyses associated with the tissue- and organ-specific expression patterns. The high intraspecific nucleotide diversity of two SUS genes was detected in the natural population, and the πnonsyn/πsyn ratio was significantly less than 1; therefore, SUS genes appear to be evolving in Populus, primarily under purifying selection. This is the first comprehensive study of the SUS gene family in woody plants; the analysis includes genome organization, gene structure, and phylogeny across land plant lineages, as well as expression profiling in Populus.

Effects of high temperature on photosynthesis and related gene expression in poplar

BackgroundHigh temperature, whether transitory or constant, causes physiological, biochemical and molecular changes that adversely affect tree growth and productivity by reducing photosynthesis. To elucidate the photosynthetic adaption response and examine the recovery capacity of trees under heat stress, we measured gas exchange, chlorophyll fluorescence, electron transport, water use efficiency, and reactive oxygen-producing enzyme activities in heat-stressed plants.ResultsWe found that photosynthesis could completely recover after less than six hours of high temperature treatment, which might be a turning point in the photosynthetic response to heat stress. Genome-wide gene expression analysis at six hours of heat stress identified 29,896 differentially expressed genes (15,670 up-regulated and 14,226 down-regulated), including multiple classes of transcription factors. These interact with each other and regulate the expression of photosynthesis-related genes in response to heat stress, controlling carbon fixation and changes in stomatal conductance. Heat stress of more than twelve hours caused reduced electron transport, damaged photosystems, activated the glycolate pathway and caused H2O2 production; as a result, photosynthetic capacity did not recover completely.ConclusionsThis study provides a systematic physiological and global gene expression profile of the poplar photosynthetic response to heat stress and identifies the main limitations and threshold of photosynthesis under heat stress. It will expand our understanding of plant thermostability and provides a robust dataset for future studies.

Genome-wide identification of novel long non-coding RNAs in Populus tomentosa tension wood, opposite wood and normal wood xylem by RNA-seq

Increasing evidence shows that long non-coding RNAs (lncRNAs) function as important regulatory factors during plant development, but few reports have examined lncRNAs in trees. Here, we report our genome-scale identification and characterization of lncRNAs differentially expressed in the xylem of tension wood, opposite wood and normal wood in Populustomentosa, by high-throughput RNA sequencing. We identified 1,377 putative lncRNAs by computational analysis, and expression and structure analyses showed that the lncRNAs had lower expression levels and shorter lengths than protein-coding transcripts in Populus. Of the 776 differently expressed (log2FC ≥1 or ≤-1, FDR ≤0.01) lncRNAs, 389 could potentially target 1,151 genes via trans-regulatory effects. Functional annotation of these target genes demonstrated that they are involved in fundamental processes, and in specific mechanisms such as response to stimuli. We also identified 16 target genes involved in wood formation, including cellulose and lignin biosynthesis, suggesting a potential role for lncRNAs in wood formation. In addition, three lncRNAs harbor precursors of four miRNAs, and 25 were potentially targeted by 44 miRNAs where a negative expression relationship between them was detected by qRT-PCR. Thus, a network of interactions among the lncRNAs, miRNAs and mRNAs was constructed, indicating widespread regulatory interactions between non-coding RNAs and mRNAs. Lastly, qRT-PCR validation confirmed the differential expression of these lncRNAs, and revealed that they have tissue-specific expression in P. tomentosa. This study presents the first global identification of lncRNAs and their potential functions in wood formation, providing a starting point for detailed dissection of the functions of lncRNAs in Populus.

Exploring the Secrets of Long Noncoding RNAs

High-throughput sequencing has revealed that the majority of RNAs have no capacity to encode protein. Among these non-coding transcripts, recent work has focused on the roles of long noncoding RNAs (lncRNAs) of >200 nucleotides. Although many of their attributes, such as patterns of expression, remain largely unknown, lncRNAs have key functions in transcriptional, post-transcriptional, and epigenetic gene regulation; Also, new work indicates their functions in scaffolding ribonuclear protein complexes. In plants, genome-wide identification of lncRNAs has been conducted in several species, including Zea mays, and recent research showed that lncRNAs regulate flowering time in the photoperiod pathway, and function in nodulation. In this review, we discuss the basic mechanisms by which lncRNAs regulate key cellular processes, using the large body of knowledge on animal and yeast lncRNAs to illustrate the significance of emerging work on lncRNAs in plants.

Genome-Wide Analysis of Gene Expression in Response to Drought Stress in Populus simonii

Drought tolerance varies considerably in different species of Populus, and is a complex trait, involving the interplay of a vast array of genes. Although Populus simonii is one of the most important commercial plantation tree species in China, genes controlling drought-stress tolerance in this organism have not yet been identified. Here, transcriptomic changes during drought stress in P. simonii were detected using an Affymetrix GeneChip. In total, 1,028 transcripts were identified as differentially expressed under drought stress, of which 496 transcripts increased and 532 decreased in abundance (two-way ANOVA, fold change >4 and P < 0.01). Expression changes of 20 candidate genes were validated by real-time PCR, indicating significant differences between the controls and water-deficit treatments. Gene annotation demonstrated that the majority of these P. simonii genes encode products involved in phytohormone metabolism, osmoregulation, and oxidative stress. Based on gene ontology (GO) classification, the increased genes were classified into six significantly enriched GO terms involved in 64 pathways, and the decreased genes were classified into 39 significantly enriched GO terms representing 42 pathways. Bioinformatics analysis showed that the differential expression of plant hormone-related, transcription factor, cytochrome P450 gene superfamily, carbohydrate metabolism, and amino acid transporter genes may contribute to drought-stress tolerance in P. simonii. Our study provides global gene expression patterns during drought stress and will be valuable for further study of the molecular mechanisms of drought tolerance in P. simonii.

Polymorphic simple sequence repeat (SSR) loci within cellulose synthase (PtoCesA) genes are associated with growth and wood properties in Populus tomentosa.

Chinese white poplar (Populus tomentosa), an important commercial tree species for timber and pulp production in northern China, has been used to examine the individual genes and allelic diversity responsible for complex traits controlling growth and lignocellulosic biosynthesis. Taking advantage of the low degree of linkage disequilibrium (LD) within P. tomentosa association populations, we examined associations between 15 cellulose synthase (PtoCesA) genes and traits including growth and wood properties. Thirty-six novel simple sequence repeat (SSR) markers within PtoCesA genes were detected by re-sequencing and genotyped in an association population (460 individuals). Single-marker and haplotype-based LD approaches were used to identify significant marker-trait associations. Family-based linkage studies and real-time PCR testing were conducted to validate the functional significance of SSR variation. Fifteen single-marker associations from seven PtoCesA genes and nine haplotype-based associations within six genes were identified in the association population (false discovery rate Q < 0.05). Next, five SSR marker-trait associations (Q < 0.05) from four PtoCesA genes were successfully validated in a linkage mapping population (1200 individuals). The results imply a functional role for these genes in mediating wood properties, demonstrating the potential of combining single-marker and haplotype-based LD approaches to detect functional allelic variation underlying quantitative traits in a low-LD population.

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.

Sexual Dimorphism Floral MicroRNA Profiling and Target Gene Expression in Andromonoecious Poplar (Populus tomentosa)

Although the molecular basis of poplar sex-specific flower development remains largely unknown, increasing evidence indicates an essential role for microRNAs (miRNAs). The specific miRNA types and precise miRNA expression patterns in dioecious plant flower development remain unclear. Here, we used andromonoecious poplar, an exceptional model system, to eliminate the confounding effects of genetic background of dioecious plants. This system, combined with high-throughput sequencing and computational analysis, allowed us to characterize sex-specific miRNAomes from female and male flowers. Comparative miRNAome analysis combined with quantitative real-time PCR revealed the expression patterns of 27 miRNAs in poplar flower and showed that the targets of these miRNAs are involved in flower organogenesis, Ca2+ transport, phytohormone synthesis and metabolism, and DNA methylation. This paper describes a complex regulatory network consisting of these miRNAs expressed in sex-specific flower development in a dioecious plant. The conserved and novel miRNA locations were annotated in the Populus trichocarpa genome. Among these, miRNA Pto-F70 and 4 targets are located in the sex-determination regions of chromosome XIX. Furthermore, two novel miRNAs, Pto-F47 and Pto-F68, were shown for the first time to be regulatory factors in phytohormone interactions. To our knowledge, this report is the first systematic investigation of sex-specific flower-related miRNAs and their targets in poplar, and it deepens our understanding of the important regulatory functions of miRNAs in female and male flower development in this dioecious plant.

Allelic Variation in PtGA20Ox Associates with Growth and Wood Properties in Populus spp

Populus tomentosa is an economically important tree crop that produces wood for lumber, pulp, paper, and biofuels. Wood quality traits are likely to be strongly affected by the plant hormone gibberellic acid (GA), which regulates growth. GA20Ox encodes one of the major regulatory enzymes of GA biosynthesis and may therefore play a large role in growth and wood quality. Here, linkage disequilibrium (LD) studies were used to identify significant associations between single nucleotide polymorphisms (SNPs) within PtGA20Ox and growth and wood-quality traits of P. tomentosa. We isolated a full-length GA20Ox cDNA from Populus tomentosa by reverse transcription (RT)-PCR; this 1401 bp cDNA clone had an open reading frame of 1158 bp and encoded a protein of 385 amino acids. PtGA20Ox transcripts were maximally expressed in the mature xylem of vascular tissues, suggesting that PtGA20Ox is highly expressed and specifically associated with secondary xylem formation. Resequencing the PtGA20Ox locus of 36 individuals identified 55 SNPs, and the frequency of SNPs was 1/31 bp. The 29 most common SNPs (frequency>0.1) were genotyped in an association population (426 individuals) that was also phenotyped for key growth and wood quality traits. LD did not extend over the entire gene (r 2<0.1, within 500 bp), demonstrating that a candidate-gene-based LD approach may the best way to understand the molecular basis underlying quantitative variation in this species. SNP- and haplotype-based association analyses indicated that four SNPs (false discovery rate Q<0.05) and 14 haplotypes (P<0.05) were significantly associated with growth and wood properties. The phenotypic variance explained by each SNP ranged from 3.44% to 14.47%. The SNP markers identified in this study can be applied to breeding programs for the improvement of growth and wood-property traits by marker-assisted selection.

Genetic architecture of growth traits in Populus revealed by integrated quantitative trait locus (QTL) analysis and association studies

Deciphering the genetic architecture underlying polygenic traits in perennial species can inform molecular marker-assisted breeding. Recent advances in high-throughput sequencing have enabled strategies that integrate linkage-linkage disequilibrium (LD) mapping in Populus. We used an integrated method of quantitative trait locus (QTL) dissection with a high-resolution linkage map and multi-gene association mapping to decipher the nature of genetic architecture (additive, dominant, and epistatic effects) of potential QTLs for growth traits in a Populus linkage population (1200 progeny) and a natural population (435 individuals). Seventeen QTLs for tree height, diameter at breast height, and stem volume mapped to 11 linkage groups (logarithm of odds (LOD) ≥ 2.5), and explained 2.7-18.5% of the phenotypic variance. After comparative mapping and transcriptome analysis, 187 expressed genes (10 046 common single nucleotide polymorphisms (SNPs)) were selected from the segmental homology regions (SHRs) of 13 QTLs. Using multi-gene association models, we observed 202 significant SNPs in 63 promising genes from 10 QTLs (P ≤ 0.0001; FDR ≤ 0.10) that exhibited reproducible associations with additive/dominant effects, and further determined 11 top-ranked genes tightly linked to the QTLs. Epistasis analysis uncovered a uniquely interconnected gene-gene network for each trait. This study opens up opportunities to uncover the causal networks of interacting genes in plants using an integrated linkage-LD mapping approach.