A transcriptomics investigation into pine reproductive organ development

Deciphering the Role of SVP-Like Genes and Their Key Regulation Networks During Reproductive Cone Development in Pinus tabuliformis

Reproductive development plays an essential role in the perpetuation of genetic material and environmental adaptation. In angiosperms, the Short Vegetative Phase (SVP) serves as a flowering repressor, influencing the development of floral organs. In this study, heterologous transformation of Arabidopsis thaliana with SVP-like genes (PtSVL1 and PtSVL2) derived from Pinus tabuliformis significantly impacted stamen formation and pollen fertility, without altering flowering time. Gene co-expression networks revealed that SVP-like and SOC1-like genes function as key coregulatory transcription factors during the initial stages of cone development in P. tabuliformis. Interestingly, the regulatory module of SOC1 regulated by SVP in angiosperms is absent in conifers and conifer SVP-like exercises its function in a form that is physically bound to SOC1-like. Furthermore, combining the yeast one-hybrid scanning with co-expression network analysis, revealed that SPLs and TPSs were the principal downstream target genes of PtSVL1. Notably, the PtSPL16 promoter is positively regulated by PtSVL1, and overexpression of PtSPL16 results in delayed flowering in Arabidopsis, suggesting that the PtSVL1-PtSPL16 module plays a crucial role in regulating reproductive development in conifers. Collectively, these findings enhance our understanding of the roles of SVP-like genes in conifers and the key regulatory networks centred on PtSVL1 during reproductive cone development.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Phytohormone biosynthesis and transcriptional analyses provide insight into the main growth stage of male and female cones Pinus koraiensis

The cone is a crucial component of the whole life cycle of gymnosperm and an organ for sexual reproduction of gymnosperms. In Pinus koraiensis, the quantity and development process of male and female cones directly influence seed production, which in turn influences the tree's economic value. There are, however, due to the lack of genetic information and genomic data, the morphological development and molecular mechanism of female and male cones of P. koraiensis have not been analyzed. Long-term phenological observations were used in this study to document the main process of the growth of both male and female cones. Transcriptome sequencing and endogenous hormone levels at three critical developmental stages were then analyzed to identify the regulatory networks that control these stages of cones development. The most significant plant hormones influencing male and female cones growth were discovered to be gibberellin and brassinosteroids, according to measurements of endogenous hormone content. Additionally, transcriptome sequencing allowed the identification of 71,097 and 31,195 DEGs in male and female cones. The synthesis and control of plant hormones during cones growth were discovered via enrichment analysis of key enrichment pathways. FT and other flowering-related genes were discovered in the coexpression network of flower growth development, which contributed to the growth development of male and female cones of P. koraiensis. The findings of this work offer a cutting-edge foundation for understanding reproductive biology and the molecular mechanisms that control the growth development of male and female cones in P. koraiensis.

Model Species to Investigate the Origin of Flowers

The angiosperms, or flowering plants, arose at least 135 million years ago (Ma) and rapidly diversified to form over 300,000 species alive today. This group appears, however, to have separated from its closest living relatives, the extant gymnosperms, much earlier: over 300 Ma. Representatives of basally-diverging angiosperm lineages are of key importance to studies aimed at reconstructing the most recent common ancestor of living angiosperms, including its morphological, anatomical, eco-physiological and molecular aspects. Furthermore, evo-devo comparisons of angiosperms with living gymnosperms may help to determine how the many novel aspects of angiosperms, including those of the flower, first came about. This chapter reviews literature on the origin of angiosperms and focusses on basally-diverging angiosperms and gymnosperms that show advantages as potential experimental models, reviewing information and protocols for the use of these species in an evo-devo context. The final section suggests a means by which data from living and fossil groups could be integrated to better elucidate evolutionary events that took place on the long stem-lineage that apparently preceded the radiation of living angiosperms.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

Involvement of PtCOL5-PtNF-YC4 in reproductive cone development and gibberellin signaling in Chinese pine

