Comparative Transcriptome Analysis of Recessive Male Sterility (RGMS) in Sterile and Fertile Brassica napus Lines

Identification of Trehalose-6-Phosphate Synthase (TPS) Genes Associated with Both Source-/Sink-Related Yield Traits and Drought Response in Rapeseed (Brassica napus L.)

Trehalose-6-phosphate synthase (TPS) is an important enzyme for the synthesis of Trehalose-6-phosphate (T6P). In addition to being a signaling regulator of carbon allocation that improves crop yields, T6P also plays essential roles in desiccation tolerance. However, comprehensive studies, such as evolutionary analysis, expression analysis, and functional classification of the TPS family in rapeseed (Brassica napus L.) are lacking. Here, we identified 35 BnTPSs, 14 BoTPSs, and 17 BrTPSs in cruciferous plants, which were classified into three subfamilies. Phylogenetic and syntenic analysis of TPS genes in four cruciferous species indicated that only gene elimination occurred during their evolution. Combined phylogenetic, protein property, and expression analysis of the 35 BnTPSs suggested that changes in gene structures might have led to changes in their expression profiles and further functional differentiation during their evolution. In addition, we analyzed one set of transcriptome data from Zhongshuang11 (ZS11) and two sets of data from extreme materials associated with source-/sink-related yield traits and the drought response. The expression levels of four BnTPSs (BnTPS6, BnTPS8, BnTPS9, and BnTPS11) increased sharply after drought stress, and three differentially expressed genes (BnTPS1, BnTPS5, and BnTPS9) exhibited variable expression patterns among source and sink tissues between yield-related materials. Our findings provide a reference for fundamental studies of TPSs in rapeseed and a framework for future functional research of the roles of BnTPSs in both yield and drought resistance.

Comparative transcriptome analysis provides insight into the important pathways and key genes related to the pollen abortion in the thermo-sensitive genic male sterile line 373S in Brassica napus L

The thermo-sensitive genic male sterility (TGMS) system plays a key role in the production of two-line hybrids in rapeseed (Brassica napus). To uncover key cellular events and genetic regulation associated with TGMS, a combined study using cytological methods and RNA-sequencing analysis was conducted for the rapeseed TGMS line 373S. Cytological studies showed that microspore cytoplasm of 373S plants was condensed, the microspore nucleus was degraded at an early stage, the exine was irregular, and the tapetum developed abnormally, eventually leading to male sterility. RNA-sequencing analysis identified 430 differentially expressed genes (298 upregulated and 132 downregulated) between the fertile and sterile samples. Gene ontology analysis demonstrated that the most highly represented biological processes included sporopollenin biosynthetic process, pollen exine formation, and extracellular matrix assembly. Kyoto encyclopedia of genes and genomes analysis indicated that the enriched pathways included amino acid metabolism, carbohydrate metabolism, and lipid metabolism. Moreover, 26 transcript factors were identified, which may be associated with abnormal tapetum degeneration and exine formation. Subsequently, 19 key genes were selected, which are considered to regulate pollen development and even participate in pollen exine formation. Our results will provide important insight into the molecular mechanisms underlying TGMS in rapeseed.

Conjunctive Analyses of BSA-Seq and BSR-Seq Unveil the Msβ-GAL and MsJMT as Key Candidate Genes for Cytoplasmic Male Sterility in Alfalfa (Medicago sativa L.)

Knowing the molecular mechanism of male sterility in alfalfa is important to utilize the heterosis more effectively. However, the molecular mechanisms of male sterility in alfalfa are still unclear. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of male sterility in alfalfa. The BSA-seq analysis was located in a candidate region on chromosome 5 containing 626 candidate genes which were associated with male sterility in alfalfa, while the BSR-seq analysis filtered seven candidate DEGs related to male sterility, and these candidate genes including EF-Tu, β-GAL, CESA, PHGDH, and JMT. The conjunctive analyses of BSR and BSA methods revealed that the genes of Msβ-GAL and MsJMT are the common detected candidate genes involved in male sterility in alfalfa. Our research provides a theory basis for further study of the molecular mechanism of male sterility in alfalfa and significant information for the genetic breeding of Medicago sativa.

Transcriptome and Small RNA Sequencing Reveal the Mechanisms Regulating Harvest Index in Brassica napus

