A fast and powerful linear mixed model approach for genotype-environment interaction tests in large-scale GWAS.

Genotype-by-environment interaction (GEI or GxE) plays an important role in understanding complex human traits. However, it is usually challenging to detect GEI signals efficiently and accurately while adjusting for population stratification and sample relatedness in large-scale genome-wide association studies (GWAS). Here we propose a fast and powerful linear mixed model-based approach, fastGWA-GE, to test for GEI effect and G + GxE joint effect. Our extensive simulations show that fastGWA-GE outperforms other existing GEI test methods by controlling genomic inflation better, providing larger power and running hundreds to thousands of times faster. We performed a fastGWA-GE analysis of ~7.27 million variants on 452 249 individuals of European ancestry for 13 quantitative traits and five environment variables in the UK Biobank GWAS data and identified 96 significant signals (72 variants across 57 loci) with GEI test P-values < 1 × 10-9, including 27 novel GEI associations, which highlights the effectiveness of fastGWA-GE in GEI signal discovery in large-scale GWAS.

Generalized multi-SNP mediation intersection-union test.

To elucidate the molecular mechanisms underlying genetic variants identified from genome-wide association studies (GWAS) for a variety of phenotypic traits encompassing binary, continuous, count, and survival outcomes, we propose a novel and flexible method to test for mediation that can simultaneously accommodate multiple genetic variants and different types of outcome variables. Specifically, we employ the intersection-union test approach combined with the likelihood ratio test to detect mediation effect of multiple genetic variants via some mediator (e.g., the expression of a neighboring gene) on outcome. We fit high-dimensional generalized linear mixed models under the mediation framework, separately under the null and alternative hypothesis. We leverage Laplace approximation to compute the marginal likelihood of outcome and use coordinate descent algorithm to estimate corresponding parameters. Our extensive simulations demonstrate the validity of our proposed methods and substantial, up to 97%, power gains over alternative methods. Applications to real data for the study of Chlamydia trachomatis infection further showcase advantages of our methods. We believe our proposed methods will be of value and general interest in this post-GWAS era to disentangle the potential causal mechanism from DNA to phenotype for new drug discovery and personalized medicine.

Host Genetic Risk Factors for Chlamydia trachomatis-Related Infertility in Women.

BACKGROUND: Chlamydia trachomatis (Ct) infection ascending to the upper genital tract can cause infertility. Direct association of genetic variants as contributors is challenging because infertility may not be diagnosed until years after infection. Investigating the intermediate trait of ascension bridges this gap. METHODS: We identified infertility genome-wide association study (GWAS) loci using deoxyribonucleic acid from Ct-seropositive cisgender women in a tubal factor infertility study and Ct-infected cisgender women from a longitudinal pelvic inflammatory disease cohort with known fertility status. Deoxyribonucleic acid and blood messenger ribonucleic acid from 2 additional female cohorts with active Ct infection and known endometrial infection status were used to investigate the impact of infertility single-nucleotide polymorphisms (SNPs) on Ct ascension. A statistical mediation test examined whether multiple infertility SNPs jointly influenced ascension risk by modulating expression of mediator genes. RESULTS: We identified 112 candidate infertility GWAS loci, and 31 associated with Ct ascension. The SNPs altered chlamydial ascension by modulating expression of 40 mediator genes. Mediator genes identified are involved in innate immune responses including type I interferon production, T-cell function, fibrosis, female reproductive tract health, and protein synthesis and degradation. CONCLUSIONS: We identified Ct-related infertility loci and their potential functional effects on Ct ascension.

Multi-SNP mediation intersection-union test.

