Hongju Jian

Identification of Rapeseed MicroRNAs Involved in Early Stage Seed Germination under Salt and Drought Stresses

Drought and salinity are severe and wide-ranging abiotic stresses that substantially affect crop germination, development and productivity, and seed germination is the first critical step in plant growth and development. To comprehensively investigate small-RNA targets and improve our understanding of miRNA-mediated post-transcriptional regulation networks during Brassica napus seed imbibition under drought and salt stresses, we constructed three small-RNA libraries from B. napus variety ZS11 embryos exposed to salt (200 mM NaCl, denoted “S”), drought (200 g L−1 PEG-6000, denoted “D”), and distilled water (denoted “CK”) during imbibition and sequenced them using an Illumina Genome Analyzer. A total of 11,528,557, 12,080,081, and 12,315,608 raw reads were obtained from the CK, D, and S libraries, respectively. Further analysis identified 85 known miRNAs belonging to 31 miRNA families and 882 novel miRNAs among the three libraries. Comparison of the D and CK libraries revealed significant down-regulation of six miRNA families, miR156, miR169, miR860, miR399, miR171, and miR395, whereas only miR172 was significantly up-regulated. In contrast, comparison of the S library with the CK library showed significant down-regulation of only two miRNA families: miRNA393 and miRNA399. Putative targets for 336, 376, and 340 novel miRNAs were successfully predicted in the CK, D, and S libraries, respectively, and 271 miRNA families and 20 target gene families [including disease resistance protein (DIRP), drought-responsive family protein (DRRP), early responsive to dehydration stress protein (ERD), stress-responsive alpha-beta barrel domain protein (SRAP), and salt tolerance homolog2 (STH2)] were confirmed as being core miRNAs and genes involved in the seed imbibition response to salt and drought stresses. The sequencing results were partially validated by quantitative RT-PCR for both conserved and novel miRNAs as well as the predicted target genes. Our data suggest that diverse and complex miRNAs are involved in seed imbibition, indicating that miRNAs are involved in plant hormone regulation, and may play important roles during seed germination under salt- or drought-stress conditions.

Genome-Wide Analysis and Expression Profiling of the SUC and SWEET Gene Families of Sucrose Transporters in Oilseed Rape (Brassica napus L.)

Sucrose is the principal transported product of photosynthesis from source leaves to sink organs. SUTs/SUCs (sucrose transporters or sucrose carriers) and SWEETs (Sugars Will Eventually be Exported Transporters) play significant central roles in phloem loading and unloading. SUTs/SUCs and SWEETs are key players in sucrose translocation and are associated with crop yields. The SUT/SUC and SWEET genes have been characterized in several plant species, but a comprehensive analysis of these two gene families in oilseed rape has not yet been reported. In our study, 22 and 68 members of the SUT/SUCs and SWEET gene families, respectively, were identified in the oilseed rape (Brassica napus) genome through homology searches. An analysis of the chromosomal distribution, phylogenetic relationships, gene structures, motifs and the cis-acting regulatory elements in the promoters of BnSUC and BnSWEET genes were analyzed. Furthermore, we examined the expression of the 18 BnSUC and 16 BnSWEET genes in different tissues of “ZS11” and the expression of 9 BnSUC and 7 BnSWEET genes in “ZS11” under various conditions, including biotic stress (Sclerotinia sclerotiorum), abiotic stresses (drought, salt and heat), and hormone treatments (abscisic acid, auxin, cytokinin, brassinolide, gibberellin, and salicylic acid). In conclusion, our study provides the first comprehensive analysis of the oilseed rape SUC and SWEET gene families. Information regarding the phylogenetic relationships, gene structure and expression profiles of the SUC and SWEET genes in the different tissues of oilseed rape helps to identify candidates with potential roles in specific developmental processes. Our study advances our understanding of the important roles of sucrose transport in oilseed rape.

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.

Identification of microRNAs Actively Involved in Fatty Acid Biosynthesis in Developing Brassica napus Seeds Using High-Throughput Sequencing

Seed development has a critical role during the spermatophyte life cycle. In Brassica napus, a major oil crop, fatty acids are synthesized and stored in specific tissues during embryogenesis, and understanding the molecular mechanism underlying fatty acid biosynthesis during seed development is an important research goal. In this study, we constructed three small RNA libraries from early seeds at 14, 21, and 28 days after flowering (DAF) and used high-throughput sequencing to examine microRNA (miRNA) expression. A total of 85 known miRNAs from 30 families and 1160 novel miRNAs were identified, of which 24, including 5 known and 19 novel miRNAs, were found to be involved in fatty acid biosynthesis.bna-miR156b, bna-miR156c, bna-miR156g, novel_mir_1706, novel_mir_1407, novel_mir_173, and novel_mir_104 were significantly down-regulated at 21 DAF and 28 DAF, whereas bna-miR159, novel_mir_1081, novel_mir_19 and novel_mir_555 were significantly up-regulated. In addition, we found that some miRNAs regulate functional genes that are directly involved in fatty acid biosynthesis and that other miRNAs regulate the process of fatty acid biosynthesis by acting on a large number of transcription factors. The miRNAs and their corresponding predicted targets were partially validated by quantitative RT-PCR. Our data suggest that diverse and complex miRNAs are involved in the seed development process and that miRNAs play important roles in fatty acid biosynthesis during seed development.

Genome-wide analysis and expression profiling of the GRF gene family in oilseed rape (Brassica napus L.)

