Genome wide association analyses to understand genetic basis of flowering and plant height under three levels of nitrogen application in Brassica juncea (L.) Czern & Coss

Physiology, genomics, and evolutionary aspects of desert plants

BACKGROUND

Despite the exposure to arid environmental conditions across the globe ultimately hampering the sustainability of the living organism, few plant species are equipped with several unique genotypic, biochemical, and physiological features to counter such harsh conditions. Physiologically, they have evolved with reduced leaf size, spines, waxy cuticles, thick leaves, succulent hydrenchyma, sclerophyll, chloroembryo, and photosynthesis in nonfoliar and other parts. At the biochemical level, they are evolved to perform efficient photosynthesis through Crassulacean acid metabolism (CAM) and C4 pathways with the formation of oxaloacetic acid (Hatch-Slack pathway) instead of the C3 pathway. Additionally, comparative genomics with existing data provides ample evidence of the xerophytic plants' positive selection to adapt to the arid environment. However, adding more high-throughput sequencing of xerophyte plant species is further required for a comparative genomic study toward trait discovery related to survival. Learning from the mechanism to survive in harsh conditions could pave the way to engineer crops for future sustainable agriculture.

AIM OF THE REVIEW

The distinct physiology of desert plants allows them to survive in harsh environments. However, the genomic composition also contributes significantly to this and requires great attention. This review emphasizes the physiological and genomic adaptation of desert plants. Other important parameters, such as desert biodiversity and photosynthetic strategy, are also discussed with recent progress in the field. Overall, this review discusses the different features of desert plants, which prepares them for harsh conditions intending to translate knowledge to engineer plant species for sustainable agriculture.

KEY SCIENTIFIC CONCEPTS OF REVIEW

This review comprehensively presents the physiology, molecular mechanism, and genomics of desert plants aimed towards engineering a sustainable crop.

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.

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.

Identifying SNP markers associated with distinctness, uniformity, and stability testing in Egyptian fenugreek genotypes

Distinctness, uniformity, and stability (DUS) test is the legal requirement in crop breeding to grant the intellectual property right for new varieties by evaluating their morphological characteristics across environments. On the other hand, molecular markers accurately identify genetic variations and validate the purity of the cultivars. Therefore, genomic DUS can improve the efficiency of traditional DUS testing. In this study, 112 Egyptian fenugreek genotypes were grown in Egypt at two locations: Wadi El-Natrun (Wadi), El-Beheira Governorate, with salty and sandy soil, and Giza, Giza governorate, with loamy clay soil. Twelve traits were measured, of which four showed a high correlation above 0.94 over the two locations. We observed significant genotype-by-location interactions (GxL) for seed yield, as it was superior in Wadi, with few overlapping genotypes with Giza. We attribute this superiority in Wadi to the maternal habitat, as most genotypes grew in governorates with newly reclaimed salty and sandy soil. As a first step toward genomic DUS, we performed an association study, and out of 38,142 SNPs, we identified 39 SNPs demonstrating conditional neutrality and four showing pleiotropic effects. Forty additional SNPs overlapped between both locations, each showing a similar impact on the associated trait. Our findings highlight the importance of GxL in validating the effect of each SNP to make better decisions about its suitability in the marker-assisted breeding program and demonstrate its potential use in registering new plant varieties.

Mapping and Functional Analysis of QTL for Kernel Number per Row in Tropical and Temperate-Tropical Introgression Lines of Maize (Zea mays L.)

Kernel number per row (KNR) is an essential component of maize (Zea mays L.) grain yield (GY), and understanding its genetic mechanism is crucial to improve GY. In this study, two F₇ recombinant inbred line (RIL) populations were created using a temperate-tropical introgression line TML418 and a tropical inbred line CML312 as female parents and a backbone maize inbred line Ye107 as the common male parent. Bi-parental quantitative trait locus (QTL) mapping and genome-wide association analysis (GWAS) were then performed on 399 lines of the two maize RIL populations for KNR in two different environments using 4118 validated single nucleotide polymorphism (SNP) markers. This study aimed to: (1) detect molecular markers and/or the genomic regions associated with KNR; (2) identify the candidate genes controlling KNR; and (3) analyze whether the candidate genes are useful in improving GY. The authors reported a total of 7 QTLs tightly linked to KNR through bi-parental QTL mapping and identified 21 SNPs significantly associated with KNR through GWAS. Among these, a highly confident locus qKNR7-1 was detected at two locations, Dehong and Baoshan, with both mapping approaches. At this locus, three novel candidate genes (Zm00001d022202, Zm00001d022168, Zm00001d022169) were identified to be associated with KNR. These candidate genes were primarily involved in the processes related to compound metabolism, biosynthesis, protein modification, degradation, and denaturation, all of which were related to the inflorescence development affecting KNR. These three candidate genes were not reported previously and are considered new candidate genes for KNR. The progeny of the hybrid Ye107 × TML418 exhibited strong heterosis for KNR, which the authors believe might be related to qKNR7-1. This study provides a theoretical foundation for future research on the genetic mechanism underlying KNR in maize and the use of heterotic patterns to develop high-yielding hybrids.

