Integrative RNA- and miRNA-Profile Analysis Reveals a Likely Role of BR and Auxin Signaling in Branch Angle Regulation of B. napus

Homoeologous exchanges contribute to branch angle variations in rapeseed: Insights from transcriptome, QTL-seq and gene functional analysis

Branch angle (BA) is a critical morphological trait that significantly influences planting density, light interception and ultimately yield in plants. Despite its importance, the regulatory mechanism governing BA in rapeseed remains poorly understood. In this study, we generated 109 transcriptome data sets for 37 rapeseed accessions with divergent BA phenotypes. Relative to adaxial branch segments, abaxial segments accumulated higher levels of auxin and exhibited lower expression of six TCP1 homologues and one GA20ox3. A co-expression network analysis identified two modules highly correlated with BA. The modules contained homologues to known BA control genes, such as FUL, YUCCA6, TCP1 and SGR3. Notably, a homoeologous exchange (HE), occurring at the telomeres of A09, was prevalent in large BA accessions, while an A02-C02 HE was common in small BA accessions. In their corresponding regions, these HEs explained the formation of hub gene hotspots in the two modules. QTL-seq analysis confirmed that the presence of a large A07-C06 HE (~8.1 Mb) was also associated with a small BA phenotype, and BnaA07.WRKY40.b within it was predicted as candidate gene. Overexpressing BnaA07.WRKY40.b in rapeseed increased BA by up to 20°, while RNAi- and CRISPR-mediated mutants (BnaA07.WRKY40.b and BnaC06.WRKY40.b) exhibited decreased BA by up to 11.4°. BnaA07.WRKY40.b was exclusively localized to the nucleus and exhibited strong expression correlations with many genes related to gravitropism and plant architecture. Taken together, our study highlights the influence of HEs on rapeseed plant architecture and confirms the role of WRKY40 homologues as novel regulators of BA.

Identification of a major QTL and candidate genes analysis for branch angle in rapeseed (Brassica napus L.) using QTL-seq and RNA-seq

Introduction

Branching angle is an essential trait in determining the planting density of rapeseed (Brassica napus L.) and hence the yield per unit area. However, the mechanism of branching angle formation in rapeseed is not well understood.

Methods

In this study, two rapeseed germplasm with extreme branching angles were used to construct an F2 segregating population; then bulked segregant analysis sequencing (BSA-seq) and quantitative trait loci (QTL) mapping were utilized to localize branching anglerelated loci and combined with transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qPCR) for candidate gene mining.

Results and discussion

A branching angle-associated quantitative trait loci (QTL) was mapped on chromosome C3 (C3: 1.54-2.65 Mb) by combining BSA-seq as well as traditional QTL mapping. A total of 54 genes had SNP/Indel variants within the QTL interval were identified. Further, RNA-seq of the two parents revealed that 12 of the 54 genes were differentially expressed between the two parents. Finally, we further validated the differentially expressed genes using qPCR and found that six of them presented consistent differential expression in all small branching angle samples and large branching angles, and thus were considered as candidate genes related to branching angles in rapeseed. Our results introduce new candidate genes for the regulation of branching angle formation in rapeseed, and provide an important reference for the subsequent exploration of its formation mechanism.

A multi-trait GWAS-based genetic association network controlling soybean architecture and seed traits

Ideal plant architecture is essential for enhancing crop yields. Ideal soybean (Glycine max) architecture encompasses an appropriate plant height, increased node number, moderate seed weight, and compact architecture with smaller branch angles for growth under high-density planting. However, the functional genes regulating plant architecture are far not fully understood in soybean. In this study, we investigated the genetic basis of 12 agronomic traits in a panel of 496 soybean accessions with a wide geographical distribution in China. Analysis of phenotypic changes in 148 historical elite soybean varieties indicated that seed-related traits have mainly been improved over the past 60 years, with targeting plant architecture traits having the potential to further improve yields in future soybean breeding programs. In a genome-wide association study (GWAS) of 12 traits, we detected 169 significantly associated loci, of which 61 overlapped with previously reported loci and 108 new loci. By integrating the GWAS loci for different traits, we constructed a genetic association network and identified 90 loci that were associated with a single trait and 79 loci with pleiotropic effects. Of these 79 loci, 7 hub-nodes were strongly linked to at least three related agronomic traits. qHub_5, containing the previously characterized Determinate 1 (Dt1) locus, was associated not only with plant height and node number (as determined previously), but also with internode length and pod range. Furthermore, we identified qHub_7, which controls three branch angle-related traits; the candidate genes in this locus may be beneficial for breeding soybean with compact architecture. These findings provide insights into the genetic relationships among 12 important agronomic traits in soybean. In addition, these studies uncover valuable loci for further functional gene studies and will facilitate molecular design breeding of soybean architecture.

