Genome-wide Association Studies Methods Research Articles

Comparison of the statistical methods for genome-wide association studies on simulated quantitative traits of domesticated goats (Capra hircus L.)

Genome-wide association studies (GWAS) are one of the essential methods used to find associations between phenotypes and genes. The main problem in GWAS methods is false positives that occur from family relationships and population structure. In this study, the performances and powers of traditional and state-of-the-art statistical methods used in genome-wide association studies were compared on nine simulated quantitative traits using genotypic data of the domesticated goats (Capra hircus L.). The traits were simulated with varying degrees of heritability from 0.05 to 0.8. The findings revealed that Blink and FarmCPU are the outperforming methods to control false positives more efficiently for the medium and high heritability traits. At the same time, Blink promises to identify true positive signals more efficiently for the low heritability traits. The other multi-loci statistical methods SUPER and MLMM perform worse than Blink and FarmCPU but better than all of the single-locus methods GLM, MLM, CMLM and ECLMM. Additionally, Bonferroni or Holm's adjustment procedures were better in multiple testing of associations at the 1 % significance level to get more accurate results for the high heritability traits.

Application of SVR-Mediated GWAS for Identification of Durable Genetic Regions Associated with Soybean Seed Quality Traits.

Soybean (Glycine max L.) is an important food-grade strategic crop worldwide because of its high seed protein and oil contents. Due to the negative correlation between seed protein and oil percentage, there is a dire need to detect reliable quantitative trait loci (QTL) underlying these traits in order to be used in marker-assisted selection (MAS) programs. Genome-wide association study (GWAS) is one of the most common genetic approaches that is regularly used for detecting QTL associated with quantitative traits. However, the current approaches are mainly focused on estimating the main effects of QTL, and, therefore, a substantial statistical improvement in GWAS is required to detect associated QTL considering their interactions with other QTL as well. This study aimed to compare the support vector regression (SVR) algorithm as a common machine learning method to fixed and random model circulating probability unification (FarmCPU), a common conventional GWAS method in detecting relevant QTL associated with soybean seed quality traits such as protein, oil, and 100-seed weight using 227 soybean genotypes. The results showed a significant negative correlation between soybean seed protein and oil concentrations, with heritability values of 0.69 and 0.67, respectively. In addition, SVR-mediated GWAS was able to identify more relevant QTL underlying the target traits than the FarmCPU method. Our findings demonstrate the potential use of machine learning algorithms in GWAS to detect durable QTL associated with soybean seed quality traits suitable for genomic-based breeding approaches. This study provides new insights into improving the accuracy and efficiency of GWAS and highlights the significance of using advanced computational methods in crop breeding research.

Using DEPendency of Association on the Number of Top Hits (DEPTH) as a Complementary Tool to Identify Novel Colorectal Cancer Susceptibility Loci.

DEPendency of association on the number of Top Hits (DEPTH) is an approach to identify candidate susceptibility regions by considering the risk signals from overlapping groups of sequential variants across the genome. We applied a DEPTH analysis using a sliding window of 200 SNPs to colorectal cancer data from the Colon Cancer Family Registry (CCFR; 5,735 cases and 3,688 controls), and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO; 8,865 cases and 10,285 controls) studies. A DEPTH score > 1 was used to identify candidate susceptibility regions common to both analyses. We compared DEPTH results against those from conventional genome-wide association study (GWAS) analyses of these two studies as well as against 132 published susceptibility regions. Initial DEPTH analysis revealed 2,622 (CCFR) and 3,686 (GECCO) candidate susceptibility regions, of which 569 were common to both studies. Bootstrapping revealed 40 and 49 candidate susceptibility regions in the CCFR and GECCO data sets, respectively. Notably, DEPTH identified at least 82 regions that would not be detected using conventional GWAS methods, nor had they been identified by previous colorectal cancer GWASs. We found four reproducible candidate susceptibility regions (2q22.2, 2q33.1, 6p21.32, 13q14.3). The highest DEPTH scores were in the human leukocyte antigen locus at 6p21 where the strongest associated SNPs were rs762216297, rs149490268, rs114741460, and rs199707618 for the CCFR data, and rs9270761 for the GECCO data. DEPTH can identify candidate susceptibility regions for colorectal cancer not identified using conventional analyses of larger datasets. DEPTH has potential as a powerful complementary tool to conventional GWAS analyses for discovering susceptibility regions within the genome.

