Transcriptome profiling of differentially expressed genes in floral buds and flowers of male sterile and fertile lines in watermelon

Identification of differentially expressed miRNAs between male sterile and fertile floral buds in watermelon (Citrullus lanatus L.) via high-throughput sequencing

This experiment used floral buds from watermelon genic male sterile dual-purpose lines as materials to explore the differentially expressed miRNAs (DEMs) between male fertile and sterile floral buds of watermelon. Paraffin sectioning technology was employed for a cytological analysis, and small RNA sequencing was used to explore miRNAs related to anther or pollen development. Cytological analysis indicated that the abnormal development of tapetal cells may cause microspore abortion. Small RNA sequencing identified a total of 314 miRNAs (29 known and 285 novel, which belonged to 12 and 61 miRNA families, respectively) in floral buds. Differential expression revealed 36 (5 known and 31 novel) DEMs between male fertile and sterile buds, 7 and 29 of which were up-regulated and down-regulated, respectively. Target genes analysis showed that the 36 DEMs were predicted to target 577 genes, and these targets might participate in various biological processes, such as response to metal ions, floral organ development, stamen development, anther development, pollen maturation, and programmed cell death. Moreover, pathway analysis indicated that these genes were mainly enriched in purine metabolism, starch and sucrose metabolism, RNA transport, and other pathways. In addition, the 55 miRNA-target modules, including 3 known and 16 novel miRNAs with 30 target genes, might be related to anther or pollen development in watermelon. Our findings provide important miRNA-target modules related to watermelon anther or pollen development and can lay the foundation for biological functional analysis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04084-6.

Comparative transcriptome analysis reveals the involvement of an MYB transcriptional activator, SmMYB108, in anther dehiscence in eggplant

Male sterility is a highly attractive agronomic trait as it effectively prevents self-fertilization and facilitates the production of high-quality hybrid seeds in plants. Timely release of mature pollen following anther dehiscence is essential for stamen development in flowering plants. Although several theories have been proposed regarding this, the specific mechanism of anther development in eggplant remains elusive. In this study, we selected an R2R3-MYB transcription factor gene, SmMYB108, that encodes a protein localized primarily in the nucleus by comparing the transcriptomics of different floral bud developmental stages of the eggplant fertile line, F142. Quantitative reverse transcription polymerase chain reaction revealed that SmMYB108 was preferentially expressed in flowers, and its expression increased significantly on the day of flowering. Overexpression of SmMYB108 in tobacco caused anther dehiscence. In addition, we found that SmMYB108 primarily functions as a transcriptional activator via C-terminal activation (amino acid 262-317). Yeast one-hybrid and dual-luciferase reporter assays revealed that genes (SmMYB21, SmARF6, and SmARF8) related to anther development targeted the SmMYB108 promoter. Overall, our results provide insights into the molecular mechanisms involved in the regulation of anther development by SmMYB108.

Virus-induced gene silencing for in planta validation of gene function in cucurbits

Virus-induced gene silencing (VIGS) is a powerful tool for high-throughput analysis of gene function. Here, we developed the VIGS vector pCF93, from which expression of the cucumber fruit mottle mosaic virus genome is driven by the cauliflower mosaic virus 35S promoter to produce viral transcripts in inoculated plants. To test the utility of the pCF93 vector, we identified candidate genes related to male sterility (MS) in watermelon (Citrullus lanatus), which is recalcitrant to genetic transformation. Specifically, we exploited previously reported reference-based and de novo transcriptome data to define 38 differentially expressed genes between a male-sterile line and its fertile near-isogenic line in the watermelon cultivar DAH. We amplified 200- to 300-bp fragments of these genes, cloned them into pCF93, and inoculated DAH with the resulting VIGS clones. The small watermelon cultivar DAH enabled high-throughput screening using a small cultivation area. We simultaneously characterized the phenotypes associated with each of the 38 candidate genes in plants grown in a greenhouse. Silencing of 8 of the 38 candidate genes produced male-sterile flowers with abnormal stamens and no pollen. We confirmed the extent of gene silencing in inoculated flowers using reverse transcription-qPCR. Histological analysis of stamens from male-fertile and male-sterile floral buds and mature flowers revealed developmental defects and shrunken pollen sacs. Based on these findings, we propose that the pCF93 vector and our VIGS system will facilitate high-throughput analysis for the study of gene function in watermelons.

Disruption of REC8 in Meiosis I led to watermelon seedless

In triploid watermelon (Citrullus lanatus), the homologous chromosomes of germ cells are disorder during meiosis, resulting in the failure of seeds formation and producing seedless fruit. Therefore, mutating the genes specifically functioning in meiosis may be an alternative way to achieve seedless watermelon. REC8, as a key component of the cohesin complex in meiosis, is dramatically essential for sister chromatid cohesion and chromosome segregation. However, the role of REC8 in meiosis has not yet been characterized in watermelon. Here, we identified ClREC8 as a member of RAD21/REC8 family with a high expression in male and female flowers of watermelon. In situ hybridization analysis showed that ClREC8 was highly expressed at the early stage of meiosis during pollen formation. Knocking out ClREC8 in watermelon led to decline of pollen vitality. After pollinating with foreign normal pollen, the ovaries of ClREC8 knockout lines could inflate normally but failed to form seeds. We further compared the meiosis chromosomes of pollen mother cells in different stages between the knockout lines and the corresponding wild type. The results indicated that ClREC8 was required for the monopolar orientation of the sister kinetochores in Meiosis I. Additionally, transcriptome sequencing (RNA-seq) analysis between WT and the knockout lines revealed that the disruption of ClREC8 caused the expression levels of mitosis-related genes and meiosis-related genes to decrease. Our results demonstrated ClREC8 has a specific role in Meiosis I of watermelon germ cells, and loss-of-function of the ClREC8 led to seedless fruit, which may provide an alternative strategy to breed cultivars with seedless watermelon.

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.

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.

Disruption of the bHLH transcription factor Abnormal Tapetum 1 causes male sterility in watermelon

