Evidence from genome-wide simple sequence repeat markers for a polyphyletic origin and secondary centers of genetic diversity of Brassica juncea in China and India

Transcriptome Dynamics of Brassica juncea Leaves in Response to Omnivorous Beet Armyworm (Spodoptera exigua, Hübner)

Cruciferous plants manufacture glucosinolates (GSLs) as special and important defense compounds against insects. However, how insect feeding induces glucosinolates in Brassica to mediate insect resistance, and how plants regulate the strength of anti-insect defense response during insect feeding, remains unclear. Here, mustard (Brassica juncea), a widely cultivated Brassica plant, and beet armyworm (Spodoptera exigua), an economically important polyphagous pest of many crops, were used to analyze the changes in GSLs and transcriptome of Brassica during insect feeding, thereby revealing the plant-insect interaction in Brassica plants. The results showed that the content of GSLs began to significantly increase after 48 h of herbivory by S. exigua, with sinigrin as the main component. Transcriptome analysis showed that a total of 8940 DEGs were identified in mustard challenged with beet armyworm larvae. The functional enrichment results revealed that the pathways related to the biosynthesis of glucosinolate and jasmonic acid were significantly enriched by upregulated DEGs, suggesting that mustard might provide a defense against herbivory by inducing JA biosynthesis and then promoting GSL accumulation. Surprisingly, genes regulating JA catabolism and inactivation were also activated, and both JA signaling repressors (JAZs and JAMs) and activators (MYCs and NACs) were upregulated during herbivory. Taken together, our results indicate that the accumulation of GSLs regulated by JA signaling, and the regulation of active and inactive JA compound conversion, as well as the activation of JA signaling repressors and activators, collectively control the anti-insect defense response and avoid over-stunted growth in mustard during insect feeding.

Exchanges of economic plants along the land silk road

BACKGROUD

The greatest contribution of the Silk Road is to communicate among different countries and nationalities, and promote two-way cultural exchanges between the East and the West. We now have clearer understanding about how material civilization and religious culture of Central Asia and West Asia spread eastward along the Land Silk Road. However, there is controversial about how crops migrate along the Land Silk Road.

RESULTS

We summarize archaeology, genetics, and genomics data to explore crop migration patterns. Of the 207 crops that were domesticated along the Land Silk Road, 19 for which genomic evidence was available were selected for discussion.

CONCLUSIONS

There were conflicting lines of evidence for the domestication of Tibetan barley, mustard, lettuce, buckwheat, and chickpea. The main reasons for the conflicting results may include incomplete early knowledge, record differences in different period, sample sizes, and data analysis techniques. There was strong evidence that Tibetan barley, barley, wheat, and jujube were introduced into China before the existence of the Land Silk Road; and mustard, lettuce, buckwheat, chickpea, alfalfa, walnut, cauliflower, grape, spinach, apple, cucumber, mulberry, and pea spread to China via trade and human migration along the Land Silk Road.

First Report of Root-knot Nematode, Meloidogyne Graminicola on Brassica Juncea in China

In southern China, the staple food rice (Oryza sativa) field is commonly rotated with brown mustard Brassica juncea. Root-knot nematodes (RKNs) are a major threat to rice production. From 2019 to 2021, B. juncea in 56 fields from 26 counties in Guangxi Province were observed with symptoms of leaf yellowing, stunting, and several hook-shaped galls on the roots. Females and egg masses of Meloidogyne sp. were found within the galls. The females, males, and second-stage juveniles (J2s) were collected, and identified with morphological and molecular characteristics and female perineal patterns. Brassica juncea was transplanted in pots and a pathogenicity test was conducted to confirm the species as Meloidogyne graminicola. In China, this is the first record of a natural infection of mustard with M. graminicola, and this finding has great importance for Chinese mustard production, since this nematode may damage mustard plants and become an additional problem for this crop.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

Expanding the genetic variation of Brassica juncea by introgression of the Brassica rapa Genome

