A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING

Targeting Induced Local Lesions in Genomes (TILLING): advances and opportunities for fast tracking crop breeding

The intensification of food production via conventional crop breeding alone is inadequate to cater for global hunger. The development of precise and expeditious high throughput reverse genetics approaches has hugely benefited modern plant breeding programs. Targeting Induced Local Lesions in Genomes (TILLING) is one such reverse genetics approach which employs chemical/physical mutagenesis to create new genetic sources and identifies superior/novel alleles. Owing to technical limitations and sectional applicability of the original TILLING protocol, it has been timely modified. Successions include: EcoTILLING, Double stranded EcoTILLING (DEcoTILLING), Self-EcoTILLING, Individualized TILLING (iTILLING), Deletion-TILLING (De-TILLING), PolyTILLING, and VeggieTILLING. This has widened its application to a variety of crops and needs. They can characterize mutations in coding as well as non-coding regions and can overcome complexities associated with the large genomes. Combining next generation sequencing tools with the existing TILLING protocols has enabled screening of huge germplasm collections and mutant populations for the target genes. In silico TILLING platforms have transformed TILLING into an exciting breeding approach. The present review outlines these multifarious TILLING modifications for precise mutation detection and their application in advance breeding programmes together with relevant case studies. Appropriate use of these protocols will open up new avenues for crop improvement in the twenty first century.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Plant Functional Genomics Based on High-Throughput CRISPR Library Knockout Screening: A Perspective

Plant biology studies in the post-genome era have been focused on annotating genome sequences' functions. The established plant mutant collections have greatly accelerated functional genomics research in the past few decades. However, most plant genome sequences' roles and the underlying regulatory networks remain substantially unknown. Clustered, regularly interspaced short palindromic repeat (CRISPR)-associated systems are robust, versatile tools for manipulating plant genomes with various targeted DNA perturbations, providing an excellent opportunity for high-throughput interrogation of DNA elements' roles. This study compares methods frequently used for plant functional genomics and then discusses different DNA multi-targeted strategies to overcome gene redundancy using the CRISPR-Cas9 system. Next, this work summarizes recent reports using CRISPR libraries for high-throughput gene knockout and function discoveries in plants. Finally, this work envisions the future perspective of optimizing and leveraging CRISPR library screening in plant genomes' other uncharacterized DNA sequences.

Is it the end of TILLING era in plant science?

Since its introduction in 2000, the TILLING strategy has been widely used in plant research to create novel genetic diversity. TILLING is based on chemical or physical mutagenesis followed by the rapid identification of mutations within genes of interest. TILLING mutants may be used for functional analysis of genes and being nontransgenic, they may be directly used in pre-breeding programs. Nevertheless, classical mutagenesis is a random process, giving rise to mutations all over the genome. Therefore TILLING mutants carry background mutations, some of which may affect the phenotype and should be eliminated, which is often time-consuming. Recently, new strategies of targeted genome editing, including CRISPR/Cas9-based methods, have been developed and optimized for many plant species. These methods precisely target only genes of interest and produce very few off-targets. Thus, the question arises: is it the end of TILLING era in plant studies? In this review, we recap the basics of the TILLING strategy, summarize the current status of plant TILLING research and present recent TILLING achievements. Based on these reports, we conclude that TILLING still plays an important role in plant research as a valuable tool for generating genetic variation for genomics and breeding projects.

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus

BACKGROUND

Increasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content.

RESULTS

Here, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC.

CONCLUSIONS

The present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

Variations in Total Protein and Amino Acids in the Sequenced Sorghum Mutant Library

Sorghum (Sorghum bicolor) is the fifth most important cereal crop worldwide; however, its utilization in food products can be limited due to reduced nutritional quality related to amino acid composition and protein digestibility in cooked products. Low essential amino acid levels and digestibility are influenced by the composition of the sorghum seed storage proteins, kafirins. In this study, we report a core collection of 206 sorghum mutant lines with altered seed storage proteins. Wet lab chemistry analysis was conducted to evaluate the total protein content and 23 amino acids, including 19 protein-bound and 4 non-protein amino acids. We identified mutant lines with diverse compositions of essential and non-essential amino acids. The highest total protein content in these lines was almost double that of the wild-type (BTx623). The mutants identified in this study can be used as a genetic resource to improve the sorghum grain quality and determine the molecular mechanisms underlying the biosynthesis of storage protein and starch in sorghum seeds.

Investigation of mutation load and rate in androgenic mutant lines of rapeseed in early generations evaluated by high-density SNP genotyping

Oilseed rape (Brassica napus) is an important oil crop distributed worldwide with a broad adaptation to different climate zones. The cultivation of rapeseed is one of the most commercially viable areas in crop production. Altogether 269,093 ha of rapeseed are cultivated in Kazakhstan. However, all rapeseed cultivars and lines cultivated in Kazakhstan on an industrial scale predominantly belong to the foreign breeding system. Therefore, the formation of a diverse genetic pool for breeding new, highly productive cultivars adopted to the environmental conditions of Kazakhstan is the most important goal in country selection programs. In this work, we have developed ethyl methanesulfonate (EMS) doubled haploid mutant lines from plant material of cultivars 'Galant' and 'Kris' to broad diversity of rapeseed in Kazakhstan. The development of mutant lines was performed via embryo callusogenesis or embryo secondary callusogenesis. Mutants were investigated by Brassica90k SNP array, and we were able to locate 24,657 SNPs from 26,256 SNPs filtered by quality control on the genome assembly (Bra_napus_v2.0). Only 18,831 SNPs were assigned to the available annotated genomic features. The most frequent combination of mutations according to reference controls was adenine with guanine (70%), followed by adenine with cytosine (28.8%), and only minor fractions were cytosine with guanine (0.54%) and adenine with thymine (0.59%). We revealed 5606.27 markers for 'Kris' and 4893.01 markers for 'Galant' by mutation occurrence. Most mutation occurrences were occupied by double mutations where progenitors and offspring were homozygous by different alleles, enabling the selection of appropriate genotypes in a short period of time. Regarding the biological impact of mutations, 861 variants were reported as having a low predicted impact, with 1042 as moderate and 121 as high; all others were reported as belonging to non-coding sequences, intergenic regions, and other features with the effect of modifiers. Protein encoding genes, such as wall-associated receptor kinase-like protein 5, TAO1-like disease resistance protein, receptor-like protein 12, and At5g42460-like F-box protein, contained more than two variable positions, with an impact on their biological activities. Nevertheless, the obtained mutant lines were able to survive and reproduce. Mutant lines, which include moderate and high impact mutations in encoding genes, are a perfect pool not only for MAS but also for the investigation of the fundamental basis of protein functions. For the first time, a collection of mutant lines was developed in our country to improve the selection of local rapeseed cultivars.