Genome-wide study for SNPS associated with average daily growth for genotyping by sequencing and identification of candidate genes in jalgin merino

To improve the productive qualities of sheep bred in sheep farms, an important factor is the selection activities aimed at selecting animals for economically valuable traits. For this purpose, methods of markerassociated selection based on genotyping data are widely used. Genome-wide association study (GWAS) is used to identify single nucleotide substitutions suitable for genotyping by sequencing and to discover new candidate genes for productive traits. This work is devoted to the search for genome loci associated with the "average daily gain" trait using the GWAS method in Jalgin Merino for further use in breeding. A genome-wide study was carried out using Illumina Ovine Infinium HD BeadChip 600K on rams aged 12 months, bred in the Ipatovsky district of the Stavropol Territory. Quality control of genotyping, as well as genetic and statistical analysis of the obtained data, was performed using PLINK V.1.07. Visualization and plotting were carried out using the QQman package in the R programming language. A genome- wide study for associations showed a significant relationship between 10 detected SNPs and the average daily gain. Analysis of the localization of single nucleotide substitutions found 5 genes next to them: BABAM2, ELL2, ENSOARG00000019376, PLA2G10, ENSOARG00000003785. Their protein products are involved in the implementation of genetic information, cell differentiation and the flow of important metabolic processes in the body. Further research is needed to study the structure of the proposed candidate genes in order to understand the mechanisms of their influence on the phenotypic manifestation of economically valuable traits in sheep of domestic breeds.

Multi-layered genetic approaches to identify approved drug targets

Drugs targeting genes linked to disease via evidence from human genetics have increased odds of approval. Approaches to prioritize such genes include genome-wide association studies (GWASs), rare variant burden tests in exome sequencing studies (Exome), or integration of a GWAS with expression/protein quantitative trait loci (eQTL/pQTL-GWAS). Here, we compare gene-prioritization approaches on 30 clinically relevant traits and benchmark their ability to recover drug targets. Across traits, prioritized genes were enriched for drug targets with odds ratios (ORs) of 2.17, 2.04, 1.81, and 1.31 for the GWAS, eQTL-GWAS, Exome, and pQTL-GWAS methods, respectively. Adjusting for differences in testable genes and sample sizes, GWAS outperforms e/pQTL-GWAS, but not the Exome approach. Furthermore, performance increased through gene network diffusion, although the node degree, being the best predictor (OR= 8.7), revealed strong bias in literature-curated networks. In conclusion, we systematically assessed strategies to prioritize drug target genes, highlighting the promises and pitfalls of current approaches.

Re-sequencing and morphological data revealed the genetics of stone shell and kernel traits in apricot.

Kernel-using apricot (Prunus armeniaca L.) is an economically important fruit tree species in arid areas owing to its hardiness and cold and drought tolerance. However, little is known about its genetic background and trait inheritances. In the present study, we first evaluated the population structure of 339 apricot accessions and the genetic diversity of kernel-using apricots using whole genome re-sequencing. Second, the phenotypic data of 222 accessions were investigated for two consecutive seasons (2019 and 2020) for 19 traits, including kernel and stone shell traits and the pistil abortion rate of flowers. Heritability and correlation coefficient of traits were also estimated. The stone shell length (94.46%) showed the highest heritability, followed by the length/width ratio (92.01%) and length/thickness ratio (92.00%) of the stone shell, whereas breaking force of the nut (17.08%) exhibited a very low heritability. A genome-wide association study (GWAS) using general linear model and generalized linear mixed model revealed 122 quantitative trait loci (QTLs). The QTLs of the kernel and stone shell traits were unevenly assigned on the eight chromosomes. Out of the 1,614 candidate genes identified in the 13 consistently reliable QTLs found using the two GWAS methods and/or in the two seasons, 1,021 were annotated. The sweet kernel trait was assigned to chromosome 5 of the genome, similar to the almond, and a new locus was also mapped at 17.34-17.51 Mb on chromosome 3, including 20 candidate genes. The loci and genes identified here will be of significant use in molecular breeding efforts, and the candidate genes could play essential roles in exploring the mechanisms of genetic regulation.

