Anatomy and transcript profiling of gynoecium development in female sterile Brassica napus mediated by one alien chromosome from Orychophragmus violaceus

The Transcriptome and Metabolome Reveal the Potential Mechanism of Lodging Resistance in Intergeneric Hybrids between Brassica napus and Capsella bursa-pastoris

Lodging is one of the main reasons for the reduction in seed yield and is the limitation of mechanized harvesting in B. napus. The dissection of the regulatory mechanism of lodging resistance is an important goal in B. napus. In this study, the lodging resistant B. napus line, YG689, derived from the hybridization between B. napus cv. Zhongyou 821 (ZY821) and Capsella bursa-pastoris, was used to dissect the regulation mechanism of hard stem formation by integrating anatomical structure, transcriptome and metabolome analyses. It was shown that the lignocellulose content of YG689 is higher than that of ZY821, and some differentially expressed genes (DEGs) involved in the lignocellulose synthesis pathway were revealed by transcriptome analyses. Meanwhile, GC–TOF–MS and UPLC–QTOF–MS identified 40, 54, and 31 differential metabolites in the bolting stage, first flower stage, and the final flower stage. The differential accumulation of these metabolites might be associated with the lignocellulose biosynthesis in B. napus. Finally, some important genes that regulate the metabolic pathway of lignocellulose biosynthesis, such as BnaA02g18920D, BnaA10g15590D, BnaC05g48040D, and NewGene_216 were identified in B. napus through the combination of transcriptomics and metabolomics data. The present results explored the potential regulatory mechanism of lignocellulose biosynthesis, which provided a new clue for the breeding of B. napus with lodging resistance in the future.

Screening of Differentially Expressed Genes and Localization Analysis of Female Gametophyte at the Free Nuclear Mitosis Stage in Pinus tabuliformis Carr

Female sterility is a common phenomenon in the plant world, and systematic research has not been carried out in gymnosperms. In this study, the ovules of No. 28 sterile line and No. 15 fertile line Pinus tabuliformis were used as materials, and a total of 18 cDNA libraries were sequenced by the HiSeqTM 4000 platform to analyze the differentially expressed genes (DEGs) and simple sequence repeats (SSRs) between the two lines. In addition, this study further analyzed the DEGs involved in the signal transduction of plant hormones, revealing that the signal pathways related to auxin, cytokinin, and gibberellin were blocked in the sterile ovule. Additionally, real-time fluorescent quantitative PCR verified that the expression trend of DEGs related to plant hormones was consistent with the results of high-throughput sequencing. Frozen sections and fluorescence in situ hybridization (FISH) were used to study the temporal and spatial expression patterns of PtRab in the ovules of P. tabuliformis. It was found that PtRab was significantly expressed in female gametophytes and rarely expressed in the surrounding diploid tissues. This study further explained the molecular regulation mechanism of female sterility in P. tabuliformis, preliminarily mining the key factors of ovule abortion in gymnosperms at the transcriptional level.

Dosage-Dependent Gynoecium Development and Gene Expression in Brassica napus-Orychophragmus violaceus Addition Lines

Distant hybridization usually leads to female sterility of the hybrid but the mechanism behind this is poorly understood. Complete pistil abortion but normal male fertility was shown by one Brassica napus-Orychophragmus violaceus monosomic alien addition line (MA, AACC + 1 IO, 2n = 39) produced previously. To study the effect of a single O. violaceus chromosome addition on pistil development in different genetic backgrounds, hybrids between the MA and B. carinata (BBCC), B. juncea (AABB), and two synthetic hexaploids (AABBCC) were firstly produced in this study which show complete female sterility. A microspore culture was further performed to produce the haploid monosomic alien addition line (HMA, AC + 1 IO, 2n = 20) and disomic addition line (DA, AACC + 2 IO, 2n = 40) together with haploid (H, AC, 2n = 19) and double haploid (DH, AACC, 2n = 38) plants of B. napus from MA to investigate the dosage effect of the alien O. violaceus chromosome on pistil development and gene expression. Compared to MA, the development of the pistils of DA and HMA was completely or partially recovered, in which the pistils could swell and elongate to a normal shape after open pollination, although no seeds were produced. Comparative RNA-seq analyses revealed that the numbers of the differentially expressed genes (DEGs) were significantly different, dosage-dependent, and consistent with the phenotypic difference in pairwise comparisons of HMA vs. H, DA vs. DH, MA vs. DH, MA vs. DA, and MA vs. HMA. The gene ontology (GO) enrichment analysis of DEGs showed that a number of genes involved in the development of the gynoecium, embryo sac, ovule, and integuments. Particularly, several common DEGs for pistil development shared in HMA vs. H and DA vs. DH showed functions in genotoxic stress response, auxin transport, and signaling and adaxial/abaxial axis specification. The results provided updated information for the molecular mechanisms behind the gynoecium development of B. napus responding to the dosage of alien O. violaceus chromosomes.

