Identification of Putative Candidate Genes for Water Stress Tolerance in Canola (Brassica napus)

Combined GWAS and eGWAS reveals the genetic basis underlying drought tolerance in emmer wheat (Triticum turgidum L.)

Drought is one of the major environmental constraints for wheat production world-wide. As the progenitor and genetic reservoir of common wheat, emmer wheat is considered as an invaluable gene pool for breeding drought-tolerant wheat. Combining GWAS and eGWAS analysis of 107 accessions, we identified 86 QTLs, 105 462 eQTLs as well as 68 eQTL hotspots associating with drought tolerance (DT) in emmer wheat. A complex regulatory network composed of 185 upstream regulator and 2432 downstream drought-responsive candidates was developed, of which TtOTS1 was found to play a negative effect in determining DT through affecting root development. This study sheds light on revealing the genetic basis underlying DT, which will provide the indispensable genes and germplasm resources for elite drought tolerance wheat improvement and breeding.

Identification of novel genes associated with herbicide tolerance in Lentil (Lens culinaris ssp. culinaris Medik.)

Weeds pose a major constraint in lentil cultivation, leading to decrease farmers' revenues by reducing the yield and increasing the management costs. The development of herbicide tolerant cultivars is essential to increase lentil yield. Even though herbicide tolerant lines have been identified in lentils, breeding efforts are still limited and lack proper validation. Marker assisted selection (MAS) can increase selection accuracy at early generations. Total 292 lentil accessions were evaluated under different dosages of two herbicides, metribuzin and imazethapyr, during two seasons at Marchouch, Morocco and Terbol, Lebanon. Highly significant differences among accessions were observed for days to flowering (DF) and maturity (DM), plant height (PH), biological yield (BY), seed yield (SY), number of pods per plant (NP), as well as the reduction indices (RI) for PH, BY, SY and NP. A total of 10,271 SNPs markers uniformly distributed along the lentil genome were assayed using Multispecies Pulse SNP chip developed at Agriculture Victoria, Melbourne. Meta-GWAS analysis was used to detect marker-trait associations, which detected 125 SNPs markers associated with different traits and clustered in 85 unique quantitative trait loci. These findings provide valuable insights for initiating MAS programs aiming to enhance herbicide tolerance in lentil crop.

Exogenous melatonin improves germination rate in buckwheat under high temperature stress by regulating seed physiological and biochemical characteristics

The germinations of three common buckwheat (Fagopyrum esculentum) varieties and two Tartary buckwheat (Fagopyrum tataricum) varieties seeds are known to be affected by high temperature. However, little is known about the physiological mechanism affecting germination and the effect of melatonin (MT) on buckwheat seed germination under high temperature. This work studied the effects of exogenous MT on buckwheat seed germination under high temperature. MT was sprayed. The parameters, including growth, and physiological factors, were examined. The results showed that exogenous MT significantly increased the germination rate (GR), germination potential (GP), radicle length (RL), and fresh weight (FW) of these buckwheat seeds under high-temperature stress and enhanced the content of osmotic adjustment substances and enzyme activity. Comprehensive analysis revealed that under high-temperature stress during germination, antioxidant enzymes play a predominant role, while osmotic adjustment substances work synergistically to reduce the extent of damage to the membrane structure, serving as the primary key indicators for studying high-temperature resistance. Consequently, our results showed that MT had a positive protective effect on buckwheat seeds exposed to high temperature stress, providing a theoretical basis for improving the ability to adapt to high temperature environments.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Alleviating sweetpotato salt tolerance through exogenous glutathione and melatonin: A profound mechanism for active oxygen detoxification and preservation of photosynthetic organs

Salt stress profoundly impacts sweetpotato production. Exogenous glutathione (GSH) and melatonin (MT) promoted plant growth under stress, but their specific roles and mechanisms in sweetpotato salt tolerance need exploration. This study investigated GSH and MT's regulatory mechanisms in sweetpotato under salt stress. Salt stress significantly reduces both growth and biomass by hindering photosynthesis, root traits, K+ content, and K+/Na+ balance, leading to oxidative stress and excessive hydrogen peroxide (H2O2), superoxide ion (O2•-), and malondialdehyde (MDA) production and Na+ accumulation. Nevertheless, GSH (2 mM) and MT (25 μM) pre-treatments effectively mitigated salt-induced oxidative damage and protected the plasma membrane. They reduced osmotic pressure by enhancing K+ uptake, K+/Na+ regulation, osmolyte accumulation, and reducing Na+ accumulation. Improved stomatal traits, chloroplast and grana lamella preservation, and maintenance of mesophyll cells, cell wall, and mitochondrial structure were observed with GSH and MT pre-treatments under salt stress, therefore boosting the photosynthetic system and enhancing plant growth and biomass. Moreover, the findings also indicate that the positive outcomes of GSH and MT pre-treatments result from elevated antioxidant levels, enhanced enzymatic activity, and upregulated expression of sodium hydrogen exchanger 2 (NHX2), K+transporter 1 (AKT1), and cation/H+exchanger (CHX), CBL-interacting protein kinase 1 (CIPK1), and antioxidant enzyme genes. These mechanisms enhance structural stability in photosynthesis and reduce salt stress. Evidently, MT pre-treatment exhibited superior effects compared to GSH. These findings provide a firm theoretical basis for employing GSH and MT to enhance salt tolerance in sweetpotato cultivation.

Role of Microbes in Alleviating Crop Drought Stress: A Review

Drought stress is an annual global phenomenon that has devastating effects on crop production, so numerous studies have been conducted to improve crop drought resistance. Plant-associated microbiota play a crucial role in crop health and growth; however, we have a limited understanding of the key processes involved in microbiome-induced crop adaptation to drought stress. In this review, we summarize the adverse effects of drought stress on crop growth in terms of germination, photosynthesis, nutrient uptake, biomass, and yield, with a focus on the response of soil microbial communities to drought stress and plant-microbe interactions under drought stress. Moreover, we review the morpho-physiological, biochemical, and molecular mechanisms underlying the mitigation effect of microbes on crop drought stress. Finally, we highlight future research directions, including the characterization of specific rhizosphere microbiome species with corresponding root exudates and the efficiency of rhizobacteria inoculants under drought conditions. Such research will advance our understanding of the complex interactions between crops and microbes and improve crop resistance to drought stress through the application of beneficial drought-adaptive microbes.

Integration of genome-wide association studies, metabolomics, and transcriptomics reveals phenolic acid- and flavonoid-associated genes and their regulatory elements under drought stress in rapeseed flowers

Introduction

Biochemical and metabolic processes help plants tolerate the adverse effects of drought. In plants accumulating bioactive compounds, understanding the genetic control of the biosynthesis of biochemical pathways helps the discovery of candidate gene (CG)-metabolite relationships.

Methods

The metabolic profile of flowers in 119 rapeseed (Brassica napus) accessions was assessed over two irrigation treatments, one a well-watered (WW) condition and the other a drought stress (DS) regime. We integrated information gained from 52,157 single-nucleotide polymorphism (SNP) markers, metabolites, and transcriptomes to identify linked SNPs and CGs responsible for the genetic control of flower phenolic compounds and regulatory elements.

Results

