Transcriptome analysis of the differentially expressed genes in the male and female shrub willows (Salix suchowensis)

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

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.

Mapping the sex determination region in the Salix F1 hybrid common parent population confirms a ZW system in six diverse species

Within the genus Salix, there are approximately 350 species native primarily to the northern hemisphere and adapted to a wide range of habitats. This diversity can be exploited to mine novel alleles conferring variation important for production as a bioenergy crop, but also to identify evolutionarily important genes, such as those involved in sex determination. To leverage this diversity, we created a mapping population by crossing 6 Salix species (Salix viminalis, Salix suchowensis, Salix integra, Salix koriyanagi, Salix udensis, and Salix alberti) to common male and female Salix purpurea parents. Each family was genotyped via genotyping-by-sequencing and assessed for kinship and population structure as well as the construction of 16 backcross linkage maps to be used as a genetic resource for breeding and selection. Analyses of population structure resolved both the parents and F₁ progeny to their respective phylogenetic section and indicated that the S. alberti parent was misidentified and was most likely S.suchowensis. Sex determining regions were identified on Salix chromosome 15 in the female-informative maps for seven of the eight families indicating that these species share a common female heterogametic ZW sex system. The eighth family, S. integra × S. purpurea, was entirely female and had a truncated chromosome 15. Beyond sex determination, the Salix F₁ hybrid common parent population (Salix F₁ HCP) introduced here will be useful in characterizing genetic factors underlying complex traits, aid in marker-assisted selection, and support genome assemblies for this promising bioenergy crop.

Distribution of simple sequence repeats, transcription factors, and differentially expressed genes in the NGS-based transcriptome of male and female seabuckthorn (Hippophae salicifolia)

Seabuckthorn (Hippophae salicifolia) is a perennial, multipurpose wonder plant, popular for its immense medicinal, nutritional, and therapeutic properties. However, due to the lack of whole-genome-based studies, the molecular mechanism governing distinct sexual phenotypes is still not clear. We employed the high-throughput NGS Illumina NovaSeq paired-end technology to generate whole transcriptome profiles of male and female plants of H. salicifolia. In total, 3.2 million raw short reads were generated with an average length of 150 bp, including 50911358 reads from the male leaf tissue samples and 45850364 reads from the female leaf tissue samples. Clustering of the high-quality reads yielded de novo short read assembly of 50259 transcripts of >100 bp length. The final transcripts were assigned Gene Ontology (GO) terms. The digital expression of genes was studied using the DESeq2 of R package that identified 7180 differentially expressed genes (DEGs) between the male and female plant samples. Further, 10,850 simple sequence repeats, and 8,351 transcription factors, distributed in more than 85 transcription families, were also mined from the final assembled transcriptome. Next, COG and KEGG pathway analyses were performed to assign biological functional terms to the DEGs. The findings of the present study will provide a valuable resource for gene expression discovery and other functional genomics studies aiming towards the selection of candidate genes for the development of sex-specific markers in seabuckthorn and other closely related species.Communicated by Ramaswamy H. Sarma.

Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Sex dimorphism and gene expression were studied in developing catkins in 159 F2 individuals from the bioenergy crop Salix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

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 analysis of microsatellite and sex-linked marker identification in Gleditsia sinensis

BACKGROUND

Gleditsia sinensis Lam. (Leguminosae), a dioecious perennial arbor, demonstrates important medicinal properties and economic value. These properties can be harnessed depending on the sex of the plant. However, the sex of the plants is difficult to identify accurately through morphological methods before the flowering.

RESULTS

We used bulked segregant analysis to screen sex-specific simple sequence repeat (SSR) markers in G. sinensis. Five male and five female plants were pooled to form the male and female bulks, respectively, and subjected to whole-genome sequencing. After high-throughput sequencing, 5,350,359 sequences were obtained, in which 2,065,210 SSRs were searched. Among them, the number of duplicated SSRs was the highest. The male plants could reach 857,874, which accounted for 60.86% of the total number of male plants. The female plants could reach 1,447,603, which accounted for 56.25% of the total model of the female plants. Among all the nucleotide repeat types, the A/T-rich motif was the most abundant. A total of 309,516 female strain-specific SSRs were selected by clustering. After designing the primers, the male and female gene pools were amplified, and five pairs of primers (i.e., 27, 34, 36, 39, and 41) were found to amplify the differential bands in the male and female gene pools. Using the five pairs of primers, we performed PCR verification on 10 individuals of known sex, which constructed the gene pool. The female plants amplified a single fragment of lengths (i.e., 186, 305, 266, 203, and 260 bp) and no male plant strip, thereby completing the identification of the male and female sexes of the G. sinensis.

