Combined transcriptome and metabolome analysis reveals the effects of light quality on maize hybrids.

BACKGROUND: Heterosis, or hybrid vigor, refers to the phenotypic superiority of an F1 hybrid relative to its parents in terms of growth rate, biomass production, grain yield, and stress tolerance. Light is an energy source and main environmental cue with marked impacts on heterosis in plants. Research into the production applications and mechanism of heterosis has been conducted for over a century and a half, but little is known about the effect of light on plant heterosis. RESULTS: In this study, an integrated transcriptome and metabolome analysis was performed using maize (Zea mays L.) inbred parents, B73 and Mo17, and their hybrids, B73 × Mo17 (BM) and Mo17 × B73 (MB), grown in darkness or under far-red, red, or blue light. Most differentially expressed genes (73.72-92.50%) and differentially accumulated metabolites (84.74-94.32%) exhibited non-additive effects in BM and MB hybrids. Gene Ontology analysis revealed that differential genes and metabolites were involved in glutathione transfer, carbohydrate transport, terpenoid biosynthesis, and photosynthesis. The darkness, far-red, red, and blue light treatments were all associated with phenylpropanoid-flavonoid biosynthesis by Weighted Gene Co-expression Network Analysis and Kyoto Encyclopedia of Genes and Genomes enrichment analysis. Five genes and seven metabolites related to phenylpropanoid-flavonoid biosynthesis pathway were identified as potential contributors to the interactions between maize heterosis and light conditions. Consistent with the strong mid-parent heterosis observed for metabolites, significant increases in both fresh and dry weights were found in the MB and BM hybrids compared with their inbred parents. Unexpectedly, increasing light intensity resulted in higher biomass heterosis in MB, but lower biomass heterosis in BM. CONCLUSIONS: The transcriptomic and metabolomic results provide unique insights into the effects of light quality on gene expression patterns and genotype-environment interactions, and have implications for gene mining of heterotic loci to improve maize production.

Genomic diversity of aquaporins across genus Oryza provides a rich genetic resource for development of climate resilient rice cultivars.

BACKGROUND: Plant aquaporins are critical genetic players performing multiple biological functions, especially climate resilience and water-use efficiency. Their genomic diversity across genus Oryza is yet to be explored. RESULTS: This study identified 369 aquaporin-encoding genes from 11 cultivated and wild rice species and further categorized these into four major subfamilies, among which small basic intrinsic proteins are speculated to be ancestral to all land plant aquaporins. Evolutionarily conserved motifs in peptides of aquaporins participate in transmembrane transport of materials and their relatively complex gene structures provide an evolutionary playground for regulation of genome structure and transcription. Duplication and evolution analyses revealed higher genetic conservation among Oryza aquaporins and strong purifying selections are assisting in conserving the climate resilience associated functions. Promoter analysis highlighted enrichment of gene upstream regions with cis-acting regulatory elements involved in diverse biological processes, whereas miRNA target site prediction analysis unveiled substantial involvement of osa-miR2102-3p, osa-miR2927 and osa-miR5075 in post-transcriptional regulation of gene expression patterns. Moreover, expression patterns of japonica aquaporins were significantly perturbed in response to different treatment levels of six phytohormones and four abiotic stresses, suggesting their multifarious roles in plants survival under stressed environments. Furthermore, superior haplotypes of seven conserved orthologous aquaporins for higher thousand-grain weight are reported from a gold mine of 3,010 sequenced rice pangenomes. CONCLUSIONS: This study unveils the complete genomic atlas of aquaporins across genus Oryza and provides a comprehensive genetic resource for genomics-assisted development of climate-resilient rice cultivars.

Utilization of natural alleles for heat adaptability QTLs at the flowering stage in rice.

