Introgressing Subgenome Components from Brassica rapa and B. carinata to B. juncea for Broadening Its Genetic Base and Exploring Intersubgenomic Heterosis

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B. juncea (AABB)

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B. carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

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.

Introgression of Heterotic Genomic Segments from Brassica carinata into Brassica juncea for Enhancing Productivity

Interspecific hybridization resulted in the creation of B. juncea introgression lines (ILs) generated from B. carinata with increased productivity and adaptability. Forty ILs were crossed with their respective B. juncea recipient parents to generate introgression line hybrids (ILHs) and the common tester (SEJ 8) was used to generate test hybrids (THs). Mid-parent heterosis in ILHs and standard heterosis in THs were calculated for eight yield and yield-related traits. Heterotic genomic regions were dissected using ten ILs with significant mid-parent heterosis in ILHs and standard heterosis in THs for seed yield. A high level of heterosis for seed yield was contributed by 1000 seed weight (13.48%) in D31_ILHs and by total siliquae/plant (14.01%) and siliqua length (10.56%) in PM30_ILHs. The heterotic ILs of DRMRIJ 31 and Pusa Mustard 30 were examined using polymorphic SNPs between the parents, and a total of 254 and 335 introgressed heterotic segments were identified, respectively. This investigation discovered potential genes, viz., PUB10, glutathione S transferase, TT4, SGT, FLA3, AP2/ERF, SANT4, MYB, and UDP-glucosyl transferase 73B3 that were previously reported to regulate yield-related traits. The heterozygosity of the FLA3 gene significantly improved siliqua length and seeds per siliqua in ILHs of Pusa Mustard 30. This research proved that interspecific hybridization is an effective means of increasing the diversity of cultivated species by introducing new genetic variants and improving the level of heterosis.

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Population genomics and subgenome evolution of the allotetraploid frog Xenopus laevis in southern Africa

Allotetraploid genomes have two distinct genomic components called subgenomes that are derived from separate diploid ancestral species. Many genomic characteristics such as gene function, expression, recombination, and transposable element mobility may differ significantly between subgenomes. To explore the possibility that subgenome population structure and gene flow may differ as well, we examined genetic variation in an allotetraploid frog-the African clawed frog (Xenopus laevis)-over the dynamic and varied habitat of its native range in southern Africa. Using reduced representation genome sequences from 91 samples from 12 localities, we found no strong evidence that population structure and gene flow differed substantially by subgenome. We then compared patterns of population structure in the nuclear genome to the mitochondrial genome using Sanger sequences from 455 samples from 183 localities. Our results provide further resolution to the geographic distribution of mitochondrial and nuclear diversity in this species and illustrate that population structure in both genomes corresponds roughly with variation in seasonal rainfall and with the topography of southern Africa.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Development of de-novo transcriptome assembly and SSRs in allohexaploid Brassica with functional annotations and identification of heat-shock proteins for thermotolerance

Crop Brassicas contain monogenomic and digenomic species, with no evidence of a trigenomic Brassica in nature. Through somatic fusion (Sinapis alba + B. juncea), a novel allohexaploid trigenomic Brassica (H1 = AABBSS; 2n = 60) was produced and used for transcriptome analysis to uncover genes for thermotolerance, annotations, and microsatellite markers for future molecular breeding. Illumina Novaseq 6000 generated a total of 76,055,546 paired-end raw reads, which were used for de-novo assembly, resulting in the development of 486,066 transcripts. A total of 133,167 coding sequences (CDSs) were predicted from transcripts with a mean length of 507.12 bp and 46.15% GC content. The BLASTX search of CDSs against public protein databases showed a maximum of 126,131 (94.72%) and a minimum of 29,810 (22.39%) positive hits. Furthermore, 953,773 gene ontology (GO) terms were found in 77,613 (58.28%) CDSs, which were divided into biological processes (49.06%), cellular components (31.67%), and molecular functions (19.27%). CDSs were assigned to 144 pathways by a pathway study using the KEGG database and 1,551 pathways by a similar analysis using the Reactome database. Further investigation led to the discovery of genes encoding over 2,000 heat shock proteins (HSPs). The discovery of a large number of HSPs in allohexaploid Brassica validated our earlier findings for heat tolerance at seed maturity. A total of 15,736 SSRs have been found in 13,595 CDSs, with an average of one SSR per 4.29 kb length and an SSR frequency of 11.82%. The first transcriptome assembly of a meiotically stable allohexaploid Brassica has been given in this article, along with functional annotations and the presence of SSRs, which could aid future genetic and genomic studies.