Combined Transcriptome Analysis Reveals the Ovule Abortion Regulatory Mechanisms in the Female Sterile Line of Pinus tabuliformis Carr

Ovule abortion is a common phenomenon in plants that has an impact on seed production. Previous studies of ovule and female gametophyte (FG) development have mainly focused on angiosperms, especially in Arabidopsis thaliana. However, because it is difficult to acquire information about ovule development in gymnosperms, this remains unclear. Here, we investigated the transcriptomic data of natural ovule abortion mutants (female sterile line, STE) and the wild type (female fertile line, FER) of Pinus tabuliformis Carr. to evaluate the mechanism of ovule abortion during the process of free nuclear mitosis (FNM). Using single-molecule real-time (SMRT) sequencing and next-generation sequencing (NGS), 18 cDNA libraries via Illumina and two normalized libraries via PacBio, with a total of almost 400,000 reads, were obtained. Our analysis showed that the numbers of isoforms and alternative splicing (AS) patterns were significantly variable between FER and STE. The functional annotation results demonstrate that genes involved in the auxin response, energy metabolism, signal transduction, cell division, and stress response were differentially expressed in different lines. In particular, AUX/IAA, ARF2, SUS, and CYCB had significantly lower expression in STE, showing that auxin might be insufficient in STE, thus hindering nuclear division and influencing metabolism. Apoptosis in STE might also have affected the expression levels of these genes. To confirm the transcriptomic analysis results, nine pairs were confirmed by quantitative real-time PCR. Taken together, these results provide new insights into ovule abortion in gymnosperms and further reveal the regulatory mechanisms of ovule development.

The SAP function in pistil development was proved by two allelic mutations in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Pistil development is a complicated process in plants, and female sterile mutants are ideal material for screening and cloning pistil development-related genes. Using the female sterile mutant (fsm1), BraA04g009730.3C was previously predicted as a candidate mutant gene encoding the STERILE APETALA (SAP) transcriptional regulator. In the current study, a parallel female sterile mutant (fsm2) was derived from EMS mutagenesis of a Chinese cabbage DH line 'FT' seeds.

RESULTS

Both fsm2 and fsm1 mutant phenotypes exhibited pistil abortion and smaller floral organs. Genetic analysis indicated that the phenotype of mutant fsm2 was also controlled by a single recessive nuclear gene. Allelism testing showed that the mutated fsm1 and fsm2 genes were allelic. A single-nucleotide mutation (G-to-A) in the first exon of BraA04g009730.3C caused a missense mutation from GAA (glutamic acid) to GGA (glycine) in mutant fsm2 plants. Both allelic mutations of BraA04g009730.3C in fsm1 and fsm2 conferred the similar pistil abortion phenotype, which verified the SAP function in pistil development. To probe the mechanism of SAP-induced pistil abortion, we compared the mutant fsm1 and wild-type 'FT' pistil transcriptomes. Among the 3855 differentially expressed genes obtained, 29 were related to ovule development and 16 were related to organ size.

CONCLUSION

Our study clarified the function of BraA04g009730.3C and revealed that it was responsible for ovule development and organ size. These results lay a foundation to elucidate the molecular mechanism of pistil development in Chinese cabbage.

Anatomy and RNA-Seq reveal important gene pathways regulating sex differentiation in a functionally Androdioecious tree, Tapiscia sinensis

BACKGROUND

Gametogenesis is a key step in the production of ovules or pollen in higher plants. The sex-determination aspects of gametogenesis have been well characterized in the model plant Arabidopsis. However, little is known about this process in androdioecious plants. Tapiscia sinensis Oliv. is a functionally androdioecious tree, with both male and hermaphroditic individuals. Hermaphroditic flowers (HFs) are female-fertile flowers that can produce functional pollen and set fruits. However, compared with male flowers (MFs), the pollen viability and number of pollen grains per flower are markedly reduced in HFs. MFs are female-sterile flowers that fail to set fruit and that eventually drop.

