The landscape of genomic imprinting across diverse adult human tissues

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

Investigating the potential of single-cell DNA methylation data to detect allele-specific methylation and imprinting

Allele-specific methylation (ASM) is an epigenetic modification whereby one parental allele becomes methylated and the other unmethylated at a specific locus. ASM is most often driven by the presence of nearby heterozygous variants that influence methylation, but also occurs somatically in the context of genomic imprinting. In this study, we investigate ASM using publicly available single-cell reduced representation bisulfite sequencing (scRRBS) data on 608 B cells sampled from six healthy B cell samples and 1,230 cells from 11 chronic lymphocytic leukemia (CLL) samples. We developed a likelihood-based criterion to test whether a CpG exhibited ASM, based on the distributions of methylated and unmethylated reads both within and across cells. Applying our likelihood ratio test, 65,998 CpG sites exhibited ASM in healthy B cell samples according to a Bonferroni criterion (p < 8.4 × 10-9), and 32,862 CpG sites exhibited ASM in CLL samples (p < 8.5 × 10-9). We also called ASM at the sample level. To evaluate the accuracy of our method, we called heterozygous variants from the scRRBS data, which enabled variant-based calls of ASM within each cell. Comparing sample-level ASM calls to the variant-based measures of ASM, we observed a positive predictive value of 76%-100% across samples. We observed high concordance of ASM across samples and an overrepresentation of ASM in previously reported imprinted genes and genes with imprinting binding motifs. Our study demonstrates that single-cell bisulfite sequencing is a potentially powerful tool to investigate ASM, especially as studies expand to increase the number of samples and cells sequenced.

A novel imprinted locus on bovine chromosome 18 homologous with human chromosome 16q24.1

Genomic imprinting is an epigenetic regulation mechanism in mammals resulting in the parentally dependent monoallelic expression of genes. Imprinting disorders in humans are associated with several congenital syndromes and cancers and remain the focus of many medical studies. Cattle is a better model organism for investigating human embryo development than mice. Imprinted genes usually cluster on chromosomes and are regulated by different methylation regions (DMRs) located in imprinting control regions that control gene expression in cis. There is an imprinted locus on human chromosome 16q24.1 associated with congenital lethal developmental lung disease in newborns. However, genomic imprinting on bovine chromosome 18, which is homologous with human chromosome 16 has not been systematically studied. The aim of this study was to analyze the allelic expressions of eight genes (CDH13, ATP2C2, TLDC1, COTL1, CRISPLD2, ZDHHC7, KIAA0513, and GSE1) on bovine chromosome 18 and to search the DMRs associated gene allelic expression. Three transcript variants of the ZDHHC7 gene (X1, X2, and X5) showed maternal imprinting in bovine placentas. In addition, the monoallelic expression of X2 and X5 was tissue-specific. Five transcripts of the KIAA0513 gene showed tissue- and isoform-specific monoallelic expression. The CDH13, ATP2C2, and TLDC1 genes exhibited tissue-specific imprinting, however, COTL1, CRISLPLD2, and GSE1 escaped imprinting. Four DMRs, established after fertilization, were found in this region. Two DMRs were located between the ZDHHC7 and KIAA0513 genes, and two were in exon 1 of the CDH13 and ATP2C2 genes, respectively. The results from this study support future studies on the molecular mechanism to regulate the imprinting of candidate genes on bovine chromosome 18.

Biallelic SOX8 Variants Associated With Novel Syndrome With Myopathy, Skeletal Deformities, Intellectual Disability, and Ovarian Dysfunction

Background and Objectives

The human genome contains ∼20,000 genes, each of which has its own set of complex regulatory systems to govern precise expression in each developmental stage and cell type. Here, we report a female patient with congenital weakness, respiratory failure, skeletal dysplasia, contractures, short stature, intellectual delay, respiratory failure, and amenorrhea who presented to Medical Genetics service with no known cause for her condition.

Methods

Whole-exome and whole-genome sequencing were conducted, as well as investigational functional studies to assess the effect of SOX8 variant.

Results

The patient was found to have biallelic SOX8 variants (NM_014587.3:c.422+5G>C; c.583dup p.(His195ProfsTer11)). SOX8 is a transcriptional regulator, which is predicted to be imprinted (expressed from only one parental allele), but this has not yet been confirmed. We provide evidence that while SOX8 was maternally expressed in adult-derived fibroblasts and lymphoblasts, it was biallelically expressed in other cell types and therefore suggest that biallelic variants are associated with this recessive condition. Functionally, we showed that the paternal variant had the capacity to affect mRNA splicing while the maternal variant resulted in low levels of a truncated protein, which showed decreased binding at and altered expression of SOX8 targets.

Discussion

Our findings associate SOX8 variants with this novel condition, highlight how complex genome regulation can complicate novel disease-gene identification, and provide insight into the molecular pathogenesis of this disease.

Imprinting disorders

Imprinting disorders (ImpDis) are congenital conditions that are characterized by disturbances of genomic imprinting. The most common individual ImpDis are Prader-Willi syndrome, Angelman syndrome and Beckwith-Wiedemann syndrome. Individual ImpDis have similar clinical features, such as growth disturbances and developmental delay, but the disorders are heterogeneous and the key clinical manifestations are often non-specific, rendering diagnosis difficult. Four types of genomic and imprinting defect (ImpDef) affecting differentially methylated regions (DMRs) can cause ImpDis. These defects affect the monoallelic and parent-of-origin-specific expression of imprinted genes. The regulation within DMRs as well as their functional consequences are mainly unknown, but functional cross-talk between imprinted genes and functional pathways has been identified, giving insight into the pathophysiology of ImpDefs. Treatment of ImpDis is symptomatic. Targeted therapies are lacking owing to the rarity of these disorders; however, personalized treatments are in development. Understanding the underlying mechanisms of ImpDis, and improving diagnosis and treatment of these disorders, requires a multidisciplinary approach with input from patient representatives.

Genomic-transcriptomic evolution in lung cancer and metastasis

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Widespread allele-specific topological domains in the human genome are not confined to imprinted gene clusters

BACKGROUND