Isogenic human trophectoderm cells demonstrate the role of NDUFA4 and associated variants in ZIKV infection

Population-based genome-wide association studies (GWAS) normally require a large sample size, which can be labor intensive and costly. Recently, we reported a human induced pluripotent stem cell (hiPSC) array-based GWAS method, identifying NDUFA4 as a host factor for Zika virus (ZIKV) infection. In this study, we extended our analysis to trophectoderm cells, which constitute one of the major routes of mother-to-fetus transmission of ZIKV during pregnancy. We differentiated hiPSCs from various donors into trophectoderm cells. We then infected cells carrying loss of function mutations in NDUFA4, harboring risk versus non-risk alleles of SNPs (rs917172 and rs12386620) or having deletions in the NDUFA4 cis-regulatory region with ZIKV. We found that loss/reduction of NDUFA4 suppressed ZIKV infection in trophectoderm cells. This study validated our published hiPSC array-based system as a useful platform for GWAS and confirmed the role of NDUFA4 as a susceptibility locus for ZIKV in disease-relevant trophectoderm cells.

Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep.

Economic losses due to copper intoxication or deficiency is a problem encountered by sheep farmers. The aim of this study was to investigate the ovine genome for genomic regions and candidate genes responsible for variability in liver copper concentration. Liver samples were collected from slaughtered lambs of the Merinoland breed from two farms, and used for measurement of copper concentration and genome-wide association study (GWAS). A total of 45,511 SNPs and 130 samples were finally used for analysis, in which single-locus and several multi-locus GWAS (SL-GWAS; ML-GWAS) methods were employed. Gene enrichment analysis was performed for identified candidate genes to detect gene ontology (GO) terms significantly associated with hepatic copper levels. The SL-GWAS and a minimum of two ML-GWAS identified two and thirteen significant SNPs, respectively. Within genomic regions surrounding identified SNPs, we observed nine promising candidate genes such as DYNC1I2, VPS35, SLC38A9 and CHMP1A. GO terms such as lysosomal membrane, mitochondrial inner membrane and sodium:proton antiporter activity were significantly enriched. Genes involved in these identified GO terms mediate multivesicular body (MVB) fusion with lysosome for degradation and control mitochondrial membrane permeability. This reveals the polygenic status of this trait and candidate genes for further studies on breeding for copper tolerance in sheep.

Kernel-based genetic association analysis for microbiome phenotypes identifies host genetic drivers of beta-diversity