RESULTS

Compared with HF, a notable cause of MF female sterility in T. sinensis is when the early gynoecium meristem is disrupted. During the early stage of HF development (stage 6), the ring meristem begins to form as a ridge around the center of the flower. At this stage, the internal fourth-whorl organ is stem-like rather than carpelloid in MF. A total of 52,945 unigenes were identified as transcribed in MF and HF. A number of differentially expressed genes (DEGs) and metabolic pathways were detected as involved in the development of the gynoecium, especially the ovule, carpel and style. At the early gynoecium development stage, DEGs were shown to function in the metabolic pathways regulating ethylene biosynthesis and signal transduction (upstream regulator), auxin, cytokinin transport and signalling, and sex determination (or flower meristem identity).

CONCLUSIONS

Pathways for the female sterility model were initially proposed to shed light on the molecular mechanisms of gynoecium development at early stages in T. sinensis.

Identification of genomic regions associated with multi-silique trait in Brassica napus

BACKGROUND

Although rapeseed (Brassica napus L.) mutant forming multiple siliques was morphologically described and considered to increase the silique number per plant, an important agronomic trait in this crop, the molecular mechanism underlying this beneficial trait remains unclear. Here, we combined bulked-segregant analysis (BSA) and whole genome re-sequencing (WGR) to map the genomic regions responsible for the multi-silique trait using two pools of DNA from the near-isogenic lines (NILs) zws-ms (multi-silique) and zws-217 (single-silique). We used the Euclidean Distance (ED) to identify genomic regions associated with this trait based on both SNPs and InDels. We also conducted transcriptome sequencing to identify differentially expressed genes (DEGs) between zws-ms and zws-217.

RESULTS

Genetic analysis using the ED algorithm identified three SNP- and two InDel-associated regions for the multi-silique trait. Two highly overlapped parts of the SNP- and InDel-associated regions were identified as important intersecting regions, which are located on chromosomes A09 and C08, respectively, including 2044 genes in 10.20-MB length totally. Transcriptome sequencing revealed 129 DEGs between zws-ms and zws-217 in buds, including 39 DEGs located in the two abovementioned associated regions. We identified candidate genes involved in multi-silique formation in rapeseed based on the results of functional annotation.

CONCLUSIONS

This study identified the genomic regions and candidate genes related to the multi-silique trait in rapeseed.

Conserved and novel responses to cytokinin treatments during flower and fruit development in Brassica napus and Arabidopsis thaliana

Hormones are an important component in the regulatory networks guiding plant development. Cytokinins are involved in different physiological and developmental processes in plants. In the model plant Arabidopsis thaliana, cytokinin application during gynoecium development produces conspicuous phenotypes. On the other hand, Brassica napus, also known as canola, is a crop plant belonging to the Brassicaceae family, as A. thaliana. This makes B. napus a good candidate to study whether the cytokinin responses observed in A. thaliana are conserved in the same plant family. Here, we observed that cytokinin treatment in B. napus affects different traits of flower and fruit development. It increases ovule and seed number, affects stamen filament elongation and anther maturation, and causes a conspicuous overgrowth of tissue in petals and gynoecia. Furthermore, cytokinin recovers replum development in both wild type B. napus and in the A. thaliana rpl ntt double mutant, in which no replum is visible. These results indicate both conserved and novel responses to cytokinin in B. napus. Moreover, in this species, some cytokinin-induced phenotypes are inherited to the next, untreated generation, suggesting that cytokinins may trigger epigenetic modifications.

Production of red-flowered oilseed rape via the ectopic expression of Orychophragmus violaceus OvPAP2

