A Genome-Wide Identification and Expression Pattern of LMCO Gene Family from Turnip (Brassica rapa L.) under Various Abiotic Stresses

Laccase-like multi-copper oxidases (LMCOs) are a group of enzymes involved in the oxidation of numerous substrates. Recently, these enzymes have become extremely popular due to their practical applications in various fields of biology. LMCOs generally oxidize various substrates by linking four-electron reduction of the final acceptor, molecular oxygen (O₂), to water. Multi-copper oxidases related to laccase are extensively distributed as multi-gene families in the genome sequences of higher plants. The current study thoroughly investigated the LMCO gene family (Br-Lac) and its expression pattern under various abiotic stresses in B. rapa L. A total of 18 Br-Lac gene family members located on five different chromosomes were identified. Phylogenetic analysis classified the documented Br-Lac genes into seven groups: Group-I (four genes), Group-II (nine genes), Group-III (eight genes), Group-IV (four genes), Group-V (six genes), and Group-VI and Group-VII (one gene each). The key features of gene structure and responsive motifs shared the utmost resemblance within the same groups. Additionally, a divergence study also assessed the evolutionary features of Br-Lac genes. The anticipated period of divergence ranged from 12.365 to 39.250 MYA (million years ago). We also identified the pivotal role of the 18 documented members of the LMCO (Br-lac) gene family using quantitative real-time qRT-PCR. Br-Lac-6, Br-Lac-7, Br-Lac-8, Br-Lac-16, Br-Lac-17, and Br-Lac-22 responded positively to abiotic stresses (i.e., drought, heat, and salinity). These findings set the stage for the functional diversity of the LMCO genes in B. rapa.

Genome-wide analysis of Catalase gene family reveal insights into abiotic stress response mechanism in Brassica juncea and B. rapa

Environmental stresses affect the yield and productivity of Brassica crops. Catalases are important antioxidant enzymes involved in reducing excess hydrogen peroxide produced by environmental stresses. In the present study, nine and seven CAT family members in two oilseed Brassica species (B. juncea and B. rapa) were identified with complete characterization based on gene and protein structure. Phylogenetic classification categorized CAT proteins into three classes and differentiated the monocot and dicot-specific CAT proteins. Further, the gene and protein characterizations revealed a high degree of conservation across the CAT family members. Differences were observed in the CAT-HEME binding affinity in CAT1, CAT2, and CAT3 isozymes, which could suggest their differential enzyme activities in different conditions. Furthermore, protein-protein interaction with other antioxidant proteins suggested their coordinated role in ROS scavenging mechanisms. Notably, the differential gene expression of BjuCATs and BraCATs and CAT enzyme activities suggested their crucial roles in major abiotic stresses faced by Brassica species. Promoter analysis in BjuCATs and BraCATs suggested the presence of abiotic-stress responsive cis-regulatory elements. Gene regulatory network analysis suggested miRNA and TF mediated stress response in BjuCATs and BraCATs. CAT family screening and characterization in Brassica sp. has established a basic ground for further functional validation in abiotic and heavy-metal stresses which can help in developing stress tolerant crops.

A novel CLAVATA1 mutation causes multilocularity in Brassica rapa

Locules are the seed-bearing structure of fruits. Multiple locules are associated with increased fruit size and seed set, and therefore, control of locule number is an important agronomic trait. Locule number is controlled in part by the CLAVATA-WUSCHEL pathway. Disruption of either the CLAVATA1 receptor-like kinase or its ligand CLAVATA3 can cause larger floral meristems and an increased number of locules. In an EMS mutagenized population of Brassica rapa, we identified a mutant allele that raises the number of locules from four to a range of from six to eight. Linkage mapping and genetic analysis support that the mutant phenotype is due to a missense mutation in a CLAVATA 1 (CLV1) homolog. In addition to increased locule number, additional internal gynoecia are formed in brclv1 individuals, suggesting a failure to terminate floral meristem development, which results in decreased seed production.

Optimization of Regeneration and Agrobacterium-Mediated Transformation Protocols for Bi and Multilocular Varieties of Brassica rapa