SUMMARY: Tens of thousands of reproducibly identified GWAS (Genome-Wide Association Studies) variants, with the vast majority falling in non-coding regions resulting in no eventual protein products, call urgently for mechanistic interpretations. Although numerous methods exist, there are few, if any methods, for simultaneously testing the mediation effects of multiple correlated SNPs via some mediator (e.g. the expression of a gene in the neighborhood) on phenotypic outcome. We propose multi-SNP mediation intersection-union test (SMUT) to fill in this methodological gap. Our extensive simulations demonstrate the validity of SMUT as well as substantial, up to 92%, power gains over alternative methods. In addition, SMUT confirmed known mediators in a real dataset of Finns for plasma adiponectin level, which were missed by many alternative methods. We believe SMUT will become a useful tool to generate mechanistic hypotheses underlying GWAS variants, facilitating functional follow-up. AVAILABILITY AND IMPLEMENTATION: The R package SMUT is publicly available from CRAN at https://CRAN.R-project.org/package=SMUT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discovery of Blood Transcriptional Endotypes in Women with Pelvic Inflammatory Disease.

Sexually transmitted infections with Chlamydia trachomatis and/or Neisseria gonorrhoeae and rates of pelvic inflammatory disease (PID) in women continue to rise, with reinfection being common because of poor adaptive immunity. Diagnosis remains imprecise, and pathogenesis data are derived primarily from monoinfection of mice with C. trachomatis or N. gonorrhoeae By comparing blood mRNA responses of women with C. trachomatis- and/or N. gonorrhoeae-induced PID and histologic endometritis with those from women with C. trachomatis and/or N. gonorrhoeae infection limited to their cervix and asymptomatic uninfected women determined via microarray, we discovered important pathogenic mechanisms in PID and response differences that provide a pathway to biomarker discovery. Women with N. gonorrhoeae- and/or C. trachomatis-induced PID exhibit overexpression of myeloid cell genes and suppression of protein synthesis, mitochondrial oxidative phosphorylation, and T cell-specific genes. Coinfected women exhibited the greatest activation of cell death pathways and suppression of responses essential for adaptive immunity. Women solely infected with C. trachomatis expressed elevated levels of type I and type II IFN genes, and enhanced type I IFN-induced chemokines in cervical secretions were associated with ascension of C. trachomatis to the endometrium. Blood microarrays reveal discrete pathobiological endotypes in women with PID that are driven by pathogen invasion of the upper genital tract.

Cervical Cytokines Associated With Chlamydia trachomatis Susceptibility and Protection.

BACKGROUND: Chlamydia trachomatis can cause reproductive morbidities after ascending to the upper genital tract of women, and repeated infection can lead to worse disease. Data related to protective immune responses at the cervical mucosa that could limit chlamydial infection to the cervix and/or prevent reinfection inform vaccine approaches and biomarkers of risk. METHODS: We measured 48 cytokines in cervical secretions from women having chlamydial cervical infection alone (n = 92) or both cervical and endometrial infection (n = 68). Univariable regression identified cytokines associated with differential odds of endometrial infection and reinfection risk, and multivariable stepwise regression identified cytokine ratios associated with differential risk. RESULTS: Elevated interleukin (IL) 15/CXCL10 (odds ratio [OR], 0.55 [95% confidence interval {CI}, .37-.78]), IL-16/tumor necrosis factor-α (OR, 0.66 [95% CI, .45-.93]), and CXCL14/IL-17A (OR, 0.73 [95% CI, .54-.97]) cytokine ratios were significantly (P ≤ .05) associated with decreased odds of endometrial infection. A higher Flt-3L/IL-14 ratio was significantly (P = .001) associated with a decreased risk of reinfection (hazard ratio, 0.71 [95% CI, .58-.88]). CONCLUSIONS: Cytokines involved in humoral, type I interferon, and T-helper (Th) 17 responses were associated with susceptibility to C. trachomatis, whereas cytokines involved in Th1 polarization, recruitment, and activation were associated with protection against ascension and reinfection.

Whole-Exome Sequencing to Identify Novel Biological Pathways Associated With Infertility After Pelvic Inflammatory Disease.