Although male sterility has been identified as a useful trait for hybrid vigor utilization and hybrid seed production, its underlying molecular mechanisms in Cucurbitaceae species are still largely unclear. Here, a spontaneous male-sterile watermelon mutant, Se18, was reported to have abnormal tapetum development, which resulted in completely aborted pollen grains. Map-based cloning demonstrated that the causal gene Citrullus lanatus Abnormal Tapetum 1 (ClATM1) encodes a basic helix-loop-helix (bHLH) transcription factor with a 10-bp deletion and produces a truncated protein without the bHLH interaction and functional (BIF) domain in Se18 plants. qRT-PCR and RNA in situ hybridization showed that ClATM1 is specifically expressed in the tapetum layer and in microsporocytes during stages 6-8a of anther development. The genetic function of ClATM1 in regulating anther development was verified by CRISPR/Cas9-mediated mutagenesis. Moreover, ClATM1 was significantly downregulated in the Se18 mutant, displaying a clear dose effect at the transcriptional level. Subsequent dual-luciferase reporter, β-glucuronidase (GUS) activity, and yeast one-hybrid assays indicated that ClATM1 could activate its own transcriptional expression through promoter binding. Collectively, ClATM1 is the first male sterility gene cloned from watermelon, and its self-regulatory activity provides new insights into the molecular mechanism underlying anther development in plants.

Global gene expression analysis of pigeonpea with male sterility conditioned by A2 cytoplasm

Cytoplasmic male sterility(CMS), a maternally inherited trait, provides a promising means to harness yield gains associated with hybrid vigor. In pigeonpea [Cajanus cajan (L.) Huth], nine types of sterility-inducing cytoplasm have been reported, of which A₂ and A₄ have been successfully deployed in hybrid breeding. Unfortunately, molecular mechanism of the CMS trait is poorly understood because of limited research invested. More recently, an association between a mitochondrial gene (nad7) and A₄ -CMS has been demonstrated in pigeonpea; however, the mechanism underlying A₂ -CMS still remains obscure. The current investigation aimed to analyze the differences in A₂ -CMS line (ICPL 88039A) and its isogenic maintainer line (ICPL 88039B) at transcriptome level using next-generation sequencing. Gene expression profiling uncovered a set of 505 genes that showed altered expression in response to CMS, of which, 412 genes were upregulated while 93 were downregulated in the fertile maintainer line vs. the CMS line. Further, gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein-protein interaction (PPI) network analyses revealed association of CMS in pigeonpea with four major pathways: glucose and lipid metabolism, ATP production, pollen development and pollen tube growth, and reactive oxygen species (ROS) scavenging. Patterns of digital gene expression were confirmed by quantitative real-time polymerase chain reaction (qRT-PCR) of six candidate genes. This study elucidates candidate genes and metabolic pathways having potential associations with pollen development and male sterility in pigeonpea A₂ -CMS. New insights on molecular mechanism of CMS trait in pigeonpea will be helpful to accelerate heterosis utilization for enhancing productivity gains in pigeonpea.

Mapping and Analysis of a Novel Genic Male Sterility Gene in Watermelon (Citrullus lanatus)

Seed production is critical for watermelon production, which mostly involves first-generation hybrid varieties. However, watermelon hybrid seed production currently requires complex procedures, including artificial isolation and pollination. Therefore, the development and use of a male-sterile system to generate watermelon hybrids can simplify the process. The scarcity of male-sterile watermelon germplasm resources necessitates the use of molecular breeding methods. Unfortunately, the genes responsible for male sterility in watermelon have not been cloned. Thus, the genetic basis of the male sterility remains unknown. In this study, two DNA pools derived from male-sterile and normal plants in the F2 population were used for whole-genome resequencing. The Illumina high-throughput sequencing resulted in 62.99 Gbp clean reads, with a Q30 of 80% after filtering. On the basis of the SNP index association algorithm, eight candidate regions (0.32 Mb) related to specific traits were detected on chromosome 6. Expression pattern analyses and watermelon transformation studies generated preliminary evidence that Cla006625 encodes a pollen-specific leucine-rich repeat protein (ClaPEX1) influencing the male sterility of watermelon. The identification and use of genic male sterility genes will promote watermelon male sterility research and lay the foundation for the efficient application of seed production technology.

Comparative transcriptome between male fertile and male sterile alfalfa (Medicago varia)

Male sterility is an important factor in improving crop quality and yield through heterosis breeding. In this study, we analyzed the transcriptomes of male fertile (MF) and male sterile (MS) alfalfa flower buds using the Illumina HiSeq™ 4000 platform. A total of 54.05 million clean reads were generated and assembled into 65,777 unigenes with an average length of 874 bp. The differentially expressed genes (DEGs) between the MF and MS flowers at three stages of pollen development were identified, and there were 3832, 5678 and 5925 DEGs respectively in stages 1, 2 and 3. GO and KEGG functional enrichment analysis revealed 12, 12, 6 and 12 key branch-point genes involved in circadian rhythm, transcription factors, pollen development and flavonoid biosynthesis. Our findings provide novel insights into the mechanism of male sterility in alfalfa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01026-x.

The impaired biosynthetic networks in defective tapetum lead to male sterility in watermelon

