Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability

Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown—as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

Redistribution of H3K27me3 upon DNA hypomethylation results in de-repression of Polycomb target genes

BACKGROUND

DNA methylation and the Polycomb repression system are epigenetic mechanisms that play important roles in maintaining transcriptional repression. Recent evidence suggests that DNA methylation can attenuate the binding of Polycomb protein components to chromatin and thus plays a role in determining their genomic targeting. However, whether this role of DNA methylation is important in the context of transcriptional regulation is unclear.

RESULTS

By genome-wide mapping of the Polycomb Repressive Complex 2-signature histone mark, H3K27me3, in severely DNA hypomethylated mouse somatic cells, we show that hypomethylation leads to widespread H3K27me3 redistribution, in a manner that reflects the local DNA methylation status in wild-type cells. Unexpectedly, we observe striking loss of H3K27me3 and Polycomb Repressive Complex 2 from Polycomb target gene promoters in DNA hypomethylated cells, including Hox gene clusters. Importantly, we show that many of these genes become ectopically expressed in DNA hypomethylated cells, consistent with loss of Polycomb-mediated repression.

CONCLUSIONS

An intact DNA methylome is required for appropriate Polycomb-mediated gene repression by constraining Polycomb Repressive Complex 2 targeting. These observations identify a previously unappreciated role for DNA methylation in gene regulation and therefore influence our understanding of how this epigenetic mechanism contributes to normal development and disease.

Genome-wide methylation profiling in Crohn's disease identifies altered epigenetic regulation of key host defense mechanisms including the Th17 pathway

BACKGROUND

Germline variation in the 71 Crohn's disease (CD) loci implicated by genome-wide association studies (GWAS) only accounts for approximately 25% of estimated heritability. The contribution of epigenetic alterations to disease pathogenesis is emerging as a research priority.

MATERIALS AND METHODS

The methylation status of 27,578 CpG sites across the genome was analyzed using the Illumina Human Methylation27 assay in DNA extracted from whole blood samples from 40 adult females (21 ileal CD, 19 healthy controls) and 16 girls with childhood-onset CD, all nonsmokers. Our primary analysis compared methylation profiles in adult cases and controls.

RESULTS

Our data define a global methylation profile characteristic of ileal CD. In all, 1117 sites were differentially methylated (corrected P < 0.01); 50 showed significantly altered methylation in cases compared with controls (uncorrected P < 10(-6), corrected P < 0.0006), including genes altering immune activation: MAPK13, FASLG, PRF1, S100A13, RIPK3, and IL-21R. Gene ontology analyses implicated immunity-related pathways as targets of epigenetic modification (immune system processes [P = 1.3 × 10(-22)], immune response [P = 8.1 × 10(-16)], defense responses to bacteria [P = 1.8 × 10(-15)]). Ingenuity canonical pathway analyses implicated dendritic cell activity (P = 2.4 × 10(-8)) and differential regulation of cytokines by interleukin (IL)-17A and IL-17F (P = 5.8 × 10(-7)). We identified a significant enrichment of methylation changes within 50 kb of CD GWAS loci (8.6-fold [P = 0.021] in adults; 2.4-fold [P = 0.009] in adults and children combined), including IL-27, IL-19, TNF, MST1, and NOD2. Methylation status was predictive of disease status (sensitivity 0.71, specificity 0.83). Disease activity, drug therapy, NOD2 and DNMT3A genotypes were not associated with methylation changes.

CONCLUSIONS

These data provide an important insight into the impact of epigenetic mechanisms in the pathogenesis of CD.

Abundant pleiotropy in human complex diseases and traits

We present a systematic review of pleiotropy among SNPs and genes reported to show genome-wide association with common complex diseases and traits. We find abundant evidence of pleiotropy; 233 (16.9%) genes and 77 (4.6%) SNPs show pleiotropic effects. SNP pleiotropic status was associated with gene location (p = 0.024; pleiotropic SNPs more often exonic [14.5% versus 4.9% for nonpleiotropic, trait-associated SNPs] and less often intergenic [15.8% versus 23.6%]), "predicted transcript consequence" (p = 0.001; pleiotropic SNPs more often predicted to be structurally deleterious [5% versus 0.4%] but not more often in regulatory sequences), and certain disease classes. We develop a method to calculate the likelihood that pleiotropic links between traits occurred more often than expected and demonstrate that this approach can identify etiological links that are already known (such as between fetal hemoglobin and malaria risk) and those that are not yet established (e.g., between plasma campesterol levels and gallstones risk; and between immunoglobulin A and juvenile idiopathic arthritis). Examples of pleiotropy will accumulate over time, but it is already clear that pleiotropy is a common property of genes and SNPs associated with disease traits, and this will have implications for identification of molecular targets for drug development, future genetic risk-profiling, and classification of diseases.