Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10-20), ER-negative BC (P=1.1 × 10-13), BRCA1-associated BC (P=7.7 × 10-16) and triple negative BC (P-diff=2 × 10-5). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10-3) and ABHD8 (P<2 × 10-3). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3′-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk

A Germline Variant at 8q24 Contributes to Familial Clustering of Prostate Cancer in Men of African Ancestry

Although men of African ancestry have a high risk of prostate cancer (PCa), no genes or mutations have been identified that contribute to familial clustering of PCa in this population. We investigated whether the African ancestry–specific PCa risk variant at 8q24, rs72725854, is enriched in men with a PCa family history in 9052 cases, 143 cases from high-risk families, and 8595 controls of African ancestry. We found the risk allele to be significantly associated with earlier age at diagnosis, more aggressive disease, and enriched in men with a PCa family history (32% of high-risk familial cases carried the variant vs 23% of cases without a family history and 12% of controls). For cases with two or more first-degree relatives with PCa who had at least one family member diagnosed at age <60 yr, the odds ratios for TA heterozygotes and TT homozygotes were 3.92 (95% confidence interval [CI] = 2.13–7.22) and 33.41 (95% CI = 10.86–102.84), respectively. Among men with a PCa family history, the absolute risk by age 60 yr reached 21% (95% CI = 17–25%) for TA heterozygotes and 38% (95% CI = 13–65%) for TT homozygotes. We estimate that in men of African ancestry, rs72725854 accounts for 32% of the total familial risk explained by all known PCa risk variants. Patient summary: We found that rs72725854, an African ancestry–specific risk variant, is more common in men with a family history of prostate cancer and in those diagnosed with prostate cancer at younger ages. Men of African ancestry may benefit from the knowledge of their carrier status for this genetic risk variant to guide decisions about prostate cancer screening. © 2020 The AuthorsThe African ancestry–specific prostate cancer risk variant at 8q24, rs72725854, is enriched in men diagnosed at younger ages and men with a prostate cancer family history. Carriers of this risk allele would benefit from regular and earlier prostate cancer screening

Genome-wide identification of directed gene networks using large-scale population genomics data

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene–gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10−10), among which transcription factors were overrepresented (Fisher’s P = 3.3 × 10−7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser

Autosomal genetic variation is associated with DNA methylation in regions variably escaping X-chromosome inactivation

X-chromosome inactivation (XCI), i.e., the inactivation of one of the female X chromosomes, restores equal expression of X-chromosomal genes between females and males. However, ~10% of genes show variable degrees of escape from XCI between females, although little is known about the causes of variable XCI. Using a discovery data-set of 1867 females and 1398 males and a replication sample of 3351 females, we show that genetic variation at three autosomal loci is associated with female-specific changes in X-chromosome methylation. Through cis-eQTL expression analysis, we map these loci to the genes SMCHD1/METTL4, TRIM6/HBG2, and ZSCAN9. Low-expression alleles of the loci are predominantly associated with mild hypomethylation of CpG islands near genes known to variably escape XCI, implicating the autosomal genes in variable XCI. Together, these results suggest a genetic basis for variable escape from XCI and highlight the potential of a population genomics approach to identify genes involved in XCI

Skewed X-inactivation is common in the general female population

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants

SARS-CoV-2-specific nasal IgA wanes 9 months after hospitalisation with COVID-19 and is not induced by subsequent vaccination

BACKGROUND: Most studies of immunity to SARS-CoV-2 focus on circulating antibody, giving limited insights into mucosal defences that prevent viral replication and onward transmission. We studied nasal and plasma antibody responses one year after hospitalisation for COVID-19, including a period when SARS-CoV-2 vaccination was introduced. METHODS: In this follow up study, plasma and nasosorption samples were prospectively collected from 446 adults hospitalised for COVID-19 between February 2020 and March 2021 via the ISARIC4C and PHOSP-COVID consortia. IgA and IgG responses to NP and S of ancestral SARS-CoV-2, Delta and Omicron (BA.1) variants were measured by electrochemiluminescence and compared with plasma neutralisation data. FINDINGS: Strong and consistent nasal anti-NP and anti-S IgA responses were demonstrated, which remained elevated for nine months (p < 0.0001). Nasal and plasma anti-S IgG remained elevated for at least 12 months (p < 0.0001) with plasma neutralising titres that were raised against all variants compared to controls (p < 0.0001). Of 323 with complete data, 307 were vaccinated between 6 and 12 months; coinciding with rises in nasal and plasma IgA and IgG anti-S titres for all SARS-CoV-2 variants, although the change in nasal IgA was minimal (1.46-fold change after 10 months, p = 0.011) and the median remained below the positive threshold determined by pre-pandemic controls. Samples 12 months after admission showed no association between nasal IgA and plasma IgG anti-S responses (R = 0.05, p = 0.18), indicating that nasal IgA responses are distinct from those in plasma and minimally boosted by vaccination. INTERPRETATION: The decline in nasal IgA responses 9 months after infection and minimal impact of subsequent vaccination may explain the lack of long-lasting nasal defence against reinfection and the limited effects of vaccination on transmission. These findings highlight the need to develop vaccines that enhance nasal immunity. FUNDING: This study has been supported by ISARIC4C and PHOSP-COVID consortia. ISARIC4C is supported by grants from the National Institute for Health and Care Research and the Medical Research Council. Liverpool Experimental Cancer Medicine Centre provided infrastructure support for this research. The PHOSP-COVD study is jointly funded by UK Research and Innovation and National Institute of Health and Care Research. The funders were not involved in the study design, interpretation of data or the writing of this manuscript

Large-scale phenotyping of patients with long COVID post-hospitalization reveals mechanistic subtypes of disease

One in ten severe acute respiratory syndrome coronavirus 2 infections result in prolonged symptoms termed long coronavirus disease (COVID), yet disease phenotypes and mechanisms are poorly understood1. Here we profiled 368 plasma proteins in 657 participants ≥3 months following hospitalization. Of these, 426 had at least one long COVID symptom and 233 had fully recovered. Elevated markers of myeloid inflammation and complement activation were associated with long COVID. IL-1R2, MATN2 and COLEC12 were associated with cardiorespiratory symptoms, fatigue and anxiety/depression; MATN2, CSF3 and C1QA were elevated in gastrointestinal symptoms and C1QA was elevated in cognitive impairment. Additional markers of alterations in nerve tissue repair (SPON-1 and NFASC) were elevated in those with cognitive impairment and SCG3, suggestive of brain–gut axis disturbance, was elevated in gastrointestinal symptoms. Severe acute respiratory syndrome coronavirus 2-specific immunoglobulin G (IgG) was persistently elevated in some individuals with long COVID, but virus was not detected in sputum. Analysis of inflammatory markers in nasal fluids showed no association with symptoms. Our study aimed to understand inflammatory processes that underlie long COVID and was not designed for biomarker discovery. Our findings suggest that specific inflammatory pathways related to tissue damage are implicated in subtypes of long COVID, which might be targeted in future therapeutic trials