Characterization of APX and APX-R gene family in Brassica juncea and B. rapa for tolerance against abiotic stresses

APX and APX-R gene families were identified and characterized in two important oilseed species of Brassica. Gene expression under abiotic stress conditions, recombinant protein expression, and analysis further divulged their drought, heat, and salt-responsive behavior. Ascorbate peroxidases (APX) are heme-dependent enzymes that rid the cells of H₂O₂ and regulate diverse biological processes. In the present study, we performed APX gene family characterization in two Brassica sp. (B. juncea and B. rapa) as these are commercially important oilseed crops and affected severely by abiotic stresses. We identified 16 BjuAPX and 9 BraAPX genes and 2 APX-R genes each in B. juncea and B. rapa genomes, respectively. Phylogenetic analysis divided the APX genes into five distinct clades, which exhibited conservation in the gene structure, motif organization, and sub-cellular location within the clade. Structural analysis of APX and APX-R proteins revealed the amino acid substitutions in conserved domains of APX-R proteins. The expression profiling of BjuAPX and BraAPX genes showed that 3 BjuAPX, 7BraAPX, and 2 BraAPX-R genes were drought and heat responsive. Notably, BjuAAPX1a, BjuAPX1d, BjuAAPX6, BraAAPX1a, BraAAPX2, and BraAAPX3b showed high expression levels in RT-qPCR. Cis-regulatory elements in APX and APX-R gene promoters supported the differential behavior of these genes. Further, two stress-responsive genes BjuAPX1d and BraAAPX2 were cloned, characterized, and their roles were validated under heat, drought, salt, and cold stress in bacterial expression system. This study for the first time reports the presence of APX activity in dimeric and LMW form of purified BraAAPX2 protein. The study may help pave way for developing abiotic stress-tolerant Brassica crops.

Deployment of Brassica carinata A. Braun Derived Brassica juncea (L.) Czern. Lines for Improving Heterosis and Water Use Efficiency Under Water Deficit Stress Conditions

Among Brassica species, Ethiopian mustard (Brassica carinata A. Braun) is known to tolerate most abiotic stresses, including drought. Drought caused by low and erratic rainfall in semi-arid regions consistently challenges rapeseed mustard productivity. Development of B. carinata-derived lines (CDLs) in Brassica juncea (L.) Czern. nuclear background, carrying genomic segments from B. carinata, are expected to tolerate moisture deficit stress conditions. The present study was, thus, aimed to establish the phenomenon "heterosis" for drought tolerance and water use efficiency by evaluating 105 hybrids developed from intermating 15 CDLs in half diallel fashion. Data on 17 seed yield and yield contributing traits were recorded under two different environments, viz., irrigated and rainfed conditions. Traits under study were found to be governed by both additive and non-additive types of gene action. Average degree of dominance was higher (>2) for yield and yield contributing traits, viz., secondary branches/plant, point to first siliqua on main shoot, total siliquae/plant, 1,000-seed weight, seed yield/plant, biological yield, harvest index, and seed yield/hectare under rainfed conditions, clearly indicating that higher productivity under drought conditions can be realised through the development of hybrids. Out of 15, highly significant general combining ability (GCA) effects for seven CDLs were observed under rainfed condition. Furthermore, nine and six hybrids expressed highly significant specific combining ability (SCA) effects and > 50% heterobeltiosis for yield contributing traits under rainfed and irrigated conditions, respectively. Water use efficiency (WUE) of parental CDLs and hybrids varied from 2.05 to 2.57 kg m-3 under rainfed, while 1.10 to 1.28 kg m-3 under irrigated conditions. Hybrids expressed higher WUE than parental lines under both water regimes. Furthermore, selection indices such as drought tolerance index (DTI) and mean relative performance (MRP) were identified to be efficient in the selection of productive CDLs and hybrids under drought conditions. Nine hybrids, identified as highly productive in the present study, can further be exploited for improving the yield of Indian mustard in drought-prone areas. Usefulness of interspecific hybridisation in the development of B. carinata-derived B. juncea lines for improving heterosis and WUE is, thus, well demonstrated through the present study.

Development of B. carinata with super-high erucic acid content through interspecific hybridization

KEY MESSAGE