Harvest index (HI), the ratio of harvested seed weight to total aboveground biomass weight, is an economically critical value reflecting the convergence of complex agronomic traits. HI values in rapeseed (Brassica napus) remain much lower than in other major crops, and the underlying regulatory network is largely unknown. In this study, we performed mRNA and small RNA sequencing to reveal the mechanisms shaping HI in B. napus during the seed-filling stage. A total of 8,410 differentially expressed genes (DEGs) between high-HI and low-HI accessions in four tissues (silique pericarp, seed, leaves, and stem) were identified. Combining with co-expression network, 72 gene modules were identified, and a key gene BnaSTY46 was found to participate in retarded establishment of photosynthetic capacity to influence HI. Further research found that the genes involved in circadian rhythms and response to stimulus may play important roles in HI and that their transcript levels were modulated by differentially expressed microRNAs (DEMs), and we identified 903 microRNAs (miRNAs), including 46 known miRNAs and 857 novel miRNAs. Furthermore, transporter activity-related genes were critical to enhancing HI in good cultivation environments. Of 903 miRNAs, we found that the bna-miR396-Bna.A06SRp34a/Bna.A01EMB3119 pair may control the seed development and the accumulation of storage compounds, thus contributing to higher HI. Our findings uncovered the underlying complex regulatory network behind HI and offer potential approaches to rapeseed improvement.

Comparative transcriptome analysis reveals differential gene expression in sterile and fertile rubber tree varieties during flower bud differentiation

Plant male sterility (MS) is an important agronomic trait that provides an efficient tool for hybridization and heterosis utilization of crops. Based on phenotypic and cytological observations, our study performed a multi-comparison transcriptome analysis strategy on multiple sterile and fertile rubber tree varieties using RNA-seq. Compared with the male-fertile varieties, a total of 1590 differentially expressed genes (DEGs) were detected in male-sterile varieties, including 970 up-regulated and 620 down-regulated transcripts in sterile varieties. Key DEGs were further assessed focusing on anther development, microsporogenesis and plant hormone metabolism. Twenty DEGs were selected randomly to validate transcriptome data using quantitative real-time PCR (qRT-PCR). Eleven key genes were subjected to expression pattern analysis using qRT-PCR and fluorescence in situ hybridization. Among them, nine genes, i.e., A6, GAI1, ACA7, TKPR1, CYP704B1, XTH26, MS1, MS35 and MYB33, that regulate callose metabolism, pollen wall formation, tapetum and microspores development were identified as candidate male-sterile genes. These findings provide insights into the molecular mechanism of male sterility in rubber tree.

The characterization and candidate gene isolation for a novel male-sterile mutant ms40 in maize

A novel genic male-sterile mutant ms40 was obtained from EMS treated RP125. The key candidate gene ZmbHLH51 located on chromosome 4 was identified by map-based cloning. This study further enriched the male sterile gene resources for both production applications and theoretical studies of abortion mechanisms. Maize male-sterile mutant 40 (ms40) was obtained from the progeny of the ethyl methanesulfonate (EMS) treated inbred line RP125. Genetic analysis indicated that the sterility was controlled by a single recessive nuclear gene. Cytological observation of anthers revealed that the cuticles of ms40 anthers were abnormal, and no Ubisch bodies were observed on the inner surface of ms40 anthers through scanning electron microscopy(SEM). Moreover, its tapetum exhibited delayed degradation and then blocked the formation of normal microspores. Using map-based cloning strategy, the ms40 locus was found to locate in a 282-kb interval on chromosome 4, and five annotated genes were predicted within this region. PCR-based sequencing detected a single non-synonymous SNP (G > A) that changed glycine (G) to arginine (A) in the seventh exon of Zm00001d053895, while no sequence difference between ms40 and RP125 was found for the other four genes. Zm00001d053895 encodes the bHLH transcription factor ZmbHLH51 which is localized in the nucleus. Phylogenetic analysis showed that ZmbHLH51 had the highest homology with Sb04g001650, a tapetum degeneration retardation (TDR) bHLH transcription factor in Sorghum bicolor. Co-expression analysis revealed a total of 1192 genes co-expressed with ZmbHLH51 in maize, 647 of which were anther-specific genes. qRT-PCR results suggested the expression levels of some known genes related to anther development were affected in ms40. In summary, these findings revealed the abortion characteristics of ms40 anthers and lay a foundation for further studies on the mechanisms of male fertility.