It is well documented that the CO/NF-YB/NF-YC trimer (NF-Y-CO) binds and regulates the FT promoter. However, the FT/TFL1-like (FLOWERING LOCUS T/TERMINALFLOWER1-like) genes in gymnosperms are all flowering suppressors, and the regulation model of NF-Y in gymnosperms is different from that in angiosperms. Here, using Chinese pine (Pinus tabuliformis), we identified a CONSTANS-LIKE gene, PtCOL5, the expression of which was strongly induced during cones development and it functioned as a repressor of flowering. PtNF-YC4, which interacted with PtCOL5, was highly correlated with PtCOL5 during growth and development, has been demonstrated. Moreover, PtNF-YC4 and PtCOL5 can bind to PtTFL2 promoter, and their interaction can enhance PtTFL2 expression. Interestingly, we found PtNF-YC4 and PtCOL5 were involved in gibberellin signaling and their interaction was inhibited by PtDELLA protein, thus affecting PtTFL2 expression. Collectively, PtCOL5-PtNF-YC4 was involved in reproductive cone development and gibberellin signaling in Chinese pine. Our findings uncovered reproductive cone development and signal transduction mechanism of COL-NF-Y in gymnosperms.

Morpho-histology, endogenous hormone dynamics, and transcriptome profiling in Dacrydium pectinatum during female cone development

Dacrydium pectinatum de Laubenfels is a perennial dioeciously gymnosperm species dominant in tropical montane rain forests. Due to deforestation, natural disasters, long infancy, and poor natural regeneration ability, the population of this species has been significantly reduced and listed as an endangered protected plant. To better understand the female cone development in D. pectinatum, we examined the morphological and anatomical changes, analyzed the endogenous hormone dynamics, and profiled gene expression. The female reproductive structures were first observed in January. The morpho-histological observations suggest that the development of the D. pectinatum megaspore can be largely divided into six stages: early flower bud differentiation, bract primordium differentiation, ovule primordium differentiation, dormancy, ovule maturity, and seed maturity. The levels of gibberellins (GA), auxin (IAA), abscisic acid (ABA), and cytokinin (CTK) fluctuate during the process of female cone development. The female cones of D. pectinatum need to maintain a low level of GA3-IAA-ABA steady state to promote seed germination. The first transcriptome database for female D. pectinatum was generated, revealing 310,621 unigenes. Differential expression analyses revealed several floral (MADS2, AGL62, and LFY) and hormone biosynthesis and signal transduction (CKX, KO, KAO, ABA4, ACO, etc.) genes that could be critical for female cone development. Our study provides new insights into the cone development in D. pectinatum and the foundation for female cone induction with hormones.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

MADS-box transcription factors MADS11 and DAL1 interact to mediate the vegetative-to-reproductive transition in pine

The reproductive transition is an important event that is crucial for plant survival and reproduction. Relative to the thorough understanding of the vegetative phase transition in angiosperms, a little is known about this process in perennial conifers. To gain insight into the molecular basis of the regulatory mechanism in conifers, we used temporal dynamic transcriptome analysis with samples from seven different ages of Pinus tabuliformis to identify a gene module substantially associated with aging. The results first demonstrated that the phase change in P. tabuliformis occurred as an unexpectedly rapid transition rather than a slow, gradual progression. The age-related gene module contains 33 transcription factors and was enriched in genes that belong to the MADS (MCMl, AGAMOUS, DEFICIENS, SRF)-box family, including six SOC1-like genes and DAL1 and DAL10. Expression analysis in P. tabuliformis and a late-cone-setting P. bungeana mutant showed a tight association between PtMADS11 and reproductive competence. We then confirmed that MADS11 and DAL1 coordinate the aging pathway through physical interaction. Overexpression of PtMADS11 and PtDAL1 partially rescued the flowering of 35S::miR156A and spl1,2,3,4,5,6 mutants in Arabidopsis (Arabidopsis thaliana), but only PtMADS11 could rescue the flowering of the ft-10 mutant, suggesting PtMADS11 and PtDAL1 play different roles in flowering regulatory networks in Arabidopsis. The PtMADS11 could not alter the flowering phenotype of soc1-1-2, indicating it may function differently from AtSOC1 in Arabidopsis. In this study, we identified the MADS11 gene in pine as a regulatory mediator of the juvenile-to-adult transition with functions differentiated from the angiosperm SOC1.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Comparative transcriptome analyses reveal two distinct transcriptional modules associated with pollen shedding time in pine