Heterosis has been widely applied in watermelon breeding, because of the higher resistance and yield of hybrid. As the basis of heterosis utilization, genic male sterility (GMS) is an important tool for facilitating hybrid seed production, while the detailed mechanism in watermelon is still largely unknown. Here, we report a spontaneous mutant Se18 exhibited complete male sterility due to the uniquely multilayered tapetum and the un-meiotic pollen mother cells during pollen development. Using TMT based quantitative proteomic analyses, a total of 348 differentially abundant proteins (DAPs) were detected with the overwhelming majority down-regulated in mutant Se18. By analyzing the putative orthologs/homologs of Arabidopsis GMS related genes, the biosynthesis and transport of sporopollenin and tryphine precursors were predictably altered in mutant compared to its sibling wild type. Moreover, the general phenylpropanoid pathway as well as its related metabolisms was also expectably impaired in mutant, coincident with the pale yellow petals. Notably, some key transcriptional factors regulating tapetum development, together with their down-regulated targets, offered potentially valuable candidates regarding of male sterility. Collectively, the disrupted regulatory networks underlying male sterility of watermelon was proposed, which provide novel insights into genetic mechanism of male reproductive process and rich gene resources for future research. SIGNIFICANCE: Watermelon is an importantly economical cucurbit crop worldwide, with high nutritional value. Although several male sterile mutants have been identified in watermelon, the underlying molecular mechanism is poorly elucidated. Comparative cytological analysis revealed that the defective development of tapetum was responsible for male sterility in mutant Se18. Combined with the morphological comparison, male floral buds at 2.0-2.5 mm in diameter were confirmed with no obvious phenotypic differences but distinct cytological defects, which were in turn sampled for TMT based proteomic analyses. Referring to functionally characterized GMS related genes, the genetic pathway DYT1-TDF1-AMS-MS188-MS1 regulating tapetum development, together with some downstream targets, were considerably altered in mutant Se18. Moreover, enrichment analyses illustrated the general phenylpropanoid related metabolisms, as well as the biosynthesis and transport of sporopollenin and tryphine precursors, were significantly disrupted in defective anther development. Collectively, the proposed regulatory networks in watermelon not only contribute to a better understanding of molecular mechanisms underlying male sterility, but also provide valuable GMS related candidates for future researches.

RNA Sequencing Reveals Rice Genes Involved in Male Reproductive Development under Temperature Alteration

Rice (Oryza sativa L.) is one of the most important food crops, providing food for nearly half of the world population. Rice grain yields are affected by temperature changes. Temperature stresses, both low and high, affect male reproductive development, resulting in yield reduction. Thermosensitive genic male sterility (TGMS) rice is sterile at high temperature and fertile at low temperature conditions, facilitating hybrid production, and is a good model to study effects of temperatures on male development. Semithin sections of the anthers of a TGMS rice line under low (fertile) and high (sterile) temperature conditions showed differences starting from the dyad stage, suggesting that genes involved in male development play a role during postmeiotic microspore development. Using RNA sequencing (RNA-Seq), transcriptional profiling of TGMS rice panicles at the dyad stage revealed 232 genes showing differential expression (DEGs) in a sterile, compared to a fertile, condition. Using qRT-PCR to study expression of 20 selected DEGs using panicles of TGMS and wild type rice plants grown under low and high temperature conditions, revealed that six out of the 20 selected genes may be unique to TGMS, while the other 14 genes showed common responses to temperatures in both TGMS and wild-type rice plants. The results presented here would be useful for further investigation into molecular mechanisms controlling TGMS and rice responses to temperature alteration.

Full-length transcriptome analysis reveals the differences between floral buds of recessive genic male-sterile line (RMS3185A) and fertile line (RMS3185B) of cabbage

Starch and sucrose metabolism and plant-pathogen interaction pathways play a dominate role in recessive genic male sterility (RGMS) of cabbage (Brassica oleracea L. var. capitata). RGMS is common in plants and has been widely applied as an effective and economic system for hybrid seed production in many crops. However, little is known regarding the molecular mechanisms of RGMS in cabbage. Hence, full-length transcriptomic and physiological analysis were performed in the spontaneous RGMS mutant RMS3185A and its near-isogenic fertile line (NIL) RMS3185B of small (< 1.6 mm in diameter), medium (~ 2.5 mm in diameter), and large floral buds (~ 3.4 mm in diameter) to identify the differentially expressed genes (DEGs) associated with RMGS. The pollen abnormalities between RMS3185B and RMS3185A appeared at the large floral bud stage. In contrast with RMS3185B, the mature anthers and stamens of RMS3185A were shorter than those of RMS3185B, and the anthers did not dehiscent. The concentrations of glucose, fructose, trehalose, starch, and cellulose in RMS3185A were all significantly lower than those in large floral buds of RMS3185B. PacBio sequencing results showed that DEGs were mainly concentrated in large floral bud stage. In combination with the KEGG enrichment analysis of DEGs in GO terms related to cell wall, pollen and anther, pentose and glucuronate interconversions (ko00040), starch and sucrose metabolism (ko00500), and plant-pathogen interaction (ko04626) were significantly enriched. Among which, cell-wall/pectin-related genes of eighteen PEI, twenty-two PEL, three PG, and fifteen PGL involved in ko00040, and one UGDH, one SPS, four CWINV, four TPP/TPS, and four EGL involved in ko00500, as well as plant-pathogen interaction genes, including sixteen calcium-dependent protein kinase (CDPK), one cyclic nucleotide-gated ion channel (CNGC), and twenty-three calcium-binding protein CML (CML), were significantly down-regulated in RMS3185A relative to that in RMS3185B. Besides, genes involved in ko04626, including two CML and one transcription factor WRKY33, were up-regulated in RMS3185A relative to that in RMS3185B. In conclusion, we hypothesized that the expression alterations of these genes were responsible for calcium signaling and sugar metabolism, thus affecting the occurrence of RGMS in cabbage.

Primary Cell Wall Modifying Proteins Regulate Wall Mechanics to Steer Plant Morphogenesis

Plant morphogenesis involves multiple biochemical and physical processes inside the cell wall. With the continuous progress in biomechanics field, extensive studies have elucidated that mechanical forces may be the most direct physical signals that control the morphology of cells and organs. The extensibility of the cell wall is the main restrictive parameter of cell expansion. The control of cell wall mechanical properties largely determines plant cell morphogenesis. Here, we summarize how cell wall modifying proteins modulate the mechanical properties of cell walls and consequently influence plant morphogenesis.

Analysis of differentially expressed genes and pathways associated with male sterility lines in watermelon via bulked segregant RNA-seq