In a genome-wide association study (GWAS), of the SNPs tested, 29,310 SNPs were qualified to assess the population structure and linkage disequilibrium (LD), of which several SNPs for radical scavenging activity (RSA) and total flavanol content (TFLC) were common between the two irrigation conditions and pleiotropic SNPs were found for chlorogenic and coumaric acids content. The principal component analysis (PCA) and stepwise regression showed that chlorogenic acid and epicatechin in WW and myricetin in DS conditions were the most important components for RSA. The hierarchical cluster analysis (HCA) showed that vanillic acid, myricetin, gallic acid, and catechin were closely associated in both irrigation conditions. Analysis of GWAS showed that 60 CGs were identified, of which 18 were involved in stress-induced pathways, phenylpropanoid pathway, and flavonoid modifications. Of the CGs, PAL1, CHI, UGT89B1, FLS3, CCR1, and CYP75B137 contributed to flavonoid biosynthetic pathways. The results of RNA sequencing (RNA-seq) revealed that the transcript levels of PAL, CHI, and CYP75B137 known as early flavonoid biosynthesis-related genes and FLS3, CCR1, and UGT89B1 related to the later stages were increased during drought conditions. The transcription factors (TFs) NAC035 and ERF119 related to flavonoids and phenolic acids were upregulated under drought conditions.

Discussion

These findings expand our knowledge on the response mechanisms to DS, particularly regarding the regulation of key phenolic biosynthetic genes in rapeseed. Our data also provided specific linked SNPs for marker-assisted selection (MAS) programs and CGs as resources toward realizing metabolomics-associated breeding of rapeseed.

Comparative transcriptome analysis of two contrasting genotypes provides new insights into the drought response mechanism in pigeon pea (Cajanus cajan L. Millsp.)

BACKGROUND

Despite plant's ability to adapt and withstand challenging environments, drought poses a severe threat to their growth and development. Although pigeon pea is already quite resistant to drought, the prolonged dehydration induced by the aberrant climate poses a serious threat to their survival and productivity.

OBJECTIVE

Comparative physiological and transcriptome analyses of drought-tolerant (CO5) and drought-sensitive (CO1) pigeon pea genotypes subjected to drought stress were carried out in order to understand the molecular basis of drought tolerance in pigeon pea.

METHODS

The transcriptomic analysis allowed us to examine how drought affects the gene expression of C. cajan. Using bioinformatics tools, the unigenes were de novo assembled, annotated, and functionally evaluated. Additionally, a homology-based sequence search against the droughtDB database was performed to identify the orthologs of the DEGs.

RESULTS

1102 potential drought-responsive genes were found to be differentially expressed genes (DEGs) between drought-tolerant and drought-sensitive genotypes. These included Abscisic acid insensitive 5 (ABI5), Nuclear transcription factor Y subunit A-7 (NF-YA7), WD40 repeat-containing protein 55 (WDR55), Anthocyanidin reductase (ANR) and Zinc-finger homeodomain protein 6 (ZF-HD6) and were highly expressed in the tolerant genotype. Further, GO analysis revealed that the most enriched classes belonged to biosynthetic and metabolic processes in the biological process category, binding and catalytic activity in the molecular function category and nucleus and protein-containing complex in the cellular component category. Results of KEGG pathway analysis revealed that the DEGs were significantly abundant in signalling pathways such as plant hormone signal transduction and MAPK signalling pathways. Consequently, in our investigation, we have identified and validated by qPCR a group of genes involved in signal reception and propagation, stress-specific TFs, and basal regulatory genes associated with drought response.

CONCLUSION

In conclusion, our comprehensive transcriptome dataset enabled the discovery of candidate genes connected to pathways involved in pigeon pea drought response. Our research uncovered a number of unidentified genes and transcription factors that could be used to understand and improve susceptibility to drought.

Genetic Diversity and Population Structure in Ethiopian Mustard (Brassica carinata A. Braun) as Revealed by Single Nucleotide Polymorphism Markers

Ethiopian mustard (Brassica carinata A. Braun) is currently one of the potential oilseeds dedicated to the production for biofuel and other bio-industrial applications. The crop is assumed to be native to Ethiopia where a number of diversified B. carinata germplasms are found and conserved ex situ. However, there is very limited information on the genetic diversity and population structure of the species. This study aimed to investigate the genetic diversity and population structure of B. carinata genotypes of different origins using high-throughput single nucleotide polymorphism (SNP) markers. We used Brassica 90K Illumina InfiniumTM SNP array for genotyping 90 B. carinata genotypes, and a total of 11,499 informative SNP markers were used for investigating the population structure and genetic diversity. The structure analysis, principal coordinate analysis (PcoA) and neighbor-joining tree analysis clustered the 90 B. carinata genotypes into two distinct subpopulations (Pop1 and Pop2). The majority of accessions (65%) were clustered in Pop1, mainly obtained from Oromia and South West Ethiopian People (SWEP) regions. Pop2 constituted dominantly of breeding lines and varieties, implying target selection contributed to the formation of distinct populations. Analysis of molecular variance (AMOVA) revealed a higher genetic variation (93%) within populations than between populations (7%), with low genetic differentiation (PhiPT = 0.07) and poor correlation between genetic and geographical distance (R = 0.02). This implies the presence of gene flow (Nm > 1) and weak geographical structure of accessions. Genetic diversity indices showed the presence of moderate genetic diversity in B. carinata populations with an average genetic diversity value (HE = 0.31) and polymorphism information content (PIC = 0.26). The findings of this study provide important and relevant information for future breeding and conservation efforts of B. carinata.

Transcriptome analysis reveals key regulatory genes for root growth related to potassium utilization efficiency in rapeseed (Brassica napus L.)

Root system architecture (RSA) is the primary predictor of nutrient intake and significantly influences potassium utilization efficiency (KUE). Uncertainty persists regarding the genetic factors governing root growth in rapeseed. The root transcriptome analysis reveals the genetic basis driving crop root growth. In this study, RNA-seq was used to profile the overall transcriptome in the root tissue of 20 Brassica napus accessions with high and low KUE. 71,437 genes in the roots displayed variable expression profiles between the two contrasting genotype groups. The 212 genes that had varied expression levels between the high and low KUE lines were found using a pairwise comparison approach. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional classification analysis revealed that the DEGs implicated in hormone and signaling pathways, as well as glucose, lipid, and amino acid metabolism, were all differently regulated in the rapeseed root system. Additionally, we discovered 33 transcription factors (TFs) that control root development were differentially expressed. By combining differential expression analysis, weighted gene co-expression network analysis (WGCNA), and recent genome-wide association study (GWAS) results, four candidate genes were identified as essential hub genes. These potential genes were located fewer than 100 kb from the peak SNPs of QTL clusters, and it was hypothesized that they regulated the formation of the root system. Three of the four hub genes' homologs-BnaC04G0560400ZS, BnaC04G0560400ZS, and BnaA03G0073500ZS-have been shown to control root development in earlier research. The information produced by our transcriptome profiling could be useful in revealing the molecular processes involved in the growth of rapeseed roots in response to KUE.

Melatonin-priming enhances maize seedling drought tolerance by regulating the antioxidant defense system