BACKGROUND: Heat stress threatens rice yield and quality at flowering stage. In this study, average relative seed setting rate under heat stress (RHSR) and genotypes of 284 varieties were used for a genome-wide association study. RESULTS: We identified eight and six QTLs distributed on chromosomes 1, 3, 4, 5, 7 and 12 in the full population and indica, respectively. qHTT4.2 was detected in both the full population and indica as an overlapping QTL. RHSR was positively correlated with the accumulation of heat-tolerant superior alleles (SA), and indica accession contained at least two heat-tolerant SA with average RHSR greater than 43%, meeting the needs of stable production and heat-tolerant QTLs were offer yield basic for chalkiness degree, amylose content, gel consistency and gelatinization temperature. Chalkiness degree, amylose content, and gelatinization temperature under heat stress increased with accumulation of heat-tolerant SA. Gel consistency under heat stress decreased with polymerization of heat-tolerant SA. The study revealed qHTT4.2 as a stable heat-tolerant QTL that can be used for breeding that was detected in the full population and indica. And the grain quality of qHTT4.2-haplotype1 (Hap1) with chalk5, wx, and alk was better than that of qHTT4.2-Hap1 with CHALK5, WX, and ALK. Twelve putative candidate genes were identified for qHTT4.2 that enhance RHSR based on gene expression data and these genes were validated in two groups. Candidate genes LOC_Os04g52830 and LOC_Os04g52870 were induced by high temperature. CONCLUSIONS: Our findings identify strong heat-tolerant cultivars and heat-tolerant QTLs with great potential value to improve rice tolerance to heat stress, and suggest a strategy for the breeding of yield-balance-quality heat-tolerant crop varieties.

Genome-wide investigation and expression analysis of APETALA-2 transcription factor subfamily reveals its evolution, expansion and regulatory role in abiotic stress responses in Indica Rice (Oryza sativa L. ssp. indica).

Rice is an important cereal crop that serves as staple food for more than half of the world population. Abiotic stresses resulting from changing climatic conditions are continuously threating its yield and production. Genes in APETALA-2 (AP2) family encode transcriptional regulators implicated during regulation of developmental processes and abiotic stress responses but their identification and characterization in indica rice was still missing. In this context, twenty-six genes distributed among eleven chromosomes in Indica rice encoding AP2 transcription-factor subfamily were identified and their diverse haplotypes were studied. Phylogenetic analysis of OsAP2 TF family-members grouped them into three clades indicating conservation of clades among cereals. Segmental duplications were observed to be principal route of evolution, supporting the higher positive selection-pressure, which were estimated to be originated about 10.57 to 56.72 million years ago (MYA). Conserved domain analysis and intron-exon distribution pattern of identified OsAP2s revealed their exclusive distribution among the specific clades of the phylogenetic tree. Moreover, the members of osa-miR172 family were also identified potentially targeting four OsAP2 genes. The real-time quantitative expression profiling of OsAP2s under heat stress conditions in contrasting indica rice genotypes revealed the differential expression pattern of OsAP2s (6 genes up-regulated and 4 genes down-regulated) in stress- and genotype-dependent manner. These findings unveiled the evolutionary pathways of AP2-TF in rice, and can help the functional characterization under developmental and stress responses.

Loci and natural alleles underlying robust roots and adaptive domestication of upland ecotype rice in aerobic conditions.

A robust (long and thick) root system is characteristic of upland japonica rice adapted to drought conditions. Using deep sequencing and large scale phenotyping data of 795 rice accessions and an integrated strategy combining results from high resolution mapping by GWAS and linkage mapping, comprehensive analyses of genomic, transcriptomic and haplotype data, we identified large numbers of QTLs affecting rice root length and thickness (RL and RT) and shortlisted relatively few candidate genes for many of the identified small-effect QTLs. Forty four and 97 QTL candidate genes for RL and RT were identified, and five of the RL QTL candidates were validated by T-DNA insertional mutation; all have diverse functions and are involved in root development. This work demonstrated a powerful strategy for highly efficient cloning of moderate- and small-effect QTLs that is difficult using the classical map-based cloning approach. Population analyses of the 795 accessions, 202 additional upland landraces, and 446 wild rice accessions based on random SNPs and SNPs within robust loci suggested that there could be much less diversity in robust-root candidate genes among upland japonica accessions than in other ecotypes. Further analysis of nucleotide diversity and allele frequency in the robust loci among different ecotypes and wild rice accessions showed that almost all alleles could be detected in wild rice, and pyramiding of robust-root alleles could be an important genetic characteristic of upland japonica. Given that geographical distribution of upland landraces, we suggest that during domestication of upland japonica, the strongest pyramiding of robust-root alleles makes it a unique ecotype adapted to aerobic conditions.

