Addressing population-specific multiple testing burdens in genetic association studies

Genetic analyses unravel the causal association of cytokine levels on lichen simplex chronicus risk: insights from a mendelian randomization study

Lichen simplex chronicus (LSC) presents a challenge in dermatology due to its elusive pathogenic mechanisms. While associations between circulating inflammatory cytokines and LSC were observed, the definitive causal dynamics remain to be elucidated. Our study used a two-sample Mendelian randomization (MR) approach to investigate causal relationships. We applied a suite of MR methodologies, including IVW, Weighted Median, MR-Egger, Weighted Mode, Simple Mode, MR-PRESSO, and the Steiger test, to ensure robust causal inference. Our analysis confirmed the causal impact of genetically determined cytokine levels on LSC risk, particularly MMP-10 (OR = 0.493, P = 0.004) and DNER (OR = 0.651, P = 0.043) in risk attenuation. We also found a positive causal correlation between GDNF levels (OR = 1.871, P = 0.007) and LSC prevalence. Notably, bidirectional causality was observed between DNER and LSC. Consistency across various MR analyses and sensitivity analyses confirmed the absence of horizontal pleiotropy, validating the causal estimates. This pioneering MR investigation unveils a novel genetically anchored causal relationship between the circulating levels of MMP-10, DNER, and GDNF and LSC risk. Although further validation is requisite, our findings augment the understanding of cytokine mediation in LSC and underscore prospective avenues for research.

A Multibreed Genome-Wide Association Study for Cattle Leukocyte Telomere Length

Telomeres are terminal DNA regions of chromosomes that prevent chromosomal fusion and degradation during cell division. In cattle, leukocyte telomere length (LTL) is associated with longevity, productive lifespan, and disease susceptibility. However, the genetic basis of LTL in this species is less studied than in humans. In this study, we utilized the whole-genome resequencing data of 239 animals from 17 cattle breeds for computational leukocyte telomere length estimation and subsequent genome-wide association study of LTL. As a result, we identified 42 significant SNPs, of which eight were found in seven genes (EXOC6B, PTPRD, RPS6KC1, NSL1, AGBL1, ENSBTAG00000052188, and GPC1) when using covariates for two major breed groups (Turano-Mongolian and European). Association analysis with covariates for breed effect detected 63 SNPs, including 13 in five genes (EXOC6B, PTPRD, RPS6KC1, ENSBTAG00000040318, and NELL1). The PTPRD gene, demonstrating the top signal in analysis with breed effect, was previously associated with leukocyte telomere length in cattle and likely is involved in the mechanism of alternative lengthening of telomeres. The single nucleotide variants found could be tested for marker-assisted selection to improve telomere-length-associated traits.

Genetic Dissection of Spike Productivity Traits in the Siberian Collection of Spring Barley

Barley (Hordeum vulgare L.) is one of the most commonly cultivated cereals worldwide. Its local varieties can represent a valuable source of unique genetic variants useful for crop improvement. The aim of this study was to reveal loci contributing to spike productivity traits in Siberian spring barley and to develop diagnostic DNA markers for marker-assisted breeding programs. For this purpose we conducted a genome-wide association study using a panel of 94 barley varieties. In total, 64 SNPs significantly associated with productivity traits were revealed. Twenty-three SNP markers were validated by genotyping in an independent sample set using competitive allele-specific PCR (KASP). Finally, fourteen markers associated with spike productivity traits on chromosomes 2H, 4H and 5H can be suggested for use in breeding programs.

Germline Cancer Gene Expression Quantitative Trait Loci Are Associated with Local and Global Tumor Mutations

Somatic mutations drive cancer development and are relevant to patient responses to treatment. Emerging evidence shows that variations in the somatic genome can be influenced by the germline genetic background. However, the mechanisms underlying these germline-somatic associations remain largely obscure. We hypothesized that germline variants can influence somatic mutations in a nearby cancer gene ("local impact") or a set of recurrently mutated cancer genes across the genome ("global impact") through their regulatory effect on gene expression. To test this hypothesis, tumor targeted sequencing data from 12,413 patients across 11 cancer types in the Dana-Farber Profile cohort were integrated with germline cancer gene expression quantitative trait loci (eQTL) from the Genotype-Tissue Expression Project. Variants that upregulate ATM expression were associated with a decreased risk of somatic ATM mutations across 8 cancer types. GLI2, WRN, and CBFB eQTL were associated with global tumor mutational burden of cancer genes in ovarian cancer, glioma, and esophagogastric carcinoma, respectively. An EPHA5 eQTL was associated with mutations in cancer genes specific to colorectal cancer, and eQTL related to expression of APC, WRN, GLI1, FANCA, and TP53 were associated with mutations in genes specific to endometrial cancer. These findings provide evidence that germline-somatic associations are mediated through expression of specific cancer genes, opening new avenues for research on the underlying biological processes.

SIGNIFICANCE

Analysis of associations between the germline genetic background and somatic mutations in patients with cancer suggests that germline variants can influence local and global tumor mutations by altering expression of cancer-related genes. See related commentary by Kar, p. 1165.

Late-onset neonatal sepsis: genetic differences by sex and involvement of the NOTCH pathway

BACKGROUND

Late-Onset Neonatal Sepsis (LOS) is a rare condition, involving widespread infection, immune disruption, organ dysfunction, and often death. Because exposure to pathogens is not completely preventable, identifying susceptibility factors is critical to characterizing the pathophysiology and developing interventions. Prior studies demonstrated both genetics and infant sex influence susceptibility. Our study was designed to identify LOS associated genetic variants.

METHODS

We performed an exploratory genome wide association study (GWAS) with 224 LOS cases and 273 controls from six European countries. LOS was defined as sepsis presenting from 3 to 90 days of age; diagnosis was established by clinical criteria consensus guidelines. We tested for association with both autosomal and X-chromosome variants in the total sample and in sex-stratified analyses.

RESULTS

In total, 71 SNPs associated with neonatal sepsis at p < 1 × 10^-4 in at least one analysis. Most importantly, sex-stratified analyses revealed associations with multiple SNPs (28 in males and 16 in females), but no variants from single-sex analyses associated with sepsis in the other sex. Pathway analyses showed NOTCH signaling is over-represented among genes linked to these SNPS.

CONCLUSION

Our results indicate genetic susceptibility to LOS is sexually dimorphic and corroborate that NOTCH signaling plays a role in determining risk.

IMPACT