The regeneration of the high-yielding multilocular types has not been attempted, although successful regeneration and transformation in brassica have been done. Here, we report efficient regeneration and transformation protocols for two B. rapa genotypes; UAF11 and Toria. The B. rapa cv UAF11 is a multilocular, non-shattering, and high-yielding genotype, while Toria is the bilocular type. For UAF11 8 shoots and for Toria 7 shoots, explants were observed on MS supplemented with 3 mg/L BAP + 0.4 mg/L NAA + 0.01 mg/L GA₃ + 5 mg/L AgNO₃ + 0.75 mg/L Potassium Iodide (KI), MS salt supplemented with 1 mg/L IBA and 0.37 mg/L KI produced an equal number of roots (3) in UAF11 and Toria. For the establishment of transformation protocols, Agrobacterium-mediated floral dip transformation was attempted using different induction media, infection time, and flower stages. The induction medium III yielded a maximum of 7.2% transformants on half-opened flowers and 5.2% transformants on fully opened flowers in UAF11 and Toria, respectively, with 15 min of inoculation. This study would provide the basis for the improvement of tissue culture and transformation protocols in multilocular and bilocular Brassica genotypes.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Effect of Beta-Lactam Antibiotics on Microspore Embryogenesis in Brassica Species

Antibiotics are widely applied for plant cultivation in vitro to eliminate bacterial contamination. However, they can have both positive and negative effects on the cells of cultivated plants, and these effects largely depend on the type antibiotic used and its concentration. The objective of the present study was to estimate the effect of β-lactam antibiotics ampicillin (Amp) and cefotaxime (Cef) on microspore embryogenesis induction in vitro in the Brassica species. The performed experiments confirmed cefotaxime inhibits microspores in B. napus and B. oleracea, even in concentrations as low as 50 mg/L. The highest embryo yield was obtained for B. napus in the NLN-13 medium with added ampicillin in concentrations of 50–100 mg/L as an antimicrobial agent. This embryo yield was significantly higher than that obtained in a medium without supplemented antibiotics and two times higher than in the medium with added cefotaxime. Analogous results were obtained for B. oleracea and B. rapa.

Gene Expression Changes During the Allo-/Deallopolyploidization Process of Brassica napus

Gene expression changes due to allopolyploidization have been extensively studied in plants over the past few decades. Nearly all these studies focused on comparing the changes before and after genome merger. In this study, we used the uniquely restituted Brassica rapa (RBR, A_eA_e, 2n = 20) obtained from Brassica napus (A_nA_nC_nC_n, 2n = 38) to analyze the gene expression changes and its potential mechanism during the process of allo-/deallopolyploidization. RNA-seq-based transcriptome profiling identified a large number of differentially expressed genes (DEGs) between RBR and natural B. rapa (A_rA_r), suggesting potential effects of allopolyploidization/domestication of AA component of B. napus at the tetrapolyploid level. Meanwhile, it was revealed that up to 20% of gene expressions were immediately altered when compared with those in the A_n-subgenome. Interestingly, one fifth of these changes are in fact indicative of the recovery of antecedent gene expression alternations occurring since the origin of B. napus and showed association with homoeologous expression bias between A_n and C_n subgenomes. Enrichment of distinct gene ontology (GO) categories of the above sets of genes further indicated potential functional cooperation of the A_n and C_n subgenome of B. napus. Whole genome methylation analysis revealed a small number of DEGs were identified in the differentially methylated regions.

Characterization and expression profiles of miRNAs in the triploid hybrids of Brassica napus and Brassica rapa

BACKGROUND

Polyploidy provides a means of interspecific genome transfer to incorporate preferable traits from progenitor to progeny. However, few studies on miRNA expression profiles of interspecific hybrids of B. napus (AnAnCnCn) and B. rapa (ArAr) have been reported.

RESULTS

Here, we apply small RNA sequencing to explore miRNA expression patterns between B. napus, B. rapa and their F1 hybrid. Bioinformatics analysis identified 376, 378, 383 conserved miRNAs and 82, 76, 82 novel miRNAs in B. napus, B. rapa and the F1 hybrid, respectively. Moreover, 213 miRNAs were found to be differentially expressed between B. napus, B. rapa and the F1 hybrid. The present study also shows 211 miRNAs, including 77 upregulated and 134 downregulated miRNAs, to be nonadditively expressed in the F1 hybrid. Furthermore, miRNA synteny analysis revealed high genomic conservation between the genomes of B. napus, B. rapa and their F1 hybrid, with some miRNA loss and gain events in the F1 hybrid.

CONCLUSIONS