Distinct nucleotide patterns among three subgenomes of bread wheat and their potential origins during domestication after allopolyploidization.

BACKGROUND: The speciation and fast global domestication of bread wheat have made a great impact on three subgenomes of bread wheat. DNA base composition is an essential genome feature, which follows the individual-strand base equality rule and [AT]-increase pattern at the genome, chromosome, and polymorphic site levels among thousands of species. Systematic analyses on base compositions of bread wheat and its wild progenitors could facilitate further understanding of the evolutionary pattern of genome/subgenome-wide base composition of allopolyploid species and its potential causes. RESULTS: Genome/subgenome-wide base-composition patterns were investigated by using the data of polymorphic site in 93 accessions from worldwide populations of bread wheat, its diploid and tetraploid progenitors, and their corresponding reference genome sequences. Individual-strand base equality rule and [AT]-increase pattern remain in recently formed hexaploid species bread wheat at the genome, subgenome, chromosome, and polymorphic site levels. However, D subgenome showed the fastest [AT]-increase across polymorphic site from Aegilops tauschii to bread wheat than that on A and B subgenomes from wild emmer to bread wheat. The fastest [AT]-increase could be detected almost all chromosome windows on D subgenome, suggesting different mechanisms between D and other two subgenomes. Interestingly, the [AT]-increase is mainly contributed by intergenic regions at non-selective sweeps, especially the fastest [AT]-increase of D subgenome. Further transition frequency and sequence context analysis indicated that three subgenomes shared same mutation type, but D subgenome owns the highest mutation rate on high-frequency mutation type. The highest mutation rate on D subgenome was further confirmed by using a bread-wheat-private SNP set. The exploration of loci/genes related to the [AT] value of D subgenome suggests the fastest [AT]-increase of D subgenome could be involved in DNA repair systems distributed on three subgenomes of bread wheat. CONCLUSIONS: The highest mutation rate is detected on D subgenome of bread wheat during domestication after allopolyploidization, leading to the fastest [AT]-increase pattern of D subgenome. The phenomenon may come from the joint action of multiple repair systems inherited from its wild progenitors.

Insights on Calcium-Dependent Protein Kinases (CPKs) Signaling for Abiotic Stress Tolerance in Plants.

Abiotic stresses are the major limiting factors influencing the growth and productivity of plants species. To combat these stresses, plants can modify numerous physiological, biochemical, and molecular processes through cellular and subcellular signaling pathways. Calcium-dependent protein kinases (CDPKs or CPKs) are the unique and key calcium-binding proteins, which act as a sensor for the increase and decrease in the calcium (Ca) concentrations. These Ca flux signals are decrypted and interpreted into the phosphorylation events, which are crucial for signal transduction processes. Several functional and expression studies of different CPKs and their encoding genes validated their versatile role for abiotic stress tolerance in plants. CPKs are indispensable for modulating abiotic stress tolerance through activation and regulation of several genes, transcription factors, enzymes, and ion channels. CPKs have been involved in supporting plant adaptation under drought, salinity, and heat and cold stress environments. Diverse functions of plant CPKs have been reported against various abiotic stresses in numerous research studies. In this review, we have described the evaluated functions of plant CPKs against various abiotic stresses and their role in stress response signaling pathways.

Nucleotide diversity, natural variation, and evolution of Flexible culm-1 and Strong culm-2 lodging resistance genes in rice.

Lodging resistance is one of the vital traits in yield improvement and sustainability. Culm wall thickness, diameter, and strength are different traits that can govern the lodging resistance in rice. The genes SCM2 and FC1 have been isolated for culm thickness, strength, and flexibility, but their functional nucleotide variations were still unknown. We used a 13× deep sequence of 795 diverse genotypes to present the functional variation and SNP diversity in SCM2 and FC1. The major functional variant for the SCM2 gene was at position 27480181 and for the FC1 gene at position 31072992. Haplotype analysis of both genes provided their various allelic differences among haplotypes. SCM2 alleles further presented the evolution of Oryza sativa L. subsp. indica and subsp. japonica genomes from common parent in different geographical zones, while the haplotypes of FC1 suggested their evolution from different strains of the common parent Oryza rufipogon. SCM2 showed purifying selection and functional associations with rare alleles, while FC1 displayed balanced selection favored by multiple heterozygous alleles. Genotypes with an allelic combination of SCM2-3 and FC1-2 in japonica background exhibited striking resistance against lodging, which can be used in further breeding programs.

