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

Combined BSA-Seq and RNA-Seq Analysis to Identify Candidate Genes Associated with Aluminum Toxicity in Rapeseed (Brassica napus L.)

Exchangeable aluminum (Al) ions released from acidic soils with pH < 5.5 inhibit root elongation of crops, ultimately leading to yield reduced. It is necessary to identify the quantitative trait locus (QTLs) and candidate genes that confer toxicity resistance to understand the mechanism and improve tolerance of rapeseed. In this study, an F2 segregating population was derived from a cross between Al-tolerance inbred line FDH188 (R178) and -sensitive inbred line FDH152 (S169), and the F2:3 were used as materials to map QTLs associated with the relative elongation of taproot (RET) under Al toxicity stress. Based on bulked segregant analysis sequencing (BSA-seq), three QTLs (qAT-A07-1, qAT-A07-2, and qAT-A09-1) were detected as significantly associated with RET, and 656 candidate genes were screened. By combined BSA and RNA-seq analysis, 55 candidate genes showed differentially expressed, including genes encoding ABC transporter G (ABCG), zinc finger protein, NAC, ethylene-responsive transcription factor (ERF), etc. These genes were probably positive factors in coping with Al toxicity stress in rapeseed. This study provides new insight into exploring the QTLs and candidate genes' response to Al toxicity stress by combined BSA-seq and RNA-seq and is helpful to further research on the mechanism of Al resistance in rapeseed.

Identification of Genomic Regions for Traits Associated with Flowering in Cassava (Manihot esculenta Crantz)

Flowering in cassava (Manihot esculenta Crantz) is crucial for the generation of botanical seed for breeding. However, genotypes preferred by most farmers are erect and poor at flowering or never flower. To elucidate the genetic basis of flowering, 293 diverse cassava accessions were evaluated for flowering-associated traits at two locations and seasons in Uganda. Genotyping using the Diversity Array Technology Pty Ltd. (DArTseq) platform identified 24,040 single-nucleotide polymorphisms (SNPs) distributed on the 18 cassava chromosomes. Population structure analysis using principal components (PCs) and kinships showed three clusters; the first five PCs accounted for 49.2% of the observed genetic variation. Linkage disequilibrium (LD) estimation averaged 0.32 at a distance of ~2850 kb (kilo base pairs). Polymorphism information content (PIC) and minor allele frequency (MAF) were 0.25 and 0.23, respectively. A genome-wide association study (GWAS) analysis uncovered 53 significant marker-trait associations (MTAs) with flowering-associated traits involving 27 loci. Two loci, SNPs S5_29309724 and S15_11747301, were associated with all the traits. Using five of the 27 SNPs with a Phenotype_Variance_Explained (PVE) ≥ 5%, 44 candidate genes were identified in the peak SNP sites located within 50 kb upstream or downstream, with most associated with branching traits. Eight of the genes, orthologous to Arabidopsis and other plant species, had known functional annotations related to flowering, e.g., eukaryotic translation initiation factor and myb family transcription factor. This study identified genomic regions associated with flowering-associated traits in cassava, and the identified SNPs can be useful in marker-assisted selection to overcome hybridization challenges, like unsynchronized flowering, and candidate gene validation.

Genome-wide identification and expression analysis of phenylalanine ammonia-lyase (PAL) family in rapeseed (Brassica napus L.)

BACKGROUND

Phenylalanine ammonia-lyase (PAL), as a key enzyme in the phenylalanine metabolism pathway in plants, plays an important role in the response to environmental stress. However, the PAL family responding to abiotic stress has not been fully characterized in rapeseed.

RESULTS

In this study, we conducted a genome-wide study of PAL family, and analyzed their gene structure, gene duplication, conserved motifs, cis-acting elements and response to stress treatment. A total of 17 PALs were identified in the rapeseed genome. Based on phylogenetic analysis, the BnPALs were divided into four clades (I, II, IV, and V). The prediction of protein structure domain presented that all BnPAL members contained a conservative PAL domain. Promoter sequence analysis showed that the BnPALs contain many cis-acting elements related to hormone and stress responses, indicating that BnPALs are widely involved in various biological regulatory processes. The expression profile showed that the BnPALs were significantly induced under different stress treatments (NaCl, Na2CO3, AlCl3, and PEG), suggesting that BnPAL family played an important role in response to abiotic stress.

CONCLUSIONS

Taken together, our research results comprehensively characterized the BnPAL family, and provided a valuable reference for revealing the role of BnPALs in the regulation of abiotic stress responses in rapeseed.

Integration of GWAS and RNA-Seq Analysis to Identify SNPs and Candidate Genes Associated with Alkali Stress Tolerance at the Germination Stage in Mung Bean

Soil salt-alkalization seriously impacts crop growth and productivity worldwide. Breeding and applying tolerant varieties is the most economical and effective way to address soil alkalization. However, genetic resources for breeders to improve alkali tolerance are limited in mung bean. Here, a genome-wide association study (GWAS) was performed to detect alkali-tolerant genetic loci and candidate genes in 277 mung bean accessions during germination. Using the relative values of two germination traits, 19 QTLs containing 32 SNPs significantly associated with alkali tolerance on nine chromosomes were identified, and they explained 3.6 to 14.6% of the phenotypic variance. Moreover, 691 candidate genes were mined within the LD intervals containing significant trait-associated SNPs. Transcriptome sequencing of alkali-tolerant accession 132-346 under alkali and control conditions after 24 h of treatment was conducted, and 2565 DEGs were identified. An integrated analysis of the GWAS and DEGs revealed six hub genes involved in alkali tolerance responses. Moreover, the expression of hub genes was further validated by qRT-PCR. These findings improve our understanding of the molecular mechanism of alkali stress tolerance and provide potential resources (SNPs and genes) for the genetic improvement of alkali tolerance in mung bean.