BACKGROUND: Ideal management of sexually transmitted infections (STI) may require risk markers for pathology or vaccine development. Previously, we identified common genetic variants associated with chlamydial pelvic inflammatory disease (PID) and reduced fecundity. As this explains only a proportion of the long-term morbidity risk, we used whole-exome sequencing to identify biological pathways that may be associated with STI-related infertility. METHODS: We obtained stored DNA from 43 non-Hispanic black women with PID from the PID Evaluation and Clinical Health Study. Infertility was assessed at a mean of 84 months. Principal component analysis revealed no population stratification. Potential covariates did not significantly differ between groups. Sequencing kernel association test was used to examine associations between aggregates of variants on a single gene and infertility. The results from the sequencing kernel association test were used to choose "focus genes" (P < 0.01; n = 150) for subsequent Ingenuity Pathway Analysis to identify "gene sets" that are enriched in biologically relevant pathways. RESULTS: Pathway analysis revealed that focus genes were enriched in canonical pathways including, IL-1 signaling, P2Y purinergic receptor signaling, and bone morphogenic protein signaling. CONCLUSIONS: Focus genes were enriched in pathways that impact innate and adaptive immunity, protein kinase A activity, cellular growth, and DNA repair. These may alter host resistance or immunopathology after infection. Targeted sequencing of biological pathways identified in this study may provide insight into STI-related infertility.

Analysis of Factors Driving Incident and Ascending Infection and the Role of Serum Antibody in Chlamydia trachomatis Genital Tract Infection.

BACKGROUND: Chlamydia trachomatis genital tract infection is a major cause of female reproductive morbidity. Risk factors for ascending infection are unknown, and the role for antibody in protection is not well established. METHODS: We recruited 225 women from urban outpatient clinics and followed them for a median of 12 months. We performed a cross-sectional analysis of serum anti-chlamydial immunoglobulin G (IgG), behavioral factors, and microbiological factors associated with endometrial infection at enrollment, and a longitudinal analysis of factors associated with incident infection. RESULTS: Oral contraceptives (adjusted relative risk [RR], 2.02 [95% confidence interval {CI}, 1.38-2.97]) and gonorrhea (adjusted RR, 1.66 [95% CI, 1.07-2.60]) were associated with endometrial infection. Gonorrhea (adjusted hazard ratio [HR], 3.09 [95% CI, 1.41-6.78]), cervical infection at enrollment (adjusted HR, 2.33 [95% CI, 1.07-5.11]), and exposure to uncircumcised partners (adjusted HR, 2.65 [95% CI, 1.21-5.82]) or infected partners (adjusted HR, 4.99 [95% CI, 2.66-9.39]) significantly increased the risk of incident infection. Seropositivity was associated with a reduced cervical burden (P < .05) but no differences in rates of ascending infection (adjusted RR, 1.24 [95% CI, .71-2.19]) or incident infection (adjusted HR, 0.94 [95% CI, .52-1.69]). CONCLUSIONS: Serum anti-chlamydial IgG is not associated with a lowered rate of ascending or repeat infection. Identification of factors associated with ascending infection and increased risk of incident infection provide guidance for targeted screening of women at increased risk for sequelae.

Identification of Chlamydia trachomatis Antigens Recognized by T Cells From Highly Exposed Women Who Limit or Resist Genital Tract Infection.

BACKGROUND: Natural infection induces partial immunity to Chlamydia trachomatis Identification of chlamydial antigens that induce interferon γ (IFN-) secretion by T cells from immune women could advance vaccine development. METHODS: IFN-γ production by CD4+ and CD8+ peripheral blood T cells from 58 high-risk women was measured after coculture with antigen-presenting cells preincubated with recombinant Escherichia coli expressing a panel of 275 chlamydial antigens. Quantile median regression analysis was used to compare frequencies of IFN-γ responses in women with only cervical infection to those in women with endometrial infection and frequencies in women who remained uninfected for over 1 year to those in women who developed incident infection. Statistical methods were then used to identify proteins associated with protection. RESULTS: A higher median frequency of CD8+ T-cell responses was detected in women with lower genital tract chlamydial infection, compared with those with upper genital tract chlamydial infection (13.8% vs 9.5%; P =04), but the CD4+ T-cell response frequencies were not different. Women who remained uninfected displayed a greater frequency of positive CD4+ T-cell responses (29% vs 18%; P < .0001), compared with women who had incident infection, while the frequencies of CD8+ T-cell responses did not differ. A subset of proteins involved in central metabolism, type III secretion, and protein synthesis were associated with protection. CONCLUSIONS: Investigations in naturally exposed women reveal protective responses and identify chlamydial vaccine candidate antigens.

