Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber

Transcriptome analysis of gall oak (Quercus infectoria): De novo assembly, functional annotation and metabolic pathways analysis

Gall oak (Quercus infectoria) is a native tree of Iran, whose gall extract is used to treat many diseases. The presence of abundant secondary metabolites with various bioactivities in this plant has made it medically important. Despite its medicinal value, due to the lack of genomic information, the biosynthetic pathways of these compounds in this species are still unknown. The current research was aimed at observing, characterizing, and investigating the biosynthetic pathways of these compounds in Q.infectoria. De novo transcriptome assembly was conducted using the RNA sequencing technique. A total of 89,335 unigenes were generated, of which 6928 unigenes showed differential expression in leaves compared to root tissue. Gene ontology examination of DEGs revealed GO-term enrichment was related to cellular processes and enzyme activity. KEGG enrichment analysis for DEGs showed that most unigenes were related to metabolic pathways and biosynthesis of secondary metabolites. Moreover, 39 families of transcription factors were identified, of which the C2H2, bZIP, bHLH, and ERF TFs had the highest frequency. In the absence of a reference genome, the overall study of transcriptome will provide a reference for future functional and comparative studies. Moreover, the data obtained from sequencing and de novo assembly can be a valuable scientific resource for Q.infectoria.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

The Patterns of Male and Female Flowers in Flowering Stage May Not Be Optimal Resource Allocation for Fruit and Seed Growth

Changes in the proportions of male and female flowers in monoecious plants in response to external environmental conditions are directly related to the reproductive fitness of plants. The monoecious cucumber (Cucumber sativus) plant was used in this study to assess the responses of sex differentiation and the breeding process to nutrient supply and the degree of artificial pollination using pollen solutions of different concentrations. We found that the nutrient supply significantly improved the number of female flowers, while pollination treatments did not obviously increase the number of male flowers. Continuous pollination changed the number of female flowers especially in the later stage of the pollination experiment. Therefore, pollination changed the ratio of male and female flowers in the flowering stage of cucumber. Pollination treatment affected the fruit growth, seed set, and fruit yield. The number of fruit, fruit set percentage, and total seeds per plant did not increase with the pollination level, but individual fruit weight and seed number in one fruit did increase. The differentiation of male and female flowers in the flowering stage of cucumber is a response to nutrient and pollination resources, but this response is not the optimal resource allocation for subsequent fruit development and seed maturity, which suggests that the response of plants to external environment resources is short-term and direct.

The Dynamics of Flower Development in Castanea sativa Mill

The sweet chestnut tree (Castanea sativa Mill.) is one of the most significant Mediterranean tree species, being an important natural resource for the wood and fruit industries. It is a monoecious species, presenting unisexual male catkins and bisexual catkins, with the latter having distinct male and female flowers. Despite the importance of the sweet chestnut tree, little is known regarding the molecular mechanisms involved in the determination of sexual organ identity. Thus, the study of how the different flowers of C. sativa develop is fundamental to understand the reproductive success of this species and the impact of flower phenology on its productivity. In this study, a C. sativa de novo transcriptome was assembled and the homologous genes to those of the ABCDE model for floral organ identity were identified. Expression analysis showed that the C. sativa B- and C-class genes are differentially expressed in the male flowers and female flowers. Yeast two-hybrid analysis also suggested that changes in the canonical ABCDE protein-protein interactions may underlie the mechanisms necessary to the development of separate male and female flowers, as reported for the monoecious Fagaceae Quercus suber. The results here depicted constitute a step towards the understanding of the molecular mechanisms involved in unisexual flower development in C. sativa, also suggesting that the ABCDE model for flower organ identity may be molecularly conserved in the predominantly monoecious Fagaceae family.

Sex biased expression of hormone related genes at early stage of sex differentiation in papaya flowers

