Margaret A. Taub

Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis

Epigenetic mechanisms integrate genetic and environmental causes of disease, but comprehensive genome-wide analyses of epigenetic modifications have not yet demonstrated robust association with common diseases. Using Illumina HumanMethylation450 arrays on 354 anti-citrullinated protein antibody–associated rheumatoid arthritis cases and 337 controls, we identified two clusters within the major histocompatibility complex (MHC) region whose differential methylation potentially mediates genetic risk for rheumatoid arthritis. To reduce confounding factors that have hampered previous epigenome-wide studies, we corrected for cellular heterogeneity by estimating and adjusting for cell-type proportions in our blood-derived DNA samples and used mediation analysis to filter out associations likely to be a consequence of disease. Four CpGs also showed an association between genotype and variance of methylation. The associations for both clusters replicated at least one CpG (P < 0.01), with the rest showing suggestive association, in monocyte cell fractions in an independent cohort of 12 cases and 12 controls. Thus, DNA methylation is a potential mediator of genetic risk.

Confirming genes influencing risk to cleft lip with/without cleft palate in a case–parent trio study

A collection of 1,108 case–parent trios ascertained through an isolated, nonsyndromic cleft lip with or without cleft palate (CL/P) was used to replicate the findings from a genome-wide association study (GWAS) conducted by Beaty et al. (Nat Genet 42:525–529, 2010), where four different genes/regions were identified as influencing risk to CL/P. Tagging SNPs for 33 different genes were genotyped (1,269 SNPs). All four of the genes originally identified as showing genome-wide significance (IRF6, ABCA4 and MAF, plus the 8q24 region) were confirmed in this independent sample of trios (who were primarily of European and Southeast Asian ancestry). In addition, eight genes classified as ‘second tier’ hits in the original study (PAX7, THADA, COL8A1/FILIP1L, DCAF4L2, GADD45G, NTN1, RBFOX3 and FOXE1) showed evidence of linkage and association in this replication sample. Meta-analysis between the original GWAS trios and these replication trios showed PAX7, COL8A1/FILIP1L and NTN1 achieved genome-wide significance. Tests for gene–environment interaction between these 33 genes and maternal smoking found evidence for interaction with two additional genes: GRID2 and ELAVL2 among European mothers (who had a higher rate of smoking than Asian mothers). Formal tests for gene–gene interaction (epistasis) failed to show evidence of statistical interaction in any simple fashion. This study confirms that many different genes influence risk to CL/P.

Overcoming bias and systematic errors in next generation sequencing data

Considerable time and effort has been spent in developing analysis and quality assessment methods to allow the use of microarrays in a clinical setting. As is the case for microarrays and other high-throughput technologies, data from new high-throughput sequencing technologies are subject to technological and biological biases and systematic errors that can impact downstream analyses. Only when these issues can be readily identified and reliably adjusted for will clinical applications of these new technologies be feasible. Although much work remains to be done in this area, we describe consistently observed biases that should be taken into account when analyzing high-throughput sequencing data. In this article, we review current knowledge about these biases, discuss their impact on analysis results, and propose solutions.

Identification of Functional Variants for Cleft Lip with or without Cleft Palate in or near PAX7, FGFR2, and NOG by Targeted Sequencing of GWAS Loci

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

Measuring cell-type specific differential methylation in human brain tissue

The behavior of epigenetic mechanisms in the brain is obscured by tissue heterogeneity and disease-related histological changes. Not accounting for these confounders leads to biased results. We develop a statistical methodology that estimates and adjusts for celltype composition by decomposing neuronal and non-neuronal differential signal. This method provides a conceptual framework for deconvolving heterogeneous epigenetic data from postmortem brain studies. We apply it to find cell-specific differentially methylated regions between prefrontal cortex and hippocampus. We demonstrate the utility of the method on both Infinium 450k and CHARM data.

Association of DNA Methylation Differences With Schizophrenia in an Epigenome-Wide Association Study