Transcriptome Profiling Reveals Molecular Changes during Flower Development between Male Sterile and Fertile Chinese Cabbage (Brassica rapa ssp. pekinensis) Lines

Male sterility exists widely in flowering plants and is used as a fascinating tool by breeders for creating hybrid varieties. Herein, stamen samples from male sterile CCR20000 and male fertile CCR20001 lines during two developmental stages were employed to elucidate the molecular changes during flower development in fertile and sterile Chinese cabbage lines. RNA-seq revealed weak transcriptional activity in the sterile line, which may have led to the abnormal stamen development. The differentially expressed genes were enriched in plant hormone, carbon metabolism, and biosynthesis of amino acid pathways. Important genes with opposite patterns of regulation between the two lines have been associated with the male sterility trait. Members of the transcription factor families such as AP2, MYB, bHLH, and WRKY were highly active in the regulation of structural genes involved in pollen fertility. This study generated important genomic information to support the exploitation of the male sterility trait in Chinese cabbage breeding programs.

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.

Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress

BACKGROUND

Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea).

RESULTS

In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content.

CONCLUSIONS

Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

Integrative analysis of transcriptomic and proteomic changes related to male sterility in Tagetes erecta

Male sterile and male fertile two-type lines are important in heterosis utilization and breeding in Tagetes erecta, but the genes and pathways involved in male sterility are poorly understood. To explore these topics, transcriptome data (by RNA-seq) and proteome data (by iTRAQ) were gathered from flower buds of the male sterile line 'MS2-2' and male fertile line 'MF2-2' and integrated for a better understanding of the underlying molecular mechanisms of male sterility in T. erecta. The RNA-seq procedure generated 285,139,740 clean reads and 63359 unigenes and 6640 differentially expressed genes (DEGs) were identified, of which 4136 were downregulated and 2504 were upregulated in 'MS2-2'. DEGs related to flower development, pollen development, pollen wall assembly, endogenous hormones and transcription factors were identified. The iTRAQ analysis identified 3950 proteins in total; 789 were differentially expressed proteins (381 upregulated, 408 downregulated), which were mainly annotated to the Ribosome, Carbon metabolism and Biosynthesis of amino acids pathways. An association analysis revealed strong correlation (r Pearson = 0.6019) between the transcriptomic and proteomic data, and 256 and 34 proteins showed the same and opposite expression patterns with regard to their transcripts, respectively. Pathways such as photosynthesis, fatty acid biosynthesis and phenylpropanoid biosynthesis which influence tapetum and pollen development in male sterile plants, were significantly enriched at the transcript and protein levels. Most genes involved in these pathways were downregulated in 'MS2-2'. The low expression of these genes or functional loss of proteins could be associated with flower development, pollen development and related to changes in fertility in T. erecta. This study provided transcriptomic and proteomic information for T. erecta that could illuminate the mechanism of male sterility.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Genome-Wide Identification and Expression Profiling of Monosaccharide Transporter Genes Associated with High Harvest Index Values in Rapeseed (Brassica napus L.)

Sugars are important throughout a plant’s lifecycle. Monosaccharide transporters (MST) are essential sugar transporters that have been identified in many plants, but little is known about the evolution or functions of MST genes in rapeseed (Brassica napus). In this study, we identified 175 MST genes in B. napus, 87 in Brassica oleracea, and 83 in Brassica rapa. These genes were separated into the sugar transport protein (STP), polyol transporter (PLT), vacuolar glucose transporter (VGT), tonoplast monosaccharide transporter (TMT), inositol transporter (INT), plastidic glucose transporter (pGlcT), and ERD6-like subfamilies, respectively. Phylogenetic and syntenic analysis indicated that gene redundancy and gene elimination have commonly occurred in Brassica species during polyploidization. Changes in exon-intron structures during evolution likely resulted in the differences in coding regions, expression patterns, and functions seen among BnMST genes. In total, 31 differentially expressed genes (DEGs) were identified through RNA-seq among materials with high and low harvest index (HI) values, which were divided into two categories based on the qRT-PCR results, expressed more highly in source or sink organs. We finally identified four genes, including BnSTP5, BnSTP13, BnPLT5, and BnERD6-like14, which might be involved in monosaccharide uptake or unloading and further affect the HI of rapeseed. These findings provide fundamental information about MST genes in Brassica and reveal the importance of BnMST genes to high HI in B. napus.

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.

