A Global View of Transcriptome Dynamics During Male Floral Bud Development in Populus tomentosa

Morphology, sucrose metabolism and gene network reveal the molecular mechanism of seed fiber development in poplar

Poplar is an important tree species for ecological protection, wood production, bioenergy and urban greening; it has been widely planted worldwide. However, the catkin fibers produced by female poplars can cause environmental pollution and safety hazards during spring. This study focused on Populus tomentosa, and revealed the sucrose metabolism regulatory mechanism of catkin fibers development from morphological, physiological and molecular aspects. Paraffin section suggested that poplar catkin fibers were not seed hairs and produced from the epidermal cells of funicle and placenta. Sucrose degradation via invertase and sucrose synthase played the dominant role during poplar catkin fibers development. The expression patterns revealed that sucrose metabolism-related genes played important roles during catkin fibers development. Y1H analysis indicated that there was a potential interaction between sucrose synthase 2 (PtoSUS2)/vacuolar invertase 3 (PtoVIN3) and trichome-regulating MYB transcription factors in poplar. Finally, the two key genes, PtoSUS2 and PtoVIN3, had roles in Arabidopsis trichome density, indicating that sucrose metabolism is important in poplar catkin fibers development. This study is not only helpful for clarifying the mechanism of sucrose regulation during trichome development in perennial woody plants, but also establishes a foundation to solve poplar catkin fibers pollution through genetic engineering methods.

Poplar coma morphogenesis and miRNA regulatory networks by combining ovary tissue sectioning and deep sequencing

Poplar coma, commonly referred to as "seed hairs", is a tuft of trichomes attached to the seed coat that helps seed dispersal. However, they can also trigger health impacts for humans, including sneezing, shortness of breath, and skin irritation. Despite efforts to study the regulatory mechanism of herbaceous trichome formation, poplar coma remains poorly understood. In this study, we showed that the epidermal cells of the funiculus and placenta are the origin of poplar coma based on observations of paraffin sections. Small RNA (sRNA) and degradome libraries were also constructed at three stages of poplar coma development, including initiation and elongation stages. Based on 7,904 miRNA-target pairs identified by small RNA and degradome sequencing, we constructed a miRNA-transcript factor and a stage-specific miRNA regulatory network. By combining paraffin section observation and deep sequencing, our research will provide greater insight into the molecular mechanisms of poplar coma development.

Verification of the Protective Effects of Poplar Phenolic Compounds Against Poplar Anthracnose

Poplar anthracnose caused by Colletotrichum gloeosporioides is one of the most important diseases widely distributed in poplar-growing areas in China, causing serious economic and ecological losses. In this study, three poplar species showed different resistance to poplar anthracnose: Populus × canadensis was resistant, Populus tomentosa was susceptible, and P. × beijingensis showed intermediate resistance. However, it remains uncertain whether phenolic compounds in poplar are involved in this resistance. Therefore, we determined the concentrations of phenolic compounds and their antifungal activity. Before and after the C. gloeosporioides inoculation, 20 phenolic compounds were detected in P. × canadensis and the number increased from 12 to 14 in P. × beijingensis but decreased from seven to four in P. tomentosa. Thus, phenolic compounds may be positively correlated with the degree of disease resistance. We selected seven phenolic compounds for further analysis, which varied considerably in content after inoculation with C. gloeosporioides. These seven compounds were salicin, arbutin, benzoic acid, salicylic acid, chlorogenic acid, ferulic acid, and naringenin, which helped poplar trees to limit the growth of C. gloeosporioides and differed in their antifungal effects, with phenolic acids having the strongest inhibitory effect. In addition, the optimal concentrations of different substances varied. We demonstrate that these phenolic compounds produced by poplar do play a certain role in the process of poplar resistance to anthracnose. These findings lay a foundation for future research into the antifungal mechanism of poplar trees and may be useful for enhancing the prevention and control of poplar anthracnose.

Genome-wide analysis and molecular dissection of the SPL gene family in Fraxinus mandshurica

BACKGROUND