Brassica juncea is an important vegetable and oil crop cultivated worldwide. To increase its genetic variation, we introgressed the A-genome of Brassica rapa into B. juncea. We used three each of heading and semi-heading B. juncea accessions as recipient parents and a B. rapa line, B9008, as the donor parent. We obtained 101 BC1S1 lines in total with expanded phenotypic variations such as leafy head shapes. We developed 132 SNP markers, which could distinguish the A genome of B. juncea from B. rapa genome, and tracked the introgression of B. rapa segments in the new B. juncea germplasm. On average, 59.2% of the B. juncea A genome in the B. juncea introgression lines was covered by the donor segments. We also identified three markers whose donor genotype frequencies were significantly lower than the theoretical value, suggesting strong selection of the recipient genotype during the introgression process. We provide an effective strategy to evaluate the diversity of the new germplasm based on the combination of parental re-sequencing data and marker genotyping results. Further genetic analysis of 1642 SNPs showed that the genetic diversity of the new B. juncea germplasm with the introgressed B. rapa genome was significantly increased. This study illustrated the potential for expanding the genetic diversity of B. juncea through the introgression of the B. rapa genome.

Genomic signatures of vegetable and oilseed allopolyploid Brassica juncea and genetic loci controlling the accumulation of glucosinolates

Allopolyploid Brassica juncea crops in Brassicaceae are becoming increasingly revitalized as vegetables and oilseeds owing to wide adaptability and significant economic values. However, the genomic differentiation of diversified vegetables and oilseed B. juncea and the genetic basis underlying glucosinolates accumulation have yet to be elucidated. To address this knowledge gap, we report the sequencing of pairwise genomes of vegetable and oilseed B. juncea at chromosome scale. Comparative genomics analysis unveils panoramic structural variation footprints, particularly the genetic loci of HSP20 and TGA1 associated with abiotic and biotic stresses responses between oilseed and vegetable subgroups. We anchored two major loci of MYB28 (HAG1) orthologues caused by copy number variations on A02 and A09 chromosomes using scored genomic SNPs-based GWAS that are responsible for seed oil quality-determining glucosinolates biosynthesis. These findings will provide valuable repertories of polyploidy genomic information enabling polyploidy genome evolution studies and precise genomic selections for crucial traits like functional components of glucosinolates in B. juncea crops and beyond.

Data on agronomic traits, biochemical composition of lipids, proteins and polysaccharides and rheological measurement in a brown mustard seed collection

The data were collected from a brown mustard seeds collection of 18 accessions during two years and in three distinct sites of production in France. The 18 accessions of mustard seeds were selected to be representative of genetic, agronomical and technological variabilities. All accessions were produced in the "Bourgogne" area. This article describes agronomical data (PMG, yield), genotyping data, global composition of mustard seeds (lipids, proteins and polysaccharides) and fine composition of the previous macronutrients potentially involved in the technological properties (fatty acids, storage proteins and osidic composition of polysaccharides). Additional data regarding the potential rheological property of each accessions were also reported. These data can be reused by food industries, breeders and geneticists in order to understand pedoclimatic effects (year and location) and the relation between mustard seed composition and the end-uses properties (paste mustard quality).

The evolution of the expansin gene family in Brassica species

Expansin gene (EXP) family plays important roles in plant growth and crop improvement. However, it has not been well studied in the Brassica genus that includes several important agricultural and horticultural crops. To get insight to the evolution and expansion of EXP family in Brassica, Brassica EXPs which are homologues of 35 known AtEXPs of Arabidopsis were comprehensively and systematically analyzed in the present study. In total, 340 Brassica EXPs were clustered into four groups that corresponded multiple alignment to four subfamilies of AtEXPs, with divergent conserved motifs and cis-acting elements among groups. To understand the expansion of EXP family, an integrated genomic block system was constructed among Arabidopsis and Brassica species based on 24 known ancestral karyotype blocks. Obvious gene loss, segmental duplication, tandem duplication and DNA sequence repeat events were found during the expansion of Brassica EXPs, of which the segmental duplication was possibly the major driving force. The divergence time was estimated in 1109 orthologs pairs of EXPs, revealing the divergence of Brassica EXPs from AtEXPs during ~30 MYA, and the divergence of EXPs among Brassica species during 13.50-17.94 MYA. Selective mode analysis revealed that the purifying selection was the major contributor to expansion of Brassica EXPs. This study provides new insights into the evolution and expansion of the EXP family in Brassica genus.