Growth regulating-factors (GRFs) are plant-specific transcription factors that help regulate plant growth and development. Genome-wide identification and evolutionary analyses of GRF gene families have been performed in Arabidopsis thaliana, Zea mays, Oryza sativa, and Brassica rapa, but a comprehensive analysis of the GRF gene family in oilseed rape (Brassica napus) has not yet been reported. In the current study, we identified 35 members of the BnGRF family in B. napus. We analyzed the chromosomal distribution, phylogenetic relationships (Bayesian Inference and Neighbor Joining method), gene structures, and motifs of the BnGRF family members, as well as the cis-acting regulatory elements in their promoters. We also analyzed the expression patterns of 15 randomly selected BnGRF genes in various tissues and in plant varieties with different harvest indices and gibberellic acid (GA) responses. The expression levels of BnGRFs under GA treatment suggested the presence of possible negative feedback regulation. The evolutionary patterns and expression profiles of BnGRFs uncovered in this study increase our understanding of the important roles played by these genes in oilseed rape.

Genome-Wide Analysis of the PYL Gene Family and Identification of PYL Genes That Respond to Abiotic Stress in Brassica napus

Abscisic acid (ABA) is an endogenous phytohormone that plays important roles in the regulation of plant growth, development, and stress responses. The pyrabactin resistance 1-like (PYR/PYL) protein is a core regulatory component of ABA signaling networks in plants. However, no details regarding this family in Brassica napus are available. Here, 46 PYLs were identified in the B. napus genome. Based on phylogenetic analysis, BnPYR1 and BnPYL1-3 belong to subfamily I, BnPYL7-10 belong to subfamily II, and BnPYL4-6 and BnPYL11-13 belong to subfamily III. Analysis of BnPYL conserved motifs showed that every subfamily contained four common motifs. By predicting cis-elements in the promoters, we found that all BnPYL members contained hormone- and stress-related elements and that expression levels of most BnPYLs were relatively higher in seeds at the germination stage than those in other organs or at other developmental stages. Gene Ontology (GO) enrichment showed that BnPYL genes mainly participate in responses to stimuli. To identify crucial PYLs mediating the response to abiotic stress in B. napus, expression changes in 14 BnPYL genes were determined by quantitative real-time RT-PCR after drought, heat, and salinity treatments, and identified BnPYR1-3, BnPYL1-2, and BnPYL7-2 in respond to abiotic stresses. The findings of this study lay a foundation for further investigations of PYL genes in B. napus.

Genetic and transcriptomic analyses of lignin- and lodging-related traits in Brassica napus

Key messageCandidate genes associated with lignin and lodging traits were identified by combining phenotypic, genotypic, and gene expression data in B. napus.AbstractBrassica napus is one of the world’s most important oilseed crops, but its yield can be dramatically reduced by lodging, bending, and falling of its vertical stems. Lignin has been shown to contribute to stem mechanical strength. In this study, we found that the syringyl/guaiacyl (S/G) monolignol ratio exhibits a significant negative correlation with disease and lodging resistance. A total of 92 and 50 SNP and SSR loci, respectively, were found to be significantly associated with five traits, breaking force, breaking strength, lodging coefficient, acid detergent lignin content, and the S/G monolignol ratio using GWAS. To identify novel genes involved in lignin biosynthesis, transcriptome sequencing of high- (H) and low (L)-ADL content accessions was performed. The up-regulated genes were mainly involved in glycoside catabolic processes (especially glucosinolate catabolism) and cell wall biogenesis, while down-regulated genes were involved in glucosinolate biosynthesis, indicating that crosstalk exists between glucosinolate metabolic processes and lignin biosynthesis. Integrating this differential expression with the GWAS analysis, we identified four candidate genes regulating lignin, including glycosyl hydrolase (BnaA01g00480D), CYT1 (BnaA04g22820D), and two encoding transcription factors, SHINE1 (ERF family) and DAR6 (LIM family). This study provides insight into the genetic control of lodging and lignin in B. napus.

Screening of Candidate Leaf Morphology Genes by Integration of QTL Mapping and RNA Sequencing Technologies in Oilseed Rape (Brassica napus L.)

Leaf size and shape play important roles in agronomic traits, such as yield, quality and stress responses. Wide variations in leaf morphological traits exist in cultivated varieties of many plant species. By now, the genetics of leaf shape and size have not been characterized in Brassica napus. In this study, a population of 172 recombinant inbred lines (RILs) was used for quantitative trait locus (QTL) analysis of leaf morphology traits. Furthermore, fresh young leaves of extreme lines with more leaf lobes (referred to as ‘A’) and extreme lines with fewer lobes (referred to as ‘B’) selected from the RIL population and leaves of dissected lines (referred to as ‘P’) were used for transcriptional analysis. A total of 31 QTLs for the leaf morphological traits tested in this study were identified on 12 chromosomes, explaining 5.32–39.34% of the phenotypic variation. There were 8, 6, 2, 5, 8, and 2 QTLs for PL (petiole length), PN (lobe number), LW (lamina width), LL (Lamina length), LL/LTL (the lamina size ratio) and LTL (leaf total length), respectively. In addition, 74, 1,166 and 1,272 differentially expressed genes (DEGs) were identified in ‘A vs B’, ‘A vs P’ and ‘B vs P’ comparisons, respectively. The Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used to predict the functions of these DEGs. Gene regulators of leaf shape and size, such as ASYMMETRIC LEAVES 2, gibberellin 20-oxidase 3, genes encoding gibberellin-regulated family protein, genes encoding growth-regulating factor and KNOTTED1-like homeobox were also detected in DEGs. After integrating the QTL mapping and RNA sequencing data, 33 genes, including a gene encoding auxin-responsive GH3 family protein and a gene encoding sphere organelles protein-related gene, were selected as candidates that may control leaf shape. Our findings should be valuable for studies of the genetic control of leaf morphological trait regulation in B. napus.

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.