Genome-Wide Association Study Statistical Models: A Review

Genome-Wide Association Study as an Efficacious Approach to Discover Candidate Genes Associated with Body Linear Type Traits in Dairy Cattle

Body shape traits are very important and play a crucial role in the economic development of dairy farming. By improving the accuracy of selection for body size traits, we can enhance economic returns across the dairy industry and on farms, contributing to the future profitability of the dairy sector. Registered body conformation traits are reliable and cost-effective tools for use in national cattle breeding selection programs. These traits are significantly related to the production, longevity, mobility, health, fertility, and environmental adaptation of dairy cows. Therefore, they can be considered indirect indicators of economically important traits in dairy cows. Utilizing efficacious genetic methods, such as genome-wide association studies (GWASs), allows for a deeper understanding of the genetic architecture of complex traits through the identification and application of genetic markers. In the current review, we summarize information on candidate genes and genomic regions associated with body conformation traits in dairy cattle worldwide. The manuscript also reviews the importance of body conformation, the relationship between body conformation traits and other traits, heritability, influencing factors, and the genetics of body conformation traits. The information on candidate genes related to body conformation traits provided in this review may be helpful in selecting potential genetic markers for the genetic improvement of body conformation traits in dairy cattle.

Exploring the Interplay between the Hologenome and Complex Traits in Bovine and Porcine Animals Using Genome-Wide Association Analysis

Genome-wide association studies (GWAS) significantly enhance our ability to identify trait-associated genomic variants by considering the host genome. Moreover, the hologenome refers to the host organism's collective genetic material and its associated microbiome. In this study, we utilized the hologenome framework, called Hologenome-wide association studies (HWAS), to dissect the architecture of complex traits, including milk yield, methane emissions, rumen physiology in cattle, and gut microbial composition in pigs. We employed four statistical models: (1) GWAS, (2) Microbial GWAS (M-GWAS), (3) HWAS-CG (hologenome interaction estimated using COvariance between Random Effects Genome-based restricted maximum likelihood (CORE-GREML)), and (4) HWAS-H (hologenome interaction estimated using the Hadamard product method). We applied Bonferroni correction to interpret the significant associations in the complex traits. The GWAS and M-GWAS detected one and sixteen significant SNPs for milk yield traits, respectively, whereas the HWAS-CG and HWAS-H each identified eight SNPs. Moreover, HWAS-CG revealed four, and the remaining models identified three SNPs each for methane emissions traits. The GWAS and HWAS-CG detected one and three SNPs for rumen physiology traits, respectively. For the pigs' gut microbial composition traits, the GWAS, M-GWAS, HWAS-CG, and HWAS-H identified 14, 16, 13, and 12 SNPs, respectively. We further explored these associations through SNP annotation and by analyzing biological processes and functional pathways. Additionally, we integrated our GWA results with expression quantitative trait locus (eQTL) data using transcriptome-wide association studies (TWAS) and summary-based Mendelian randomization (SMR) methods for a more comprehensive understanding of SNP-trait associations. Our study revealed hologenomic variability in agriculturally important traits, enhancing our understanding of host-microbiome interactions.

Genome-wide association study and expression of candidate genes for Fe and Zn concentration in sorghum grains

Sorghum germplasm showed grain Fe and Zn genetic variability, but a few varieties were biofortified with these minerals. This work contributes to narrowing this gap. Fe and Zn concentrations along with 55,068 high-quality GBS SNP data from 140 sorghum accessions were used in this study. Both micronutrients exhibited good variability with respective ranges of 22.09-52.55 ppm and 17.92-43.16 ppm. Significant marker-trait associations were identified on chromosomes 1, 3, and 5. Two major effect SNPs (S01_72265728 and S05_58213541) explained 35% and 32% of Fe and Zn phenotypic variance, respectively. The SNP S01_72265728 was identified in the cytochrome P450 gene and showed a positive effect on Fe accumulation in the kernel, while S05_58213541 was intergenic near Sobic.005G134800 (zinc-binding ribosomal protein) and showed negative effect on Zn. Tissue-specific in silico expression analysis resulted in higher levels of Sobic.003G350800 gene product in several tissues such as leaf, root, flower, panicle, and stem. Sobic.005G188300 and Sobic.001G463800 were expressed moderately at grain maturity and anthesis in leaf, root, panicle, and seed tissues. The candidate genes expressed in leaves, stems, and grains will be targeted to improve grain and stover quality. The haplotypes identified will be useful in forward genetics breeding.