SQUAMOSA promoter binding protein-like (SPL) is a unique family of transcription factors in plants, which is engaged in regulating plant growth and development, physiological and biochemical processes. Fraxinus mandshurica is an excellent timber species with a wide range of uses in northeastern China and enjoys a high reputation in the international market. SPL family analysis has been reported in some plants while SPL family analysis of Fraxinus mandshurica has not been reported.

RESULTS

We used phylogeny, conserved motifs, gene structure, secondary structure prediction, miR156 binding sites, promoter cis elements and GO annotation to systematically analyze the FmSPLs family. This was followed by expression analysis by subcellular localization, expression patterns at various tissue sites, abiotic stress and hormone induction. Because FmSPL2 is highly expressed in flowers it was selected to describe the SPL gene family of Fraxinus mandshurica by ectopic expression. Among them, 10 FmSPL genes that were highly expressed at different loci were selected for expression analysis under abiotic stress (NaCl and Cold) and hormone induction (IAA and ABA). These 10 FmSPL genes showed corresponding trends in response to both abiotic stress and hormone induction. We showed that overexpression of FmSPL2 in transgenic Nicotiana tabacum L. resulted in taller plants, shorter root length, increased root number, rounded leaves, and earlier flowering time.

CONCLUSIONS

We identified 36 SPL genes, which were classified into seven subfamilies based on sequence analysis. FmSPL2 was selected for subsequent heterologous expression by analysis of expression patterns in various tissues and under abiotic stress and hormone induction, and significant phenotypic changes were observed in the transgenic Nicotiana tabacum L. These results provide insight into the evolutionary origin and biological significance of plant SPL. The aim of this study was to lay the foundation for the genetic improvement of Fraxinus mandshurica and the subsequent functional analysis of FmSPL2.

Time-series RNA-Seq transcriptome profiling reveals novel insights about cold acclimation and de-acclimation processes in an evergreen shrub of high altitude

The high-altitude alpine regions are characterized by highly variable and harsh environmental conditions. However, relatively little is known about the diverse mechanisms adopted by alpine plants to adapt to these stressful conditions. Here, we studied variation in transcriptome and physiological adjustments occurring across the year at high elevation environments in the leaf tissue of Rhododendron anthopogon, an evergreen shrub of Himalaya. The samples were collected at 12 different time-points, from August until snowfall in November 2017, and then from June to September 2018. It was observed that with a drop in both ambient air temperature and photoperiod towards onset of winter, the freezing resistance of plants increased, resulting in 'cold acclimation'. Further, 'de-acclimation' was associated with a decrease in freezing resistance and increase in photosynthetic efficiency of leaves during spring. A considerable amount of variation was observed in the transcriptome in a time-dependent sequential manner, with a total of 9,881 differentially expressed genes. Based on gene expression profiles, the time-points could be segregated into four clusters directly correlating with the distinct phases of acclimation: non-acclimation (22-August-2017, 14-August-2018, 31-August-2018), early cold acclimation (12-September-2017, 29-September-2017), late cold acclimation (11-October-2017, 23-October-2017, 04-November-2017, 18-September-2018) and de-acclimation (15-June-2018, 28-June-2018, 14-July-2018). Cold acclimation was a gradual process, as indicated by presence of an intermediate stage (early acclimation). However, the plants can by-pass this stage when sudden decrease in temperature is encountered. The maximum variation in expression levels of genes occurred during the transition to de-acclimation, hence was 'transcriptionally' the most active phase. The similar or higher expression levels of genes during de-acclimation in comparison to non-acclimation suggested that molecular functionality is re-initiated after passing through the harsh winter conditions.

ARR17 controls dioecy in Populus by repressing B-class MADS-box gene expression