Drought stress (DS) challenges sustainable agriculture production by limiting crop growth and development. The objective of the study was to evaluate the effect of melatonin-priming on enzymatic and non-enzymatic antioxidant defense mechanisms and its relation with leaf ultrastructure and stomatal traits in maize (Zea mays L) seedlings under DS (PEG-6000). DS drastically decreased seed germination, plant growth, and leaf chlorophyll content due to excessive reactive oxygen species (ROS) production. Melatonin-priming significantly (P < 0.05) increased seed germination, root length, shoot length, fresh seedling weight, proline content, total soluble protein content, sugar content, chlorophyll content, and stomatal aperture size by 101%, 30%, 133%, 51%, 22%, 59%, 54%, 20%, and 424%, compared to no priming (NP) under DS, respectively. Similarly, priming improved leaf ultrastructure and reduced the amount of chlorophyll loss and oxidative damage in maize seedlings. Melatonin seed priming with 500 µM melatonin (M2) greatly increased superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), glutathione, and ascorbate (AsA) activity, by 65%, 63%, 94%, 41%, and 55% compared to NP under DS and by 0.26%, 8%, 33%, 42%, and 15% under no-stress (NS), respectively. Melatonin-priming also reduced malondialdehyde content, electrolyte leakage, hydrogen peroxide (H2O2) content, and superoxide anion (O2-) content by 26%, 31%, 31%, and 33% compared to NP under DS and by 8%, 18%, 10%, and 11% under NS, respectively. In response to DS, melatonin-priming also stabilized the chloroplast structure, sustained cell expansion, protected cell walls, and greatly improved stomatal traits, including stomatal number, length, and width. Our results suggest that melatonin-priming improves drought tolerance in maize seedlings by alleviating the negative effect of ROS.

Differential responses of contrasting low phosphorus tolerant cotton genotypes under low phosphorus and drought stress

BACKGROUND

Drought is one of the main reasons for low phosphorus (P) solubility and availability.

AIMS

The use of low P tolerant cotton genotypes might be a possible option to grow in drought conditions.

METHODS

This study investigates the tolerance to drought stress in contrasting low P-tolerant cotton genotypes (Jimian169; strong tolerant to low P and DES926; weak tolerant to low P). In hydroponic culture, the drought was artificially induced with 10% PEG in both cotton genotypes followed by low (0.01 mM KH₂PO₄) and normal (1 mM KH₂PO₄) P application.

RESULTS

The results showed that under low P, PEG-induced drought greatly inhibited growth, dry matter production, photosynthesis, P use efficiency, and led to oxidative stress from excessive malondialdehyde (MDA) and higher accumulation of reactive oxygen species (ROS), and these effects were more in DES926 than Jimian169. Moreover, Jimian169 alleviated oxidative damage by improving the antioxidant system, photosynthetic activities, and an increase in the levels of osmoprotectants like free amino acids, total soluble proteins, total soluble sugars, and proline.

CONCLUSIONS

The present study suggests that the low P-tolerant cotton genotype can tolerate drought conditions through high photosynthesis, antioxidant capacity, and osmotic adjustment.

Genome-Wide association analysis of phenotypic traits in Bambara groundnut under drought-stressed and non-stressed conditions based on DArTseq SNP

INTRODUCTION

Bambara groundnut (BG) (Vigna subterranea [L.] Verdc) is an indigenous, resilient, but underutilized leguminous crop that occurs mostly as genetically heterogeneous landraces with limited information on the drought tolerant attributes. This study elucidates the associations between sequencing-based diversity array technology (DArTseq) and phenotypic character as well as differing indices related to drought tolerance in one hundred accessions of Bambara groundnut.

METHODS

The field experiments were conducted at IITA research stations in Kano and Ibadan between 2016 and 2018 planting seasons. The experiments were arranged in randomised complete block design with three replications, under the different water regimes. The phenotypic traits evaluated was further to construct the dendrogram. Genome-wide association mapping was conducted based on 5927 DArTs loci with < 20% missing data.

RESULTS AND DISCUSSIONS

The genome wide association study predicted drought tolerance in Bambara accessions for geometric mean productivity (GMP) and stress tolerance index (STI). TVSu-423 had the highest GMP and STI values (28.50, 2.40), while TVSu-2017 had the lowest at GMP (1.74) and STI (0.01) respectively. The relative water content (%) was significantly higher for accessions; TVSu-266 (60.35, 61.49), TVSu-2 (58.29, 53.94), and TVSu-411 (55.17, 58.92) in 2016/2017 and 2017/2018, respectively. The phenotypic characters studied delineated the accessions into two major clusters and five distinct sub-clusters, indicating variations across all the geographical locations. The 5,927 DArTseq genomic markers in association with STI further grouped the 100 accessions into two main clusters. TVSu-1897 from Botswana (Southern Africa) was in the first cluster, while the remaining 99 accessions from Western, Central, and Eastern Africa made up the second cluster. The eight significant Quantitative Trait Loci (QTLs) (24346377|F|0-22:A>G-22:A>G, 24384105|F|0-56:A>G33 :A> G, 24385643|F|0-53:G>C-53:G>C, 24385696|F|0-43:A>G-43:A>G, 4177257|F|0-44:A>T-44:A>T, 4182070|F|0-66:G>A-66:G>A, 4183483|F|0-24:G>A-24:G>A, 4183904|F|0-11:C>T-11:C>T) identified with Bonferroni threshold was in association with STI, indicative of variations under the drought-stressed condition. The observation of consistent SNPs in the 2016 and 2017 planting seasons, as well as in combination with the 2016 and 2017 planting seasons, led to the designation of these QTLs as significant. The drought selected accessions could form basis for hybridization breeding. The identified quantitative trait loci could be useful in marker-assisted selection in drought molecular breeding programs.

Identification of candidate genes for LepR1 resistance against Leptosphaeria maculans in Brassica napus

Utilising resistance (R) genes, such as LepR1, against Leptosphaeria maculans, the causal agent of blackleg in canola (Brassica napus), could help manage the disease in the field and increase crop yield. Here we present a genome wide association study (GWAS) in B. napus to identify LepR1 candidate genes. Disease phenotyping of 104 B. napus genotypes revealed 30 resistant and 74 susceptible lines. Whole genome re-sequencing of these cultivars yielded over 3 million high quality single nucleotide polymorphisms (SNPs). GWAS in mixed linear model (MLM) revealed a total of 2,166 significant SNPs associated with LepR1 resistance. Of these SNPs, 2108 (97%) were found on chromosome A02 of B. napus cv. Darmor bzh v9 with a delineated LepR1_mlm1 QTL at 15.11-26.08 Mb. In LepR1_mlm1, there are 30 resistance gene analogs (RGAs) (13 nucleotide-binding site-leucine rich repeats (NLRs), 12 receptor-like kinases (RLKs), and 5 transmembrane-coiled-coil (TM-CCs)). Sequence analysis of alleles in resistant and susceptible lines was undertaken to identify candidate genes. This research provides insights into blackleg resistance in B. napus and assists identification of the functional LepR1 blackleg resistance gene.

Differential alternative splicing genes and isoform co-expression networks of Brassica napus under multiple abiotic stresses

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.

Genome-wide association, RNA-seq and iTRAQ analyses identify candidate genes controlling radicle length of wheat