Genic male sterility (GMS) is a common and important trait, which is widely used for the production of hybrid seeds. However, the molecular mechanism of GMS in watermelon remains poorly understood. In this study, we comparatively analyzed the transcriptome profiles of sterile and fertile floral buds using the bulked segregant analysis (BSA) and transcriptome sequencing (RNA-seq). A total of 2507 differentially expressed genes (DEGs) including 593 up-regulated and 1914 down-regulated, were identified to be related to male sterility in watermelon line Se18. Gene ontology (GO) analysis showed that 57 GO terms were significantly enriched, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed plant hormone signal transduction, glycolysis/gluconeogenesis, starch and sucrose metabolism, plant-pathogen interaction, phenylpropanoid biosynthesis pathways were obviously enriched. Furthermore, the efficiency of the RNA-seq analysis was validated by quantitative real-time PCR (qRT-PCR). Among the DEGs, some valuable candidate genes involved in pollen development were identified, such as gene Cla000029, a bHLH transcription factor and homologous to MS1 in Arabidopsis. Moreover, other DEGs including MYB gene Cla012590 (MYB26), Cla017100 (MYB21), etc., also provide useful information for further understanding the function of key genes involved in pollen development. This study provides new insights into the global network of male sterility in watermelon.

Transcriptome profiling of differentially expressed genes in cytoplasmic male-sterile line and its fertility restorer line in pigeon pea (Cajanus cajan L.)

BACKGROUND

Pigeon pea (Cajanus cajan L.) is the sixth major legume crop widely cultivated in the Indian sub-continent, Africa, and South-east Asia. Cytoplasmic male-sterility (CMS) is the incompetence of flowering plants to produce viable pollens during anther development. CMS has been extensively utilized for commercial hybrid seeds production in pigeon pea. However, the molecular basis governing CMS in pigeon pea remains unclear and undetermined. In this study transcriptome analysis for exploring differentially expressed genes (DEGs) between cytoplasmic male-sterile line (AKCMS11) and its fertility restorer line (AKPR303) was performed using Illumina paired-end sequencing.

RESULTS

A total of 3167 DEGs were identified, of which 1432 were up-regulated and 1390 were down-regulated in AKCMS11 in comparison to AKPR303. By querying, all the 3167 DEGs against TAIR database, 34 pigeon pea homologous genes were identified, few involved in pollen development (EMS1, MS1, ARF17) and encoding MYB and bHLH transcription factors with lower expression in the sterile buds, implying their possible role in pollen sterility. Many of these DEGs implicated in carbon metabolism, tricarboxylic acid cycle (TCA), oxidative phosphorylation and elimination of reactive oxygen species (ROS) showed reduced expression in the AKCMS11 (sterile) buds.

CONCLUSION

The comparative transcriptome findings suggest the potential role of these DEGs in pollen development or abortion, pointing towards their involvement in cytoplasmic male-sterility in pigeon pea. The candidate DEGs identified in this investigation will be highly significant for further research, as they could lend a comprehensive basis in unravelling the molecular mechanism governing CMS in pigeon pea.

Energy metabolism involved in fertility of the wheat TCMS line YS3038

In the wheat TCMS line YS3038, the anther development is inhibited from late uninucleate stage to the binucleate stage. The disruption of energy metabolism pathways by aberrant transcriptional regulation causes the male sterility under low temperatures. The utilization of thermosensitive male sterile (TMS) lines provides a basis for two-line breeding. Previous work, including morphological and cytological observations, has shown that the development process of the TMS line YS3038 is inhibited from the late uninucleate stage to the binucleate stage. Transcriptomics studies could now help to elucidate the overall expression of related genes in a specific reproductive process, revealing the metabolic network and its regulatory mechanism of the reproductive process from the transcription level. Considering the fertility characteristics of YS3038, three important stages for transcriptome analysis were determined to be the early uninucleate, late uninucleate and binucleate stages. The number of differentially expressed genes (DEGs) was found to be highest in the binucleate stage, and most were related to energy metabolism. Quantitative PCR analysis of selected genes related to energy metabolism revealed that their expression patterns were consistent with the sequencing results. Analysis of the fertility mechanism of YS3038 showed that although the tapetum of anthers was degraded in advance of the tetrad stage, the development of microspores did not result in obvious abnormalities until the binucleate stage, because the genes involved in energy metabolism pathways, including starch and sucrose metabolism (SSM), glycolysis, the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and respiration electron transport chain are differentially expressed under sterile and fertile conditions. Therefore, the pollen in YS3038 was sterile.

Pectate lyases: Their role in plants and importance in fruit ripening

Plant cell walls are complex structures that are modified throughout development. They are a major contributor to the properties of plant structure and act as barriers against pathogens. The primary cell walls of plants are composed of polysaccharides and proteins. The polysaccharide fraction is divided into components cellulose, hemicelluloses and pectin, are all modified during fruit ripening. Pectin plays an important role in intercellular adhesion and controlling the porosity of the wall. A large number of pectin degrading enzymes have been characterised from plants and they are involved in numerous aspects of plant development. The role of pectate lyases in plant development has received little attention, probably because they are normally associated with the action of plant pathogenic organisms. However their importance in plant development and ripening is now becoming well established and new information about the role of pectate lyases in plant development forms the focus of this review.

Current understanding of male sterility systems in vegetable Brassicas and their exploitation in hybrid breeding

Overview of the current status of GMS and CMS systems available in Brassica vegetables, their molecular mechanism, wild sources of sterile cytoplasm and exploitation of male sterility in hybrid breeding. The predominantly herbaceous family Brassicaceae (crucifers or mustard family) encompasses over 3700 species, and many of them are scientifically and economically important. The genus Brassica is an economically important genus within the tribe Brassicaceae that comprises important vegetable, oilseed and fodder crops. Brassica vegetables display strong hybrid vigor, and heterosis breeding is the integral part in their improvement. Commercial production of F₁ hybrid seeds in Brassica vegetables requires an effective male sterility system. Among the available male sterility systems, cytoplasmic male sterility (CMS) is the most widely exploited in Brassica vegetables. This system is maternally inherited and studied intensively. A limited number of reports about the genic male sterility (GMS) are available in Brassica vegetables. The GMS system is reported to be dominant, recessive and trirecessive in nature in different species. In this review, we discuss the available male sterility systems in Brassica vegetables and their potential use in hybrid breeding. The molecular mechanism of mt-CMS and causal mitochondrial genes of CMS has been discussed in detail. Finally, the exploitation of male sterility system in heterosis breeding of Brassica vegetables, future prospects and need for further understanding of these systems are highlighted.