BACKGROUND

Seasonal flowering time is an ecologically and economically important trait in temperate trees. Previous studies have shown that temperature in many tree species plays a pivotal role in regulating flowering time. However, genetic control of flowering time is not synchronised in different individual trees under comparable temperature conditions, the underlying molecular mechanism is mainly to be investigated.

RESULTS

In the present study, we analysed the transcript abundance in male cones and needles from six early pollen-shedding trees (EPs) and six neighbouring late pollen-shedding trees (LPs) in Pinus tabuliformis at three consecutive time points in early spring. We found that the EPs and LPs had distinct preferred transcriptional modules in their male cones and, interestingly, the expression pattern was also consistently maintained in needles even during the winter dormancy period. Additionally, the preferred pattern in EPs was also adopted by other fast-growing tissues, such as elongating new shoots. Enhancement of nucleic acid synthesis and stress resistance pathways under cold conditions can facilitate rapid growth and maintain higher transcriptional activity.

CONCLUSIONS

During the cold winter and early spring seasons, the EPs were more sensitive to relatively warmer temperatures and showed higher transcriptomic activity than the LPs, indicating that EPs required less heat accumulation for pollen shedding than LPs. These results provided a transcriptomic-wide understanding of the temporal regulation of pollen shedding in pines.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

The sex expression and sex determining mechanism in Pistacia species

Generally, Pistacia species are dioecious, but monoecious strains in several populations have been found, providing excellent models for studying sex differentiation and sex determination mechanisms. Although the mechanisms of sex determination and sex evolution have been extensively studied, related research on heterozygous woody plants is limited. Here, we discuss the expressions of various sex types, which showed broad diversity and complex instability. We have also reviewed the sex determination systems in the plant kingdom and the morphological, cytological, physiological, and molecular aspects of the sex-linked markers in Pistacia trees. Moreover, hypotheses to explain the origin of monoecy are discussed, which is more likely to be the interaction between sex-related genes and environment factors in female plants. Besides, further prospects for the utilization of monoecious resources and the research directions of sex determination mechanism are proposed. This study provides important information on sex expression and provides more insights into sex differentiation and determination.

Pinus massoniana Introgression Hybrids Display Differential Expression of Reproductive Genes

Pinus massoniana and P. hwangshanensis are two conifer species located in southern China, which are of both economic and ornamental value. Around the middle and lower reaches of the Yangtze River, P. massoniana occurs mainly at altitudes below 700 m, while P. hwangshanensis can be found above 900 m. At altitudes where the distribution of both pines overlaps, a natural introgression hybrid exists, which we will further refer to as the Z pine. This pine has a morphological character that shares attributes of both P. massoniana and P. hwangshanensis. However, compared to the other two pines, its reproductive structure, the pinecone, has an ultra-low ripening rate with seeds that germinate poorly. In this study, we aimed to find the reason for the impaired cone maturation by comparing transcriptome libraries of P. massoniana and Z pine cones at seven successive growth stages. After sequencing and assembly, we obtained unigenes and then annotated them against NCBI’s non-redundant nucleotide and protein sequences, Swiss-Prot, Clusters of Orthologous Groups, Gene Ontology and KEGG Orthology databases. Gene expression levels were estimated and differentially expressed genes (DEGs) of the two pines were mined and analyzed. We found that several of them indeed relate to reproductive process. At every growth stage, these genes are expressed at a higher level in P. massoniana than in the Z pine. These data provide insight into understanding which molecular mechanisms are altered between P. massoniana and the Z pine that might cause changes in the reproductive process.

Partial correlation analysis of transcriptomes helps detangle the growth and defense network in spruce