The radicle, present in the embryo of a seed, is the first root to emerge at germination, and its rapid growth is essential for establishment and survival of the seedling. However, there are few studies on the critical mechanisms underlying radicle and then radicle length in wheat seedlings, despite its importance as a food crop throughout the world. In the present study, 196 wheat accessions from the Huanghuai Wheat Region were screened to measure radicle length under 4 hydroponic culture environments over 3 years. Different expression genes and proteins (DEGs/DEPs) between accessions with extremely long [Yunong 949 (WRL1), Zhongyu 9,302 (WRL2)] and short roots [Yunong 201 (WRS1), Beijing 841 (WRS2)] were identified in 12 sets of root tissue samples by RNA-seq and iTRAQ (Isobaric tags for relative and absolute quantification). Phenotypic results showed that the elongation zone was significantly longer in root accessions with long roots compared to the short-rooted accessions. A genome-wide association study (GWAS) identified four stable chromosomal regions significantly associated with radicle length, among which 1A, 4A, and 7A chromosomes regions explained 7.17% to12.93% of the phenotypic variation. The omics studies identified the expression patterns of 24 DEGs/DEPs changed at both the transcriptional and protein levels. These DEGs/DEPs were mainly involved in carbon fixation in photosynthetic organisms, photosynthesis and phenylpropanoid biosynthesis pathways. TraesCS1A02G104100 and TraesCS2B02G519100 were involved in the biosynthesis of tricin-lignins in cell walls and may affect the extension of cell walls in the radicle elongation zone. A combination of GWAS and RNA-seq analyses revealed 19 DEGs with expression changes in the four accessions, among which, TraesCS1A02G422700 (a cysteine-rich receptor-like protein kinase 6, CRK6) also showed upregulation in the comparison group by RNA-seq, iTRAQ, and qRT-PCR. BSMV-mediated gene silencing also showed that TaCRK6 improves root development in wheat. Our data suggest that TaCRK6 is a candidate gene regulating radicle length in wheat.

Gene Differential Co-Expression Networks Based on RNA-Seq: Construction and Its Applications

Gene co-expression network (GCN) becomes an increasingly important tool in omics data analysis. A great challenge for GCN construction is that the sample size is far lower than the number of genes. Traditional methods rely on considerable samples. Moreover, association signals are likely weak, nonlinear and stochastic, which are difficult to be identified among thousands of candidates. In this paper, the gray correlation coefficient (GCC) is introduced, and a novel method to construct gene differential co-expression networks (GDCNs) is proposed. Based on the GDCNs, three measures are proposed to explore informative genes. The proposed method can make full use of the information provided by a handful of samples and overcome the shortages of GCNs, which can evaluate the changes of co-expression relationships that are possibly triggered by treatments. Based on RNA-seq data of Brassica napus, GDCNs under multiple experimental conditions are constructed and investigated. It is found that the GCC-based method is very robust to data processing. The GDCNs facilitate the inference of gene functions and the identification of informative genes that are responsible for stress responsiveness. The GDCN-based approaches integrate the 'guilt by association' and the 'guilt by rewiring' rules, which provide alternative tools for omics data analysis.

Integration of GWAS and transcriptome analyses to identify SNPs and candidate genes for aluminum tolerance in rapeseed (Brassica napus L.)

BACKGROUND

The exchangeable aluminum (Al), released from the acid soils, is another addition to the environmental stress factors in the form of Al toxicity stress. Al stress affects the normal crop development and reduces the overall yield of rapeseed (Brassica napus L.). The response mechanism of plants to Al toxicity is complicated and difficult to understand with few QTL related studies in rapeseed under Al toxicity stress.

RESULT

Using 200,510 SNPs developed by SLAF-seq (specific-locus amplified fragment sequencing) technology, we carried out the genome-wide association analysis (GWAS) in a population of 254 inbred lines of B. napus with large genetic variation and Al-tolerance differences. There were 43 SNPs significantly associated with eight Al-tolerance traits in the seedling stage were detected on 14 chromosomes, and 777 candidate genes were screened at the flanking 100 kb region of these SNPs. Moreover, RNA-seq detected 8291 and 5341 DEGs (the differentially expressed gene) in the Al -tolerant line (ATL) and -sensitive line (ASL), respectively. Based on integration of GWAS and RNA-seq analysis, 64 candidate genes from GWAS analysis differentially expressed at least once in 6 h vs 0 h or 24 h vs 0 h conditions in ATL or ASL. Moreover, four out of sixty-four candidate genes (BnaA03g30320D, BnaA10g11500D, BnaC03g38360D and BnaC06g30030D) were differentially expressed in both 6 h and 24 h compared to 0 h (control) conditions in both lines. The proposed model based on the candidate genes excavated in this study highlighted that Al stress disturb the oxidation-redox balance, causing abnormal synthesis and repair of cell wall and ABA signal transduction, ultimately resulting in inhibition of root elongation.

CONCLUSIONS

The integration of GWAS and transcriptome analysis provide an effective strategy to explore the SNPs and candidate genes, which has a potential to develop molecular markers for breeding Al tolerant rapeseed varieties along with theoretical basis of molecular mechanisms for Al toxicity response of Brassica napus plants.

Transcriptomic, Metabolomic and Ionomic Analyses Reveal Early Modulation of Leaf Mineral Content in Brassica napus under Mild or Severe Drought

While it is generally acknowledged that drought is one of the main abiotic factors affecting plant growth, how mineral nutrition is specifically and negatively affected by water deficit has received very little attention, other than being analyzed as a consequence of reduced growth. Therefore, Brassica napus plants were subjected to a gradual onset of water deficits (mild, severe, or severe extended), and leaves were analyzed at the ionomic, transcriptomic and metabolic levels. The number of Differentially Expressed Genes (DEGs) and of the most differentially accumulated metabolites increased from mild (525 DEGs, 57 metabolites) to severe (5454 DEGs, 78 metabolites) and severe extended (9346 DEGs, 95 metabolites) water deficit. Gene ontology enrichment analysis of the 11,747 DEGs identified revealed that ion transport was one of the most significant processes affected, even under mild water deficit, and this was also confirmed by the shift in ionomic composition (mostly micronutrients with a strong decrease in Mo, Fe, Zn, and Mn in leaves) that occurred well before growth reduction. The metabolomic data and most of the transcriptomic data suggested that well-known early leaf responses to drought such as phytohormone metabolism (ABA and JA), proline accumulation, and oxidative stress defense were induced later than repression of genes related to nutrient transport.

Genomic background selection to reduce the mutation load after random mutagenesis

Random mutagenesis is a standard procedure to increase allelic variation in a crop species, especially in countries where the use of genetically modified crops is limited due to legal constraints. The chemical mutagen EMS is used in many species to induce random mutations throughout the genome with high mutation density. The major drawback for functional analysis is a high background mutation load in a single plant that must be eliminated by subsequent backcrossing, a time and resource-intensive activity. Here, we demonstrate that genomic background selection combined with marker-assisted selection is an efficient way to select individuals with reduced background mutations within a short period. We identified BC1 plants with a significantly higher share of the recurrent parent genome, thus saving one backcross generation. Furthermore, spring rapeseed as the recurrent parent in a backcrossing program could accelerate breeding by reducing the generation cycle. Our study depicts the potential for reducing the background mutation load while accelerating the generation cycle in EMS-induced winter oilseed rape populations by integrating genomic background selection.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

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

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

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

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