Reduced glucosinolate content in oilseed rape (Brassica napus L.) by random mutagenesis of BnMYB28 and BnCYP79F1 genes

The presence of anti-nutritive compounds like glucosinolates (GSLs) in the rapeseed meal severely restricts its utilization as animal feed. Therefore, reducing the GSL content to < 18 µmol/g dry weight in the seeds is a major breeding target. While candidate genes involved in the biosynthesis of GSLs have been described in rapeseed, comprehensive functional analyses are missing. By knocking out the aliphatic GSL biosynthesis genes BnMYB28 and BnCYP79F1 encoding an R2R3 MYB transcription factor and a cytochrome P450 enzyme, respectively, we aimed to reduce the seed GSL content in rapeseed. After expression analyses on single paralogs, we used an ethyl methanesulfonate (EMS) treated population of the inbred winter rapeseed 'Express617' to detect functional mutations in the two gene families. Our results provide the first functional analysis by knock-out for the two GSL biosynthesis genes in winter rapeseed. We demonstrate that independent knock-out mutants of the two genes possessed significantly reduced seed aliphatic GSLs, primarily progoitrin. Compared to the wildtype Express617 control plants (36.3 µmol/g DW), progoitrin levels were decreased by 55.3% and 32.4% in functional mutants of BnMYB28 (16.20 µmol/g DW) and BnCYP79F1 (24.5 µmol/g DW), respectively. Our study provides a strong basis for breeding rapeseed with improved meal quality in the future.

Direct access to millions of mutations by whole genome sequencing of an oilseed rape mutant population

