Associative transcriptomics of traits in the polyploid crop species Brassica napus

Genomic and transcriptome analyses reveal potential contributors to erucic acid biosynthesis in seeds of rapeseed (Brassica napus)

KEY MESSAGE

Through comprehensive genomic and transcriptomic analyses, we identified a set of 23 genes that act up- or downstream of erucic acid content (EAC) production in rapeseed seeds. We selected example genes to showcase the distribution of single nucleotide polymorphisms, haplotypes associated with EAC phenotypes, and the creation of molecular markers differentiating low EAC and high EAC genotypes. Erucic acid content (EAC) is a crucial trait in rapeseed, with low LEAC oil recognized for its health benefits and high EA oil holding industrial value. Despite its significance, the genomic consequences of intensive LEAC-cultivar selection and the genetic basis underlying EA regulation remain largely unexplored. To address this knowledge gap, we conducted selective signal analyses, genome-wide association studies (GWAS), and transcriptome analyses. Our investigation unveiled the genetic footprints resulting from LEAC selection in germplasm populations, drawing attention to specific loci that contribute to enriching diversity. By integrating GWAS and transcriptome analyses, we identified a set of 23 genes that play a significant role in determining EAC in seeds or are downstream consequences of EA-level alterations. These genes have emerged as promising candidates for elucidating the potential mechanisms governing EAC in rapeseed. To exemplify the findings, we selected specific genes to demonstrate the distribution of single nucleotide polymorphisms and haplotypes associated with different EAC phenotypes. Additionally, we showcased to develop molecular markers distinguishing between LEAC and high EAC genotypes.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

Pathogen lifestyle determines host genetic signature of quantitative disease resistance loci in oilseed rape (Brassica napus)

KEY MESSAGE

Using associative transcriptomics, our study identifies genes conferring resistance to four diverse fungal pathogens in crops, emphasizing key genetic determinants of multi-pathogen resistance. Crops are affected by several pathogens, but these are rarely studied in parallel to identify common and unique genetic factors controlling diseases. Broad-spectrum quantitative disease resistance (QDR) is desirable for crop breeding as it confers resistance to several pathogen species. Here, we use associative transcriptomics (AT) to identify candidate gene loci associated with Brassica napus constitutive QDR to four contrasting fungal pathogens: Alternaria brassicicola, Botrytis cinerea, Pyrenopeziza brassicae, and Verticillium longisporum. We did not identify any shared loci associated with broad-spectrum QDR to fungal pathogens with contrasting lifestyles. Instead, we observed QDR dependent on the lifestyle of the pathogen-hemibiotrophic and necrotrophic pathogens had distinct QDR responses and associated loci, including some loci associated with early immunity. Furthermore, we identify a genomic deletion associated with resistance to V. longisporum and potentially broad-spectrum QDR. This is the first time AT has been used for several pathosystems simultaneously to identify host genetic loci involved in broad-spectrum QDR. We highlight constitutive expressed candidate loci for broad-spectrum QDR with no antagonistic effects on susceptibility to the other pathogens studies as candidates for crop breeding. In conclusion, this study represents an advancement in our understanding of broad-spectrum QDR in B. napus and is a significant resource for the scientific community.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Toxic effects of lead on plants: integrating multi-omics with bioinformatics to develop Pb-tolerant crops

MAIN CONCLUSION

Lead disrupts plant metabolic homeostasis and key structural elements. Utilizing modern biotechnology tools, it's feasible to develop Pb-tolerant varieties by discovering biological players regulating plant metabolic pathways under stress. Lead (Pb) has been used for a variety of purposes since antiquity despite its toxic nature. After arsenic, lead is the most hazardous heavy metal without any known beneficial role in the biological system. It is a crucial inorganic pollutant that affects plant biochemical and morpho-physiological attributes. Lead toxicity harms plants throughout their life cycle and the extent of damage depends on the concentration and duration of exposure. Higher levels of lead exposure disrupt numerous key metabolic activities of plants including oxygen-evolving complex, organelles integrity, photosystem II connectivity, and electron transport chain. This review summarizes the detrimental effects of lead toxicity on seed germination, crop growth, and yield, oxidative and ultra-structural alterations, as well as nutrient absorption, transport, and assimilation. Further, it discusses the Pb-induced toxic modulation of stomatal conductance, photosynthesis, respiration, metabolic-enzymatic activity, osmolytes accumulation, and antioxidant activity. It is a comprehensive review that reports on omics-based studies along with morpho-physiological and biochemical modifications caused by lead stress. With advances in DNA sequencing technologies, genomics and transcriptomics are gradually becoming popular for studying Pb stress effects in plants. Proteomics and metabolomics are still underrated and there is a scarcity of published data, and this review highlights both their technical and research gaps. Besides, there is also a discussion on how the integration of omics with bioinformatics and the use of the latest biotechnological tools can aid in developing Pb-tolerant crops. The review concludes with core challenges and research directions that need to be addressed soon.

Transcriptome-wide association study coupled with eQTL analysis reveals the genetic connection between gene expression and flowering time in Arabidopsis

Genome-wide association study (GWAS) has improved our understanding of complex traits, but challenges exist in distinguishing causation versus association caused by linkage disequilibrium. Instead, transcriptome-wide association studies (TWAS) detect direct associations between expression levels and phenotypic variations, providing an opportunity to better prioritize candidate genes. To assess the feasibility of TWAS, we investigated the association between transcriptomes, genomes, and various traits in Arabidopsis, including flowering time. The associated genes formerly known to regulate growth allometry or metabolite production were first identified by TWAS. Next, for flowering time, six TWAS-newly identified genes were functionally validated. Analysis of the expression quantitative trait locus (eQTL) further revealed a trans-regulatory hotspot affecting the expression of several TWAS-identified genes. The hotspot covers the FRIGIDA (FRI) gene body, which possesses multiple haplotypes differentially affecting the expression of downstream genes, such as FLOWERING LOCUS C (FLC) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1). We also revealed multiple independent paths towards the loss of function of FRI in natural accessions. Altogether, this study demonstrates the potential of combining TWAS with eQTL analysis to identify important regulatory modules of FRI-FLC-SOC1 for quantitative traits in natural populations.

A causal inference and Bayesian optimisation framework for modelling multi-trait relationships-Proof-of-concept using Brassica napus seed yield under controlled conditions