Comparative Analysis of the Metabolic Profiles of Yellow- versus Black-Seeded Rapeseed Using UPLC-HESI-MS/MS and Transcriptome Analysis

The high levels of secondary metabolites in rapeseed play important roles in determining the oil quality and feeding value. Here, we characterized the metabolic profiles in seeds of various yellow- and black-seeded rapeseed accessions. Two hundred and forty-eight features were characterized, including 31 phenolic acids, 54 flavonoids, 24 glucosinolates, 65 lipid compounds, and 74 other polar compounds. The most abundant phenolic acids and various flavonoids (epicatechin, isorhamnetin, kaempferol, quercetin, and their derivatives) were widely detected and showed significant differences in distribution between the yellow- and black-seeded rapeseed. Furthermore, the related genes (e.g., BnTT3, BnTT18, BnTT10, BnTT12, and BnBAN) involved in the proanthocyanidin pathway had lower expression levels in yellow-seeded rapeseed, strongly suggesting that the seed coat color could be mainly determined by the levels of epicatechin and their derivatives. These results improve our understanding of the primary constituents of rapeseed and lay the foundation for breeding novel varieties with a high nutritional value.

Investigation of the genes associated with a male sterility mutant (msm) in Chinese cabbage (Brassica campestris ssp. pekinensis) using RNA-Seq

In Chinese cabbage, hybrid seed production is performed using male sterility lines, an important approach to heterosis utilization. In this study, a stably inherited male sterile mutant msm was obtained from the 'FT'-doubled haploid line of Chinese cabbage using isolated microspore culture combined with ⁶⁰Co γ-ray mutagenesis. The genetic backgrounds of 'FT' and msm were highly consistent; however, compared with wild-type 'FT', msm exhibited completely degenerated stamens and no pollen phenotype. Other characters showed no significant differences. Cytological observations revealed that stamen abortion in msm begins during the tetrad period and that tapetum cells were abnormally expanded and highly vacuolated, leading to microspore abortion. Genetic analysis indicated that the msm mutant phenotype is controlled by a single recessive nuclear gene. Comparative transcriptome analysis of 'FT' and msm flower buds using RNA-Seq technology revealed 1653 differentially expressed genes, among which, a large number associated with male sterility were detected, including 64 pollen development- and pollen tube growth-related genes, 94 pollen wall development-related genes, 11 phytohormone-related genes, and 16 transcription factor-related genes. An overwhelming majority of these genes were down-regulated in msm compared with 'FT'. Furthermore, KEGG pathway analysis indicated that a variety of carbohydrate metabolic and lipid metabolic pathways were significantly enriched, which may be related to pollen abortion. The expression patterns of 24 male sterility-related genes were analyzed using qRT-PCR. In addition, 24,476 single-nucleotide polymorphisms and 413,073 insertion-deletion events were specifically detected in msm. These results will facilitate elucidation of the regulatory mechanisms underlying male sterility in Chinese cabbage.

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.

Transcriptomic analysis reveals the mechanism of thermosensitive genic male sterility (TGMS) of Brassica napus under the high temperature inducement

Background

The thermo-sensitive genic male sterility (TGMS) of Brassica napus facilitates reproductive researches and hybrid seed production. Considering the complexity and little information about the molecular mechanism involved in B. napus TGMS, comparative transcriptomic analyses were peroformed for the sterile (160S-MS) and fertile (160S-MF) flowers to identify potential crucial genes and pathways associated with TGMS.

Results

In total, RNA-seq analysis showed that 2202 genes (561 up-regulated and 1641 down-regulated) were significantly differentially expressed in the fertile flowers of 160S-MF at 25 °C when compared the sterile flower of 160S-MS at 15 °C. Detailed analysis revealed that expression changes in genes encoding heat shock proteins, antioxidant, skeleton protein, GTPase and calmodulin might be involved in TGMS of B. napus. Moreover, gene expression of some key members in plant hormone signaling pathways, such as auxin, gibberellins, jasmonic acid, abscisic acid, brassinosteroid signalings, were significantly surppressed in the flowers of 160S, suggesting that these genes might be involved in the regulation in B. napus TGMS. Here, we also found that transcription factor MADS, NFY, HSF, MYB/C and WRKY might play a crucial role in male fertility under the high temperature condition.