Disomic alien chromosome addition Brassica carinata lines with super-high erucic acid content were developed through interspecific hybridization with B. juncea and characterized using molecular, cytological and biochemical techniques. Brassica carinata [A.] Braun (BBCC, 2n = 34) is a climate-resilient oilseed. Its seed oil is high in erucic acid (> 40%), rendering it well suited for the production of biofuel and other bio-based applications. To enhance the competitiveness of B. carinata with high erucic B. napus (HEAR), lines with super-high erucic acid content were developed through interspecific hybridization. To this end, a fad2B null allele from Brassica juncea (AABB, 2n = 36) was introgressed into B. carinata, resulting in a B. carinata fad2B mutant with erucic acid levels of over 50%. Subsequently, the FAE allele from B. rapa spp. yellow sarson (AA, 2n = 20) was transferred to the fad2B B. carinata line, yielding lines with erucic acid contents of up to 57.9%. Molecular analysis using the Brassica 90 K Illumina Infinium™ SNP genotyping array identified these lines as disomic alien chromosome addition lines, with two extra A08 chromosomes containing the BrFAE gene. The alien chromosomes from B. rapa were clearly distinguished by molecular cytogenetics in one of the addition lines. Analysis of microspore-derived offspring and hybrids from crosses with a CMS B. carinata line showed that the transfer rate of the A08 chromosome into male gametes was over 98%, resulting in almost completely stable transmission of an A08 chromosome copy into the progeny. The increase in erucic acid levels was accompanied by changes in the proportions of other fatty acids depending on the genetic changes that were introduced in the interspecific hybrids, providing valuable insights into erucic acid metabolism in Brassica.

Current Status of the Disease-Resistant Gene(s)/QTLs, and Strategies for Improvement in Brassica juncea

Brassica juncea is a major oilseed crop in tropical and subtropical countries, especially in south-east Asia like India, China, Bangladesh, and Pakistan. The widespread cultivation of genetically similar varieties tends to attract fungal pathogens which cause heavy yield losses in the absence of resistant sources. The conventional disease management techniques are often expensive, have limited efficacy, and cause additional harm to the environment. A substantial approach is to identify and use of resistance sources within the Brassica hosts and other non-hosts to ensure sustainable oilseed crop production. In the present review, we discuss six major fungal pathogens of B. juncea: Sclerotinia stem rot (Sclerotinia sclerotiorum), Alternaria blight (Alternaria brassicae), White rust (Albugo candida), Downy mildew (Hyaloperonospora parasitica), Powdery mildew (Erysiphe cruciferarum), and Blackleg (Leptoshaeria maculans). From discussing studies on pathogen prevalence in B. juncea, the review then focuses on highlighting the resistance sources and quantitative trait loci/gene identified so far from Brassicaceae and non-filial sources against these fungal pathogens. The problems in the identification of resistance sources for B. juncea concerning genome complexity in host subpopulation and pathotypes were addressed. Emphasis has been laid on more elaborate and coordinated research to identify and deploy R genes, robust techniques, and research materials. Examples of fully characterized genes conferring resistance have been discussed that can be transformed into B. juncea using advanced genomics tools. Lastly, effective strategies for B. juncea improvement through introgression of novel R genes, development of pre-breeding resistant lines, characterization of pathotypes, and defense-related secondary metabolites have been provided suggesting the plan for the development of resistant B. juncea.

Development of a Yellow-Seeded Stable Allohexaploid Brassica Through Inter-Generic Somatic Hybridization With a High Degree of Fertility and Resistance to Sclerotinia sclerotiorum

The Brassica coenospeceis have treasure troves of genes that could be beneficial if introgressed into cultivated Brassicas to combat the current conditions of climate change. Introducing genetic variability through plant speciation with polyploidization is well documented, where ploidy augmentation of inter-generic allohexaploids using somatic hybridization has significantly contributed to genetic base broadening. Sinapis alba is a member of the Brassicaceae family that possesses valuable genes, including genes conferring resistance to Sclerotinia sclerotiorum, Alternaria brassicae, pod shattering, heat, and drought stress. This work aimed to synthesize stable allohexaploid (AABBSS) Brassica while incorporating the yellow-seed trait and resistance to S. sclerotiorum stem rot. The two fertile and stable allohexaploids were developed by polyethylene glycol mediated protoplast fusions between Brassica juncea (AABB) and S. alba (SS) and named as JS1 and JS2. These symmetric hybrids (2n = 60) were validated using morphological and molecular cytology techniques and were found to be stable over consecutive generations. The complete chromosome constitution of the three genomes was determined through genomic in situ hybridization of mitotic cells probed with S. alba genomic DNA labeled with fluorescein isothiocyanate. These two allohexaploids showed 24 hybridization signals demonstrating the presence of complete diploid chromosomes from S. alba and 36 chromosomes from B. juncea. The meiotic pollen mother cell showed 30 bivalent sets of all the 60 chromosomes and none of univalent or trivalent observed during meiosis. Moreover, the backcross progeny 1 plant revealed 12 hybridization signals out of a total of 48 chromosome counts. Proper pairing and separation were recorded at the meiotic metaphase and anaphase, which proved the stability of the allohexaploid and their backcross progeny. When screening, the allohexaploid (JS2) of B. juncea and S. alba displayed a high degree of resistance to S. sclerotiorum rot along with a half-yellow and half-brown (mosaic) seed coat color, while the B. juncea and S. alba allohexaplopid1 (JS1) displayed a yellow seed coat color with the same degree of resistance to Sclerotinia rot.