Genome-Wide Association Studies of Salt Tolerance at the Seed Germination Stage and Yield-Related Traits in Brassica napus L

Salt stress severely affects crop growth and development and reduces the yield of Brassica napus. Exploring natural genetic variations for high salt tolerance in B. napus seedlings is an effective approach to improve productivity under salt stress. Using 10,658 high-quality single nucleotide polymorphic (SNP) markers developed by specific-locus amplified fragment sequencing (SLAF-seq) technology, genome-wide association studies (GWAS) were performed to investigate the genetic basis of salt tolerance and yield-related traits of B. napus. The results revealed that 77 and 497 SNPs were significantly associated with salt tolerance and yield-related traits, of which 40 and 58 SNPs were located in previously reported QTLs/SNPs, respectively. We identified nineteen candidate genes orthologous with Arabidopsis genes known to be associated with salt tolerance and seven potential candidates controlling both salt tolerance and yield. Our study provides a novel genetic resource for the breeding of high-yield cultivars resistant to salt stress.

Finger millet RNA-seq reveals differential gene expression associated with tolerance to aluminum toxicity and provides novel genomic resources

Eleusine coracana, finger millet, is a multipurpose crop cultivated in arid and semi-arid regions of Africa and Asia. RNA sequencing (RNA-seq) was used in this study to obtain valuable genomic resources and identify genes differentially expressed between Al-tolerant and Al-susceptible genotypes. Two groups of finger millet genotypes were used: Al-tolerant (215836, 215845, and 229722) and Al-susceptible (212462, 215804 and 238323). The analysis of the RNA-seq data resulted in 198,546 unigenes, 56.5% of which were annotated with significant hits in one or more of the following six databases: NR (48.8%), GO (29.7%), KEGG (45%), PlantTFDB (19.0%), Uniprot (49.2%), and NT (46.2%). It is noteworthy that only 220 unigenes in the NR database had significant hits against finger millet sequences suggesting that finger millet's genomic resources are scarce. The gene expression analysis revealed that 322 genes were significantly differentially expressed between the Al-tolerant and Al-susceptible genotypes, of which 40.7% were upregulated while 59.3% were downregulated in Al-tolerant genotypes. Among the significant DEGs, 54.7% were annotated in the GO database with the top hits being ATP binding (GO:0005524) and DNA binding (GO:0003677) in the molecular function, DNA integration (GO:0015074) and cell redox homeostasis in the biological process, as well as cellular anatomical entity and intracellular component in the cellular component GO classes. Several of the annotated DEGs were significantly enriched for their corresponding GO terms. The KEGG pathway analysis resulted in 60 DEGs that were annotated with different pathway classes, of which carbohydrate metabolism and signal transduction were the most prominent. The homologs of a number of significant DEGs have been previously reported as being associated with Al or other abiotic stress responses in various crops, including carboxypeptidase SOL1, HMA3, AP2, bZIP, C3H, and WRKY TF genes. A more detailed investigation of these and other DEGs will enable genomic-led breeding for Al tolerance in finger millet. RNA-seq data analysis also yielded 119,073 SNP markers, the majority of which had PIC values above 0.3, indicating that they are highly informative. Additionally, 3,553 single-copy SSR markers were identified, of which trinucleotide SSRs were the most prevalent. These genomic resources contribute substantially to the enrichment of genomic databases for finger millet, and facilitate future research on this crop.

Reprisal of Schima superba to Mn stress and exploration of its defense mechanism through transcriptomic analysis

One of the most diverse protein families, ATP-binding cassette (ABC) transporters, play a role in disease resistance, heavy metal tolerance, and food absorption.Differentially expressed genes contribute in the investigation of plant defense mechanisms under varying stress conditions. To elucidate the molecular mechanisms involved in Mn metal stress, we performed a transcriptomic analysis to explore the differential gene expression in Schima superba with the comparison of control. A total of 79.84 G clean data was generated and 6558 DEGs were identified in response to Mn metal stress. Differentially expressed genes were found to be involved in defense, signaling pathways, oxidative burst, transcription factors and stress responses. Genes important in metal transport were more expressive in Mn stress than control plants. The investigation of cis-acting regions in the ABC family indicated that these genes might be targeted by a large variety of trans-acting elements to control a variety of stress circumstances. Moreover, genes involved in defense responses, the mitogen-activated protein kinase (MAPK) signaling and signal transduction in S. superba were highly induced in Mn stress. Twenty ABC transporters were variably expressed on 1st, 5th, and 10th day of Mn treatment, according to the qRT PCR data. Inclusively, our findings provide an indispensable foundation for an advanced understanding of the metal resistance mechanisms. Our study will enrich the sequence information of S. superba in a public database and would provide a new understanding of the molecular mechanisms of heavy metal tolerance and detoxification.