Genes associate with late onset neonatal sepsis. Notch pathway genes are overrepresented in associations with sepsis. Genes associating with sepsis do not overlap between males and females. Sexual dimorphism can lead to sex specific treatment of sepsis.

Novel genetic loci associated with osteoarthritis in multi-ancestry analyses in the Million Veteran Program and UK Biobank

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain.

Genome-wide study on 72,298 individuals in Korean biobank data for 76 traits

Genome-wide association studies (GWAS) on diverse ancestry groups are lacking, resulting in deficits of genetic discoveries and polygenic scores. We conducted GWAS for 76 phenotypes in Korean biobank data, namely the Korean Genome and Epidemiology Study (KoGES) (n = 72,298). Our analysis discovered 2,242 associated loci, including 122 novel associations, many of which were replicated in Biobank Japan (BBJ) GWAS. We also applied several up-to-date methods for genetic association tests to increase the power, discovering additional associations that are not identified in simple case-control GWAS. We evaluated genetic pleiotropy to investigate genes associated with multiple traits. Following meta-analysis of 32 phenotypes between KoGES and BBJ, we further identified 379 novel associations and demonstrated the improved predictive performance of polygenic risk scores by using the meta-analysis results. The summary statistics of 76 KoGES GWAS phenotypes are publicly available, contributing to a better comprehension of the genetic architecture of the East Asian population.

Genetic and environmental drivers of large-scale epigenetic variation in Thlaspi arvense

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

Causal Effects of Circulating Cytokines on the Risk of Psoriasis Vulgaris: A Mendelian Randomization Study

Background: Psoriasis vulgaris is an inflammatory skin disease. Observational studies have shown associations between circulating cytokine levels and psoriasis vulgaris. But the causal relationship between circulating cytokine and psoriasis vulgaris remains elusive.

Methods: To assess the causal effects of cytokine levels on the risk of psoriasis vulgaris and vice versa, we performed a two-sample Mendelian randomization (MR) study by using the inverse-variance weighted (IVW), weighted median, and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) in genome-wide association summary statistics of 41 circulating cytokines in up to 8,293 individuals and psoriasis vulgaris in 399,883 individuals.

Results: We identified that increasing RANTES level induced an elevated risk of psoriasis vulgaris in IVW (β = 0.33, S.E. = 0.12, p = 0.006). This causal effect showed consistency across the weighted median (β = 0.35, S.E. = 0.15, p = 0.022) and MR-PRESSO method (β = 0.33, S.E. = 0.11, p = 0.028).

Conclusions: Our results suggest a potential causal effect of elevated RANTES concentration on the increased risk of psoriasis vulgaris.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Leveraging deep phenotyping from health check-up cohort with 10,000 Korean individuals for phenome-wide association study of 136 traits

The expanding use of the phenome-wide association study (PheWAS) faces challenges in the context of using International Classification of Diseases billing codes for phenotype definition, imbalanced study population ethnicity, and constrained application of the results in research. We performed a PheWAS utilizing 136 deep phenotypes corroborated by comprehensive health check-ups in a Korean population, along with trans-ethnic comparisons through using the UK Biobank and Biobank Japan Project. Meta-analysis with Korean and Japanese population was done. The PheWAS associated 65 phenotypes with 14,101 significant variants (P < 4.92 × 10-10). Network analysis, visualization of cross-phenotype mapping, and causal inference mapping with Mendelian randomization were conducted. Among phenotype pairs from the genotype-driven cross-phenotype associations, we evaluated penetrance in correlation analysis using a clinical database. We focused on the application of PheWAS in order to make it robust and to aid the derivation of biological meaning post-PheWAS. This comprehensive analysis of PheWAS results based on a health check-up database will provide researchers and clinicians with a panoramic overview of the networks among multiple phenotypes and genetic variants, laying groundwork for the practical application of precision medicine.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma

BACKGROUND

Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent.

METHODS

We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 μg in children, 100 μg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV1. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 μg twice daily in children and 500 μg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 μg twice daily in children, 250 μg twice daily in adolescents and adults), and 5 × ICS versus 100 μg fluticasone plus a LABA (salmeterol 50 μg twice daily). We used a genome-wide significance threshold of p<1·6 × 10-4, and tested for replication using independent cohorts of individuals of African descent with asthma.

FINDINGS