Genetic study for seed germination and shattering in Euphorbia lagascae in response to different seed treatments

Euphorbia lagascae Spreng is a promising emerging oilseed crop, with its seed oil accounting for approximately 50% of the seed weight. Euphorbia oil contains a significant amount of vernolic acid, comprising two-thirds of its composition, which boasts various industrial applications, including acting as a stabilizer-plasticizer and natural dye. However, this species was known to have a high degree of seed-shattering and a low germination rate, which act as two important barriers to large-scale production and exploitation. Therefore, the objective of this study is to determine the genetic control of seed germination and seed-shattering traits in order to develop a reliable pipeline that would be applicable for industries and breeders to select superior E. lagascae lines and design a robust breeding scheme in a short time at reduced labor costs. For this objective, five different wild-type genotypes of E. lagascae that demonstrated high germination potential were crossed with an ethyl methanesulfonate (EMS) mutant genotype that produces non-shattering capsules. The F2 populations from two successful crosses (A and B) were separated into three different treated groups for seed germination evaluation and to study the segregation of 200 individuals per F2 population. The three treatments were: light, gibberellic acid (GA3), and control treatment. Consequently, plants treated with approximately 250 μmol/m2/s of light showed significant improvement in germination up to 75% in cross A and 82.4 % in cross B compared with the control plants and the group treated with 0.05% GA3. According to the chi-square test results, the inheritance pattern of seed germination in response to light treatment follows a 3:1 segregation ratio between germinated and non-germinated seeds, indicating a dominant gene action in the F2 generation. The same conclusion was followed for the shattering trait in the group treated with light, which was also simply inherited as a 3:1 ratio for shattering vs. non-shattering capsules. Our results emphasize the importance and significance of light treatment in producing uniform populations through acceptable germination and shattering resistance of the mutant genotypes of E. lagascae. This is the first report of light treatment that significantly improved seed germination of E. lagascae, which may enhance efforts in the development of this new industrial crop as a feedstock for vernolic acid production.

Genome-wide association analysis identifies a candidate gene controlling seed size and yield in Xanthoceras sorbifolium Bunge

Yellow horn (Xanthoceras sorbifolium Bunge) is a woody oilseed tree species whose seed oil is rich in unsaturated fatty acids and rare neuronic acids, and can be used as a high-grade edible oil or as a feedstock for biodiesel production. However, the genetic mechanisms related to seed yield in yellow horn are not well elucidated. This study identified 2 164 863 SNP loci based on 222 genome-wide resequencing data of yellow horn germplasm. We conducted genome-wide association study (GWAS) analysis on three core traits (hundred-grain weight, single-fruit seed mass, and single-fruit seed number) that influence seed yield for the years 2022 and 2020, and identified 399 significant SNP loci. Among these loci, the Chr10_24013014 and Chr10_24012613 loci caught our attention due to their consistent associations across multiple analyses. Through Sanger sequencing, we validated the genotypes of these two loci across 16 germplasms, confirming their consistency with the GWAS analysis results. Downstream of these two significant loci, we identified a candidate gene encoding an AP2 transcription factor protein, which we named XsAP2. RT-qPCR analysis revealed high expression of the XsAP2 gene in seeds, and a significant negative correlation between its expression levels and seed hundred-grain weight, as well as single-fruit seed mass, suggesting its potential role in the normal seed development process. Transgenic Arabidopsis lines with the overexpressed XsAP2 gene exhibited varying degrees of reduction in seed size, number of seeds per silique, and number of siliques per plant compared with wild-type Arabidopsis. Combining these results, we hypothesize that the XsAP2 gene may have a negative regulatory effect on seed yield of yellow horn. These results provide a reference for the molecular breeding of high-yielding yellow horn.

Editorial for the Special Issue "Genetics Studies on Wheat"