CONCLUSIONS

This study provides accurate sex identification strategies between female and male plants, thus improving the utilization rate of G. sinensis resources.

Towards sex identification of Asian Palmyra palm (Borassus flabellifer L.) by DNA fingerprinting, suppression subtractive hybridization and de novo transcriptome sequencing

BACKGROUND

Asian Palmyra palm, the source of palm-sugar, is dioecious with a long juvenile period requiring at least 12 years to reach its maturity. To date, there is no reliable molecular marker for identifying sexes before the first bloom, limiting crop designs and utilization. We aimed to identify sex-linked markers for this palm using PCR-based DNA fingerprinting, suppression subtractive hybridization (SSH) and transcriptome sequencing.

METHODS

DNA fingerprints were generated between males and females based on RAPD, AFLP, SCoT, modified SCoT, ILP, and SSR techniques. Large-scale cloning and screening of SSH libraries and de novo transcriptome sequencing of male and female cDNA from inflorescences were performed to identify sex-specific genes for developing sex-linked markers.

RESULTS

Through extensive screening and re-testing of the DNA fingerprints (up to 1,204 primer pairs) and transcripts from SSH (>10,000 clones) and transcriptome data, however, no sex-linked marker was identified. Although de novo transcriptome sequencing of male and female inflorescences provided ∼32 million reads and 187,083 assembled transcripts, PCR analysis of selected sex-highly represented transcripts did not yield any sex-linked marker. This result may suggest the complexity and small sex-determining region of the Asian Palmyra palm. To this end, we provide the first global transcripts of male and female inflorescences of Asian Palmyra palm. Interestingly, sequence annotation revealed a large proportion of transcripts related to sucrose metabolism, which corresponds to the sucrose-rich sap produced in the inflorescences, and these transcripts will be useful for further understanding of sucrose production in sugar crop plants. Provided lists of sex-specific and differential-expressed transcripts would be beneficial to the further study of sexual development and sex-linked markers in palms and related species.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

Comparative genomics and transcriptomics analysis reveals evolution patterns of selection in the Salix phylogeny

BACKGROUND

Willows are widely distributed in the northern hemisphere and have good adaptability to different living environment. The increasing of genome and transcriptome data provides a chance for comparative analysis to study the evolution patterns with the different origin and geographical distributions in the Salix phylogeny.

RESULTS

Transcript sequences of 10 Salicaceae species were downloaded from public databases. All pairwise of orthologues were identified by comparative analysis in these species, from which we constructed a phylogenetic tree and estimated the rate of diverse. Divergence times were estimated in the 10 Salicaceae using comparative transcriptomic analysis. All of the fast-evolving positive selection sequences were identified, and some cold-, drought-, light-, universal-, and heat- resistance genes were discovered.

CONCLUSIONS

The divergence time of subgenus Vetrix and Salix was about 17.6-16.0 Mya during the period of Middle Miocene Climate Transition (21-14 Mya). Subgenus Vetrix diverged to migratory and resident groups when the climate changed to the cool and dry trend by 14 Mya. Cold- and light- stress genes were involved in positive selection among the resident Vetrix, and which would help them to adapt the cooling stage. Universal- stress genes exhibited positive selection among the migratory group and subgenus Salix. These data are useful for comprehending the adaptive evolution and speciation in the Salix lineage.

Sex-biased gene expression in flowers, but not leaves, reveals secondary sexual dimorphism in Populus balsamifera