We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 μg fluticasone plus salmeterol on chromosome 12 (odds ratio [ORlocal African] 3·95, 95% CI 2·02-7·72, p=6·1 × 10-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, ORallele dose 0·17, 95% CI 0·07-0·42, p=8·4 × 10-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (ORlocal African 3·35, 1·98-5·67, p=6·8 × 10-6) containing a locus adjacent to TPST2 (rs5752429, ORallele dose 0·21, 0·09-0·52, p=5·7 × 10-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts.

INTERPRETATION

BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma.

FUNDING

National Institutes of Health, National Heart, Lung, and Blood Institute.

SNP Data Science for Classification of Bipolar Disorder I and Bipolar Disorder II

Bipolar disorder I (BD-I) and bipolar disorder II (BD-II) have specific characteristics and clear diagnostic criteria, but quite different treatment guidelines. In clinical practice, BD-II is commonly mistaken as a mild form of BD-I. This study uses data science technique to identify the important Single Nucleotide Polymorphisms (SNPs) significantly affecting the classifications of BD-I and BD-II, and develops a set of complementary diagnostic classifiers to enhance the diagnostic process. Screening assessments and SNP genotypes of 316 Han Chinese were performed with the Affymetrix Axiom Genome-Wide TWB Array Plate. The results show that the classifier constructed by 23 SNPs reached the area under curve of ROC (AUC) level of 0.939, while the classifier constructed by 42 SNPs reached the AUC level of 0.9574, which is a mere addition of 1.84 percent. The accuracy rate of classification increased by 3.46 percent. This study also uses Gene Ontology (GO) and Pathway to conduct a functional analysis and identify significant items, including calcium ion binding, GABA-A receptor activity, Rap1 signaling pathway, ECM proteoglycans, IL12-mediated signaling events, Nicotine addiction), and PI3K-Akt signaling pathway. The study can address time-consuming SNPs identification and also quantify the effect of SNP-SNP interactions.

GWAS identifies an ortholog of the rice D11 gene as a candidate gene for grain size in an international collection of hexaploid wheat

Grain size is a key agronomic trait that contributes to grain yield in hexaploid wheat. Grain length and width were evaluated in an international collection of 157 wheat accessions. These accessions were genetically characterized using a genotyping-by-sequencing (GBS) protocol that produced 73,784 single nucleotide polymorphism (SNP) markers. GBS-derived genotype calls obtained on Chinese Spring proved extremely accurate when compared to the reference (> 99.9%) and showed > 95% agreement with calls made at SNP loci shared with the 90 K SNP array on a subset of 71 Canadian wheat accessions for which both types of data were available. This indicates that GBS can yield a large amount of highly accurate SNP data in hexaploid wheat. The genetic diversity analysis performed using this set of SNP markers revealed the presence of six distinct groups within this collection. A GWAS was conducted to uncover genomic regions controlling variation for grain length and width. In total, seven SNPs were found to be associated with one or both traits, identifying three quantitative trait loci (QTLs) located on chromosomes 1D, 2D and 4A. In the vicinity of the peak SNP on chromosome 2D, we found a promising candidate gene (TraesCS2D01G331100), whose rice ortholog (D11) had previously been reported to be involved in the regulation of grain size. These markers will be useful in breeding for enhanced wheat productivity.

A high-resolution genome-wide association study of the grain ionome and agronomic traits in rice Oryza sativa subsp. indica

This study presents a comprehensive study of the genetic bases controlling variation in the rice ionome employing genome-wide association studies (GWAS) with a diverse panel of indica accessions, each genotyped with 5.2 million markers. GWAS was performed for twelve elements including B, Ca, Co, Cu, Fe, K, Mg, Mn, Mo, Na, P, and Zn and four agronomic traits including days to 50% flowering, grain yield, plant height and thousand grain weight. GWAS identified 128 loci associated with the grain elements and 57 associated with the agronomic traits. There were sixteen co-localization regions containing QTL for two or more traits. Fourteen grain element quantitative trait loci were stable across growing environments, which can be strong candidates to be used in marker-assisted selection to improve the concentrations of nutritive elements in rice grain. Potential candidate genes were revealed including OsNAS3 linked to the locus that controls the variation of Zn and Co concentrations. The effects of starch synthesis and grain filling on multiple grain elements were elucidated through the likely involvement of OsSUS1 and OsGSSB1 genes. Overall, our study provides crucial insights into the genetic basis of ionomic variations in rice and will facilitate improvement in breeding for trace mineral content.

A Multiancestry Sex-Stratified Genome-Wide Association Study of Spontaneous Clearance of Hepatitis C Virus

BACKGROUND

Spontaneous clearance of acute hepatitis C virus (HCV) infection is more common in women than in men, independent of known risk factors.

METHODS

To identify sex-specific genetic loci, we studied 4423 HCV-infected individuals (2903 male, 1520 female) of European, African, and Hispanic ancestry. We performed autosomal, and X chromosome sex-stratified and combined association analyses in each ancestry group.

RESULTS

A male-specific region near the adenosine diphosphate-ribosylation factor-like 5B (ARL5B) gene was identified. Individuals with the C allele of rs76398191 were about 30% more likely to have chronic HCV infection than individuals with the T allele (OR, 0.69; P = 1.98 × 10-07), and this was not seen in females. The ARL5B gene encodes an interferon-stimulated gene that inhibits immune response to double-stranded RNA viruses. We also identified suggestive associations near septin 6 and ribosomal protein L39 genes on the X chromosome. In box sexes, allele G of rs12852885 was associated with a 40% increase in HCV clearance compared with the A allele (OR, 1.4; P = 2.46 × 10-06). Septin 6 facilitates HCV replication via interaction with the HCV NS5b protein, and ribosomal protein L39 acts as an HCV core interactor.

CONCLUSIONS

These novel gene associations support differential mechanisms of HCV clearance between the sexes and provide biological targets for treatment or vaccine development.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

A deoxyribonuclease 1-like 3 genetic variant associates with asthma exacerbations

The analysis of genetic variants from six genes whose expression is predictive of asthma exacerbations revealed a novel association of a regulatory polymorphism from DNASE1L3 in African-Americans and Latinos.

Recent shifts in the genomic ancestry of Mexican Americans may alter the genetic architecture of biomedical traits

People in the Americas represent a diverse continuum of populations with varying degrees of admixture among African, European, and Amerindigenous ancestries. In the United States, populations with non-European ancestry remain understudied, and thus little is known about the genetic architecture of phenotypic variation in these populations. Using genotype data from the Hispanic Community Health Study/Study of Latinos, we find that Amerindigenous ancestry increased by an average of ~20% spanning 1940s-1990s in Mexican Americans. These patterns result from complex interactions between several population and cultural factors which shaped patterns of genetic variation and influenced the genetic architecture of complex traits in Mexican Americans. We show for height how polygenic risk scores based on summary statistics from a European-based genome-wide association study perform poorly in Mexican Americans. Our findings reveal temporal changes in population structure within Hispanics/Latinos that may influence biomedical traits, demonstrating a need to improve our understanding of admixed populations.

The influence of anthropogenic habitat fragmentation on the genetic structure and diversity of the malaria vector Anopheles cruzii (Diptera: Culicidae)

Fragmentation of natural environments as a result of human interference has been associated with a decrease in species richness and increase in abundance of a few species that have adapted to these environments. The Brazilian Atlantic Forest, which has been undergoing an intense process of fragmentation and deforestation caused by human-made changes to the environment, is an important hotspot for malaria transmission. The main vector of simian and human malaria in this biome is the mosquito Anopheles cruzii. Anthropogenic processes reduce the availability of natural resources at the tree canopies, An. cruzii primary habitat. As a consequence, An. cruzii moves to the border of the Atlantic Forest nearing urban areas seeking resources, increasing their contact with humans in the process. We hypothesized that different levels of anthropogenic changes to the environment can be an important factor in driving the genetic structure and diversity in An. cruzii populations. Five different hypotheses using a cross-sectional and a longitudinal design were tested to assess genetic structure in sympatric An. cruzii populations and microevolutionary processes driving these populations. Single nucleotide polymorphisms were used to assess microgeographic genetic structure in An. cruzii populations in a low-endemicity area in the city of São Paulo, Brazil. Our results show an overall weak genetic structure among the populations, indicating a high gene flow system. However, our results also pointed to the presence of significant genetic structure between sympatric An. cruzii populations collected at ground and tree-canopy habitats in the urban environment and higher genetic variation in the ground-level population. This indicates that anthropogenic modifications leading to habitat fragmentation and a higher genetic diversity and structure in ground-level populations could be driving the behavior of An. cruzii, ultimately increasing its contact with humans. Understanding how anthropogenic changes in natural areas affect An. cruzii is essential for the development of more effective mosquito control strategies and, on a broader scale, for malaria-elimination efforts in the Brazilian Atlantic Forest.

Identification of pleiotropic genes between risk factors of stroke by multivariate metaCCA analysis

Genome-wide association studies (GWASs) have identified more than 20 genetic loci as risk predictors associated with stroke. However, these studies were generally performed for single-trait and failed to consider the pleiotropic effects of these risk genes among the multiple risk factors for stroke. In this study, we applied a novel metaCCA method followed by gene-based VEGAS2 analysis to identify the risk genes for stroke that may overlap between seven correlated risk factors (including atrial fibrillation, hypertension, coronary artery disease, heart failure, diabetes, body mass index, and total cholesterol level) by integrating seven corresponding GWAS data. We detected 20 potential pleiotropic genes that may be associated with multiple risk factors of stroke. Furthermore, using gene-to-trait pathway analysis, we suggested six potential risk genes (FUT8, GMIP, PLA2G6, PDE3A, SMARCA4, SKAPT) that may affect ischemic or hemorrhage stroke through multiple intermediate factors such as MAPK family. These findings provide novel insight into the genetic determinants contributing to the concurrent development of biological conditions that may influence stroke susceptibility, and also indicate some potential therapeutic targets that can be further studied for the prevention of cerebrovascular disease.

Genome-wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward-related ventral striatum activity in African- and European-Americans

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome-wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European-Americans (EA; 2927 cases) and 3132 African-Americans (AA: 1315 cases) participating in the family-based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h 2 in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome-wide significant (GWS; P < 5E-08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion-deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans-ancestral meta-analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward-related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non-European samples with distinct patterns of substance use may lead to the identification of novel ancestry-specific genetic markers of risk.

Multi-Ancestry Genome-Wide Association Study of Spontaneous Clearance of Hepatitis C Virus

BACKGROUND & AIMS

Spontaneous clearance of hepatitis C virus (HCV) occurs in approximately 30% of infected persons and less often in populations of African ancestry. Variants in major histocompatibility complex (MHC) and in interferon lambda genes are associated with spontaneous HCV clearance, but there have been few studies of these variants in persons of African ancestry. We performed a dense multi-ancestry genome-wide association study of spontaneous clearance of HCV, focusing on individuals of African ancestry.

METHODS

We performed genotype analyses of 4423 people from 3 ancestry groups: 2201 persons of African ancestry (445 with HCV clearance and 1756 with HCV persistence), 1739 persons of European ancestry (701 with HCV clearance and 1036 with HCV persistence), and 486 multi-ancestry Hispanic persons (173 with HCV clearance and 313 with HCV persistence). Samples were genotyped using Illumina (San Diego, CA) arrays and statistically imputed to the 1000 Genomes Project. For each ancestry group, the association of single-nucleotide polymorphisms with HCV clearance was tested by log-additive analysis, and then a meta-analysis was performed.

RESULTS

In the meta-analysis, significant associations with HCV clearance were confirmed at the interferon lambda gene locus IFNL4-IFNL3 (19q13.2) (P = 5.99 × 10-50) and the MHC locus 6p21.32 (P = 1.15 × 10-21). We also associated HCV clearance with polymorphisms in the G-protein-coupled receptor 158 gene (GPR158) at 10p12.1 (P = 1.80 × 10-07). These 3 loci had independent, additive effects of HCV clearance, and account for 6.8% and 5.9% of the variance of HCV clearance in persons of European and African ancestry, respectively. Persons of African or European ancestry carrying all 6 variants were 24-fold and 11-fold, respectively, more likely to clear HCV infection compared with individuals carrying none or 1 of the clearance-associated variants.

CONCLUSIONS

In a meta-analysis of data from 3 studies, we found variants in MHC genes, IFNL4-IFNL3, and GPR158 to increase odds of HCV clearance in patients of European and African ancestry. These findings could increase our understanding of immune response to and clearance of HCV infection.

A simple and accurate method to determine genomewide significance for association tests in sequencing studies

Whole-exome sequencing (WES) and whole-genome sequencing (WGS) studies are underway to investigate the impact of genetic variants on complex diseases and traits. It is customary to perform single-variant association tests for common variants and region-based association tests for rare variants. The latter may target variants with similar or opposite effects, interrogate variants with different frequencies or different functional annotations, and examine a variety of regions. The large number of tests that are performed necessitates adjustment for multiple testing. The conventional Bonferroni correction is overly conservative as the test statistics are correlated. To address this challenge, we propose a simple and accurate method based on parametric bootstrap to assess genomewide significance. We show that the correlations of the test statistics are determined primarily by the genotypes, such that the same significance threshold can be used in different studies that share a common sequencing platform. We demonstrate the usefulness of the proposed method with WES data from the National Heart, Lung, and Blood Institute Exome Sequencing Project and WGS data from the 1000 Genomes Project. We recommend the p value of 5 × 1 0 - 9 as the genomewide significance threshold for testing all common and low-frequency variants (MAFs ≥ 0.1%) in the human genome.

Variation in DNA methylation of human blood over a 1-year period using the Illumina MethylationEPIC array

Assessing DNA methylation profiles in human blood has become a major focus of epidemiologic inquiry. Understanding variability in CpG-specific DNA methylation over moderate periods of time is a critical first step in identifying CpG sites that are candidates for DNA methylation-based etiologic, diagnostic and prognostic predictors of pathogenesis. Using the Illumina MethylationEPIC [850K] BeadArray, DNA methylation was profiled in paired whole blood samples collected approximately 1 year apart from 35 healthy women enrolled in the Nurses Study II cohort. The median intraclass correlation coefficient (ICC) across all CpG loci was 0.19 [Interquartile Range (IQR) 0.00-0.50]; 74.8% of ICCs were in the low range (0-0.5), 16.9% in the mid-range of ICCs (0.5-0.8), and 8.3% in the high-range of ICCs (0.8-1). ICCs were similar for CpG probes on the 450K Illumina array (median 0.17) and the new probes added to the 850K array (median 0.21). ICCs for CpG loci on the sex chromosomes and known metastable epialleles were high (median 0.71, 0.97, respectively), and ICCs among methylation quantitative trait loci (mQTL) CpGs were significantly higher as compared to non-mQTL CpGs (median 0.73, 0.16, respectively, P < 2 × 10-16). We observed wide variation in DNA methylation stability over a 1-year period. Probes considered non-stable, due to substantial variation over a moderate period of time and with minimal variability across individuals could be removed in large epidemiological studies. Moreover, adjusting for technical variation that arises from using high-dimensional arrays is critical.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma

Short-acting β₂-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the U.S^{1, 2}. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma³. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1,830 Latinos (Total=2,779). Lastly, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p=7.69×10⁻⁹). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958, and rs7081864, p≤5×10⁻⁸). Our results failed to replicate in three additional populations of 416 Latinos and 1,615 African Americans. Our findings indicate that both population specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Genetic Determinants of Telomere Length in African American Youth