Oilseed rape (Brassica napus L.), which has yellow flowers, is both an important oil crop and a traditional tourism resource in China, whereas the Orychophragmus violaceus, which has purple flowers, likely possesses a candidate gene or genes to alter the flower colour of oilseed rape. A previously established B. napus line has a particular pair of O. violaceus chromosomes (M4) and exhibits slightly red petals. In this study, the transcriptomic analysis of M4, B. napus (H3), and O. violaceus with purple petals (OvP) and with white petals (OvW) revealed that most anthocyanin biosynthesis genes were up-regulated in both M4 and OvP. Read assembly and sequence alignment identified a homolog of AtPAP2 in M4, which produced the O. violaceus transcript (OvPAP2). The overexpression of OvPAP2 via the CaMV35S promoter in Arabidopsis thaliana led to different levels of anthocyanin accumulation in most organs, including the petals. However, the B. napus overexpression plants showed anthocyanin accumulation primarily in the anthers, but not the petals. However, when OvPAP2 was driven by the petal-specific promoter XY355, the transgenic B. napus plants produced red anthers and red petals. The results of metabolomic experiments showed that specific anthocyanins accumulated to high levels in the red petals. This study illustrates the feasibility of producing red-flowered oilseed rape, thereby enhancing its ornamental value, via the ectopic expression of the OvPAP2 gene. Moreover, the practical application of this study for insect pest management in the crop is discussed.

RNA-sequencing-based transcriptome and biochemical analyses of steroidal saponin pathway in a complete set of Allium fistulosum-A. cepa monosomic addition lines

The genus Allium is a rich source of steroidal saponins, and its medicinal properties have been attributed to these bioactive compounds. The saponin compounds with diverse structures play a pivotal role in Allium’s defense mechanism. Despite numerous studies on the occurrence and chemical structure of steroidal saponins, their biosynthetic pathway in Allium species is poorly understood. The monosomic addition lines (MALs) of the Japanese bunching onion (A. fistulosum, FF) with an extra chromosome from the shallot (A. cepa Aggregatum group, AA) are powerful genetic resources that enable us to understand many physiological traits of Allium. In the present study, we were able to isolate and identify Alliospiroside A saponin compound in A. fistulosum with extra chromosome 2A from shallot (FF2A) and its role in the defense mechanism against Fusarium pathogens. Furthermore, to gain molecular insight into the Allium saponin biosynthesis pathway, high-throughput RNA-Seq of the root, bulb, and leaf of AA, MALs, and FF was carried out using Illumina's HiSeq 2500 platform. An open access Allium Transcript Database (Allium TDB, http://alliumtdb.kazusa.or.jp) was generated based on RNA-Seq data. The resulting assembled transcripts were functionally annotated, revealing 50 unigenes involved in saponin biosynthesis. Differential gene expression (DGE) analyses of AA and MALs as compared with FF (as a control) revealed a strong up-regulation of the saponin downstream pathway, including cytochrome P450, glycosyltransferase, and beta-glucosidase in chromosome 2A. An understanding of the saponin compounds and biosynthesis-related genes would facilitate the development of plants with unique saponin content and, subsequently, improved disease resistance.

Mapping and Identifying a Candidate Gene (Bnmfs) for Female-Male Sterility through Whole-Genome Resequencing and RNA-Seq in Rapeseed (Brassica napus L.)

In oilseed crops, carpel and stamen development play vital roles in pollination and rapeseed yield, but the genetic mechanisms underlying carpel and stamen development remain unclear. Herein, a male- and female-sterile mutant was obtained in offspring of a (Brassica napus cv. Qingyou 14) × (Qingyou 14 × B. rapa landrace Dahuang) cross. Subsequently, F₂-F₉ populations were generated through selfing of the heterozygote plants among the progeny of each generation. The male- and female-sterility exhibited stable inheritance in successive generations and was controlled by a recessive gene. The mutant kept the same chromosome number (2n = 38) as B. napus parent but showed abnormal meiosis for male and female. One candidate gene for the sterility was identified by simple sequence repeat (SSR) and insertion deletion length polymorphism (InDel) markers in F₇-F₉ plants, and whole-genome resequencing with F₈ pools and RNA sequencing with F₉ pools. Whole-genome resequencing found three candidate intervals (35.40-35.68, 35.74-35.75, and 45.34-46.45 Mb) on chromosome C3 in B. napus and candidate region for Bnmfs was narrowed to approximately 1.11-Mb (45.34-46.45 M) by combining SSR and InDel marker analyses with whole-genome resequencing. From transcriptome profiling in 0-2 mm buds, all of the genes in the candidate interval were detected, and only two genes with significant differences (BnaC03g56670D and BnaC03g56870D) were revealed. BnaC03g56870D was a candidate gene that shared homology with the CYP86C4 gene of Arabidopsis thaliana. Quantitative reverse transcription (qRT)-PCR analysis showed that Bnmfs primarily functioned in flower buds. Thus, sequencing and expression analyses provided evidence that BnaC03g56870D was the candidate gene for male and female sterility in the B. napus mutant.