Tissue-Specific Transcriptome and Metabolome Analysis Reveals the Response Mechanism of Brassica napus to Waterlogging Stress

During the growth period of rapeseed, if there is continuous rainfall, it will easily lead to waterlogging stress, which will seriously affect the growth of rapeseed. Currently, the mechanisms of rapeseed resistance to waterlogging stress are largely unknown. In this study, the rapeseed (Brassica napus) inbred lines G230 and G218 were identified as waterlogging-tolerant rapeseed and waterlogging-sensitive rapeseed, respectively, through a potted waterlogging stress simulation and field waterlogging stress experiments. After six days of waterlogging stress at the seedling stage, the degree of leaf aging and root damage of the waterlogging-tolerant rapeseed G230 were lower than those of the waterlogging-sensitive rapeseed G218. A physiological analysis showed that waterlogging stress significantly increased the contents of malondialdehyde, soluble sugar, and hydrogen peroxide in rape leaves and roots. The transcriptomic and metabolomic analysis showed that the differential genes and the differential metabolites of waterlogging-tolerant rapeseed G230 were mainly enriched in the metabolic pathways, biosynthesis of secondary metabolites, flavonoid biosynthesis, and vitamin B6 metabolism. Compared to G218, the expression levels of some genes associated with flavonoid biosynthesis and vitamin B metabolism were higher in G230, such as CHI, DRF, LDOX, PDX1.1, and PDX2. Furthermore, some metabolites involved in flavonoid biosynthesis and vitamin B6 metabolism, such as naringenin and epiafzelechin, were significantly up-regulated in leaves of G230, while pyridoxine phosphate was only significantly down-regulated in roots and leaves of G218. Furthermore, foliar spraying of vitamin B6 can effectively improve the tolerance to waterlogging of G218 in the short term. These results indicate that flavonoid biosynthesis and vitamin B6 metabolism pathways play a key role in the waterlogging tolerance and hypoxia stress resistance of Brassica napus and provide new insights for improving the waterlogging tolerance and cultivating waterlogging-tolerant rapeseed varieties.

Auxin Transporter OsPIN1b, a Novel Regulator of Leaf Inclination in Rice (Oryza sativa L.)

Leaf inclination is one of the most important components of the ideal architecture, which effects yield gain. Leaf inclination was shown that is mainly regulated by brassinosteroid (BR) and auxin signaling. Here, we reveal a novel regulator of leaf inclination, auxin transporter OsPIN1b. Two CRISPR-Cas9 homozygous mutants, ospin1b-1 and ospin1b-2, with smaller leaf inclination compared to the wild-type, Nipponbare (WT/NIP), while overexpression lines, OE-OsPIN1b-1 and OE-OsPIN1b-2 have opposite phenotype. Further cell biological observation showed that in the adaxial region, OE-OsPIN1b-1 has significant bulge compared to WT/NIP and ospin1b-1, indicating that the increase in the adaxial cell division results in the enlarging of the leaf inclination in OE-OsPIN1b-1. The OsPIN1b was localized on the plasma membrane, and the free IAA contents in the lamina joint of ospin1b mutants were significantly increased while they were decreased in OE-OsPIN1b lines, suggesting that OsPIN1b might action an auxin transporter such as AtPIN1 to alter IAA content and leaf inclination. Furthermore, the OsPIN1b expression was induced by exogenous epibrassinolide (24-eBL) and IAA, and ospin1b mutants are insensitive to BR or IAA treatment, indicating that the effecting leaf inclination is regulated by OsPIN1b. This study contributes a new gene resource for molecular design breeding of rice architecture.