Understanding Regulatory Mechanisms of Brain Function and Disease through 3D Genome Organization.

The human genome has a complex and dynamic three-dimensional (3D) organization, which plays a critical role for gene regulation and genome function. The importance of 3D genome organization in brain development and function has been well characterized in a region- and cell-type-specific fashion. Recent technological advances in chromosome conformation capture (3C)-based techniques, imaging approaches, and ligation-free methods, along with computational methods to analyze the data generated, have revealed 3D genome features at different scales in the brain that contribute to our understanding of genetic mechanisms underlying neuropsychiatric diseases and other brain-related traits. In this review, we discuss how these advances aid in the genetic dissection of brain-related traits.

Understanding the function of regulatory DNA interactions in the interpretation of non-coding GWAS variants.

Genome-wide association studies (GWAS) have identified a vast number of variants associated with various complex human diseases and traits. However, most of these GWAS variants reside in non-coding regions producing no proteins, making the interpretation of these variants a daunting challenge. Prior evidence indicates that a subset of non-coding variants detected within or near cis-regulatory elements (e.g., promoters, enhancers, silencers, and insulators) might play a key role in disease etiology by regulating gene expression. Advanced sequencing- and imaging-based technologies, together with powerful computational methods, enabling comprehensive characterization of regulatory DNA interactions, have substantially improved our understanding of the three-dimensional (3D) genome architecture. Recent literature witnesses plenty of examples where using chromosome conformation capture (3C)-based technologies successfully links non-coding variants to their target genes and prioritizes relevant tissues or cell types. These examples illustrate the critical capability of 3D genome organization in annotating non-coding GWAS variants. This review discusses how 3D genome organization information contributes to elucidating the potential roles of non-coding GWAS variants in disease etiology.

Genetic susceptibility loci for Chlamydia trachomatis endometrial infection influence expression of genes involved in T cell function, tryptophan metabolism and epithelial integrity.

Objectives: Identify genetic loci of enhanced susceptibility to Chlamydial trachomatis (Ct) upper genital tract infection in women. Methods: We performed an integrated analysis of DNA genotypes and blood-derived mRNA profiles from 200 Ct-exposed women to identify expression quantitative trait loci (eQTL) and determine their association with endometrial chlamydial infection using a mediation test. We further evaluated the effect of a lead eQTL on the expression of CD151 by immune cells from women with genotypes associated with low and high whole blood expression of CD151, respectively. Results: We identified cis-eQTLs modulating mRNA expression of 81 genes (eGenes) associated with altered risk of ascending infection. In women with endometrial infection, eGenes involved in proinflammatory signaling were upregulated. Downregulated eGenes included genes involved in T cell functions pivotal for chlamydial control. eGenes encoding molecules linked to metabolism of tryptophan, an essential chlamydial nutrient, and formation of epithelial tight junctions were also downregulated in women with endometrial infection. A lead eSNP rs10902226 was identified regulating CD151, a tetrospanin molecule important for immune cell adhesion and migration and T cell proliferation. Further in vitro experiments showed that women with a CC genotype at rs10902226 had reduced rates of endometrial infection with increased CD151 expression in whole blood and T cells when compared to women with a GG genotype. Conclusions: We discovered genetic variants associated with altered risk for Ct ascension. A lead eSNP for CD151 is a candidate genetic marker for enhanced CD4 T cell function and reduced susceptibility.

Inferring Regulatory Networks From Mixed Observational Data Using Directed Acyclic Graphs.

