A powerful approach to identify replicable variants in genome-wide association studies.

Replicability is the cornerstone of modern scientific research. Reliable identifications of genotype-phenotype associations that are significant in multiple genome-wide association studies (GWASs) provide stronger evidence for the findings. Current replicability analysis relies on the independence assumption among single-nucleotide polymorphisms (SNPs) and ignores the linkage disequilibrium (LD) structure. We show that such a strategy may produce either overly liberal or overly conservative results in practice. We develop an efficient method, ReAD, to detect replicable SNPs associated with the phenotype from two GWASs accounting for the LD structure. The local dependence structure of SNPs across two heterogeneous studies is captured by a four-state hidden Markov model (HMM) built on two sequences of p values. By incorporating information from adjacent locations via the HMM, our approach provides more accurate SNP significance rankings. ReAD is scalable, platform independent, and more powerful than existing replicability analysis methods with effective false discovery rate control. Through analysis of datasets from two asthma GWASs and two ulcerative colitis GWASs, we show that ReAD can identify replicable genetic loci that existing methods might otherwise miss.

Effect of serum metabolites on the risk of iridocyclitis: a bidirectional Mendelian randomization study.

Previous research has linked serum metabolite levels to iridocyclitis, yet their causal relationship remains unexplored. This study investigated this potential causality by analyzing pooled data from 7824 iridocyclitis patients in a Genome-Wide Association Study (GWAS) using Mendelian randomization (MR) and linkage disequilibrium score regression (LDSC). Employing rigorous quality control and comprehensive statistical methods, including sensitivity analyses, we examined the influence of 486 serum metabolites on iridocyclitis. Our MR analysis identified 23 metabolites with significant causal effects on iridocyclitis, comprising 17 known and 6 unidentified metabolites. Further refinement using Cochran's Q test and MR-PRESSO indicated 16 metabolites significantly associated with iridocyclitis risk. LDSC highlighted the heritability of certain metabolites, underscoring genetic influences on their levels. Notably, tryptophan, proline, theobromine, and 7-methylxanthine emerged as risk factors, while 3,4-dihydroxybutyrate appeared protective. These findings enhance our understanding of the metabolic interactions in iridocyclitis, offering insights for diagnosis, unraveling pathophysiological mechanisms, and informing potential avenues for prevention and personalized treatment.

Exploring the interconnected between type 2 diabetes mellitus and nonalcoholic fatty liver disease: Genetic correlation and Mendelian randomization analysis.

Epidemiological and clinical studies have indicated a higher risk of nonalcoholic fatty liver disease (NAFLD) and type 2 diabetes mellitus (T2DM), implying a potentially shared genetic etiology, which is still less explored. Genetic links between T2DM and NAFLD were assessed using linkage disequilibrium score regression and pleiotropic analysis under composite null hypothesis. European GWAS data have identified shared genes, whereas SNP-level pleiotropic analysis under composite null hypothesis has explored pleiotropic loci. generalized gene-set analysis of GWAS data determines pleiotropic pathways and tissue enrichment using eQTL mapping to identify associated genes. Mendelian randomization analysis was used to investigate the causal relationship between NAFLD and T2DM. Linkage disequilibrium score regression analysis revealed a strong genetic correlation between T2DM and NAFLD, and identified 24 pleiotropic loci. These single-nucleotide polymorphisms are primarily involved in biosynthetic regulation, RNA biosynthesis, and pancreatic development. generalized gene-set analysis of GWAS data analysis revealed significant enrichment in multiple brain tissues. Gene mapping using these 3 methods led to the identification of numerous pleiotropic genes, with differences observed in liver and kidney tissues. These genes were mainly enriched in pancreas, brain, and liver tissues. The Mendelian randomization method indicated a significantly positive unidirectional causal relationship between T2DM and NAFLD. Our study identified a shared genetic structure between NAFLD and T2DM, providing new insights into the genetic pathogenesis and mechanisms of NAFLD and T2DM comorbidities.

Dissecting the shared genetic architecture between endometriosis and polycystic ovary syndrome.