Telomere length (TL) is associated with numerous disease states and is affected by genetic and environmental factors. However, TL has been mostly studied in adult populations of European or Asian ancestry. These studies have identified 34 TL-associated genetic variants recently used as genetic proxies for TL. The generalizability of these associations to pediatric populations and racially diverse populations, specifically of African ancestry, remains unclear. Furthermore, six novel variants associated with TL in a population of European children have been identified but not validated. We measured TL from whole blood samples of 492 healthy African American youth (children and adolescents between 8 and 20 years old) and performed the first genome-wide association study of TL in this population. We were unable to replicate neither the 34 reported genetic associations found in adults nor the six genetic associations found in European children. However, we discovered a novel genome-wide significant association between TL and rs1483898 on chromosome 14. Our results underscore the importance of examining genetic associations with TL in diverse pediatric populations such as African Americans.

Fine-mapping analysis of a chromosome 2 region linked to resistance to Mycobacterium tuberculosis infection in Uganda reveals potential regulatory variants

Tuberculosis (TB) is a major public health burden worldwide, and more effective treatment is sorely needed. Consequently, uncovering causes of resistance to Mycobacterium tuberculosis (Mtb) infection is of special importance for vaccine design. Resistance to Mtb infection can be defined by a persistently negative tuberculin skin test (PTST-) despite living in close and sustained exposure to an active TB case. While susceptibility to Mtb is, in part, genetically determined, relatively little work has been done to uncover genetic factors underlying resistance to Mtb infection. We examined a region on chromosome 2q previously implicated in our genomewide linkage scan by a targeted, high-density association scan for genetic variants enhancing PTST- in two independent Ugandan TB household cohorts (n = 747 and 471). We found association with SNPs in neighboring genes ZEB2 and GTDC1 (peak meta p = 1.9 × 10^-5) supported by both samples. Bioinformatic analysis suggests these variants may affect PTST- by regulating the histone deacetylase (HDAC) pathway, supporting previous results from transcriptomic analyses. An apparent protective effect of PTST- against body-mass wasting suggests a link between resistance to Mtb infection and healthy body composition. Our results provide insight into how humans may escape latent Mtb infection despite heavy exposure.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma

RATIONALE

Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response.

OBJECTIVES

To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children.

METHODS

We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR.

MEASUREMENTS AND MAIN RESULTS

We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings.

CONCLUSIONS

The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

Identification of novel loci associated with maturity and yield traits in early maturity soybean plant introduction lines

BACKGROUND

To continue to meet the increasing demands of soybean worldwide, it is crucial to identify key genes regulating flowering and maturity to expand the cultivated regions into short season areas. Although four soybean genes have been successfully utilized in early maturity breeding programs, new genes governing maturity are continuously being identified suggesting that there remains as yet undiscovered loci governing agronomic traits of interest. The objective of this study was to identify novel loci and genes involved in a diverse set of early soybean maturity using genome-wide association (GWA) analyses to identify loci governing days to maturity (DTM), flowering (DTF) and pod filling (DTPF), as well as yield and 100 seed weight in Canadian environments. To do so, soybean plant introduction lines varying significantly for maturity, but classified as early varieties, were used. Plants were phenotyped for the five agronomic traits for five site-years and GWA approaches used to identify candidate loci and genes affecting each trait.

RESULTS

Genotyping using genotyping-by-sequencing and microarray methods identified 67,594 single nucleotide polymorphisms, of which 31,283 had a linkage disequilibrium < 1 and minor allele frequency > 0.05 and were used for GWA analyses. A total of 9, 6, 4, 5 and 2 loci were detected for GWA analyses for DTM, DTF, DTPF, 100 seed weight and yield, respectively. Regions of interest, including a region surrounding the E1 gene for flowering and maturity, and several novel loci, were identified, with several loci having pleiotropic effects. Novel loci affecting maturity were identified on chromosomes five and 13 and reduced maturity by 7.2 and 3.3 days, respectively. Novel loci for maturity and flowering contained genes orthologous to known Arabidopsis flowering genes, while loci affecting yield and 100 seed weight contained genes known to cause dwarfism.

CONCLUSIONS

This study demonstrated substantial variation in soybean agronomic traits of interest, including maturity and flowering dates as well as yield, and the utility of GWA analyses in identifying novel genetic factors underlying important agronomic traits. The loci and candidate genes identified serve as promising targets for future studies examining the mechanisms underlying the related soybean traits.

Sphingolipid Metabolic Pathway Impacts Thiazide Diuretics Blood Pressure Response: Insights From Genomics, Metabolomics, and Lipidomics

BACKGROUND

Although hydrochlorothiazide (HCTZ) is a well-established first-line antihypertensive in the United States, <50% of HCTZ treated patients achieve blood pressure (BP) control. Thus, identifying biomarkers that could predict the BP response to HCTZ is critically important. In this study, we utilized metabolomics, genomics, and lipidomics to identify novel pathways and biomarkers associated with HCTZ BP response.

METHODS AND RESULTS

First, we conducted a pathway analysis for 13 metabolites we recently identified to be significantly associated with HCTZ BP response. From this analysis, we found the sphingolipid metabolic pathway as the most significant pathway (P=5.8E-05). Testing 78 variants, within 14 genes involved in the sphingolipid metabolic canonical pathway, with the BP response to HCTZ identified variant rs6078905, within the SPTLC3 gene, as a novel biomarker significantly associated with the BP response to HCTZ in whites (n=228). We found that rs6078905 C-allele carriers had a better BP response to HCTZ versus noncarriers (∆SBP/∆DBP: -11.4/-6.9 versus -6.8/-3.5 mm Hg; ∆SBP P=6.7E-04; ∆DBP P=4.8E-04). Additionally, in blacks (n=148), we found genetic signals in the SPTLC3 genomic region significantly associated with the BP response to HCTZ (P<0.05). Last, we observed that rs6078905 significantly affects the baseline level of 4 sphingomyelins (N24:2, N24:3, N16:1, and N22:1; false discovery rate <0.05), from which N24:2 sphingomyelin has a significant correlation with both HCTZ DBP-response (r=-0.42; P=7E-03) and SBP-response (r=-0.36; P=2E-02).