Functional Characterization of the Cystine-Rich-Receptor-like Kinases (CRKs) and Their Expression Response to Sclerotinia sclerotiorum and Abiotic Stresses in Brassica napus

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that bind to the calcium ion to regulate stress-signaling and plant development-related pathways, as indicated by several pieces of evidence. However, the CRK gene family hasn’t been inadequately examined in Brassica napus. In our study, 27 members of the CRK gene family were identified in Brassica napus, which are categorized into three phylogenetic groups and display synteny relationship to the Arabidopsis thaliana orthologs. All the CRK genes contain highly conserved N-terminal PKINASE domain; however, the distribution of motifs and gene structure were variable conserved. The functional divergence analysis between BnaCRK groups indicates a shift in evolutionary rate after duplication events, demonstrating that BnaCRKs might direct a specific function. RNA-Seq datasets and quantitative real-time PCR (qRT-PCR) exhibit the complex expression profile of the BnaCRKs in plant tissues under multiple stresses. Nevertheless, BnaA06CRK6-1 and BnaA08CRK8 from group B were perceived to play a predominant role in the Brassica napus stress signaling pathway in response to drought, salinity, and Sclerotinia sclerotiorum infection. Insights gained from this study improve our knowledge about the Brassica napus CRK gene family and provide a basis for enhancing the quality of rapeseed.

Multiple Functions of MiRNAs in Brassica napus L

The worldwide climate changes every year due to global warming, waterlogging, drought, salinity, pests, and pathogens, impeding crop productivity. Brassica napus is one of the most important oil crops in the world, and rapeseed oil is considered one of the most health-beneficial edible vegetable oils. Recently, miRNAs have been found and confirmed to control the expression of targets under disruptive environmental conditions. The mechanism is through the formation of the silencing complex that mediates post-transcriptional gene silencing, which pairs the target mRNA and target cleavage and/or translation inhibition. However, the functional role of miRNAs and targets in B. napus is still not clarified. This review focuses on the current knowledge of miRNAs concerning development regulation and biotic and abiotic stress responses in B. napus. Moreover, more strategies for miRNA manipulation in plants are discussed, along with future perspectives, and the enormous amount of transcriptome data available provides cues for miRNA functions in B. napus. Finally, the construction of the miRNA regulatory network can lead to the significant development of climate change-tolerant B. napus through miRNA manipulation.

Genetic dissection of branch architecture in oilseed rape (Brassica napus L.) germplasm

Branch architecture is an important factor influencing rapeseed planting density, mechanized harvest, and yield. However, its related genes and regulatory mechanisms remain largely unknown. In this study, branch angle (BA) and branch dispersion degree (BD) were used to evaluate branch architecture. Branch angle exhibited a dynamic change from an increase in the early stage to a gradual decrease until reaching a stable state. Cytological analysis showed that BA variation was mainly due to xylem size differences in the vascular bundle of the branch junction. The phenotypic analysis of 327 natural accessions revealed that BA in six environments ranged from 24.3° to 67.9°, and that BD in three environments varied from 4.20 cm to 21.4 cm, respectively. A total of 115 significant loci were detected through association mapping in three models (MLM, mrMLM, and FarmCPU), which explained 0.53%-19.4% of the phenotypic variations. Of them, 10 loci were repeatedly detected in different environments and models, one of which qBAD.A03-2 was verified as a stable QTL using a secondary segregation population. Totally, 1066 differentially expressed genes (DEGs) were identified between branch adaxial- and abaxial- sides from four extremely large or small BA/BD accessions through RNA sequencing. These DEGs were significantly enriched in the pathways related to auxin biosynthesis and transport as well as cell extension such as indole alkaloid biosynthesis, other glycan degradation, and fatty acid elongation. Four known candidate genes BnaA02g16500D (PIN1), BnaA03g10430D (PIN2), BnaC03g06250D (LAZY1), and BnaC06g20640D (ARF17) were identified by both GWAS and RNA-seq, all of which were involved in regulating the asymmetric distribution of auxins. Our identified association loci and candidate genes provide a theoretical basis for further study of gene cloning and genetic improvement of branch architecture.