Advancing designer crops for climate resilience through an integrated genomics approach

Climate change and exponential population growth are exposing an immediate need for developing future crops that are highly resilient and adaptable to changing environments to maintain global food security in the next decade. Rigorous selection from long domestication history has rendered cultivated crops genetically disadvantaged, raising concerns in their ability to adapt to these new challenges and limiting their usefulness in breeding programmes. As a result, future crop improvement efforts must rely on integrating various genomic strategies ranging from high-throughput sequencing to machine learning, in order to exploit germplasm diversity and overcome bottlenecks created by domestication, expansive multi-dimensional phenotypes, arduous breeding processes, complex traits and big data.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.)

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.

Genome-Wide Comparative Analysis of Flowering-Time Genes; Insights on the Gene Family Expansion and Evolutionary Perspective

In polyploids, whole genome duplication (WGD) played a significant role in genome expansion, evolution and diversification. Many gene families are expanded following polyploidization, with the duplicated genes functionally diversified by neofunctionalization or subfunctionalization. These mechanisms may support adaptation and have likely contributed plant survival during evolution. Flowering time is an important trait in plants, which affects critical features, such as crop yields. The flowering-time gene family is one of the largest expanded gene families in plants, with its members playing various roles in plant development. Here, we performed genome-wide identification and comparative analysis of flowering-time genes in three palnt families i.e., Malvaceae, Brassicaceae, and Solanaceae, which indicate these genes were expanded following the event/s of polyploidization. Duplicated genes have been retained during evolution, although genome reorganization occurred in their flanking regions. Further investigation of sequence conservation and similarity network analyses provide evidence for functional diversification of duplicated genes during evolution. These functionally diversified genes play important roles in plant development and provide advantages to plants for adaptation and survival in response to environmental changes encountered during evolution. Collectively, we show that flowering-time genes were expanded following polyploidization and retained as large gene family by providing advantages from functional diversification during evolution.

Unraveling the Relationship Between Seed Yield and Yield-Related Traits in a Diversity Panel of Brassica juncea Using Multi-Traits Mixed Model

The response to selection in any crop improvement program depends on the degree of variance and heritability. The objective of the current study was to explain variance and heritability components in Indian mustard Brassica juncea (L). Czern & Coss to recognize promising genotypes for effective breeding. Two hundred and eighty-nine diverse accessions of Indian mustard belonging to four continents were analyzed for yield and yield-related traits (20 traits) over two seasons (2017-2018 and 2018-2019) using an alpha lattice design. The genetic variance was found to be significant (P ≤ 0.01) for the individual and under pooled analysis for all of the evaluated traits, demonstrating the presence of significant genetic variability in the diversity panel, which bids greater opportunities for utilizing these traits in future breeding programs. High heritability combined with high genetic advance as percent of mean and genotypic coefficient of variation was observed for flowering traits, plant height traits, seed size, and seed yield/plant; hence, a better genetic gain is expected upon the selection of these traits over subsequent generations. Both correlation and stepwise regression analysis indicated that the main shoot length, biological yield, total seed yield, plant height up to the first primary branch, seed size, total siliqua count, days to flowering initiation, plant height at maturity, siliquae on the main shoot, main shoot length, and siliqua length were the most significant contributory traits for seed yield/plant. Also, promising genotypes were identified among the diversity panel, which can be utilized as a donor to improve Indian mustard further. These results indicated a greater scope for improving seed yield per plant directly through a selection of genotypes having the parsimonious combination of these nine traits.

Nitrogen Signaling Genes and SOC1 Determine the Flowering Time in a Reciprocal Negative Feedback Loop in Chinese Cabbage (Brassica rapa L.) Based on CRISPR/Cas9-Mediated Mutagenesis of Multiple BrSOC1 Homologs

Precise flowering timing is critical for the plant life cycle. Here, we examined the molecular mechanisms and regulatory network associated with flowering in Chinese cabbage (Brassica rapa L.) by comparative transcriptome profiling of two Chinese cabbage inbred lines, “4004” (early bolting) and “50” (late bolting). RNA-Seq and quantitative reverse transcription PCR (qPCR) analyses showed that two positive nitric oxide (NO) signaling regulator genes, nitrite reductase (BrNIR) and nitrate reductase (BrNIA), were up-regulated in line “50” with or without vernalization. In agreement with the transcription analysis, the shoots in line “50” had substantially higher nitrogen levels than those in “4004”. Upon vernalization, the flowering repressor gene Circadian 1 (BrCIR1) was significantly up-regulated in line “50”, whereas the flowering enhancer genes named SUPPRESSOR OF OVEREXPRESSION OF CONSTANCE 1 homologs (BrSOC1s) were substantially up-regulated in line “4004”. CRISPR/Cas9-mediated mutagenesis in Chinese cabbage demonstrated that the BrSOC1-1/1-2/1-3 genes were involved in late flowering, and their expression was mutually exclusive with that of the nitrogen signaling genes. Thus, we identified two flowering mechanisms in Chinese cabbage: a reciprocal negative feedback loop between nitrogen signaling genes (BrNIA1 and BrNIR1) and BrSOC1s to control flowering time and positive feedback control of the expression of BrSOC1s.