Comparative transcriptome analysis reveals new molecular pathways for cucumber genes related to sex determination

Transcriptome data and qPCR analysis revealed new insight into genes regulatory mechanism related to cucumber sex determination. Cucumber (Cucumis sativus L.) is an economically important crop cultivated worldwide. Enhancing the genomic resources for cucumber may enable the regulation of traits relevant to crop productivity and quality. Sequencing technologies and bioinformatics tools provide opportunities for the development of such resources. The aims of this study were to identify and characterize the genes involved in sex determination and flower morphogenesis in cucumber isogenic lines that differed regarding flower sex type. We obtained transcripts for 933 genes related to shoot apex development, among which 310 were differentially expressed genes (DEGs) among the male, female, and hermaphroditic lines. We performed gene ontology and molecular network analyses and explored the DEGs related to already known processes like: hormone synthesis and signaling, lipid and sugar metabolism; and also newly discovered processes related to cell wall, membrane, and cytoskeleton modifications; ion homeostasis which appears to be important for ethylene perception and signaling, and genes expression mediated by transcription factors related to floral organ identities. We proposed a new model of regulatory mechanism network of sex development in cucumber. Our results may be useful for clarifying the molecular genetics and the functional mechanisms underlying the sex determination processes.

An evolutionarily conserved non-synonymous SNP in a leucine-rich repeat domain determines anthracnose resistance in watermelon

A non-synonymous SNP of CC-NBS-LRR was firstly mapped to confer resistance to anthracnose in watermelon. Newly proposed LRR domain harboring the SNP is evolutionary conserved in the Cucurbitaceae and Fabaceae. Anthracnose disease caused by Colletotrichum devastates many plants. Despite the importance of the disease, the mechanisms of resistance against it are poorly understood. Here, we identified a non-synonymous single-nucleotide polymorphism (SNP) located in a leucine-rich repeat domain as a marker for resistance to anthracnose race 1 in watermelon, using a combination of genetic analyses. We validated this SNP in segregating populations and 59 watermelon accessions using high-resolution melting assays and Sanger sequencing. We demonstrated that the resulting arginine-to-lysine substitution is particularly conserved among the Cucurbitaceae and Fabaceae. We identified a conserved motif, IxxLPxSxxxLYNLQTLxL, found in 1007 orthologues/paralogues from 89 plant species, and discovered that residue 18 of this motif could determine resistance to disease caused by external invaders. This study provides a step forward in understanding anthracnose resistance in watermelon, as well as functional and evolutionary insight into leucine-rich repeat proteins.

Seedless mutant 'Wuzi Ougan' (Citrus suavissima Hort. ex Tanaka 'seedless') and the wild type were compared by iTRAQ-based quantitative proteomics and integratedly analyzed with transcriptome to improve understanding of male sterility

Background

Bud mutation is a vital method of citrus. ‘Wuzi Ougan’ (mutant type, MT) as a bud variant of ‘Ougan’ (wild type, WT) was first found in 1996 and has become popular because of its male sterility and seedless character. Previous analysis of its cytological sections and transcriptome revealed that the abnormal microsporogenesis that occurs before the tetrad stage of anther development might be the result of down-regulated oxidation-reduction biological processes in MT. To reveal the mechanism behind the male sterility in MT at the post-transcriptional stage, proteome profiling and integrative analysis on previously obtained transcriptome and proteome data were performed in two strains.

Results

The proteome profiling was performed by iTRAQ (isobaric Tags for relative and absolute quantitation) analysis and 6201 high-confidence proteins were identified, among which there were 487 differentially expressed proteins (DEPs) in one or more developmental stages of anthers between MT and WT. The main functional subcategories associated with the main category biological process into which the DEPs were classified were sporopollenin biosynthesis process and pollen exine formation. The enriched pathways were phenylpropanoid biosynthesis, flavonoid biosynthesis, and phenylalanine metabolism. Moreover, there were eight pathways linked in terms of being related to phenylpropanoid metabolism. Eighteen important genes related to phenylpropanoid metabolism were also analysized by qRT-PCR (quantitative real time PCR). An integrative analysis of the fold change at the transcript (log2 FPKM ratios) and protein (log1.2 iTRAQ ratios) levels was performed to reveal the consistency of gene expression at transcriptional and proteomic level. In general, the expression of genes and proteins tended to be positively correlated, in which the correlation coefficients were 0.3414 (all genes and all proteins) and 0.5686 (DEPs and according genes).

Conclusion

This study is the first to offer a comprehensive understanding of the gene regulation in ‘Wuzi Ougan’ and its wild type, especially during the microsporocyte to meiosis stage. Specifically, the involved genes include those in phenylpropanoid biosynthesis, flavonoid biosynthesis, and phenylalanine metabolism, as determined by integrative transcriptome and proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0693-9) contains supplementary material, which is available to authorized users.

Comparative analysis of the male inflorescence transcriptome profiles of an ms22 mutant of maize