There is widespread interest in the three-dimensional chromatin conformation of the genome and its impact on gene expression. However, these studies frequently do not consider parent-of-origin differences, such as genomic imprinting, which result in monoallelic expression. In addition, genome-wide allele-specific chromatin conformation associations have not been extensively explored. There are few accessible bioinformatic workflows for investigating allelic conformation differences and these require pre-phased haplotypes which are not widely available.

RESULTS

We developed a bioinformatic pipeline, "HiCFlow," that performs haplotype assembly and visualization of parental chromatin architecture. We benchmarked the pipeline using prototype haplotype phased Hi-C data from GM12878 cells at three disease-associated imprinted gene clusters. Using Region Capture Hi-C and Hi-C data from human cell lines (1-7HB2, IMR-90, and H1-hESCs), we can robustly identify the known stable allele-specific interactions at the IGF2-H19 locus. Other imprinted loci (DLK1 and SNRPN) are more variable and there is no "canonical imprinted 3D structure," but we could detect allele-specific differences in A/B compartmentalization. Genome-wide, when topologically associating domains (TADs) are unbiasedly ranked according to their allele-specific contact frequencies, a set of allele-specific TADs could be defined. These occur in genomic regions of high sequence variation. In addition to imprinted genes, allele-specific TADs are also enriched for allele-specific expressed genes. We find loci that have not previously been identified as allele-specific expressed genes such as the bitter taste receptors (TAS2Rs).

CONCLUSIONS

This study highlights the widespread differences in chromatin conformation between heterozygous loci and provides a new framework for understanding allele-specific expressed genes.

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our toolchain, which consists of additions to the VG toolkit and a standalone tool, RPVG, can construct spliced pangenome graphs, map RNA sequencing data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. We show that this workflow improves accuracy over state-of-the-art RNA sequencing mapping methods, and that it can efficiently quantify haplotype-specific transcript expression without needing to characterize the haplotypes of a sample beforehand.

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our toolchain, which consists of additions to the VG toolkit and a standalone tool, RPVG, can construct spliced pangenome graphs, map RNA sequencing data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. We show that this workflow improves accuracy over state-of-the-art RNA sequencing mapping methods, and that it can efficiently quantify haplotype-specific transcript expression without needing to characterize the haplotypes of a sample beforehand.

Random allelic expression in the adult human body

Genes are typically assumed to express both parental alleles similarly, yet cell lines show random allelic expression (RAE) for many autosomal genes that could shape genetic effects. Thus, understanding RAE in human tissues could improve our understanding of phenotypic variation. Here, we develop a methodology to perform genome-wide profiling of RAE and biallelic expression in GTEx datasets for 832 people and 54 tissues. We report 2,762 autosomal genes with some RAE properties similar to randomly inactivated X-linked genes. We found that RAE is associated with rapidly evolving regions in the human genome, adaptive signaling processes, and genes linked to age-related diseases such as neurodegeneration and cancer. We define putative mechanistic subtypes of RAE distinguished by gene overlaps on sense and antisense DNA strands, aggregation in clusters near telomeres, and increased regulatory complexity and inputs compared with biallelic genes. We provide foundations to study RAE in human phenotypes, evolution, and disease.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Beckwith-Wiedemann syndrome with long QT caused by a deletion involving KCNQ1 but not KCNQ1OT1:TSS-DMR

Beckwith-Wiedemann syndrome (BWS) is an imprinting disorder with characteristic features, such as overgrowth, macroglossia, and exomphalos. Hypomethylation of the KCNQ1OT1:TSS-differentially methylated region (DMR) on the 11p15.5 imprinted region is the most common etiology of BWS. KCNQ1 on 11p15.5 is expressed from the maternally inherited allele in most tissues, but is biparentally expressed in the heart, and maternal KCNQ1 transcription is required to establish the maternal DNA imprint in the KCNQ1OT1:TSS-DMR. Loss of function variants in KCNQ1 result in long QT syndrome type 1 (LQT1). To date, eight patients with BWS due to KCNQ1 splice variants or structural abnormalities involving KCNQ1 but not the KCNQ1OT1:TSS-DMR have been reported (KCNQ1-BWS), and four of them had LQT1. We report a Japanese boy with BWS and LQT1 presenting with extreme hypomethylation of the KCNQ1OT1:TSS-DMR caused by a de novo 215-kb deletion including KCNQ1 but not the KCNQ1OT1:TSS-DMR on the maternal allele. He was born by emergency cesarean section due to suspicion of placental abruption at 30 weeks of gestation. His birth weight and length were +1.6 SD and +1.0 SD, respectively. His placental weight was +3.9 SD, and histological examination of his placenta was consistent with mesenchymal dysplasia. He had BWS clinical features, including macroglossia, ear creases and pits, body asymmetry, and rectus abdominis muscle dehiscence, and BWS was therefore diagnosed. LQT1 was first noticed at three months in a preoperative examination for lingual frenectomy. The summarized data of our patient and the previously reported eight patients in KCNQ1-BWS showed more frequent and earlier preterm births and smaller sized birth weight in KCNQ1-BWS cases than those with BWS caused by epimutation of the KCNQ1OT1:TSS-DMR. In addition, in five of nine patients with KCNQ1-BWS, LQT1 was detected, and two of them were identified at school age. In our patient and in another single case with LQT1, the LQT1 was not detected early despite neonatal ECG monitoring. For BWS patients with extreme hypomethylation of the KCNQ1OT1:TSS-DMR, searching for CNVs involving KCNQ1 and mutation screening for KCNQ1 should be considered together with periodic ECG monitoring. (338/500 words).

Increased copy number of imprinted genes in the chromosomal region 20q11-q13.32 is associated with resistance to antitumor agents in cancer cell lines

BACKGROUND

Parent of origin-specific allelic expression of imprinted genes is epigenetically controlled. In cancer, imprinted genes undergo both genomic and epigenomic alterations, including frequent copy number changes. We investigated whether copy number loss or gain of imprinted genes in cancer cell lines is associated with response to chemotherapy treatment.

RESULTS