The number of dioecious species for which the genetic basis of sex determination has been resolved is rapidly increasing. Nevertheless, the molecular mechanisms downstream of the sex determinants remain largely elusive. Here, by RNA-sequencing early-flowering isogenic aspen (Populus tremula) lines differing exclusively for the sex switch gene ARR17, we show that a narrowly defined genetic network controls differential development of female and male flowers. Although ARR17 encodes a type-A response regulator supposedly involved in cytokinin (CK) hormone signalling, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated arr17 knockout only affected the expression of a strikingly small number of genes, indicating a specific role in the regulation of floral development rather than a generic function in hormone signalling. Notably, the UNUSUAL FLORAL ORGANS (UFO) gene, encoding an F-box protein acting as a transcriptional cofactor with LEAFY (LFY) to activate B-class MADS-box gene expression, and the B-class gene PISTILLATA (PI), necessary for male floral organ development, were strongly de-repressed in the arr17 CRISPR mutants. Our data highlight a CK-independent role of the poplar response regulator ARR17 and further emphasize the minimal differences between female and male individuals. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Comprehensive Analyses of Four PtoNF-YC Genes from Populus tomentosa and Impacts on Flowering Timing

Flowering is an important link in the life process of angiosperms, and it is also an important sign of the transformation of plants from vegetative to reproductive growth. Although the flowering regulation network of Arabidopsis is well-understood, there has been little research on the molecular mechanisms of perennial woody plant flower development regulation. Populus tomentosa is a unique Chinese poplar species with fast growth, strong ecological adaptability, and a long lifecycle. However, it has a long juvenile phase, which seriously affects its breeding process. Nuclear factor-Y (NF-Y) is an important type of transcription factor involved in the regulation of plant flowering. However, there are few reports on PtoNF-Y gene flowering regulation, and the members of the PtNF-YC subfamily are unknown. In this study, four key genes were cloned and analyzed for sequence characteristics, gene structure, genetic evolution, expression patterns, and subcellular localization. The plant expression vector was further constructed, and transgenic Arabidopsis and P. tomentosa plants were obtained through genetic transformation and a series of molecular tests. The flowering time and other growth characteristics were analyzed. Finally, the expression level of flowering genes was detected by quantitative PCR, the interaction between PtoNF-YC and PtoCOL proteins was measured using the yeast two-hybrid system to further explain the flowering regulation mechanism, and the molecular mechanisms by which PtNF-YC6 and PtNF-YC8 regulate poplar flowering were discussed. These results lay the foundation for elucidating the molecular regulation mechanism of PtoNF-YC in flowering and furthering the molecular design and breeding of poplar, while providing a reference for other flowering woody plants.

A MADS-box gene with similarity to FLC is induced by cold and correlated with epigenetic changes to control budbreak in kiwifruit

Temperate perennials require exposure to chilling temperatures to resume growth in the following spring. Growth and dormancy cycles are controlled by complex genetic regulatory networks and are governed by epigenetic mechanisms, but the specific genes and mechanisms remain poorly understood. To understand how seasonal changes and chilling regulate dormancy and growth in the woody perennial vine kiwifruit (Ac, Actinidia chinensis), a transcriptome study of kiwifruit buds in the field and controlled conditions was performed. A MADS-box gene with homology to Arabidopsis FLOWERING LOCUS C (FLC) was identified and characterized. Elevated expression of AcFLC-like (AcFLCL) was detected during bud dormancy and chilling. A long noncoding (lnc) antisense transcript with an expression pattern opposite to AcFLCL and shorter sense noncoding RNAs were identified. Chilling induced an increase in trimethylation of lysine-4 of histone H3 (H3K4me3) in the 5' end of the gene, indicating multiple layers of epigenetic regulation in response to cold. Overexpression of AcFLCL in kiwifruit gave rise to plants with earlier budbreak, whilst gene editing using CRISPR-Cas9 resulted in transgenic lines with substantially delayed budbreak, suggesting a role in activation of growth. These results have implications for the future management and breeding of perennials for resilience to changing climate.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

Comparative Transcriptome Reveals the Genes' Adaption to Herkogamy of Lumnitzera littorea (Jack) Voigt

