Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

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

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.     

Genome-Wide Association Study Dissects the Genetic Architecture of Seed Weight and Seed Quality in Rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium^® SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using ‘pseudomolecules’ representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.     

Isozyme analysis as a tool for introgression of Sinapis alba germ plasm into Brassica napus

Summary Isozyme analysis of Brassica napus cv Topas and CGRC5006 as well as of Sinapis alba cv Emergo revealed significant polymorphism between the two species for the isozymes, aconitate hydratase, glucose phosphate isomerase, and diaphorase. F₁ hybrids between B. napus 5006 and S. alba cv Emergo were backcrossed to B. napus cv Topas, and the S₁ progeny of the first two backcrosses were studied isozymically. At the backcross one level the frequency of S. alba or S. alba plus B. napus patterns observed ranged from 18% to 87% across the four lines studied. There were differences between lines for the frequency of S. alba patterns, which could have an impact on the efficiency of selection for subsequent backcrossing. By the backcross two generation in one of the two lines studied, GR86-24, the S. alba patterns for GPI and DIA had been lost, while in the other line, GR86-28, the S. alba pattern for ACO had been lost, resulting in lost opportunity for S. alba gene transfer. In a wide cross such as S. alba x B. napus, which requires an intensive effort to accomplish, the isozymes ACO, GPI, and DIA may serve as useful markers to ensure gene transfer between the two species has occurred. In addition, the identification of lines with divergent isozyme patterns from B. napus will provide the basis for establishing linkages between S. alba traits of interest and isozyme markers.     

A triallelic genetic male sterility locus in Brassica napus: an integrative strategy for its physical mapping and possible local chromosome evolution around it

Background and Aims

Spontaneous male sterility is an advantageous trait for both constructing efficient pollination control systems and for understanding the developmental process of the male reproductive unit in many crops. A triallelic genetic male-sterile locus (BnMs5) has been identified in Brassica napus; however, its complicated genome structure has greatly hampered the isolation of this locus. The aim of this study was to physically map BnMs5 through an integrated map-based cloning strategy and analyse the local chromosomal evolution around BnMs5.

Methods

A large F₂ population was used to integrate the existing genetic maps around BnMs5. A map-based cloning strategy in combination with comparative mapping among B. napus, Arabidopsis, Brassica rapa and Brassica oleracea was employed to facilitate the identification of a target bacterial artificial chromosome (BAC) clone covering the BnMs5 locus. The genomic sequences from the Brassica species were analysed to reveal the regional chromosomal evolution around BnMs5.

Key Results

BnMs5 was finally delimited to a 0·3-cM genetic fragment from an integrated local genetic map, and was anchored on the B. napus A8 chromosome. Screening of a B. napus BAC clone library and identification of the positive clones validated that JBnB034L06 was the target BAC clone. The closest flanking markers restrict BnMs5 to a 21-kb region on JBnB034L06 containing six predicted functional genes. Good collinearity relationship around BnMs5 between several Brassica species was observed, while violent chromosomal evolutionary events including insertions/deletions, duplications and single nucleotide mutations were also found to have extensively occurred during their divergence.

Conclusions

This work represents major progress towards the molecular cloning of BnMs5, as well as presenting a powerful, integrative method to mapping loci in plants with complex genomic architecture, such as the amphidiploid B. napus.     

Genetic investigation of the origination of allopolyploid with virtually synthesized lines: application to the C subgenome of Brassica napus

Although there are a number of different allopolyploids in the plant kingdom, the exact ancestral parents of some allopolyploids have not been well characterized. We propose a strategy in which virtual allopolyploid lines derived from different types of parental species are used to investigate the progenitors of an allopolyploid. The genotypes of the parental lines and the natural allopolyploid were established using a set of DNA molecular markers. The genotypes of the virtual lines were then derived from those of the parental lines, and compared extensively with that of the natural allopolyploid. We applied this strategy to investigate the progenitors of the C subgenome of Brassica napus (rapeseed, AACC). A total of 39 accessions from 10 wild and 7 cultivated types of the B. oleracea cytodeme (CC), and 4 accessions of B. rapa (AA) were used to construct 156 virtual rapeseed lines. Genetic structure was compared among natural rapeseed, virtual rapeseed lines, and their parental lines by principal component analysis and analysis of ancestry. Our data showed that the C subgenome of natural rapeseed was related closely to the genome of cultivated B. oleracea and its related wild types, such as B. incana, B. bourgeaui, B. montana, B. oleracea ssp. oleracea and B. cretica. This finding indicated that these types or their progeny might be ancestral donors of the C subgenome of rapeseed. The successful application of the strategy of virtual allopolyploidy in rapeseed demonstrates that it can possibly be used to identify the progenitors of an allopolyploid species.     

