Genome-Wide Association Study Reveals Novel QTLs and Candidate Genes for Grain Number in Rice

Grain number per panicle (GNPP), determined mainly by panicle branching, is vital for rice yield. The dissection of the genetic basis underlying GNPP could help to improve rice yield. However, genetic resources, including quantitative trait loci (QTL) or genes for breeders to enhance rice GNPP, are still limited. Here, we conducted the genome-wide association study (GWAS) on the GNPP, primary branch number (PBN), and secondary branch number (SBN) of 468 rice accessions. We detected a total of 18 QTLs, including six for GNPP, six for PBN, and six for SBN, in the whole panel and the indica and japonica subpanels of 468 accessions. More importantly, qPSG1 was a common QTL for GNPP, PBN, and SBN and was demonstrated by chromosome segment substitution lines (CSSLs). Considering gene annotation, expression, and haplotype analysis, seven novel and strong GNPP-related candidate genes were mined from qPSG1. Our results provide clues to elucidate the molecular regulatory network of GNPP. The identified QTLs and candidate genes will contribute to the improvement of GNPP and rice yield via molecular marker-assisted selection (MAS) breeding and genetic engineering techniques.     

RNA-Seq and Genome-Wide Association Studies Reveal Potential Genes for Rice Seed Shattering

The loss of the shattering ability is one of the key events in rice domestication. The strength of the seed shattering ability is closely related to the harvest yield and the adaptability of modern mechanical harvesting methods. In this study, using a population of 587 natural rice cultivars, quantitative trait loci associated with seed shattering were detected by genome-wide association studies (GWASs). We consider the quantitative trait loci (QTLs) qBTS1 and qBTS3 to be the key loci for seed shattering in rice. Additionally, the abscission zone (AZ) and nonabscission zone (NAZ) of materials with a loss of shattering (DZ129) and easy shattering (W517) were subjected to RNA-Seq, and high-quality differential expression profiles were obtained. The AZ-specific differentially expressed genes (DEGs) of W517 were significantly enriched in plant hormone signal transduction, while the AZ-specific DEGs of DZ129 were enriched in phenylpropanoid biosynthesis. We identified candidate genes for the lignin-associated laccase precursor protein (LOC_Os01g63180) and the glycoside hydrolase family (LOC_Os03g14210) in the QTLs qBTS1 (chromosome 1) and qBTS3 (chromosome 3), respectively. In summary, our findings lay the foundation for the further cloning of qBTS1 and qBTS3, which would provide new insights into seed shattering in rice.     

Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein

Rice blast is a worldwide fungal disease that seriously affects the yield and quality of rice. Identification of resistance genes against rice blast disease is one of the effective ways to control this disease. However, panicle blast resistance genes, which are useful in the fields, have rarely been studied due to the difficulty in phenotypic identification and the environmental influences. Here, panicle blast resistance-3 (Pb3) was identified by a genome-wide association study (GWAS) based on the panicle blast resistance phenotypes of 230 Rice Diversity Panel I (RDP-I) accessions with 700,000 single-nucleotide polymorphism (SNP) markers. A total of 16 panicle blast resistance loci (PBRLs) within three years including one repeated locus PBRL3 located in chromosome 11 were identified. In addition, 7 genes in PBRL3 were identified as candidate genes by haplotype analysis, which showed significant differences between resistant and susceptible varieties. Among them, one nucleotide-binding domain and Leucine-rich Repeat (NLR) gene Pb3 was highly conserved in multiple resistant rice cultivars, and its expression was significantly induced after rice blast inoculation. Evolutionary analysis showed that Pb3 was a typical disease resistance gene containing coiled-coil, NB-ARC, and LRR domains. T-DNA insertion mutants and CRISPR lines of Pb3 showed significantly reduced panicle blast resistance. These results indicate that Pb3 is a panicle blast resistance gene and GWAS is a rapid method for identifying panicle blast resistance in rice.     

The Integration of Genome-Wide Association Study and Homology Analysis to Explore the Genomic Regions and Candidate Genes for Panicle-Related Traits in Foxtail Millet