Lumnitzera littorea (Jack) Voigt is among the most endangered mangrove species in China. The morphology and evolution of L. littorea flowers have received substantial attention for their crucial reproductive functions. However, little is known about the genomic regulation of flower development in L. littorea. In this study, we characterized the morphology of two kinds of L. littorea flowers and performed comparative analyses of transcriptome profiles of the two different flowers. Morphological observation showed that some flowers have a column embedded in the petals while others produce a stretched flower style during petal unfolding in flowering. By using RNA-seq, we obtained 138,857 transcripts that were assembled into 82,833 unigenes with a mean length of 1055.48 bp. 82,834 and 34,997 unigenes were assigned to 52 gene ontology (GO) functional groups and 364 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 4,267 differentially expressed genes (DEGs), including 1,794 transcription factors (TFs), were identified between two types of flowers. These TFs are mainly involved in bHLH, B3, bZIP, MYB-related, and NAC family members. We further validated that 12 MADS-box genes, including 4 MIKC-type and 8 M-type TFs, were associated with the pollinate of L. littorea by herkogamy. Our current results provide valuable information for genetic analysis of L. littorea flowering and may be useful for illuminating its adaptive evolutionary mechanisms.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

EST-SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives

BACKGROUND

The species Olea europaea includes cultivated varieties (subsp. europaea var. europaea), wild plants (subsp. europaea var. sylvestris), and five other subspecies spread over almost all continents. Single nucleotide polymorphisms in the expressed sequence tag able to underline intra-species differentiation are not yet identified, beyond a few plastidial markers.

METHODS

In the present work, more than 1000 transcript-specific SNP markers obtained by the genotyping of 260 individuals were studied. These genotypes included cultivated, oleasters, and samples of subspecies guanchica, and were analyzed in silico, in order to identify polymorphisms on key genes distinguishing different Olea europaea forms.

RESULTS

Phylogeny inference and principal coordinate analysis allowed to detect two distinct clusters, clearly separating wilds and guanchica samples from cultivated olives, meanwhile the structure analysis made possible to differentiate these three groups. Sequences carrying the polymorphisms that distinguished wild and cultivated olives were analyzed and annotated, allowing to identify 124 candidate genes that have a functional role in flower development, stress response, or involvement in important metabolic pathways. Signatures of selection that occurred during olive domestication, were detected and reported.

CONCLUSION

This deep EST-SNP analysis provided important information on the genetic and genomic diversity of the olive complex, opening new opportunities to detect gene polymorphisms with potential functional and evolutionary roles, and to apply them in genomics-assisted breeding, highlighting the importance of olive germplasm conservation.

Gene expression trajectories during male and female reproductive development in balsam poplar (Populus balsamifera L.)

Plant reproductive development from the first appearance of reproductively committed axes through to floral maturation requires massive and rapid remarshalling of gene expression. In dioecious species such as poplar this is further complicated by divergent male and female developmental programs. We used seven time points in male and female balsam poplar (Populus balsamifera L.) buds and catkins representing the full annual flowering cycle, to elucidate the effects of time and sex on gene expression during reproductive development. Time (developmental stage) is dominant in patterning gene expression with the effect of sex nested within this. Here, we find (1) evidence for five successive waves of alterations to the chromatin landscape which may be important in setting the overall reproductive trajectory, regardless of sex. (2) Each individual developmental stage is further characterized by marked sex-differential gene expression. (3) Consistent sexually differentiated gene expression regardless of developmental stage reveal candidates for high-level regulators of sex and include the female-specific poplar ARR17 homologue. There is also consistent male-biased expression of the MADS-box genes PISTILLATA and APETALA3. Our work provides insights into expression trajectories shaping reproductive development, its potential underlying mechanisms, and sex-specific translation of the genome information into reproductive structures in balsam poplar.

Altered sucrose metabolism and plant growth in transgenic Populus tomentosa with altered sucrose synthase PtSS3