Background: Previous study suggested evidence for coexistence and similarities between endometriosis and polycystic ovary syndrome (PCOS), but it is unclear regarding the shared genetic architecture and causality underlying the phenotypic similarities observed for endometriosis and PCOS. Methods: By leveraging summary statistics from public genome-wide association studies regarding endometriosis (European-based: N=470,866) and PCOS (European-based: N=210,870), we explored the genetic correlation that shared between endometriosis and PCOS using linkage disequilibrium score regression. Shared risk SNPs were derived using PLACO (Pleiotropic analysis under composite null hypothesis) and FUMA (Functional Mapping and Annotation of Genetic Associations). The potential causal association between endometriosis and PCOS was investigated using two-sample Mendelian randomization (MR). Linkage disequilibrium score for the specific expression of genes analysis (LDSC-SEG) were performed for tissue enrichment analysis. The expression profiles of the risk gene in tissues were further examined. Results: A positive genetic association was observed between endometriosis and PCOS. 12 significant pleiotropic loci shared between endometriosis and PCOS were identified. Genetic associations between endometriosis and PCOS were particularly enriched in uterus, endometrium and fallopian tube. Two-sample MR analysis further indicated a potential causative effect of endometriosis on PCOS, and vice versa. Microarray and RNA-seq verified the expressions of SYNE1 and DNM3 were significantly altered in the endometrium of patients with endometriosis or PCOS compared to those of control subjects. Conclusion: Our study indicates the genetic correlation and shared risk genes between PCOS and endometriosis. These findings provide insights into the potential mechanisms behind their comorbidity and the future development of therapeutics.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech.

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization.

BACKGROUND: As genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype. RESULTS: In this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation . CONCLUSIONS: GaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

Investigating the shared genetic architecture between attention-deficit/hyperactivity disorder and risk taking behavior: A large-scale genomewide cross-trait analysis.

BACKGROUND: This study aims to explore the genetic architecture shared between Attention-Deficit/Hyperactivity Disorder (ADHD) and risk behavior. METHODS: Based on the latest large-scale Genome-wide association studies (GWAS), we firstly employed Linkage disequilibrium score regression (LDSC) and Local Analysis of Variant Association (LAVA) to investigate the genetic correlation between risk behavior and ADHD. Then, we conducted cross-trait analysis to identified the Pleiotropic loci. Finally, bidirectional Mendelian randomization analysis (MR) was applied to examine the causal relationship. RESULTS: We found a significant positive genetic correlation between ADHD and risk-taking behavior (rg = 0.351, p = 6.50E-37). The cross-trait meta-analysis identified 27 significant SNPs shared between ADHD and risk behavior. The most significant locus, located near the CADM2 gene on chromosome 3, had been identified associated with this two trait (pADHD = 3.07E-05 and prisk-taking behavior = 2.47E-30). The same situation can also be observed near the FOXP2 gene on chromosome 7 (rs8180817, pmeta = 5.72E-21). We found CCDC171 gene and other genes played a significant role in ADHD and risk behavior in mRNA level. Bidirectional MR analysis found a causal relationship between them. LIMITATION: The majority of our data sources were of European origin, which may limit the generalizability of our findings to other ethnic populations. CONCLUSION: This article reveals in depth the shared genetic structure between ADHD and risk-taking behavior, finding a significant positive genetic correlation between ADHD and risk-taking behavior. Providing insights for the future treatment and management of these two traits.

Genetic correlation and causal associations between psychiatric disorders and lung cancer risk.

BACKGROUND: Patients with certain psychiatric disorders have increased lung cancer incidence. However, establishing a causal relationship through traditional epidemiological methods poses challenges. METHODS: Available summary statistics of genome-wide association studies of cigarette smoking, lung cancer, and eight psychiatric disorders, including attention deficit/hyperactivity disorder (ADHD), autism, depression, major depressive disorder, bipolar disorder, insomnia, neuroticism, and schizophrenia (range N: 46,350-1,331,010) were leveraged to estimate genetic correlations using Linkage Disequilibrium Score Regression and assess causal effect of each psychiatric disorder on lung cancer using two-sample Mendelian randomization (MR) models, comprising inverse-variance weighted (IVW), weighted median, MR-Egger, pleiotropy residual sum and outlier testing (MR-PRESSO), and a constrained maximum likelihood approach (cML-MR). RESULTS: Significant positive correlations were observed between each psychiatric disorder and both smoking and lung cancer (all FDR < 0.05), except for the correlation between autism and lung cancer. Both univariable and the cML-MA MR analyses demonstrated that liability to schizophrenia, depression, ADHD, or insomnia was associated with an increased risk of overall lung cancer. Genetic liability to insomnia was linked specifically to squamous cell carcinoma (SCC), while genetic liability to ADHD was associated with an elevated risk of both SCC and small cell lung cancer (all P < 0.05). The later was further supported by multivariable MR analyses, which accounted for smoking. LIMITATIONS: Participants were constrained to European ancestry populations. Causal estimates from binary psychiatric disorders may be biased. CONCLUSION: Our findings suggest appropriate management of several psychiatric disorders, particularly ADHD, may potentially reduce the risk of developing lung cancer.