In modern agricultural production, maize is the most successful crop utilizing heterosis. 712C-ms22 is an important male sterile material in maize. In this study, we performed transcriptome sequencing analysis of the V10 stage of male inflorescence. Through this analysis, 27.63 million raw reads were obtained, and trimming of the raw data revealed 26.63 million clean reads, with an average match rate of 94.64%. Using Tophat software, we matched these clean reads to the maize reference genome. The abundance of 39,622 genes was measured, and 35,399 genes remained after filtering out the non-expressed genes across all the samples. These genes were classified into 19 categories by clusters of orthologous groups of protein annotation. Transcriptome sequencing analysis of the male sterile and fertile 712C-ms22 maize revealed some key DEGs that may be related to metabolic pathways. qRT-PCR analysis validated the gene expression patterns identified by RNA-seq. This analysis revealed some of the essential genes responsible for pollen development and for pollen tube elongation. Our findings provide useful markers of male sterility and new insights into the global mechanisms mediating male sterility in maize.

iTRAQ-Based Proteomics Analyses of Sterile/Fertile Anthers from a Thermo-Sensitive Cytoplasmic Male-Sterile Wheat with Aegilops kotschyi Cytoplasm

A “two-line hybrid system” was developed, previously based on thermo-sensitive cytoplasmic male sterility in Aegilops kotschyi (K-TCMS), which can be used in wheat breeding. The K-TCMS line exhibits complete male sterility and it can be used to produce hybrid wheat seeds during the normal wheat-growing season; it propagates via self-pollination at high temperatures. Isobaric tags for relative and absolute quantification-based quantitative proteome and bioinformatics analyses of the TCMS line KTM3315A were conducted under different fertility conditions to understand the mechanisms of fertility conversion in the pollen development stages. In total, 4639 proteins were identified, the differentially abundant proteins that increased/decreased in plants with differences in fertility were mainly involved with energy metabolism, starch and sucrose metabolism, phenylpropanoid biosynthesis, protein synthesis, translation, folding, and degradation. Compared with the sterile condition, many of the proteins that related to energy and phenylpropanoid metabolism increased during the anther development stage. Thus, we suggest that energy and phenylpropanoid metabolism pathways are important for fertility conversion in K-TCMS wheat. These findings provide valuable insights into the proteins involved with anther and pollen development, thereby, helping to further understand the mechanism of TCMS in wheat.

De novo-based transcriptome profiling of male-sterile and fertile watermelon lines

The whole-genome sequence of watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai), a valuable horticultural crop worldwide, was released in 2013. Here, we compared a de novo-based approach (DBA) to a reference-based approach (RBA) using RNA-seq data, to aid in efforts to improve the annotation of the watermelon reference genome and to obtain biological insight into male-sterility in watermelon. We applied these techniques to available data from two watermelon lines: the male-sterile line DAH3615-MS and the male-fertile line DAH3615. Using DBA, we newly annotated 855 watermelon transcripts, and found gene functional clusters predicted to be related to stimulus responses, nucleic acid binding, transmembrane transport, homeostasis, and Golgi/vesicles. Among the DBA-annotated transcripts, 138 de novo-exclusive differentially-expressed genes (DEDEGs) related to male sterility were detected. Out of 33 randomly selected newly annotated transcripts and DEDEGs, 32 were validated by RT-qPCR. This study demonstrates the usefulness and reliability of the de novo transcriptome assembly in watermelon, and provides new insights for researchers exploring transcriptional blueprints with regard to the male sterility.

miRNAs involved in the development and differentiation of fertile and sterile flowers in Viburnum macrocephalum f. keteleeri

BACKGROUND

Sterile and fertile flowers are important evolutionary developmental phenotypes in angiosperm flowers. The development of floral organs, critical in angiosperm reproduction, is regulated by microRNAs (miRNAs). However, the mechanisms underpinning the miRNA regulation of the differentiation and development of sterile and fertile flowers remain unclear.

RESULTS

Here, based on investigations of the morphological differences between fertile and sterile flowers, we used high-throughput sequencing to characterize the miRNAs in the differentiated floral organs of Viburnum macrocephalum f. keteleeri. We identified 49 known miRNAs and 67 novel miRNAs by small RNA (sRNA) sequencing and bioinformatics analysis, and 17 of these known and novel miRNA precursors were validated by polymerase chain reaction (PCR) and Sanger sequencing. Furthermore, by comparing the sequencing results of two sRNA libraries, we found that 30 known and 39 novel miRNA sequences were differentially expressed, and 35 were upregulated and 34 downregulated in sterile compared with fertile flowers. Combined with their predicted targets, the potential roles of miRNAs in V. macrocephalum f. keteleeri flowers include involvement in floral organogenesis, cell proliferation, hormonal pathways, and stress responses. miRNA precursors and targets were further validated by quantitative real-time PCR (qRT-PCR). Specifically, miR156a-5p, miR156g, and miR156j expression levels were significantly higher in fertile flowers than in sterile flowers, while SPL genes displayed the opposite expression pattern. Considering that the targets of miR156 are predicted to be SPL genes, we propose that miR156 may be involved in the regulation of stamen development in V. macrocephalum f. keteleeri.

CONCLUSIONS

We identified miRNAs differentially expressed between fertile and sterile flowers in V. macrocephalum f. keteleeri and provided new insights into the important regulatory roles of miRNAs in the differentiation and development of fertile and sterile flowers.

Insights into the immune manipulation mechanisms of pollen allergens by protein domain profiling

Plant pollens are airborne allergens, as their inhalation causes immune activation, leading to rhinitis, conjunctivitis, sinusitis and oral allergy syndrome. A myriad of pollen proteins belonging to profilin, expansin, polygalacturonase, glucan endoglucosidase, pectin esterase, and lipid transfer protein class have been identified. In the present in silico study, the protein domains of fifteen pollen sequences were extracted from the UniProt database and submitted to the interactive web tool SMART (Simple Modular Architecture Research Tool), for finding the protein domain profiles. Analysis of the data based on custom-made scripts revealed the conservation of pathogenic domains such as OmpH, PROF, PreSET, Bet_v_1, Cpl-7 and GAS2. Further, the retention of critical domains like CHASE2, Galanin, Dak2, DALR_1, HAMP, PWI, EFh, Excalibur, CT, PbH1, HELICc, and Kelch in pollen proteins, much like cockroach allergens and lethal viruses (such as HIV, HCV, Ebola, Dengue and Zika) was observed. Based on the shared motifs in proteins of taxonomicall-ydispersed organisms, it can be hypothesized that allergens and pathogens manipulate the human immune system in a similar manner. Allergens, being inanimate, cannot replicate in human body, and are neutralized by immune system. But, when the allergens are unremitting, the immune system becomes persistently hyper-sensitized, creating an inflammatory milieu. This study is expected to contribute to the understanding of pollen allergenicity and pathogenicity.