Physiological and Biochemical Characterization of the GABA Shunt Pathway in Pea (Pisum sativum L.) Seedlings under Drought Stress

The physiological and biochemical role of the γ-aminobutyric acid (GABA) shunt pathway in green pea seedlings (Pisum sativum L.) was studied in response to soil water holding capacity levels: 80%, 60%, 40%, 20%, and 10% grown under continuous light at 25 °C for 7 days and 14 days, separately. Characterization of seeds germination pattern, seedlings growth (plant height, fresh and dry weight, and chlorophyll contents), GABA shunt metabolite (GABA, glutamate, and alanine) levels, total protein and carbohydrate levels, and oxidative damage (MDA level) were examined. Data showed a significant effect of drought stress on seed germination, plant growth, GABA shunt metabolites level, total protein and carbohydrate contents, and MDA level. A significant decline in seed germination percentage was recorded at a 20% drought level, which indicated that 20% of soil water holding capacity is the threshold value of water availability for normal germination after 14 days. Seedling fresh weight, dry weight, and plant height were significantly reduced with a positive correlation as water availability was decreased. There was a significant decrease with a positive correlation in Chl a and Chl b contents in response to 7 days and 14 days of drought. GABA shunt metabolites were significantly increased with a negative correlation as water availability decreased. Pea seedlings showed a significant increase in protein content as drought stress was increased. Total carbohydrate levels increased significantly when the amount of water availability decreased. MDA content increased slightly but significantly after 7 days and sharply after 14 days under all water stress levels. The maximum increase in MDA content was observed at 20% and 10% water levels. Overall, the significant increases in GABA, protein and carbohydrate contents were to cope with the physiological impact of drought stress on Pisum sativum L. seedlings by maintaining cellular osmotic adjustment, protecting plants from oxidative stress, balancing carbon and nitrogen (C:N) metabolism, and maintaining cell metabolic homeostasis and cell turgor. The results presented in this study indicated that severe (less than 40% water content of the holding capacity) and long-term drought stress should be avoided during the germination stage to ensure proper seedling growth and metabolism in Pisum sativum L.

Genome-Wide Association Mapping Unravels the Genetic Control of Seed Vigor under Low-Temperature Conditions in Rapeseed (Brassica napus L.)

Low temperature inhibits rapid germination and successful seedling establishment of rapeseed (Brassica napus L.), leading to significant productivity losses. Little is known about the genetic diversity for seed vigor under low-temperature conditions in rapeseed, which motivated our investigation of 13 seed germination- and emergence-related traits under normal and low-temperature conditions for 442 diverse rapeseed accessions. The stress tolerance index was calculated for each trait based on performance under non-stress and low-temperature stress conditions. Principal component analysis of the low-temperature stress tolerance indices identified five principal components that captured 100% of the seedling response to low temperature. A genome-wide association study using ~8 million SNP (single-nucleotide polymorphism) markers identified from genome resequencing was undertaken to uncover the genetic basis of seed vigor related traits in rapeseed. We detected 22 quantitative trait loci (QTLs) significantly associated with stress tolerance indices regarding seed vigor under low-temperature stress. Scrutiny of the genes in these QTL regions identified 62 candidate genes related to specific stress tolerance indices of seed vigor, and the majority were involved in DNA repair, RNA translation, mitochondrial activation and energy generation, ubiquitination and degradation of protein reserve, antioxidant system, and plant hormone and signal transduction. The high effect variation and haplotype-based effect of these candidate genes were evaluated, and high priority could be given to the candidate genes BnaA03g40290D, BnaA06g07530D, BnaA09g06240D, BnaA09g06250D, and BnaC02g10720D in further study. These findings should be useful for marker-assisted breeding and genomic selection of rapeseed to increase seed vigor under low-temperature stress.

QTL Mapping for Seed Germination Response to Drought Stress in Brassica napus

Drought stress is one of the most environmental abiotic stresses affecting seed germination and crop growth. In the present study, the genetic characteristics of seed germination under drought stress in a Brassica napus double haploid population were analyzed. Five germination-related indexes, including germination percentage (GP), root length (RL), shoot length (SL), fresh weight (FW), and root-to-shoot length ratio (R/S) under control and drought stress, were calculated, and the drought stress index (DSI), including DSI-GP, DSI-RL, DSI-SL, DSI-FW, and DSI-R/S, was determined using the quantitative trait loci (QTLs) analysis based on high-density genetic linkage map. The phenotypic analysis indicated that the R/S is an effective morphological trait in the determination of drought tolerance in the seedling stage. Thirty-nine identified QTLs were observed for these traits and then integrated into 36 consensus QTLs, in which 18 QTLs were found to affect the DSI of four traits (GP, RL, SL, and R/S). Based on the co-linearity between genetic and physical maps of B. napus, 256 candidate genes were detected, and 128 genes have single-nucleotidepolymorphisms/insertion-deletion (SNP/InDel) variations between two parents, some of which were associated with the drought stress tolerance (for example, BnaC03g32780D, BnaC03g37030D, and BnaC09g27300D). The present results laid insights into drought tolerance and its genetic bases in B. napus.

Transcriptomics reveal high regulatory diversity of drought tolerance strategies in a biennial oil crop

Spring droughts are expected to become more frequent in Central Europe as a result of climate change. Their coincidence with flowering of biennial crops like winter oilseed rape (Brassica napus) can cause major impact for yield development. However, no data is available on the diversity of genetic regulation of drought tolerance during this stage under realistic conditions. Here, we assessed the phenotypic plasticity of drought response for eight diverse B. napus accessions under field-like conditions and linked their stress response to gene and miRNA expression during early and late stress. We observed highly diverse responses, both on the phenotypic and on the gene expression level. Our data suggest that drought tolerant accessions have more effective molecular protection mechanisms like ROS scavenging, source/sink ratio and regulation of developmental timing, compared to otherwise phenotypically similar accessions. Bna.MAP3K13.C05 expression was found to be protective independently of the tolerance mechanism, indicating cross-talk to nitrogen signaling. Moreover, we identified putative miRNA genes in the B. napus genome which respond to stress and may also be involved in protective mechanisms, representing possible breeding targets.

Melatonin improves the germination rate of cotton seeds under drought stress by opening pores in the seed coat

The germination of cotton (Gossypium hirsutum L.) seeds is affected by drought stress; however, little is known about the physiological mechanism affecting germination and the effect of melatonin (MT) on cotton seed germination under drought stress. Therefore, we studied the effects of exogenous MT on the antioxidant capacity and epidermal microstructure of cotton under drought stress. The results demonstrated a retarded water absorption capacity of testa under drought stress, significantly inhibiting germination and growth in cotton seeds. Drought stress led to the accumulation of reactive oxygen species (ROS), malondialdehyde (MDA), and osmoregulatory substances (e.g., proline, soluble protein, and soluble sugars); it also decreased the activity of antioxidant enzymes and α-amylase. Drought stress inhibited gibberellin acid (GA3) synthesis and increased abscisic acid (ABA) content, seriously affecting seed germination. However, seeds pre-soaked with MT (100 µM) showed a positive regulation in the number and opening of stomata in cotton testa. The exogenous application of MT increased the germination rate, germination potential, radical length, and fresh weight, as well as the activities of superoxide dismutase (SOD), peroxidase (POD), and α-amylase. In addition, MT application increased the contents of organic osmotic substances by decreasing the hydrogen peroxide (H2O2), superoxide anion (O2 -), and MDA levels under drought stress. Further analysis demonstrated that seeds pre-soaked with MT alleviated drought stress by affecting the ABA and GA3 contents. Our findings show that MT plays a positive role in protecting cotton seeds from drought stress.