Because sexual dimorphism in plants is often less morphologically conspicuous than in animals, studies of sex-biased gene expression may provide a quantitative metric to better address their commonality, molecular pathways, consistency across tissues and taxa, and evolution. The presence of sex-biased gene expression in tissues other than the androecium or gynoecium, termed secondary sexual characters, suggests that these traits arose after the initial evolution of dioecy. Patterns of sequence evolution may provide evidence of positive selection that drove sexual specialization. We compared gene expression in male and female flowers and leaves of Populus balsamifera to assess the extent of sex-biased expression, and tested whether sex-biased genes exhibit elevated rates of protein evolution. Sex-biased expression was pervasive in floral tissue, but nearly absent in leaf tissue. Female-biased genes in flowers were associated with photosynthesis, whereas male-biased genes were associated with mitochondrial function. Sex-biased genes did not exhibit elevated rates of protein evolution, contrary to results from other studies in animals and plants. Our results suggest that the ecological and physiological constraints associated with the energetics of flowering, rather than sexual conflict, have probably shaped the differences in male and female gene expression in P. balsamifera.

Comparative de novo transcriptome analysis of male and female Sea buckthorn

Sea buckthorn is a dioecious medicinal plant found at high altitude. The plant has both male and female reproductive organs in separate individuals. In this article, whole transcriptome de novo assemblies of male and female flower bud samples were carried out using Illumina NextSeq 500 platform to determine the role of the genes involved in sex determination. Moreover, genes with differential expression in male and female transcriptomes were identified to understand the underlying sex determination mechanism. The current study showed 63,904 and 62,272 coding sequences (CDS) in female and male transcriptome data sets, respectively. 16,831 common CDS were screened out from both transcriptomes, out of which 625 were upregulated and 491 were found to be downregulated. To understand the potential regulatory roles of differentially expressed genes in metabolic networks and biosynthetic pathways: KEGG mapping, gene ontology, and co-expression network analysis were performed. Comparison with Flowering Interactive Database (FLOR-ID) resulted in eight differentially expressed genes viz. CHD3-type chromatin-remodeling factor PICKLE (PKL), phytochrome-associated serine/threonine-protein phosphatase (FYPP), protein TOPLESS (TPL), sensitive to freezing 6 (SFR6), lysine-specific histone demethylase 1 homolog 1 (LDL1), pre-mRNA-processing-splicing factor 8A (PRP8A), sucrose synthase 4 (SUS4), ubiquitin carboxyl-terminal hydrolase 12 (UBP12), known to be broadly involved in flowering, photoperiodism, embryo development, and cold response pathways. Male and female flower bud transcriptome data of Sea buckthorn may provide comprehensive information at genomic level for the identification of genetic regulation involved in sex determination.

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.

Identification of the Sex-Biased Gene Expression and Putative Sex-Associated Genes in Eucommia ulmoides Oliver Using Comparative Transcriptome Analyses

Eucommia ulmoides is a model representative of the dioecious plants with sex differentiation at initiation. Nevertheless, the genetic mechanisms of sexual dimorphism and sex determination in E. ulmoides remain poorly understood. In this study de novo transcriptome sequencing on Illumina platform generated >45 billion high-quality bases from fresh leaves of six male and female individuals of E. ulmoides. A total of 148,595 unigenes with an average length of 801 base-pairs (bp) were assembled. Through comparative transcriptome analyses, 116 differentially expressed genes (DEGs) between the males and the females were detected, including 73 male-biased genes and 43 female-biased genes. Of these DEGs, three female-biased genes were annotated to be related with the sexually dimorphic gutta content in E. ulmoides. One male-biased DEG was identified as putative MADS box gene APETALA3, a B class floral organ identity gene in the flowering plants. SNPs calling analyses further confirmed that the APETALA3-like gene was probably involved in the sex determination in E. ulmoides. Four other male-biased DEGs were potential sex-associated genes as well with segregated SNPs in accord with sex type. In addition, the SNPs density was 1.02 per kilobase (kb) in the expressed genes of E. ulmoides, implying a relatively high genetic diversity.

Dominance and Sexual Dimorphism Pervade the Salix purpurea L. Transcriptome