The improvement of crop yield is a major breeding target and there is a long history of research that has focussed on unravelling the mechanisms and processes that contribute to yield. Quantitative prediction of the interplay between morphological traits, and the effects of these trait-trait relationships on seed production remains, however, a challenge. Consequently, the extent to which crop varieties optimise their morphology for a given environment is largely unknown. This work presents a new combination of existing methodologies by framing crop breeding as an optimisation problem and evaluates the extent to which existing varieties exhibit optimal morphologies under the test conditions. In this proof-of-concept study using spring and winter oilseed rape plants grown under greenhouse conditions, we employ causal inference to model the hierarchically structured effects of 27 morphological yield traits on each other. We perform Bayesian optimisation of seed yield, to identify and quantify the morphologies of ideotype plants, which are expected to be higher yielding than the varieties in the studied panels. Under the tested growth conditions, we find that existing spring varieties occupy the optimal regions of trait-space, but that potentially high yielding strategies are unexplored in extant winter varieties. The same approach can be used to evaluate trait (morphology) space for any environment.

Exome-based new allele-specific PCR markers and transferability for sodicity tolerance in bread wheat (Triticum aestivum L.)

Targeted exome-based genotype by sequencing (t-GBS), a sequencing technology that tags SNPs and haplotypes in gene-rich regions was used in previous genome-wide association studies (GWAS) for sodicity tolerance in bread wheat. Thirty-nine novel SNPs including 18 haplotypes for yield and yield-components were identified. The present study aimed at developing SNP-derived markers by precisely locating new SNPs on ~180 bp allelic sequence of t-GBS, marker validation, and SNP functional characterization based on its exonic location. We identified unknown locations of significant SNPs/haplotypes by aligning allelic sequences on to IWGSC RefSeqv1.0 on respective chromosomes. Eighteen out of the target 39 SNP locations fulfilled the criteria for producing PCR markers, among which only eight produced polymorphic signals. These eight markers associated with yield, plants m-2, heads m-2, and harvest index, including a pleiotropic marker for yield, harvest index, and grains/head were validated for its amplification efficiency and phenotypic effects in focused identification germplasm strategy (FIGS) wheat set and a doubled haploid (DH) population (Scepter/IG107116). The phenotypic variation explained by these markers are in the range of 4.1-37.6 in the FIGS population. High throughput PCR-based genotyping using new markers and association with phenotypes in FIGS wheat set and DH population validated the effect of functional SNP on closely associated genes-calcineurin B-like- and dirigent protein, basic helix-loop-helix (BHLH-), plant homeodomain (PHD-) and helix-turn-helix myeloblastosis (HTH myb) type -transcription factor. Further, genome-wide SNP annotation using SnpEff tool confirmed that these SNPs are in gene regulatory regions (upstream, 3'-UTR, and intron) modifying gene expression and protein-coding. This integrated approach of marker design for t-GBS alleles, SNP functional annotation, and high-throughput genotyping of functional SNP offers translation solutions across crops and complex traits in crop improvement programs.

Characterization of novel loci controlling seed oil content in Brassica napus by marker metabolite-based multi-omics analysis

BACKGROUND

Seed oil content is an important agronomic trait of Brassica napus (B. napus), and metabolites are considered as the bridge between genotype and phenotype for physical traits.

RESULTS

Using a widely targeted metabolomics analysis in a natural population of 388 B. napus inbred lines, we quantify 2172 metabolites in mature seeds by liquid chromatography mass spectrometry, in which 131 marker metabolites are identified to be correlated with seed oil content. These metabolites are then selected for further metabolite genome-wide association study and metabolite transcriptome-wide association study. Combined with weighted correlation network analysis, we construct a triple relationship network, which includes 21,000 edges and 4384 nodes among metabolites, metabolite quantitative trait loci, genes, and co-expression modules. We validate the function of BnaA03.TT4, BnaC02.TT4, and BnaC05.UK, three candidate genes predicted by multi-omics analysis, which show significant impacts on seed oil content through regulating flavonoid metabolism in B. napus.

CONCLUSIONS

This study demonstrates the advantage of utilizing marker metabolites integrated with multi-omics analysis to dissect the genetic basis of agronomic traits in crops.

BnIR: A multi-omics database with various tools for Brassica napus research and breeding

In the post-genome-wide association study era, multi-omics techniques have shown great power and potential for candidate gene mining and functional genomics research. However, due to the lack of effective data integration and multi-omics analysis platforms, such techniques have not still been applied widely in rapeseed, an important oil crop worldwide. Here, we report a rapeseed multi-omics database (BnIR; http://yanglab.hzau.edu.cn/BnIR), which provides datasets of six omics including genomics, transcriptomics, variomics, epigenetics, phenomics, and metabolomics, as well as numerous "variation-gene expression-phenotype" associations by using multiple statistical methods. In addition, a series of multi-omics search and analysis tools are integrated to facilitate the browsing and application of these datasets. BnIR is the most comprehensive multi-omics database for rapeseed so far, and two case studies demonstrated its power to mine candidate genes associated with specific traits and analyze their potential regulatory mechanisms.

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

Quantitative disease resistance (QDR) controls the association of the light leaf spot pathogen with Brassica napus; four QDR loci that were in linkage disequilibrium and eight gene expression markers were identified. Quantitative disease resistance (QDR) can provide durable control of pathogens in crops in contrast to resistance (R) gene-mediated resistance which can break down due to pathogen evolution. QDR is therefore a desirable trait in crop improvement, but little is known about the causative genes, and so it is difficult to incorporate into breeding programmes. Light leaf spot, caused by Pyrenopeziza brassicae, is an important disease of oilseed rape (canola, Brassica napus). To identify new QDR gene loci, we used a high-throughput screening pathosystem with P. brassicae on 195 lines of B. napus combined with an association transcriptomics platform. We show that all resistance against P. brassicae was associated with QDR and not R gene-mediated. We used genome-wide association analysis with an improved B. napus population structure to reveal four gene loci significantly (P = 0.0001) associated with QDR in regions showing linkage disequilibrium. On chromosome A09, enhanced resistance was associated with heterozygosity for a cytochrome P450 gene co-localising with a previously described locus for seed glucosinolate content. In addition, eight significant gene expression markers with a false discovery rate of 0.001 were associated with QDR against P. brassicae. For seven of these, expression was positively correlated with resistance, whereas for one, a HXXXD-type acyl-transferase, negative correlation indicated a potential susceptibility gene. The study identifies novel QDR loci for susceptibility and resistance, including novel cryptic QDR genes associated with heterozygosity, that will inform future crop improvement.

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.

High-Throughput Association Mapping in Brassica napus L.: Methods and Applications

Oil seed rape (Braasica napus L.) is ranked second among oil seed crops cultivated globally for edible oil for human, and seed cake for animal consumption. Recent genetic and genomics advancements highlighted the diversity that exists within B. napus, which is largely discovered using the most promising genetic markers called single nucleotide polymorphism (SNP). Their calling rate is also enhanced to ~100 folds after the continuous advancements in the next generation sequencing (NGS) technologies. As the high throughput of NGS resulted in multi-Giga bases data, the detailed quality control (QC) prior to downstream analyses is a pre-requisite. It mainly involved the removal of false positives, missing proportions, filtering of low-quality SNPs, and adjustments of minor-allele frequency and heterozygosity. After marker-trait association, for conformation of target SNPs, validations of SNPs can be performed using various methods, especially allele-specific PCR assay-based methods have been utilized for SNP genotyping of genes targeting agronomic traits and somaclonal variations occurred during transgenic studies. In the present study, the authors mainly argue on the genotypic progress, and pipelines/methods that are being used for detection, calling, filtering, and validation of SNPs. Also, insight is provided into the application of SNPs in linkage and association mapping, including QTL mapping and genome-wide association studies targeting mainly developmental traits related to the root system and plant architecture, flowering time, silique, and oil quality. Briefly, the present study provides the recent information and recommendations on the SNP genotyping methods and its applications, which can be useful for marker-assisted breeding in B. napus and other crops.

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.

Natural Variation in OASC Gene for Mitochondrial O-Acetylserine Thiollyase Affects Sulfate Levels in Arabidopsis

Sulfur plays a vital role in the primary and secondary metabolism of plants, and carries an important function in a large number of different compounds. Despite this importance, compared to other mineral nutrients, relatively little is known about sulfur sensing and signalling, as well as about the mechanisms controlling sulfur metabolism and homeostasis. Sulfur contents in plants vary largely not only among different species, but also among accessions of the same species. We previously used associative transcriptomics to identify several genes potentially controlling variation in sulfate content in the leaves of Brassica napus, including an OASC gene for mitochondrial O-acetylserine thiollyase (OAS-TL), an enzyme involved in cysteine synthesis. Here, we show that loss of OASC in Arabidopsis thaliana lowers not only sulfate, but also glutathione levels in the leaves. The reduced accumulation is caused by lower sulfate uptake and translocation to the shoots; however, the flux through the pathway is not affected. In addition, we identified a single nucleotide polymorphism in the OASC gene among A. thaliana accessions that is linked to variation in sulfate content. Both genetic and transgenic complementation confirmed that the exchange of arginine at position 81 for lysine in numerous accessions resulted in a less active OASC and a lower sulfate content in the leaves. The mitochondrial isoform of OAS-TL is, thus, after the ATPS1 isoform of sulfurylase and the APR2 form of APS reductase 2, the next metabolic enzyme with a role in regulation of sulfate content in Arabidopsis.

CRISPR/Cas9-Mediated Targeted Mutagenesis of BnaCOL9 Advances the Flowering Time of Brassica napus L

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. The planting area and output of rapeseed are affected by the flowering time, which is a critical agronomic feature. COL9 controls growth and development in many different plant species as a member of the zinc finger transcription factor family. However, BnaCOL9 in rapeseed has not been documented. The aim of this study was to apply CRISPR/Cas9 technology to create an early-flowering germplasm resource to provide useful material for improving the early-maturing breeding of rapeseed. We identified four COL9 homologs in rapeseed that were distributed on chromosomes A05, C05, A03, and C03. We successfully created quadruple BnaCOL9 mutations in rapeseed using the CRISPR/Cas9 platform. The quadruple mutants of BnaCOL9 flowered earlier than the wild-type. On the other hand, the flowering time of the BnaCOL9 overexpression lines was delayed. An analysis of the expression patterns revealed that these genes were substantially expressed in the leaves and flowers. A subcellular localization experiment demonstrated that BnaCOL9 was in the nucleus. Furthermore, we discovered that two key flowering-related genes, BnaCO and BnaFT, were highly elevated in the BnaCOL9 mutants, but dramatically downregulated in the BnaCOL9 overexpression lines. Our findings demonstrate that BnaCOL9 is a significant flowering inhibitor in rapeseed and may be employed as a crucial gene for early-maturing breeding.

Winter warming post floral initiation delays flowering via bud dormancy activation and affects yield in a winter annual crop

In temperate climates many plant species use long-term detection of winter chilling as a seasonal cue. Previously the timing of flowering in winter annual plants has been shown to be controlled by the promotion of the floral transition by chilling, known as vernalization. In contrast, many temperate perennial species produce flower buds prior to winter and require winter chilling to break bud dormancy to enable bud break and flowering in the following spring. Here we show that flowering time in winter annuals can be controlled by bud dormancy and that in winter oilseed rape–reduced chilling during flower bud dormancy is associated with yield declines.

Winter annual life history is conferred by the requirement for vernalization to promote the floral transition and control the timing of flowering. Here we show using winter oilseed rape that flowering time is controlled by inflorescence bud dormancy in addition to vernalization. Winter warming treatments given to plants in the laboratory and field increase flower bud abscisic acid levels and delay flowering in spring. We show that the promotive effect of chilling reproductive tissues on flowering time is associated with the activity of two FLC genes specifically silenced in response to winter temperatures in developing inflorescences, coupled with activation of a BRANCHED1-dependent bud dormancy transcriptional module. We show that adequate winter chilling is required for normal inflorescence development and high yields in addition to the control of flowering time. Because warming during winter flower development is associated with yield losses at the landscape scale, our work suggests that bud dormancy activation may be important for effects of climate change on winter arable crop yields.

Further insight into decreases in seed glucosinolate content based on QTL mapping and RNA-seq in Brassica napus L

The QTL hotspots determining seed glucosinolate content instead of only four HAG1 loci and elucidation of a potential regulatory model for rapeseed SGC variation. Glucosinolates (GSLs) are amino acid-derived, sulfur-rich secondary metabolites that function as biopesticides and flavor compounds, but the high seed glucosinolate content (SGC) reduces seed quality for rapeseed meal. To dissect the genetic mechanism and further reduce SGC in rapeseed, QTL mapping was performed using an updated high-density genetic map based on a doubled haploid (DH) population derived from two parents that showed significant differences in SGC. In 15 environments, a total of 162 significant QTLs were identified for SGC and then integrated into 59 consensus QTLs, of which 32 were novel QTLs. Four QTL hotspot regions (QTL-HRs) for SGC variation were discovered on chromosomes A09, C02, C07 and C09, including seven major QTLs that have previously been reported and four novel major QTLs in addition to HAG1 loci. SGC was largely determined by superimposition of advantage allele in the four QTL-HRs. Important candidate genes directly related to GSL pathways were identified underlying the four QTL-HRs, including BnaC09.MYB28, BnaA09.APK1, BnaC09.SUR1 and BnaC02.GTR2a. Related differentially expressed candidates identified in the minor but environment stable QTLs indicated that sulfur assimilation plays an important rather than dominant role in SGC variation. A potential regulatory model for rapeseed SGC variation constructed by combining candidate GSL gene identification and differentially expressed gene analysis based on RNA-seq contributed to a better understanding of the GSL accumulation mechanism. This study provides insights to further understand the genetic regulatory mechanism of GSLs, as well as the potential loci and a new route to further diminish the SGC in rapeseed.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Transcriptomic analysis of rapeseed (Brassica napus. L.) seed development in Xiangride, Qinghai Plateau, reveals how its special eco-environment results in high yield in high-altitude areas

As one of the most important oil crops, rapeseed (Brassica napus) is cultivated worldwide to produce vegetable oil, animal feed, and biodiesel. As the population grows and the need for renewable energy increases, the breeding and cultivation of high-yield rapeseed varieties have become top priorities. The formation of a high rapeseed yield is so complex because it is influenced not only by genetic mechanisms but also by many environmental conditions, such as climatic conditions and different farming practices. Interestingly, many high-yield areas are located in special eco-environments, for example, in the high-altitude Xiangride area of the Qinghai Plateau. However, the molecular mechanisms underlying the formation of high yields in such a special eco-environment area remain largely unknown. Here, we conducted field yield analysis and transcriptome analysis in the Xiangride area. Compared with the yield and environmental factors in the Xinning area (a low-yielding area), we found that the relatively longer daylight length is the key to high rapeseed yield in the Xiangride area, which leads up to a 52.1% increase in rapeseed yield, especially the increase in thousand seed weight and silique number (SN). Combined with transcriptome H-cluster analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses, we can assume that the grain development of rapeseed in the Xiangride area is ahead of schedule and lasts for a long time, leading to the high-yield results in the Xiangride area, confirmed by the expression analysis by quantitative real-time polymerase chain reaction (qRT-PCR) of yield-related genes. Our results provide valuable information for further exploring the molecular mechanism underlying high yield in special ecological environments and provide a helpful reference for studying seed development characteristics in special-producing regions for Brassica napus.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus

Background

Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus.

Results

Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus.

Conclusions

This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02647-5.

Identification of Chilling Accumulation-Associated Genes for Litchi Flowering by Transcriptome-Based Genome-Wide Association Studies

Litchi is an important Sapindaceae fruit tree. Flowering in litchi is triggered by low temperatures in autumn and winter. It can be divided into early-, medium-, and late-flowering phenotypes according to the time for floral induction. Early-flowering varieties need low chilling accumulation level for floral induction, whereas the late-flowering varieties require high chilling accumulation level. In the present study, RNA-Seq of 87 accessions was performed and transcriptome-based genome-wide association studies (GWAS) was used to identify candidate genes involved in chilling accumulation underlying the time for floral induction. A total of 98,155 high-quality single-nucleotide polymorphism (SNP) sites were obtained. A total of 1,411 significantly associated SNPs and 1,115 associated genes (AGs) were identified, of which 31 were flowering-related, 23 were hormone synthesis-related, and 27 were hormone signal transduction-related. Association analysis between the gene expression of the AGs and the flowering phenotypic data was carried out, and differentially expressed genes (DEGs) in a temperature-controlled experiment were obtained. As a result, 15 flowering-related candidate AGs (CAGs), 13 hormone synthesis-related CAGs, and 11 hormone signal transduction-related CAGs were further screened. The expression levels of the CAGs in the early-flowering accessions were different from those in the late-flowering ones, and also between the flowering trees and non-flowering trees. In a gradient chilling treatment, flowering rates of the trees and the CAGs expression were affected by the treatment. Our present work for the first time provided candidate genes for genetic regulation of flowering in litchi using transcriptome-based GWAS.

A transcriptome-based association study of growth, wood quality, and oleoresin traits in a slash pine breeding population

Slash pine (Pinus elliottii Engelm.) is an important timber and resin species in the United States, China, Brazil and other countries. Understanding the genetic basis of these traits will accelerate its breeding progress. We carried out a genome-wide association study (GWAS), transcriptome-wide association study (TWAS) and weighted gene co-expression network analysis (WGCNA) for growth, wood quality, and oleoresin traits using 240 unrelated individuals from a Chinese slash pine breeding population. We developed high quality 53,229 single nucleotide polymorphisms (SNPs). Our analysis reveals three main results: (1) the Chinese breeding population can be divided into three genetic groups with a mean inbreeding coefficient of 0.137; (2) 32 SNPs significantly were associated with growth and oleoresin traits, accounting for the phenotypic variance ranging from 12.3% to 21.8% and from 10.6% to 16.7%, respectively; and (3) six genes encoding PeTLP, PeAP2/ERF, PePUP9, PeSLP, PeHSP, and PeOCT1 proteins were identified and validated by quantitative real time polymerase chain reaction for their association with growth and oleoresin traits. These results could be useful for tree breeding and functional studies in advanced slash pine breeding program.

Slash pine is an important source of original timber and resin production on commercial forest plantations. It is necessary to implement precise breeding strategies to improve timber quality and resin yield. However, little is known about the species’ molecular genetic basis. Using a transcriptome dataset with sequencing from 240 individuals in the slash pine population, we combined multiple approaches (based on gene variation, expression variation and co-expression network) to dissect the genetic structure for slash pine major breeding traits. We found that the research population could be divided into three genetic groups with a mean heterozygosity of 0.2246. We also found that six genes with important functions in slash pine resin synthesis and timber formation through association studies. Four new SNPs associatation with the average ring width were also discovered. Our results provide new insights into the molecular genetic basis of important traits in slash pine and provide a comprehensive method for association analyses of conifer tree species with large genome.

Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus

A high content of seed glucosinolates and their degradation products imposes anti-nutritional effects on livestock; therefore, persistent efforts are made to reduce the seed GSL content to increase the commercial value of rapeseed meal. Here, we dissected the genetic structure of SGC by genome-wide association studies (GWAS) combined with transcriptome-wide association studies (TWAS). Fifteen reliable quantitative trait loci (QTLs) were identified to be associated with the reduced SGC in modern B. napus cultivars by GWAS. Analysis of the selection strength and haplotypes at these QTLs revealed that low SGC was predominantly generated by the co-selection of qGSL.A02.2, qGSL.C02.1, qGSL.A09.2, and qGSL.C09.1. Integration of the results from TWAS, comprehensive bioinformatics, and POCKET algorithm analyses indicated that BnaC02.GTR2 (BnaC02g42260D) is a candidate gene underlying qGSL.C02.1. Using CRISPR/Cas9-derived Bna.gtr2s knockout mutants, we experimentally verified that both BnaC02.GTR2 and its three paralogs positively regulate seed GSL accumulation but negatively regulated vegetative tissue GSL contents. In addition, we observed smaller seeds with higher seed oil content in these Bna.gtr2 mutants. Furthermore, both RNA-seq and correlation analyses suggested that Bna.GTR2s might play a comprehensive role in seed development, such as amino acid accumulation, GSL synthesis, sugar assimilation, and oil accumulation. This study unravels the breeding selection history of low-SGC improvement and provides new insights into the molecular function of Bna.GTR2s in both seed GSL accumulation and seed development in B. napus.

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

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

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

Uncovering Trait Associations Resulting in Maximal Seed Yield in Winter and Spring Oilseed Rape

Seed yield is a complex trait for many crop species including oilseed rape (OSR) (Brassica napus), the second most important oilseed crop worldwide. Studies have focused on the contribution of distinct factors in seed yield such as environmental cues, agronomical practices, growth conditions, or specific phenotypic traits at the whole plant level, such as number of pods in a plant. However, how female reproductive traits contribute to whole plant level traits, and hence to seed yield, has been largely ignored. Here, we describe the combined contribution of 33 phenotypic traits within a B. napus diversity set population and their trade-offs at the whole plant and organ level, along with their interaction with plant level traits. Our results revealed that both Winter OSR (WOSR) and Spring OSR (SOSR); the two more economically important OSR groups in terms of oil production; share a common dominant reproductive strategy for seed yield. In this strategy, the main inflorescence is the principal source of seed yield, producing a good number of ovules, a large number of long pods with a concomitantly high number of seeds per pod. Moreover, we observed that WOSR opted for additional reproductive strategies than SOSR, presenting more plasticity to maximise seed yield. Overall, we conclude that OSR adopts a key strategy to ensure maximal seed yield and propose an ideal ideotype highlighting crucial phenotypic traits that could be potential targets for breeding.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Comparative transcriptomics of leaves of five mulberry accessions and cataloguing structural and expression variants for future prospects

Bombyx mori, a monophagous insect, prefers leaves of the certain species of Morus more than others. The preference has been attributed to morphological and anatomical features and biochemical compounds. In the present manuscript a comparison has been made among the transcriptome of leaves of the two preferred cultivated varieties and three wild types species. While assembling, high quality transcriptomes of five genotypes were constructed with a total of 100930, 151245, 89724, 181761 and 102908 transcripts from ML, MN, MS, K2 and V1 samples respectively. Further, to compare them, orthologs were identified from these assembled transcriptome. A total of 22462, 23413, 23685, 24371, 18362, 22326, 20058, 18049, 17567 and 20518 clusters of orthologs were found in one to one comparison in KvsN, KvsL, KvsS, KvsV, LvsN, LvsS, LvsV, NvsS, NvsV, and SvsV respectively. 4236 orthologs with algebraic connectivity of 1.0 were then used to compare and to find out differentially expressed transcripts from all the genotypes. A total of 1037 transcripts expressed that include some of the important morphology, anatomy and biochemical pathways regulating transcription factors (AP2/ERFs and C2H2 Zinc fingers) and signalling components were identified to express differentially. Further, these transcriptomes were used find out markers (SSR) and variants and a total of 1101013, 537245, 970877, 310437, 675772, 338400, 581189, 751477, 514999 and 257107 variants including SNP, MNP, Insertions and deletions were found in one to one comparisons. Taken together, our data could be highly useful for mulberry community worldwide as it could be utilized in mulberry breeding programs.

Validation of a novel associative transcriptomics pipeline in Brassica oleracea: identifying candidates for vernalisation response

BACKGROUND

Associative transcriptomics has been used extensively in Brassica napus to enable the rapid identification of markers correlated with traits of interest. However, within the important vegetable crop species, Brassica oleracea, the use of associative transcriptomics has been limited due to a lack of fixed genetic resources and the difficulties in generating material due to self-incompatibility. Within Brassica vegetables, the harvestable product can be vegetative or floral tissues and therefore synchronisation of the floral transition is an important goal for growers and breeders. Vernalisation is known to be a key determinant of the floral transition, yet how different vernalisation treatments influence flowering in B. oleracea is not well understood.

RESULTS

Here, we present results from phenotyping a diverse set of 69 B. oleracea accessions for heading and flowering traits under different environmental conditions. We developed a new associative transcriptomics pipeline, and inferred and validated a population structure, for the phenotyped accessions. A genome-wide association study identified miR172D as a candidate for the vernalisation response. Gene expression marker association identified variation in expression of BoFLC.C2 as a further candidate for vernalisation response.

CONCLUSIONS

This study describes a new pipeline for performing associative transcriptomics studies in B. oleracea. Using flowering time as an example trait, it provides insights into the genetic basis of vernalisation response in B. oleracea through associative transcriptomics and confirms its characterisation as a complex G x E trait. Candidate leads were identified in miR172D and BoFLC.C2. These results could facilitate marker-based breeding efforts to produce B. oleracea lines with more synchronous heading dates, potentially leading to improved yields.

Creation and judicious application of a wheat resistance gene atlas

Disease-resistance (R) gene cloning in wheat (Triticum aestivum) has been accelerated by the recent surge of genomic resources, facilitated by advances in sequencing technologies and bioinformatics. However, with the challenges of population growth and climate change, it is vital not only to clone and functionally characterize a few handfuls of R genes, but also to do so at a scale that would facilitate the breeding and deployment of crops that can recognize the wide range of pathogen effectors that threaten agroecosystems. Pathogen populations are continually changing, and breeders must have tools and resources available to rapidly respond to those changes if we are to safeguard our daily bread. To meet this challenge, we propose the creation of a wheat R-gene atlas by an international community of researchers and breeders. The atlas would consist of an online directory from which sources of resistance could be identified and deployed to achieve more durable resistance to the major wheat pathogens, such as wheat rusts, blotch diseases, powdery mildew, and wheat blast. We present a costed proposal detailing how the interacting molecular components governing disease resistance could be captured from both the host and the pathogen through biparental mapping, mutational genomics, and whole-genome association genetics. We explore options for the configuration and genotyping of diversity panels of hexaploid and tetraploid wheat, as well as their wild relatives and major pathogens, and discuss how the atlas could inform a dynamic, durable approach to R-gene deployment. Set against the current magnitude of wheat yield losses worldwide, recently estimated at 21%, this endeavor presents one route for bringing R genes from the lab to the field at a considerable speed and quantity.

Genome-wide comparative analyses of GATA transcription factors among 19 Arabidopsis ecotype genomes: Intraspecific characteristics of GATA transcription factors

GATA transcription factors (TFs) are widespread eukaryotic regulators whose DNA-binding domain is a class IV zinc finger motif (CX2CX17-20CX2C) followed by a basic region. Due to the low cost of genome sequencing, multiple strains of specific species have been sequenced: e.g., number of plant genomes in the Plant Genome Database (http://www.plantgenome.info/) is 2,174 originated from 713 plant species. Thus, we investigated GATA TFs of 19 Arabidopsis thaliana genome-widely to understand intraspecific features of Arabidopsis GATA TFs with the pipeline of GATA database (http://gata.genefamily.info/). Numbers of GATA genes and GATA TFs of each A. thaliana genome range from 29 to 30 and from 39 to 42, respectively. Four cases of different pattern of alternative splicing forms of GATA genes among 19 A. thaliana genomes are identified. 22 of 2,195 amino acids (1.002%) from the alignment of GATA domain amino acid sequences display variations across 19 ecotype genomes. In addition, maximally four different amino acid sequences per each GATA domain identified in this study indicate that these position-specific amino acid variations may invoke intraspecific functional variations. Among 15 functionally characterized GATA genes, only five GATA genes display variations of amino acids across ecotypes of A. thaliana, implying variations of their biological roles across natural isolates of A. thaliana. PCA results from 28 characteristics of GATA genes display the four groups, same to those defined by the number of GATA genes. Topologies of bootstrapped phylogenetic trees of Arabidopsis chloroplasts and common GATA genes are mostly incongruent. Moreover, no relationship between geographical distribution and their phylogenetic relationships was found. Our results present that intraspecific variations of GATA TFs in A. thaliana are conserved and evolutionarily neutral along with 19 ecotypes, which is congruent to the fact that GATA TFs are one of the main regulators for controlling essential mechanisms, such as seed germination and hypocotyl elongation.

Character changes and Transcriptomic analysis of a cassava sexual Tetraploid

BACKGROUND

Cassava (Manihot esculenta Crantz) is an important food crop known for its high starch content. Polyploid breeding is effective in its genetic improvement, and use of 2n gametes in sexual polyploid breeding is one of the potential methods for cassava breeding and improvement. In our study, the cassava sexual tetraploid (ST), which carries numerous valuable traits, was successfully generated by hybridizing 2n female gametes SC5 (♀) and 2n male gametes SC10 (♂). However, the molecular mechanisms remain unclear. To understand these underlying molecular mechanisms behind the phenotypic alterations and heterosis in ST plants, we investigated the differences in gene expression between polyploids and diploids by determining the transcriptomes of the ST plant and its parents during the tuber root enlargement period. We also compared the characters and transcriptomes of the ST plant with its parents.

RESULTS

The ST plant was superior in plant height, stem diameter, leaf area, petiole length, plant weight, and root weight than the parent plants, except the leaf number, which was lower. The number of starch granules was higher in the roots of ST plants than those in the parent plants after five months (tuber root enlargement period), which could be due to a higher leaf net photosynthetic rate leading to early filling of starch granules. Based on transcriptome analysis, we identified 2934 and 3171 differentially expressed genes (DEGs) in the ST plant as compared to its female and male parents, respectively. Pathway enrichment analyses revealed that flavonoid biosynthesis and glycolysis/gluconeogenesis were significantly enriched in the ST plants, which might contribute to the colors of petiole (purple-red), root epidermis (dark brown), and tuber starch accumulation, respectively.

CONCLUSIONS

After sexual polyploidization, the phenotype of ST has changed significantly in comparison to their diploid parents, mainly manifest as enlarged biomass, yield, early starch filling, deep colored petiole and root epidermis. The tetraploid plants were also mature early due to early starch grain filling. Owing to enriched flavonoid biosynthesis and glycolysis/gluconeogenesis, they are possibly resistant to adversity stresses and provide better yield, respectively.

A 598-bp InDel Variation in the Promoter Region of Bna.SOC1.A05 Is Predominantly Present in Winter Type Rapeseeds

During rapeseed domestication and breeding, genetic diversity allowed to adapt it to different eco-geographical regions and to shape its useful traits. Structural variations (SVs), including presence/absence variations (PAVs), are thought to play a major role in the genetic diversity and phenotypic plasticity of rapeseed. In this study, we detected a 598-bp PAV within the promoter region of an Arabidopsis ortholog of a major flowering time gene and a downstream target of FLC, SOC1, which is one of the first genes that are upregulated in rapeseed during vernalization. Further analysis showed that the insertion is present predominantly in winter types while absent in spring types. The 589-bp sequence is present only in the A sub-genome indicating that it originated from Brassica rapa. Since the genomic region around Bna.SOC1.A05 showed a strong reduction in nucleotide diversity, the insertion might represent a larger selected sweep for rapeseed adaptation. Cis-element analysis showed that the insertion contains an ACGTG box, which is the strongest binding motif for the HY5 transcription factor in Arabidopsis. In addition, expression analyses showed that mRNA levels of Bna.SOC1.A05 were lower in accessions carrying the insertion compared to the ones that had no insertion.

Parallel analysis of global garlic gene expression and alliin content following leaf wounding

BACKGROUND

Allium sativum (garlic) is an economically important food source and medicinal plant rich in sulfides and other protective substances such as alliin, the precursor of allicin biosynthesis. Cysteine, serine and sulfur is the precursor of alliin biosynthesis. However, little is known about the alliin content under abiotic stress or the mechanism by which it is synthesized.

RESULTS

The findings revealed that the content of alliin was lowest in the garlic roots, and highest in the buds. Furthermore, alliin levels decreased in mature leaves following wounding. Transcriptome data generated over time after wounding further revealed significant up-regulation of genes integral to the biosynthetic pathways of cysteine and serine in mature garlic leaves.

CONCLUSIONS

The findings suggest that differential expression of cysteine, serine and sulfide-related genes underlies the accumulation of alliin and its precursors in garlic, providing a basis for further analyses of alliin biosynthesis.

Using natural variation to understand plant responses to iron availability

Iron bioavailability varies dramatically between soil types across the globe. This has given rise to high levels of natural variation in plant iron responses, allowing members of even a single species to thrive across a wide range of soil types. In recent years we have seen the use of genome-wide association analysis to identify natural variants underlying plant responses to changes in iron availability in both Arabidopsis and important crop species. These studies have provided insights into which genes have been important in shaping local adaptation to iron availability in different plant species and have allowed the discovery of novel regulators and mechanisms, not previously identified using mutagenesis approaches. Furthermore, these studies have allowed the identification of markers that can be used to accelerate breeding of future elite varieties with increased resilience to iron stress and improved nutritional quality. The studies highlighted here show that, in addition to studying plant responses to iron alone, it is important to consider these responses within the context of plant nutrition more broadly and to also consider iron regulation in relation to additional traits of agronomic importance such as yield and disease resistance.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

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.

Unravelling the Complex Interplay of Transcription Factors Orchestrating Seed Oil Content in Brassica napus L

Transcription factors (TFs) and their complex interplay are essential for directing specific genetic programs, such as responses to environmental stresses, tissue development, or cell differentiation by regulating gene expression. Knowledge regarding TF-TF cooperations could be promising in gaining insight into the developmental switches between the cultivars of Brassica napus L., namely Zhongshuang11 (ZS11), a double-low accession with high-oil- content, and Zhongyou821 (ZY821), a double-high accession with low-oil-content. In this regard, we analysed a time series RNA-seq data set of seed tissue from both of the cultivars by mainly focusing on the monotonically expressed genes (MEGs). The consideration of the MEGs enables the capturing of multi-stage progression processes that are orchestrated by the cooperative TFs and, thus, facilitates the understanding of the molecular mechanisms determining seed oil content. Our findings show that TF families, such as NAC, MYB, DOF, GATA, and HD-ZIP are highly involved in the seed developmental process. Particularly, their preferential partner choices as well as changes in their gene expression profiles seem to be strongly associated with the differentiation of the oil content between the two cultivars. These findings are essential in enhancing our understanding of the genetic programs in both cultivars and developing novel hypotheses for further experimental studies.

Omics: The way forward to enhance abiotic stress tolerance in Brassica napus L

Plant abiotic stresses negative affects growth and development, causing a massive reduction in global agricultural production. Rapeseed (Brassica napus L.) is a major oilseed crop because of its economic value and oilseed production. However, its productivity has been reduced by many environmental adversities. Therefore, it is a prime need to grow rapeseed cultivars, which can withstand numerous abiotic stresses. To understand the various molecular and cellular mechanisms underlying the abiotic stress tolerance and improvement in rapeseed, omics approaches have been extensively employed in recent years. This review summarized the recent advancement in genomics, transcriptomics, proteomics, metabolomics, and their imploration in abiotic stress regulation in rapeseed. Some persisting bottlenecks have been highlighted, demanding proper attention to fully explore the omics tools. Further, the potential prospects of the CRISPR/Cas9 system for genome editing to assist molecular breeding in developing abiotic stress-tolerant rapeseed genotypes have also been explained. In short, the combination of integrated omics, genome editing, and speed breeding can alter rapeseed production worldwide.

Using single-plant-omics in the field to link maize genes to functions and phenotypes

Most of our current knowledge on plant molecular biology is based on experiments in controlled laboratory environments. However, translating this knowledge from the laboratory to the field is often not straightforward, in part because field growth conditions are very different from laboratory conditions. Here, we test a new experimental design to unravel the molecular wiring of plants and study gene-phenotype relationships directly in the field. We molecularly profiled a set of individual maize plants of the same inbred background grown in the same field and used the resulting data to predict the phenotypes of individual plants and the function of maize genes. We show that the field transcriptomes of individual plants contain as much information on maize gene function as traditional laboratory-generated transcriptomes of pooled plant samples subject to controlled perturbations. Moreover, we show that field-generated transcriptome and metabolome data can be used to quantitatively predict individual plant phenotypes. Our results show that profiling individual plants in the field is a promising experimental design that could help narrow the lab-field gap.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

SeedGerm: a cost-effective phenotyping platform for automated seed imaging and machine-learning based phenotypic analysis of crop seed germination

Efficient seed germination and establishment are important traits for field and glasshouse crops. Large-scale germination experiments are laborious and prone to observer errors, leading to the necessity for automated methods. We experimented with five crop species, including tomato, pepper, Brassica, barley, and maize, and concluded an approach for large-scale germination scoring. Here, we present the SeedGerm system, which combines cost-effective hardware and open-source software for seed germination experiments, automated seed imaging, and machine-learning based phenotypic analysis. The software can process multiple image series simultaneously and produce reliable analysis of germination- and establishment-related traits, in both comma-separated values (CSV) and processed images (PNG) formats. In this article, we describe the hardware and software design in detail. We also demonstrate that SeedGerm could match specialists' scoring of radicle emergence. Germination curves were produced based on seed-level germination timing and rates rather than a fitted curve. In particular, by scoring germination across a diverse panel of Brassica napus varieties, SeedGerm implicates a gene important in abscisic acid (ABA) signalling in seeds. We compared SeedGerm with existing methods and concluded that it could have wide utilities in large-scale seed phenotyping and testing, for both research and routine seed technology applications.

Insights on SNP types, detection methods and their utilization in Brassica species: Recent progress and future perspectives

The genus Brassica, family Brassicaceae (Cruciferae), comprises many important species of oil crops, vegetables and medicinal plants including B. rapa, B. oleracea, B. nigra, B. napus, B. juncea, B. carinata. Genomic researches in Brassica species is constrained by polyploidization, mainly due to its complicated genomic structure. However, rapid development of methods for detecting single nucleotide polymorphisms (SNP), such as next generation sequencing and SNP microarray, has accelerated release of reference Brassica species genomes as well as discovery of large numbers and genome-wide SNPs, thus intensifying forward genetics in this genus. In this review, we summarize biological characteristics, classification and various methods for detecting SNPs, focusing on high-throughput techniques. Moreover, we describe the pivotal roles of SNPs in genetic diversity, linkage map construction and QTL mapping, comparative genomics, linkage disequilibrium and genome-wide association studies. These insights are expected to deepen our understanding and guide further advancements in Brassica species research.

Extended TAQing system for large-scale plant genome reorganization

We previously developed a large-scale genome restructuring technology called the TAQing system. It can induce genomic rearrangements by introducing transient and conditional formation of DNA double-strand breaks (DSBs) via heat activation of a restriction enzyme TaqI, which can cleave DNA at 5'-TCGA-3' sequences in the genome at higher temperatures (37-42°C). Such heat treatment sometimes confers lethal damage in certain plant species and TaqI cannot induce rearrangements in AT-rich regions. To overcome such problems we developed an extended TAQing (Ex-TAQing) system, which enables the use of a wider range of restriction enzymes active at standard plant-growing temperatures. We established the Ex-TAQing system using MseI that can efficiently cleave DNA at room temperature (at temperatures ranging from 22 to 25°C) and the 5'-TTAA-3' sequence which is highly abundant in the Arabidopsis genome. A synthetic intron-spanning MseI gene, which was placed downstream of a heat-shock-inducible promoter, was conditionally expressed upon milder heat treatment (33°C) to generate DSBs in Arabidopsis chromosomes. Genome resequencing revealed various types of genomic rearrangements, including copy number variations, translocation and direct end-joining at MseI cleavage sites. The Ex-TAQing system could induce large-scale rearrangements in diploids more frequently (17.4%, n = 23) than the standard TAQing system. The application of this system to tetraploids generated several strains with chromosomal rearrangements associated with beneficial phenotypes, such as high salinity stress tolerance and hypersensitivity to abscisic acid. We have developed the Ex-TAQing system, allowing more diverse patterns of genomic rearrangements, by employing various types of endonucleases and have opened a way to expand the capacity for artificial genome reorganization.

Association Mapping of Seed Quality Traits Under Varying Conditions of Nitrogen Application in Brassica juncea L. Czern & Coss

Indian mustard (Brassica juncea) is a major source of vegetable oil in the Indian subcontinent. The seed cake left after the oil extraction is used as livestock feed. We examined the genetic architecture of oil, protein, and glucosinolates by conducting a genome-wide association study (GWAS), using an association panel comprising 92 diverse genotypes. We conducted trait phenotyping over 2 years at two levels of nitrogen (N) application. Genotyping by sequencing was used to identify 66,835 loci, covering 18 chromosomes. Genetic diversity and phenotypic variations were high for the studied traits. Trait performances were stable when averaged over years and N levels. However, individual performances differed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits. Population structure, principal components (PCs) analysis, and discriminant analysis of principal components (DAPCs) were included as covariates to overcome the bias due to the population stratification. We identified 16, 23, and 27 loci associated with oil, protein, and glucosinolates, respectively. We also established LD patterns and haplotype structures for the candidate genes. The average block sizes were larger on A-genome chromosomes as compared to the B- genome chromosomes. Genetic associations differed over N levels. However, meta-analysis of GWAS datasets not only improved the power to recognize associations but also helped to identify common SNPs for oil and protein contents. Annotation of the genomic region around the identified SNPs led to the prediction of 21 orthologs of the functional candidate genes related to the biosynthesis of oil, protein, and glucosinolates. Notable among these are: LACS5 (A09), FAD6 (B05), ASN1 (A06), GTR2 (A06), CYP81G1 (B06), and MYB44 (B06). The identified loci will be very useful for marker-aided breeding for seed quality modifications in B. juncea.

A transcriptomic variation map provides insights into the genetic basis of Pinus massoniana Lamb. evolution and the association with oleoresin yield

BACKGROUND

Masson pine (Pinus massoniana Lamb.), the dominant native coniferous species in southern China, is commercially important for supplying timber and oleoresin. However, knowledge of the genetic variability of masson pine germplasm is still limited. In this study, the genetic diversity and population structure of masson pine germplasm were assessed using 204 wild accessions from 10 main distribution regions using 94,194 core single-nucleotide polymorphisms (SNPs) obtained from transcriptome sequencing data.

RESULTS

The average expected heterozygosity was 0.2724, implying abundant genetic diversity within masson pine germplasm. Analysis of molecular variance (AMOVA) revealed that 3.29% of the variation was sourced from genetic differentiation. Structure analysis identified two geographically distinct groups. Discriminant analysis of principal components (DAPC) showed that one of those groups was further divided into two clusters. Sichuan and Chongqing provenance is the geographical origin, which diffused outward along two different lines. Oleoresin yield is reflected in the evolution of the two groups, and exhibits two different trends along the two lines of diffusion. The oleoresin yield may be associated with the genes of chitinase, CYP720B, cytochrome P450, ABC transporter, and AP2/ethylene-responsive transcription factor (ERF) based on SNPs and expression.

CONCLUSIONS

SNP markers from transcriptome sequencing are highly capable of evaluating genetic diversity within different species, as well as the genetic control of objective traits. The functions of these genes will be verified in future studies, and those genes strongly associated with oleoresin yield will be used to improve yields by means of early genotype selection and genetic engineering.

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

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