Bulk RNA-Seq Analysis Reveals Differentially Expressed Genes Associated with Lateral Branch Angle in Peanut

Lateral branch angle (LBA), or branch habit, is one of the most important agronomic traits in peanut. To date, the underlying molecular mechanisms of LBA have not been elucidated in peanut. To acquire the differentially expressed genes (DEGs) related to LBA, a TI population was constructed through the hybridization of a bunch-type peanut variety Tifrunner and prostrate-type Ipadur. We report the identification of DEGs related to LBA by sequencing two RNA pools, which were composed of 45 F₃ lines showing an extreme opposite bunch and prostrate phenotype. We propose to name this approach Bulk RNA-sequencing (BR-seq) as applied to several plant species. Through BR-seq analysis, a total of 3083 differentially expressed genes (DEGs) were identified, including 13 gravitropism-related DEGs, 22 plant hormone-related DEGs, and 55 transcription factors-encoding DEGs. Furthermore, we also identified commonly expressed alternatively spliced (AS) transcripts, of which skipped exon (SE) and retained intron (RI) were most abundant in the prostrate and bunch-type peanut. AS isoforms between prostrate and bunch peanut highlighted important clues to further understand the post-transcriptional regulatory mechanisms of branch angle regulation. Our findings provide not only important insights into the landscape of the regulatory pathway involved in branch angle formation but also present practical information for peanut molecular breeding in the future.

Genome-Wide Differential DNA Methylation and miRNA Expression Profiling Reveals Epigenetic Regulatory Mechanisms Underlying Nitrogen-Limitation-Triggered Adaptation and Use Efficiency Enhancement in Allotetraploid Rapeseed

Improving crop nitrogen (N) limitation adaptation (NLA) is a core approach to enhance N use efficiency (NUE) and reduce N fertilizer application. Rapeseed has a high demand for N nutrients for optimal plant growth and seed production, but it exhibits low NUE. Epigenetic modification, such as DNA methylation and modification from small RNAs, is key to plant adaptive responses to various stresses. However, epigenetic regulatory mechanisms underlying NLA and NUE remain elusive in allotetraploid B. napus. In this study, we identified overaccumulated carbohydrate, and improved primary and lateral roots in rapeseed plants under N limitation, which resulted in decreased plant nitrate concentrations, enhanced root-to-shoot N translocation, and increased NUE. Transcriptomics and RT-qPCR assays revealed that N limitation induced the expression of NRT1.1, NRT1.5, NRT1.7, NRT2.1/NAR2.1, and Gln1;1, and repressed the transcriptional levels of CLCa, NRT1.8, and NIA1. High-resolution whole genome bisulfite sequencing characterized 5094 differentially methylated genes involving ubiquitin-mediated proteolysis, N recycling, and phytohormone metabolism under N limitation. Hypermethylation/hypomethylation in promoter regions or gene bodies of some key N-metabolism genes might be involved in their transcriptional regulation by N limitation. Genome-wide miRNA sequencing identified 224 N limitation-responsive differentially expressed miRNAs regulating leaf development, amino acid metabolism, and plant hormone signal transduction. Furthermore, degradome sequencing and RT-qPCR assays revealed the miR827-NLA pathway regulating limited N-induced leaf senescence as well as the miR171-SCL6 and miR160-ARF17 pathways regulating root growth under N deficiency. Our study provides a comprehensive insight into the epigenetic regulatory mechanisms underlying rapeseed NLA, and it will be helpful for genetic engineering of NUE in crop species through epigenetic modification of some N metabolism-associated genes.

Identification of Bna.IAA7.C05 as allelic gene for dwarf mutant generated from tissue culture in oilseed rape

BACKGROUND