Genomic insights into the origin, domestication and diversification of Brassica juncea

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000-14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

Investigating genetic relationship of Brassica juncea with B. nigra via virtual allopolyploidy and hexaploidy strategy

Brassica juncea is an important economic crop of the world; however, the narrow genetic base of this crop has tremendously decreased its crop productivity. As an ancestral species of B. juncea, B. nigra is of great importance in widening the genetic diversity of B. juncea. In the present study, 42 SSR markers were employed to screen the genetic diversity among 83 B. nigra, 16 B. juncea, and other Brassica accessions. The molecular characteristics of 498 virtual B. juncea lines were deduced based on the bands of B. nigra and B. rapa via a virtual allopolyploid strategy, and then compared with natural B. juncea accessions. It was found that B. nigra had rich genetic diversity and could be classified into four subgroups, of which subgroup B-III and subgroup B-IV exhibited the closest and the most distant genetic relationship with B. juncea, respectively. To verify this, a hexaploidy strategy was applied to generated synthetic B. juncea from 20 B. nigra accessions, resulting in 45 new-type B. juncea genotypes. The genetic analyses detected that synthetic B. juncea derived from B. nigra in subgroup B-III was close to natural B. juncea, while B. juncea synthesized with B. nigra from subgroup B-IV exhibited wide genetic diversity and was most distant with current B. juncea. This study revealed a great potential of B. nigra in widening genetic diversity of B. juncea particularly using B. nigra in subgroup B-IV, and is helpful in better understanding of the genetic relationship between B. nigra and B. juncea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-020-01197-7.

Validation of an Associative Transcriptomics platform in the polyploid crop species Brassica juncea by dissection of the genetic architecture of agronomic and quality traits

The development of more productive crops will be key to addressing the challenges that climate change, population growth and diminishing resources pose to global food security. Advanced 'omics techniques can help to accelerate breeding by facilitating the identification of genetic markers for use in marker-assisted selection. Here, we present the validation of a new Associative Transcriptomics platform in the important oilseed crop Brassica juncea. To develop this platform, we established a pan-transcriptome reference for B. juncea, to which we mapped transcriptome data from a diverse panel of B. juncea accessions. From this panel, we identified 355 050 single nucleotide polymorphism variants and quantified the abundance of 93 963 transcripts. Subsequent association analysis of functional genotypes against a number of important agronomic and quality traits revealed a promising candidate gene for seed weight, BjA.TTL, as well as additional markers linked to seed colour and vitamin E content. The establishment of the first full-scale Associative Transcriptomics platform for B. juncea enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Genetic and Comparative Transcriptome Analysis Revealed DEGs Involved in the Purple Leaf Formation in Brassica juncea