The heritability of gene expression is critical in understanding heterosis and is dependent on allele-specific regulation by local and remote factors in the genome. We used RNA-Seq to test whether variation in gene expression among F1 and F2 intraspecific Salix purpurea progeny is attributable to cis- and trans-regulatory divergence. We assessed the mode of inheritance based on gene expression levels and allele-specific expression for F1 and F2 intraspecific progeny in two distinct tissue types: shoot tip and stem internode. In addition, we explored sexually dimorphic patterns of inheritance and regulatory divergence among F1 progeny individuals. We show that in S. purpurea intraspecific crosses, gene expression inheritance largely exhibits a maternal dominant pattern, regardless of tissue type or pedigree. A significantly greater number of cis- and trans-regulated genes coincided with upregulation of the maternal parent allele in the progeny, irrespective of the magnitude, whereas the paternal allele was higher expressed for genes showing cis × trans or compensatory regulation. Importantly, consistent with previous genetic mapping results for sex in shrub willow, we have delimited sex-biased gene expression to a 2 Mb pericentromeric region on S. purpurea chr15 and further refined the sex determination region. Altogether, our results offer insight into the inheritance of gene expression in S. purpurea as well as evidence of sexually dimorphic expression which may have contributed to the evolution of dioecy in Salix.

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.

Detecting the Candidate Gender Determinants by Bioinformatic Prediction of miRNAs and Their Targets from Transcriptome Sequences of the Male and Female Flowers in Salix suchowensis

MicroRNAs (miRNAs) belong to a class of small, noncoding, and endogenous single-stranded RNAs that negatively regulate gene expression at the posttranscriptional level. Potential miRNAs can be identified based on sequence homology since miRNAs are highly conserved in plants. In this study, we aligned the expressed sequence tags derived from flower buds of male and female S. suchowensis to miRNAs in the miRBase, which enable us to identify 34 potential miRNAs from flower buds of the alternate sexes. Among them, 11 were from the female and 23 were from the male. Analyzing sequence complementarity led to identification of 124 and 55 miRNA targets in the male and female flower buds, respectively. By mapping the target genes of the predicted miRNAs to the sequence assemblies of S. suchowensis, a miR156 mediated gene was detected at the gender locus of willow, which was a transcription factor involved in flower development. It is noteworthy that this target is not expressed in male flower, while it is expressed fairly highly in female flower based on the transcriptome data derived from the alternate sexes of willows. This study provides new bioinformatic clue for further exploring the genetic mechanism underlying gender determination in willows.

Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis

Willow is a widely used dioecious woody plant of Salicaceae family in China. Due to their high biomass yields, willows are promising sources for bioenergy crops. In this study, we assembled the complete mitochondrial (mt) genome sequence of S. suchowensis with the length of 644,437 bp using Roche-454 GS FLX Titanium sequencing technologies. Base composition of the S. suchowensis mt genome is A (27.43%), T (27.59%), C (22.34%), and G (22.64%), which shows a prevalent GC content with that of other angiosperms. This long circular mt genome encodes 58 unique genes (32 protein-coding genes, 23 tRNA genes and 3 rRNA genes), and 9 of the 32 protein-coding genes contain 17 introns. Through the phylogenetic analysis of 35 species based on 23 protein-coding genes, it is supported that Salix as a sister to Populus. With the detailed phylogenetic information and the identification of phylogenetic position, some ribosomal protein genes and succinate dehydrogenase genes are found usually lost during evolution. As a native shrub willow species, this worthwhile research of S. suchowensis mt genome will provide more desirable information for better understanding the genomic breeding and missing pieces of sex determination evolution in the future.

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.

Confirmation of Single-Locus Sex Determination and Female Heterogamety in Willow Based on Linkage Analysis

In this study, we constructed high-density genetic maps of Salix suchowensis and mapped the gender locus with an F₁ pedigree. Genetic maps were separately constructed for the maternal and paternal parents by using amplified fragment length polymorphism (AFLP) markers and the pseudo-testcross strategy. The maternal map consisted of 20 linkage groups that spanned a genetic distance of 2333.3 cM; whereas the paternal map contained 21 linkage groups that covered 2260 cM. Based on the established genetic maps, it was found that the gender of willow was determined by a single locus on linkage group LG_03, and the female was the heterogametic gender. Aligned with mapped SSR markers, linkage group LG_03 was found to be associated with chromosome XV in willow. It is noteworthy that marker density in the vicinity of the gender locus was significantly higher than that expected by chance alone, which indicates severe recombination suppression around the gender locus. In conclusion, this study confirmed the findings on the single-locus sex determination and female heterogamety in willow. It also provided additional evidence that validated the previous studies, which found that different autosomes evolved into sex chromosomes between the sister genera of Salix (willow) and Populus (poplar).