MPH: fast REML for large-scale genome partitioning of quantitative genetic variation.

MOTIVATION: Genome partitioning of quantitative genetic variation is useful for dissecting the genetic architecture of complex traits. However, existing methods, such as Haseman-Elston regression and linkage disequilibrium score regression, often face limitations when handling extensive farm animal datasets, as demonstrated in this study. RESULTS: To overcome this challenge, we present MPH, a novel software tool designed for efficient genome partitioning analyses using restricted maximum likelihood. The computational efficiency of MPH primarily stems from two key factors: the utilization of stochastic trace estimators and the comprehensive implementation of parallel computation. Evaluations with simulated and real datasets demonstrate that MPH achieves comparable accuracy and significantly enhances convergence, speed, and memory efficiency compared to widely used tools like GCTA and LDAK. These advancements facilitate large-scale, comprehensive analyses of complex genetic architectures in farm animals. AVAILABILITY AND IMPLEMENTATION: The MPH software is available at https://jiang18.github.io/mph/.

Autoimmune diseases and their genetic link to bronchiectasis: insights from a genetic correlation and Mendelian randomization study.

Background: Previous studies have demonstrated that autoimmune diseases are closely associated with bronchiectasis (BE). However, the causal effects between autoimmune diseases and BE remain elusive. Methods: All summary-level data were obtained from large-scale Genome-Wide Association Studies (GWAS). The univariate Mendelian randomization (UVMR) was utilized to investigate the genetic causal correlation (rg) of 12 autoimmune diseases and bronchiectasis, The Multivariable Mendelian Randomization (MVMR) method was used to explore the effects of the confounding factors. Further investigation was conducted to identify potential intermediate factors using mediation analysis. Finally, the linkage disequilibrium score regression (LDSC) method was used to identify genetic correlations among complex traits. A series of sensitivity analyses was performed to validate the robustness of the results. Results: The LDSC analysis revealed significant genetic correlations between BE and Crohn's disease (CD) (rg = 0.220, P = 0.037), rheumatoid arthritis (RA) (rg = 0.210, P = 0.021), and ulcerative colitis (UC) (rg = 0.247, P = 0.023). However, no genetic correlation was found with other autoimmune diseases (P > 0.05). The results of the primary IVW analysis suggested that for every SD increase in RA, there was a 10.3% increase in the incidence of BE (odds ratio [OR] = 1.103, 95% confidence interval [CI] 1.055-1.154, P = 1.75×10-5, FDR = 5.25×10-5). Furthermore, for every standard deviation (SD) increase in celiac disease (CeD), the incidence of BE reduced by 5.1% (OR = 0.949, 95% CI 0.902-0.999, P = 0.044, FDR = 0.044). We also observed suggestive evidence corresponding to a 3% increase in BE incidence with T1DM (OR = 1.033, 95% CI 1.001-1.066, P = 0.042, FDR = 0.063). Furthermore, MVMR analysis showed that RA was an independent risk factor for BE, whereas mediator MR analysis did not identify any mediating factors. The sensitivity analyses corroborated the robustness of these findings. Conclusion: LDSC analysis revealed significant genetic correlations between several autoimmune diseases and BE, and further MVMR analysis showed that RA is an independent risk factor for BE.

Accelerated biological aging as a potential risk factor for rheumatoid arthritis.

OBJECTS: Previous studies have suggested a potential correlation between rheumatoid arthritis (RA) and biological aging, but the intricate connections and mechanisms remain elusive. METHODS: In our study, we focused on two specific measures of biological age (PhenoAge and BioAge), which are derived from clinical biomarkers. The residuals of these measures, when compared to chronological age, are defined as biological age accelerations (BAAs). Utilizing the extensive UK Biobank dataset along with various genetic datasets, we conducted a thorough assessment of the relationship between BAAs and RA at both the individual and aggregate levels. RESULTS: Our observational studies revealed positive correlations between the two BAAs and the risk of developing both RA and seropositive RA. Furthermore, the genetic risk score (GRS) for PhenoAgeAccel was associated with an increased risk of RA and seropositive RA. Linkage disequilibrium score regression (LDSC) analysis further supported these findings, revealing a positive genetic correlation between PhenoAgeAccel and RA. PLACO analysis identified 38 lead pleiotropic single nucleotide polymorphisms linked to 301 genes, providing valuable insights into the potential mechanisms connecting PhenoAgeAccel and RA. CONCLUSION: In summary, our study has successfully revealed a positive correlation between accelerated biological aging, as measured by BAAs, and the susceptibility to RA.