甘蓝型油菜雌性不育突变体FS-M1乳突细胞的细胞学观察

雌性不育突变体FS-M₁是从甘蓝型油菜(Brassica napus)品种宁油10号中发现的。为了从细胞学角度研究FS-M₁的雌性不育机理, 利用荧光显微镜、扫描和透射电子显微镜观察分析了FS-M₁柱头乳突细胞的授粉行为和超微结构。结果表明: 花粉粒能在FS-M₁乳突细胞上附着和萌发形成花粉管, 但花粉管无法穿越柱头乳突细胞; 开花后的FS-M₁乳突细胞迅速衰退而呈干瘪萎蔫状, 在衰退过程中, FS-M₁柱头乳突细胞的细胞器数量减少, 细胞液泡化明显, 高尔基体、内质网和线粒体等一些细胞器结构被逐渐破坏。因此, 推测FS-M₁的雌性不育性是由于柱头乳突细胞发育异常造成的。     

Unravelling the complex trait of harvest index in rapeseed (Brassica napus L.) with association mapping

Background

Harvest index (HI), the ratio of grain yield to total biomass, is considered as a measure of biological success in partitioning assimilated photosynthate to the harvestable product. While crop production can be dramatically improved by increasing HI, the underlying molecular genetic mechanism of HI in rapeseed remains to be shown.

Results

In this study, we examined the genetic architecture of HI using 35,791 high-throughput single nucleotide polymorphisms (SNPs) genotyped by the Illumina BrassicaSNP60 Bead Chip in an association panel with 155 accessions. Five traits including plant height (PH), branch number (BN), biomass yield per plant (BY), harvest index (HI) and seed yield per plant (SY), were phenotyped in four environments. HI was found to be strongly positively correlated with SY, but negatively or not strongly correlated with PH. Model comparisons revealed that the A–D test (ADGWAS model) could perfectly balance false positives and statistical power for HI and associated traits. A total of nine SNPs on the C genome were identified to be significantly associated with HI, and five of them were identified to be simultaneously associated with HI and SY. These nine SNPs explained 3.42 % of the phenotypic variance in HI.

Conclusions

Our results showed that HI is a complex polygenic phenomenon that is strongly influenced by both environmental and genotype factors. The implications of these results are that HI can be increased by decreasing PH or reducing inefficient transport from pods to seeds in rapeseed. The results from this association mapping study can contribute to a better understanding of natural variations of HI, and facilitate marker-based breeding for HI.

Electronic supplementary material

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

A comparison of anther and microspore culture as a breeding tool in Brassica napus

Summary A direct comparison of microspore culture and anther culture was made in Brassica napus using F₁ crosses of Regent (canola) by Golden (rapeseed), and their reciprocals, as well as a hybrid between Reston and a highly embryogenic, canola-quality breeding line (G231) as donor plants. The study confirmed that microspore culture can be ten times more efficient than anther culture for embryo production. Embryo yields from cultures initiated from the Reston x G231 were four-fold greater than those initiated from the Regent x Golden crosses, and significant differences were also detected among cultures initiated from the different Regent x Golden crosses. These results illustrate the influence that donor plant genotype has on embryo production. However, superior embryogenic potential among donor material was not always coincident with superior plant production. The average haploid-todiploid ratio in microspore-derived regenerates was 21 for the population obtained from the Regent x Golden crosses but 11 for the Reston x G231 cross. For both types of material, the frequency of diploids increased upon repeated cycles of explanting. A field study showed that there were no differences between the populations of anther-derived and microspore-derived spontaneous diploid and doubled haploid lines, with respect to the days required for them to flower or to mature. The information is valuable for canola breeding programs considering the use of haploidy.     

Changes in fitness-associated traits due to the stacking of transgenic glyphosate resistance and insect resistance in Brassica napus L