Sex types of papaya are controlled by a pair of nascent sex chromosomes, but molecular genetic mechanisms of sex determination and sex differentiation in papaya are still unclear. We performed comparative analysis of transcriptomic profiles of male and female floral buds at the early development stage before the initiation of reproductive organ primordia at which there is no morphological difference between male and female flowers. A total of 1734 differentially expressed genes (DEGs) were identified, of which 923 showed female-biased expression and 811 showed male-biased expression. Functional annotation revealed that genes related to plant hormone biosynthesis and signaling pathways, especially in abscisic acid and auxin pathways, were overrepresented in the DEGs. Transcription factor binding motifs, such as MYB2, GAMYB, and AP2/EREBP, were enriched in the promoters of the hormone-related DEGs, and transcription factors with those motifs also exhibited differential expression between sex types. Among these DEGs, we also identified 11 genes in the non-recombining region of the papaya sex chromosomes and 9 genes involved in stamen and carpel development. Our results suggested that sex differentiation in papaya may be regulated by multiple layers of regulation and coordination and involved transcriptional, epigenetic, and phytohormone regulation. Hormones, especially ABA and auxin, transcription factors, and genes in the non-recombination region of the sex chromosome could be involved in this process. Our findings may facilitate the elucidation of signal transduction and gene interaction in sex differentiation of unisexual flowers in papaya.

Sex-biased genes and metabolites explain morphologically sexual dimorphism and reproductive costs in Salix paraplesia catkins

Dioecious species evolved from species with monomorphic sex systems in order to achieve overall fitness gains by separating male and female functions. As reproductive organs, unisexual flowers have different reproductive roles and exhibit conspicuous sexual dimorphism. To date, little is known about the temporal variations in and molecular mechanisms underlying the morphology and reproductive costs of dioecious flowers. We investigated male and female flowers of Salix paraplesia in three flowering stages before pollination (the early, blooming and late stages) via transcriptional sequencing as well as metabolite content and phenotypic analysis. We found that a large number of sex-biased genes, rather than sex-limited genes, were responsible for sexual dimorphism in S. paraplesia flowers and that the variation in gene expression in male flowers intensified this situation throughout flower development. The temporal dynamics of sex-biased genes derived from changes in reproductive function during the different flowering stages. Sexually differentiated metabolites related to respiration and flavonoid biosynthesis exhibited the same bias directions as the sex-biased genes. These sex-biased genes were involved mainly in signal transduction, photosynthesis, respiration, cell proliferation, phytochrome biosynthesis, and phenol metabolism; therefore, they resulted in more biomass accumulation and higher energy consumption in male catkins. Our results indicated that sex-biased gene expression in S. paraplesia flowers is associated with different reproductive investments in unisexual flowers; male flowers require a greater reproductive investment to meet their higher biomass accumulation and energy consumption needs.

CorkOakDB-The Cork Oak Genome Database Portal

Quercus suber (cork oak) is an evergreen tree native to the Mediterranean basin, which plays a key role in the ecology and economy of this area. Over the last decades, this species has gone through an observable decline, mostly due to environmental factors. Deciphering the mechanisms of cork oak's response to the environment and getting a deep insight into its biology are crucial to counteract biotic and abiotic stresses compromising the stability of a unique ecosystem. In the light of these setbacks, the publication of the genome in 2018 was a major step towards understanding the genetic make-up of this species. In an effort to integrate this information in a comprehensive, accessible and intuitive format, we have developed The Cork Oak Genome Database Portal (CorkOakDB). The CorkOakDB is supported by the BioData.pt e-infrastructure, the Portuguese ELIXIR node for biological data. The portal gives public access to search and explore the curated genomic and transcriptomic data on this species. Moreover, CorkOakDB provides a user-friendly interface and functional tools to help the research community take advantage of the increased accessibility to genomic information. A study case is provided to highlight the functionalities of the portal. CorkOakDB guarantees the update, curation and data collection, aiming to collect data besides the genetic/genomic information, in order to become the main repository in cork oak research. Database URL: http://corkoakdb.org/.

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.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

BACKGROUND

Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms.

RESULTS

We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development.

CONCLUSIONS

Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

Unisexual flower initiation in the monoecious Quercus suber L.: a molecular approach