Genome-wide transcriptome analysis of female-sterile rice ovule shed light on its abortive mechanism

The comprehensive transcriptome analysis of rice female-sterile line and wild-type line ovule provides an important clue for exploring the regulatory network of the formation of rice fertile female gametophyte. Ovules are the female reproductive tissues of rice (Oryza sativa L.) and play a major role in sexual reproduction. To investigate the potential mechanism of rice female gametophyte fertility, we used RNA sequencing, combined with genetic subtraction, to compare the transcriptome of the ovules of a high-frequency female-sterile line (fsv1) and a rice wild-type line (Gui 99) during ovule development. Ovules were harvested at three developmental stages: ovule containing megaspore mother cell in meiosis process (stage 1), ovule containing functional megaspore in mitosis process (stage 2), and ovule containing mature female gametophyte (stage 3). Six cDNA libraries generated a total of 42.2 million high-quality clean reads that aligned with 30,204 genes. The comparison between the fsv1 and Gui 99 ovules identified a large number of differentially expressed genes (DEGs), i.e., 45, 495, and 932 DEGs at the three ovule developmental stages, respectively. From the comparison of the two rice lines, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and MapMan analyses indicated that a large number of DEGs associated with starch and sucrose metabolism, plant hormone signal transduction, protein modification and degradation, oxidative phosphorylation, and receptor kinase. These DEGs might play roles in ovule development and fertile female gametophyte formation. Many transcription factor genes and epigenetic-related genes also exhibit different expression patterns and significantly different expression levels in two rice lines during ovule development, which might provide important information regarding the abortive mechanism of the female gametophyte in rice.

Transcriptome Analysis of Stem and Globally Comparison with Other Tissues in Brassica napus

Brassica napus is one of the most important oilseed crops in the world. However, there is currently no enough stem transcriptome information and comparative transcriptome analysis of different tissues, which impedes further functional genomics research on B. napus. In this study, the stem transcriptome of B. napus was characterized by RNA-seq technology. Approximately 13.4 Gb high-quality clean reads with an average length of 100 bp were generated and used for comparative transcriptome analysis with the existing transcriptome sequencing data of roots, leaves, flower buds, and immature embryos of B. napus. All the transcripts were annotated against GO and KEGG databases. The common genes in five tissues, differentially expressed genes (DEGs) of the common genes between stems and other tissues, and tissue-specific genes were detected, and the main biochemical activities and pathways implying the common genes, DEGs and tissue-specific genes were investigated. Accordingly, the common transcription factors (TFs) in the five tissues and tissue-specific TFs were identified, and a TFs-based regulation network between TFs and the target genes involved in 'Phenylpropanoid biosynthesis' pathway were constructed to show several important TFs and key nodes in the regulation process. Collectively, this study not only provided an available stem transcriptome resource in B. napus, but also revealed valuable comparative transcriptome information of five tissues of B. napus for future investigation on specific processes, functions and pathways.

Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq

The widespread and improper use of pyrethroid insecticides, such as deltamethrin, has resulted in the evolution of resistance in many mosquito species, including Culex pipiens pallens. With the development of high-throughput sequencing, it is possible to massively screen pyrethroid resistance-associated gene. In this study, we used Illumina-Solexa transcriptome sequencing to identify genes that are expressed differently in deltamethrin-susceptible and -resistant strains of Culex pipiens pallens as a critical knowledge base for further studies. A total of 4,961,197,620 base pairs and 55,124,418 reads were sequenced, mapped to the Culex quinquefasciatus genome and assembled into 17,679 known genes. We recorded 1826 significantly differentially expressed genes (DEGs). Among them, 1078 genes were up-regulated and 748 genes were down-regulated in the deltamethrin-resistant strain compared to -susceptible strain. These DEGs contained cytochrome P450 s, cuticle proteins, UDP-glucuronosyltransferases, lipases, serine proteases, heat shock proteins, esterases and others. Among the 1826 DEGs, we found that the transcriptional levels of CYP6AA9 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Moreover, the expression levels of the CYP6AA9 were significantly higher in the resistant strains than the susceptible strains in three different field populations. We further confirmed the association between the CYP6AA9 gene and deltamethrin resistance in mosquitoes by RNA interfering (RNAi). Altogether, we explored massive potential pyrethroid resistance-associated genes and demonstrated that CYP6AA9 participated in the pyrethroid resistance in mosquitoes.