Differentially evolved drought stress indices determine the genetic variation of Brassica napus at seedling traits by genome-wide association mapping

Drought seriously curtails growth, physiology and productivity in rapeseed (Brassica napus). Although drought tolerance is a complex trait, efficient phenotyping and genotyping has led to the identification of novel marker-trait associations underlying drought tolerance. A diverse panel of 228 Brassica accessions was phenotyped under normal (without stress) and water-stress conditions, simulated by polyethylene glycol (PEG-6000) (15% PEG stress) at the seedling stage; stress tolerance index (STI) and stress susceptibility index (SSI) values were acquired. Genome-wide association studies (GWAS) using 201 817 high quality SNPs identified 314 marker-trait associations strongly linked with drought indices and distributed across all nineteen chromosomes in both the A and C genomes. None of these quantitative trait loci (QTL) had been previously identified by other studies. We identified 85 genes underlying these QTL (most within 100 kb of associated SNPs) which were orthologous to Arabidopsis genes known to be associated with drought tolerance. Our study provides a novel resource for breeding drought-tolerant Brassica crops.

Root transcriptome profiling of contrasting wheat genotypes provides an insight to their adaptive strategies to water deficit

Water deficit limits plant growth and productivity in wheat. The effect of water deficit varies considerably in the contrasting genotypes. This study attempted comparative transcriptome profiling of the tolerant (Abura) and susceptible (AUS12671) genotypes under PEG-simulated water stress via genome-wide RNA-seq technology to understand the dynamics of tolerance mechanism. Morphological and physiological analyses indicated that the tolerant genotype Abura had a higher root growth and net photosynthesis, which accounted for its higher root biomass than AUS12671 under stress. Transcriptomic analysis revealed a total of 924 differentially expressed genes (DEGs) that were unique in the contrasting genotypes under stress across time points. The susceptible genotype AUS12671 had slightly more abundant DEGs (505) than the tolerant genotype Abura (419). Gene ontology enrichment and pathway analyses of these DEGs suggested that the two genotypes differed significantly in terms of adaptive mechanism. Predominant upregulation of genes involved in various metabolic pathways was the key adaptive feature of the susceptive genotype AUS12671 indicating its energy-consuming approach in adaptation to water deficit. In contrast, downregulation the expression of genes of key pathways, such as global and overview maps, carbohydrate metabolism, and genetic information processing was the main strategy for the tolerant genotype Abura. Besides, significantly higher number of genes encoding transcription factors (TF) families like MYB and NAC, which were reported to be associated with stress defense, were differentially expressed in the tolerant genotype Abura. Gene encoding transcription factors TIFY were only differentially expressed between stressed and non-stressed conditions in the sensitive genotype. The identified DEGs and the suggested differential adaptive strategies of the contrasting genotypes provided an insight for improving water deficit tolerance in wheat.

Genome-wide association studies for inflorescence type and remontancy in Hydrangea macrophylla

Inflorescence type and remontancy are two valuable traits in bigleaf hydrangea (Hydrangea macrophylla L.) and both are recessively inherited. Molecular marker-assisted selection (MAS) can greatly reduce the time necessary to breed cultivars with desired traits. In this study, a genome-wide association study (GWAS) using 5803 single-nucleotide polymorphisms (SNPs) was performed using a panel of 82 bigleaf hydrangea cultivars. One SNP locus (Hy_CAPS_Inflo) associated with inflorescence type was identified with general linear model (GLM) and mixed linear model (MLM) methods that explained 65.5% and 36.1% of the phenotypic variations, respectively. Twenty-three SNPs associated with remontancy were detected in GLM whereas no SNP was detected in MLM. The SNP locus (Hy_CAPS_Inflo) was converted to a cleaved amplified polymorphic sequence (CAPS) marker that showed absolute identification accuracy (100%) of inflorescence type in a validation panel consisting of eighteen H. macrophylla cultivars. The SNP was investigated in 341 F1 progenies using genotyping by sequencing (GBS) and co-segregated with inflorescence type (χ 2 = 0.12; P = 0.73). The SNP was subsequently used for breeding selection using kompetitive allele specific PCR (KASP) technology. Future directions for the use of genomics and MAS in hydrangea breeding improvement are discussed. The results presented in this study provide insights for further research on understanding genetic mechanisms behind inflorescence type and remontancy in H. macrophylla. The CAPS and KASP markers developed here will be immediately useful for applying MAS to accelerate breeding improvement in hydrangea.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering

BACKGROUND

Drought stress is a major abiotic factor that affects rapeseed (Brassica napus L.) productivity. Though previous studies indicated that long non-coding RNAs (lncRNAs) play a key role in response to drought stress, a scheme for genome-wide identification and characterization of lncRNAs' response to drought stress is still lacking, especially in the case of B. napus. In order to further understand the molecular mechanism of the response of B. napus to drought stress, we compared changes in the transcriptome between Q2 (a drought-tolerant genotype) and Qinyou8 (a drought-sensitive genotype) responding drought stress and rehydration treatment at the seedling stage.

RESULTS

A total of 5546 down-regulated and 6997 up-regulated mRNAs were detected in Q2 compared with 7824 and 10,251 in Qinyou8, respectively; 369 down-regulated and 108 up- regulated lncRNAs were detected in Q2 compared with 449 and 257 in Qinyou8, respectively. LncRNA-mRNA interaction network analysis indicated that the co-expression network of Q2 was composed of 145 network nodes and 5175 connections, while the co-expression network of Qinyou8 was composed of 305 network nodes and 22,327 connections. We further identified 34 transcription factors (TFs) corresponding to 126 differentially expressed lncRNAs in Q2, and 45 TFs corresponding to 359 differentially expressed lncRNAs in Qinyou8. Differential expression analysis of lncRNAs indicated that up- and down-regulated mRNAs co-expressed with lncRNAs participated in different metabolic pathways and were involved in different regulatory mechanisms in the two genotypes. Notably, some lncRNAs were co-expressed with BnaC07g44670D, which are associated with plant hormone signal transduction. Additionally, some mRNAs co-located with XLOC_052298, XLOC_094954 and XLOC_012868 were mainly categorized as signal transport and defense/stress response.

CONCLUSIONS

The results of this study increased our understanding of expression characterization of rapeseed lncRNAs in response to drought stress and re-watering, which would be useful to provide a reference for the further study of the function and action mechanisms of lncRNAs under drought stress and re-watering.

High Nitrogen Enhance Drought Tolerance in Cotton through Antioxidant Enzymatic Activities, Nitrogen Metabolism and Osmotic Adjustment