Genetic insights into agronomic and morphological traits of drug-type cannabis revealed by genome-wide association studies.

Cannabis sativa L., previously concealed by prohibition, is now a versatile and promising plant, thanks to recent legalization, opening doors for medical research and industry growth. However, years of prohibition have left the Cannabis research community lagging behind in understanding Cannabis genetics and trait inheritance compared to other major crops. To address this gap, we conducted a comprehensive genome-wide association study (GWAS) of nine key agronomic and morphological traits, using a panel of 176 drug-type Cannabis accessions from the Canadian legal market. Utilizing high-density genotyping-by-sequencing (HD-GBS), we successfully generated dense genotyping data in Cannabis, resulting in a catalog of 800 K genetic variants, of which 282 K common variants were retained for GWAS analysis. Through GWAS analysis, we identified 18 markers significantly associated with agronomic and morphological traits. Several identified markers exert a substantial phenotypic impact, guided us to putative candidate genes that reside in high linkage-disequilibrium (LD) with the markers. These findings lay a solid foundation for an innovative cannabis research, leveraging genetic markers to inform breeding programs aimed at meeting diverse needs in the industry.

Genetic causal relationship between immune diseases and migraine: a Mendelian randomization study.

Background: Migraine has an increased prevalence in several immune disorders, but genetic cause-effect relationships remain unclear. Mendelian randomization (MR) was used in this study to explore whether immune diseases are causally associated with migraine and its subtypes. Methods: We conducted a two-sample bidirectional multivariate Mendelian randomization study. Single-nucleotide polymorphisms (SNP) for six immune diseases, including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), type 1 diabetes mellitus (T1D), allergic rhinitis (AR), asthma and psoriasis, were used as genetic instrumental variables. Summary statistics for migraine were obtained from 3 databases: the International Headache Genetics Consortium (IHGC), UK Biobank, and FinnGen study. MR analyses were performed per outcome database for each exposure and subsequently meta-analyzed. Reverse MR analysis was performed to determine whether migraine were risk factors for immune diseases. In addition, we conducted a genetic correlation to identify shared genetic variants for these two associations. Results: No significant causal relationship was found between immune diseases and migraine and its subtypes. These results were robust with a series of sensitivity analyses. Using the linkage disequilibrium score regression method (LDSC), we detected no genetic correlation between migraine and immune diseases. Conclusion: The evidence from our study does not support a causal relationship between immune diseases and migraine. The mechanisms underlying the frequent comorbidity of migraine and several immune diseases need to be further elucidated.

GAUSS: a summary-statistics-based R package for accurate estimation of linkage disequilibrium for variants, Gaussian imputation, and TWAS analysis of cosmopolitan cohorts.

MOTIVATION: As the availability of larger and more ethnically diverse reference panels grows, there is an increase in demand for ancestry-informed imputation of genome-wide association studies (GWAS), and other downstream analyses, e.g. fine-mapping. Performing such analyses at the genotype level is computationally challenging and necessitates, at best, a laborious process to access individual-level genotype and phenotype data. Summary-statistics-based tools, not requiring individual-level data, provide an efficient alternative that streamlines computational requirements and promotes open science by simplifying the re-analysis and downstream analysis of existing GWAS summary data. However, existing tools perform only disparate parts of needed analysis, have only command-line interfaces, and are difficult to extend/link by applied researchers. RESULTS: To address these challenges, we present Genome Analysis Using Summary Statistics (GAUSS)-a comprehensive and user-friendly R package designed to facilitate the re-analysis/downstream analysis of GWAS summary statistics. GAUSS offers an integrated toolkit for a range of functionalities, including (i) estimating ancestry proportion of study cohorts, (ii) calculating ancestry-informed linkage disequilibrium, (iii) imputing summary statistics of unobserved variants, (iv) conducting transcriptome-wide association studies, and (v) correcting for "Winner's Curse" biases. Notably, GAUSS utilizes an expansive, multi-ethnic reference panel consisting of 32 953 genomes from 29 ethnic groups. This panel enhances the range and accuracy of imputable variants, including the ability to impute summary statistics of rarer variants. As a result, GAUSS elevates the quality and applicability of existing GWAS analyses without requiring access to subject-level genotypic and phenotypic information. AVAILABILITY AND IMPLEMENTATION: The GAUSS R package, complete with its source code, is readily accessible to the public via our GitHub repository at https://github.com/statsleelab/gauss. To further assist users, we provided illustrative use-case scenarios that are conveniently found at https://statsleelab.github.io/gauss/, along with a comprehensive user guide detailed in Supplementary Text S1.