We analyzed 198 human imprinted genes including protein-coding genes and noncoding RNA genes using data from tumor cell lines from the Cancer Cell Line Encyclopedia and Genomics of Drug Sensitivity in Cancer datasets. We examined whether copy number of the imprinted genes in 35 different genome locations was associated with response to cancer drug treatment. We also analyzed associations of pretreatment expression and DNA methylation of imprinted genes with drug response. Higher copy number of BLCAP, GNAS, NNAT, GNAS-AS1, HM13, MIR296, MIR298, and PSIMCT-1 in the chromosomal region 20q11-q13.32 was associated with resistance to multiple antitumor agents. Increased expression of BLCAP and HM13 was also associated with drug resistance, whereas higher methylation of gene regions of BLCAP, NNAT, SGK2, and GNAS was associated with drug sensitivity. While expression and methylation of imprinted genes in several other chromosomal regions was also associated with drug response and many imprinted genes in different chromosomal locations showed a considerable copy number variation, only imprinted genes at 20q11-q13.32 had a consistent association of their copy number with drug response. Copy number values among the imprinted genes in the 20q11-q13.32 region were strongly correlated. They were also correlated with the copy number of cancer-related non-imprinted genes MYBL2, AURKA, and ZNF217 in that chromosomal region. Expression of genes at 20q11-q13.32 was associated with ex vivo drug response in primary tumor samples from the Beat AML 1.0 acute myeloid leukemia patient cohort. Association of the increased copy number of the 20q11-q13.32 region with drug resistance may be complex and could involve multiple genes.

CONCLUSIONS

Copy number of imprinted and non-imprinted genes in the chromosomal region 20q11-q13.32 was associated with cancer drug resistance. The genes in this chromosomal region may have a modulating effect on tumor response to chemotherapy.

Further understanding of paternal uniparental disomy in Beckwith-Wiedemann syndrome

INTRODUCTION

Paternal uniparental disomy of chromosome 11 (upd(11)pat) accounts for up to 20% of molecularly confirmed Beckwith-Wiedemann spectrum (BWSp) cases. It belongs to the BWSp subgroup with the second highest tumor risk, and therefore needs particular awareness in research, diagnostics and clinical management.

AREAS COVERED

We overview the contribution of paternal (mosaic) uniparental disomy of chromosome 11 (UPD, upd(11)pat) and mosaic paternal uniparental diploidy in patients with Beckwith-Wiedemann features. The review comprises the current knowledge on their formation and their molecular and clinical consequences. Accordingly, the consequences for diagnostic testing and clinical monitoring are compiled.

EXPERT OPINION

The necessity to diagnostically identify and thus discriminate genome-wide paternal uniparental disomy, and upd(11)pat becomes obvious, due to the differences in the clinical course, disease prognosis, and treatment. In particular, monitoring of tumor development by liquid biopsy might be a promising option in the future. From the research point of view, it should be addressed why 11p is prone to mitotic recombination and thus also provide to the role of upd(11) as second hit in tumorigenesis.

Allele-specific differential regulation of monoallelically expressed autosomal genes in the cardiac lineage

Each mammalian autosomal gene is represented by two alleles in diploid cells. To our knowledge, no insights have been made in regard to allele-specific regulatory mechanisms of autosomes. Here we use allele-specific single cell transcriptomic analysis to elucidate the establishment of monoallelic gene expression in the cardiac lineage. We find that monoallelically expressed autosomal genes in mESCs and mouse blastocyst cells are differentially regulated based on the genetic background of the parental alleles. However, the genetic background of the allele does not affect the establishment of monoallelic genes in differentiated cardiomyocytes. Additionally, we observe epigenetic differences between deterministic and random autosomal monoallelic genes. Moreover, we also find a greater contribution of the maternal versus paternal allele to the development and homeostasis of cardiac tissue and in cardiac health, highlighting the importance of maternal influence in male cardiac tissue homeostasis. Our findings emphasize the significance of allele-specific insights into gene regulation in development, homeostasis and disease.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Genome-wide detection of imprinted differentially methylated regions using nanopore sequencing

Imprinting is a critical part of normal embryonic development in mammals, controlled by defined parent-of-origin (PofO) differentially methylated regions (DMRs) known as imprinting control regions. Direct nanopore sequencing of DNA provides a means to detect allelic methylation and to overcome the drawbacks of methylation array and short-read technologies. Here, we used publicly available nanopore sequencing data for 12 standard B-lymphocyte cell lines to acquire the genome-wide mapping of imprinted intervals in humans. Using the sequencing data, we were able to phase 95% of the human methylome and detect 94% of the previously well-characterized, imprinted DMRs. In addition, we found 42 novel imprinted DMRs (16 germline and 26 somatic), which were confirmed using whole-genome bisulfite sequencing (WGBS) data. Analysis of WGBS data in mouse (Mus musculus), rhesus monkey (Macaca mulatta), and chimpanzee (Pan troglodytes) suggested that 17 of these imprinted DMRs are conserved. Some of the novel imprinted intervals are within or close to imprinted genes without a known DMR. We also detected subtle parental methylation bias, spanning several kilobases at seven known imprinted clusters. At these blocks, hypermethylation occurs at the gene body of expressed allele(s) with mutually exclusive H3K36me3 and H3K27me3 allelic histone marks. These results expand upon our current knowledge of imprinting and the potential of nanopore sequencing to identify imprinting regions using only parent-offspring trios, as opposed to the large multi-generational pedigrees that have previously been required.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics

BACKGROUND

Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics.

METHODS

We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage.

RESULTS

We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions.

CONCLUSION

Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

Genome-wide DNA methylation patterns reveal clinically relevant predictive and prognostic subtypes in human osteosarcoma

Aberrant methylation of genomic DNA has been reported in many cancers. Specific DNA methylation patterns have been shown to provide clinically useful prognostic information and define molecular disease subtypes with different response to therapy and long-term outcome. Osteosarcoma is an aggressive malignancy for which approximately half of tumors recur following standard combined surgical resection and chemotherapy. No accepted prognostic factor save tumor necrosis in response to adjuvant therapy currently exists, and traditional genomic studies have thus far failed to identify meaningful clinical associations. We studied the genome-wide methylation state of primary tumors and tested how they predict patient outcomes. We discovered relative genomic hypomethylation to be strongly predictive of response to standard chemotherapy. Recurrence and survival were also associated with genomic methylation, but through more site-specific patterns. Furthermore, the methylation patterns were reproducible in three small independent clinical datasets. Downstream transcriptional, in vitro, and pharmacogenomic analysis provides insight into the clinical translation of the methylation patterns. Our findings suggest the assessment of genomic methylation may represent a strategy for stratifying patients for the application of alternative therapies.