Brassica juncea is an important dietary vegetable cultivated and consumed in China for its edible stalks and leaves. The purple leaf mustard, which is rich in anthocyanins, is eye-catching and delivers valuable nutrition. However, the molecular mechanism involved in anthocyanin biosynthesis has not been well studied in B. juncea. Here, histological and transcriptome analyses were used to characterize the purple leaf color and gene expression profiles. Free-hand section analysis showed that the anthocyanin was mainly accumulated in the adaxial epidermal leaf cells. The anthocyanin content in the purple leaves was significantly higher than that in the green leaves. To investigate the critical genes and pathways involved in anthocyanin biosynthesis and accumulation, the transcriptome analysis was used to identify the differentially expressed genes (DEGs) between the purple and green leaves from the backcrossed BC3 segregation population in B. juncea. A total of 2,286 different expressed genes were identified between the purple and green leaves. Among them, 1,593 DEGs were up-regulated and 693 DEGs were down-regulated. There were 213 differently expressed transcription factors among them. The MYB and bHLH transcription factors, which may regulate anthocyanin biosynthesis, were up-regulated in the purple leaves. Interestingly, most of the genes involved in plant-pathogen interaction pathway were also up-regulated in the purple leaves. The late biosynthetic genes involved in anthocyanin biosynthesis were highly up-regulated in the purple leaves of B. juncea. The up regulation of BjTT8 and BjMYC2 and anthocyanin biosynthetic genes (BjC4H, BjDFR, and BjANS) may activate the purple leaf formation in B. juncea. This study may help to understand the transcriptional regulation of anthocyanin biosynthesis in B. juncea.

A Mapped Locus on LG A6 of Brassica juncea Line Tumida Conferring Resistance to White Rust Contains a CNL Type R Gene

White rust, causal agent oomycete Albugo candida, is a significant disease of the cultivated Brassica species. The Indian gene pool lines of oilseed mustard, Brassica juncea, are highly susceptible to the pathogen. Resistance to A. candida has been reported in the east European gene pool lines of mustard and mapped to LG A4 in line Heera and LG A5 in line Donskaja-IV. A new resistance-conferring locus to A. candida isolate AcB1 has been mapped to LG A6 of B. juncea line Tumida-a Chinese vegetable type mustard using an F1DH mapping population that has been developed from a Tumida × Varuna (susceptible Indian gene pool line) cross. A molecular map containing 8,303 genic and GBS markers was used to map the resistance trait to an interval of 63.0 cM-70.8 cM on LG A6. Genome assemblies of Tumida and Varuna were used to find the genes present within the flanking markers discerned by genetic mapping. The most likely candidate gene in the mapped interval is BjuA046215, a CC-NBS-LRR (CNL) type R gene that encodes a protein with all the specific subdomains of the proteins encoded by such genes. Alleles of BjuA046215 in Varuna and other lines of the Indian and the east European gene pools encode proteins that have truncated LRR domains. Analysis of the syntenic regions in some of the Brassicaceae genomes and phylogenetic analysis of CNL type R genes showed BjuA046215 to be closely related to a recently described white rust resistance-conferring R gene BjuWRR1 in B. juncea Donskaja-IV, both belonging to the CNL-D group of R genes. Related R genes in Arabidopsis thaliana confer resistance to another oomycete, Peronospora parasitica.

Maternal Inheritance of U's Triangle and Evolutionary Process of Brassica Mitochondrial Genomes

The sequences and genomic structures of plant mitochondrial (mt) genomes provide unique material for phylogenetic studies. The nature of uniparental inheritance renders an advantage when utilizing mt genomes for determining the parental sources of hybridized taxa. In this study, a concatenated matrix of mt genes was used to infer the phylogenetic relationships of six cultivated Brassica taxa and explore the maternal origins of three allotetraploids. The well-resolved sister relationships between two pairs of diploid and allotetraploid taxa suggest that Brassica carinata (car) possessed a maternal origin from Brassica nigra, while Brassica juncea (jun) was maternally derived from Brassica rapa (cam). Another allotetraploid taxon, Brassica napus (cv. Wester) may have been maternally derived from the common ancestor of B. rapa and Brassica oleracea (ole), and/or have undergone (an) extra hybridization event(s) along its evolutionary history. The characteristics of Brassica mt genomic structures also supported the phylogenetic results. Sinapis arvensis was nested inside the Brassica species, sister to the B. nigra-B. carinata lineage, and possessed an mt genome structure that mostly resembled B. nigra. Collectively, the evidence supported a systematic revision that placed S. arvensis within Brassica. Finally, ancestral mt genomes at each evolutionary node of Brassica were reconstructed, and the detailed and dynamic evolution of Brassica mt genomes was successfully reproduced. The mt genome of B. nigra structurally resembled that of the Brassica ancestor the most, with only one reversion of a block, and the Brassica oleracea underwent the most drastic changes. These findings suggested that repeat-mediated recombinations were largely responsible for the observed structural variations in the evolutionary history of Brassica mt genomes.