Molecular Breeding for Improved Second Generation Bioenergy Crops

There is increasing urgency to develop and deploy sustainable sources of energy to reduce our global dependency on finite, high-carbon fossil fuels. Lignocellulosic feedstocks, used in power and liquid fuel generation, are valuable sources of non-food plant biomass. They are cultivated with minimal inputs on marginal or degraded lands to prevent competition with arable agriculture and offer significant potential for sustainable intensification (the improvement of yield without the necessity for additional inputs) through advanced molecular breeding. This article explores progress made in next generation sequencing, advanced genotyping, association genetics, and genetic modification in second generation bioenergy production. Using poplar as an exemplar where most progress has been made, a suite of target traits is also identified giving insight into possible routes for crop improvement and deployment in the immediate future.

Transcriptome Profile Analysis from Different Sex Types of Ginkgo biloba L

In plants, sex determination is a comprehensive process of correlated events, which involves genes that are differentially and/or specifically expressed in distinct developmental phases. Exploring gene expression profiles from different sex types will contribute to fully understanding sex determination in plants. In this study, we conducted RNA-sequencing of female and male buds (FB and MB) as well as ovulate strobilus and staminate strobilus (OS and SS) of Ginkgo biloba to gain insights into the genes potentially related to sex determination in this species. Approximately 60 Gb of clean reads were obtained from eight cDNA libraries. De novo assembly of the clean reads generated 108,307 unigenes with an average length of 796 bp. Among these unigenes, 51,953 (47.97%) had at least one significant match with a gene sequence in the public databases searched. A total of 4709 and 9802 differentially expressed genes (DEGs) were identified in MB vs. FB and SS vs. OS, respectively. Genes involved in plant hormone signal and transduction as well as those encoding DNA methyltransferase were found to be differentially expressed between different sex types. Their potential roles in sex determination of G. biloba were discussed. Pistil-related genes were expressed in male buds while anther-specific genes were identified in female buds, suggesting that dioecism in G. biloba was resulted from the selective arrest of reproductive primordia. High correlation of expression level was found between the RNA-Seq and quantitative real-time PCR results. The transcriptome resources that we generated allowed us to characterize gene expression profiles and examine differential expression profiles, which provided foundations for identifying functional genes associated with sex determination in G. biloba.

Characterization of early transcriptional responses to cadmium in the root and leaf of Cd-resistant Salix matsudana Koidz

BACKGROUND

Salix matsudana Koidz. is a fast growing tree species. It has a high cadmium (Cd) tolerance capacity, making it potentially suitable for phytoremediation. Presently, transcriptomic and physiological Cd response mechanisms are poorly understood. Transcriptomic analysis in early response to high (50 μM) Cd levels was investigated in leaf and root of Cd-resistant S. matsudana Koidz..

RESULTS

Analysis of the response profiles demonstrate the existence of a complex transcriptional network in the root and leaf when exposed to Cd. The main response in the root involved up-regulation of genes associated with defence response via callose deposition in the cell wall and cell wall thickening. In the leaf, transcripts related to biotic stress signalling and secondary metabolism were activated. Additionally, many lignin and brassinosteroids synthesis pathway genes were induced mainly in the leaf, indicating that gene response to Cd was tissue-specific. The Cd transcriptome results were consistent with observed physiological changes.

CONCLUSION

The sub-localization, transcriptional network, and physiological regulation demonstrate the tissue-specific manner of Cd response, and provide a novel insight into in early response of tree species to Cd exposure.

Sex-Biased Temporal Gene Expression in Male and Female Floral Buds of Seabuckthorn (Hippophae rhamnoides)

Seabuckthorn is an economically important dioecious plant in which mechanism of sex determination is unknown. The study was conducted to identify seabuckthorn homologous genes involved in floral development which may have role in sex determination. Forty four putative Genes involved in sex determination (GISD) reported in model plants were shortlisted from literature survey, and twenty nine seabuckthorn homologous sequences were identified from available seabuckthorn genomic resources. Of these, 21 genes were found to differentially express in either male or female flower bud stages. HrCRY2 was significantly expressed in female flower buds only while HrCO had significant expression in male flowers only. Among the three male and female floral development stages (FDS), male stage II had significant expression of most of the GISD. Information on these sex-specific expressed genes will help in elucidating sex determination mechanism in seabuckthorn.