Paternally expressed gene 3 (Pw1/Peg3) promotes sexual dimorphism in metabolism and behavior

The paternally expressed gene 3 (Pw1/Peg3) is a mammalian-specific parentally imprinted gene expressed in stem/progenitor cells of the brain and endocrine tissues. Here, we compared phenotypic characteristics in Pw1/Peg3 deficient male and female mice. Our findings indicate that Pw1/Peg3 is a key player for the determination of sexual dimorphism in metabolism and behavior. Mice carrying a paternally inherited Pw1/Peg3 mutant allele manifested postnatal deficits in GH/IGF dependent growth before weaning, sex steroid dependent masculinization during puberty, and insulin dependent fat accumulation in adulthood. As a result, Pw1/Peg3 deficient mice develop a sex-dependent global shift of body metabolism towards accelerated adiposity, diabetic-like insulin resistance, and fatty liver. Furthermore, Pw1/Peg3 deficient males displayed reduced social dominance and competitiveness concomitant with alterations in the vasopressinergic architecture in the brain. This study demonstrates that Pw1/Peg3 provides an epigenetic context that promotes male-specific characteristics through sex steroid pathways during postnatal development.

Pw1/Peg3 is under parental specific epigenetic regulation. We propose that Pw1/Peg3 confers a selective advantage in mammals by regulating sexual dimorphism. To address this question, we examined the consequences of Pw1/Peg3 loss of function in mice in an age- and sex-dependent context and found that Pw1/Peg3 mutants display reduced sexual dimorphism in growth, metabolism and behaviors. Our findings support the intralocus sexual conflict model of genomic imprinting where it contributes in sexual differentiation. Furthermore, our observations provide a unifying role of sex steroid signaling as a common property of Pw1/Peg3 expressing stem/progenitor cells and differentiated endocrine cells, both of which remain proliferative in response to gonadal hormones in adult life.

Functional organization of the maternal and paternal human 4D Nucleome

Every human somatic cell inherits a maternal and a paternal genome, which work together to give rise to cellular phenotypes. However, the allele-specific relationship between gene expression and genome structure through the cell cycle is largely unknown. By integrating haplotype-resolved genome-wide chromosome conformation capture, mature and nascent mRNA, and protein binding data from a B lymphoblastoid cell line, we investigate this relationship both globally and locally. We introduce the maternal and paternal 4D Nucleome, enabling detailed analysis of the mechanisms and dynamics of genome structure and gene function for diploid organisms. Our analyses find significant coordination between allelic expression biases and local genome conformation, and notably absent expression bias in universally essential cell cycle and glycolysis genes. We propose a model in which coordinated biallelic expression reflects prioritized preservation of essential gene sets.

Gaining insight into metabolic diseases from human genetic discoveries

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

Power calculator for detecting allelic imbalance using hierarchical Bayesian model

OBJECTIVE

Allelic imbalance (AI) is the differential expression of the two alleles in a diploid. AI can vary between tissues, treatments, and environments. Methods for testing AI exist, but methods are needed to estimate type I error and power for detecting AI and difference of AI between conditions. As the costs of the technology plummet, what is more important: reads or replicates?

RESULTS

We find that a minimum of 2400, 480, and 240 allele specific reads divided equally among 12, 5, and 3 replicates is needed to detect a 10, 20, and 30%, respectively, deviation from allelic balance in a condition with power > 80%. A minimum of 960 and 240 allele specific reads divided equally among 8 replicates is needed to detect a 20 or 30% difference in AI between conditions with comparable power. Higher numbers of replicates increase power more than adding coverage without affecting type I error. We provide a Python package that enables simulation of AI scenarios and enables individuals to estimate type I error and power in detecting AI and differences in AI between conditions.

Heterogeneity is not always a source of noise: Stochastic gene expression in regulatory heterozygote

Zygosity of diploid genome (i.e., degree to which two parental alleles of a gene have varied genetic sequences) adds another dimension to stochastic gene expression. The allelic imbalance in chromatin accessibility or divergence in regulatory sequences leads to fitness effects but the quantitative aspects thereof are largely left unexplored. We investigate diploid gene expression systems with homozygous (the same) and heterozygous (varied) combination of alleles in cis-regulatory sequences, not in structural gene loci, and characterize the zygosity-associated stochastic fluctuations in protein abundance. An emerging feat of heterozygosity is its counterintuitive capacity for genetic noise control. Especially when the noise is dominantly contributed to by the fluctuations in duty cycle ("reliability") rather than in process speed ("productivity") of gene expression machinery, its interallelic discrepancy acts to reduce the gene expression noise. These findings offer a novel insight into the rich repertoire of balancing selection enriched in the regulatory elements of immune response genes.

The Role of Long Non-coding RNAs in Human Imprinting Disorders: Prospective Therapeutic Targets

Genomic imprinting is a term used for an intergenerational epigenetic inheritance and involves a subset of genes expressed in a parent-of-origin-dependent way. Imprinted genes are expressed preferentially from either the paternally or maternally inherited allele. Long non-coding RNAs play essential roles in regulating this allele-specific expression. In several well-studied imprinting clusters, long non-coding RNAs have been found to be essential in regulating temporal- and spatial-specific establishment and maintenance of imprinting patterns. Furthermore, recent insights into the epigenetic pathological mechanisms underlying human genomic imprinting disorders suggest that allele-specific expressed imprinted long non-coding RNAs serve as an upstream regulator of the expression of other protein-coding or non-coding imprinted genes in the same cluster. Aberrantly expressed long non-coding RNAs result in bi-allelic expression or silencing of neighboring imprinted genes. Here, we review the emerging roles of long non-coding RNAs in regulating the expression of imprinted genes, especially in human imprinting disorders, and discuss three strategies targeting the central long non-coding RNA UBE3A-ATS for the purpose of developing therapies for the imprinting disorders Prader-Willi syndrome and Angelman syndrome. In summary, a better understanding of long non-coding RNA-related mechanisms is key to the development of potential therapeutic targets for human imprinting disorders.