Several plant species display a temporal separation of the male and female flower organ development to enhance outbreeding; however, little is known regarding the genetic mechanisms controlling this temporal separation. Quercus suber is a monoecious oak tree with accentuated protandry: in late winter, unisexual male flowers emerge adjacent to the swollen buds, whereas unisexual female flowers emerge in the axils of newly formed leaves formed during spring (4-8 weeks after male flowering). Here, a phylogenetic profiling has led to the identification of cork oak homologs of key floral regulatory genes. The role of these cork oak homologs during flower development was identified with functional studies in Arabidopsis thaliana. The expression profile throughout the year of flower regulators (inducers and repressors), in leaves and buds, suggests that the development of male and female flowers may be preceded by separated induction events. Female flowers are most likely induced during the vegetative flush occurring in spring, whereas male flowers may be induced in early summer. Male flowers stay enclosed within the pre-dormant buds, but complete their development before the vegetative flush of the following year, displaying a long period of anthesis that spans the dormant period. Our results portray a genetic mechanism that may explain similar reproductive habits in other monoecious tree species.

Transcriptomic differences between male and female Trachycarpus fortunei

Trachycarpus fortunei (Hook.) is a typical dioecious plant, which has important economic value. There is currently no sex identification method for the early stages of T. fortunei growth. The aim of this study was to obtain expression and site differences between male and female T. fortunei transcriptomes. Using the Illumina sequencing platform, the transcriptomes of T. fortunei male and female plants were sequenced. By analyzing transcriptomic differences, the chromosomal helical binding protein (CHD1), serine/threonine protein kinase (STPK), cytochrome P450 716B1, and UPF0136 were found to be specifically expressed in T. fortunei males. After single nucleotide polymorphism (SNP) detection, a total of 12 male specific sites were found and the THUMP domain protein homologs were found to be male-biased expressed. Cytokinin dehydrogenase 6 (CKX6) was upregulated in male flowers and the lower concentrations of cytokinin (CTK) may be more conducive to male flower development. During new leaf growth, flavonoid and flavonol biosynthesis were initiated. Additionally, the flavonoids, 3',5'-hydroxylase (F3'5'H), flavonoids 3'-hydroxylase, were upregulated, which may cause the pale yellow phenotype. Based on these data, it can be concluded that inter-sex differentially expressed genes (DEGs) and specific SNP loci may be associated with sex determination in T. fortunei.

Genome-Wide Identification of Epigenetic Regulators in Quercus suber L

Modifications of DNA and histones, including methylation and acetylation, are critical for the epigenetic regulation of gene expression during plant development, particularly during environmental adaptation processes. However, information on the enzymes catalyzing all these modifications in trees, such as Quercus suber L., is still not available. In this study, eight DNA methyltransferases (DNA Mtases) and three DNA demethylases (DDMEs) were identified in Q. suber. Histone modifiers involved in methylation (35), demethylation (26), acetylation (8), and deacetylation (22) were also identified in Q. suber. In silico analysis showed that some Q. suber DNA Mtases, DDMEs and histone modifiers have the typical domains found in the plant model Arabidopsis, which might suggest a conserved functional role. Additional phylogenetic analyses of the DNA and histone modifier proteins were performed using several plant species homologs, enabling the classification of the Q. suber proteins. A link between the expression levels of each gene in different Q. suber tissues (buds, flowers, acorns, embryos, cork, and roots) with the functions already known for their closest homologs in other species was also established. Therefore, the data generated here will be important for future studies exploring the role of epigenetic regulators in this economically important species.

Flowering Phenology and Characteristics of Pollen Aeroparticles of Quercus Species in Korea

In recent decades, airborne allergens for allergic respiratory diseases have been found to increase significantly by a process of converting coniferous forests into broad-leaved forests in Korea. This study was conducted to evaluate factors, including airborne pollen counts, micromorphology, and flowering phenology, that can affect oak pollen-related allergic symptoms. The catkin of Mongolian oak (Quercus mongolica Fisch. ex Ledeb.) showed the most rapidly blooming catkin on Julian day 104 in flower development. Among six species, the last flowering was observed on Julian day 119 in Korean oak (Quercus dentata Thunb.). The pollen dispersal was persisted for about 32 days from Julian day 104 to Julian day 136. Airborne pollen was observed about 2 weeks after flowering phase H, the senescence phase. Pollen size varied by species, with the largest from Q. mongolica (polar axis length, PL = 31.72 µm, equatorial axis length, EL = 39.05 µm) and the smallest from Jolcham oak (Quercus serrata Murray) (PL = 26.47 µm, EL = 32.32 µm). Regarding pollen wall structure, endexine of Q. dentata was coarsely laminated or fragmented. The endexine thicknesses of Sawtooth oak (Quercus acutissima Carruth.) and Q. serrata were thick and stable, whereas Galcham oak (Quercus aliena Blume), Q. mongolica, and Oriental cork oak (Quercus variabilis Blume) had thinner endexines. The area occupied by pollenkitt of Q. variabilis was significantly larger than that of Q. acutissima. Importantly, Q. variabilis had a distinctly thick 17 kDa protein band, a presumed major allergen. Oak species differ in pollen protein composition, and thus there is a possibility that the allergenic activity of pollen proteins vary depending on oak species. This study highlights the fact that native oak species in Korea differ in flowering pattern of male flowers, pollen morphology, and pollen chemical constituents. These discrepancies in flowering and pollen properties imply variable allergic responses to oak pollen from different species.