Conclusion

High temperature can significant affect gene expression in the flowers. The findings in the current study improve our understanding of B. napus TGMS at the molecular level and also provide an effective foundation for male fertility researches in other important economic crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6008-3) contains supplementary material, which is available to authorized users.

Analysis of metabolic pathways related to fertility restoration and identification of fertility candidate genes associated with Aegilops kotschyi cytoplasm in wheat (Triticum aestivum L.)

BACKGROUND

Thermo-sensitive male-sterility based on Aegilops kotschyi cytoplasm (K-TCMS) plays an important role in hybrid wheat breeding. This has important possible applications in two-line hybrid wheat breeding but the genetic basis and molecular regulation mechanism related to fertility restoration are poorly understood. In this study, comparative transcriptome profiling based on RNA sequencing was conducted for two near-isogenic lines comprising KTM3315R and its sterile counterpart KTM3315A, a total of six samples (3 repetitions per group), in order to identify fertility restoration genes and their metabolic pathways.

RESULTS

In total, 2642 significant differentially expressed genes (DEGs) were detected, among which 1238 were down-regulated and 1404 were up-regulated in fertile anthers. Functional annotation enrichment analysis identified important pathways related to fertility restoration, such as carbohydrate metabolism, phenylpropanoid metabolism and biosynthesis, as well as candidate genes encoding pectin methylesterase and flavanone 3-hydroxylase. Moreover, transcription factor analysis showed that a large number of DEGs were mainly involved with the WRKY, bHLH, and MYB transcription factor families. Determination of total soluble sugar and flavonoid contents demonstrated that important metabolic pathways and candidate genes are associated with fertility restoration. Twelve DEGs were selected and detected by quantitative reverse-transcribed PCR, and the results indicated that the transcriptome sequencing results were reliable.

CONCLUSIONS

Our results indicate that identified DEGs were related to the fertility restoration and they proved to be crucial in Aegilops kotschyi cytoplasm. These findings also provide a basis for exploring the molecular regulation mechanism associated with wheat fertility restoration as well as screening and cloning related genes.

Comparative transcriptome analysis and ChIP-sequencing reveals stage-specific gene expression and regulation profiles associated with pollen wall formation in Brassica rapa

Background

Genic male sterility (GMS) line is an important approach to utilize heterosis in Brassica rapa, one of the most widely cultivated vegetable crops in Northeast Asia. However, the molecular genetic mechanisms of GMS remain to be largely unknown.

Results

Detailed phenotypic observation of ‘Bcajh97-01A/B’, a B. rapa genic male sterile AB line in this study revealed that the aberrant meiotic cytokinesis and premature tapetal programmed cell death occurring in the sterile line ultimately resulted in microspore degeneration and pollen wall defect. Further gene expression profile of the sterile and fertile floral buds of ‘Bcajh97-01A/B’ at five typical developmental stages during pollen development supported the result of phenotypic observation and identified stage-specific genes associated with the main events associated with pollen wall development, including tapetum development or functioning, callose metabolism, pollen exine formation and cell wall modification. Additionally, by using ChIP-sequencing, the genomic and gene-level distribution of trimethylated histone H3 lysine 4 (H3K4) and H3K27 were mapped on the fertile floral buds, and a great deal of pollen development-associated genes that were covalently modified by H3K4me³ and H3K27me³ were identified.

Conclusions

Our study provids a deeper understanding into the gene expression and regulation network during pollen development and pollen wall formation in B. rapa, and enabled the identification of a set of candidate genes for further functional annotation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5637-x) contains supplementary material, which is available to authorized users.

BrSKS13, a multiple-allele-inherited male sterility-related gene in Chinese cabbage (Brassica rapa L. ssp. pekinensis), affects pollen development and pollination/fertilization process