Identification of Nephrogenic Therapeutic Biomarkers of Wilms Tumor Using Machine Learning

Wilms tumor is the most common renal malignancy in children, with a survival rate of more than 90%; however, treatment outcomes for certain patient subgroups, such as those with bilateral and recurrent diseases, remain significantly below this survival rate. Therefore, it remains essential to identify new biomarkers and develop effective therapeutic strategies. Based on the Therapeutically Applicable Research to Generate Effective Treatments and Gene Expression Omnibus RNA microarray datasets, we have identified eight differentially expressed genes in Wilms tumors as renal-specific in 33 randomly selected adult tumors. The risk model, constructed using survival forest and multivariate Cox regression, can effectively predict the prognosis; the risk score is an independent prognostic factor in Wilms tumor. Gene set enrichment analysis showed that most of the signature genes were involved in regulating human development-related pathways. At the same time, patients in the high-risk group exhibited more sensitive immunological and chemotherapeutic properties than those in the low-risk group. These results provide new insights into personalized and precise Wilms tumor treatment strategies.

Targeted RNA-seq improves efficiency, resolution, and accuracy of allele specific expression for human term placentas

Genomic imprinting is an epigenetic mechanism that results in allele-specific expression (ASE) based on the parent of origin. It is known to play a role in the prenatal and postnatal allocation of maternal resources in mammals. ASE detected by whole transcriptome RNA-seq (wht-RNAseq) has been widely used to analyze imprinted genes using reciprocal crosses in mice to generate large numbers of informative SNPs. Studies in humans are more challenging due to the paucity of SNPs and the poor preservation of RNA in term placentas and other tissues. Targeted RNA-seq (tar-RNAseq) can potentially mitigate these challenges by focusing sequencing resources on the regions of interest in the transcriptome. Here, we compared tar-RNAseq and wht-RNAseq in a study of ASE in known imprinted genes in placental tissue collected from a healthy human cohort in Mali, West Africa. As expected, tar-RNAseq substantially improved the coverage of SNPs. Compared to wht-RNAseq, tar-RNAseq produced on average four times more SNPs in twice as many genes per sample and read depth at the SNPs increased fourfold. In previous research on humans, discordant ASE values for SNPs of the same gene have limited the ability to accurately quantify ASE. We show that tar-RNAseq reduces this limitation as it unexpectedly increased the concordance of ASE between SNPs of the same gene, even in cases of degraded RNA. Studies aimed at discovering associations between individual variation in ASE and phenotypes in mammals and flowering plants will benefit from the improved power and accuracy of tar-RNAseq.

Conservation of Imprinting and Methylation of MKRN3, MAGEL2 and NDN Genes in Cattle

Genomic imprinting is the epigenetic mechanism of transcriptional regulation that involves differential DNA methylation modification. Comparative analysis of imprinted genes between species can help us to investigate the biological significance and regulatory mechanisms of genomic imprinting. MKRN3, MAGEL2 and NDN are three maternally imprinted genes identified in the human PWS/AS imprinted locus. This study aimed to assess the allelic expression of MKRN3, MAGEL2 and NDN and to examine the differentially methylated regions (DMRs) of bovine PWS/AS imprinted domains. An expressed single-nucleotide polymorphism (SNP)-based approach was used to investigate the allelic expression of MKRN3, MAGEL2 and NDN genes in bovine adult tissues and placenta. Consistent with the expression in humans and mice, we found that the MKRN3, MAGEL2 and NDN genes exhibit monoallelic expression in bovine somatic tissues and the paternal allele expressed in the bovine placenta. Three DMRs, PWS-IC, MKRN3 and NDN DMR, were identified in the bovine PWS/AS imprinted region by analysis of the DNA methylation status in bovine tissues using the bisulfite sequencing method and were located in the promoter and exon 1 of the SNRPN gene, NDN promoter and 5' untranslated region (5'UTR) of MKRN3 gene, respectively. The PWS-IC DMR is a primary DMR inherited from the male or female gamete, but NDN and MKRN3 DMR are secondary DMRs that occurred after fertilization by examining the methylation status in gametes.

New susceptible locus, rs9428555, is associated with pediatric-onset immunoglobulin A nephropathy and immunoglobulin A vasculitis in Koreans

BACKGROUND

Immunoglobulin A nephropathy (IgAN) is one of the most common primary forms of glomerulonephritis, while IgA vasculitis (IgAV) is the most common systemic vasculitis in children.

OBJECTIVE

Herein we aimed to uncover single nucleotide polymorphism (SNP) markers associated with these two related diseases by applying association tests and Sanger sequencing.

METHODS

Within the discovery stage, genomic DNA in blood samples from 101 enrolled patients were genotyped by the Korean Biobank Array. Association tests were performed with 397 Korean reference genomes. In the validation stage, 26 independent samples were genotyped by Sanger sequencing.

RESULTS

Four SNPs were identified (P < 5 × 10^-8) in the discovery stage. To determine whether the genotypes determined by SNP array were accurate, additional genotyping via Sanger sequencing was performed. As a result, only one SNP, rs9428555, was properly genotyped. In the validation stage, the minor allele (A > G) was found in as many as 15 out of 26 samples (minor allele frequency = 0.288), even though this minor allele is rare in East Asians (< 3%).

CONCLUSIONS

We found rs9428555 as a novel susceptible locus associated with the development of both IgAN and IgAV in Koreans. Though we cannot conclude rs9428555 is the unique susceptible locus of IgAN and IgAV, it is likely a good marker as the minor allele of this SNP occurred much more often in the patient group here versus within East Asians as a whole.

The RNA Atlas expands the catalog of human non-coding RNAs

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

Genome-wide detection of cytosine methylation by single molecule real-time sequencing