CONCLUSIONS

This study provides insight into potential pharmacometabolomic and genetic mechanisms underlying HCTZ BP response and suggests that SPTLC3 is a potential determinant of the BP response to HCTZ.

CLINICAL TRIAL REGISTRATION

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00246519.

Current Scope and Challenges in Phenome-Wide Association Studies

PURPOSE OF REVIEW

Over many decades, researchers have been designing studies to investigate the relationship between genotypes and phenotypes to gain an understanding about the effect of genetics on disease. Recently, a high-throughput approach called phenome-wide associations studies (PheWAS) have been extensively used to identify associations between genetic variants and many diseases and traits simultaneously. In this review, we describe the value of PheWAS along with methodological issues and challenges in interpretation for current applications of PheWAS.

RECENT FINDINGS

PheWAS have uncovered a paradigm to identify new associations for genetic loci across many diseases. The application of PheWAS have been effective with phenotype data from electronic health records, epidemiological studies, and clinical trials data.

SUMMARY

The key strength of a PheWAS is to identify the association of one or more genetic variants with multiple phenotypes, which can showcase interconnections among the phenotypes due to shared genetic associations. While the PheWAS approach appears promising, there are a number of challenges that need to be addressed to provide additional robustness to PheWAS findings.

An eQTL variant of ZXDC is associated with IFN-γ production following Mycobacterium tuberculosis antigen-specific stimulation

There is a large inter-individual variability in the response to Mycobacterium tuberculosis infection. In previous linkage analyses, we identified a major locus on chromosome region 8q controlling IFN-γ production after stimulation with live BCG (Bacillus Calmette-Guérin), and a second locus on chromosome region 3q affecting IFN-γ production triggered by the 6-kDa early secretory antigen target (ESAT-6), taking into account the IFN-γ production induced by BCG (IFNγ-ESAT6_BCG). High-density genotyping and imputation identified ~100,000 variants within each linkage region, which we tested for association with the corresponding IFN-γ phenotype in families from a tuberculosis household contact study in France. Significant associations were replicated in a South African familial sample. The most convincing association observed was that between the IFNγ-ESAT6_BCG phenotype and rs9828868 on chromosome 3q (p = 9.8 × 10⁻⁶ in the French sample). This variant made a significant contribution to the linkage signal (p < 0.001), and a trend towards the same association was observed in the South African sample. This variant was reported to be an eQTL of the ZXDC gene, biologically linked to monocyte IL-12 production through CCL2/MCP1. The identification of rs9828868 as a genetic driver of IFNγ production in response to mycobacterial antigens provides new insights into human anti-tuberculosis immunity.

A genome-wide quantitative trait locus (QTL) linkage scan of NEO personality factors in Latino families segregating bipolar disorder

Personality traits have been suggested as potential endophenotypes for Bipolar Disorder (BP), as they can be quantitatively measured and show correlations with BP. The present study utilized data from 2,745 individuals from 686 extended pedigrees originally ascertained for having multiplex cases of BP (963 cases of BPI or schizoaffective BP). Subjects were assessed with the NEO Personality Inventory, Revised (NEO PI-R) and genotyped using the Illumina HumanLinkage-24 Bead Chip, with an average genetic coverage of 0.67 cM. Two point linkage scores were calculated for each trait as a quantitative variable using SOLAR (Sequential Oligogenic Linkage Analysis Routines). Suggestive evidence for linkage was found for neuroticism at 1q32.1 (LOD = 2.52), 6q23.3 (2.32), 16p12 (2.79), extraversion at 4p15.3 (2.33), agreeableness at 4q31.1 (2.37), 5q34 (2.80), 7q31.1 (2.56), 16q22 (2.52), and conscientiousness at 4q31.1 (2.50). Each of the above traits have been shown to be correlated with the broad BP phenotype in this same sample. In addition, for the trait of openness, we found significant evidence of linkage to chromosome 3p24.3 (rs336610, LOD = 4.75) and suggestive evidence at 1q43 (2.74), 5q35.1 (3.03), 11q14.3 (2.61), 11q21 (2.30), and 19q13.1 (2.52). These findings support previous linkage findings of the openness trait to chromosome 19q13 and the agreeableness trait to 4q31 and identify a number of new loci for personality endophenotypes related to bipolar disorder.

Methods to analyze big data in pharmacogenomics research

The scale and scope of pharmacogenomics research continues to expand as the cost and efficiency of molecular data generation techniques advance. These new technologies give rise to enormous opportunity for the identification of important genetic and genomic factors important for drug treatment response. With this opportunity come significant challenges. Most of these can be categorized as 'big data' issues, facing not only pharmacogenomics, but other fields in the life sciences as well. In this review, we describe some of the analysis techniques and tools being implemented for genetic/genomic discovery in pharmacogenomics.

Gene Set Enrichment Analyses: lessons learned from the heart failure phenotype

BACKGROUND

Genetic studies for complex diseases have predominantly discovered main effects at individual loci, but have not focused on genomic and environmental contexts important for a phenotype. Gene Set Enrichment Analysis (GSEA) aims to address this by identifying sets of genes or biological pathways contributing to a phenotype, through gene-gene interactions or other mechanisms, which are not the focus of conventional association methods.

RESULTS

Approaches that utilize GSEA can now take input from array chips, either gene-centric or genome-wide, but are highly sensitive to study design, SNP selection and pruning strategies, SNP-to-gene mapping, and pathway definitions. Here, we present lessons learned from our experience with GSEA of heart failure, a particularly challenging phenotype due to its underlying heterogeneous etiology.

CONCLUSIONS

This case study shows that proper data handling is essential to avoid false-positive results. Well-defined pipelines for quality control are needed to avoid reporting spurious results using GSEA.

Genetic Effects on the Correlation Structure of CVD Risk Factors: Exome-Wide Data From a Ghanaian Population