A multiplex PCR for rapid identification of Brassica species in the triangle of U

BACKGROUND

Within the Brassicaceae, six species from the genus Brassica are widely cultivated throughout the world as oilseed, condiment, fodder or vegetable crops. The genetic relationships among the six Brassica species are described by U's triangle model. Extensive shared traits and diverse morphotypes among Brassica species make identification and classification based on phenotypic data alone challenging and unreliable, especially when dealing with large germplasm collections. Consequently, a major issue for genebank collections is ensuring the correct identification of species. Molecular genotyping based on simple sequence repeat (SSR) marker sequencing or the Illumina Infinium Brassica napus 60K single nucleotide polymorphism (SNP) array has been used to identify species and assess genetic diversity of Brassica collections. However, these methods are technically challenging, expensive and time-consuming, making them unsuitable for routine or rapid screening of Brassica accessions for germplasm management. A cheaper, faster and simpler method for Brassica species identification is described here.

RESULTS

A multiplex polymerase chain reaction (MPCR) consisting of new and existing primers specific to the Brassica A, B and C genomes was able to reliably distinguish all six Brassica species in the triangle of U with 16 control samples of known species identity. Further validation against 120 Brassica accessions previously genotyped showed that the MPCR is highly accurate and comparable to more advanced techniques such as SSR marker sequencing or the Illumina Infinium B. napus 60K SNP array. In addition, the MPCR was sensitive enough to detect seed contaminations in pooled seed samples of Brassica accessions.

CONCLUSION

A cheap and fast multiplex PCR assay for identification of Brassica species in the triangle of U was developed and validated in this study. The MPCR assay can be readily implemented in any basic molecular laboratory and should prove useful for the management of Brassica germplasm collections in genebanks.

Polyploidy and the relationship between leaf structure and function: implications for correlated evolution of anatomy, morphology, and physiology in Brassica

BACKGROUND

Polyploidy is well studied from a genetic and genomic perspective, but the morphological, anatomical, and physiological consequences of polyploidy remain relatively uncharacterized. Whether these potential changes bear on functional integration or are idiosyncratic remains an open question. Repeated allotetraploid events and multiple genomic combinations as well as overlapping targets of artificial selection make the Brassica triangle an excellent system for exploring variation in the connection between plant structure (anatomy and morphology) and function (physiology). We examine phenotypic integration among structural aspects of leaves including external morphology and internal anatomy with leaf-level physiology among several species of Brassica. We compare diploid and allotetraploid species to ascertain patterns of phenotypic correlations among structural and functional traits and test the hypothesis that allotetraploidy results in trait disintegration allowing for transgressive phenotypes and additional evolutionary and crop improvement potential.

RESULTS

Among six Brassica species, we found significant effects of species and ploidy level for morphological, anatomical and physiological traits. We identified three suites of intercorrelated traits in both diploid parents and allotetraploids: Morphological traits (such as leaf area and perimeter) anatomic traits (including ab- and ad- axial epidermis) and aspects of physiology. In general, there were more correlations between structural and functional traits for allotetraploid hybrids than diploid parents. Parents and hybrids did not have any significant structure-function correlations in common. Of particular note, there were no significant correlations between morphological structure and physiological function in the diploid parents. Increased phenotypic integration in the allotetraploid hybrids may be due, in part, to increased trait ranges or simply different structure-function relationships.

CONCLUSIONS