BackgroundUnderstanding human genetic influences on the gut microbiota helps elucidate the mechanisms by which genetics may influence health outcomes. Typical microbiome genome-wide association studies (GWAS) marginally assess the association between individual genetic variants and individual microbial taxa. We propose a novel approach, the covariate-adjusted kernel RV (KRV) framework, to map genetic variants associated with microbiome beta-diversity, which focuses on overall shifts in the microbiota. The KRV framework evaluates the association between genetics and microbes by comparing similarity in genetic profiles, based on groups of variants at the gene level, to similarity in microbiome profiles, based on the overall microbiome composition, across all pairs of individuals. By reducing the multiple-testing burden and capturing intrinsic structure within the genetic and microbiome data, the KRV framework has the potential of improving statistical power in microbiome GWAS.ResultsWe apply the covariate-adjusted KRV to the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) in a two-stage (first gene-level, then variant-level) genome-wide association analysis for gut microbiome beta-diversity. We have identified an immunity-related gene, IL23R, reported in a previous microbiome genetic association study and discovered 3 other novel genes, 2 of which are involved in immune functions or autoimmune disorders. In addition, simulation studies show that the covariate-adjusted KRV has a greater power than other microbiome GWAS methods that rely on univariate microbiome phenotypes across a range of scenarios.ConclusionsOur findings highlight the value of the covariate-adjusted KRV as a powerful microbiome GWAS approach and support an important role of immunity-related genes in shaping the gut microbiome composition.3iRoYfA_uxPN1-JktGz5GLVideo

Genome-wide association study in two-row spring barley landraces identifies QTL associated with plantlets root system architecture traits in well-watered and osmotic stress conditions.

Water availability is undoubtedly one of the most important environmental factors affecting crop production. Drought causes a gradual deprivation of water in the soil from top to deep layers and can occur at diverse stages of plant development. Roots are the first organs that perceive water deficit in soil and their adaptive development contributes to drought adaptation. Domestication has contributed to a bottleneck in genetic diversity. Wild species or landraces represent a pool of genetic diversity that has not been exploited yet in breeding program. In this study, we used a collection of 230 two-row spring barley landraces to detect phenotypic variation in root system plasticity in response to drought and to identify new quantitative trait loci (QTL) involved in root system architecture under diverse growth conditions. For this purpose, young seedlings grown for 21 days in pouches under control and osmotic-stress conditions were phenotyped and genotyped using the barley 50k iSelect SNP array, and genome-wide association studies (GWAS) were conducted using three different GWAS methods (MLM GAPIT, FarmCPU, and BLINK) to detect genotype/phenotype associations. In total, 276 significant marker-trait associations (MTAs; p-value (FDR)< 0.05) were identified for root (14 and 12 traits under osmotic-stress and control conditions, respectively) and for three shoot traits under both conditions. In total, 52 QTL (multi-trait or identified by at least two different GWAS approaches) were investigated to identify genes representing promising candidates with a role in root development and adaptation to drought stress.

Revealing the correlations between brain cortical characteristics and susceptibility genes for Alzheimer disease: a cross-sectional study.

Alzheimer disease (AD) is a progressive neurodegenerative disease closely related to genes and characterized by the atrophy of the cerebral cortex. Correlations between imaging phenotypes and the susceptibility genes for AD, as demonstrated in the findings of genome-wide association studies (GWASs), still need to be addressed due to the complicated structure of the human cortex. In our study, an improved GWAS method, whole cortex characteristics GWAS (WCC-GWAS), was proposed. The WCC-GWAS uses multiple cortex characteristics of gray-matter volume (GMV), cortical thickness (CT), cortical surface area (CSA), and local gyrification index (LGI). A cohort of 496 participants was enrolled and divided into 4 groups: normal control (NC; n=122), early mild cognitive impairment (EMCI; n=196), late mild cognitive impairment (LMCI; n=62), and AD (n=116). Based on the Desikan-Killiany atlas, the brain was parcellated into 68 brain regions, and the WCC of each brain region was individually calculated. Four cortex characteristics of GMV, CT, CSA, and LGI across the 4 groups optimized with multiple comparisons and the ReliefF algorithm were taken as magnetic resonance imaging (MRI) brain phenotypes. Under the model of multiple linear additive genetic regression, the correlations between the MRI brain phenotypes and single-nucleotide polymorphisms (SNPs) were deduced. The findings identified 2 prominent correlations. First, rs7309929 of neuron navigator 3 (NAV3) located on chromosome 12 correlated with the decreased GMV for the left middle temporal gyrus (P=0.0074). Second, rs11250992 of long intergenic non-protein-coding RNA 700 (LINC00700) located on chromosome 10 correlated with the decreased CT for the left supramarginal gyrus (P=0.0019). The findings suggested that the correlations between phenotypes and genotypes could be effectively evaluated. The strategy of extracting MRI phenotypes as endophenotypes provided valuable indications in AD GWAS.