Different expression pattern of flowering pathway genes contribute to male or female organ development during floral transition in the monoecious weed Ambrosia artemisiifolia L. (Asteraceae)

The highly allergenic and invasive weed Ambrosia artemisiifolia L. is a monoecius plant with separated male and female flowers. The genetic regulation of floral morphogenesis is a less understood field in the reproduction biology of this species. Therefore the objective of this work was to investigate the genetic control of sex determination during floral organogenesis. To this end, we performed a genome-wide transcriptional profiling of vegetative and generative tissues during the plant development comparing wild-growing and in vitro cultivated plants. RNA-seq on Illumina NextSeq 500 platform with an integrative bioinformatics analysis indicated differences in 80 floral gene expressions depending on photoperiodic and endogenous initial signals. Sex specificity of genes was validated based on RT-qPCR experiments. We found 11 and 16 uniquely expressed genes in female and male transcriptomes that were responsible particularly to maintain fertility and against abiotic stress. High gene expression of homologous such as FD, FT, TFL1 and CAL, SOC1, AP1 were characteristic to male and female floral meristems during organogenesis. Homologues transcripts of LFY and FLC were not found in the investigated generative and vegetative tissues. The repression of AP1 by TFL1 homolog was demonstrated in male flowers resulting exclusive expression of AP2 and PI that controlled stamen and carpel formation in the generative phase. Alterations of male and female floral meristem differentiation were demonstrated under photoperiodic and hormonal condition changes by applying in vitro treatments.

New Insights of Salicylic Acid Into Stamen Abortion of Female Flowers in Tung Tree (Vernicia fordii)

Tung tree (Vernicia fordii), an economically important woody oil plant, is a monoecious and diclinous species with male and female flowers on the same inflorescence. The extremely low proportion of female flowers leads to low fruit yield in tung orchards. The female flower normally develops along with stamen abortion; otherwise sterile ovules will be produced. However, little knowledge is known about the molecular basis of the female flower development in tung tree. In this study, integrated analyses of morphological and cytological observations, endogenous phytohormone assay and RNA-seq were conducted to understand the molecular mechanism of the female flower development in tung tree. Cytological observation suggested that the abortion of stamens in female flowers (SFFs) belongs to the type of programmed cell death (PCD), which was caused by tapetum degeneration at microspore mother cell stage. A total of 1,366 differentially expressed genes (DEGs) were identified in female flowers by RNA-seq analysis, of which 279 (20.42%) DEGs were significantly enriched in phenylpropanoid biosynthesis, phenylalanine metabolism, flavonoid biosynthesis, starch and sucrose metabolism, and plant hormone signal transduction. Stage-specific transcript identification detected dynamically expressed genes of important transcription regulators in female flowers that may be involved in PCD and floral organ development. Gene expression patterns revealed that 17 anther and pollen development genes and 37 PCD-related genes might be involved in the abortion of SFF. Further analyses of phytohormone levels and co-expression networks suggested that salicylic acid (SA) accumulation could trigger PCD and inhibit the development of SFF in tung tree. This study provides new insights into the role of SA in regulating the abortion of SFF to develop normal female flowers.

The draft genome sequence of cork oak