Drought is one of the most important abiotic stresses and hampers many plant physiological processes under suboptimal nitrogen (N) concentration. Seedling tolerance to drought stress is very important for optimum growth and development, however, the enhancement of plant stress tolerance through N application in cotton is not fully understood. Therefore, this study investigates the role of high N concentration in enhancing drought stress tolerance in cotton. A hydroponic experiment supplying low (0.25 mM) and high (5 mM) N concentrations, followed by 150 g L-1 polyethylene glycol (PEG)-induced stress was conducted in a growth chamber. PEG-induced drought stress inhibited seedling growth, led to oxidative stress from excessive malondialdehyde (MDA) generation, and reduced N metabolism. High N concentrations alleviated oxidative damage and stomatal limitation by increasing antioxidant enzymatic activities, leaf relative water content, and photosynthesis in cotton seedlings under drought stress. The results revealed that the ameliorative effects of high N concentration may be ascribed to the enhancement of N metabolizing enzymes and an increase in the amounts of osmoprotectants like free amino acids and total soluble protein. The present data suggest that relatively high N concentrations may contribute to drought stress tolerance in cotton through N metabolism, antioxidant capacity, and osmotic adjustment.

Transcriptome, physiological and biochemical analysis of Triarrhena sacchariflora in response to flooding stress

Background

In recent decades, the frequency of flooding is increasing with the change of global climate. Flooding has become one of the major abiotic stresses that seriously affect growth and development of plants. Triarrhena sacchariflora Nakai has been considered a promising energy crop for utilization in ethanol production. Flooding stress is among the most severe abiotic stressors in the production of Nakai. However, the physiological and molecular biological mechanisms of Nakai response to flooding is still unclear. In the present study, in order to understand the molecular mechanisms of Nakai in response to flooding stress, the transcriptome, physiological and biochemical were investigated.

Results

The results demonstrated that significant physiological changes were observed in photosynthetic system, antioxidative enzyme activity, chlorophyll, carotenoid, proline, lipid peroxidation and soluble sugar content under normal and flooding treatments. Such as, the chlorophyll, carotenoid contents and photosynthetic system were significantly decreased. Whereas, the antioxidative enzyme activity, proline, lipid peroxidation and soluble sugar has increased first and then decreased under treatments compared with the normal plants. Additionally, a total of 8832, 6608 and 3649 unigenes were validated to be differentially expressed under different treatments, respectively. Besides, gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of the different expression levels of genes also presented processes, which involved in photosynthesis, sucrose catabolism, glycolysis, stress response and defense, phytohormone biosynthesis and signal transduction.

Conclusions

The results provide a comprehensive view of the complex molecular events involved in the response to flooding stress of Nakai leaves, which also will promote the research in the development of flood-resistant crops and provide new tools for Nakai breeders.

Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus

BACKGROUND

Increasing seed oil content is one of the most important targets for rapeseed (Brassica napus) breeding. However, genetic mechanisms of mature seed oil content in Brassica napus (B. napus) remain little known. To identify oil content-related genes, a genome-wide association study (GWAS) was performed using 588 accessions.

RESULTS

High-throughput genome resequencing resulted in 385,692 high-quality single nucleotide polymorphism (SNPs) with a minor allele frequency (MAF) > 0.05. We identified 17 loci that were significantly associated with seed oil content, among which 12 SNPs were distributed on the A3 (11 loci) and A1 (one loci) chromosomes, and five novel significant SNPs on the C5 (one loci) and C7 (four loci) chromosomes, respectively. Subsequently, we characterized differentially expressed genes (DEGs) between the seeds and silique pericarps on main florescences and primary branches of extremely high- and low-oil content accessions (HO and LO). A total of 64 lipid metabolism-related DEGs were identified, 14 of which are involved in triacylglycerols (TAGs) biosynthesis and assembly. Additionally, we analyzed differences in transcription levels of key genes involved in de novo fatty acid biosynthesis in the plastid, TAGs assembly and lipid droplet packaging in the endoplasmic reticulum (ER) between high- and low-oil content B. napus accessions.

CONCLUSIONS

The combination of GWAS and transcriptome analyses revealed seven candidate genes located within the confidence intervals of significant SNPs. Current findings provide valuable information for facilitating marker-based breeding for higher seed oil content in B. napus.

Transcriptome Profile Analysis of Winter Rapeseed (Brassica napus L.) in Response to Freezing Stress, Reveal Potentially Connected Events to Freezing Stress

Winter rapeseed is not only an important oilseed crop, but also a winter cover crop in Northern China, where its production was severely limited by freezing stress. As an overwinter crop, the production is severely limited by freezing stress. Therefore, understanding the physiological and molecular mechanism of winter rapeseed (Brassica napus L.) in freezing stress responses becomes essential for the improvement and development of freezing-tolerant varieties of Brassica napus. In this study, morphological, physiological, ultrastructure and transcriptome changes in the Brassica napus line "2016TS(G)10" (freezing-tolerance line) that was exposed to -2 °C for 0 h, 1 h, 3 h and 24 h were characterized. The results showed that freezing stress caused seedling dehydration, and chloroplast dilation and degradation. The content of malondialdehyde (MDA), proline, soluble protein and soluble sugars were increased, as well as the relative electrolyte leakage (REL) which was significantly increased at frozen 24 h. Subsequently, RNA-seq analysis revealed a total of 98,672 UniGenes that were annotated in Brassica napus and 3905 UniGenes were identified as differentially expressed genes after being exposed to freezing stress. Among these genes, 2312 (59.21%) were up-regulated and 1593 (40.79%) were down-regulated. Most of these DEGs were significantly annotated in the carbohydrates and energy metabolism, signal transduction, amino acid metabolism and translation. Most of the up-regulated DEGs were especially enriched in plant hormone signal transduction, starch and sucrose metabolism pathways. Transcription factor enrichment analysis showed that the AP2/ERF, WRKY and MYB families were also significantly changed. Furthermore, 20 DEGs were selected to validate the transcriptome profiles via quantitative real-time PCR (qRT-PCR). In conclusion, the results provide an overall view of the dynamic changes in physiology and insights into the molecular regulation mechanisms of winter Brassica napus in response to freezing treatment, expanding our understanding on the complex molecular mechanism in plant response to freezing stress.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

Associating transcriptional regulation for rapid germination of rapeseed (Brassica napus L.) under low temperature stress through weighted gene co-expression network analysis