Single molecule real-time (SMRT) sequencing theoretically offers the opportunity to directly assess certain base modifications of native DNA molecules without any prior chemical/enzymatic conversions and PCR amplification, using kinetic signals of a DNA polymerase. However, the kinetic signal changes caused by 5mC modification are extremely subtle. Hence, the robust genome-wide measurement of 5mC modification has not been achieved. We enhanced 5mC detection using SMRT sequencing by holistically analyzing kinetic signals of a DNA polymerase and sequence context for every base within a measurement window. We employed a convolutional neural network to train a methylation classification model, leading to genome-wide 5mC detection. The sensitivity and specificity reached 90% and 94%, with a 99% correlation of overall methylation level with bisulfite sequencing.

5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human–mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Features and mechanisms of canonical and noncanonical genomic imprinting

Genomic imprinting is the monoallelic expression of a gene based on parent of origin and is a consequence of differential epigenetic marking between the male and female germlines. Canonically, genomic imprinting is mediated by allelic DNA methylation. However, recently it has been shown that maternal H3K27me3 can result in DNA methylation-independent imprinting, termed "noncanonical imprinting." In this review, we compare and contrast what is currently known about the underlying mechanisms, the role of endogenous retroviral elements, and the conservation of canonical and noncanonical genomic imprinting.

Dynamic Expression of Imprinted Genes in the Developing and Postnatal Pituitary Gland

In mammals, imprinted genes regulate many critical endocrine processes such as growth, the onset of puberty and maternal reproductive behaviour. Human imprinting disorders (IDs) are caused by genetic and epigenetic mechanisms that alter the expression dosage of imprinted genes. Due to improvements in diagnosis, increasing numbers of patients with IDs are now identified and monitored across their lifetimes. Seminal work has revealed that IDs have a strong endocrine component, yet the contribution of imprinted gene products in the development and function of the hypothalamo-pituitary axis are not well defined. Postnatal endocrine processes are dependent upon the production of hormones from the pituitary gland. While the actions of a few imprinted genes in pituitary development and function have been described, to date there has been no attempt to link the expression of these genes as a class to the formation and function of this essential organ. This is important because IDs show considerable overlap, and imprinted genes are known to define a transcriptional network related to organ growth. This knowledge deficit is partly due to technical difficulties in obtaining useful transcriptomic data from the pituitary gland, namely, its small size during development and cellular complexity in maturity. Here we utilise high-sensitivity RNA sequencing at the embryonic stages, and single-cell RNA sequencing data to describe the imprinted transcriptome of the pituitary gland. In concert, we provide a comprehensive literature review of the current knowledge of the role of imprinted genes in pituitary hormonal pathways and how these relate to IDs. We present new data that implicate imprinted gene networks in the development of the gland and in the stem cell compartment. Furthermore, we suggest novel roles for individual imprinted genes in the aetiology of IDs. Finally, we describe the dynamic regulation of imprinted genes in the pituitary gland of the pregnant mother, with implications for the regulation of maternal metabolic adaptations to pregnancy.

Inducible uniparental chromosome disomy to probe genomic imprinting at single-cell level in brain and beyond

Genomic imprinting is an epigenetic mechanism that results in parental allele-specific expression of ~1% of all genes in mouse and human. Imprinted genes are key developmental regulators and play pivotal roles in many biological processes such as nutrient transfer from the mother to offspring and neuronal development. Imprinted genes are also involved in human disease, including neurodevelopmental disorders, and often occur in clusters that are regulated by a common imprint control region (ICR). In extra-embryonic tissues ICRs can act over large distances, with the largest surrounding Igf2r spanning over 10 million base-pairs. Besides classical imprinted expression that shows near exclusive maternal or paternal expression, widespread biased imprinted expression has been identified mainly in brain. In this review we discuss recent developments mapping cell type specific imprinted expression in extra-embryonic tissues and neocortex in the mouse. We highlight the advantages of using an inducible uniparental chromosome disomy (UPD) system to generate cells carrying either two maternal or two paternal copies of a specific chromosome to analyze the functional consequences of genomic imprinting. Mosaic Analysis with Double Markers (MADM) allows fluorescent labeling and concomitant induction of UPD sparsely in specific cell types, and thus to over-express or suppress all imprinted genes on that chromosome. To illustrate the utility of this technique, we explain how MADM-induced UPD revealed new insights about the function of the well-studied Cdkn1c imprinted gene, and how MADM-induced UPDs led to identification of highly cell type specific phenotypes related to perturbed imprinted expression in the mouse neocortex. Finally, we give an outlook on how MADM could be used to probe cell type specific imprinted expression in other tissues in mouse, particularly in extra-embryonic tissues.

Insulin-Like Growth Factor 2 As a Possible Neuroprotective Agent and Memory Enhancer-Its Comparative Expression, Processing and Signaling in Mammalian CNS

A number of studies performed on rodents suggest that insulin-like growth factor 2 (IGF-2) or its analogs may possibly be used for treating some conditions like Alzheimer’s disease, Huntington’s disease, autistic spectrum disorders or aging-related cognitive impairment. Still, for translational research a comparative knowledge about the function of IGF-2 and related molecules in model organisms (rats and mice) and humans is necessary. There is a number of important differences in IGF-2 signaling between species. In the present review we emphasize species-specific patterns of IGF-2 expression in rodents, humans and some other mammals, using, among other sources, publicly available transcriptomic data. We provide a detailed description of Igf2 mRNA expression regulation and pre-pro-IGF-2 protein processing in different species. We also summarize the function of IGF-binding proteins. We describe three different receptors able to bind IGF-2 and discuss the role of IGF-2 signaling in learning and memory, as well as in neuroprotection. We hope that comprehensive understanding of similarities and differences in IGF-2 signaling between model organisms and humans will be useful for development of more effective medicines targeting IGF-2 receptors.

Imprinted Genes and Multiple Sclerosis: What Do We Know?

Multiple sclerosis (MS) is a chronic autoimmune neurodegenerative disease of the central nervous system that arises from interplay between non-genetic and genetic risk factors. The epigenetics functions as a link between these factors, affecting gene expression in response to external influence, and therefore should be extensively studied to improve the knowledge of MS molecular mechanisms. Among others, the epigenetic mechanisms underlie the establishment of parent-of-origin effects that appear as phenotypic differences depending on whether the allele was inherited from the mother or father. The most well described manifestation of parent-of-origin effects is genomic imprinting that causes monoallelic gene expression. It becomes more obvious that disturbances in imprinted genes at the least affecting their expression do occur in MS and may be involved in its pathogenesis. In this review we will focus on the potential role of imprinted genes in MS pathogenesis.