Panicle traits are important factors affecting yield, and their improvement has long been a critical goal in foxtail millet breeding. In order to understand the genetic basis of panicle formation, a large-scale genome-wide association study (GWAS) was performed in this study for six panicle-related traits based on 706,646 high-polymorphism SNP loci in 407 accessions. As a result, 87 quantitative trait loci (QTL) regions with a physical distance of less than 100 kb were detected to be associated with these traits in three environments. Among them, 27 core regions were stably detected in at least two environments. Based on rice–foxtail millet homologous comparison, expression, and haplotype analysis, 27 high-confidence candidate genes in the QTL regions, such as Si3g11200 (OsDER1), Si1g27910 (OsMADS6), Si7g27560 (GS5), etc., affected panicle-related traits by involving multiple plant growth regulator pathways, a photoperiod response, as well as panicle and grain development. Most of these genes showed multiple effects on different panicle-related traits, such as Si3g11200 affecting all six traits. In summary, this study clarified a strategy based on the integration of GWAS, a homologous comparison, and haplotype analysis to discover the genomic regions and candidate genes for important traits in foxtail millet. The detected QTL regions and candidate genes could be further used for gene clone and marker-assisted selection in foxtail millet breeding.     

Identification of a Negative Regulator for Salt Tolerance at Seedling Stage via a Genome-Wide Association Study of Thai Rice Populations

Salt stress is a major limiting factor in crop production and yield in many regions of the world. The objective of this study was to identify the genes responsible for salt tolerance in Thai rice populations. We performed a genome-wide association study with growth traits, relative water content, and cell membrane stability at the seedling stage, and predicted 25 putative genes. Eleven of them were located within previously reported salt-tolerant QTLs (ST-QTLs). OsCRN, located outside the ST-QTLs, was selected for gene characterization using the Arabidopsis mutant line with T-DNA insertion in the orthologous gene. Mutations in the AtCRN gene led to the enhancement of salt tolerance by increasing the ability to maintain photosynthetic pigment content and relative water content, while the complemented lines with ectopic expression of OsCRN showed more susceptibility to salt stress detected by photosynthesis performance. Moreover, the salt-tolerant rice varieties showed lower expression of this gene than the susceptible rice varieties under salt stress conditions. The study concludes that by acting as a negative regulator, OsCRN plays an important role in salt tolerance in rice.     

Genome-Wide Association Analysis of Senescence-Related Traits in Maize

Senescence is a programmed process that involves the destruction of the photosynthesis apparatus and the relocation of nutrients to the grain. Identifying senescence-associated genes is essential to adapting varieties for the duration of the cultivation cycle. A genome-wide association study (GWAS) was performed using 400 inbred maize lines with 156,164 SNPs to study the genetic architecture of senescence-related traits and their relationship with agronomic traits. We estimated the timing of senescence to be 45 days after anthesis in the whole plant and specifically in the husks. A list of genes identified in a previous RNAseq experiment as involved in senescence (core senescence genes) was used to propose candidate genes in the vicinity of the significant SNPs. Forty-six QTLs of moderate to high effect were found for senescence traits, including specific QTLs for husk senescence. The allele that delayed senescence primarily increased grain yield and moisture. Seven and one significant SNPs were found in the coding and promoter regions of eight core senescence genes, respectively. These genes could be potential candidates for generating a new variation by genome editing for functional analysis and breeding purposes, particularly Zm00001d014796, which could be responsible for a QTL of senescence found in multiple studies.     

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H² = 0.5–0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.     

Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene,Pb2, Encoding NLR Protein

Rice blast is one of the main diseases in rice and can occur in different rice growth stages. Due to the complicated procedure of panicle blast identification and instability of panicle blast infection influenced by the environment, most cloned rice resistance genes are associated with leaf blast. In this study, a rice panicle blast resistance gene, Pb2, was identified by genome-wide association mapping based on the panicle blast resistance phenotypes of 230 Rice Diversity Panel 1 (RDP1) accessions with 700,000 single-nucleotide polymorphism (SNP) markers. A genome-wide association study identified 18 panicle blast resistance loci (PBRL) within two years, including 9 reported loci and 2 repeated loci (PBRL2 and PBRL13, PBRL10 and PBRL18). Among them, the repeated locus (PBRL10 and PBRL18) was located in chromosome 11. By haplotype and expression analysis, one of the Nucleotide-binding domain and Leucine-rich Repeat (NLR) Pb2 genes was highly conserved in multiple resistant rice cultivars, and its expression was significantly upregulated after rice blast infection. Pb2 encodes a typical NBS-LRR protein with NB-ARC domain and LRR domain. Compared with wild type plants, the transgenic rice of Pb2 showed enhanced resistance to panicle and leaf blast with reduced lesion number. Subcellular localization of Pb2 showed that it is located on plasma membrane, and GUS tissue-staining observation found that Pb2 is highly expressed in grains, leaf tips and stem nodes. The Pb2 transgenic plants showed no difference in agronomic traits with wild type plants. It indicated that Pb2 could be useful for breeding of rice blast resistance.     