Slow germination speed caused by low temperature stress intensifies the risk posed by adverse environmental factors, contributing to low germination rate and reduced production of rapeseed. The purpose of this study was to understand the transcriptional regulation mechanism for rapid germination of rapeseed. The results showed that seed components and size do not determine the seed germination speed. Different temporal transcriptomic profiles were generated under normal and low temperature conditions in genotypes with fast and slow germination speeds. Using weight gene co-expression network analysis, 37 823 genes were clustered into 15 modules with different expression patterns. There were 10 233 and 9111 differentially expressed genes found to follow persistent tendency of up- and down-regulation, respectively, which provided the conditions necessary for germination. Hub genes in the continuous up-regulation module were associated with phytohormone regulation, signal transduction, the pentose phosphate pathway, and lipolytic metabolism. Hub genes in the continuous down-regulation module were involved in ubiquitin-mediated proteolysis. Through pairwise comparisons, 1551 specific upregulated DEGs were identified for the fast germination speed genotype under low temperature stress. These DEGs were mainly enriched in RNA synthesis and degradation metabolisms, signal transduction, and defense systems. Transcription factors, including WRKY, bZIP, EFR, MYB, B3, DREB, NAC, and ERF, are associated with low temperature stress in the fast germination genotype. The aquaporin NIP5 and late embryogenesis abundant (LEA) protein genes contributed to the water uptake and transport under low temperature stress during seed germination. The ethylene/H₂O₂-mediated signal pathway plays an important role in cell wall loosening and embryo extension during germination. The ROS-scavenging system, including catalase, aldehyde dehydrogenase, and glutathione S-transferase, was also upregulated to alleviate ROS toxicity in the fast germinating genotype under low temperature stress. These findings should be useful for molecular assisted screening and breeding of fast germination speed genotypes for rapeseed.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Multiple genomic regions influence root morphology and seedling growth in cultivated sunflower (Helianthus annuus L.) under well-watered and water-limited conditions

With climate change and an ever-increasing human population threatening food security, developing a better understanding of the genetic basis of crop performance under stressful conditions has become increasingly important. Here, we used genome-wide association studies to genetically dissect variation in seedling growth traits in cultivated sunflower (Helianthus annuus L.) under well-watered and water-limited (i.e., osmotic stress) conditions, with a particular focus on root morphology. Water limitation reduced seedling size and produced a shift toward deeper rooting. These effects varied across genotypes, and we identified 13 genomic regions that were associated with traits of interest across the two environments. These regions varied in size from a single marker to 186.2 Mbp and harbored numerous genes, some of which are known to be involved in the plant growth/development as well as the response to osmotic stress. In many cases, these associations corresponded to growth traits where the common allele outperformed the rare variant, suggesting that selection for increased vigor during the evolution of cultivated sunflower might be responsible for the relatively high frequency of these alleles. We also found evidence of pleiotropy across multiple traits, as well as numerous environmentally independent genetic effects. Overall, our results indicate the existence of genetic variation in root morphology and allocation and further suggest that the majority of alleles associated with these traits have consistent effects across environments.

Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L

BACKGROUND

Brassica napus (B. napus) encompasses diverse transcription factors (TFs), but thorough identification and characterization of TF families, as well as their transcriptional responsiveness to multifarious stresses are still not clear.

RESULTS

Totally 2167 TFs belonging to five families were genome-widely identified in B. napus, including 518 BnAP2/EREBPs, 252 BnbZIPs, 721 BnMYBs, 398 BnNACs and 278 BnWRKYs, which contained some novel members in comparison with existing results. Sub-genome distributions of BnAP2/EREBPs and BnMYBs indicated that the two families might have suffered from duplication and divergence during evolution. Synteny analysis revealed strong co-linearity between B. napus and its two ancestors, although chromosomal rearrangements have occurred and 85 TFs were lost. About 7.6% and 9.4% TFs of the five families in B. napus were novel genes and conserved genes, which both showed preference on the C sub-genome. RNA-Seq revealed that more than 80% TFs were abiotic stress inducible and 315 crucial differentially expressed genes (DEGs) were screened out. Network analysis revealed that the 315 DEGs are highly co-expressed. The homologous gene network in A. thaliana revealed that a considerable amount of TFs could trigger the differential expression of targeted genes, resulting in a complex clustered network with clusters of genes responsible for targeted stress responsiveness.

CONCLUSIONS

We identified and characterized five TF families in B. napus. Some crucial members and regulatory networks involved in different abiotic stresses have been explored. The investigations deepen our understanding of TFs for stress tolerance in B. napus.

Mitochondrial Biogenesis in Diverse Cauliflower Cultivars under Mild and Severe Drought. Impaired Coordination of Selected Transcript and Proteomic Responses, and Regulation of Various Multifunctional Proteins

Mitochondrial responses under drought within Brassica genus are poorly understood. The main goal of this study was to investigate mitochondrial biogenesis of three cauliflower (Brassica oleracea var. botrytis) cultivars with varying drought tolerance. Diverse quantitative changes (decreases in abundance mostly) in the mitochondrial proteome were assessed by two-dimensional gel electrophoresis (2D PAGE) coupled with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Respiratory (e.g., complex II, IV (CII, CIV) and ATP synthase subunits), transporter (including diverse porin isoforms) and matrix multifunctional proteins (e.g., components of RNA editing machinery) were diversely affected in their abundance under two drought levels. Western immunoassays showed additional cultivar-specific responses of selected mitochondrial proteins. Dehydrin-related tryptic peptides (found in several 2D spots) immunopositive with dehydrin-specific antisera highlighted the relevance of mitochondrial dehydrin-like proteins for the drought response. The abundance of selected mRNAs participating in drought response was also determined. We conclude that mitochondrial biogenesis was strongly, but diversely affected in various cauliflower cultivars, and associated with drought tolerance at the proteomic and functional levels. However, discussed alternative oxidase (AOX) regulation at the RNA and protein level were largely uncoordinated due to the altered availability of transcripts for translation, mRNA/ribosome interactions, and/or miRNA impact on transcript abundance and translation.

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

A user guide to the Brassica 60K Illumina Infinium™ SNP genotyping array

The Brassica napus 60K Illumina Infinium™ SNP array has had huge international uptake in the rapeseed community due to the revolutionary speed of acquisition and ease of analysis of this high-throughput genotyping data, particularly when coupled with the newly available reference genome sequence. However, further utilization of this valuable resource can be optimized by better understanding the promises and pitfalls of SNP arrays. We outline how best to analyze Brassica SNP marker array data for diverse applications, including linkage and association mapping, genetic diversity and genomic introgression studies. We present data on which SNPs are locus-specific in winter, semi-winter and spring B. napus germplasm pools, rather than amplifying both an A-genome and a C-genome locus or multiple loci. Common issues that arise when analyzing array data will be discussed, particularly those unique to SNP markers and how to deal with these for practical applications in Brassica breeding applications.

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.

Transcriptomic basis for drought-resistance in Brassica napus L

Based on transcriptomic data from four experimental settings with drought-resistant and drought-sensitive cultivars under drought and well-watered conditions, statistical analysis revealed three categories encompassing 169 highly differentially expressed genes (DEGs) in response to drought in Brassica napus L., including 37 drought-resistant cultivar-related genes, 35 drought-sensitive cultivar-related genes and 97 cultivar non-specific ones. We provide evidence that the identified DEGs were fairly uniformly distributed on different chromosomes and their expression patterns are variety specific. Except commonly enriched in response to various stimuli or stresses, different categories of DEGs show specific enrichment in certain biological processes or pathways, which indicated the possibility of functional differences among the three categories. Network analysis revealed relationships among the 169 DEGs, annotated biological processes and pathways. The 169 DEGs can be classified into different functional categories via preferred pathways or biological processes. Some pathways might simultaneously involve a large number of shared DEGs, and these pathways are likely to cross-talk and have overlapping biological functions. Several members of the identified DEGs fit to drought stress signal transduction pathway in Arabidopsis thaliana. Finally, quantitative real-time PCR validations confirmed the reproducibility of the RNA-seq data. These investigations are profitable for the improvement of crop varieties through transgenic engineering.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.