Multiple-allele-inherited male sterility (MAMS) is important in Chinese cabbage (Brassica rapa L. ssp. pekinensis) breeding, but the molecular mechanisms leading to male sterility are poorly understood. In this study, we cloned a novel gene, BrSKS13, that is differentially expressed in fertile and sterile flower buds of Chinese cabbage. BrSKS13 is most similar to Arabidopsis thaliana AT3G13400 (SKS13) and encodes a predicted 61.87 kDa protein with three cupredoxin superfamily conserved domains in the multicopper oxidase family. Semi-quantitative reverse-transcription PCR (sqRT-PCR) showed that expression of BrSKS13 is higher in fertile buds than in sterile buds. Quantitative RT-PCR (qRT-PCR) and in situ hybridization showed that BrSKS13 is highly expressed in fertile anthers, peaking at pollen-maturation stage VI, but is weakly expressed in other tissues and floral organs. Expression patterns of BrSKS13 promoter::GUS reporter fusions in Arabidopsis showed that the BrSKS13 promoter drives expression of the GUS gene only in anthers. The relative expression of Brsks13 in fertile buds was higher than in sterile buds for all other MAMS lines of Chinese cabbage examined. These results suggest that BrSKS13 affects pollen development. In situ hybridization analysis of flower stigmas at different times after pollination showed that BrSKS13 expression was first observed in stigmas and immature seeds at 1 h after pollination, and the signal intensity in seeds increased with increasing maturity. Compared to Col-0, A. thaliana sks13 mutant plants have shorter and fewer siliques, shriveled pollen grains, pollen tube abnormalities, and reduced seed number. The phenotype of sks13 mutant was recovered by over-expressing BrSKS13. Our results suggest that BrSKS13 affects pollen development and the pollination/fertilization process, and will enable further study of the genetic mechanisms underlying MAMS in Chinese cabbage.

Fine mapping of the major QTL for seed coat color in Brassica rapa var. Yellow Sarson by use of NIL populations and transcriptome sequencing for identification of the candidate genes

Yellow seed is a desirable trait in Brassica oilseed crops. The B. rapa var. Yellow Sarson carry unique yellow seed color genes which are not only important for the development of yellow-seeded oilseed B. rapa cultivars but this variant can also be used to develop yellow-seeded B. napus. In this study, we developed near-isogenic lines (NILs) of Yellow Sarson for the major seed coat color QTL SCA9-2 of the chromosome A9 and used the NILs to fine map this QTL region and to identify the candidate genes through linkage mapping and transcriptome sequencing of the developing seeds. From the 18.4 to 22.79 Mb region of SCA9-2, six SSR markers showing 0.63 to 5.65% recombination were developed through linkage analysis and physical mapping. A total of 55 differentially expressed genes (DEGs) were identified in the SCA9-2 region through transcriptome analysis; these included three transcription factors, Bra028039 (NAC), Bra023223 (C2H2 type zinc finger), Bra032362 (TIFY), and several other genes which encode unknown or nucleic acid binding protein; these genes might be the candidates and involved in the regulation of seed coat color in the materials used in this study. Several biosynthetic pathways, including the flavonoid, phenylpropanoid and suberin biosynthetic pathways were significantly enriched through GO and KEGG enrichment analysis of the DEGs. This is the first comprehensive study to understand the yellow seed trait of Yellow Sarson through employing linkage mapping and global transcriptome analysis approaches.

Transcriptome and Hormone Comparison of Three Cytoplasmic Male Sterile Systems in Brassica napus

The interaction between plant mitochondria and the nucleus markedly influences stress responses and morphological features, including growth and development. An important example of this interaction is cytoplasmic male sterility (CMS), which results in plants producing non-functional pollen. In current research work, we compared the phenotypic differences in floral buds of different Brassica napus CMS (Polima, Ogura, Nsa) lines with their corresponding maintainer lines. By comparing anther developmental stages between CMS and maintainer lines, we identified that in the Nsa CMS line abnormality occurred at the tetrad stage of pollen development. Phytohormone assays demonstrated that IAA content decreased in sterile lines as compared to maintainer lines, while the total hormone content was increased two-fold in the S₂ stage compared with the S₁ stage. ABA content was higher in the S₁ stage and exhibited a two-fold decreasing trend in S₂ stage. Sterile lines however, had increased ABA content at both stages compared with the corresponding maintainer lines. Through transcriptome sequencing, we compared differentially expressed unigenes in sterile and maintainer lines at both (S₁ and S₂) developmental stages. We also explored the co-expressed genes of the three sterile lines in the two stages and classified these genes by gene function. By analyzing transcriptome data and validating by RT-PCR, it was shown that some transcription factors (TFs) and hormone-related genes were weakly or not expressed in the sterile lines. This research work provides preliminary identification of the pollen abortion stage in Nsa CMS line. Our focus on genes specifically expressed in sterile lines may be useful to understand the regulation of CMS.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.