Genomic and chromosomal instability in early generation allotetraploids may allow Brassica species to explore new trait space and potentially reach higher adaptive peaks than their progenitor species could, despite temporary fitness costs associated with unstable genomes. The trait correlations that disappear after hybridization as well as the novel trait correlations observed in allotetraploid hybrids may represent relatively evolutionarily labile associations and therefore could be ideal targets for artificial selection and crop improvement.

The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection

The Brassica genus encompasses three diploid and three allopolyploid genomes, but a clear understanding of the evolution of agriculturally important traits via polyploidy is lacking. We assembled an allopolyploid Brassica juncea genome by shotgun and single-molecule reads integrated to genomic and genetic maps. We discovered that the A subgenomes of B. juncea and Brassica napus each had independent origins. Results suggested that A subgenomes of B. juncea were of monophyletic origin and evolved into vegetable-use and oil-use subvarieties. Homoeolog expression dominance occurs between subgenomes of allopolyploid B. juncea, in which differentially expressed genes display more selection potential than neutral genes. Homoeolog expression dominance in B. juncea has facilitated selection of glucosinolate and lipid metabolism genes in subvarieties used as vegetables and for oil production. These homoeolog expression dominance relationships among Brassicaceae genomes have contributed to selection response, predicting the directional effects of selection in a polyploid crop genome.

The first genetic map of a synthesized allohexaploid Brassica with A, B and C genomes based on simple sequence repeat markers

We present the first genetic map of an allohexaploid Brassica species, based on segregating microsatellite markers in a doubled haploid mapping population generated from a hybrid between two hexaploid parents. This study reports the first genetic map of trigenomic Brassica. A doubled haploid mapping population consisting of 189 lines was obtained via microspore culture from a hybrid H16-1 derived from a cross between two allohexaploid Brassica lines (7H170-1 and Y54-2). Simple sequence repeat primer pairs specific to the A genome (107), B genome (44) and C genome (109) were used to construct a genetic linkage map of the population. Twenty-seven linkage groups were resolved from 274 polymorphic loci on the A genome (109), B genome (49) and C genome (116) covering a total genetic distance of 3178.8 cM with an average distance between markers of 11.60 cM. This is the first genetic framework map for the artificially synthesized Brassica allohexaploids. The linkage groups represent the expected complement of chromosomes in the A, B and C genomes from the original diploid and tetraploid parents. This framework linkage map will be valuable for QTL analysis and future genetic improvement of a new allohexaploid Brassica species, and in improving our understanding of the genetic control of meiosis in new polyploids.

Co-linearity and divergence of the A subgenome of Brassica juncea compared with other Brassica species carrying different A subgenomes

Background

There are three basic Brassica genomes (A, B, and C) and three parallel sets of subgenomes distinguished in the diploid Brassica (i.e.: B. rapa, A^rA^r; B. nigra, BⁿⁱBⁿⁱ; B. oleracea, C^oC^o) and the derived allotetraploid species (i.e.: B. juncea, A^jA^jB^jB^j; B. napus, AⁿAⁿCⁿCⁿ; B. carinata, B^cB^cC^cC^c). To understand subgenome differentiation in B. juncea in comparison to other A genome-carrying Brassica species (B. rapa and B. napus), we constructed a dense genetic linkage map of B. juncea, and conducted population genetic analysis on diverse lines of the three A-genome carrying Brassica species using a genotyping-by-sequencing approach (DArT-seq).

Results

A dense genetic linkage map of B. juncea was constructed using an F₂ population derived from Sichuan Yellow/Purple Mustard. The map included 3329 DArT-seq markers on 18 linkage groups and covered 1579 cM with an average density of two markers per cM. Based on this map and the alignment of the marker sequences with the physical genome of Arabidopsis thaliana, we observed strong co-linearity of the ancestral blocks among the different A subgenomes but also considerable block variation. Comparative analyses at the level of genome sequences of B. rapa and B. napus, and marker sequence anchored on the genetic map of B. juncea, revealed a total of 30 potential inversion events across large segments and 20 potential translocation events among the three A subgenomes. Population genetic analysis on 26 accessions of the three A genome-carrying Brassica species showed that the highest genetic distance were estimated when comparing A^j-Aⁿ than between Aⁿ-A^r and A^j-A^r subgenome pairs.