全基因组关联分析在植物中的应用

全基因组关联分析(GWAS)已被广泛应用到植物遗传学和育种相关的重要性状研究中。关联分析是一种基于连锁不平衡来识别分子标记之间或候选基因与性状之间关系的方法。简单介绍了GWAS的发展背景、研究原理及研究策略,对GWAS在重要的植物基因位点方面的应用研究进展进行了综述。     

Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.)

The root tissues play important roles in water and nutrient acquisition, environmental adaptation, and plant development. In this study, a diversity panel of 388 wheat accessions was collected to investigate nine root system architecture (RSA) traits at the three-leaf stage under two growing environments: outdoor pot culture (OPC) and indoor pot culture (IPC). Phenotypic analysis revealed that root development was faster under OPC than that under IPC and a significant correlation was observed between the nine RSA traits. The 660K single-nucleotide polymorphism (SNP) chip was used for a genome-wide association study (GWAS). Significant SNPs with a threshold of −log₁₀ (p-value) ≥ 4 were considered. Thus, 36 quantitative trait loci (QTLs), including 13 QTL clusters that were associated with more than one trait, were detected, and 31 QTLs were first identified. The QTL clusters on chromosomes 3D and 5B were associated with four and five RSA traits, respectively. Two candidate genes, TraesCS2A01G516200 and TraesCS7B01G036900, were found to be associated with more than one RSA trait using haplotype analysis, and preferentially expressed in the root tissues. These favourable alleles for RSA traits identified in this study may be useful to optimise the root system in wheat.     

Grain Size Associated Genes and the Molecular Regulatory Mechanism in Rice

Grain size is a quantitative trait that is controlled by multiple genes. It is not only a yield trait, but also an important appearance quality of rice. In addition, grain size is easy to be selected in evolution, which is also a significant trait for studying rice evolution. In recent years, many quantitative trait loci (QTL)/genes for rice grain size were isolated by map-based cloning or genome-wide association studies, which revealed the genetic and molecular mechanism of grain size regulation in part. Here, we summarized the QTL/genes cloned for grain size and the regulation mechanism with a view to provide the theoretical basis for improving rice yield and breeding superior varieties.     

Genome-Wide Identification and Analysis of the Metallothionein Genes in Oryza Genus

Metallothionein (MT) proteins are low molecular mass, cysteine-rich, and metal-binding proteins that play an important role in maintaining metal homeostasis and stress response. However, the evolutionary relationships and functional differentiation of MT in the Oryza genus remain unclear. Here we identified 53 MT genes from six Oryza genera, including O. sativa ssp. japonica, O. rufipogon, O. sativa ssp. indica, O. nivara, O. glumaepatula, and O. barthii. The MT genes were clustered into four groups based on phylogenetic analysis. MT genes are unevenly distributed on chromosomes; almost half of the MT genes were clustered on chromosome 12, which may result from a fragment duplication containing the MT genes on chromosome 12. Five pairs of segmental duplication events and ten pairs of tandem duplication events were found in the rice MT family. The Ka/Ks values of the fifteen duplicated MT genes indicated that the duplicated MT genes were under a strong negative selection during evolution. Next, combining the promoter activity assay with gene expression analysis revealed different expression patterns of MT genes. In addition, the expression of OsMT genes was induced under different stresses, including NaCl, CdCl₂, ABA, and MeJ treatments. Additionally, we found that OsMT genes were mainly located in chloroplasts. These results imply that OsMT genes play different roles in response to these stresses. All results provide important insights into the evolution of the MT gene family in the Oryza genus, and will be helpful to further study the function of MT genes.     

The Analysis of a Genome-Wide Association Study (GWAS) of Overweight and Obesity in Psoriasis

There is evidence that the concomitance of psoriasis and obesity may originate from the interplay between multiple genetic pathways and involve gene–gene interactions. The aim of this study was to compare the genetic background related to obesity among psoriatic patients versus healthy controls by means of a Genome-Wide Association Study (GWAS). A total of 972 psoriatic patients and a total of 5878 healthy donors were enrolled in this study. DNA samples were genotyped for over 500,000 single nucleotide polymorphisms (SNPs) using Infinium CoreExome BeadChips (Illumina, San Diego, CA, USA). Statistical analysis identified eleven signals (p < 1 × 10⁻⁵) associated with BMI across the study groups and revealed a varying effect size in each sub-cohort. Seven of the alternative alleles (rs1558902 in the FTO gene, rs696574 in the CALCRL gene, as well as rs10968110, rs4551082, rs4609724, rs9320269, and rs2338833,) are associated with increased BMI among all psoriatic patients and four (rs1556519 in the ITLN2 gene, rs12972098 in the AC003006.7 gene, rs12676670 in the PAG1 gene, and rs1321529) are associated with lower BMI. The results of our study may lead to further insights into the understanding of the pathogenesis of obesity among psoriatic patients.     