Genome-wide identification of pistil-specific genes expressed during fruit set initiation in tomato (Solanum lycopersicum)

Fruit set involves the developmental transition of an unfertilized quiescent ovary in the pistil into a fruit. While fruit set is known to involve the activation of signals (including various plant hormones) in the ovary, many biological aspects of this process remain elusive. To further expand our understanding of this process, we identified genes that are specifically expressed in tomato (Solanum lycopersicum L.) pistils during fruit set through comprehensive RNA-seq-based transcriptome analysis using 17 different tissues including pistils at six different developmental stages. First, we identified 532 candidate genes that are preferentially expressed in the pistil based on their tissue-specific expression profiles. Next, we compared our RNA-seq data with publically available transcriptome data, further refining the candidate genes that are specifically expressed within the pistil. As a result, 108 pistil-specific genes were identified, including several transcription factor genes that function in reproductive development. We also identified genes encoding hormone-like peptides with a secretion signal and cysteine-rich residues that are conserved among some Solanaceae species, suggesting that peptide hormones may function as signaling molecules during fruit set initiation. This study provides important information about pistil-specific genes, which may play specific roles in regulating pistil development in relation to fruit set.

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

Background

Kiwifruit is a perennial, deciduous and functionally dioecious plant. However, very little is known about the whole-genome molecular mechanisms contributing to distinct sexual phenotypes. To gain a global view of genes differentially expressed between male and female flowers, we analyzed genome-wide gene expression profiles in the flowers of male and female plants using high-throughput RNA sequencing.

Results

A total of 53.5 million reads were generated. Based on the alignments of unigenes to kiwifruit genome predicted genes, a total of 39,040 unique genes with a mean length of 970 bp were identified. There were 2,503 UniGenes differentially expressed between female and male flowers, with 1,793 up-regulated and 710 down-regulated in the female flowers. Moreover, the gene expression pattern of 17 out of 19 unigenes differentially expressed between male and female flowers revealed by RNA-Seq was confirmed by real-time quantitative PCR (qRT-PCR).

Conclusions

Here, we obtained a large number of EST sequences from female and male flowers of kiwifruit. This comparative transcriptome analysis provides an invaluable resource for gene expression, genomics, and functional genomic studies in A. chinensis and its related species. This study also represents a first step toward the investigation of genes involved in kiwifruit sex determination.

IRREGULAR POLLEN EXINE1 Is a Novel Factor in Anther Cuticle and Pollen Exine Formation

Anther cuticle and pollen exine are protective barriers for pollen development and fertilization. Despite that several regulators have been identified for anther cuticle and pollen exine development in rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), few genes have been characterized in maize (Zea mays) and the underlying regulatory mechanism remains elusive. Here, we report a novel male-sterile mutant in maize, irregular pollen exine1 (ipe1), which exhibited a glossy outer anther surface, abnormal Ubisch bodies, and defective pollen exine. Using map-based cloning, the IPE1 gene was isolated as a putative glucose-methanol-choline oxidoreductase targeted to the endoplasmic reticulum. Transcripts of IPE1 were preferentially accumulated in the tapetum during the tetrad and early uninucleate microspore stage. A biochemical assay indicated that ipe1 anthers had altered constituents of wax and a significant reduction of cutin monomers and fatty acids. RNA sequencing data revealed that genes implicated in wax and flavonoid metabolism, fatty acid synthesis, and elongation were differentially expressed in ipe1 mutant anthers. In addition, the analysis of transfer DNA insertional lines of the orthologous gene in Arabidopsis suggested that IPE1 and their orthologs have a partially conserved function in male organ development. Our results showed that IPE1 participates in the putative oxidative pathway of C16/C18 ω-hydroxy fatty acids and controls anther cuticle and pollen exine development together with MALE STERILITY26 and MALE STERILITY45 in maize.

Identification of Candidate Genes and Biosynthesis Pathways Related to Fertility Conversion by Wheat KTM3315A Transcriptome Profiling

The Aegilops kotschyi thermo-sensitive cytoplasmic male sterility (K-TCMS) system may facilitate hybrid wheat (Triticum aestivum L.) seed multiplication and production. The K-TCMS line is completely male sterile during the normal wheat-growing season, whereas its fertility can be restored in a high-temperature environment. To elucidate the molecular mechanisms responsible for male sterility/fertility conversion and candidate genes involved with pollen development in K-TCMS, we employed RNA-seq to sequence the transcriptomes of anthers from K-TCMS line KTM3315A during development under sterile and fertile conditions. We identified 16840 differentially expressed genes (DEGs) in different stages including15157 known genes (15135 nuclear genes and 22 plasmagenes) and 1683 novel genes. Bioinformatics analysis identified possible metabolic pathways involved with fertility based on KEGG pathway enrichment of the DEGs expressed in fertile and sterile plants. We found that most of the genes encoding key enzyme in the phenylpropanoid biosynthesis and jasmonate biosynthesis pathways were significant upregulated in uninucleate, binuclate or trinucleate stage, which both interact with MYB transcription factors, and that link between all play essential roles in fertility conversion. The relevant DEGs were verified by quantitative RT-PCR. Thus, we suggested that phenylpropanoid biosynthesis and jasmonate biosynthesis pathways were involved in fertility conversion of K-TCMS wheat. This will provide a new perspective and an effective foundation for the research of molecular mechanisms of fertility conversion of CMS wheat. Fertility conversion mechanism in thermo-sensitive cytoplasmic male sterile/fertile wheat involves the phenylpropanoid biosynthesis pathway, jasmonate biosynthesis pathway, and MYB transcription factors.

