Periklis Makrythanasis

Passive and active DNA methylation and the interplay with genetic variation in gene regulation

DNA methylation is an essential epigenetic mark whose role in gene regulation and its dependency on genomic sequence and environment are not fully understood. In this study we provide novel insights into the mechanistic relationships between genetic variation, DNA methylation and transcriptome sequencing data in three different cell-types of the GenCord human population cohort. We find that the association between DNA methylation and gene expression variation among individuals are likely due to different mechanisms from those establishing methylation-expression patterns during differentiation. Furthermore, cell-type differential DNA methylation may delineate a platform in which local inter-individual changes may respond to or act in gene regulation. We show that unlike genetic regulatory variation, DNA methylation alone does not significantly drive allele specific expression. Finally, inferred mechanistic relationships using genetic variation as well as correlations with TF abundance reveal both a passive and active role of DNA methylation to regulatory interactions influencing gene expression. DOI: http://dx.doi.org/10.7554/eLife.00523.001

Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing.

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans’ lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

Mutations in ZMYND10, a Gene Essential for Proper Axonemal Assembly of Inner and Outer Dynein Arms in Humans and Flies, Cause Primary Ciliary Dyskinesia

Primary ciliary dyskinesia (PCD) is a ciliopathy characterized by airway disease, infertility, and laterality defects, often caused by dual loss of the inner dynein arms (IDAs) and outer dynein arms (ODAs), which power cilia and flagella beating. Using whole-exome and candidate-gene Sanger resequencing in PCD-affected families afflicted with combined IDA and ODA defects, we found that 6/38 (16%) carried biallelic mutations in the conserved zinc-finger gene BLU (ZMYND10). ZMYND10 mutations conferred dynein-arm loss seen at the ultrastructural and immunofluorescence level and complete cilia immotility, except in hypomorphic p.Val16Gly (c.47T>G) homozygote individuals, whose cilia retained a stiff and slowed beat. In mice, Zmynd10 mRNA is restricted to regions containing motile cilia. In a Drosophila model of PCD, Zmynd10 is exclusively expressed in cells with motile cilia: chordotonal sensory neurons and sperm. In these cells, P-element-mediated gene silencing caused IDA and ODA defects, proprioception deficits, and sterility due to immotile sperm. Drosophila Zmynd10 with an equivalent c.47T>G (p.Val16Gly) missense change rescued mutant male sterility less than the wild-type did. Tagged Drosophila ZMYND10 is localized primarily to the cytoplasm, and human ZMYND10 interacts with LRRC6, another cytoplasmically localized protein altered in PCD. Using a fly model of PCD, we conclude that ZMYND10 is a cytoplasmic protein required for IDA and ODA assembly and that its variants cause ciliary dysmotility and PCD with laterality defects.

DEPDC5 mutations in families presenting as autosomal dominant nocturnal frontal lobe epilepsy

Objective: To study the prevalence of DEPDC5 mutations in a series of 30 small European families with a phenotype compatible with autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE). Methods: Thirty unrelated families referred with ADNFLE were recruited in France, Italy, Germany, Belgium, and Norway. Whole-exome sequencing was performed in 10 probands and direct sequencing of the DEPDC5 coding sequence in 20 probands. Testing for nonsense-mediated messenger RNA decay (NMD) was performed in lymphoblastic cells. Results: Exome sequencing revealed a splice acceptor mutation (c.2355-2A>G) in DEPDC5 in the proband of a German family. In addition, 3 nonsense DEPDC5 mutations (p.Arg487*, p.Arg1087*, and p.Trp1369*) were detected in the probands of 2 French and one Belgian family. The nonsense mutations p.Arg487* and p.Arg1087* were targeted by NMD, leading to the degradation of the mutated transcripts. At the clinical level, 78% of the patients with DEPDC5 mutations were drug resistant. Conclusions: DEPDC5 loss-of-function mutations were found in 13% of the families with a presentation of ADNFLE. The rate of drug resistance was high in patients with DEPDC5 mutations. Small ADNFLE pedigrees with DEPDC5 mutations might actually represent a part of the broader familial focal epilepsy with variable foci phenotype.