Increasingly, genetically modified crops are being developed to express multiple 'stacked' traits for different types of transgenes, for example, herbicide resistance, insect resistance, crop quality and tolerance to environmental stresses. The release of crops that express multiple traits could result in ecological changes in weedy environments if feral crop plants or hybrids formed with compatible weeds results in more competitive plants outside of agriculture. To examine the effects of combining transgenes, we developed a stacked line of canola (Brassica napus L.) from a segregating F(2) population that expresses both transgenic glyphosate resistance (CP4 EPSPS) and lepidopteran insect resistance (Cry1Ac). Fitness-associated traits were evaluated between this stacked genotype and five other Brassica genotypes in constructed mesocosm plant communities exposed to insect herbivores (Plutella xylostella L.) or glyphosate-drift. Vegetative biomass, seed production and relative fecundity were all reduced in stacked trait plants when compared with non-transgenic plants in control treatments, indicating potential costs of expressing multiple transgenes without selection pressure. Although costs of the transgenes were offset by selective treatment, the stacked genotype continued to produce fewer seeds than either single transgenic line. However, the increase in fitness of the stacked genotype under selective pressure contributed to an increased number of seeds within the mesocosm community carrying unselected, hitchhiking transgenes. These results demonstrate that the stacking of these transgenes in canola results in fitness costs and benefits that are dependent on the type and strength of selection pressure, and could also contribute to changes in plant communities through hitchhiking of unselected traits.     

Reduced Susceptibility of Brassica napus to Pratylenchus neglectus in Plants with Elevated Root Levels of 2-Phenylethyl Glucosinolate

The effect of canola (Brassica napus) as a crop suppressive to Pratylenchus neglectus is in part due to the release of nematicidal isothiocyanates, particularly 2-phenylethyl isothiocyanate, from degrading root tissues. However, many cultivars of canola are relatively susceptible to P. neglectus and will fail to reduce soil populations of the nematode. A survey of B. napus accessions and closely related species revealed limited scope to decrease the susceptibility of canola through conventional intercrossing. Susceptibility to P. neglectus was not related to the total glucosinolate levels, but there were signifi- cant, negative correlations (r = -0.619, -0.517; P < 0.001) between root levels of 2-phenylethyl glucosinolate (isothiocyanate precursor) and plant susceptibility to P. neglectus: plants containing more than a certain threshold level of 2-phenylethyl glucosinolate showed reduced susceptibility to the nematode. Selection for high root levels of 2-phenylethyl glucosinolate should reduce the susceptibility of the plants during the growing season while also increasing the nematicidal impact of the degrading root tissues, thereby improving the suppressive benefits of the crop when used in rotation with cereals.     

Mining expressed sequence tags of rapeseed (Brassica napus L.) to predict the drought responsive regulatory network

Physiology and Molecular Biology of Plants - It is of great significance to understand the regulatory mechanisms by which plants deal with drought stress. Two EST libraries derived from rapeseed...     

EST-based in silico identification and in vitro test of antimicrobial peptides in Brassica napus

Background

Brassica napus is the third leading source of vegetable oil in the world after soybean and oil palm. The accumulation of gene sequences, especially expressed sequence tags (ESTs) from plant cDNA libraries, has provided a rich resource for genes discovery including potential antimicrobial peptides (AMPs). In this study, we used ESTs including those generated from B. napus cDNA libraries of seeds, pathogen-challenged leaves and deposited in the public databases, as a model, to perform in silico identification and consequently in vitro confirmation of putative AMP activities through a highly efficient system of recombinant AMP prokaryotic expression.

Results