Genetic insight into a polygenic trait using a novel genome-wide association approach in a wild amphibian population.

Body size variation is central in the evolution of life-history traits in amphibians, but the underlying genetic architecture of this complex trait is still largely unknown. Herein, we studied the genetic basis of body size and fecundity of the alternative morphotypes in a wild population of the Greek smooth newt (Lissotriton graecus). By combining a genome-wide association approach with linkage disequilibrium network analysis, we were able to identify clusters of highly correlated loci thus maximizing sequence data for downstream analysis. The putatively associated variants explained 12.8% to 44.5% of the total phenotypic variation in body size and were mapped to genes with functional roles in the regulation of gene expression and cell cycle processes. Our study is the first to provide insights into the genetic basis of complex traits in newts and provides a useful tool to identify loci potentially involved in fitness-related traits in small data sets from natural populations in non-model species.

Bayesian approach to assessing population differences in genetic risk of disease with application to prostate cancer.

Population differences in risk of disease are common, but the potential genetic basis for these differences is not well understood. A standard approach is to compare genetic risk across populations by testing for mean differences in polygenic scores, but existing studies that use this approach do not account for statistical noise in effect estimates (i.e., the GWAS betas) that arise due to the finite sample size of GWAS training data. Here, we show using Bayesian polygenic score methods that the level of uncertainty in estimates of genetic risk differences across populations is highly dependent on the GWAS training sample size, the polygenicity (number of causal variants), and genetic distance (FST) between the populations considered. We derive a Wald test for formally assessing the difference in genetic risk across populations, which we show to have calibrated type 1 error rates under a simplified assumption that all SNPs are independent, which we achieve in practise using linkage disequilibrium (LD) pruning. We further provide closed-form expressions for assessing the uncertainty in estimates of relative genetic risk across populations under the special case of an infinitesimal genetic architecture. We suggest that for many complex traits and diseases, particularly those with more polygenic architectures, current GWAS sample sizes are insufficient to detect moderate differences in genetic risk across populations, though more substantial differences in relative genetic risk (relative risk > 1.5) can be detected. We show that conventional approaches that do not account for sampling error from the training sample, such as using a simple t-test, have very high type 1 error rates. When applying our approach to prostate cancer, we demonstrate a higher genetic risk in African Ancestry men, with lower risk in men of European followed by East Asian ancestry.

Poor statistical power in population-based association study of gene interaction.

BACKGROUND: Statistical epistasis, or "gene-gene interaction" in genetic association studies, means the nonadditive effects between the polymorphic sites on two different genes affecting the same phenotype. In the genetic association analysis of complex traits, nevertheless, the researchers haven't found enough clues of statistical epistasis so far. METHODS: We developed a statistical model where the statistical epistasis was presented as an extra linkage disequilibrium between the polymorphic sites of different risk genes. The power of statistical test for identifying the gene-gene interaction was calculated and then compared in different hypothesis scenarios. RESULTS: Our results show the statistical power increases with the increasing of interaction coefficient, relative risk, and linkage disequilibrium with genetic markers. However, the power of interaction discovery is much lower than that of regular single-site association test. When rigorous criteria were employed in statistical tests, the identification of gene-gene interaction became a very difficult task. Since the criterion of significance was given to be p-value ≤ 5.0 × 10-8, the same as that of many genome-wide association studies, there is little chance to identify the gene-gene interaction in all kind of circumstances. CONCLUSIONS: The lack of epistasis tends to be an inevitable result caused by the statistical principles of methods in the genetic association studies and therefore is the inherent characteristic of the research itself.

GWAS determined genetic loci associated with callus induction in oil palm tissue culture.

KEY MESSAGE: GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.