Induced mutations are an essential source of genetic variation in plant breeding. Ethyl methanesulfonate (EMS) mutagenesis has been frequently applied, and mutants have been detected by phenotypic or genotypic screening of large populations. In the present study, a rapeseed M₂ population was derived from M₁ parent cultivar 'Express' treated with EMS. Whole genomes were sequenced from fourfold (4×) pools of 1988 M₂ plants representing 497 M₂ families. Detected mutations were not evenly distributed and displayed distinct patterns across the 19 chromosomes with lower mutation rates towards the ends. Mutation frequencies ranged from 32/Mb to 48/Mb. On average, 284 442 single nucleotide polymorphisms (SNPs) per M₂ DNA pool were found resulting from EMS mutagenesis. 55% of the SNPs were C → T and G → A transitions, characteristic for EMS induced ('canonical') mutations, whereas the remaining SNPs were 'non-canonical' transitions (15%) or transversions (30%). Additionally, we detected 88 725 high confidence insertions and deletions per pool. On average, each M₂ plant carried 39 120 canonical mutations, corresponding to a frequency of one mutation per 23.6 kb. Approximately 82% of such mutations were located either 5 kb upstream or downstream (56%) of gene coding regions or within intergenic regions (26%). The remaining 18% were located within regions coding for genes. All mutations detected by whole genome sequencing could be verified by comparison with known mutations. Furthermore, all sequences are accessible via the online tool 'EMSBrassica' (http://www.emsbrassica.plantbreeding.uni-kiel.de), which enables direct identification of mutations in any target sequence. The sequence resource described here will further add value for functional gene studies in rapeseed breeding.

Acetylome reprograming participates in the establishment of fruit metabolism during polyploidization in citrus

Polyploidization leads to novel phenotypes and is a major force in evolution. However, the relationship between the evolution of new traits and variations in the post-translational modifications (PTM) of proteins during polyploidization has not been studied. Acetylation of lysine residues is a common protein PTM that plays a critical regulatory role in central metabolism. To test whether changes in metabolism in citrus fruit is associated with the reprogramming of lysine acetylation (Kac) in non-histone proteins during allotetraploidization, we performed a global acetylome analysis of fruits from a synthetic allotetraploid citrus and its diploid parents. A total of 4,175 Kac sites were identified on 1,640 proteins involved in a wide range of fruit traits. In the allotetraploid, parental dominance (i.e. resemblance to one of the two parents) in specific fruit traits, such as fruit acidity and flavonol metabolism, was highly associated with parental Kac level dominance in pertinent enzymes. This association is due to Kac-mediated regulation of enzyme activity. Moreover, protein Kac probably contributes to the discordance between the transcriptomic and proteomic variations during allotetraploidization. The acetylome reprogramming can be partially explained by the expression pattern of several lysine deacetylases (KDACs). Overexpression of silent information regulator 2 (CgSRT2) and histone deacetylase 8 (CgHDA8) diverted metabolic flux from primary metabolism to secondary metabolism and partially restored a metabolic status to the allotetraploid, which expressed attenuated levels of CgSRT2 and CgHDA8. Additionally, KDAC inhibitor treatment greatly altered metabolism in citrus fruit. Collectively, these findings reveal the important role of acetylome reprogramming in trait evolution during polyploidization.

Perspectives for integrated insect pest protection in oilseed rape breeding

In the past, breeding for incorporation of insect pest resistance or tolerance into cultivars for use in integrated pest management schemes in oilseed rape/canola (Brassica napus) production has hardly ever been approached. This has been largely due to the broad availability of insecticides and the complexity of dealing with high-throughput phenotyping of insect performance and plant damage parameters. However, recent changes in the political framework in many countries demand future sustainable crop protection which makes breeding approaches for crop protection as a measure for pest insect control attractive again. At the same time, new camera-based tracking technologies, new knowledge-based genomic technologies and new scientific insights into the ecology of insect-Brassica interactions are becoming available. Here we discuss and prioritise promising breeding strategies and direct and indirect breeding targets, and their time-perspective for future realisation in integrated insect pest protection of oilseed rape. In conclusion, researchers and oilseed rape breeders can nowadays benefit from an array of new technologies which in combination will accelerate the development of improved oilseed rape cultivars with multiple insect pest resistances/tolerances in the near future.

A Review of Erucic Acid Production in Brassicaceae Oilseeds: Progress and Prospects for the Genetic Engineering of High and Low-Erucic Acid Rapeseeds (Brassica napus)

Erucic acid (C22:1, ω-9, EA) is a very-long-chain monounsaturated fatty acid (FA) that is an important oleochemical product with a wide range of uses in metallurgy, machinery, rubber, the chemical industry, and other fields because of its hydrophobicity and water resistance. EA is not easily digested and absorbed in the human body, and high-EA rapeseed (HEAR) oil often contains glucosinolates. Both glucosinolates and EA are detrimental to health and can lead to disease, which has resulted in strict guidelines by regulatory bodies on maximum EA contents in oils. Increasingly, researchers have attempted to enhance the EA content in Brassicaceae oilseeds to serve industrial applications while conversely reducing the EA content to ensure food safety. For the production of both LEAR and HEAR, biotechnology is likely to play a fundamental role. Elucidating the metabolic pathways of EA can help inform the improvement of Brassicaceae oilseeds through transgenic technology. In this paper, we introduce the industrial applications of HEAR oil and health benefits of low-EA rapeseed (LEAR) oil first, following which we review the biosynthetic pathways of EA, introduce the EA resources from plants, and focus on research related to the genetic engineering of EA in Brassicaceae oilseeds. In addition, the effects of the environment on EA production are addressed, and the safe cultivation of HEAR and LEAR is discussed. This paper supports further research into improving FAs in Brassicaceae oilseeds through transgenic technologies and molecular breeding techniques, thereby advancing the commercialization of transgenic products for better application in various fields.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

Comparison between Germinated Seed and Isolated Microspore EMS Mutagenesis in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Mutagenesis is an important tool for breeding and genomic research. In this study, the germinated seeds and isolated microspores of a double haploid line ‘FT’ were treated with EMS, respectively, with the aim of comparing the effects of the two approaches on generating mutants in Chinese cabbage. For microspore EMS mutagenesis, the isolated microspores were treated with 0.12% EMS for 20 min, a total of 1268 plantlets were obtained, and 15 M1 mutants were screened with a mutation frequency of 1.2%. For seed EMS mutagenesis, 7800 germinated seeds were treated with 0.8% EMS for 12 h, and a total of 701 M2 mutants were screened, with a mutation frequency of 18.78%. In total, 716 mutants with heritable morphological variation including leaf color, leaf shape, leafy head, bolting, and fertility, were obtained from the EMS mutagenesis experiments. Homozygous mutant plants could be screened from M1 lines by microspore mutagenesis, and M2 lines by seed mutagenesis. The mutation frequency was higher in seed mutagenesis than in microspore mutagenesis. Based on these results, we propose that seed EMS mutagenesis is more suitable to generate a large-scale mutant library, and the microspore EMS mutagenesis is conducive to rapidly obtaining homozygous mutants.

CRISPR/Cas9-Targeted Mutagenesis of BnaFAE1 Genes Confers Low-Erucic Acid in Brassica napus

Rapeseed (Brassica napus) is an important oilseed crop widely planted in the world, providing substantial edible oil and other nutrients for mankind. The composition of fatty acids affects the edible and processing quality of vegetable oils, among which erucic acid (EA) is potentially to cause health problems. Therefore, low erucic acid (LEA) has always been a breeding trait of B. napus. Fatty acid elongase 1 (FAE1) plays a decisive role in the synthesis of EA. There are two functional homologous copies of FAE1 on the A08 and C03 chromosomes in B. napus. In this study, we used CRISPR/Cas9 technology to create targeted mutations on these two homologous copies of BnaFAE1 in three B. napus germplasms with high EA (>30%) and high oil (>50%). Our results show that the EA content was significantly reduced by more than 10 percentage points in the mutant of BnaC03.FAE1 (c03), while the double mutation of BnaA08.FAE1 and BnaC03.FAE1 (a08c03) resulted in nearly zero EA in three BnaFAE1-edited germplasms, and the oleic acid content was increased in different degrees. In addition, knockout of BnaA08.FAE1 or/and BnaC03.FAE1 mildly decreased seed oil content, but had no significant effect on other agronomic traits. In general, we successfully created low EA germplasms of B. napus, which provides a feasible way for future low EA breeding.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Gene duplication and stress genomics in Brassicas: Current understanding and future prospects

Polyploidy or whole genome duplication (WGD) is an evolutionary phenomenon that happened in all angiosperms multiple times over millions of years. Extensive studies on the model plant Arabidopsis thaliana genome have revealed that it has undergone five rounds of WGDs followed, in the Brassicaceae tribe, by a characteristic whole genome triplication (WGT). In addition, small-scale events such as tandem or segmental duplications and retrotransposition also enable plants to reshape their genomes. Over the decades, extensive research efforts have been undertaken to understand the evolutionary significance of polyploidy. On the other hand, much less attention has been paid to understanding the impact of gene duplication on the diversification of important stress response genes. The main objective of this review is to discuss key aspects of gene and genome duplications with a focus on genes primarily regulated by osmotic stresses. The focal family is the Brassicaceae, since it (i) underwent multiple rounds of WGDs plus WGTs, (ii) hosts many economically important crops and wild relatives that are tolerant to a range of stresses, and (iii) comprises many species that have already been sequenced. Diverse molecular mechanisms that lead to structural and regulatory alterations of duplicated genes are discussed. Examples are drawn from recent literature to elucidate expanded, stress responsive gene families identified from different Brassica crops. A combined bioinformatic and transcriptomic method has been proposed and tested on a known stress-responsive gene pair to prove that stress-responsive duplicated allelic variants can be identified by this method. Finally, future prospects for engineering these genes into crops to enhance stress tolerance are discussed, and important resources for Brassica genome research are provided.

Development and screening of EMS mutants with altered seed oil content or fatty acid composition in Brassica napus

Brassica napus is an important oilseed crop in the world, and the mechanism of seed oil biosynthesis in B. napus remains unclear. In order to study the mechanism of oil biosynthesis and generate germplasms for breeding, an ethyl methanesulfonate (EMS) mutant population with ~100 000 M2 lines was generated using Zhongshuang 11 as the parent line. The EMS-induced genome-wide mutations in M2-M4 plants were assessed. The average number of mutations including single nucleotide polymorphisms and insertion/deletion in M2-M4 was 21 177, 28 675 and 17 915, respectively. The effects of the mutations on gene function were predicted in M2-M4 mutants, respectively. We screened the seeds from 98 113 M2 lines, and 9415 seed oil content and fatty acid mutants were identified. We further confirmed 686 mutants with altered seed oil content and fatty acid in advanced generation (M4 seeds). Five representative M4 mutants with increased oleic acid were re-sequenced, and the potential causal variations in FAD2 and ROD1 genes were identified. This study generated and screened a large scale of B. napus EMS mutant population, and the identified mutants could provide useful genetic resources for the study of oil biosynthesis and genetic improvement of seed oil content and fatty acid composition of B. napus in the future.

Molecular and genetic bases of heat stress responses in crop plants and breeding for increased resilience and productivity

To ensure the food security of future generations and to address the challenge of the 'no hunger zone' proposed by the FAO (Food and Agriculture Organization), crop production must be doubled by 2050, but environmental stresses are counteracting this goal. Heat stress in particular is affecting agricultural crops more frequently and more severely. Since the discovery of the physiological, molecular, and genetic bases of heat stress responses, cultivated plants have become the subject of intense research on how they may avoid or tolerate heat stress by either using natural genetic variation or creating new variation with DNA technologies, mutational breeding, or genome editing. This review reports current understanding of the genetic and molecular bases of heat stress in crops together with recent approaches to creating heat-tolerant varieties. Research is close to a breakthrough of global relevance, breeding plants fitter to face the biggest challenge of our time.

Combining QTL Mapping with Genome Resequencing Identifies an Indel in an R Gene that is Associated with Variation in Leaf Rust Disease Resistance in Poplar

Poplar trees (Populus spp.) are important and are widely grown worldwide. However, the extensive occurrence of leaf rust disease caused by Melampsora spp. seriously inhibits their growth and reduces their biomass. In our previous study, a high-quality genetic map was constructed for the poplar F1 population I-69 × XYY by using next-generation sequencing-based genotyping-by-sequencing. Here, we collected phenotypic data on leaf rust disease resistance on three different dates for all 300 progenies of the F1 population. Combining a high-quality genetic map and phenotypic data, we were able to detect 11 major quantitative trait loci (QTLs) for leaf rust disease resistance. Among these 11 QTLs, two pairs were detected on at least two dates. In the corresponding genomic sequence, we found that resistance (R) gene clusters were located in these two QTL regions. By using genome resequencing, PCR confirmation and statistical analysis, a 611-bp deletion within an R gene in one QTL region was found to be associated with variation in leaf rust disease resistance. A PCR-based examination of this 611-bp deletion was performed. This 611-bp deletion was also found to affect mRNA splicing and form a new protein with the loss of some key protein domains. Based on this study, we were able to determine the genetic architecture of variation in poplar leaf rust disease resistance, and the 611-bp deletion in the R gene could be used as a diagnostic marker for future poplar molecular breeding.

Identification of phytic acid mutants in oilseed rape (Brassica napus) by large-scale screening of mutant populations through amplicon sequencing

Brassica napus (oilseed rape) is an important oil crop in temperate regions, which originated from hybridization of Brassica oleracea and Brassica rapa. Owing to its polyploidy, the functional study of single genes is cumbersome. Phytic acid is considered as an antinutritive compound, and we aimed to knock out the underlying synthesis and transporter genes to identify low phytic acid mutants. We implemented a high-throughput next-generation sequencing screening protocol for an ethylmethane sulfonate population of 7680 plants in six gene families (BnMIPS, BnMIK, Bn2-PGK, BnIPK1, BnIPK2, and BnMRP5) with two paralogues for each gene. A total of 1487 mutations were revealed, and the vast majority (96%) were confirmed by Sanger sequencing. Furthermore, the characterization of double mutants of Bn.2-PGK2 showed a significant reduction of phytic acid contents. We propose to use three-dimensional pooling combined with amplicon stacking and next-generation sequencing to identify mutations in polyploid oilseed rape in a fast and cost-effective manner for complex metabolic pathways. Furthermore, the mutants identified in Bn2-PGK2 might be a very valuable resource for industrial production of oilseed rape protein for human consumption.

An efficient TILLING platform for cultivated tobacco

Targeting-induced local lesions in genomes (TILLING) is a powerful reverse-genetics tool that enables high-throughput screening of genomic variations in plants. Although TILLING has been developed for many diploid plants, the technology has been used in very few polyploid species due to their genomic complexity. Here, we established an efficient capillary electrophoresis-based TILLING platform for allotetraploid cultivated tobacco (Nicotiana tabacum L.) using an ethyl methanesulfonate (EMS)-mutagenized population of 1,536 individuals. We optimized the procedures for endonuclease preparation, leaf tissue sampling, DNA extraction, normalization, pooling, PCR amplification, heteroduplex formation, and capillary electrophoresis. In a test screen using seven target genes with eight PCR fragments, we obtained 118 mutants. The mutation density was estimated to be approximately one mutation per 106 kb on average. Phenotypic analyses showed that mutations in two heavy metal transporter genes, HMA2S and HMA4T, led to reduced accumulation of cadmium and zinc, which was confirmed independently using CRISPR/Cas9 to generate knockout mutants. Our results demonstrate that this powerful TILLING platform (available at http://www.croptilling.org) can be used in tobacco to facilitate functional genomics applications.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

Phenotypic, biochemical and genomic variability in generations of the rapeseed (Brassica napus L.) mutant lines obtained via chemical mutagenesis

The phenotypic, biochemical and genetic variability was studied in M2-M5 generations of ethyl methansulfonat (EMS, 0.2%) mutagenized rapeseed lines generated from canola, '00', B. napus cv. Vikros. EMS mutagenesis induced extensive diversity in morphological and agronomic traits among mutant progeny resulted in selection of EMS populations of B. napus- and B. rapa-morphotypes. The seeds of the obtained mutant lines were high-protein, low in oil and stabilized in contents of main fatty acids which make them useful for feed production. Despite the increased level of various meiotic abnormalities revealed in EMS populations, comparative karyotype analysis and FISH-based visualization of 45S and 5S rDNA indicated a high level of karyotypic stability in M2-M5 plants, and therefore, the obtained mutant lines could be useful in further rapeseed improvement. The revealed structural chromosomal reorganizations in karyotypes of several plants of B. rapa-type indicate that rapeseed breeding by chemical mutagenesis can result in cytogenetic instability in the mutant progeny, and therefore, it should include the karyotype examination. Our findings demonstrate that EMS at low concentrations has great potential in rapeseed improvement.

Identification of SNP loci and candidate genes related to four important fatty acid composition in Brassica napus using genome wide association study

Rapeseed oil (canola, Brassica napus L.) is an important healthy vegetable oil throughout the world, the nutritional and economical value of which largely depends on its seed fatty acid composition. In this study, based on 201,187 SNP markers developed from the SLAF-seq (specific locus amplified fragment sequencing), a genome wide association study of four important fatty acid content traits (erucic acid, oleic acid, linoleic acid and linolenic acid) in a panel of 300 inbred lines of rapeseed in two environments (JXAU and JXRIS) was carried out. A total of 148 SNP loci significantly associated with these traits were detected by MLM model analysis respectively, and 30 SNP loci on A08 and C03 chromosomes were detected in three traits of erucic acid, oleic acid and linoleic acid contents simultaneously. Furthermore, 108 highly favorable alleles for increasing oleic acid and linoleic acid content, also for decreasing erucic acid content simultaneously were observed. By a basic local alignment search tool (BLAST) search with in a distance of 100 Kb around these significantly SNP-trait associations, we identified 20 orthologs of the functional candidate genes related to fatty acid biosynthesis, including the known vital fatty acid biosynthesis genes of BnaA.FAE1 and BnaC. FAE1 on the A08 and C03 chromosomes, and other potential candidate genes involving in the fatty acid biosynthesis pathway, such as the orthologs genes of FAD2, LACS09, KCS17, CER4, TT16 and ACBP5. This study lays a basis for uncovering the genetic variations and the improvement of fatty acid composition in B. napus.

A missense mutation of STERILE APETALA leads to female sterility in Chinese cabbage (Brassica campestris ssp. pekinensis)

Flower development is essential for the sexual reproduction and crop yield of plants. Thus, a better understanding of plant sterility from the perspective of morphological and molecular genetics is imperative. In our previous study, a recessive female-sterile Chinese cabbage mutant fsm was obtained from a doubled haploid line 'FT' via an isolated microspore culture combined with EMS mutagenesis. Pistil aniline blue staining and stigma scanning observation showed that the growth of the stigma papillar cells and pollen tubes of the mutant fsm were normal. Therefore, the female sterility was due to abnormal development of the ovules. To map the mutant fsm, 3108 F₂ individuals were selected for linkage analysis. Two closely linked markers, Indel-I2 and Indel-I7, were localized on the flanking region of fsm at distances of 0.05 cM and 0.06 cM, respectively. The physical distance between Indel-I2 and Indel-I7 was ~ 1376 kb, with 107 genes remaining in the target region. This region was located on the chromosome A04 centromere, on which low recombination rates and a high frequency of repetitive sequences were present. Whole-genome re-sequencing detected a single-nucleotide (C-to-A) transition (TCG/TAG) on the exon of BraA04001030, resulting in a premature stop codon. Genotyping revealed that the female-sterile phenotype was fully cosegregated with this SNP. BraA04001030 encodes a homologue of STERILE APETALA (SAP) transcriptional regulator, which plays vital roles in floral development. The results of the present study suggest that BraA04001030 is a strong candidate gene for fsm and provide the basis for exploring the molecular mechanism underlying female sterility in Chinese cabbage.

An APETALA1 ortholog affects plant architecture and seed yield component in oilseed rape (Brassica napus L.)

BACKGROUND

Increasing the productivity of rapeseed as one of the widely cultivated oil crops in the world is of upmost importance. As flowering time and plant architecture play a key role in the regulation of rapeseed yield, understanding the genetic mechanism underlying these traits can boost the rapeseed breeding. Meristem identity genes are known to have pleiotropic effects on plant architecture and seed yield in various crops. To understand the function of one of the meristem identity genes, APETALA1 (AP1) in rapeseed, we performed phenotypic analysis of TILLING mutants under greenhouse conditions. Three stop codon mutant families carrying a mutation in Bna.AP1.A02 paralog were analyzed for different plant architecture and seed yield-related traits.

RESULTS

It was evident that stop codon mutation in the K domain of Bna.AP1.A02 paralog caused significant changes in flower morphology as well as plant architecture related traits like plant height, branch height, and branch number. Furthermore, yield-related traits like seed yield per plant and number of seeds per plants were also significantly altered in the same mutant family. Apart from phenotypic changes, stop codon mutation in K domain of Bna.AP1.A02 paralog also altered the expression of putative downstream target genes like Bna.TFL1 and Bna.FUL in shoot apical meristem (SAM) of rapeseed. Mutant plants carrying stop codon mutations in the COOH domain of Bna.AP1.A02 paralog did not have a significant effect on plant architecture, yield-related traits or the expression of the downstream targets.

CONCLUSIONS

We found that Bna.AP1.A02 paralog has pleiotropic effect on plant architecture and yield-related traits in rapeseed. The allele we found in the current study with a beneficial effect on seed yield can be incorporated into rapeseed breeding pool to develop new varieties.

Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar

High-quality and dense genetic maps were constructed, and leaf shape variation was dissected by QTL mapping in poplar. Species in the genus Populus, also known as poplars, are important woody species and considered model plants for perennial trees. High-density genetic maps are valuable genomic resources for population genetics. Here, we generated a high-quality and dense genetic map for an F1 poplar population using high-throughput NGS-based genotyping. A total of 92,097 high-quality SNP markers were developed by stringent filtering and identification. In total, 889 and 1650 SNPs formed the female and male genetic maps, respectively. To test the application of the genetic maps, QTL mapping of leaf shape was conducted for this F1 population. A total of nine parameters were scored for leaf shape variation in three different environments. Combining genetic maps and measurements of the nine leaf shape parameters, we mapped a total of 42 significant QTLs. The highest LOD score of all QTLs was 9.2, and that QTL explained the most (15.13%) trait variation. A total of nine QTLs could be detected in at least two environments, and they were located in two genomic regions. Within these two QTL regions, some candidate genes for regulating leaf shape were predicted through functional annotation. The successful mapping of leaf shape QTLs demonstrated the utility of our genetic maps. According to the performance of this study, we were able to provide high-quality and dense genetic maps and dissect the leaf shape variation in poplar.

Integrated QTL and eQTL Mapping Provides Insights and Candidate Genes for Fatty Acid Composition, Flowering Time, and Growth Traits in a F2 Population of a Novel Synthetic Allopolyploid Brassica napus

Brassica napus (B. napus, AACC), is an economically important allotetraploid crop species that resulted from hybridization between two diploid species, Brassica rapa (AA) and Brassica olereacea (CC). We have created one new synthetic B. napus genotype Da-Ae (AACC) and one introgression line Da-Ol-1 (AACC), which were used to generate an F₂ mapping population. Plants in this F₂ mapping population varied in fatty acid content, flowering time, and growth-related traits. Using quantitative trait locus (QTL) mapping, we aimed to determine if Da-Ae and Da-Ol-1 provided novel genetic variation beyond what has already been found in B. napus. Making use of the genotyping information generated from RNA-seq data of these two lines and their F₂ mapping population of 166 plants, we constructed a genetic map consisting of 2,021 single nucleotide polymorphism markers that spans 2,929 cM across 19 linkage groups. Besides the known major QTL identified, our high resolution genetic map facilitated the identification of several new QTL contributing to the different fatty acid levels, flowering time, and growth-related trait values. These new QTL probably represent novel genetic variation that existed in our new synthetic B. napus strain. By conducting genome-wide expression variation analysis in our F₂ mapping population, genetic regions that potentially regulate many genes across the genome were revealed. A FLOWERING LOCUS C gene homolog, which was identified as a candidate regulating flowering time and multiple growth-related traits, was found underlying one of these regions. Integrated QTL and expression QTL analyses also helped us identified candidate causative genes associated with various biological traits through expression level change and/or possible protein function modification.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.

Translational Research: Exploring and Creating Genetic Diversity

The crop selection process has created a genetic bottleneck ultimately restricting breeding output. Wild relatives of major crops as well as the so-called 'neglected plant' species represent a reservoir of genetic diversity that remains underutilized. These species could be used as a tool to discover new alleles of agronomic interest or could be the target of breeding programs. Targeted induced local lesions in the genome (TILLING) can be used to translate in neglected crops what has been discovered in major crops and reciprocally. However, random mutagenesis, used in TILLING approaches, provides only a limited density of mutational events at a defined target locus. Alternatively, clustered regularly interspaced short palindromic repeats (CRISPR) associated 9 (Cas9) fused to a cytidine deaminase could serve as a localized mutagenic agent to produce high-density mutant populations. Artificial evolution is at hand.

Phenomic and Genomic Characterization of a Mutant Platform in Cucurbita pepo

The Cucurbita pepo genome comprises 263 Mb and 34,240 gene models organized in 20 different chromosomes. To improve our understanding of gene function we have generated an EMS mutant platform, consisting of 3,751 independent M2 families. The quality of the collection has been evaluated based on phenotyping and whole-genome re-sequencing (WGS) results. The phenotypic evaluation of the whole platform at seedling stage has demonstrated that the rate of variation for easily observable traits is more than 10%. The percentage of families with albino or chlorotic seedlings exceeded 3%, similar or higher to that found in other EMS collections of cucurbit crops. A rapid screening of the library for triple ethylene response in etiolated seedlings allowed the identification of four ethylene-insensitive mutants, that were found to be semidominant (ein1, ein2, and ein3) or dominant (EIN4). By evaluating 4 adult plants from 300 independent families more than 28% of apparent mutations were found for vegetative and reproductive traits, including plant vigor, leaf size and shape, sex expression and sex determination, and fruit set and development. Two pools of genomic DNA derived from 20 plants of two mutant families were subjected to WGS by using NGS methodology, estimating the density, spectrum, distribution and impact of EMS induced mutation. The number of EMS mutations in the genomes of families L1 and L2 was 1,704 and 859, respectively, which represents a density of 11.8 and 6 mutations per Mb, respectively. As expected, the predominant EMS induced mutations were C > T and G > A transitions (80.3% in L1, and 61% L2), that were found to be randomly distributed along the 20 chromosomes of C. pepo. The mutations were mostly affecting intergenic regions, but 7.9 and 6% of the identified EMS mutations in L1 and L2, respectively, were located in the exome, and 0.4 and 0.2% had a moderate and high putative impact on gene functions. These results provide information regarding the potential use of the obtained mutant platform in the discovery of novel alleles for both functional genomics and Cucurbita breeding by using direct- or reverse-genetic approaches.

Identification of Light-Independent Anthocyanin Biosynthesis Mutants Induced by Ethyl Methane Sulfonate in Turnip "Tsuda" (Brassica rapa)

The epidermis of swollen storage roots in purple cultivars of turnip “Tsuda” (Brassica rapa) accumulates anthocyanin in a light-dependent manner, especially in response to UV-A light, of which the mechanism is unclear. In this study, we mutagenized 15,000 seeds by 0.5% (v/v) ethyl methane sulfonate (EMS) and obtained 14 mutants with abnormal anthocyanin production in their epidermis of swollen storage roots. These mutants were classified into two groups: the red mutants with constitutive anthocyanin accumulation in their epidermis of storage roots even in underground parts in darkness and the white mutants without anthocyanin accumulation in the epidermis of storage roots in aboveground parts exposed to sunlight. Test cross analysis demonstrated that w9, w68, w204, r15, r21, r30 and r57 contained different mutations responsible for their phenotypic variations. Further genetic analysis of four target mutants (w9, w68, w204 and r15) indicated that each of them was controlled by a different recessive gene. Intriguingly, the expression profiles of anthocyanin biosynthesis genes, including structural and regulatory genes, coincided with their anthocyanin levels in the epidermis of storage roots in the four target mutants. We proposed that potential genes responsible for the mutations should be upstream factors of the anthocyanin biosynthesis pathway in turnips, which provided resources to further investigate the mechanisms of light-induced anthocyanin accumulation.

Breeding histories and selection criteria for oilseed rape in Europe and China identified by genome wide pedigree dissection

Selection breeding has played a key role in the improvement of seed yield and quality in oilseed rape (Brassica napus L.). We genotyped Tapidor (European), Ningyou7 (Chinese) and their progenitors with the Brassica 60 K Illumina Infinium SNP array and mapped a total of 29,347 SNP markers onto the reference genome of Darmor-bzh. Identity by descent (IBD) refers to a haplotype segment of a chromosome inherited from a shared common ancestor. IBDs identified on the C subgenome were larger than those on the A subgenome within both the Tapidor and Ningyou7 pedigrees. IBD number and length were greater in the Ningyou7 pedigree than in the Tapidor pedigree. Seventy nine QTLs for flowering time, seed quality and root morphology traits were identified in the IBDs of Tapidor and Ningyou7. Many more candidate genes had been selected within the Ningyou7 pedigree than within the Tapidor pedigree. These results highlight differences in the transfer of favorable gene clusters controlling key traits during selection breeding in Europe and China.

Genome-wide association mapping and Identification of candidate genes for fatty acid composition in Brassica napus L. using SNP markers

BACKGROUND

B. napus (oilseed) is an important source of edible vegetable oil, and its nutritional and economic value is determined by its fatty acid composition and content.

RESULTS

Using the Brassica 60 K SNP array, we performed a genome-wide association study of fatty acid composition in a population of 520 genetically diverse oilseed accessions. Using the PCA + K model in TASSEL 5.2.1, we identified 62 genomic regions that were significantly associated with the composition of seven fatty acids, and five consensus regions that mapped to the A2, A8, A9, C1, and C3 chromosomes, respectively, of the Brassica napus Darmor-bzh genome. We then identified 24 orthologs of the functional candidate genes involved in fatty acid biosynthesis, excluding BnaA.FAE1 and BnaC.FAE1 on the A8 and C3 homologous genome blocks, which are known to have critical roles in the fatty acid biosynthesis pathway, and potential orthologs of these genes (e.g., LACS9, KCR1, FAB1, LPAT4, KCS17, CER4, TT16, and ACBP5).

CONCLUSIONS

Our results demonstrate the power of association mapping in identifying genes of interest in B. napus and provide insight into the genetic basis of fatty acid biosynthesis in B. napus. Furthermore, our findings may facilitate marker-based breeding efforts aimed at improving fatty acid composition and quality in B. napus.

Progress in Genome Editing Technology and Its Application in Plants

Genome editing technology (GET) is a versatile approach that has progressed rapidly as a mechanism to alter the genotype and phenotype of organisms. However, conventional genome modification using GET cannot satisfy current demand for high-efficiency and site-directed mutagenesis, retrofitting of artificial nucleases has developed into a new avenue within this field. Based on mechanisms to recognize target genes, newly-developed GETs can generally be subdivided into three cleavage systems, protein-dependent DNA cleavage systems (i.e., zinc-finger nucleases, ZFN, and transcription activator-like effector nucleases, TALEN), RNA-dependent DNA cleavage systems (i.e., clustered regularly interspaced short palindromic repeats-CRISPR associated proteins, CRISPR-Cas9, CRISPR-Cpf1, and CRISPR-C2c1), and RNA-dependent RNA cleavage systems (i.e., RNA interference, RNAi, and CRISPR-C2c2). All these techniques can lead to double-stranded (DSB) or single-stranded breaks (SSB), and result in either random mutations via non-homologous end-joining (NHEJ) or targeted mutation via homologous recombination (HR). Thus, site-directed mutagenesis can be induced via targeted gene knock-out, knock-in, or replacement to modify specific characteristics including morphology-modification, resistance-enhancement, and physiological mechanism-improvement along with plant growth and development. In this paper, an non-comprehensive review on the development of different GETs as applied to plants is presented.

Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus

Seed oil content is an important agricultural trait in rapeseed breeding. Although numerous quantitative trait locus (QTL) have been identified, most of them cannot be applied in practical breeding mainly due to environmental instability or large confidence intervals. The purpose of this study was to identify and validate high quality and more stable QTLs by combining linkage mapping and genome-wide association study (GWAS). For linkage mapping, we constructed two F₂ populations from crosses of high-oil content (∼50%) lines 6F313 and 61616 with a low-oil content (∼40%) line 51070. Two high density linkage maps spanned 1987cM (1659 bins) and 1856cM (1746 bins), respectively. For GWAS, we developed more than 34,000 high-quality SNP markers based on 227 accessions. Finally, 40 QTLs and 29 associations were established by linkage and association mapping in different environments. After merging the results, 32 consensus QTLs were obtained and 7 of them were identified by both mapping methods. Seven overlapping QTLs covered an average confidence interval of 183kb and explained the phenotypic variation of 10.23 to 24.45%. We further developed allele-specific PCR primers to identify each of the seven QTLs. These stable QTLs should be useful in gene cloning and practical breeding application.

Characteristics Analysis of F1 Hybrids between Genetically Modified Brassica napus and B. rapa

A number of studies have been conducted on hybridization between transgenic Brassica napus and B. rapa or backcross of F1 hybrid to their parents. However, trait changes must be analyzed to evaluate hybrid sustainability in nature. In the present study, B. rapa and transgenic (BrAGL20) B. napus were hybridized to verify the early flowering phenomenon of F1 hybrids, and F1 hybrid traits were analyzed to predict their impact on sustainability. Flowering of F1 hybrid has been induced slightly later than that of the transgenic B. napus, but flowering was available in the greenhouse without low temperature treatment to young plant, similar to the transgenic B. napus. It is because the BrAGL20 gene has been transferred from transgenic B. napus to F1 hybrid. The size of F1 hybrid seeds was intermediate between those of B. rapa and transgenic B. napus, and ~40% of F1 pollen exhibited abnormal size and morphology. The form of the F1 stomata was also intermediate between that of B. rapa and transgenic B. napus, and the number of stomata was close to the parental mean. Among various fatty acids, the content of erucic acid exhibited the greatest change, owing to the polymorphism of parental FATTY ACID ELONGASE 1 alleles. Furthermore, F2 hybrids could not be obtained. However, BC1 progeny were obtained by hand pollination of B. rapa with F1 hybrid pollen, with an outcrossing rate of 50%.

Microspore Induced Doubled Haploids Production from Ethyl Methanesulfonate (EMS) Soaked Flower Buds Is an Efficient Strategy for Mutagenesis in Chinese Cabbage

Chinese cabbage buds were soaked with Ethyl methanesulfonate (EMS) to induce mutagenesis. The influence of different EMS concentrations and treatment durations on microspore development, embryo production rate and seedling rate were evaluated in five Chinese cabbage genotypes. Mutations in four color-related genes were identified using high resolution melting (HRM) curves of their PCR products. The greatest embryo production and seedling rates were observed when buds were treated with 0.03 to 0.1% EMS for 5 to 10 min, while EMS concentrations greater than 0.1% were lethal to the microspores. In total, 142 mutants with distinct variations in leaf shape, leaf color, corolla size, flower color, bolting time and downy mildew resistance were identified from 475 microspore culture derived Doubled Haploids. Our results demonstrate that microspore derived Doubled Haploids from EMS soaked buds represents an efficient approach to rapidly generate homozygous Chinese cabbage mutants.

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.

TILLING by sequencing to identify induced mutations in stress resistance genes of peanut (Arachis hypogaea)

Background

Targeting Induced Local Lesions in Genomes (TILLING) is a powerful reverse genetics approach for functional genomics studies. We used high-throughput sequencing, combined with a two-dimensional pooling strategy, with either minimum read percentage with non-reference nucleotide or minimum variance multiplier as mutation prediction parameters, to detect genes related to abiotic and biotic stress resistances. In peanut, lipoxygenase genes were reported to be highly induced in mature seeds infected with Aspergillus spp., indicating their importance in plant-fungus interactions. Recent studies showed that phospholipase D (PLD) expression was elevated more quickly in drought sensitive lines than in drought tolerant lines of peanut. A newly discovered lipoxygenase (LOX) gene in peanut, along with two peanut PLD genes from previous publications were selected for TILLING. Additionally, two major allergen genes Ara h 1 and Ara h 2, and fatty acid desaturase AhFAD2, a gene which controls the ratio of oleic to linoleic acid in the seed, were also used in our study. The objectives of this research were to develop a suitable TILLING by sequencing method for this allotetraploid, and use this method to identify mutations induced in stress related genes.

Results

We screened a peanut root cDNA library and identified three candidate LOX genes. The gene AhLOX7 was selected for TILLING due to its high expression in seeds and roots. By screening 768 M2 lines from the TILLING population, four missense mutations were identified for AhLOX7, three missense mutations were identified for AhPLD, one missense and two silent mutations were identified for Ara h 1.01, three silent and five missense mutations were identified for Ara h 1.02, one missense mutation was identified for AhFAD2B, and one silent mutation was identified for Ara h 2.02. The overall mutation frequency was 1 SNP/1,066 kb. The SNP detection frequency for single copy genes was 1 SNP/344 kb and 1 SNP/3,028 kb for multiple copy genes.

Conclusions

Our TILLING by sequencing approach is efficient to identify mutations in single and multi-copy genes. The mutations identified in our study can be used to further study gene function and have potential usefulness in breeding programs.

Electronic supplementary material

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

New approaches to facilitate rapid domestication of a wild plant to an oilseed crop: example pennycress (Thlaspi arvense L.)

Oilseed crops are sources of oils and seed meal having a multitude of uses. While the domestication of soybean and rapeseed took extended periods of time, new genome-based techniques have ushered in an era where crop domestication can occur rapidly. One attractive target for rapid domestication is the winter annual plant Field Pennycress (Thlaspi arvense L.; pennycress; Brassicaceae). Pennycress grows widespread throughout temperate regions of the world and could serve as a winter oilseed-producing cover crop. If grown throughout the USA Midwest Corn Belt, for example, pennycress could produce as much as 840L/ha oils and 1470kg/ha press-cake annually on 16 million hectares of farmland currently left fallow during the fall through spring months. However, wild pennycress strains have inconsistent germination and stand establishment, un-optimized maturity for a given growth zone, suboptimal oils and meal quality for biofuels and food production, and significant harvest loss due to pod shatter. In this review, we describe the virtues and current shortcomings of pennycress and discuss how knowledge from studying Arabidopsis thaliana and other Brassicas, in combination with the advent of affordable next generation sequencing, can bring about the rapid domestication and improvement of pennycress and other crops.

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.

Genome-wide investigation of genetic changes during modern breeding of Brassica napus

Considerable genome variation had been incorporated within rapeseed breeding programs over past decades. In past decades, there have been substantial changes in phenotypic properties of rapeseed as a result of extensive breeding effort. Uncovering the underlying patterns of allelic variation in the context of genome organisation would provide knowledge to guide future genetic improvement. We assessed genome-wide genetic changes, including population structure, genetic relatedness, the extent of linkage disequilibrium, nucleotide diversity and genetic differentiation based on F ST outlier detection, for a panel of 472 Brassica napus inbred accessions using a 60 k Brassica Infinium® SNP array. We found genetic diversity varied in different sub-groups. Moreover, the genetic diversity increased from 1950 to 1980 and then remained at a similar level in China and Europe. We also found ~6-10 % genomic regions revealed high F ST values. Some QTLs previously associated with important agronomic traits overlapped with these regions. Overall, the B. napus C genome was found to have more high F ST signals than the A genome, and we concluded that the C genome may contribute more valuable alleles to generate elite traits. The results of this study indicate that considerable genome variation had been incorporated within rapeseed breeding programs over past decades. These results also contribute to understanding the impact of rapeseed improvement on available genome variation and the potential for dissecting complex agronomic traits.

Development of a Cucumis sativus TILLinG platform for forward and reverse genetics

Background

Cucumber (Cucumis sativus) belongs to the Cucurbitaceae family that includes more than 800 species. The cucumber genome has been recently sequenced and annotated. Transcriptomics and genome sequencing of many plant genomes are providing information on candidate genes potentially related to agronomically important traits. To accelerate functional characterization of these genes in cucumber we have generated an EMS mutant population that can be used as a TILLinG platform for reverse genetics.

Principal Findings

A population of 3,331 M2 mutant seed families was generated using two EMS concentrations (0.5% and 0.75%). Genomic DNA was extracted from M2 families and eight-fold pooled for mutation detection by ENDO1 nuclease. To assess the quality of the mutant collection, we screened for induced mutations in five genes and identified 26 mutations. The average mutation rate was calculated as 1/1147 Kb giving rise to approximately 320 mutations per genome. We focused our characterization on three missense mutations, G33C, S238F and S249F identified in the CsACS2 sex determination gene. Protein modeling and crystallography studies predicted that mutation at G³³ may affect the protein function, whereas mutations at S²³⁸ and S²⁴⁹ may not impair the protein function. As predicted, detailed phenotypic evaluation showed that the S238F and the S249F mutant lines had no sexual phenotype. In contrast, plants homozygous for the G33C mutation showed a complete sexual transition from monoecy to andromonoecy. This result demonstrates that TILLinG is a valuable tool for functional validation of gene function in crops recalcitrant to transgenic transformation.

Conclusions

We have developed a cucumber mutant population that can be used as an efficient reverse genetics tool. The cucumber TILLinG collection as well as the previously described melon TILLinG collection will prove to be a valuable resource for both fundamental research and the identification of agronomically-important genes for crop improvement in cucurbits in general.

Progress in TILLING as a tool for functional genomics and improvement of crops

Food security is a global concern and substantial yield increases in crops are required to feed the growing world population. Mutagenesis is an important tool in crop improvement and is free of the regulatory restrictions imposed on genetically modified organisms. Targeting Induced Local Lesions in Genomes (TILLING), which combines traditional chemical mutagenesis with high-throughput genome-wide screening for point mutations in desired genes, offers a powerful way to create novel mutant alleles for both functional genomics and improvement of crops. TILLING is generally applicable to genomes whether small or large, diploid or even allohexaploid, and shows great potential to address the major challenge of linking sequence information to the function of genes and to modulate key traits for plant breeding. TILLING has been successfully applied in many crop species and recent progress in TILLING is summarized below, especially on the developments in mutation detection technology, application of TILLING in gene functional studies and crop breeding. The potential of TILLING/EcoTILLING for functional genetics and crop improvement is also discussed. Furthermore, a small-scale forward strategy including backcross and selfing was conducted to release the potential mutant phenotypes masked in M2 (or M3) plants.

Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.)

This report describes an integrative analysis of seed-oil-content quantitative trait loci (QTL) in Brassica napus , using a high-density genetic map to align QTL among different populations. Rapeseed (Brassica napus) is an important source of edible oil and sustainable energy. Given the challenge involved in using only a few genes to substantially increase the oil content of rapeseed without affecting the fatty acid composition, exploitation of a greater number of genetic loci that regulate the oil content variation among rapeseed germplasm is of fundamental importance. In this study, we investigated variation in the seed-oil content among two related genetic populations of Brassica napus, the TN double-haploid population and its derivative reconstructed-F2 population. Each population was grown in multiple experiments under different environmental conditions. Mapping of quantitative trait loci (QTL) identified 41 QTL in the TN populations. Furthermore, of the 20 pairs of epistatic interaction loci detected, approximately one-third were located within the QTL intervals. The use of common markers on different genetic maps and the TN genetic map as a reference enabled us to project QTL from an additional three genetic populations onto the TN genetic map. In summary, we used the TN genetic map of the B. napus genome to identify 46 distinct QTL regions that control seed-oil content on 16 of the 19 linkage groups of B. napus. Of these, 18 were each detected in multiple populations. The present results are of value for ongoing efforts to breed rapeseed with high oil content, and alignment of the QTL makes an important contribution to the development of an integrative system for genetic studies of rapeseed.