IMPORTANCE DNA methylation may play an important role in schizophrenia (SZ), either directly as a mechanism of pathogenesis or as a biomarker of risk. OBJECTIVE To scan genome-wide DNA methylation data to identify differentially methylated CpGs between SZ cases and controls. DESIGN, SETTING, AND PARTICIPANTS Epigenome-wide association study begun in 2008 using DNA methylation levels of 456 513 CpG loci measured on the Infinium HumanMethylation450 array (Illumina) in a consortium of case-control studies for initial discovery and in an independent replication set. Primary analyses used general linear regression, adjusting for age, sex, race/ethnicity, smoking, batch, and cell type heterogeneity. The discovery set contained 689 SZ cases and 645 controls (n = 1334), from 3 multisite consortia: the Consortium on the Genetics of Endophenotypes in Schizophrenia, the Project among African-Americans To Explore Risks for Schizophrenia, and the Multiplex Multigenerational Family Study of Schizophrenia. The replication set contained 247 SZ cases and 250 controls (n = 497) from the Genomic Psychiatry Cohort. MAIN OUTCOMES AND MEASURES Identification of differentially methylated positions across the genome in SZ cases compared with controls. RESULTS Of the 689 case participants in the discovery set, 477 (69%) were men and 258 (37%) were non-African American; of the 645 controls, 273 (42%) were men and 419 (65%) were non-African American. In our replication set, cases/controls were 76% male and 100% non-African American. We identified SZ-associated methylation differences at 923 CpGs in the discovery set (false discovery rate, <0.2). Of these, 625 showed changes in the same direction including 172 with P < .05 in the replication set. Some replicated differentially methylated positions are located in a top-ranked SZ region from genome-wide association study analyses. CONCLUSIONS AND RELEVANCE This analysis identified 172 replicated new associations with SZ after careful correction for cell type heterogeneity and other potential confounders. The overlap with previous genome-wide association study data can provide potential insights into the functional relevance of genetic signals for SZ.

Reproducible RNA-seq analysis using recount2

c 16. Köster, J. & Rahmann, S. Bioinformatics 28, 2520– 2522 (2012). 17. Di Tommaso, P. et al. PeerJ 3, e1273 (2015). 18. Goecks, J., Nekrutenko, A. & Taylor, J. Genome Biol. 11, R86 (2010). 19. Blankenberg, D. et al. Genome Biol. 15, 403 (2014). 20. Vivian, J. et al. Preprint at bioRxiv http://biorxiv.org/ content/early/2016/07/07/062497 (2016). 21. Stamatakis, A. Bioinformatics 22, 2688–2690 (2006). 22. Byron, S.A., Van Keuren-Jensen, K.R., Engelthaler, D.M., Carpten, J.D. & Craig, D.W. Nat. Rev. Genet. 17, 257–271 (2016).

Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.

A continuum of admixture in the Western Hemisphere revealed by the African Diaspora genome

The African Diaspora in the Western Hemisphere represents one of the largest forced migrations in history and had a profound impact on genetic diversity in modern populations. To date, the fine-scale population structure of descendants of the African Diaspora remains largely uncharacterized. Here we present genetic variation from deeply sequenced genomes of 642 individuals from North and South American, Caribbean and West African populations, substantially increasing the lexicon of human genomic variation and suggesting much variation remains to be discovered in African-admixed populations in the Americas. We summarize genetic variation in these populations, quantifying the postcolonial sex-biased European gene flow across multiple regions. Moreover, we refine estimates on the burden of deleterious variants carried across populations and how this varies with African ancestry. Our data are an important resource for empowering disease mapping studies in African-admixed individuals and will facilitate gene discovery for diseases disproportionately affecting individuals of African ancestry.

Cooperation between Referees and Authors Increases Peer Review Accuracy

Peer review is fundamentally a cooperative process between scientists in a community who agree to review each others work in an unbiased fashion. Peer review is the foundation for decisions concerning publication in journals, awarding of grants, and academic promotion. Here we perform a laboratory study of open and closed peer review based on an online game. We show that when reviewer behavior was made public under open review, reviewers were rewarded for refereeing and formed significantly more cooperative interactions (13% increase in cooperation, P = 0.018). We also show that referees and authors who participated in cooperative interactions had an 11% higher reviewing accuracy rate (P = 0.016). Our results suggest that increasing cooperation in the peer review process can lead to a decreased risk of reviewing errors.