Identification of novel inhibitors against VP40 protein of Marburg virus by integrating molecular modeling and dynamics approaches.

Marburg virus (MV) is a highly etiological agent of haemorrhagic fever in humans and has spread across the world. Its outbreaks caused a 23-90% human death rate. However, there are currently no authorized preventive or curative measures yet. VP40 is the MV matrix protein, which builds protein shell underneath the viral envelope and confers hallmark filamentous. VP40 alone is able to induce assembly and budding of filamentous virus-like particles (VLPs), which resemble authentic virions. As a result, this research is credited with clarifying the function of VP40 and leading to the discovery of new therapeutic targets effective in combating MV disease (MVD). Virtual screening, molecular docking and molecular dynamics (MD) simulation were used to find the putative active chemicals based on a 3D pharmacophore model of the protein's active site cavity. Initially, andrographidine-C, a potent inhibitor was selected for the development of the pharmacophore model. Later, a library of 30,000 compounds along with the andrographidine-C was docked against VP40 protein. Three best hits including avanafil, diuvaretin and macrourone were subjected to further MD simulation analysis, as these compounds had better binding affinities as compared to andrographidine-C. Furthermore, throughout the 100 ns simulations, the back bone of VP40 protein in presence of avanafil, diuvaretin and macrourone remained stable which was further validated by MM-PBSA analysis. Additionally, all of these compounds depict maximum drug-like properties. The predicted drugs based on the ligand, avanafil, diuvaretin and macrourone could be exploited and developed as an alternative or complementary therapy for the treatment of MVD.Communicated by Ramaswamy H. Sarma.

Microsatellite markers-aided dissection of iron, zinc and cadmium accumulation potential in Triticum aestivum.

Background: Wheat is a staple cereal food around the globe. It provides a significant source of proteins, carbohydrates, and other micronutrients to humans. When grown on cadmium (Cd) contaminated soils, the uptake of trace elements e.g., iron (Fe) and zinc (Zn) has also been affected drastically that in turn affected the wheat grain. Methods: In this study, wheat accessions were used to investigate the impact of soil application of Zn (5 mg/kg, 20 mg/kg) and Cd (0 mg/kg, 10 mg/kg) on accumulation of these elements in wheat grains. A total of 45 Fe, Zn, and Cd transporter-related genes were used to design 101 gene-specific SSR (simple sequence repeat) markers. Results: In response to Cd stress, application of 20 mg/Kg Zn improved Fe (64.6 ug/g) and Zn (48.3 ug/g) accumulation in wheat grains as well as agronomic traits. Marker trait association revealed that SSR markers based on NAM-B1 gene (PR01 and PR02) were associated with Zn accumulation. Similarly, SSR markers based on TaVTL5-2B_5 (PR19 PR20), TaVTL5-2B_2 (PR25, PR26), TaVTL5-2D_3 (PR30), TaVTL2-2A (PR31), TaVTL1-6A (PR32), and TaVTL2-2D_1 (PR37) were significantly associated with Fe accumulation, while HMA3-5B1 (PR62) and TaNRAMP3-7D (PR89) were linked to Cd accumulation in grains. The highly associated markers may be used in marker-assisted selection of suitable wheat genotypes for breeding bio-fortified varieties with low Cd accumulation.

Contribution of Duplicated Nucleotide-Binding Leucine-Rich Repeat (NLR) Genes to Wheat Disease Resistance.