In plants, there can be a trade-off between resource allocations to growth vs defense. Here, we use partial correlation analysis of gene expression to make inferences about the nature of this interaction. We studied segregating progenies of Interior spruce subject to weevil attack. In a controlled experiment, we measured pre-attack plant growth and post-attack damage with several morphological measures, and profiled transcriptomes of 188 progeny. We used partial correlations of individual transcripts (expressed sequence tags, ESTs) with pairs of growth/defense traits to identify important nodes and edges in the inferred underlying gene network, for example, those pairs of growth/defense traits with high mutual correlation with a single EST transcript. We give a method to identify such ESTs. A terpenoid ABC transporter gene showed strongest correlations (P = 0.019); its transcript represented a hub within the compact 166-member gene-gene interaction network (P = 0.004) of the negative genetic correlations between growth and subsequent pest attack. A small 21-member interaction network (P = 0.004) represented the uncovered positive correlations. Our study demonstrates partial correlation analysis identifies important gene networks underlying growth and susceptibility to the weevil in spruce. In particular, we found transcripts that strongly modify the trade-off between growth and defense, and allow identification of networks more central to the trade-off.

Transcriptome Analysis of Pinus Massoniana Lamb. Microstrobili During Sexual Reversal

The normal megastrobilli and microstrobilli before and after the sexual reversal in Pinus massoniana Lamb. were studied and classified using a transcriptomic approach. In the analysis, a total of 190,023 unigenes were obtained with an average length of 595 bp. The annotated unigenes were divided into 56 functional groups and 130 metabolic pathways involved in the physiological and biochemical processes related to ribosome biogenesis, carbon metabolism, and amino acid biosynthesis. Analysis revealed 4,758 differentially expressed genes (DEGs) between the mega- and microstrobili from the polycone twig. The DEGs between the mega- and microstrobili from the normal twig were 5,550. In the polycone twig, 1,188 DEGs were identified between the microstrobili and the sexually reversed megastrobili. Concerning plant hormone signal transduction pathways, the DEGs from both the normal and polycone twigs displayed distinct male or female associated expression patterns. There were 36 common hormone-related DEGs from the two types of twigs of P. massoniana. Interestingly, expression of these DEGs was up-regulated in the bisexual strobili, which underwent the sexual reversal. A portion of MADS-box genes in the bisexual strobili were up-regulated relative to expression in microstrobili.

Comparative Transcriptome Analysis of Male and Female Conelets and Development of Microsatellite Markers in Pinus bungeana, an Endemic Conifer in China

The sex determination in gymnosperms is still poorly characterized due to the lack of genomic/transcriptome resources and useful molecular genetic markers. To enhance our understanding of the molecular mechanisms of the determination of sexual recognition of reproductive structures in conifers, the transcriptome of male and female conelets were characterized in a Chinese endemic conifer species, Pinus bungeana Zucc. ex Endl. The 39.62 Gb high-throughput sequencing reads were obtained from two kinds of sexual conelets. After de novo assembly of the obtained reads, 85,305 unigenes were identified, 53,944 (63.23%) of which were annotated with public databases. A total of 12,073 differentially expressed genes were detected between the two types of sexes in P. bungeana, and 5766 (47.76%) of them were up-regulated in females. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enriched analysis suggested that some of the genes were significantly associated with the sex determination process of P. bungeana, such as those involved in tryptophan metabolism, zeatin biosynthesis, and cysteine and methionine metabolism, and the phenylpropanoid biosynthesis pathways. Meanwhile, some important plant hormone pathways (e.g., the gibberellin (GA) pathway, carotenoid biosynthesis, and brassinosteroid biosynthesis (BR) pathway) that affected sexual determination were also induced in P. bungeana. In addition, 8791 expressed sequence tag-simple sequence repeats (EST-SSRs) from 7859 unigenes were detected in P. bungeana. The most abundant repeat types were dinucleotides (1926), followed by trinucleotides (1711). The dominant classes of the sequence repeat were A/T (4942) in mononucleotides and AT/AT (1283) in dinucleotides. Among these EST-SSRs, 84 pairs of primers were randomly selected for the characterization of potential molecular genetic markers. Finally, 19 polymorphic EST-SSR primers were characterized. We found low to moderate levels of genetic diversity (N_A = 1.754; H_O = 0.206; H_E = 0.205) across natural populations of P. bungeana. The cluster analysis revealed two distinct genetic groups for the six populations that were sampled in this endemic species, which might be caused by the fragmentation of habitats and long-term geographic isolation among different populations. Taken together, this work provides important insights into the molecular mechanisms of sexual identity in the reproductive organs of P. bungeana. The molecular genetic resources that were identified in this study will also facilitate further studies in functional genomics and population genetics in the Pinus species.