De novo transcriptome reconstruction of a thermo-sensitive male sterility mutant in rapeseed (Brassica napus; Brassicaceae)

PREMISE OF THE STUDY

SP2S is a spontaneous thermo-sensitive genic male sterility (TGMS) mutation that facilitates two-line hybrid breeding in Brassica napus (Brassicaceae). De novo assembly of the floral bud transcriptome of SP2S can provide a foundation for deciphering the transcriptional regulation of SP2S in response to temperature change.

METHODS

mRNAs of the young floral buds of SP2S and its near-isogenic line SP2F grown under cool (16°C)/warm (22°C) conditions were sequenced on an Illumina Solexa platform, producing 239.7 million short reads with a total length of 19.95 Gbp.

RESULTS

The reads were assembled de novo using the Trinity program, resulting in 135,702 transcripts with an average length of 784 bp, an N50 value of 1221 bp, and a total length of 107 Mbp. We identified 24,157 cDNA-derived simple sequence repeats in the assembly. We found 137 and 195 single-nucleotide polymorphisms and 49 and 51 differentially regulated KEGG orthology groups when comparing sample SP2S at 22°C vs. SP2S at 16°C and sample SP2S at 22°C vs. SP2F at 22°C, respectively.

DISCUSSION

The numerous differentially expressed genes and the derived single-nucleotide polymorphisms show abnormal transcriptional regulation in the TGMS system. These results outline an intricate transcriptional regulation that occurred in the rapeseed TGMS SP2S when the temperature changed.

Cloning, molecular evolution and functional characterization of ZmbHLH16, the maize ortholog of OsTIP2 (OsbHLH142)

The transcription factor ZmbHLH16, the maize ortholog of OsTIP2 (OsbHLH142), was isolated in the present study. Tissue expression analysis showed that ZmbHLH16 is preferentially expressed in male reproductive organs. Subcellular location analysis of ZmbHLH16 via rice protoplast indicated that it is located in the nucleus. Through nucleotide variation analysis, 36 polymorphic sites in ZmbHLH16, including 23 single nucleotide polymorphisms and 13 InDels, were detected among 78 maize inbred lines. Neutrality tests and linkage disequilibrium analysis showed that ZmbHLH16 experienced no significant evolutionary pressure. Yeast one-hybrid experiment showed that the first 80 residues in the N-terminus of ZmbHLH16 had transactivation activity, whereas the full length did not. Genome-wide coexpression analysis showed that 395 genes were coexpressed with ZmbHLH16. Among these genes, the transcription factor ZmbHLH51 had similar expression pattern and identical subcellular localization to those of ZmbHLH16. Subsequently, the interaction between ZmbHLH51 and ZmbHLH16 was verified by yeast two-hybrid experiment. Through yeast two-hybrid analysis of series truncated ZmbHLH16 fragments, we found not only the typical bHLH domain [175-221 amino acids (a.a.)], but also that the 81-160 a.a. and 241-365 a.a. of ZmbHLH16 could interact with ZmbHLH51. All these results lay the foundation for further understanding the functions of ZmbHLH16.

Genome-Wide Identification and Structural Analysis of bZIP Transcription Factor Genes in Brassica napus

The basic region/leucine zipper motif (bZIP) transcription factor family is one of the largest families of transcriptional regulators in plants. bZIP genes have been systematically characterized in some plants, but not in rapeseed (Brassica napus). In this study, we identified 247 BnbZIP genes in the rapeseed genome, which we classified into 10 subfamilies based on phylogenetic analysis of their deduced protein sequences. The BnbZIP genes were grouped into functional clades with Arabidopsis genes with similar putative functions, indicating functional conservation. Genome mapping analysis revealed that the BnbZIPs are distributed unevenly across all 19 chromosomes, and that some of these genes arose through whole-genome duplication and dispersed duplication events. All expression profiles of 247 bZIP genes were extracted from RNA-sequencing data obtained from 17 different B. napus ZS11 tissues with 42 various developmental stages. These genes exhibited different expression patterns in various tissues, revealing that these genes are differentially regulated. Our results provide a valuable foundation for functional dissection of the different BnbZIP homologs in B. napus and its parental lines and for molecular breeding studies of bZIP genes in B. napus.