Plant height is one of the most important agronomic traits in many crops due to its influence on lodging resistance and yield performance. Although progress has been made in the use of dwarfing genes in crop improvement, identification of new dwarf germplasm is still of significant interest for breeding varieties with increased yield.

RESULTS

Here we describe a dominant, dwarf mutant G7 of Brassica napus with down-curved leaves derived from tissue culture. To explore the genetic variation responsible for the dwarf phenotype, the mutant was crossed to a conventional line to develop a segregating F₂ population. Bulks were formed from plants with either dwarf or conventional plant height and subjected to high throughput sequencing analysis via mutation mapping (MutMap). The dwarf mutation was mapped to a 0.6 Mb interval of B. napus chromosome C05. Candidate gene analysis revealed that one SNP causing an amino acid change in the domain II of Bna.IAA7.C05 may contribute to the dwarf phenotype. This is consistent with the phenotype of a gain-of-function indole-3-acetic acid (iaa) mutant in Bna.IAA7.C05 reported recently. GO and KEGG analysis of RNA-seq data revealed the down-regulation of auxin related genes, including many other IAA and small up regulated response (SAUR) genes, in the dwarf mutant.

CONCLUSION

Our studies characterize a new allele of Bna.IAA7.C05 responsible for the dwarf mutant generated from tissue culture. This may provide a valuable genetic resource for breeding for lodging resistance and compact plant stature in B. napus.

The miRNA-mRNA Networks Involving Abnormal Energy and Hormone Metabolisms Restrict Tillering in a Wheat Mutant dmc

Tillers not only determine plant architecture but also influence crop yield. To explore the miRNA regulatory network restraining tiller development in a dwarf-monoculm wheat mutant (dmc) derived from Guomai 301 (wild type, WT), we employed miRNome and transcriptome integrative analysis, real-time qRT-PCR, histochemistry, and determinations of the key metabolites and photosynthesis parameters. A total of 91 differentially expressed miRNAs (DEMs) were identified between dmc and WT. Among them, 40 key DEMs targeted 45 differentially expressed genes (DEGs) including the key DEGs encode growth-regulating factors (GRF), auxin response factors (ARF), and other proteins involved in the metabolisms of hormones and carbohydrates, etc. Compared with WT, both the chlorophyll contents and the photosynthesis rate were lower in dmc. The contents of glucose, sucrose, fructose, and maltose were lower in dmc. The contents of auxin (IAA) and zeatin (ZA) were significantly lower, but gibberellin (GA) was significantly higher in the tiller tissues of dmc. This research demonstrated that the DEMs regulating hormone and carbohydrate metabolisms were important causes for dmc to not tiller. A primary miRNA-mRNA regulatory model for dmc tillering was established. The lower photosynthesis rate, insufficient energy, and abnormal hormone metabolisms restrict tillering in dmc.

Joint RNA-Seq and miRNA Profiling Analyses to Reveal Molecular Mechanisms in Regulating Thickness of Pod Canopy in Brassica napus

Oilseed rape (Brassica napus) is the second largest oilseed crop worldwide. As an architecture component of B. napus, thickness of pod canopy (TPC) plays an important role in yield formation, especially under high-density cultivation conditions. However, the mechanisms underlying the regulation of TPC remain unclear. RNA and microRNA (miRNA) profiling of two groups of B. napus lines with significantly different TPC at the bolting with a tiny bud stage revealed differential expressions of numerous genes involved in nitrogen-related pathways. Expression of several nitrogen-related response genes, including ASP5, ASP2, ASN3, ATCYSC1, PAL2, APT2, CRTISO, and COX15, was dramatically changed in the thick TPC lines compared to those in the thin TPC lines. Differentially expressed miRNAs also included many involved in nitrogen-related pathways. Expression of most target genes was negatively associated with corresponding miRNAs, such as miR159, miR6029, and miR827. In addition, 12 (including miR319, miR845, and miR158) differentially expressed miRNAs between two plant tissues sampled (stem apex and flower bud) were identified, implying that they might have roles in determining overall plant architecture. These results suggest that nitrogen signaling may play a pivotal role in regulating TPC in B. napus.