The gibberellin GID1-DELLA signalling module exists in evolutionarily ancient conifers

Gibberellins (GAs) participate in controlling various aspects of basic plant growth responses. With the exception of bryophytes, GA signalling in land plants, such as lycophytes, ferns and angiosperms, is mediated via GIBBERELLIN-INSENSITIVE DWARF1 (GID1) and DELLA proteins. To explore whether this GID1-DELLA mechanism is present in pines, we cloned an orthologue (PtGID1) of Arabidopsis AtGID1a and two putative DELLA proteins (PtDPL; PtRGA) from Pinus tabuliformis, a widespread indigenous conifer species in China, and studied their recombinant proteins. PtGID1 shares with AtGID1a the conserved HSL motifs for GA binding and an N-terminal feature that are essential for interaction with DELLA proteins. Indeed, A. thaliana 35S:PtGID1 overexpressors showed a strong GA-hypersensitive phenotype compared to the wild type. Interactions between PtGID1 and PtDELLAs, but also interactions between the conifer-angiosperm counterparts (i.e. between AtGID1 and PtDELLAs and between PtGID1 and AtDELLA), were detected in vivo. This demonstrates that pine has functional GID1-DELLA components. The Δ17-domains within PtDPL and PtRGA were identified as potential interaction sites within PtDELLAs. Our results show that PtGID1 has the ability to interact with DELLA and functions as a GA receptor. Thus, a GA-GID1-DELLA signalling module also operates in evolutionarily ancient conifers.

Transcriptome profile analysis reveals the regulation mechanism of floral sex differentiation in Jatropha curcas L

The seeds of Jatropha curcas contain a high percentage of biodiesel. However, low seed yield which was limited by its poor female flowers was a bottleneck for its utilization. Here, we compared the transcriptomic profiles of five different samples during floral sex differentiation stages using Illumina Hiseq 4000. Our results showed that hundreds of differentially expressed genes (DEGs) were detected in floral sex initiation period, but thousands of DEGs were involved in the stamens and ovules development process. Moreover, the DEGs were mainly shown up-regulation in male floral initiation, but mainly down-regulation in female floral initiation. Male floral initiation was associated with the flavonoid biosynthesis pathway while female floral initiation was related to the phytohormone signal transduction pathway. Cytokinin (CTK) signaling triggered the initiation of female floral primordium, thereafter other phytohormones co-promoted the female floral development. In addition, the floral organ identity genes played important roles in floral sex differentiation process and displayed a general conservation of the ABCDE model in J. curcas. To the best of our knowledge, this data is the first comprehensive analysis of the underlying regulatory mechanism and the related genes during floral sex differentiation in J. curcas, which help in engineering high-yielding varieties of J. curcas.

Revisiting the phosphatidylethanolamine-binding protein (PEBP) gene family reveals cryptic FLOWERING LOCUS T gene homologs in gymnosperms and sheds new light on functional evolution

Angiosperms and gymnosperms are two major groups of extant seed plants. It has been suggested that gymnosperms lack FLOWERING LOCUS T (FT), a key integrator at the core of flowering pathways in angiosperms. Taking advantage of newly released gymnosperm genomes, we revisited the evolutionary history of the plant phosphatidylethanolamine-binding protein (PEBP) gene family through phylogenetic reconstruction. Expression patterns in three gymnosperm taxa and heterologous expression in Arabidopsis were studied to investigate the functions of gymnosperm FT-like and TERMINAL FLOWER 1 (TFL1)-like genes. Phylogenetic reconstruction suggests that an ancient gene duplication predating the divergence of seed plants gave rise to the FT and TFL1 genes. Expression patterns indicate that gymnosperm TFL1-like genes play a role in the reproductive development process, while GymFT1 and GymFT2, the FT-like genes resulting from a duplication event in the common ancestor of gymnosperms, function in both growth rhythm and sexual development pathways. When expressed in Arabidopsis, both spruce FT-like and TFL1-like genes repressed flowering. Our study demonstrates that gymnosperms do have FT-like and TFL1-like genes. Frequent gene and genome duplications contributed significantly to the expansion of the plant PEBP gene family. The expression patterns of gymnosperm PEBP genes provide novel insight into the functional evolution of this gene family.