Multivariate GWAS analysis reveals loci associated with liver functions in continental African populations.

Liver disease is any condition that causes liver damage and inflammation and may likely affect the function of the liver. Vital biochemical screening tools that can be used to evaluate the health of the liver and help diagnose, prevent, monitor, and control the development of liver disease are known as liver function tests (LFT). LFTs are performed to estimate the level of liver biomarkers in the blood. Several factors are associated with differences in concentration levels of LFTs in individuals, such as genetic and environmental factors. The aim of our study was to identify genetic loci associated with liver biomarker levels with a shared genetic basis in continental Africans, using a multivariate genome-wide association study (GWAS) approach. We used two distinct African populations, the Ugandan Genome Resource (UGR = 6,407) and South African Zulu cohort (SZC = 2,598). The six LFTs used in our analysis were: aspartate transaminase (AST), alanine transaminase (ALT), alkaline phosphatase (ALP), gamma-glutamyl transferase (GGT), total bilirubin, and albumin. A multivariate GWAS of LFTs was conducted using the exact linear mixed model (mvLMM) approach implemented in GEMMA and the resulting P-values were presented in Manhattan and quantile-quantile (QQ) plots. First, we attempted to replicate the findings of the UGR cohort in SZC. Secondly, given that the genetic architecture of UGR is different from that of SZC, we further undertook similar analysis in the SZC and discussed the results separately. A total of 59 SNPs reached genome-wide significance (P = 5x10-8) in the UGR cohort and with 13 SNPs successfully replicated in SZC. These included a novel lead SNP near the RHPN1 locus (lead SNP rs374279268, P-value = 4.79x10-9, Effect Allele Frequency (EAF) = 0.989) and a lead SNP at the RGS11 locus (lead SNP rs148110594, P-value = 2.34x10-8, EAF = 0.928). 17 SNPs were significant in the SZC, while all the SNPs fall within a signal on chromosome 2, rs1976391 mapped to UGT1A was identified as the lead SNP within this region. Using multivariate GWAS method improves the power to detect novel genotype-phenotype associations for liver functions not found with the standard univariate GWAS in the same dataset.

Utilizing Variants Identified with Multiple Genome-Wide Association Study Methods Optimizes Genomic Selection for Growth Traits in Pigs

Improving the prediction accuracies of economically important traits in genomic selection (GS) is a main objective for researchers and breeders in the livestock industry. This study aims at utilizing potentially functional SNPs and QTLs identified with various genome-wide association study (GWAS) models in GS of pig growth traits. We used three well-established GWAS methods, including the mixed linear model, Bayesian model and meta-analysis, as well as 60K SNP-chip and whole genome sequence (WGS) data from 1734 Yorkshire and 1123 Landrace pigs to detect SNPs related to four growth traits: average daily gain, backfat thickness, body weight and birth weight. A total of 1485 significant loci and 24 candidate genes which are involved in skeletal muscle development, fatty deposition, lipid metabolism and insulin resistance were identified. Compared with using all SNP-chip data, GS with the pre-selected functional SNPs in the standard genomic best linear unbiased prediction (GBLUP), and a two-kernel based GBLUP model yielded average gains in accuracy by 4 to 46% (from 0.19 ± 0.07 to 0.56 ± 0.07) and 5 to 27% (from 0.16 ± 0.06 to 0.57 ± 0.05) for the four traits, respectively, suggesting that the prioritization of preselected functional markers in GS models had the potential to improve prediction accuracies for certain traits in livestock breeding.

Genome-Wide Association Studies (GWAS).