Cork oak (Quercus suber) is native to southwest Europe and northwest Africa where it plays a crucial environmental and economical role. To tackle the cork oak production and industrial challenges, advanced research is imperative but dependent on the availability of a sequenced genome. To address this, we produced the first draft version of the cork oak genome. We followed a de novo assembly strategy based on high-throughput sequence data, which generated a draft genome comprising 23,347 scaffolds and 953.3 Mb in size. A total of 79,752 genes and 83,814 transcripts were predicted, including 33,658 high-confidence genes. An InterPro signature assignment was detected for 69,218 transcripts, which represented 82.6% of the total. Validation studies demonstrated the genome assembly and annotation completeness and highlighted the usefulness of the draft genome for read mapping of high-throughput sequence data generated using different protocols. All data generated is available through the public databases where it was deposited, being therefore ready to use by the academic and industry communities working on cork oak and/or related species.

De novo transcriptome assembly and comparative analysis between male and benzyladenine-induced female inflorescence buds of Plukenetia volubilis

Plukenetia volubilis is a promising oilseed crop due to its seeds being rich in unsaturated fatty acids, especially alpha-linolenic acid. P. volubilis is monoecious, with separate male and female flowers on the same inflorescence. We previously reported that male flowers were converted to female flowers by exogenous cytokinin (6-benzyladenine, 6-BA) treatment in P. volubilis. To identify candidate genes associated with floral sex differentiation of P. volubilis, we performed de novo transcriptome assembly and comparative analysis on control male inflorescence buds (MIB) and female inflorescence buds (FIB) induced by 6-BA using Illumina sequencing technology. A total of 57,664 unigenes with an average length of 979 bp were assembled from 104.1 million clean reads, and 45,235 (78.45%) unigenes were successfully annotated in the public databases. Notably, Gene Ontology analyses revealed that 4193 and 3880 unigenes were enriched in the categories of reproduction and reproductive processes, respectively. Differential expression analysis identified 1385 differentially expressed unigenes between MIB and FIB, of which six unigenes related to cytokinin and auxin signaling pathways and 16 important transcription factor (TF) genes including MADS-box family members were identified. In particular, several unigenes encoding important TFs, such as homologs of CRABS CLAW, RADIALIS-like 1, RADIALIS-like 2, HECATE 2, WUSCHEL-related homeobox 9, and SUPERMAN, were expressed at higher levels in FIB than in MIB. The expression patterns of the 36 selected unigenes revealed by transcriptome analysis were successfully validated by quantitative real-time PCR. This study not only provides comprehensive gene expression profiles of P. volubilis inflorescence buds, but also lays the foundation for research on the molecular mechanism of floral sex determination in P. volubilis and other monoecious plants.

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.

Role of floral organ identity genes in the development of unisexual flowers of Quercus suber L

Monoecious species provide an excellent system to study the specific determinants that underlie male and female flower development. Quercus suber is a monoecious species with unisexual flowers at inception. Despite the overall importance of this and other tree species with a similar reproductive habit, little is known regarding the mechanisms involved in the development of their male and female flowers. Here, we have characterised members of the ABCDE MADS-box gene family of Q. suber. The temporal expression of these genes was found to be sex-biased. The B-class genes, in particular, are predominantly, or exclusively (in the case of QsPISTILLATA), expressed in the male flowers. Functional analysis in Arabidopsis suggests that the B-class genes have their function conserved. The identification of sex-biased gene expression plus the identification of unusual protein-protein interactions suggest that the floral organ identity of Q. suber may be under control of specific changes in the dynamics of the ABCDE model. This study constitutes a major step towards the characterisation of the mechanisms involved in reproductive organ identity in a monoecious tree with a potential contribution towards the knowledge of conserved developmental mechanisms in other species with a similar sex habit.

Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis)

BACKGROUND

Garden asparagus (Asparagus officinalis) is a highly valuable vegetable crop of commercial and nutritional interest. It is also commonly used to investigate the mechanisms of sex determination and differentiation in plants. However, the sex expression mechanisms in asparagus remain poorly understood.

RESULTS