Construction of regulatory networks using cross-sectional expression profiling of genes is desired, but challenging. The Directed Acyclic Graph (DAG) provides a general framework to infer causal effects from observational data. However, most existing DAG methods assume that all nodes follow the same type of distribution, which prohibit a joint modeling of continuous gene expression and categorical variables. We present a new mixed DAG (mDAG) algorithm to infer the regulatory pathway from mixed observational data containing both continuous variables (e.g. expression of genes) and categorical variables (e.g. categorical phenotypes or single nucleotide polymorphisms). Our method can identify upstream causal factors and downstream effectors closely linked to a variable and generate hypotheses for causal direction of regulatory pathways. We propose a new permutation method to test the conditional independence of variables of mixed types, which is the key for mDAG. We also utilize an L 1 regularization in mDAG to ensure it can recover a large sparse DAG with limited sample size. We demonstrate through extensive simulations that mDAG outperforms two well-known methods in recovering the true underlying DAG. We apply mDAG to a cross-sectional immunological study of Chlamydia trachomatis infection and successfully infer the regularity network of cytokines. We also apply mDAG to a large cohort study, generating sensible mechanistic hypotheses underlying plasma adiponectin level. The R package mDAG is publicly available from CRAN at https://CRAN.R-project.org/package=mDAG.

Anti-chlamydia IgG and IgA are insufficient to prevent endometrial chlamydia infection in women, and increased anti-chlamydia IgG is associated with enhanced risk for incident infection.

PROBLEM: Chlamydia infections in women can ascend to the upper genital tract, and repeated infections are common, placing women at risk for sequelae. The protective role of anti-chlamydia antibodies to surface exposed antigens in ascending and incident infection is unclear. METHOD OF STUDY: A whole-bacterial ELISA was used to quantify chlamydia-specific IgG and IgA in serum and cervical secretions of 151 high-risk women followed longitudinally. Correlations were determined between antibody and cervical burden, and causal mediation analysis investigated the effect of antibody on ascension. We examined the relationship of antibody to incident infection using the marginal Cox model. RESULTS: Serum and cervical anti-chlamydia IgG and cervical IgA levels correlated inversely with cervical burden. While lower burden was associated with reduced ascension, causal mediation analysis revealed that the indirect effects of antibody mediated through reductions in bacterial burden were insufficient to prevent ascension. Analysis of women uninfected at enrollment revealed that serum and cervical anti-chlamydia IgG were associated with increased risk of incident infection; hazard ratio increased 3.6-fold (95% CI, 1.3-10.3), and 22.6-fold (95% CI, 3.1-165.2) with each unit of serum and cervical IgG, respectively. CONCLUSION: Although anti-chlamydia IgG and IgA correlated with reduced cervical chlamydia burden, they failed to prevent ascension and increased levels of anti-chlamydia IgG were associated with increased risk for incident infection.

Gene Expression Signatures Can Aid Diagnosis of Sexually Transmitted Infection-Induced Endometritis in Women.

Sexually transmitted infection (STI) of the upper reproductive tract can result in inflammation and infertility. A biomarker of STI-induced upper tract inflammation would be significant as many women are asymptomatic and delayed treatment increases risk of sequelae. Blood mRNA from 111 women from three cohorts was profiled using microarray. Unsupervised analysis revealed a transcriptional profile that distinguished 9 cases of STI-induced endometritis from 18 with cervical STI or uninfected controls. Using a hybrid feature selection algorithm we identified 21 genes that yielded maximal classification accuracy within our training dataset. Predictive accuracy was evaluated using an independent testing dataset of 5 cases and 10 controls. Sensitivity was evaluated in a separate test set of 12 women with asymptomatic STI-induced endometritis in whom cervical burden was determined by PCR; and specificity in an additional test set of 15 uninfected women with pelvic pain due to unknown cause. Disease module preservation was assessed in 42 women with a clinical diagnosis of pelvic inflammatory disease (PID). We also tested the ability of the biomarker to discriminate STI-induced endometritis from other diseases. The biomarker was 86.7% (13/15) accurate in correctly distinguishing cases from controls in the testing dataset. Sensitivity was 83.3% (5/6) in women with high cervical Chlamydia trachomatis burden and asymptomatic endometritis, but 0% (0/6) in women with low burden. Specificity in patients with non-STI-induced pelvic pain was 86.7% (13/15). Disease modules were preserved in all 8 biomarker predicted cases. The 21-gene biomarker was highly discriminatory for systemic infections, lupus, and appendicitis, but wrongly predicted tuberculosis as STI-induced endometritis in 52.4%. A 21-gene biomarker can identify asymptomatic women with STI-induced endometritis that places them at risk for chronic disease development and discriminate STI-induced endometritis from non-STI pelvic pain and other diseases.