Clustering suicidal phenotypes and genetic associations with brain-derived neurotrophic factor in patients with substance use disorders

Suicide attempts (SA), especially recurrent SA or serious SA, are common in substance use disorders (SUD). However, the genetic component of SA in SUD samples remains unclear. Brain-derived neurotrophic factor (BDNF) alleles and levels have been repeatedly involved in stress-related psychopathology. This investigation uses a within-cases study of BDNF and associated factors in three suicidal phenotypes ('any', 'recurrent', and 'serious') of outpatients seeking treatment for opiate and/or cocaine use disorder. Phenotypic characterization was ascertained using a semi-structured interview. After thorough quality control, 98 SNPs of BDNF and associated factors (the BDNF pathway) were extracted from whole-genome data, leaving 411 patients of Caucasian ancestry, who had reliable data regarding their SA history. Binary and multinomial regression with the three suicidal phenotypes were further performed to adjust for possible confounders, along with hierarchical clustering and compared to controls (N = 2504). Bayesian analyses were conducted to detect pleiotropy across the suicidal phenotypes. Among 154 (37%) ever suicide attempters, 104 (68%) reported at least one serious SA and 96 (57%) two SA or more. The median number of non-tobacco SUDs was three. The BDNF gene remained associated with lifetime SA in SNP-based (rs7934165, rs10835210) and gene-based tests within the clinical sample. rs10835210 clustered with serious SA. Bayesian analysis identified genetic correlation between 'any' and 'serious' SA regarding rs7934165. Despite limitations, 'serious' SA was shown to share both clinical and genetic risk factors of SA-not otherwise specified, suggesting a shared BDNF-related pathophysiology of SA in this population with multiple SUDs.

DNA methylation abnormalities of imprinted genes in congenital heart disease: a pilot study

Background

Congenital heart disease (CHD) is resulted from the interaction of genetic aberration and environmental factors. Imprinted genes, which are regulated by epigenetic modifications, are essential for the normal embryonic development. However, the role of imprinted genes in the etiology of CHD remains unclear.

Methods

After the samples were treated with bisulfate salt, imprinted genes methylation were measured by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. T test and One-way ANOVA were performed to evaluate the differences among groups. Odds ratios (ORs) were performed to evaluate the incidence risk of CHD in relation to methylation levels.

Results

We investigated the alterations of imprinted gene germline differential methylation regions (gDMRs) methylation in patients with CHD. Eighteen imprinted genes that are known to affect early embryonic development were selected and the methylation modification genes were detected by massarray in 27 CHD children and 28 healthy children. Altered gDMR methylation level of 8 imprinted genes was found, including 2 imprinted genes with hypermethylation of GRB10 and MEST and 6 genes with hypomethylation of PEG10, NAP1L5, INPP5F, PLAGL1, NESP and MEG3. Stratified analysis showed that the methylation degree of imprinted genes was different in different types of CHD. Risk analysis showed that 6 imprinted genes, except MEST and NAP1L5, within a specific methylation level range were the risk factors for CHD

Conclusion

Altered methylation of imprinted genes is associated with CHD and varies in different types of CHD. Further experiments are warranted to identify the methylation characteristics of imprinted genes in different types of CHD and clarify the etiologies of imprinted genes in CHD.

Dynamics of paternal contributions to early embryo development in large animals

This review focuses on current knowledge of paternal contributions to preimplantation embryonic development with particular emphasis on large animals. Specifically, the included content aims to summarize genomic and epigenomic contributions of paternally expressed genes, their regulation, and chromatin structure that are indispensable for early embryo development. The accumulation of current knowledge will summarize conserved allelic function among species to include functional molecular and genomic studies across large domestic animals in context with reference to founding experimental models.

Cell-Type Specificity of Genomic Imprinting in Cerebral Cortex

In mammalian genomes, a subset of genes is regulated by genomic imprinting, resulting in silencing of one parental allele. Imprinting is essential for cerebral cortex development, but prevalence and functional impact in individual cells is unclear. Here, we determined allelic expression in cortical cell types and established a quantitative platform to interrogate imprinting in single cells. We created cells with uniparental chromosome disomy (UPD) containing two copies of either the maternal or the paternal chromosome; hence, imprinted genes will be 2-fold overexpressed or not expressed. By genetic labeling of UPD, we determined cellular phenotypes and transcriptional responses to deregulated imprinted gene expression at unprecedented single-cell resolution. We discovered an unexpected degree of cell-type specificity and a novel function of imprinting in the regulation of cortical astrocyte survival. More generally, our results suggest functional relevance of imprinted gene expression in glial astrocyte lineage and thus for generating cortical cell-type diversity.

A vast resource of allelic expression data spanning human tissues

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

Alterations in nicotinic receptor alpha5 subunit gene differentially impact early and later stages of cocaine addiction: a translational study in transgenic rats and patients

Cocaine addiction is a chronic and relapsing disorder with an important genetic component. Human candidate gene association studies showed that the single nucleotide polymorphism (SNP) rs16969968 in the α5 subunit (α5SNP) of nicotinic acetylcholine receptors (nAChRs), previously associated with increased tobacco dependence, was linked to a lower prevalence of cocaine use disorder (CUD). Three additional SNPs in the α5 subunit, previously shown to modify α5 mRNA levels, were also associated with CUD, suggesting an important role of the subunit in this pathology. To investigate the link between this subunit and CUD, we submitted rats knockout for the α5 subunit gene (α5KO), or carrying the α5SNP, to cocaine self-administration (SA) and showed that the acquisition of cocaine-SA was impaired in α5SNP rats while α5KO rats exhibited enhanced cocaine-induced relapse associated with altered neuronal activity in the nucleus accumbens. In addition, we observed in a human cohort of patients with CUD that the α5SNP was associated with a slower transition from first cocaine use to CUD. We also identified a novel SNP in the β4 nAChR subunit, part of the same gene cluster in the human genome and potentially altering CHRNA5 expression, associated with shorter time to relapse to cocaine use in patients. In conclusion, the α5SNP is protective against CUD by influencing early stages of cocaine exposure while CHRNA5 expression levels may represent a biomarker for the risk to relapse to cocaine use. Drugs modulating α5 containing nAChR activity may thus represent a novel therapeutic strategy against CUD.