Auxin Metabolism Controls Developmental Decisions in Land Plants

Unlike animals, whose body plans are set during embryo development, plants maintain the ability to initiate new organs throughout their life cycle. Auxin is a key regulator of almost all aspects of plant development, including morphogenesis and adaptive responses. Cellular auxin concentrations influence whether a cell will divide, grow, or differentiate, thereby contributing to organ formation, growth, and ultimately plant shape. Auxin gradients are established and maintained by a tightly regulated interplay between metabolism, signalling, and transport. Auxin is synthesised, stored, and inactivated by a multitude of parallel pathways that are all tightly regulated. Here we summarise the remarkable progress that has been achieved in identifying some key components of these pathways and the genetic complexity underlying their precise regulation.

The Gene Structure and Expression Level Changes of the GH3 Gene Family in Brassica napus Relative to Its Diploid Ancestors

The GH3 gene family plays a vital role in the phytohormone-related growth and developmental processes. The effects of allopolyploidization on GH3 gene structures and expression levels have not been reported. In this study, a total of 38, 25, and 66 GH3 genes were identified in Brassica rapa (A_rA_r), Brassica oleracea (C_oC_o), and Brassica napus (A_nAC_nC_n), respectively. BnaGH3 genes were unevenly distributed on chromosomes with 39 on A_n and 27 on C_n, in which six BnaGH3 genes may appear as new genes. The whole genome triplication allowed the GH3 gene family to expand in diploid ancestors, and allopolyploidization made the GH3 gene family re-expand in B. napus. For most BnaGH3 genes, the exon-intron compositions were similar to diploid ancestors, while the cis-element distributions were obviously different from its ancestors. After allopolyploidization, the expression patterns of GH3 genes from ancestor species changed greatly in B. napus, and the orthologous gene pairs between A_n/A_r and C_n/C_o had diverged expression patterns across four tissues. Our study provides a comprehensive analysis of the GH3 gene family in B. napus, and these results could contribute to identifying genes with vital roles in phytohormone-related growth and developmental processes.

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.

Integrated physiologic, genomic and transcriptomic strategies involving the adaptation of allotetraploid rapeseed to nitrogen limitation

BACKGROUND

Nitrogen (N) is a macronutrient that is essential for optimal plant growth and seed yield. Allotetraploid rapeseed (A_nA_nC_nC_n, 2n = 4x = 38) has a higher requirement for N fertilizers whereas exhibiting a lower N use efficiency (NUE) than cereal crops. N limitation adaptation (NLA) is pivotal for enhancing crop NUE and reducing N fertilizer use in yield production. Therefore, revealing the genetic and molecular mechanisms underlying NLA is urgent for the genetic improvement of NUE in rapeseed and other crop species with complex genomes.

RESULTS

In this study, we integrated physiologic, genomic and transcriptomic analyses to comprehensively characterize the adaptive strategies of oilseed rape to N limitation stresses. Under N limitations, we detected accumulated anthocyanin, reduced nitrate (NO₃^-) and total N concentrations, and enhanced glutamine synthetase activity in the N-starved rapeseed plants. High-throughput transcriptomics revealed that the pathways associated with N metabolism and carbon fixation were highly over-represented. The expression of the genes that were involved in efficient N uptake, translocation, remobilization and assimilation was significantly altered. Genome-wide identification and molecular characterization of the microR827-NLA1-NRT1.7 regulatory circuit indicated the crucial role of the ubiquitin-mediated post-translational pathway in the regulation of rapeseed NLA. Transcriptional analysis of the module genes revealed their significant functional divergence in response to N limitations between allotetraploid rapeseed and the model Arabidopsis. Association analysis in a rapeseed panel comprising 102 genotypes revealed that BnaC5.NLA1 expression was closely correlated with the rapeseed low-N tolerance.

CONCLUSIONS