TGMS in Rapeseed (Brassica napus) Resulted in Aberrant Transcriptional Regulation, Asynchronous Microsporocyte Meiosis, Defective Tapetum, and Fused Sexine

The thermo-sensitive genic male sterility (TGMS) line SP2S is a spontaneous rapeseed mutation with several traits that are favorable for the production of two-line hybrids. To uncover the key cellular events and genetic regulation associated with TGMS expression, a combined study using cytological observation, transcriptome profiling, and gene expression analysis was conducted for SP2S and its near-isogenic line SP2F grown under warm conditions. Asynchronous microsporocyte meiosis and abnormal tapetal plastids and elaioplasts were demonstrated in the anther of SP2S. The tetrad microspore did not undergo mitosis before the cytoplasm degenerated. Delayed degradation of the tetrad wall, which led to tetrad microspore aggregation, resulted in postponement of sexine (outer layer of pollen exine) formation and sexine fusion in the tetrad. The nexine (foot layer of exine) was also absent. The delay of tetrad wall degradation and abnormality of the exine structure suggested that the defective tapetum lost important functions. Based on transcriptomic comparisons between young flower buds of SP2S and SP2F plants, a total of 465 differentially expressed transcripts (DETs) were identified, including 303 up-regulated DETs and 162 down-regulated DETs in SP2S. Several genes encoding small RNA degrading nuclease 2, small RNA 2'-O-methyltransferase, thioredoxin reductase 2, regulatory subunit A alpha isoform of serine/threonine-protein phosphatase 2A, glycine rich protein 1A, transcription factor bHLH25, leucine-rich repeat receptor kinase At3g14840 like, and fasciclin-like arabinogalactan proteins FLA19 and FLA20 were greatly depressed in SP2S. Interestingly, a POLLENLESS3-LIKE 2 gene encoding the Arabidopsis MS5 homologous protein, which is necessary for microsporocyte meiosis, was down-regulated in SP2S. Other genes that were up-regulated in SP2S encoded glucanase A6, ethylene-responsive transcription factor 1A-like, pollen-specific SF3, stress-associated endoplasmic reticulum protein 2, WRKY transcription factors and pentatricopeptide repeat (PPR) protein At1g07590. The tapetum-development-related genes, including BnEMS1, BnDYT1, and BnAMS, were slightly up-regulated in 3-mm-long flower buds or their anthers, and their downstream genes, BnMS1 and BnMYB80, which affect callose dissolution and exine formation, were greatly up-regulated in SP2S. This aberrant genetic regulation corresponded well with the cytological abnormalities. The results suggested that expression of TGMS associates with complex transcriptional regulation.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

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.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.

Transcriptomic Analysis of Differentially Expressed Genes during Flower Organ Development in Genetic Male Sterile and Male Fertile Tagetes erecta by Digital Gene-Expression Profiling

Tagetes erecta is an important commercial plant of Asteraceae family. The male sterile (MS) and male fertile (MF) two-type lines of T. erecta have been utilized in F1 hybrid production for many years, but no report has been made to identify the genes that specify its male sterility that is caused by homeotic conversion of floral organs. In this study, transcriptome assembly and digital gene expression profiling were performed to generate expression profiles of MS and MF plants. A cDNA library was generated from an equal mixture of RNA isolated from MS and MF flower buds (1 mm and 4 mm in diameter). Totally, 87,473,431 clean tags were obtained and assembled into 128,937 transcripts among which 65,857 unigenes were identified with an average length of 1,188 bp. About 52% of unigenes (34,176) were annotated in Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KO (KEGG Ortholog database) and/or GO. Taking the above transcriptome as reference, 125 differentially expressed genes were detected in both developmental stages of MS and MF flower buds. MADS-box genes were presumed to be highly related to male sterility in T. erecta based on histological and cytological observations. Twelve MADS-box genes showed significantly different expression levels in flower buds 4 mm in diameter, whereas only one gene expressed significantly different in flower buds 1 mm in diameter between MS and MF plants. This is the first transcriptome analysis in T. erecta and will provide a valuable resource for future genomic studies, especially in flower organ development and/or differentiation.