Wheat has a large and diverse repertoire of NLRs involved in disease resistance, with over 1500 NLRs detected in some studies. These NLR genes occur as singletons or clusters containing copies of NLRs from different phylogenetic clades. The number of NLRs and cluster size can differ drastically among ecotypes and cultivars. Primarily, duplication has led to the evolution and diversification of NLR genes. Among the various mechanisms, whole genome duplication (WGD) is the most intense and leading cause, contributing to the complex evolutionary history and abundant gene set of hexaploid wheat. Tandem duplication or recombination is another major mechanism of NLR gene expansion in wheat. The diversity and divergence of duplicate NLR genes are responsible for the broad-spectrum resistance of most plant species with limited R genes. Understanding the mechanisms underlying the rapid evolution and diversification of wheat NLR genes will help improve disease resistance in crops. The present review focuses on the diversity and divergence of duplicate NLR genes and their contribution to wheat disease resistance. Moreover, we provide an overview of disease resistance-associated gene duplication and the underlying strategies in wheat.

Genome wide association study on development and evolution of glutinous rice.

BACKGROUND: Glutinous rice as a special endosperm type is consumed as a staple food in East Asian countries by consumers' preference. Genetic studies on glutinous rice could be conducive to improve rice quality and understand its development and evolution. Therefor, we sought to explore more genes related to glutinous by genome wide association study and research the formation history for glutinous. RESULTS: Here, genome-wide association study was performed to explore the associated loci/genes underlying glutinous rice by using 2108 rice accessions. Combining the expression patterns analysis, 127, 81, and 48 candidate genes were identified to be associated with endosperm type in whole rice panel, indica, and japonica sub-populations. There were 32 genes, including three starch synthesis-related genes Wx, SSG6, and OsSSIIa, detected simultaneously in the whole rice panel and subpopulations, playing important role in determining glutinous rice. The combined haplotype analyses revealed that the waxy haplotypes combination of three genes mainly distributed in Southeast Asia (SEA), SEA islands (SER) and East Asia islands (EAR). Through population structure and genetic differentiation, we suggest that waxy haplotypes of the three genes firstly evolved or were directly inherited from wild rice in japonica, and then introgressed into indica in SER, SEA and EAR. CONCLUSIONS: The cloning and natural variation analysis of waxy-related genes are of great significance for the genetic improvement of quality breeding and comprehend the history in glutinous rice. This work provides valuable information for further gene discovery and understanding the evolution and formation for glutinous rice in SEA, SER and EAR.

Computational analysis and expression profiling of potassium transport-related gene families in mango (Mangifera indica) indicate their role in stress response and fruit development.

Mango (Mangifera indica) fruit is known for its taste, health benefits, and drought tolerance. Potassium (K+) is one of the most abundant ions in a plant cell. It is important for various biological functions related to plant growth, development, and flowering/fruiting. It significantly contributes to fruit yield, quality, and drought tolerance in plants. However, molecular mechanisms comprising K+ transport in mango are least known. In the present study, 37 members of K+ transport-related genes (PTGs) were identified in mango, which include 22 K+ transporters (16 HAKs, 1 HKT, and 6 KEAs) and 15 K+ channels (6 TPKs and 8 Shakers). All PTGs were predicted to be expressed at the plasma membrane and possess characteristic motifs and domains. Phylogenetic analysis identified a strong kinship of PTGs among Oryza sativa, Arabidopsis thaliana, Cicer arietinum, Malus domestica, and M. indica. The promoter analysis identified 60 types of cis-elements related to various biological processes. RNA-seq-based expression profiling identified that MiTPK1.2, MiHAK1, MiHAK2.1, HAK6.1, and MiAKT1.1 were most upregulated in roots and that MiKEA2, MiAKT2, and MiAKT1 were upregulated in leaves. Moreover, MiAKT6, MiHAK1.1, MiKAT2, MiKAT2.1, MiHKT1, MiTPK1.1, MiHAK7, and MiHAK12 were highly expressed during the five growth stages of mango fruit. The current study is the first comprehensive report on K+ transport system in tropical fruits. Therefore, it will provide the foundation knowledge for the functional characterization of K+ genes in mango and related plants.

Interaction estimation of pathogenicity determinant protein ßC1 encoded by Cotton leaf curl Multan Betasatellite with Nicotiana benthamiana Nuclear Transport Factor 2.