De novo transcriptome sequencing via Illumina paired-end sequencing revealed more than 26 billion bases of high-quality sequence data from male and female asparagus flower buds. A total of 72,626 unigenes with an average length of 979 bp were assembled. In comparative transcriptome analysis, 4876 differentially expressed genes (DEGs) were identified in the possible sex-determining stage of female and male/supermale flower buds. Of these DEGs, 433, including 285 male/supermale-biased and 149 female-biased genes, were annotated as flower related. Of the male/supermale-biased flower-related genes, 102 were probably involved in anther development. In addition, 43 DEGs implicated in hormone response and biosynthesis putatively associated with sex expression and reproduction were discovered. Moreover, 128 transcription factor (TF)-related genes belonging to various families were found to be differentially expressed, and this finding implied the essential roles of TF in sex determination or differentiation in asparagus. Correlation analysis indicated that miRNA-DEG pairs were also implicated in asparagus sexual development.

CONCLUSIONS

Our study identified a large number of DEGs involved in the sex expression and reproduction of asparagus, including known genes participating in plant reproduction, plant hormone signaling, TF encoding, and genes with unclear functions. We also found that miRNAs might be involved in the sex differentiation process. Our study could provide a valuable basis for further investigations on the regulatory networks of sex determination and differentiation in asparagus and facilitate further genetic and genomic studies on this dioecious species.

The quest for epigenetic regulation underlying unisexual flower development in Cucumis melo

BACKGROUND

Melon (Cucumis melo) is an important vegetable crop from the Cucurbitaceae family and a reference model specie for sex determination, fruit ripening and vascular fluxes studies. Nevertheless, the nature and role of its epigenome in gene expression regulation and more specifically in sex determination remains largely unknown.

RESULTS

We have investigated genome wide H3K27me3 and H3K9ac histone modifications and gene expression dynamics, in five melon organs. H3K9ac and H3K27me3 were mainly distributed along gene-rich regions and constrained to gene bodies. H3K9ac was preferentially located at the TSS, whereas H3K27me3 distributed uniformly from TSS to TES. As observed in other species, H3K9ac and H3K27me3 correlated with high and low gene expression levels, respectively. Comparative analyses of unisexual flowers pointed out sex-specific epigenetic states of TFs involved in ethylene response and flower development. Chip-qPCR analysis of laser dissected carpel and stamina primordia, revealed sex-specific histone modification of MADS-box genes. Using sex transition mutants, we demonstrated that the female promoting gene, CmACS11, represses the expression of the male promoting gene CmWIP1 via deposition of H3K27me3.

CONCLUSIONS

Our findings reveal the organ-specific landscapes of H3K9ac and H3K27me3 in melon. Our results also provide evidence that the sex determination genes recruit histone modifiers to orchestrate unisexual flower development in monoecious species.

Tissue-specific transcriptomic profiling provides new insights into the reproductive ecology and biology of the iconic seagrass species Posidonia oceanica

Seagrasses form extensive meadows in shallow coastal waters and are among the world's most productive ecosystems. Seagrasses can produce both clonally and sexually, and flowering has long been considered infrequent, but important for maintaining genetically diverse stands. Here we investigate the molecular mechanisms involved in flowering of the seagrass Posidonia oceanica, an iconic species endemic to the Mediterranean. We generated a de novo transcriptome of this non-model species for leaf, male and female flower tissue of three individuals, and present molecular evidence for genes that may be involved in the flowering process and on the reproductive biology of the species. We present evidence that suggests that P. oceanica exhibits a strategy of protogyny, where the female part of the hermaphroditic flower develops before the male part, in order to avoid self-fertilization. We found photosynthetic genes to be up-regulated in the female flower tissues, indicating that this may be capable of photosynthesis. Finally, we detected a number of interesting genes, previously known to be involved in flowering pathways responding to light and temperature cues and in pathways involved in anthocyanin and exine synthesis. This first comparative transcriptomic approach of leaf, male and female tissue provides a basis for functional genomics research on flower development in P. oceanica and other seagrass species.

Transcriptome Analysis of a Female-sterile Mutant (fsm) in Chinese Cabbage (Brassica campestris ssp. pekinensis)