Transcriptomic Analysis Reveals Mechanisms of Sterile and Fertile Flower Differentiation and Development in Viburnum macrocephalum f. keteleeri

Sterile and fertile flowers are an important evolutionary developmental (evo-devo) phenotype in angiosperm flowers, playing important roles in pollinator attraction and sexual reproductive success. However, the gene regulatory mechanisms underlying fertile and sterile flower differentiation and development remain largely unknown. Viburnum macrocephalum f. keteleeri, which possesses fertile and sterile flowers in a single inflorescence, is a useful candidate species for investigating the regulatory networks in differentiation and development. We developed a de novo-assembled flower reference transcriptome. Using RNA sequencing (RNA-seq), we compared the expression patterns of fertile and sterile flowers isolated from the same inflorescence over its rapid developmental stages. The flower reference transcriptome consisted of 105,683 non-redundant transcripts, of which 5,675 transcripts showed significant differential expression between fertile and sterile flowers. Combined with morphological and cytological changes between fertile and sterile flowers, we identified expression changes of many genes potentially involved in reproductive processes, phytohormone signaling, and cell proliferation and expansion using RNA-seq and qRT-PCR. In particular, many transcription factors (TFs), including MADS-box family members and ABCDE-class genes, were identified, and expression changes in TFs involved in multiple functions were analyzed and highlighted to determine their roles in regulating fertile and sterile flower differentiation and development. Our large-scale transcriptional analysis of fertile and sterile flowers revealed the dynamics of transcriptional networks and potentially key components in regulating differentiation and development of fertile and sterile flowers in Viburnum macrocephalum f. keteleeri. Our data provide a useful resource for Viburnum transcriptional research and offer insights into gene regulation of differentiation of diverse evo-devo processes in flowers.

Comparative transcriptome profiling of the fertile and sterile flower buds of a dominant genic male sterile line in sesame (Sesamum indicum L.)

BACKGROUND

Sesame (Sesamum indicum L.) is a globally important oilseed crop with highly-valued oil. Strong hybrid vigor is frequently observed within this crop, which can be exploited by the means of genic male sterility (GMS). We have previously developed a dominant GMS (DGMS) line W1098A that has great potential for the breeding of F₁ hybrids. Although it has been genetically and anatomically characterized, the underlying molecular mechanism for male sterility remains unclear and therefore limits the full utilization of such GMS line. In this study, RNA-seq based transcriptome profiling was carried out in two near-isogenic DGMS lines (W1098A and its fertile counterpart, W1098B) to identify differentially expressed genes (DEGs) related to male sterility.

RESULTS

A total of 1,502 significant DEGs were detected, among which 751 were up-regulated and 751 were down-regulated in sterile flower buds. A number of DEGs were implicated in both ethylene and JA synthesis & signaling pathway; the expression of which were either up- or down-regulated in the sterile buds, respectively. Moreover, the majority of NAC and WRKY transcription factors implicated from the DEGs were up-regulated in sterile buds. By querying the Plant Male Reproduction Database, 49 sesame homologous genes were obtained; several of these encode transcription factors (bHLH089, MYB99, and AMS) that showed reduced expression in sterile buds, thus implying the possible role in specifying or determining tapetal fate and development. The predicted effect of allelic variants on the function of their corresponding DEGs highlighted several Insertions/Deletions (InDels), which might be responsible for the phenotype of sterility/fertility in DGMS lines.

CONCLUSION

The present comparative transcriptome study suggested that both hormone signaling pathway and transcription factors control the male sterility of DGMS in sesame. The results also revealed that several InDels located in DEGs prone to cause loss of function, which might contribute to male sterility. These findings provide valuable genomic resources for a deeper insight into the molecular mechanism underlying DGMS.

RNA-Seq Library Construction Methods for Transcriptome Analysis

Next-generation sequencing (NGS) technologies have revolutionized the study of genomics with an ever-expanding list of applications. RNA-Seq has emerged as a powerful method, applying transcriptome analysis to a wider range of organisms-most significantly, non-model organisms lacking prior genomic sequencing. Whereas an initial concern of NGS datasets was the potential limitation of short read lengths, short read sequences have been successfully employed in creation of de novo transcriptome assemblies that allow for subsequent mapping of reads for expression analysis. Prior genomic sequence knowledge is no longer a requirement for identification of functional transcriptional elements and for global gene expression characterization. Significant cost reductions in generating RNA-Seq data, and improvements in de novo assemblers, has allowed the analysis of transcriptomes in heretofore unsequenced plant species. These protocols describe standard methods for constructing RNA-Seq libraries to be sequenced on Illumina sequencing platforms for comprehensive transcriptome analysis. © 2016 by John Wiley & Sons, Inc.

Comprehensive identification of sexually dimorphic genes in diverse cattle tissues using RNA-seq

BACKGROUND

Molecular mechanisms associated with sexual dimorphism in cattle have not been well elucidated. Furthermore, as recent studies have implied that gene expression patterns are highly tissue specific, it is essential to investigate gene expression in a variety of tissues using RNA-seq. Here, we employed and compared two statistical methods, a simple two group test and Analysis of deviance (ANODEV), in order to investigate bovine sexually dimorphic genes in 40 RNA-seq samples distributed across two factors: sex and tissue.

RESULTS

As a result, we detected 752 sexually dimorphic genes across tissues from two statistical approaches and identified strong tissue-specific patterns of gene expression. Additionally, significantly detected sex-related genes shared between two mammal species (cattle and rat) were identified using qRT-PCR.

CONCLUSIONS

Results of our analyses reveal that sexual dimorphism of metabolic tissues and pituitary gland in cattle involves various biological processes. Several differentially expressed genes between sexes in cattle and rat species are shared, but show tissue-specific patterns. Finally, we concluded that two distinct statistical approaches have their advantages and disadvantages in RNA-seq studies investigating multiple tissues.