Plasma concentration of plasminogen activator inhibitor-1 (PAI-1) is highly correlated with several cardiovascular disease (CVD) risk factors. It also plays a direct role in CVD, including myocardial infarction and stroke, by impeding the dissolution of thrombi in the blood. Insofar as PAI-1 links CVD's risk factors to its endpoints, genetic variants modulating the relationship between PAI-1 and risk factors may be of particular clinical and biological interest. The high heritability of PAI-1, which has not been explained by genetic association studies, may also, in large part, be due to this relationship with CVD risk factors. Using exome-wide data from 1,032 Ghanaian study participants, we tested for heterogeneity of correlation by genotype between PAI-1 and 4 CVD risk factors (body mass index, triglycerides, mean arterial pressure, and fasting glucose) under the hypothesis that loci involved in the relationship between PAI-1 and other risk factors will also modify their correlational structure. We found more significant heterogeneities of correlation by genotype than we found marginal effects, with no evidence of type I inflation. The most significant result among all univariate and multivariate tests performed in this study was the heterogeneity of correlation between PAI-1 and mean arterial pressure at rs10738554, near SLC24A2, a gene previously associated with high blood pressure in African Americans.

NFAT5 and SLC4A10 Loci Associate with Plasma Osmolality

Disorders of water balance, an excess or deficit of total body water relative to body electrolyte content, are common and ascertained by plasma hypo- or hypernatremia, respectively. We performed a two-stage genome-wide association study meta-analysis on plasma sodium concentration in 45,889 individuals of European descent (stage 1 discovery) and 17,637 additional individuals of European descent (stage 2 replication), and a transethnic meta-analysis of replicated single-nucleotide polymorphisms in 79,506 individuals (63,526 individuals of European descent, 8765 individuals of Asian Indian descent, and 7215 individuals of African descent). In stage 1, we identified eight loci associated with plasma sodium concentration at P<5.0 × 10^-6 Of these, rs9980 at NFAT5 replicated in stage 2 meta-analysis (P=3.1 × 10^-5), with combined stages 1 and 2 genome-wide significance of P=5.6 × 10^-10 Transethnic meta-analysis further supported the association at rs9980 (P=5.9 × 10^-12). Additionally, rs16846053 at SLC4A10 showed nominally, but not genome-wide, significant association in combined stages 1 and 2 meta-analysis (P=6.7 × 10^-8). NFAT5 encodes a ubiquitously expressed transcription factor that coordinates the intracellular response to hypertonic stress but was not previously implicated in the regulation of systemic water balance. SLC4A10 encodes a sodium bicarbonate transporter with a brain-restricted expression pattern, and variant rs16846053 affects a putative intronic NFAT5 DNA binding motif. The lead variants for NFAT5 and SLC4A10 are cis expression quantitative trait loci in tissues of the central nervous system and relevant to transcriptional regulation. Thus, genetic variation in NFAT5 and SLC4A10 expression and function in the central nervous system may affect the regulation of systemic water balance.

Genetic architecture differences between pediatric and adult-onset inflammatory bowel diseases in the Polish population

Most inflammatory bowel diseases (IBDs) are classic complex disorders represented by common alleles. Here we aimed to define the genetic architecture of pediatric and adult-onset IBDs for the Polish population. A total of 1495 patients were recruited, including 761 patients with Crohn’s disease (CD; 424 pediatric), 734 patients with ulcerative colitis (UC; 390 pediatric), and 934 healthy controls. Allelotyping employed a pooled-DNA genome-wide association study (GWAS) and was validated by individual genotyping. Whole exome sequencing (WES) was performed on 44 IBD patients diagnosed before 6 years of age, 45 patients diagnosed after 40 years of age, and 18 healthy controls. Altogether, out of 88 selected SNPs, 31 SNPs were replicated for association with IBD. A novel BRD2 (rs1049526) association reached significance of P = 5.2 × 10⁻¹¹ and odds ratio (OR) = 2.43. Twenty SNPs were shared between pediatric and adult patients; 1 and 7 were unique to adult-onset and pediatric-onset IBD, respectively. WES identified numerous rare and potentially deleterious variants in IBD-associated or innate immunity-associated genes. Deleterious alleles in both groups were over-represented among rare variants in affected children. Our GWAS revealed differences in the polygenic architecture of pediatric- and adult-onset IBD. A significant accumulation of rare and deleterious variants in affected children suggests a contribution by yet unexplained genetic components.

Resetting the bar: Statistical significance in whole-genome sequencing-based association studies of global populations

Genome-wide association studies (GWAS) of common disease have been hugely successful in implicating loci that modify disease risk. The bulk of these associations have proven robust and reproducible, in part due to community adoption of statistical criteria for claiming significant genotype-phenotype associations. As the cost of sequencing continues to drop, assembling large samples in global populations is becoming increasingly feasible. Sequencing studies interrogate not only common variants, as was true for genotyping-based GWAS, but variation across the full allele frequency spectrum, yielding many more (independent) statistical tests. We sought to empirically determine genome-wide significance thresholds for various analysis scenarios. Using whole-genome sequence data, we simulated sequencing-based disease studies of varying sample size and ancestry. We determined that future sequencing efforts in >2,000 samples of European, Asian, or admixed ancestry should set genome-wide significance at approximately P = 5 × 10^-9 , and studies of African samples should apply a more stringent genome-wide significance threshold of P = 1 × 10^-9 . Adoption of a revised multiple test correction will be crucial in avoiding irreproducible claims of association.

Phenome-Wide Association Study to Explore Relationships between Immune System Related Genetic Loci and Complex Traits and Diseases

We performed a Phenome-Wide Association Study (PheWAS) to identify interrelationships between the immune system genetic architecture and a wide array of phenotypes from two de-identified electronic health record (EHR) biorepositories. We selected variants within genes encoding critical factors in the immune system and variants with known associations with autoimmunity. To define case/control status for EHR diagnoses, we used International Classification of Diseases, Ninth Revision (ICD-9) diagnosis codes from 3,024 Geisinger Clinic MyCode^® subjects (470 diagnoses) and 2,899 Vanderbilt University Medical Center BioVU biorepository subjects (380 diagnoses). A pooled-analysis was also carried out for the replicating results of the two data sets. We identified new associations with potential biological relevance including SNPs in tumor necrosis factor (TNF) and ankyrin-related genes associated with acute and chronic sinusitis and acute respiratory tract infection. The two most significant associations identified were for the C6orf10 SNP rs6910071 and “rheumatoid arthritis” (ICD-9 code category 714) (p_METAL = 2.58 x 10⁻⁹) and the ATN1 SNP rs2239167 and “diabetes mellitus, type 2” (ICD-9 code category 250) (p_METAL = 6.39 x 10⁻⁹). This study highlights the utility of using PheWAS in conjunction with EHRs to discover new genotypic-phenotypic associations for immune-system related genetic loci.