Most of the breeding targets are quantitative traits. In exploring the quantitative trait locus (QTL) system of a trait, linkage mapping was established using sparse polymerase chain reaction (PCR) markers. With the genome-wide sequencing technology advanced, genome-wide association study (GWAS) was developed for natural (germplasm) populations using dense genomic markers, which facilitates the identification of the complete QTL system with their multiple alleles on genomic locations. GWAS makes use of the linkage disequilibrium (LD) due to historical saturate recombination and high-density genomic markers to detect QTLs through statistical test for the association between molecular markers and phenotypes. However, due to inbreeding and mixture of source populations, the germplasm population often has complex and unknown structure, which leads to false positives/negatives in GWAS. Various GWAS methods have been proposed to reduce false positives/negatives, including those of the general linear model and the mixed linear model, which focused mainly on finding a handful of major QTLs under single-locus model for major gene cloning and could not detect directly the multiple alleles using bi-allelic single-nucleotide polymorphism (SNP) marker. As a relatively thorough detection of QTLs with their multiple alleles is required for germplasm population, the restricted two-stage multi-locus multi-allele GWAS (RTM-GWAS) procedure was proposed for identifying the QTL system with varying multiple alleles. From the RTM-GWAS results, a QTL-allele matrix is constructed as a compact form of the population genetic constitution, which can be further used for crop genetic and breeding studies, including major gene mining, population evolution, and breeding by genetic design.

Genome-wide association studies on collagen contents trait for meat quality in Hanwoo.

Beef consumers valued meat quality traits such as texture, tenderness, juiciness, flavor, and meat color that determining consumers' purchasing decision. Most research on meat quality has focused on marbling, a key characteristic related to meat eating quality. However, other important traits such as meat texture, tenderness, and color have not much studied in cattle. Among these traits, meat tenderness and texture of cattle are among the most important factors affecting quality evaluation of consumers. Collagen is the main component of connective tissues.It greatly affects meat tenderness. The objective of this study was to determine significant variants and candidate genes associated with collagen contents trait (total collagen) through genome-wide association studies (GWAS). Phenotypic and genomic data from 135 Hanwoo were used. The BLUPF90 family program and GRAMMAR method for GWAS were applied in this study. A total of 73 potential single nucleotide polymorphisms (SNPs) showed significant associations with collagen content. They were located in or near 108 candidate genes. TMEM135 and ME3 genes were identified to have the most significant SNPs associated with collagen contents trait. Data indicated that these genes were related to collagen. Biological processes and pathways for the prediction of biological functions of candidate genes were confirmed. We found that candidate genes were involved in positive regulation of CREB transcription factor activity and actin cytoskeleton related to tenderness and texture of beef. Three genes (CRTC3, MYO1C and MYLK4) belonging to these biological functions were related to tenderness. These results provide a basis for improving genomic characteristics of Hanwoo for the production of tender beef. Furthermore, they could be used they could be used as an index to select desired traits for consumers.

Identifying the Genetic Basis of Mineral Elements in Rice Grain Using Genome-Wide Association Mapping.

Mineral malnutrition is a major problem in many rice-consuming countries. It is essential to know the genetic mechanisms of accumulation of mineral elements in the rice grain to provide future solutions for this issue. This study was conducted to identify the genetic basis of six mineral elements (Cu, Fe, K, Mg, Mn, and Zn) by using three models for single-locus and six models for multi-locus analysis of a genome-wide association study (GWAS) using 174 diverse rice accessions and 6565 SNP markers. To declare a SNP as significant, -log10(P) ≥ 3.0 and 15% FDR significance cut-off values were used for single-locus models, while LOD ≥ 3.0 was used for multi-locus models. Using these criteria, 147 SNPs were detected by one or two GWAS methods at -log10(P) ≥ 3.0, 48 of which met the 15% FDR significance cut-off value. Single-locus models outperformed multi-locus models before applying multi-test correction, but once applied, multi-locus models performed better. While 14 (~29%) of the identified quantitative trait loci (QTLs) after multiple test correction co-located with previously reported genes/QTLs and marker associations, another 34 trait-associated SNPs were novel. After mining genes within 250 kb of the 48 significant SNP loci, in silico and gene enrichment analyses were conducted to predict their potential functions. These shortlisted genes with their functions could guide future experimental validation, helping us to understand the complex molecular mechanisms controlling rice grain mineral elements.