This study not only provides useful resources for exploring global miRNA expression patterns and genome structure but also facilitates genetic research on the roles of miRNAs in genomic interactions of Brassica allopolyploids.

Genome-wide identification and characterization of abiotic-stress responsive SOD (superoxide dismutase) gene family in Brassica juncea and B. rapa

Background

Abiotic stresses like drought, heat, cold and salinity cause major productivity loss in the rapeseed-mustard crops (Brassica). Major efforts have been made in the past to identify genes that provide resistance against such stresses. Superoxide dismutase (SOD) proteins, member of the metallo-enzyme family play vital role in protecting plants against abiotic stresses. In the present study, genome-wide analysis of abiotic stress responsive SOD gene family has been done in B. juncea and B. rapa.

Results

A total of 29 and 18 SOD genes were identified in B. juncea and B. rapa respectively and chromosome location mapping indicated their wide distribution across genome. On the basis of domain composition, the SODs were phylogenetically classified into sub-groups which was also substantiated by the gene structure and sub-cellular locations of SOD proteins. Functional annotation of SODs was also done by Gene Ontology (GO) mapping and the result was corroborated by the identified cis-regulatory elements in the promoter region of SOD genes. Based on FPKM analysis of SRA data available for drought, heat and salt stress, we identified 14 and 10 abiotic stress responsive SOD genes in B. rapa and B. juncea respectively. The differential expression analysis under drought and heat stress of identified abiotic-stress responsive SOD genes was done through quantitative Real Time PCR.

Conclusion

We identified abiotic-stress responsive genes that could help in improving the plant tolerance against abiotic stresses. This was the first study to describe the genome-wide analysis of SOD gene family in B. rapa and B. juncea, and the results will help in laying basic ground for future work of cloning and functional validation of SOD genes during abiotic stresses leading to Brassica crop improvement.

Electronic supplementary material

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

Genome-wide identification and expression analysis of chitinase gene family in Brassica rapa reveals its role in clubroot resistance

Chitinases, a category of pathogenesis-related proteins, are responsible for catalyzing the hydrolysis of chitin into the N-acetyl-d-glucosamine. Therefore, chitinases are believed to function as a guardian against chitin-containing pathogens. Here, we examined the role of the Brassica rapa chitinase family genes in clubroot disease. A total of 33 chitinase genes were identified and grouped into five classes based on their conserved domain. They were distributed unevenly across eight chromosomes in B. rapa, and 31 of them contained few introns (≤2). In addition, the expression of these genes was organ-specific, and 14 genes were expressed differentially in response to Plasmodiophora brassicae challenge of clubroot-susceptible (CS NIL) and resistant (CR NIL) lines. Furthermore, reduced pathogen DNA content and clubroot symptoms were observed in the CS NILs after their treatment with chitin oligosaccharides 24 h prior to inoculation with P. brassicae. The findings indicate that chitinases play a crucial role in pathogen resistance of the host plants. The results offer an insight into the role of chitinase in B. rapa-P. brassicae interaction.

Genome-Wide Gene Expression Disturbance by Single A1/C1 Chromosome Substitution in Brassica rapa Restituted From Natural B. napus

Alien chromosome substitution (CS) lines are treated as vital germplasms for breeding and genetic mapping. Previously, a whole set of nine Brassica rapa-oleracea monosonic alien addition lines (MAALs, C1-C9) was established in the background of natural B. napus genotype "Oro," after the restituted B. rapa (RBR) for Oro was realized. Herein, a monosomic substitution line with one alien C1 chromosome (Cs1) in the RBR complement was selected in the progenies of MAAL C1 and RBR, by the PCR amplification of specific gene markers and fluorescence in situ hybridization. Cs1 exhibited the whole plant morphology similar to RBR except for the defective stamens without fertile pollen grains, but it produced some seeds and progeny plants carrying the C1 chromosome at high rate besides those without the alien chromosome after pollinated by RBR. The viability of the substitution and its progeny for the RBR diploid further elucidated the functional compensation between the chromosome pairs with high homoeology. To reveal the impact of such aneuploidy on genome-wide gene expression, the transcriptomes of MAAL C1, Cs1 and euploid RBR were analyzed. Compared to RBR, Cs1 had sharply reduced gene expression level across chromosome A1, demonstrating the loss of one copy of A1 chromosome. Both additional chromosome C1 in MAAL and substitutional chromosome C1 in Cs1 caused not only cis-effect but also prevalent trans-effect differentially expressed genes. A dominant gene dosage effects prevailed among low expressed genes across chromosome A1 in Cs1, and moreover, dosage effects for some genes potentially contributed to the phenotype deviations. Our results provided novel insights into the transcriptomic perturbation and gene dosage effects on phenotype in CS related to one naturally evolved allopolyploid.