Improvement of wood quality is an important focus of forest genetics and breeding research. Sucrose synthase (SS) catalyzes the reaction of sucrose and uridine diphosphate into uridine diphosphate glucose and fructose. It is a key enzyme involved in cell wall formation during secondary growth by providing the UDP-Glucose substrate for cellulose biosynthesis. In this study, we isolated the single-copy gene PtSS3 from the SS gene family of Populus tomentosa and analyzed its structure. To identify its function in secondary growth, we generated 19 transgenic lines of P. tomentosa using PtSS3 overexpression (OE) and artificial microRNA (amiRNA) constructs. We also performed comprehensive analyses of the transgenic P. tomentosa plants, including phenotypic analyses, quantitative real-time PCR, enzyme activity assays and sugar metabolism. We found significantly higher PtSS3 enzyme activity, fructose, and glucose levels and significantly lower sucrose levels in the stems and leaves of OE-PtSS3 plants. The opposite trend was observed in the amiRNA-PtSS3 lines. Gene expression analyses showed that PtSS3 transcript levels in stems and leaves were up-regulated in the OE-PtSS3 lines and down-regulated in the amiRNA-PtSS3 lines, and the OE-PtSS3 plants grew taller than the wild-type and amiRNA-PtSS3 plants. These findings indicate that PtSS3 plays an important role in sucrose metabolism and growth of trees.

Generation of a Comprehensive Transcriptome Atlas and Transcriptome Dynamics in Medicinal Cannabis

Cannabinoids are the main medicinal compounds of interest in the plant Cannabis sativa, that are primarily synthesised in the glandular trichomes; found on female floral buds. The content, composition and yield of secondary metabolites (cannabinoids and terpenoids) is influenced by the plant's genetics and environment. Some initial gene expression experiments have been performed from strains of this plant species that contrasted in cannabinoid production, however the present knowledge about detailed trichome transcriptomics in this species is limited. An extensive transcriptome atlas was generated by RNA sequencing using root, shoot, flower and trichome tissues from a female plant strain (Cannbio-2) and was enhanced with the addition of vegetative and reproductive tissues from a male cannabis plant. Differential gene expression analysis identified genes preferentially expressed in different tissues. Detailed trichomics was performed from extractions specifically from glandular trichomes as well as female floral tissues at varying developmental stages, to identify stage-specific differentially expressed genes. Candidate genes involved in terpene and cannabinoid synthesis were identified and the majority were found to have an abundant expression in trichomes. The comprehensive transcriptome is a significant resource in cannabis for further research of functional genomics to improve the yield of specialised metabolites with high pharmacological value.

Constitutive Expression of an Apple FLC3-like Gene Promotes Flowering in Transgenic Blueberry under Nonchilling Conditions

MADS-box transcription factors FLOWERING LOCUS C (FLC) and APETALA1 (AP1)/CAULIFLOWER (CAL) have an opposite effect in vernalization-regulated flowering in Arabidopsis. In woody plants, a functional FLC-like gene has not been verified through reverse genetics. To reveal chilling-regulated flowering mechanisms in woody fruit crops, we conducted phylogenetic analysis of the annotated FLC-like proteins of apple and found that these proteins are grouped more closely to Arabidopsis AP1 than the FLC group. An FLC3-like MADS-box gene from columnar apple trees (Malus domestica) (MdFLC3-like) was cloned for functional analysis through a constitutive transgenic expression. The MdFLC3-like shows 88% identity to pear's FLC-like genes and 82% identity to blueberry's CAL1 gene (VcCAL1). When constitutively expressed in a highbush blueberry (Vaccinium corymbosum L.) cultivar 'Legacy', the MdFLC3-like induced expressions of orthologues of three MADS-box genes, including APETALA1, SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1, and CAL1. As a consequence, in contrast to the anticipated late flowering associated with an overexpressed FLC-like, the MdFLC3-like promoted flowering of transgenic blueberry plants under nonchilling conditions where nontransgenic 'Legacy' plants could not flower. Thus, the constitutively expressed MdFLC3-like in transgenic blueberries functioned likely as a blueberry's VcCAL1. The results are anticipated to facilitate future studies for revealing chilling-mediated flowering mechanisms in woody plants.