Female-sterile mutants are ideal materials for studying pistil development in plants. Here, we identified a female-sterile mutant fsm in Chinese cabbage. This mutant, which exhibited stable inheritance, was derived from Chinese cabbage DH line 'FT' using a combination of isolated microspore culture and ethyl methanesulfonate mutagenesis. Compared with the wild-type line 'FT,' the fsm plants exhibited pistil abortion, and floral organs were also relatively smaller. Genetic analysis indicated that the phenotype of fsm is controlled by a single recessive nuclear gene. Morphological observations revealed that the presence of abnormal ovules in fsm likely influenced normal fertilization process, ultimately leading to female sterility. Comparative transcriptome analysis on the flower buds of 'FT' and fsm using RNA-Seq revealed a total of 1,872 differentially expressed genes (DEGs). Of these, a number of genes involved in pistil development were identified, such as PRETTY FEW SEEDS 2 (PFS2), temperature-induced lipocalin (TIL), AGAMOUS-LIKE (AGL), and HECATE (HEC). Furthermore, GO and KEGG pathway enrichment analyses of the DEGs suggested that a variety of biological processes and metabolic pathways are significantly enriched during pistil development. In addition, the expression patterns of 16 DEGs, including four pistil development-related genes and 12 floral organ development-related genes, were analyzed using qRT-PCR. A total of 31,272 single nucleotide polymorphisms were specifically detected in fsm. These results contribute to shed light on the regulatory mechanisms underlying pistil development in Chinese cabbage.

Bursaphelenchus xylophilus and B. mucronatus secretomes: a comparative proteomic analysis

The pinewood nematode, Bursaphelenchus xylophilus, recognized as a worldwide major forest pest, is a migratory endoparasitic nematode with capacity to feed on pine tissues and also on fungi colonizing the trees. Bursaphelenchus mucronatus, the closest related species, differs from B. xylophilus on its pathogenicity, making this nematode a good candidate for comparative analyses. Secretome profiles of B. xylophilus and B. mucronatus were obtained and proteomic differences were evaluated by quantitative SWATH-MS. From the 681 proteins initially identified, 422 were quantified and compared between B. xylophilus and B. mucronatus secretomes and from these, 243 proteins were found differentially regulated: 158 and 85 proteins were increased in B. xylophilus and B. mucronatus secretomes, respectively. While increased proteins in B. xylophilus secretome revealed a strong enrichment in proteins with peptidase activity, the increased proteins in B. mucronatus secretome were mainly related to oxidative stress responses. The changes in peptidases were evaluated at the transcription level by RT-qPCR, revealing a correlation between the mRNA levels of four cysteine peptidases with secretion levels. The analysis presented expands our knowledge about molecular basis of B. xylophilus and B. mucronatus hosts interaction and supports the hypothesis of a key role of secreted peptidases in B. xylophilus pathogenicity.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.

Characterization of the cork oak transcriptome dynamics during acorn development

BACKGROUND

Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water.

RESULTS

A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability.

CONCLUSIONS

To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.

RNA-Seq and Gene Network Analysis Uncover Activation of an ABA-Dependent Signalosome During the Cork Oak Root Response to Drought

Quercus suber (cork oak) is a West Mediterranean species of key economic interest, being extensively explored for its ability to generate cork. Like other Mediterranean plants, Q. suber is significantly threatened by climatic changes, imposing the need to quickly understand its physiological and molecular adaptability to drought stress imposition. In the present report, we uncovered the differential transcriptome of Q. suber roots exposed to long-term drought, using an RNA-Seq approach. 454-sequencing reads were used to de novo assemble a reference transcriptome, and mapping of reads allowed the identification of 546 differentially expressed unigenes. These were enriched in both effector genes (e.g., LEA, chaperones, transporters) as well as regulatory genes, including transcription factors (TFs) belonging to various different classes, and genes associated with protein turnover. To further extend functional characterization, we identified the orthologs of differentially expressed unigenes in the model species Arabidopsis thaliana, which then allowed us to perform in silico functional inference, including gene network analysis for protein function, protein subcellular localization and gene co-expression, and in silico enrichment analysis for TFs and cis-elements. Results indicated the existence of extensive transcriptional regulatory events, including activation of ABA-responsive genes and ABF-dependent signaling. We were then able to establish that a core ABA-signaling pathway involving PP2C-SnRK2-ABF components was induced in stressed Q. suber roots, identifying a key mechanism in this species' response to drought.