Genome-Wide Association Mapping Identifies New Candidate Genes for Cold Stress and Chilling Acclimation at Seedling Stage in Rice (Oryza sativa L.)

Rice (Oryza sativa L.) is a chilling-sensitive staple food crop, and thus, low temperature significantly affects rice growth and yield. Many studies have focused on the cold shock of rice although chilling acclimation is more likely to happen in the field. In this paper, a genome-wide association study (GWAS) was used to identify the genes that participated in cold stress and chilling accumulation. A total of 235 significantly associated single-nucleotide polymorphisms (SNPs) were identified. Among them, we detected 120 and 88 SNPs for the relative shoot fresh weight under cold stress and chilling acclimation, respectively. Furthermore, 11 and 12 quantitative trait loci (QTLs) were identified for cold stress and chilling acclimation, respectively, by integrating the co-localized SNPs. Interestingly, we identified 10 and 15 candidate genes in 11 and 12 QTLs involved in cold stress and chilling acclimation, respectively, and two new candidate genes (LOC_Os01g62410, LOC_Os12g24490) were obviously up-regulated under chilling acclimation. Furthermore, OsMYB3R-2 (LOC_Os01g62410) that encodes a R1R2R3 MYB gene was associated with cold tolerance, while a new C3HC4-type zinc finger protein-encoding gene LOC_Os12g24490 was found to function as a putative E3 ubiquitin-protein ligase in rice. Moreover, haplotype, distribution, and Wright’s fixation index (FST) of both genes showed that haplotype 3 of LOC_Os12g24490 is more stable in chilling acclimation, and the SNP (A > T) showed a difference in latitudinal distribution. FST analysis of SNPs in OsMYB3R-2 (LOC_Os01g62410) and LOC_Os12g24490 indicated that several SNPs were under selection in rice indica and japonica subspecies. This study provided new candidate genes in genetic improvement of chilling acclimation response in rice.     

Genome-Wide Identification of Barley Long Noncoding RNAs and Analysis of Their Regulatory Interactions during Shoot and Grain Development

Long noncoding RNAs (lncRNAs) are a class of RNA molecules with gene regulatory functions in plant development and the stress response. Although the number of lncRNAs identified in plants is rapidly increasing, very little is known about their role in barley development. In this study, we performed global identification of barley lncRNAs based on 53 RNAseq libraries derived from nine different barley tissues and organs. In total, 17,250 lncRNAs derived from 10,883 loci were identified, including 8954 novel lncRNAs. Differential expression of lncRNAs was observed in the developing shoot apices and grains, the two organs that have a direct influence on the final yield. The regulatory interaction of differentially expressed lncRNAs with the potential target genes was evaluated. We identified 176 cis-acting lncRNAs in shoot apices and 424 in grains, while the number of trans-acting lncRNAs in these organs was 1736 and 540, respectively. The potential target protein-coding genes were identified, and their biological function was annotated using MapMan ontology. This is the first insight into the roles of lncRNAs in barley development on the genome-wide scale, and our results provide a solid background for future functional studies.     

Genome-Wide Association Analysis and Genomic Prediction of Thyroglobulin Plasma Levels

Thyroglobulin (Tg) is an iodoglycoprotein produced by thyroid follicular cells which acts as an essential substrate for thyroid hormone synthesis. To date, only one genome-wide association study (GWAS) of plasma Tg levels has been performed by our research group. Utilizing recent advancements in computation and modeling, we apply a Bayesian approach to the probabilistic inference of the genetic architecture of Tg. We fitted a Bayesian sparse linear mixed model (BSLMM) and a frequentist linear mixed model (LMM) of 7,289,083 variants in 1096 healthy European-ancestry participants of the Croatian Biobank. Meta-analysis with two independent cohorts (total n = 2109) identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) within the ST6GAL1 gene (p<5×10−8). BSLMM revealed additional association signals on chromosomes 1, 8, 10, and 14. For ST6GAL1 and the newly uncovered genes, we provide physiological and pathophysiological explanations of how their expression could be associated with variations in plasma Tg levels. We found that the SNP-heritability of Tg is 17% and that 52% of this variation is due to a small number of 16 variants that have a major effect on Tg levels. Our results suggest that the genetic architecture of plasma Tg is not polygenic, but influenced by a few genes with major effects.