In total, 35,788 were generated from cDNA libraries of pathogen-challenged leaves and 187,272 ESTs from seeds of B. napus, and the 644,998 ESTs of B. napus were downloaded from the EST database of PlantGDB. They formed 201,200 unigenes. First, all the known AMPs from the AMP databank (APD2 database) were individually queried against all the unigenes using the BLASTX program. A total of 972 unigenes that matched the 27 known AMP sequences in APD2 database were extracted and annotated using Blast2GO program. Among these unigenes, 237 unigenes from B. napus pathogen-challenged leaves had the highest ratio (1.15 %) in this unigene dataset, which is 13 times that of the unigene datasets of B. napus seeds (0.09 %) and 2.3 times that of the public EST dataset. About 87 % of each EST library was lipid-transfer protein (LTP) (32 % of total unigenes), defensin, histone, endochitinase, and gibberellin-regulated proteins. The most abundant unigenes in the leaf library were endochitinase and defensin, and LTP and histone in the pub EST library. After masking of the repeat sequence, 606 peptides that were orthologous matched to different AMP families were found. The phylogeny and conserved structural motifs of seven AMPs families were also analysed. To investigate the antimicrobial activities of the predicted peptides, 31 potential AMP genes belonging to different AMP families were selected to test their antimicrobial activities after bioinformatics identification. The AMP genes were all optimized according to Escherichia coli codon usage and synthetized through one-step polymerase chain reaction method. The results showed that 28 recombinant AMPs displayed expected antimicrobial activities against E. coli and Micrococcus luteus and Sclerotinia sclerotiorum strains.

Conclusion

The study not only significantly expanded the number of known/predicted peptides, but also contributed to long-term plant genetic improvement for increased resistance to diverse pathogens of B.napus. These results proved that the high-throughput method developed that combined an in silico procedure with a recombinant AMP prokaryotic expression system is considerably efficient for identification of new AMPs from genome or EST sequence databases.

Electronic supplementary material

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

Construction of a bacterial artificial chromosome library for the Rongchang pig breed and its use for the identification of genes involved in intramuscular fat deposition

In a search for genes affecting intramuscular fat deposition, we constructed a bacterial artificial chromosome (BAC) library for the whole genome of Rongchang pig, a domestic Chinese swine breed. The library consisted of approximately 192,000 clones, with an averaged insert size of 116 kb. Frequency of non-insert clone of the BAC library was no higher than 1.8%, based on estimation of 220 BAC clones randomly selected. We estimated the coverage of the library to be more than seven porcine genome equivalents. Subsequent screening of the BAC library with a three-step PCR procedure resulted in identification of seven candidate genes that were potentially involved in intramuscular fat deposition. The number of positive BAC clones ranged from 2 to 4 for each of the seven genes. One positive clone, containing the lipin1 gene, was fully sequenced by shotgun method to generate 118,041 bp porcine genomic sequences. The BAC clone contained complete DNA sequence of porcine lipin1 gene including all the exons and introns. Our results indicate that this BAC library is a useful tool for gene identification and help to serve as an important resource for future porcine genomic study.     

Quantitative trait loci for seed yield and yield-related traits, and their responses to reduced phosphorus supply in Brassica napus

Background and Aims

One of the key targets of breeding programmes in rapeseed (Brassica napus) is to develop high-yield varieties. However, the lack of available phosphorus (P) in soils seriously limits rapeseed production. The aim of this study was to dissect the genetic control of seed yield and yield-related traits in B. napus grown with contrasting P supplies.

Methods

Two-year field trials were conducted at one site with normal and low P treatments using a population of 124 recombinant inbred lines derived from a cross between ‘B104-2’ and ‘Eyou Changjia’. Seed yield, seed weight, seed number, pod number, plant height, branch number and P efficiency coefficient (PEC) were investigated. Quantitative trait locus (QTL) analysis was performed by composite interval mapping.

Key Results

The phenotypic values of most of the tested traits were reduced under the low P conditions. In total, 74 putative QTLs were identified, contributing 7·3–25·4 % of the phenotypic variation. Of these QTLs, 16 (21·6 %) were detected in two seasons and in the mean value of two seasons, and eight QTLs for two traits were conserved across P levels. Low-P-specific QTLs were clustered on chromosomes A1, A6 and A8. By comparative mapping between Arabidopsis and B. napus, 161 orthologues of 146 genes involved in Arabidopsis P homeostasis and/or yield-related trait control were associated with 45 QTLs corresponding to 23 chromosomal regions. Four gene-based markers developed from genes involved in Arabidopsis P homeostasis were mapped to QTL intervals.

Conclusions

Different genetic determinants were involved in controlling seed yield and yield-related traits in B. napus under normal and low P conditions. The QTLs detected under reduced P supply may provide useful information for improving the seed yield of B. napus in soils with low P availability in marker-assisted selection.     

Role of six single nucleotide polymorphisms,risk factors in coronary disease,in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.