DNA Methylation Profiling of Uniparental Disomy Subjects Provides a Map of Parental Epigenetic Bias in the Human Genome

Identification of rare de novo epigenetic variations in congenital disorders

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.

XX Disorder of Sex Development is associated with an insertion on chromosome 9 and downregulation of RSPO1 in dogs (Canis lupus familiaris)

Remarkable progress has been achieved in understanding the mechanisms controlling sex determination, yet the cause for many Disorders of Sex Development (DSD) remains unknown. Of particular interest is a rare XX DSD subtype in which individuals are negative for SRY, the testis determining factor on the Y chromosome, yet develop testes or ovotestes, and both of these phenotypes occur in the same family. This is a naturally occurring disorder in humans (Homo sapiens) and dogs (C. familiaris). Phenotypes in the canine XX DSD model are strikingly similar to those of the human XX DSD subtype. The purposes of this study were to identify 1) a variant associated with XX DSD in the canine model and 2) gene expression alterations in canine embryonic gonads that could be informative to causation. Using a genome wide association study (GWAS) and whole genome sequencing (WGS), we identified a variant on C. familiaris autosome 9 (CFA9) that is associated with XX DSD in the canine model and in affected purebred dogs. This is the first marker identified for inherited canine XX DSD. It lies upstream of SOX9 within the canine ortholog for the human disorder, which resides on 17q24. Inheritance of this variant indicates that XX DSD is a complex trait in which breed genetic background affects penetrance. Furthermore, the homozygous variant genotype is associated with embryonic lethality in at least one breed. Our analysis of gene expression studies (RNA-seq and PRO-seq) in embryonic gonads at risk of XX DSD from the canine model identified significant RSPO1 downregulation in comparison to XX controls, without significant upregulation of SOX9 or other known testis pathway genes. Based on these data, a novel mechanism is proposed in which molecular lesions acting upstream of RSPO1 induce epigenomic gonadal mosaicism.

Detection of Imprinted Genes by Single-Cell Allele-Specific Gene Expression

Genomic imprinting results in parental-specific gene expression. Imprinted genes are involved in the etiology of rare syndromes and have been associated with common diseases such as diabetes and cancer. Standard RNA bulk cell sequencing applied to whole-tissue samples has been used to detect imprinted genes in human and mouse models. However, lowly expressed genes cannot be detected by using RNA bulk approaches. Here, we report an original and robust method that combines single-cell RNA-seq and whole-genome sequencing into an optimized statistical framework to analyze genomic imprinting in specific cell types and in different individuals. Using samples from the probands of 2 family trios and 3 unrelated individuals, 1,084 individual primary fibroblasts were RNA sequenced and more than 700,000 informative heterozygous single-nucleotide variations (SNVs) were genotyped. The allele-specific coverage per gene of each SNV in each single cell was used to fit a beta-binomial distribution to model the likelihood of a gene being expressed from one and the same allele. Genes presenting a significant aggregate allelic ratio (between 0.9 and 1) were retained to identify of the allelic parent of origin. Our approach allowed us to validate the imprinting status of all of the known imprinted genes expressed in fibroblasts and the discovery of nine putative imprinted genes, thereby demonstrating the advantages of single-cell over bulk RNA-seq to identify imprinted genes. The proposed single-cell methodology is a powerful tool for establishing a cell type-specific map of genomic imprinting.

Loss of heterozygosity detected at three short tandem repeat locus commonly used for human DNA identification in a case of paternity testing

Short tandem repeat (STR) is widely used for DNA profiling in forensic sciences for its stable inheritance. Genomic variations in STR loci may affect the results of the genotyping. In this study, using STR profiling and genome-wide chromosomal microarray assay, we detected the incidence of uniparental disomy or copy-neutral loss of heterozygosity (LOH) in a case of a parental testing, which altered the genotype of three commonly used STR markers including D2S1338, D2S441 and D2S1776. To the best of our knowledge, this is the first time found that LOH affect the genotyping of STR markers commonly used for paternity testing. Our findings demonstrated that the incidence of LOH in the genome may dramatically alter the results of DNA identification, and suggested that genomic structure variation need to be taking into consideration in the DNA identification using STR markers.

Cumulative Impact of Polychlorinated Biphenyl and Large Chromosomal Duplications on DNA Methylation, Chromatin, and Expression of Autism Candidate Genes

Rare variants enriched for functions in chromatin regulation and neuronal synapses have been linked to autism. How chromatin and DNA methylation interact with environmental exposures at synaptic genes in autism etiologies is currently unclear. Using whole-genome bisulfite sequencing in brain tissue and a neuronal cell culture model carrying a 15q11.2-q13.3 maternal duplication, we find that significant global DNA hypomethylation is enriched over autism candidate genes and affects gene expression. The cumulative effect of multiple chromosomal duplications and exposure to the pervasive persistent organic pollutant PCB 95 altered methylation of more than 1,000 genes. Hypomethylated genes were enriched for H2A.Z, increased maternal UBE3A in Dup15q corresponded to reduced levels of RING1B, and bivalently modified H2A.Z was altered by PCB 95 and duplication. These results demonstrate the compounding effects of genetic and environmental insults on the neuronal methylome that converge upon dysregulation of chromatin and synaptic genes.

Narrowing the FOXF1 distant enhancer region on 16q24.1 critical for ACDMPV

Background

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disorder caused by heterozygous point mutations or genomic deletions involving FOXF1 or its 60-kb tissue-specific enhancer region mapping 270 kb upstream of FOXF1 and involving fetal lung-expressed long non-coding RNA genes and CpG-enriched sites. Recently, we have proposed that the FOXF1 locus at 16q24.1 may be a subject of genomic imprinting.

Findings

Using custom-designed aCGH and Sanger sequencing, we have identified a novel de novo 104 kb genomic deletion upstream of FOXF1 in a patient with histopathologically verified full phenotype of ACDMPV. This deletion allowed us to further narrow the FOXF1 enhancer region and identify its critical 15-kb core interval, essential for lung development. This interval harbors binding sites for lung-expressed transcription factors, including GATA3, ESR1, and YY1, and is flanked by the lncRNA genes and CpG islands. Bisulfite sequencing of one of these CpG islands on the non-deleted allele showed that it is predominantly methylated on the maternal chromosome 16.

Conclusions

Substantial narrowing and bisulfite sequencing of the FOXF1 enhancer region on 16q24.1 provided new insights into its regulatory function and genomic imprinting.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0278-2) contains supplementary material, which is available to authorized users.