MLL2 mutation detection in 86 patients with Kabuki syndrome: a genotype-phenotype study

Recently, pathogenic variants in the MLL2 gene were identified as the most common cause of Kabuki (Niikawa–Kuroki) syndrome (MIM#147920). To further elucidate the genotype–phenotype correlation, we studied a large cohort of 86 clinically defined patients with Kabuki syndrome (KS) for mutations in MLL2. All patients were assessed using a standardized phenotype list and all were scored using a newly developed clinical score list for KS (MLL2‐Kabuki score 0–10). Sequencing of the full coding region and intron–exon boundaries of MLL2 identified a total of 45 likely pathogenic mutations (52%): 31 nonsense, 10 missense and four splice‐site mutations, 34 of which were novel. In five additional patients, novel, i.e. non‐dbSNP132 variants of clinically unknown relevance, were identified. Patients with likely pathogenic nonsense or missense MLL2 mutations were usually more severely affected (median ‘MLL2‐Kabuki score’ of 6) as compared to the patients without MLL2 mutations (median ‘MLL2‐Kabuki score’ of 5), a significant difference (p < 0.0014). Several typical facial features such as large dysplastic ears, arched eyebrows with sparse lateral third, blue sclerae, a flat nasal tip with a broad nasal root, and a thin upper and a full lower lip were observed more often in mutation positive patients.

Diagnostic exome sequencing to elucidate the genetic basis of likely recessive disorders in consanguineous families.

Rare, atypical, and undiagnosed autosomal‐recessive disorders frequently occur in the offspring of consanguineous couples. Current routine diagnostic genetic tests fail to establish a diagnosis in many cases. We employed exome sequencing to identify the underlying molecular defects in patients with unresolved but putatively autosomal‐recessive disorders in consanguineous families and postulated that the pathogenic variants would reside within homozygous regions. Fifty consanguineous families participated in the study, with a wide spectrum of clinical phenotypes suggestive of autosomal‐recessive inheritance, but with no definitive molecular diagnosis. DNA samples from the patient(s), unaffected sibling(s), and the parents were genotyped with a 720K SNP array. Exome sequencing and array CGH (comparative genomic hybridization) were then performed on one affected individual per family. High‐confidence pathogenic variants were found in homozygosity in known disease‐causing genes in 18 families (36%) (one by array CGH and 17 by exome sequencing), accounting for the clinical phenotype in whole or in part. In the remainder of the families, no causative variant in a known pathogenic gene was identified. Our study shows that exome sequencing, in addition to being a powerful diagnostic tool, promises to rapidly expand our knowledge of rare genetic Mendelian disorders and can be used to establish more detailed causative links between mutant genotypes and clinical phenotypes.

Efficient targeted transcript discovery via array-based normalization of RACE libraries

Rapid amplification of cDNA ends (RACE) is a widely used approach for transcript identification. Random clone selection from the RACE mixture, however, is an ineffective sampling strategy if the dynamic range of transcript abundances is large. To improve sampling efficiency of human transcripts, we hybridized the products of the RACE reaction onto tiling arrays and used the detected exons to delineate a series of reverse-transcriptase (RT)-PCRs, through which the original RACE transcript population was segregated into simpler transcript populations. We independently cloned the products and sequenced randomly selected clones. This approach, RACEarray, is superior to direct cloning and sequencing of RACE products because it specifically targets new transcripts and often results in overall normalization of transcript abundance. We show theoretically and experimentally that this strategy leads indeed to efficient sampling of new transcripts, and we investigated multiplexing the strategy by pooling RACE reactions from multiple interrogated loci before hybridization.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.

EXOSC3 mutations in pontocerebellar hypoplasia type 1: novel mutations and genotype-phenotype correlations

BackgroundPontocerebellar hypoplasia (PCH) represents a group of neurodegenerative disorders with prenatal onset. Eight subtypes have been described thus far (PCH1-8) based on clinical and genetic features. Common characteristics include hypoplasia and atrophy of the cerebellum, variable pontine atrophy, and severe mental and motor impairments. PCH1 is distinctly characterized by the combination with degeneration of spinal motor neurons. Recently, mutations in the exosome component 3 gene (EXOSC3) have been identified in approximately half of the patients with PCH subtype 1.MethodsWe selected a cohort of 99 PCH patients (90 families) tested negative for mutations in the TSEN genes, RARS2, VRK1 and CASK. Patients in this cohort were referred with a tentative diagnose PCH type 1, 2, 4, 7 or unclassified PCH. Genetic analysis of the EXOSC3 gene was performed using Sanger sequencing. Clinical data, MR images and autopsy reports of patients positive for EXOSC3 mutations were analyzed.ResultsEXOSC3 mutations were found in twelve families with PCH subtype 1, and were not found in patients with other PCH subtypes. Identified mutations included a large deletion, nonsense and missense mutations. Examination of clinical data reveals a prolonged disease course in patients with a homozygous p.D132A mutation. MRI shows variable pontine hypoplasia in EXOSC3 mediated PCH, where the pons is largely preserved in patients with a homozygous p.D132A mutation, but attenuated in patients with other mutations. Additionally, bilateral cerebellar cysts were found in patients compound heterozygous for a p.D132A mutation and a nonsense allele.ConclusionsEXOSC3 mediated PCH shows clear genotype-phenotype correlations. A homozygous p.D132A mutation leads to PCH with possible survival into early puberty, and preservation of the pons. Compound heterozygosity for a p.D132A mutation and a nonsense or p.Y109N allele, a homozygous p.G31A mutation or a p.G135E mutation causes a more rapidly progressive course leading to death in infancy and attenuation of the ventral pons.Our findings imply a clear correlation between genetic mutation and clinical outcome in EXOSC3 mediated PCH, including variable involvement of the pons.

Pathogenic Variants in PIGG Cause Intellectual Disability with Seizures and Hypotonia

Glycosylphosphatidylinositol (GPI) is a glycolipid that anchors >150 various proteins to the cell surface. At least 27 genes are involved in biosynthesis and transport of GPI-anchored proteins (GPI-APs). To date, mutations in 13 of these genes are known to cause inherited GPI deficiencies (IGDs), and all are inherited as recessive traits. IGDs mainly manifest as intellectual disability, epilepsy, coarse facial features, and multiple organ anomalies. These symptoms are caused by the decreased surface expression of GPI-APs or by structural abnormalities of GPI. Here, we present five affected individuals (from two consanguineous families from Egypt and Pakistan and one non-consanguineous family from Japan) who show intellectual disability, hypotonia, and early-onset seizures. We identified pathogenic variants in PIGG, a gene in the GPI pathway. In the consanguineous families, homozygous variants c.928C>T (p.Gln310(∗)) and c.2261+1G>C were found, whereas the Japanese individual was compound heterozygous for c.2005C>T (p.Arg669Cys) and a 2.4 Mb deletion involving PIGG. PIGG is the enzyme that modifies the second mannose with ethanolamine phosphate, which is removed soon after GPI is attached to the protein. Physiological significance of this transient modification has been unclear. Using B lymphoblasts from affected individuals of the Egyptian and Japanese families, we revealed that PIGG activity was almost completely abolished; however, the GPI-APs had normal surface levels and normal structure, indicating that the pathogenesis of PIGG deficiency is not yet fully understood. The discovery of pathogenic variants in PIGG expands the spectrum of IGDs and further enhances our understanding of this etiopathogenic class of intellectual disability.