Conclusions

The development of the dense genetic linkage map of B. juncea with informative DArT-seq marker sequences and availability of the reference sequences of the A^r, and AⁿCⁿ genomes allowed us to compare the A subgenome structure of B. juncea (A^j) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (A^r, Aⁿ and A^j) but with apparent subgenomic variation. Population genetic analysis on three A-genome carrying Brassica species support the idea that B. juncea has distinct genomic diversity, and/or evolved from a different A genome progenitor of B. napus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2343-1) contains supplementary material, which is available to authorized users.

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

Brassica juncea (A^jA^jB^jB^j), is an allotetraploid that arose from two diploid species, B. rapa (A^rA^r) and B. nigra (BⁿBⁿ). It is an old oilseed crop with unique favorable traits, but the genetic improvement on this species is limited. We developed an approach to broaden its genetic base within several generations by intensive selection. The A^r subgenome from the Asian oil crop B. rapa (A^rA^r) and the B^c subgenome from the African oil crop B. carinata (B^cB^cC^cC^c) were combined in a synthesized allohexaploid (A^rA^rB^cB^cC^cC^c), which was crossed with traditional B. juncea to generate pentaploid F₁ hybrids (A^rA^jB^cB^jC^c), with subsequent self-pollination to obtain newly synthesized B. juncea (A^r/jA^r/^jB^c/jB^c/^j). After intensive cytological screening and phenotypic selection of fertility and agronomic traits, a population of new-type B. juncea was obtained and was found to be genetically stable at the F₆ generation. The new-type B. juncea possesses good fertility and rich genetic diversity and is distinctly divergent but not isolated from traditional B. juncea, as revealed by population genetic analysis with molecular markers. More than half of its genome was modified, showing exotic introgression and novel variation. In addition to the improvement in some traits of the new-type B. juncea lines, a considerable potential for heterosis was observed in inter-subgenomic hybrids between new-type B. juncea lines and traditional B. juncea accessions. The new-type B. juncea exhibited a stable chromosome number and a novel genome composition through multiple generations, providing insight into how to significantly broaden the genetic base of crops with subgenome introgression from their related species and the potential of exploring inter-subgenomic heterosis for hybrid breeding.

Microspore culture reveals complex meiotic behaviour in a trigenomic Brassica hybrid

BACKGROUND

Development of synthetic allohexaploid Brassica (2n = AABBCC) would be beneficial for agriculture, as allelic contributions from three genomes could increase hybrid vigour and broaden adaptation. Microspore culture of a near-allohexaploid hybrid derived from the cross (B. napus × B. carinata) × B. juncea was undertaken in order to assess the frequency and distribution of homologous and homoeologous crossovers in this trigenomic hybrid. SNP and SSR molecular markers were used to detect inheritance of A, B and C genome alleles in microspore-derived (MD) progeny. SNP allele copy number was also assessed. The MD progeny were also compared to progeny derived by self-pollination and open-pollination for fertility (estimated by self-pollinated seed set and pollen viability) and DNA ploidy (measured by flow cytometry).

RESULTS