Background: Begomovirus is one of the most devastating pathogens that can cause more than 90% yield loss in various crop plants. The pathogenicity determinant ßC1, located on the betasatellite associated with monopartite begomoviruses, alters the host signaling mechanism to enhance the viral disease phenotype by undermining the host immunity. The understanding of its interacting proteins in host plants to develop disease symptoms such as curly leaves, enations, vein swelling, and chlorosis is crucial to enhance the disease resistance in crop plants. The current study was designed to reveal the contribution of ßC1 in disease pathogenicity and to unveil potential interacting partners of ßC1 protein in the model plant Nicotiana benthamiana. Methods: The ßC1 gene was cloned in pGKBT7 and used as bait against the cDNA library of N. benthamiana and its pathogenesis was tested against the healthy plant and the plants infiltrated with empty vectors. The yeast two-hybrid-based screening was performed to find the interacting factors. Successful interacting proteins were screened and evaluated in various steps and confirmed by sequence analysis. The three-dimensional structure of the Nuclear Transport Factor 2 (NTF2) protein was predicted, and in-silico protein-protein interaction was evaluated. Furthermore, protein sequence alignment and molecular phylogenetic analysis were carried out to identify its homologues in other related families. In-silico analyses were performed to validate the binding affinity of ßC1 protein with NTF2. The 3D model was predicted by using I-TASSER and then analyzed by SWISS MODEL-Workspace, RAMPAGE, and Verify 3D. The interacting amino acid residues of ßC1 protein with NTF2 were identified by using PyMOL and Chimera. Results: The agroinfiltrated leaf samples developed severe phenotypic symptoms of virus infection. The yeast-two-hybrid study identified the NTF2 as a strong interacting partner of the ßC1. The NTF2 in Solanaceae and Nicotiana was found to be evolved from the Brassica and Gossypium species. The in-silico interaction studies showed a strong binding affinity with releasing energy value of -730.6 KJ/mol, and the involvement of 10 amino acids from the middle portion towards the C-terminus and five amino acid residues from the middle portion of ßC1 to interact with six amino acids of NTF2. The study not only provided an insight into the molecular mechanism of pathogenicity but also put the foundation stone to develop the resistance genotypes for commercial purposes and food security.

Genome-wide association mapping for high temperature tolerance in wheat through 90k SNP array using physiological and yield traits.

Dissecting the genetic basis of physiological and yield traits against tolerance to heat stress is an essential in wheat breeding programs to boost up the wheat yield for sustainable food security. Herein, a genome-wide association study (GWAS) was performed to reveal the genetic basis of heat tolerance using high-density Illumina 90K Infinium SNPs array through physiological and yield indices. These indices were phenotyped on a diverse panel of foreign and domestic genotypes of Pakistan, grown in normal and heat-stressed environments. Based on STRUCTURE analysis, the studied germplasm clustered into four sub-population. Highly significant variations with a range of moderate (58.3%) to high (77.8%) heritability was observed under both conditions. Strong positive correlation existed among physiological and yield related attributes. A total of 320 significant (-log10 P ≥ 3) marker-trait associations (MTAs) were identified for the observed characters. Out of them 169 and 151 MTAs were recorded in normal and heat stress environments, respectively. Among the MTA loci, three (RAC875_c103017_302, Tdurum_contig42087_1199, and Tdurum_contig46877_488 on chromosomes 4B, 6B, and 7B respectively), two (BobWhite_c836_422 and BS00010616_51) and three (Kukri_rep_c87210_361, D_GA8KES401BNLTU_253 and Tdurum_contig1015_131) on chromosomes 5A, 1B, and 3D at the positions 243.59cM, 77.82cM and 292.51cM) showed pleiotropic effects in studied traits under normal, heat-stressed and both conditions respectively. The present study not only authenticated the numerous previously reported MTAs for examined attributes but also revealed novel MTAs under heat-stressed conditions. Identified SNPs will be beneficial in determining the novel genes in wheat to develop the heat tolerant and best yielded genotypes to fulfill the wheat requirement for the growing population.

Genome-wide association mapping of bread wheat genotypes using yield and grain morphology-related traits under different environments.

