Comparative transcript profiling explores differentially expressed genes associated with sexual phenotype in kiwifruit

Genome-wide identification, characterization and gene expression of BES1 transcription factor family in grapevine (Vitis vinifera L.)

BES1, as the most important transcription factor responsible for brassinolide (BR) signaling, has been confirmed to play a significant role in regulating plant growth and the improvement of stress resistance. The transcriptional regulatory mechanism of BES1 has been well elucidated in several plants, such as Arabidopsis thaliana (A. thaliana), Triticum aestivum L. (T. aestivum), and Oryza sativa L. (O. sativa). Nevertheless, the genome-wide analysis of the BES1 family in Vitis vinifera L. (V. vinifera). has not been comprehensively carried out. Thus, we have conducted a detailed analysis and identification of the BES1 transcription factors family in V. vinifera; a total of eight VvBES1 genes was predicted, and the phylogenetic relationships, gene structures, and Cis-acting element in their promoters were also analyzed. BES1 genes have been divided into three groups (I, II and III) based on phylogenetic relationship analysis, and most of VvBES1 genes were in group III. Also, we found that VvBES1 genes was located at seven of the total nineteen chromosomes, whereas VvBES1-2 (Vitvi04g01234) and VvBES1-5 (Vitvi18g00924) had a collinearity relationship, and their three copies are well preserved. In addition, the intron-exon model of VvBES1 genes were mostly conserved, and there existed several Cis-acting elements related to stress resistance responsive and phytohormones responsive in BES1s genes promoter. Moreover, the BES1 expressions were different in different V. vinifera organs, and BES1 expressions were different in different V. vinifera varieties under saline-alkali stress and heat stress, the expression of VvBES1 also changed with the prolongation of saline-alkali stress treatment time. The above findings could not only lay a primary foundation for the further validation of VvBES1 function, but could also provide a reference for molecular breeding in V. vinifera.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

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.

Transcriptomic Analysis of Sex-Associated DEGs in Female and Male Flowers of Kiwifruit (Actinidia deliciosa [A. Chev] C. F. Liang & A. R. Ferguson)

Kiwifruit (Actinidia deliciosa [A. Chev.], C.V. Liang & A. R. Ferguson, 1984) is a perennial plant, with morphologically hermaphroditic and functionally dioecious flowers. Fruits of this species are berries of great commercial and nutritional importance. Nevertheless, few studies have analyzed the molecular mechanisms involved in sexual differentiation in this species. To determine these mechanisms, we performed RNA-seq in floral tissue at stage 60 on the BBCH scale in cultivar ‘Hayward’ (H, female) and a seedling from ‘Green Light’ × ‘Tomuri’ (G × T, male). From these analyses, we obtained expression profiles of 24,888 (H) and 27,027 (G × T) genes, of which 6413 showed differential transcript abundance. Genetic ontology (GO) and KEGG analysis revealed activation of pathways associated with the translation of hormonal signals, plant-pathogen interaction, metabolism of hormones, sugars, and nucleotides. The analysis of the protein-protein interaction network showed that the genes ERL1, AG, AGL8, LFY, WUS, AP2, WRKY, and CO, are crucial elements in the regulation of the hormonal response for the formation and development of anatomical reproductive structures and gametophytes. On the other hand, genes encoding four Putative S-adenosyl-L-methionine-dependent methyltransferases (Achn201401, Achn281971, Achn047771 and Achn231981) were identified, which were up-regulated mainly in the male flowers. Moreover, the expression profiles of 15 selected genes through RT-qPCR were consistent with the results of RNA-seq. Finally, this work provides gene expression-based interactions between transcription factors and effector genes from hormonal signaling pathways, development of floral organs, biological and metabolic processes or even epigenetic mechanisms which could be involved in the kiwi sex-determination. Thus, in order to decode the nature of these interactions, it could be helpful to propose new models of flower development and sex determination in the Actinidia genus.

Induction, identification and genetics analysis of tetraploid Actinidia chinensis

Actinidia chinensis is a commercially important fruit, and tetraploid breeding of A. chinensis is of great significance for economic benefit. In order to obtain elite tetraploid cultivars, tetraploid plants were induced by colchicine treatment with leaves of diploid A. chinensis 'SWFU03'. The results showed that the best treatment was dipping leaves 30 h in 60 mg l^-1 colchicine solutions, with induction rate reaching 26%. Four methods, including external morphology comparison, stomatal guard cell observation, chromosome number observation and flow cytometry analysis were used to identify the tetraploid of A. chinensis. Using the induction system and flow cytometry analysis methods, 187 tetraploid plants were identified. Three randomly selected tetraploid plants and their starting diploid plants were further subjected to transcriptome analysis, real-time quantitative polymerase chain reaction (RT-qPCR) and methylation-sensitive amplification polymorphism (MSAP) analysis. The transcriptome analysis results showed that there were a total of 2230 differentially expressed genes (DEG) between the diploid and tetraploid plants, of which 660 were downregulated and 1570 upregulated. The DEGs were mainly the genes involved in growth and development, stress resistance and antibacterial ability in plants. RT-qPCR results showed that the gene expression levels of the growth and stress resistance of tetraploid plants were higher than those of diploid ones at the transcriptome level. MSAP analysis of DNA methylation results showed that tetraploid plants had lower methylation ratio than diploid ones. The present results were valuable to further explore the epigenetics of diploid and tetraploid kiwifruit plants.

Sex change in kiwifruit (Actinidia chinensis Planch.): a developmental framework for the bisexual to unisexual floral transition

The developmental morphology of male and female kiwifruit flowers is tracked to delimit a framework of events to aid the study of divergence in floral gene expression. The transition from hermaphrodite to unisexual development of kiwifruit (Actinidia chinensis Planch) flowers has been reported previously, but differences in gene expression controlling sexual development for this species have not been associated with the major developmental changes occurring within pistils. We investigated the key stages in male and female flower development to define the point at which meristematic activities diverge in the two sexes. A combination of scanning electron microscopy and light microscopy was used to investigate pistil development from the earliest stages. We identified seven distinct stages characterized by differences in ovary size and shape, macrosporogenesis, ovule primordium development, anther locule lengthening, microspore wall thickening, and pollen degeneration. Sex differences were evident from the initial stage of development, with a laterally compacted gynoecium in male flowers. However, the key developmental stage, at which tissue differentiation clearly deviated between the two sexes, was stage 3, when flowers were 3.5 to 4.5 mm in length at approximately 10 d from initiation of stamen development. At this stage, male flowers lacked evident carpel meristem development as denoted by a lack of ovule primordium formation. Pollen degeneration in female flowers, probably driven by programmed cell death, occurred at the late stage 6, while the final stage 7 was represented by pollen release. As the seven developmental stages are associated with specific morphological differences, including flower size, the scheme suggested here can provide the required framework for the future study of gene expression during the regulation of flower development in this crop 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.