Sex and gender: modifiers of health, disease, and medicine

Clinicians can encounter sex and gender disparities in diagnostic and therapeutic responses. These disparities are noted in epidemiology, pathophysiology, clinical manifestations, disease progression, and response to treatment. This Review discusses the fundamental influences of sex and gender as modifiers of the major causes of death and morbidity. We articulate how the genetic, epigenetic, and hormonal influences of biological sex influence physiology and disease, and how the social constructs of gender affect the behaviour of the community, clinicians, and patients in the health-care system and interact with pathobiology. We aim to guide clinicians and researchers to consider sex and gender in their approach to diagnosis, prevention, and treatment of diseases as a necessary and fundamental step towards precision medicine, which will benefit men's and women's health.

Genetic control of CCL24, POR, and IL23R contributes to the pathogenesis of sarcoidosis

Sarcoidosis is a genetically complex systemic inflammatory disease that affects multiple organs. We present a GWAS of a Japanese cohort (700 sarcoidosis cases and 886 controls) with replication in independent samples from Japan (931 cases and 1,042 controls) and the Czech Republic (265 cases and 264 controls). We identified three loci outside the HLA complex, CCL24, STYXL1-SRRM3, and C1orf141-IL23R, which showed genome-wide significant associations (P < 5.0 × 10⁻⁸) with sarcoidosis; CCL24 and STYXL1-SRRM3 were novel. The disease-risk alleles in CCL24 and IL23R were associated with reduced CCL24 and IL23R expression, respectively. The disease-risk allele in STYXL1-SRRM3 was associated with elevated POR expression. These results suggest that genetic control of CCL24, POR, and IL23R expression contribute to the pathogenesis of sarcoidosis. We speculate that the CCL24 risk allele might be involved in a polarized Th1 response in sarcoidosis, and that POR and IL23R risk alleles may lead to diminished host defense against sarcoidosis pathogens.

Akira Meguro et al. report a genome-wide association study for sarcoidosis—a systemic inflammatory disease—in the Japanese population. They identify 3 non-HLA loci with genome-wide significance, 2 of which have not been previously associated with sarcoidosis in any population.

Genetic susceptibility to nicotine addiction: Advances and shortcomings in our understanding of the CHRNA5/A3/B4 gene cluster contribution

Over the last decade, robust human genetic findings have been instrumental in elucidating the heritable basis of nicotine addiction (NA). They highlight coding and synonymous polymorphisms in a cluster on chromosome 15, encompassing the CHRNA5, CHRNA3 and CHRNB4 genes, coding for three subunits of the nicotinic acetylcholine receptor (nAChR). They have inspired an important number of preclinical studies, and will hopefully lead to the definition of novel drug targets for treating NA. Here, we review these candidate gene and genome-wide association studies (GWAS) and their direct implication in human brain function and NA-related phenotypes. We continue with a description of preclinical work in transgenic rodents that has led to a mechanistic understanding of several of the genetic hits. We also highlight important issues with regards to CHRNA3 and CHRNB4 where we are still lacking a dissection of their role in NA, including even in preclinical models. We further emphasize the use of human induced pluripotent stem cell-derived models for the analysis of synonymous and intronic variants on a human genomic background. Finally, we indicate potential avenues to further our understanding of the role of this human genetic variation. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Loss of imprinting of the human-specific imprinted gene ZNF597 causes prenatal growth retardation and dysmorphic features: implications for phenotypic overlap with Silver-Russell syndrome

BACKGROUND

ZNF597, encoding a zinc-finger protein, is the human-specific maternally expressed imprinted gene located on 16p13.3. The parent-of-origin expression of ZNF597 is regulated by the ZNF597:TSS-DMR, of which only the paternal allele acquires methylation during postimplantation period. Overexpression of ZNF597 may contribute to some of the phenotypes associated with maternal uniparental disomy of chromosome 16 (UPD(16)mat), and some patients with UPD(16)mat presenting with Silver-Russell syndrome (SRS) phenotype have recently been reported.

METHODS

A 6-year-old boy presented with prenatal growth restriction, macrocephaly at birth, forehead protrusion in infancy and clinodactyly of the fifth finger. Methylation, expression, microsatellite marker, single nucleotide polymorphism array and trio whole-exome sequencing analyses were conducted.

RESULTS

Isolated hypomethylation of the ZNF597:TSS-DMR and subsequent loss of imprinting and overexpression of ZNF597 were confirmed in the patient. Epigenetic alterations, such as UPD including UPD(16)mat and other methylation defects, were excluded. Pathogenic sequence or copy number variants affecting his phenotypes were not identified, indicating that primary epimutation occurred postzygotically.

CONCLUSION

We report the first case of isolated ZNF597 imprinting defect, showing phenotypic overlap with SRS despite not satisfying the clinical SRS criteria. A novel imprinting disorder entity involving the ZNF597 imprinted domain can be speculated.

The influence of DNA methylation on monoallelic expression

Monoallelic gene expression occurs in diploid cells when only one of the two alleles of a gene is active. There are three main classes of genes that display monoallelic expression in mammalian genomes: (1) imprinted genes that are monoallelically expressed in a parent-of-origin dependent manner; (2) X-linked genes that undergo random X-chromosome inactivation in female cells; (3) random monoallelically expressed single and clustered genes located on autosomes. The heritability of monoallelic expression patterns during cell divisions implies that epigenetic mechanisms are involved in the cellular memory of these expression states. Among these, methylation of CpG sites on DNA is one of the best described modification to explain somatic inheritance. Here, we discuss the relevance of DNA methylation for the establishment and maintenance of monoallelic expression patterns among these three groups of genes, and how this is intrinsically linked to development and cellular states.