Genome-wide association studies of five free amino acid levels in rice.

Rice (Oryza sativa L.) is one of the important staple foods for human consumption and livestock use. As a complex quality trait, free amino acid (FAA) content in rice is of nutritional importance. To dissect the genetic mechanism of FAA level, five amino acids' (Val, Leu, Ile, Arg, and Trp) content and 4,325,832 high-quality SNPs of 448 rice accessions were used to conduct genome-wide association studies (GWAS) with nine different methods. Of these methods, one single-locus method (GEMMA), seven multi-locus methods (mrMLM, pLARmEB, FASTmrEMMA, pKWmEB, FASTmrMLM, ISIS EM-BLASSO, and FarmCPU), and the recent released 3VmrMLM were adopted for methodological comparison of quantitative trait nucleotide (QTN) detection and identification of stable quantitative trait nucleotide loci (QTLs). As a result, 987 QTNs were identified by eight multi-locus GWAS methods; FASTmrEMMA detected the most QTNs (245), followed by 3VmrMLM (160), and GEMMA detected the least QTNs (0). Among 88 stable QTLs identified by the above methods, 3VmrMLM has some advantages, such as the most common QTNs, the highest LOD score, and the highest proportion of all detected stable QTLs. Around these stable QTLs, candidate genes were found in the GO classification to be involved in the primary metabolic process, biosynthetic process, and catalytic activity, and shown in KEGG analysis to have participated in metabolic pathways, biosynthesis of amino acids, and tryptophan metabolism. Natural variations of candidate genes resulting in the content alteration of five FAAs were identified in this association panel. In addition, 95 QTN-by-environment interactions (QEIs) of five FAA levels were detected by 3VmrMLM only. GO classification showed that the candidate genes got involved in the primary metabolic process, transport, and catalytic activity. Candidate genes of QEIs played important roles in valine, leucine, and isoleucine degradation (QEI_09_03978551 and candidate gene LOC_Os09g07830 in the Leu dataset), tryptophan metabolism (QEI_01_00617184 and candidate gene LOC_Os01g02020 in the Trp dataset), and glutathione metabolism (QEI_12_09153839 and candidate gene LOC_Os12g16200 in the Arg dataset) pathways through KEGG analysis. As an alternative of the multi-locus GWAS method, these findings suggested that the application of 3VmrMLM may provide new insights into better understanding FAA accumulation and facilitate the molecular breeding of rice with high FAA level.

Prioritized candidate causal haplotype blocks in plant genome-wide association studies.

Genome wide association studies (GWAS) can play an essential role in understanding genetic basis of complex traits in plants and animals. Conventional SNP-based linear mixed models (LMM) that marginally test single nucleotide polymorphisms (SNPs) have successfully identified many loci with major and minor effects in many GWAS. In plant, the relatively small population size in GWAS and the high genetic diversity found in many plant species can impede mapping efforts on complex traits. Here we present a novel haplotype-based trait fine-mapping framework, HapFM, to supplement current GWAS methods. HapFM uses genotype data to partition the genome into haplotype blocks, identifies haplotype clusters within each block, and then performs genome-wide haplotype fine-mapping to prioritize the candidate causal haplotype blocks of trait. We benchmarked HapFM, GEMMA, BSLMM, GMMAT, and BLINK in both simulated and real plant GWAS datasets. HapFM consistently resulted in higher mapping power than the other GWAS methods in high polygenicity simulation setting. Moreover, it resulted in smaller mapping intervals, especially in regions of high LD, achieved by prioritizing small candidate causal blocks in the larger haplotype blocks. In the Arabidopsis flowering time (FT10) datasets, HapFM identified four novel loci compared to GEMMA's results, and the average mapping interval of HapFM was 9.6 times smaller than that of GEMMA. In conclusion, HapFM is tailored for plant GWAS to result in high mapping power on complex traits and improved on mapping resolution to facilitate crop improvement.