Empirical estimation of genome-wide significance thresholds based on the 1000 Genomes Project data set

To assess the statistical significance of associations between variants and traits, genome-wide association studies (GWAS) should employ an appropriate threshold that accounts for the massive burden of multiple testing in the study. Although most studies in the current literature commonly set a genome-wide significance threshold at the level of P=5.0 × 10^-8, the adequacy of this value for respective populations has not been fully investigated. To empirically estimate thresholds for different ancestral populations, we conducted GWAS simulations using the 1000 Genomes Phase 3 data set for Africans (AFR), Europeans (EUR), Admixed Americans (AMR), East Asians (EAS) and South Asians (SAS). The estimated empirical genome-wide significance thresholds were P_sig=3.24 × 10^-8 (AFR), 9.26 × 10^-8 (EUR), 1.83 × 10^-7 (AMR), 1.61 × 10^-7 (EAS) and 9.46 × 10^-8 (SAS). We additionally conducted trans-ethnic meta-analyses across all populations (ALL) and all populations except for AFR (ΔAFR), which yielded P_sig=3.25 × 10^-8 (ALL) and 4.20 × 10^-8 (ΔAFR). Our results indicate that the current threshold (P=5.0 × 10^-8) is overly stringent for all ancestral populations except for Africans; however, we should employ a more stringent threshold when conducting a meta-analysis, regardless of the presence of African samples.

Novel genetic risk factors for asthma in African American children: Precision Medicine and the SAGE II Study

Asthma, an inflammatory disorder of the airways, is the most common chronic disease of children worldwide. There are significant racial/ethnic disparities in asthma prevalence, morbidity, and mortality among US children. This trend is mirrored in obesity, which may share genetic and environmental risk factors with asthma. The majority of asthma biomedical research has been performed in populations of European decent. We sought to identify genetic risk factors for asthma in African American children. We also assessed the generalizability of genetic variants associated with asthma in European and Asian populations to African American children. Our study population consisted of 1227 (812 asthma cases, 415 controls) African American children with genome-wide single nucleotide polymorphism (SNP) data. Logistic regression was used to identify associations between SNP genotype and asthma status. We identified a novel variant in the PTCHD3 gene that is significantly associated with asthma (rs660498, p = 2.2 × 10(-7)) independent of obesity status. Approximately 5 % of previously reported asthma genetic associations identified in European populations replicated in African Americans. Our identification of novel variants associated with asthma in African American children, coupled with our inability to replicate the majority of findings reported in European Americans, underscores the necessity for including diverse populations in biomedical studies of asthma.

A Locus at 5q33.3 Confers Resistance to Tuberculosis in Highly Susceptible Individuals

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

INTEGRATING CLINICAL LABORATORY MEASURES AND ICD-9 CODE DIAGNOSES IN PHENOME-WIDE ASSOCIATION STUDIES

Electronic health records (EHR) provide a comprehensive resource for discovery, allowing unprecedented exploration of the impact of genetic architecture on health and disease. The data of EHRs also allow for exploration of the complex interactions between health measures across health and disease. The discoveries arising from EHR based research provide important information for the identification of genetic variation for clinical decision-making. Due to the breadth of information collected within the EHR, a challenge for discovery using EHR based data is the development of high-throughput tools that expose important areas of further research, from genetic variants to phenotypes. Phenome-Wide Association studies (PheWAS) provide a way to explore the association between genetic variants and comprehensive phenotypic measurements, generating new hypotheses and also exposing the complex relationships between genetic architecture and outcomes, including pleiotropy. EHR based PheWAS have mainly evaluated associations with case/control status from International Classification of Disease, Ninth Edition (ICD-9) codes. While these studies have highlighted discovery through PheWAS, the rich resource of clinical lab measures collected within the EHR can be better utilized for high-throughput PheWAS analyses and discovery. To better use these resources and enrich PheWAS association results we have developed a sound methodology for extracting a wide range of clinical lab measures from EHR data. We have extracted a first set of 21 clinical lab measures from the de-identified EHR of participants of the Geisinger MyCodeTM biorepository, and calculated the median of these lab measures for 12,039 subjects. Next we evaluated the association between these 21 clinical lab median values and 635,525 genetic variants, performing a genome-wide association study (GWAS) for each of 21 clinical lab measures. We then calculated the association between SNPs from these GWAS passing our Bonferroni defined p-value cutoff and 165 ICD-9 codes. Through the GWAS we found a series of results replicating known associations, and also some potentially novel associations with less studied clinical lab measures. We found the majority of the PheWAS ICD-9 diagnoses highly related to the clinical lab measures associated with same SNPs. Moving forward, we will be evaluating further phenotypes and expanding the methodology for successful extraction of clinical lab measurements for research and PheWAS use. These developments are important for expanding the PheWAS approach for improved EHR based discovery.

Phenome-Wide Association Studies: Leveraging Comprehensive Phenotypic and Genotypic Data for Discovery

With the large volume of clinical and epidemiological data being collected, increasingly linked to extensive genotypic data, coupled with expanding high-performance computational resources, there are considerable opportunities for comprehensively exploring the networks of connections that exist between the phenome and the genome. These networks can be identified through Phenome-Wide Association Studies (PheWAS) where the association between a collection of genetic variants, or in some cases a particular clinical lab variable, and a wide and diverse range of phenotypes, diagnoses, traits, and/or outcomes are evaluated. This is a departure from the more familiar genome-wide association study (GWAS) approach, which has been used to identify single nucleotide polymorphisms (SNPs) associated with one outcome or a very limited phenotypic domain. In addition to highlighting novel connections between multiple phenotypes and elucidating more of the phenotype-genotype landscape, PheWAS can generate new hypotheses for further exploration, and can also be used to narrow the search space for research using comprehensive data collections. The complex results of PheWAS also have the potential for uncovering new mechanistic insights. We review here how the PheWAS approach has been used with data from epidemiological studies, clinical trials, and de-identified electronic health record data. We also review methodologies for the analyses underlying PheWAS, and emerging methods developed for evaluating the comprehensive results of PheWAS including genotype-phenotype networks. This review also highlights PheWAS as an important tool for identifying new biomarkers, elucidating the genetic architecture of complex traits, and uncovering pleiotropy. There are many directions and new methodologies for the future of PheWAS analyses, from the phenotypic data to the genetic data, and herein we also discuss some of these important future PheWAS developments.