We identified the physiologic and genome-wide transcriptional responses of oilseed rape to N limitation stresses, and characterized the global members of the BnamiR827-BnaNLA1s-BnaNRT1.7s regulatory circuit. The transcriptomics-assisted gene co-expression network analysis accelerates the rapid identification of central members within large gene families of plant species with complex genomes. These findings would enhance our comprehensive understanding of the physiologic responses, genomic adaptation and transcriptomic alterations of oilseed rape to N limitations and provide central gene resources for the genetic improvement of crop NLA and NUE.

Physiological studies and genome-wide microRNA profiling of cold-stressed Brassica napus

Temperature extremes, including cold, adversely impact plant growth and development. Plant responses to cold stress (CS) are regulated at both transcriptional and post-transcriptional levels. MicroRNAs (miRNAs), small non-coding RNAs, are known to be involved in post-transcriptional regulation of various developmental processes and metal stress in Brassica napus L. (canola), however, their role in response to CS is largely unknown. In this study, changes in various physiological parameters and endogenous abundance of miRNAs were characterized in spring canola seedlings (DH12075) exposed to 4 °C for 0-48 h. Cold stress induced electrolyte leakage, increased the levels of malondialdheyde and antioxidant enzymes and reduced photosynthetic efficiency. Using small RNA sequencing, 70 known and 126 novel miRNAs were identified in CS leaf tissues and among these, 25 known and 104 novel miRNAs were differentially expressed. Quantitative real-time (qRT) PCR analysis of eight selected miRNAs confirmed their CS responsiveness. Furthermore, the expression of six out of eight miRNAs exhibited an opposite trend in a winter variety of canola, 'Mendel', when compared to 'DH12075'. This first study on the B. napus miRNAome provides a framework for further functional analysis of these miRNAs and their targets in response to CS which may contribute towards the future development of cold resilient crops.

Genome-Wide Identification of MicroRNAs in Response to Cadmium Stress in Oilseed Rape (Brassica napus L.) Using High-Throughput Sequencing

MicroRNAs (miRNAs) have important roles in regulating stress-response genes in plants. However, identification of miRNAs and the corresponding target genes that are induced in response to cadmium (Cd) stress in Brassica napus remains limited. In the current study, we sequenced three small-RNA libraries from B. napus after 0 days, 1 days, and 3 days of Cd treatment. In total, 44 known miRNAs (belonging to 27 families) and 103 novel miRNAs were identified. A comprehensive analysis of miRNA expression profiles found 39 differentially expressed miRNAs between control and Cd-treated plants; 13 differentially expressed miRNAs were confirmed by qRT-PCR. Characterization of the corresponding target genes indicated functions in processes including transcription factor regulation, biotic stress response, ion homeostasis, and secondary metabolism. Furthermore, we propose a hypothetical model of the Cd-response mechanism in B. napus. Combined with qRT-PCR confirmation, our data suggested that miRNAs were involved in the regulations of TFs, biotic stress defense, ion homeostasis and secondary metabolism synthesis to respond Cd stress in B. napus.

Comparative transcriptome analysis provides new insights into erect and prostrate growth in bermudagrass (Cynodon dactylon L.)

Bermudagrass (Cynodon dactylon L.) is a prominent warm-season turf and forage grass species with multiple applications. In most C. dactylon cultivars and accessions, erect-growing stems (shoot) and prostrate-growing stems (stolon) often coexist. These two types of stems are both formed through tillering but grow in two directions with different tiller angles. Elucidating the mechanism of tiller angle regulation in bermudagrass could provide important clues to breed cultivars with different plant architectural features for diverse usage. In this study, we compared the stem internode transcriptome of two bermudagrass wild accessions with extremely different tiller angles and stem growth directions. A total of 2088 and 12,141 unigenes were preferentially expressed in prostrate-growing wild accession C792 and erect-growing wild accession C793, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology-based Annotation System (KOBAS) analyses further indicated that light- and gravity-responsive genes were enriched in accession C792, whereas lignin synthesis-related genes were enriched in accession C793, which well explains the difference in lignification of vascular bundles and mechanical tissues in the two accessions. These results not only expand our understanding of the genetic control of tiller angle and stem growth direction in bermudagrass but also provide insight for future molecular breeding of C. dactylon and other turfgrass species with different plant architectures.