In the trigenomic hybrid, homologous chromosome pairs A(j)-A(n), B(j)-B(c) and C(n)-C(c) had similar meiotic crossover frequencies and segregation to that previously observed in established Brassica species, as demonstrated by marker haplotype analysis of the MD population. Homoeologous pairing between chromosomes A1-C1, A2-C2 and A7-C6 was detected at frequencies of 12-18 %, with other homoeologous chromosome regions associating from 8 % (A3-C3) to 0-1 % (A8-C8, A8-C9) of the time. Copy number analysis revealed eight instances of additional chromosomes and 20 instances of chromosomes present in one copy in somatically doubled MD progeny. Presence of chromosome A6 was positively correlated with self-pollinated seed set and pollen viability in the MD population. Many MD progeny were unable to produce self-pollinated seed (76 %) or viable pollen (53 %), although one MD plant produced 198 self-pollinated seeds. Average fertility was significantly lower in progeny obtained by microspore culture than progeny obtained by self-pollination or open-pollination, after excluding MD progeny which had not undergone chromosome doubling.

CONCLUSIONS

Based on SNP data analysis of the microspore-derived progeny, crossover frequency per chromosome in the allohexaploid hybrid was found to be similar to that in established Brassica species, suggesting that the higher chromosome number did not significantly disrupt cellular regulation of meiosis. SNP allele copy number analysis revealed the occurrence not only of homoeologous duplication/deletion events but also other cryptic duplications and deletions that may have been the result of mitotic instability. Microspore culture simplified the assessment of chromosome behaviour in the allohexaploid hybrid but yielded progeny with lower fertility and a greater range of ploidy levels compared to progeny obtained by self- or open-pollination.

Expected Shannon Entropy and Shannon Differentiation between Subpopulations for Neutral Genes under the Finite Island Model

Shannon entropy H and related measures are increasingly used in molecular ecology and population genetics because (1) unlike measures based on heterozygosity or allele number, these measures weigh alleles in proportion to their population fraction, thus capturing a previously-ignored aspect of allele frequency distributions that may be important in many applications; (2) these measures connect directly to the rich predictive mathematics of information theory; (3) Shannon entropy is completely additive and has an explicitly hierarchical nature; and (4) Shannon entropy-based differentiation measures obey strong monotonicity properties that heterozygosity-based measures lack. We derive simple new expressions for the expected values of the Shannon entropy of the equilibrium allele distribution at a neutral locus in a single isolated population under two models of mutation: the infinite allele model and the stepwise mutation model. Surprisingly, this complex stochastic system for each model has an entropy expressable as a simple combination of well-known mathematical functions. Moreover, entropy- and heterozygosity-based measures for each model are linked by simple relationships that are shown by simulations to be approximately valid even far from equilibrium. We also identify a bridge between the two models of mutation. We apply our approach to subdivided populations which follow the finite island model, obtaining the Shannon entropy of the equilibrium allele distributions of the subpopulations and of the total population. We also derive the expected mutual information and normalized mutual information ("Shannon differentiation") between subpopulations at equilibrium, and identify the model parameters that determine them. We apply our measures to data from the common starling (Sturnus vulgaris) in Australia. Our measures provide a test for neutrality that is robust to violations of equilibrium assumptions, as verified on real world data from starlings.

Molecular characterization and population structure of the macaw palm, Acrocomia aculeata (Arecaceae), ex situ germplasm collection using microsatellites markers

The Acrocomia aculeata is one of the most promising plants for sustainable production of renewable energy. In order to understand patterns of the distribution of the allelic diversity of A. aculeata ex situ germplasm collection, the present study investigated the hypothesis that the genetic variability of the accessions may match their geographical origin. A genotypic analysis of 77 A. aculeata accessions was conducted with 6 simple sequence repeat markers. A high degree of molecular diversity among the accessions was found, with an average of 9 alleles per locus and a polymorphic information content with a mean of 0.76. A total of 4 clusters was identified by the Bayesian analysis of population structure. The highest subpopulation diversity was identified in Pop1, mainly formed by accessions from State of Mato Grosso do Sul. The populations Pop2A, Pop2B, and Pop2C, all from the State of Minas Gerais, showed high genetic variability as determined by a higher F st, and a wide genetic variance, which were identified within and among the population by analysis of molecular variance. Based on our results and on Vavilov's theory on crop origins, one possible diversity center for A. aculeata is proposed to be in a region in southeast Brazil.