Depleting water resources and increasing global temperature due to climate change are major challenges to agriculture and food security worldwide. Deciphering the underlying mechanisms of traits contributing to grain development and yield is essential for the development of climate-resilient cultivars. Therefore, this study assessed 105 bread wheat genotypes grown under control, drought, and heat-stress conditions for two crop seasons and performed a genome-wide association study (GWAS) using a 90k SNP array. The genotypes showed significant trait differences under all environmental conditions. Highly significant variation was observed, with moderate (50.09%) to high (76.19%) heritability in the studied germplasms. The studied traits were all also significantly positively correlated. A total of 541 significant associations (p ≤ 10-3) between marker and trait (MTAs) were observed after crossing the FDR <0.05 threshold for all traits. Among these, 195, 179, and 167 significant MTAs were detected under control, drought, and heat-stress conditions, respectively. Under the control and drought conditions, pleiotropic loci BS00010616_51 and BS00010868_51 were observed on chromosomes 7B and 1B situated at 186.24 cM and 35.47 cM, respectively. Pleiotropic loci BS00010868_51, Kukri_c11154_1723, and Ex_c10068_1509 were identified on chromosomes 1B, 5B, and 2A, respectively, under control and heat stress conditions. A stable and consistent locus (Excalibur_c20796_395) on chromosome 7A, located at 372.34 cM, was also linked to grain morphology and yield-related attributes in control, drought, and heat-stress conditions. The results of the current study confirmed several previously reported MTAs for the traits under consideration and identified new MTAs under harsh climatic conditions. These SNPs will aid in the discovery of novel genes in wheat. SNPs showing significant associations may be used in marker-assisted selection and allow the development of drought- and heat-tolerant genotypes with high yields to address global food security concerns.

Genome-wide association mapping of bread wheat genotypes for sustainable food security and yield potential under limited water conditions.

Determining the genetic basis of yield and water deficient tolerance in wheat is vital for wheat breeding programs. Herein, a genome-wide association study (GWAS) was performed for water deficient and yield-related attributes on wheat genotypes with high-density Illumina 90K Infinium SNP array. Major yield and drought-related attributes were phenotyped on a panel of Pakistani and foreign accessions grown in non-stressed and water deficient stressed environments during two crop cycles. Among all accessions, highly significant variations were shown in studied environments for examined characters. Water deficient conditions, reduced the wheat yield and had strong and positive correlation among relative water content and grain yield per plant. Population structure analyses based on 90,000 SNP data, classify the accessions into 4 sub-populations. Marker-trait association analyses (MTA) revealed that 134 significant SNPs were linked with yield and drought tolerance attributes. Pleotropic loci RAC875_s117925_244 and RAC875_c16333_340 located on chromosome 5A and 2A respectively, were significantly linked with relative water contents (RWC), cell membrane thermo-stability (CMT), grain per spike (GPS), spikelet per spike (SPS) and grain yield per plant (GYP). The markers Ra_c58279_684, BobWhite_c23828_341 and IAAV3414 located on chromosomes 2A, 6B and 7B respectively, showed pleotropic effects for RWC, GPS and GYP under both environments. The current experiment not only validated several MTAs reported in other studies but also discovered novel MTAs which significant under drought-stressed conditions. A total of 171 candidate genes were recognized that could be cloned and functionally characterized for the respective associated traits. For RWC and CMT, total 11 and 3 associated SNPs were mapped on coding DNA sequence (CDS) of the identified candidate genes. Isolation and characterization of the candidate genes herein mapped SNPs will be useful in discovering novel genes underpinning drought tolerance in bread wheat to fulfill the wheat demand and sustainable food security under limited water conditions.

Unveiling the genetic architecture for lodging resistance in rice (Oryza sativa. L) by genome-wide association analyses.