Transcriptome analysis of differentially expressed genes relevant to variegation in peach flowers

BACKGROUND

Variegation in flower color is commonly observed in many plant species and also occurs on ornamental peaches (Prunus persica f. versicolor [Sieb.] Voss). Variegated plants are highly valuable in the floricultural market. To gain a global perspective on genes differentially expressed in variegated peach flowers, we performed large-scale transcriptome sequencing of white and red petals separately collected from a variegated peach tree.

RESULTS

A total of 1,556,597 high-quality reads were obtained, with an average read length of 445 bp. The ESTs were assembled into 16,530 contigs and 42,050 singletons. The resulting unigenes covered about 60% of total predicted genes in the peach genome. These unigenes were further subjected to functional annotation and biochemical pathway analysis. Digital expression analysis identified a total of 514 genes differentially expressed between red and white flower petals. Since peach flower coloration is determined by the expression and regulation of structural genes relevant to flavonoid biosynthesis, a detailed examination detected four key structural genes, including C4H, CHS, CHI and F3H, expressed at a significantly higher level in red than in white petal. Except for the structural genes, we also detected 11 differentially expressed regulatory genes relating to flavonoid biosynthesis. Using the differentially expressed structural genes as the test objects, we validated the digital expression results by using quantitative real-time PCR, and the differential expression of C4H, CHS and F3H were confirmed.

CONCLUSION

In this study, we generated a large EST collection from flower petals of a variegated peach. By digital expression analysis, we identified an informative list of candidate genes associated with variegation in peach flowers, which offered a unique opportunity to uncover the genetic mechanisms underlying flower color variegation.

De novo assembly and characterization of Sophora japonica transcriptome using RNA-seq

Sophora japonica Linn (Chinese Scholar Tree) is a shrub species belonging to the subfamily Faboideae of the pea family Fabaceae. In this study, RNA sequencing of S. japonica transcriptome was performed to produce large expression datasets for functional genomic analysis. Approximate 86.1 million high-quality clean reads were generated and assembled de novo into 143010 unique transcripts and 57614 unigenes. The average length of unigenes was 901 bps with an N50 of 545 bps. Four public databases, including the NCBI nonredundant protein (NR), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Cluster of Orthologous Groups (COG), were used to annotate unigenes through NCBI BLAST procedure. A total of 27541 of 57614 unigenes (47.8%) were annotated for gene descriptions, conserved protein domains, or gene ontology. Moreover, an interaction network of unigenes in S. japonica was predicted based on known protein-protein interactions of putative orthologs of well-studied plant genomes. The transcriptome data of S. japonica reported here represents first genome-scale investigation of gene expressions in Faboideae plants. We expect that our study will provide a useful resource for further studies on gene expression, genomics, functional genomics, and protein-protein interaction in S. japonica.

Transcriptome comparison reveals key candidate genes responsible for the unusual reblooming trait in tree peonies

Tree peonies are important ornamental plants worldwide, but growing them can be frustrating due to their short and concentrated flowering period. Certain cultivars exhibit a reblooming trait that provides a valuable alternative for extending the flowering period. However, the genetic control of reblooming in tree peonies is not well understood. In this study, we compared the molecular properties and morphology of reblooming and non-reblooming tree peonies during the floral initiation and developmental processes. Using transcriptome sequencing technology, we generated 59,275 and 63,962 unigenes with a mean size of 698 bp and 699 bp from the two types of tree peonies, respectively, and identified eight differentially expressed genes that are involved in the floral pathways of Arabidopsis thaliana. These differentially regulated genes were verified through a detailed analysis of their expression pattern during the floral process by real time RT-PCR. From this combined analysis, we identified four genes, PsFT, PsVIN3, PsCO and PsGA20OX, which likely play important roles in the regulation of the reblooming process in tree peonies. These data constitute a valuable resource for the discovery of genes involved in flowering time and insights into the molecular mechanism of flowering to further accelerate the breeding of tree peonies and other perennial woody plants.