Wheat (Triticum aestivum L [...].

Genome-wide analysis of cold imbibition stress in soybean, Glycine max

In Canada, the length of the frost-free season necessitates planting crops as early as possible to ensure that the plants have enough time to reach full maturity before they are harvested. Early planting carries inherent risks of cold water imbibition (specifically less than 4°C) affecting seed germination. A marker dataset developed for a previously identified Canadian soybean GWAS panel was leveraged to investigate the effect of cold water imbibition on germination. Seed from a panel of 137 soybean elite cultivars, grown in the field at Ottawa, ON, over three years, were placed on filter paper in petri dishes and allowed to imbibe water for 16 hours at either 4°C or 20°C prior to being transferred to a constant 20°C. Observations on seed germination, defined as the presence of a 1 cm radicle, were done from day two to seven. A three-parameter exponential rise to a maximum equation (3PERM) was fitted to estimate germination, time to the one-half maximum germination, and germination uniformity for each cultivar. Genotype-by-sequencing was used to identify SNPs in 137 soybean lines, and using genome-wide association studies (GWAS - rMVP R package, with GLM, MLM, and FarmCPU as methods), haplotype block analysis, and assumed linkage blocks of ±100 kbp, a threshold for significance was established using the qvalue package in R, and five significant SNPs were identified on chromosomes 1, 3, 4, 6, and 13 for maximum germination after cold water imbibition. Percent of phenotypic variance explained (PVE) and allele substitution effect (ASE) eliminated two of the five candidate SNPs, leaving three QTL regions on chromosomes 3, 6, and 13 (Chr3-3419152, Chr6-5098454, and Chr13-29649544). Based on the gene ontology (GO) enrichment analysis, 14 candidate genes whose function is predicted to include germination and cold tolerance related pathways were identified as candidate genes. The identified QTLs can be used to select future soybean cultivars tolerant to cold water imbibition and mitigate risks associated with early soybean planting.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

Unveiling the Mysteries of Non-Mendelian Heredity in Plant Breeding

Mendelian heredity is the cornerstone of plant breeding and has been used to develop new varieties of plants since the 19th century. However, there are several breeding cases, such as cytoplasmic inheritance, methylation, epigenetics, hybrid vigor, and loss of heterozygosity (LOH), where Mendelian heredity is not applicable, known as non-Mendelian heredity. This type of inheritance can be influenced by several factors besides the genetic architecture of the plant and its breeding potential. Therefore, exploring various non-Mendelian heredity mechanisms, their prevalence in plants, and the implications for plant breeding is of paramount importance to accelerate the pace of crop improvement. In this review, we examine the current understanding of non-Mendelian heredity in plants, including the mechanisms, inheritance patterns, and applications in plant breeding, provide an overview of the various forms of non-Mendelian inheritance (including epigenetic inheritance, cytoplasmic inheritance, hybrid vigor, and LOH), explore insight into the implications of non-Mendelian heredity in plant breeding, and the potential it holds for future research.

Multi-locus genome-wide association studies reveal genomic regions and putative candidate genes associated with leaf spot diseases in African groundnut (Arachis hypogaea L.) germplasm

Early leaf spot (ELS) and late leaf spot (LLS) diseases are the two most destructive groundnut diseases in Ghana resulting in ≤ 70% yield losses which is controlled largely by chemical method. To develop leaf spot resistant varieties, the present study was undertaken to identify single nucleotide polymorphism (SNP) markers and putative candidate genes underlying both ELS and LLS. In this study, six multi-locus models of genome-wide association study were conducted with the best linear unbiased predictor obtained from 294 African groundnut germplasm screened for ELS and LLS as well as image-based indices of leaf spot diseases severity in 2020 and 2021 and 8,772 high-quality SNPs from a 48 K SNP array Axiom platform. Ninety-seven SNPs associated with ELS, LLS and five image-based indices across the chromosomes in the 2 two sub-genomes. From these, twenty-nine unique SNPs were detected by at least two models for one or more traits across 16 chromosomes with explained phenotypic variation ranging from 0.01 - 62.76%, with exception of chromosome (Chr) 08 (Chr08), Chr10, Chr11, and Chr19. Seventeen potential candidate genes were predicted at ± 300 kbp of the stable/prominent SNP positions (12 and 5, down- and upstream, respectively). The results from this study provide a basis for understanding the genetic architecture of ELS and LLS diseases in African groundnut germplasm, and the associated SNPs and predicted candidate genes would be valuable for breeding leaf spot diseases resistant varieties upon further validation.