A scalable Bayesian functional GWAS method accounting for multivariate quantitative functional annotations with applications for studying Alzheimer disease.

SummaryExisting methods for integrating functional annotations in genome-wide association studies (GWASs) to fine-map and prioritize potential causal variants are limited to using non-overlapped categorical annotations or limited by the computation burden of modeling genome-wide variants. To overcome these limitations, we propose a scalable Bayesian functional GWAS method to account for multivariate quantitative functional annotations (BFGWAS_QUANT), accompanied by a scalable computation algorithm enabling joint modeling of genome-wide variants. Simulation studies validated the performance of BFGWAS_QUANT for accurately quantifying annotation enrichment and improving GWAS power. Applying BFGWAS_QUANT to study five Alzheimer disease (AD)-related phenotypes using individual-level GWAS data (n = ∼1,000), we found that histone modification annotations have higher enrichment than expression quantitative trait locus (eQTL) annotations for all considered phenotypes, with the highest enrichment in H3K27me3 (polycomb regression). We also found that cis-eQTLs in microglia had higher enrichment than eQTLs of bulk brain frontal cortex tissue for all considered phenotypes. A similar enrichment pattern was also identified using the International Genomics of Alzheimer’s Project (IGAP) summary-level GWAS data of AD (n = ∼54,000). The strongest known APOE E4 risk allele was identified for all five phenotypes, and the APOE locus was validated using the IGAP data. BFGWAS_QUANT fine-mapped 32 significant variants from 1,073 genome-wide significant variants in the IGAP data. We also demonstrated that the polygenic risk scores (PRSs) using effect size estimates by BFGWAS_QUANT had a similar prediction accuracy as other methods assuming a sparse causal model. Overall, BFGWAS_QUANT is a useful GWAS tool for quantifying annotation enrichment and prioritizing potential causal variants.

First single-step genomic prediction and genome-wide association for body weight in Russian sturgeon (Acipenser gueldenstaedtii)

Russian sturgeon (Acipenser gueldenstaedtii)—a species of sturgeon—has a high economic value as a source of caviar and high-quality meat products. Genomic selection (GS) has the potential to shorten generation intervals and improve genetic gain for the key target traits of a breeding program. However, it has not been applied on sturgeon. We evaluated the performance of GS in breeding value estimation and carried out genome-wide association studies (GWAS) to dissect the genetic architecture of economically vital growth trait (body weight) in Russian sturgeon. Of a total of 3430 Russian sturgeons with phenotype records, 236 were sequenced. Through 10 replicates of fivefold cross-validation, we compared the accuracy of prediction obtained from a pedigree-based best linear unbiased prediction (BLUP), a single-step genomic BLUP (ssGBLUP), and a weighted ssGBLUP (WssGBLUP). We used single-step GWAS (ssGWAS) method to perform GWAS. The results showed that predictions made using ssGBLUP and WssGBLUP were 6.9% and 12.4% more accurate, respectively, compared with pedigree-based BLUP. Among the single-step methods, the performance of WssGBLUP was the best with the smallest mean squared (absolute) error. In addition, ssGWAS identified 3 potential SNPs with suggestive significance that were associated with body weight. Moreover, NOS1AP, CHMP4A, and RBMS1 genes were highlighted as functionally plausible candidate genes underlying the genetic architecture of body weight. These findings illustrate that the using GS in Russian sturgeon breeding could help in better selection of breeding candidates to facilitate complex economic growth trait.