Potential of the C Genome of the Different Variants of Brassica oleracea for Heterosis in Spring B. napus Canola

The genetic base of Brassica napus canola need to be broadened for exploitation of heterosis at a greater level in the breeding of F1 hybrid canola cultivars. In this study, we evaluated 228 inbred B. napus canola lines derived from six B. napus × B. oleracea interspecific crosses and following two breeding methods (F2- and BC1-derived lines) to understand the effect of the B. oleracea alleles on heterosis for different agronomic and seed quality traits. Test hybrids of the inbreds derived from crosses involving vars. botrytis (cauliflower), alboglabra (Chinese kale) and capitata (cabbage) cv. Badger Shipper, on an average, gave about 10% mid-parent heterosis (MPH), and about 67% of the test hybrids gave higher seed yield than the common B. napus parent indicating that B. oleracea alleles can contribute to heterosis for seed yield in spring B. napus canola hybrids. This was also evident from a positive correlation of the genetic distance of the inbred lines from the common B. napus parent with MPH for seed yield (r = 0.31) as well as with hybrid yield (r = 0.26). Almost no correlation was found between genetic distance and MPH for seed oil and protein content as well as with the performance of the test hybrids for these two traits. The occurrence of positive correlation between seed yield of the inbred lines and test hybrids suggested the importance of the genes exerting additive effect for high seed yield in the hybrids. Very little or almost no heterosis was found for the other agronomic traits as well as for seed oil and protein content. While comparing the two breeding methods, no significant difference was found for seed yield of the test hybrids or the level of MPH; however, the BC1-derived inbred and test hybrid populations flowered and matured earlier and had longer grain-filling period than the F2-derived population. Thus, the results suggested that the B. oleracea gene pool can be used in the breeding of spring B. napus canola to improve seed yield in hybrid cultivars.

Reconstituting the genome of a young allopolyploid crop, Brassica napus, with its related species

Brassica napus (Aⁿ Aⁿ Cⁿ Cⁿ ) is an important worldwide oilseed crop, but it is a young allotetraploid with a short evolutionary history and limited genetic diversity. To significantly broaden its genetic diversity and create a novel heterotic population for sustainable rapeseed breeding, this study reconstituted the genome of B. napus by replacing it with the subgenomes from 122 accessions of Brassica rapa (A^r A^r ) and 74 accessions of Brassica carinata (B^c B^c C^c C^c ) and developing a novel gene pool of B. napus through five rounds of extensive recurrent selection. When compared with traditional B. napus using SSR markers and high-throughput SNP/Indel markers through genotyping by sequencing, the newly developed gene pool and its homozygous progenies exhibited a large genetic distance, rich allelic diversity, new alleles and exotic allelic introgression across all 19 AC chromosomes. In addition to the abundant genomic variation detected in the AC genome, we also detected considerable introgression from the eight chromosomes of the B genome. Extensive trait variation and some genetic improvements were present from the early recurrent selection to later generations. This novel gene pool produced equally rich phenotypic variation and should be valuable for rapeseed genetic improvement. By reconstituting the genome of B. napus by introducing subgenomic variation within and between the related species using intense selection and recombination, the whole genome could be substantially reorganized. These results serve as an example of the manipulation of the genome of a young allopolyploid and provide insights into its rapid genome evolution affected by interspecific and intraspecific crosses.

Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea

Heterosis refers to the superior performance of F1 hybrids over their respective parental inbred lines. Although the genetic and expression basis of heterosis have been previously investigated, the metabolic basis for this phenomenon is poorly understood. In a preliminary morphological study in Brassica juncea, we observed significant heterosis at the 50% flowering stage, wherein both the growth and reproduction of F1 reciprocal hybrids were greater than that of their parents. To identify the possible metabolic causes or consequences of this heterosis, we carried out targeted LC-MS analysis of 48 primary (amino acids and sugars) and secondary metabolites (phytohormones, glucosinolates, flavonoids, and phenolic esters) in five developmental tissues at 50% flowering in hybrids and inbred parents. Principal component analysis (PCA) of metabolites clearly separated inbred lines from their hybrids, particularly in the bud tissues. In general, secondary metabolites displayed more negative heterosis values in comparison to primary metabolites. The tested primary and secondary metabolites displayed both additive and non-additive modes of inheritance in F1 hybrids, wherein the number of metabolites showing an additive mode of inheritance were higher in buds and siliques (52.77-97.14%) compared to leaf tissues (47.37-80%). Partial least regression (PLS) analysis further showed that primary metabolites, in general, displayed higher association with morphological parameters in F1 hybrids. Overall, our results are consistent with a resource-cost model for heterosis in B. juncea, where metabolite allocation in hybrids appears to favor growth, at the expense of secondary metabolism.

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.

Investigation of the Genetic Diversity and Quantitative Trait Loci Accounting for Important Agronomic and Seed Quality Traits in Brassica carinata

Brassica carinata (BBCC) is an allotetraploid in Brassicas with unique alleles for agronomic traits and has huge potential as source for biodiesel production. To investigate the genome-wide molecular diversity, population structure and linkage disequilibrium (LD) pattern in this species, we genotyped a panel of 81 accessions of B. carinata with genotyping by sequencing approach DArTseq, generating a total of 54,510 polymorphic markers. Two subpopulations were exhibited in the B. carinata accessions. The average distance of LD decay (r2 = 0.1) in B subgenome (0.25 Mb) was shorter than that of C subgenome (0.40 Mb). Genome-wide association analysis (GWAS) identified a total of seven markers significantly associated with five seed quality traits in two experiments. To further identify the quantitative trait loci (QTL) for important agronomic and seed quality traits, we phenotyped a doubled haploid (DH) mapping population derived from the "YW" cross between two parents (Y-BcDH64 and W-BcDH76) representing from the two subpopulations. The YW DH population and its parents were grown in three contrasting environments; spring (Hezheng and Xining, China), semi-winter (Wuhan, China), and spring (Wagga Wagga, Australia) across 5 years for QTL mapping. Genetic bases of phenotypic variation in seed yield and its seven related traits, and six seed quality traits were determined. A total of 282 consensus QTL accounting for these traits were identified including nine major QTL for flowering time, oleic acid, linolenic acid, pod number of main inflorescence, and seed weight. Of these, 109 and 134 QTL were specific to spring and semi-winter environment, respectively, while 39 consensus QTL were identified in both contrasting environments. Two QTL identified for linolenic acid (B3) and erucic acid (C7) were validated in the diverse lines used for GWAS. A total of 25 QTL accounting for flowering time, erucic acid, and oleic acid were aligned to the homologous QTL or candidate gene regions in the C genome of B. napus. These results would not only provide insights for genetic improvement of this species, but will also identify useful genetic variation hidden in the Cc subgenome of B. carinata to improve canola cultivars.

Molecular Diversity Analysis and Genetic Mapping of Pod Shatter Resistance Loci in Brassica carinata L

Seed lost due to easy pod dehiscence at maturity (pod shatter) is a major problem in several members of Brassicaceae family. We investigated the level of pod shatter resistance in Ethiopian mustard (Brassica carinata) and identified quantitative trait loci (QTL) for targeted introgression of this trait in Ethiopian mustard and its close relatives of the genus Brassica. A set of 83 accessions of B. carinata, collected from the Australian Grains Genebank, was evaluated for pod shatter resistance based on pod rupture energy (RE). In comparison to B. napus (RE = 2.16 mJ), B. carinata accessions had higher RE values (2.53 to 20.82 mJ). A genetic linkage map of an F₂ population from two contrasting B. carinata selections, BC73526 (shatter resistant with high RE) and BC73524 (shatter prone with low RE) comprising 300 individuals, was constructed using a set of 6,464 high quality DArTseq markers and subsequently used for QTL analysis. Genetic analysis of the F₂ and F_2:3 derived lines revealed five statistically significant QTL (LOD ≥ 3) that are linked with pod shatter resistance on chromosomes B1, B3, B8, and C5. Herein, we report for the first time, identification of genetic loci associated with pod shatter resistance in B. carinata. These characterized accessions would be useful in Brassica breeding programs for introgression of pod shatter resistance alleles in to elite breeding lines. Molecular markers would assist marker-assisted selection for tracing the introgression of resistant alleles. Our results suggest that the value of the germplasm collections can be harnessed through genetic and genomics tools.