Extraction of the Constituent Subgenomes of the Natural Allopolyploid Rapeseed (Brassica napus L.)

As the dynamic nature of progenitor genomes accompanies the speciation by interspecific hybridization, the extraction of the constituent subgenome(s) from a natural allopolyploid species of long history and then restitution of the progenitor(s) provides the unique opportunity to study the genome evolution and interplay. Herein, the A subgenome from the allotetraploid oilseed rape (Brassica napus L., AACC) was extracted through inducing the preferential elimination of C-subgenome chromosomes in intertribal crosses and the progenitor B. rapa was restituted (RBR). Then by crossing and backcrossing RBR with B. napus donor, the C subgenome was in situ dissected by adding each of its nine chromosomes to the extracted A subgenome and establishing the whole set of monosonic alien addition lines (MAALs). RBR from spring-type B. napus genotype "Oro" expressed a phenotype resembling some type of B. rapa never observed before, but showed a winter-type flowering habit. This RBR had weaker growth vigor and suffered more seriously from biotic and abiotic stresses compared with Oro. The phenotypes specific for these MAALs showed the location of the related genes on the particular C-subgenome chromosomes. These MAALs exhibited obviously different frequencies in homeologous pairing and transmission of additional C-subgenome chromosomes, which were associated with the distinct degrees of their relatedness, and even with the possible genetic regulation for meiotic pairing evolved in B. napus Finally, large scaffolds undetermined for sequence assembly of B. napus were anchored to specific C-subgenome chromosomes using MAALs.

The elucidation of stress memory inheritance in Brassica rapa plants

Plants are able to maintain the memory of stress exposure throughout their ontogenesis and faithfully propagate it into the next generation. Recent evidence argues for the epigenetic nature of this phenomenon. Small RNAs (smRNAs) are one of the vital epigenetic factors because they can both affect gene expression at the place of their generation and maintain non-cell-autonomous gene regulation. Here, we have made an attempt to decipher the contribution of smRNAs to the heat-shock-induced transgenerational inheritance in Brassica rapa plants using sequencing technology. To do this, we have generated comprehensive profiles of a transcriptome and a small RNAome (smRNAome) from somatic and reproductive tissues of stressed plants and their untreated progeny. We have demonstrated that the highest tissue-specific alterations in the transcriptome and smRNAome profile are detected in tissues that were not directly exposed to stress, namely, in the endosperm and pollen. Importantly, we have revealed that the progeny of stressed plants exhibit the highest fluctuations at the smRNAome level but not at the transcriptome level. Additionally, we have uncovered the existence of heat-inducible and transgenerationally transmitted tRNA-derived small RNA fragments in plants. Finally, we suggest that miR168 and braAGO1 are involved in the stress-induced transgenerational inheritance in plants.

Molecular Characterization and Phylogeny of a Phytoplasma Associated with Phyllody Disease of toria (Brassica rapa L. subsp. dichotoma (Roxb.)) in India

Samples from toria plants (Brassica rapa L. subsp. dichotoma (Roxb.)) exhibiting phyllody, virescence, witches broom, extensive malformation of floral parts, formation of bladder like siliquae and flower sterility were collected from four different locations in India. Sequencing and phylogenetic analysis of the 16S rRNA, a part of 23S rRNA, partial sec A genes, rp gene and 16S-23S intergenic spacer region indicated that the phytoplasmas associated with toria phyllody (TP) symptoms were identical and belonged to 16SrIX phytoplasma Pigeon pea witches'-broom (PPWB) group. The iPhyClassifier generated virtual RFLP pattern of 1.25 kb 16S rDNA sequences indicated that TP phytoplasma belongs to 16SrIX-C phytoplasma subgroup. Complete 23S rRNA gene of TP phytoplasma had 2,787 nucleotides and is the first sequence of 16SrIX phytoplasma group. Restriction digestion of 16S rDNA and 23S rDNA PCR products has also shown that TP phytoplasmas from all the four locations in India were identical. Toria is a previously unreported host for a phytoplasma in16SrIX-C subgroup.