Lodging is one of the major abiotic stresses, affecting the total crop yield and quality. The improved lodging resistance and its component traits potentially reduce the yield losses. The section modulus (SM), bending moment at breaking (M), pushing resistance (PR), and coefficient of lodging resistance (cLr) are the key elements to estimate the lodging resistance. Understanding the genetic architecture of lodging resistance-related traits will help to improve the culm strength and overall yield potential. In this study, a natural population of 795 globally diverse genotypes was further divided into two (indica and japonica) subpopulations and was used to evaluate the lodging resistance and culm strength-related traits. Significant diversity was observed among the studied traits. We carried out the genome-wide association evaluation of four lodging resistance traits with 3.3 million deep resolution single-nucleotide polymorphic (SNP) markers. The general linear model (GLM) and compressed mixed linear model (MLM) were used for the whole population and two subpopulation genome-wide association studies (GWAS), and a 1000-time permutation test was performed to remove the false positives. A total of 375 nonredundant QTLs were observed for four culm strength traits on 12 chromosomes of the rice genome. Then, 33 pleiotropic loci governing more than one trait were mined. A total of 4031 annotated genes were detected within the candidate genomic region of 33 pleiotropic loci. The functional annotations and metabolic pathway enrichment analysis showed cellular localization and transmembrane transport as the top gene ontological terms. The in silico and in vitro expression analyses were conducted to validate the three candidate genes in a pleiotropic QTL on chromosome 7. It validated OsFBA2 as a candidate gene to contribute to lodging resistance in rice. The haplotype analysis for the candidate gene revealed a significant functional variation in the promoter region. Validation and introgression of alleles that are beneficial to induce culm strength may be used in rice breeding for lodging resistance.

Exploration of the genomic atlas of Dof transcription factor family across genus Oryza provides novel insights on rice breeding in changing climate.

DNA-binding with one finger (Dof) transcription factors have been demonstrated to regulate various stresses and developmental processes in plants. Their identification and comparative evolutionary analyses in cultivated and wild species of genus oryza were yet to be explored. In this context, we report a comprehensive genomics atlas of DNA-binding with one finger (Dof) family genes in 13 diverse rice genomes (five cultivated and eight rice wild-relatives) through a genome-wide scanning approach. A galore of 238 Dof genes, identified across the genus Oryza, are categorized into seven distinct subgroups by comparative phylogenetic analysis with the model plant Arabidopsis. Conserved motifs and gene structure analyses unveiled the prevalence of species- and subgroups-specific structural and functional diversity that is expediating with the evolutionary period. Our results indicate that Dof genes might have undergone strong purifying selections and segmental duplications to expand their gene family members in corresponding Oryza genomes. We speculate that miR2927 potentially targets the Dof domain to regulate gene expression under different climatic conditions, which are supported by in-silico and wet-lab experiments-based expression profiles. In a nutshell, we report several superior haplotypes significantly associated with early flowering in a treasure trove of 3,010 sequenced rice accessions and have validated these haplotypes with two years of field evaluation-based flowering data of a representative subpanel. Finally, we have provided some insights on the resolution of Oryza species phylogeny discordance and divergence highlighting the mosaic evolutionary history of the genus Oryza. Overall, this study reports a complete genomic landscape of the Dof family in cultivated and wild Oryza species that could greatly facilitate in fast-track development of early maturing and climate-resilient rice cultivars through modern haplotype-led breeding.

Systematic Analysis of NB-ARC Gene Family in Rice and Functional Characterization of GNP12.

The NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) gene family plays a critical role in plant development. However, our understanding of the mechanisms of how NB-ARC genes regulate plant development in the plant panicle is still limited. Here, we subjected 258 NB-ARC genes in rice to genome-wide analysis to characterize their structure, function, and expression patterns. The NB-ARC genes were classified into three major groups, and group II included nine subgroups. Evolutionary analysis of NB-ARC genes in a dicotyledon plant (Arabidopsis thaliana) and two monocotyledonous plants (Oryza sativa L. and Triticum aestivum) indicated that homologous genome segments were conserved in monocotyledons and subjected to weak positive selective pressure during evolution. Dispersed and proximal replication events were detected. Expression analysis showed expression of most NB-ARC genes in roots, panicles, and leaves, and regulation at the panicle development stage in rice Ce253. The GNP12 gene encodes RGH1A protein, which regulates rice yield according to panicle length, grain number of panicle, and grain length, with eight major haplotypes. Most members of NB-ARC protein family are predicted to contain P-loop conserved domains and localize on the membrane. The results of this study will provide insight into the characteristics and evolution of NB-ARC family and suggest that GNP12 positively regulates panicle development.