Monoallelic Gene Expression in Mammals.

Monoallelic expression not due to cis-regulatory sequence polymorphism poses an intriguing problem in epigenetics because it requires the unequal treatment of two segments of DNA that are present in the same nucleus and that can indeed have absolutely identical sequences. Here, I focus on a few recent developments in the field of monoallelic expression that are of particular interest and raise interesting questions for future work. One development is regarding analyses of imprinted genes, in which recent work suggests the possibility that intriguing networks of imprinted genes exist and are important for genetic and physiological studies. Another issue that has been raised in recent years by a number of publications is the question of how skewed allelic expression should be for it to be designated as monoallelic expression and, further, what methods are appropriate or inappropriate for analyzing genomic data to examine allele-specific expression. Perhaps the most exciting recent development in mammalian monoallelic expression is a clever and carefully executed analysis of genetic diversity of autosomal genes subject to random monoallelic expression (RMAE), which provides compelling evidence for distinct evolutionary forces acting on random monoallelically expressed genes.

Mapping and characterizing N6-methyladenine in eukaryotic genomes using single-molecule real-time sequencing.

N6-Methyladenine (m6dA) has been discovered as a novel form of DNA methylation prevalent in eukaryotes; however, methods for high-resolution mapping of m6dA events are still lacking. Single-molecule real-time (SMRT) sequencing has enabled the detection of m6dA events at single-nucleotide resolution in prokaryotic genomes, but its application to detecting m6dA in eukaryotic genomes has not been rigorously examined. Herein, we identified unique characteristics of eukaryotic m6dA methylomes that fundamentally differ from those of prokaryotes. Based on these differences, we describe the first approach for mapping m6dA events using SMRT sequencing specifically designed for the study of eukaryotic genomes and provide appropriate strategies for designing experiments and carrying out sequencing in future studies. We apply the novel approach to study two eukaryotic genomes. For green algae, we construct the first complete genome-wide map of m6dA at single-nucleotide and single-molecule resolution. For human lymphoblastoid cells (hLCLs), it was necessary to integrate SMRT sequencing data with independent sequencing data. The joint analyses suggest putative m6dA events are enriched in the promoters of young full-length LINE-1 elements (L1s), but call for validation by additional methods. These analyses demonstrate a general method for rigorous mapping and characterization of m6dA events in eukaryotic genomes.

Mechanisms and consequences of widespread random monoallelic expression.

Although random monoallelic expression has been known for decades to affect genes on the X chromosome in female placental mammals, until a few years ago it was thought that there were few autosomal genes that were regulated in this manner. New tools for assaying gene expression genome-wide are now revealing that there are perhaps more genes that are subject to random monoallelic expression on mammalian autosomes than there are on the X chromosome and that these expression properties are achieved by diverse molecular mechanisms. This mode of expression has the potential to have an impact on natural selection and on the evolution of gene families.

Monoallelic expression of the human FOXP2 speech gene.

The recent descriptions of widespread random monoallelic expression (RMAE) of genes distributed throughout the autosomal genome indicate that there are more genes subject to RMAE on autosomes than the number of genes on the X chromosome where X-inactivation dictates RMAE of X-linked genes. Several of the autosomal genes that undergo RMAE have independently been implicated in human Mendelian disorders. Thus, parsing the relationship between allele-specific expression of these genes and disease is of interest. Mutations in the human forkhead box P2 gene, FOXP2, cause developmental verbal dyspraxia with profound speech and language deficits. Here, we show that the human FOXP2 gene undergoes RMAE. Studying an individual with developmental verbal dyspraxia, we identify a deletion 3 Mb away from the FOXP2 gene, which impacts FOXP2 gene expression in cis. Together these data suggest the intriguing possibility that RMAE impacts the haploinsufficiency phenotypes observed for FOXP2 mutations.

An architectural role for a nuclear noncoding RNA: NEAT1 RNA is essential for the structure of paraspeckles.

NEAT1 RNA, a highly abundant 4 kb ncRNA, is retained in nuclei in approximately 10 to 20 large foci that we show are completely coincident with paraspeckles, nuclear domains implicated in mRNA nuclear retention. Depletion of NEAT1 RNA via RNAi eradicates paraspeckles, suggesting that it controls sequestration of the paraspeckle proteins PSP1 and p54, factors linked to A-I editing. Unlike overexpression of PSP1, NEAT1 overexpression increases paraspeckle number, and paraspeckles emanate exclusively from the NEAT1 transcription site. The PSP-1 RNA binding domain is required for its colocalization with NEAT1 RNA in paraspeckles, and biochemical analyses support that NEAT1 RNA binds with paraspeckle proteins. Unlike other nuclear-retained RNAs, NEAT1 RNA is not A-I edited, consistent with a structural role in paraspeckles. Collectively, results demonstrate that NEAT1 functions as an essential structural determinant of paraspeckles, providing a precedent for a ncRNA as the foundation of a nuclear domain.

Modeling kinetic rate variation in third generation DNA sequencing data to detect putative modifications to DNA bases.

Current generation DNA sequencing instruments are moving closer to seamlessly sequencing genomes of entire populations as a routine part of scientific investigation. However, while significant inroads have been made identifying small nucleotide variation and structural variations in DNA that impact phenotypes of interest, progress has not been as dramatic regarding epigenetic changes and base-level damage to DNA, largely due to technological limitations in assaying all known and unknown types of modifications at genome scale. Recently, single-molecule real time (SMRT) sequencing has been reported to identify kinetic variation (KV) events that have been demonstrated to reflect epigenetic changes of every known type, providing a path forward for detecting base modifications as a routine part of sequencing. However, to date no statistical framework has been proposed to enhance the power to detect these events while also controlling for false-positive events. By modeling enzyme kinetics in the neighborhood of an arbitrary location in a genomic region of interest as a conditional random field, we provide a statistical framework for incorporating kinetic information at a test position of interest as well as at neighboring sites that help enhance the power to detect KV events. The performance of this and related models is explored, with the best-performing model applied to plasmid DNA isolated from Escherichia coli and mitochondrial DNA isolated from human brain tissue. We highlight widespread kinetic variation events, some of which strongly associate with known modification events, while others represent putative chemically modified sites of unknown types.

(Epi)Genetic analyses of age-related macular degeneration: case-control and discordant twin studies.

BACKGROUND/AIMS: Phenotypic discordance in monozygotic (MZ) twin pairs can have an epigenetic or genetic basis. Although age-related macular degeneration (AMD) has a strong genetic component, few studies have addressed its epigenetic basis. METHODS: Using SNP arrays, we evaluated differences in copy number variation (CNV) and allele-specific methylation (ASM) patterns (via methyl-sensitive restriction enzyme digestion of DNA) in MZ twin pairs from the US Twin Study of AMD. Further analyses examined the relationship between ASM and CNVs with AMD by both case/control analysis of ASM at candidate regions and by analysis of ASM and CNVs in twins discordant for AMD. RESULTS: The frequency of ASM sites differs between cases and controls in regions surrounding the AMD candidate genes CFH, C2 and CFB. While ASM patterns show a substantial dependence on local sequence polymorphisms, we observed dissimilar patterns of ASM between MZ twins. The genes closest to the sites where discordant MZ twins have dissimilar patterns of ASM are enriched for genes implicated in gliosis, a process associated with neovascular AMD. Similar twin-based analyses revealed no AMD-associated CNVs. CONCLUSIONS: Our results provide evidence of epigenetic influences beyond the known genetic susceptibility and implicate inflammatory responses and gliosis in the etiology of AMD.

Somatic mosaicism in schizophrenia brains reveals prenatal mutational processes.

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

Stochastic yet biased expression of multiple Dscam splice variants by individual cells.

The Drosophila melanogaster gene Dscam is essential for axon guidance and has 38,016 possible alternative splice forms. This diversity can potentially be used to distinguish cells. We analyzed the Dscam mRNA isoforms expressed by different cell types and individual cells. The choice of splice variants expressed is regulated both spatially and temporally. Different subtypes of photoreceptors express broad yet distinctive spectra of Dscam isoforms. Single-cell RT-PCR documented that individual cells express several different Dscam isoforms and allowed an estimation of the diversity that is present. For example, we estimate that each R3/R4 photoreceptor cell expresses 14-50 distinct mRNAs chosen from the spectrum of thousands of splice variants distinctive of its cell type. Thus, the Dscam repertoire of each cell is different from those of its neighbors, providing a potential mechanism for generating unique cell identity in the nervous system and elsewhere.

Coordination of the random asynchronous replication of autosomal loci.

Random monoallelic expression and asynchronous replication define an unusual class of autosomal mammalian genes. We show that every cell has randomly chosen either the maternal or paternal copy of each given autosome pair, such that alleles of these genes scattered across the chosen chromosome replicate earlier than the alleles on the homologous chromosome. Thus, chromosome-pair non-equivalence, rather than being limited to X-chromosome inactivation, is a fundamental property of mouse chromosomes.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Sex Differences in the Human Brain Transcriptome of Cases With Schizophrenia.

BACKGROUND: While schizophrenia differs between males and females in the age of onset, symptomatology, and disease course, the molecular mechanisms underlying these differences remain uncharacterized. METHODS: To address questions about the sex-specific effects of schizophrenia, we performed a large-scale transcriptome analysis of RNA sequencing data from 437 controls and 341 cases from two distinct cohorts from the CommonMind Consortium. RESULTS: Analysis across the cohorts identified a reproducible gene expression signature of schizophrenia that was highly concordant with previous work. Differential expression across sex was reproducible across cohorts and identified X- and Y-linked genes, as well as those involved in dosage compensation. Intriguingly, the sex expression signature was also enriched for genes involved in neurexin family protein binding and synaptic organization. Differential expression analysis testing a sex-by-diagnosis interaction effect did not identify any genome-wide signature after multiple testing corrections. Gene coexpression network analysis was performed to reduce dimensionality from thousands of genes to dozens of modules and elucidate interactions among genes. We found enrichment of coexpression modules for sex-by-diagnosis differential expression signatures, which were highly reproducible across the two cohorts and involved a number of diverse pathways, including neural nucleus development, neuron projection morphogenesis, and regulation of neural precursor cell proliferation. CONCLUSIONS: Overall, our results indicate that the effect size of sex differences in schizophrenia gene expression signatures is small and underscore the challenge of identifying robust sex-by-diagnosis signatures, which will require future analyses in larger cohorts.

EKLF/KLF1 expression defines a unique macrophage subset during mouse erythropoiesis.

Erythroblastic islands are a specialized niche that contain a central macrophage surrounded by erythroid cells at various stages of maturation. However, identifying the precise genetic and transcriptional control mechanisms in the island macrophage remains difficult due to macrophage heterogeneity. Using unbiased global sequencing and directed genetic approaches focused on early mammalian development, we find that fetal liver macrophages exhibit a unique expression signature that differentiates them from erythroid and adult macrophage cells. The importance of erythroid Krüppel-like factor (EKLF)/KLF1 in this identity is shown by expression analyses in EKLF-/- and in EKLF-marked macrophage cells. Single-cell sequence analysis simplifies heterogeneity and identifies clusters of genes important for EKLF-dependent macrophage function and novel cell surface biomarkers. Remarkably, this singular set of macrophage island cells appears transiently during embryogenesis. Together, these studies provide a detailed perspective on the importance of EKLF in the establishment of the dynamic gene expression network within erythroblastic islands in the developing embryo and provide the means for their efficient isolation.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Biallelic germline transcription at the kappa immunoglobulin locus.

Rearrangement of antigen receptor genes generates a vast array of antigen receptors on lymphocytes. The establishment of allelic exclusion in immunoglobulin genes requires differential treatment of the two sequence identical alleles. In the case of the kappa immunoglobulin locus, changes in chromatin structure, methylation, and replication timing of the two alleles are all potentially involved in regulating rearrangement. Additionally, germline transcription of the kappa locus which precedes rearrangement has been proposed to reflect an opening of the chromatin structure rendering it available for rearrangement. As the initial restriction of rearrangement to one allele is critical to the establishment of allelic exclusion, a key question is whether or not germline transcription at the kappa locus is monoallelic or biallelic. We have used a sensitive reverse transcription-polymerase chain reaction (RT-PCR) assay and an RNA-fluorescence in situ hybridization (FISH) to show that germline transcription of the kappa locus is biallelic in wild-type immature B cells and in recombination activating gene (RAG)-/-, mu+ B cells. Therefore, germline transcription is unlikely to dictate which allele will be rearranged first and rather reflects a general opening on both alleles that must be accompanied by a mechanism allowing one of the two alleles to be rearranged first.

Mice cloned from olfactory sensory neurons.

Cloning by nuclear transplantation has been successfully carried out in various mammals, including mice. Until now mice have not been cloned from post-mitotic cells such as neurons. Here, we have generated fertile mouse clones derived by transferring the nuclei of post-mitotic, olfactory sensory neurons into oocytes. These results indicate that the genome of a post-mitotic, terminally differentiated neuron can re-enter the cell cycle and be reprogrammed to a state of totipotency after nuclear transfer. Moreover, the pattern of odorant receptor gene expression and the organization of odorant receptor genes in cloned mice was indistinguishable from wild-type animals, indicating that irreversible changes to the DNA of olfactory neurons do not accompany receptor gene choice.

A novel hypomorphic MECP2 point mutation is associated with a neuropsychiatric phenotype.

The MECP2 gene on Xq28 encodes a transcriptional repressor, which binds to and modulates expression of active genes. Mutations in MECP2 cause classic or preserved speech variant Rett syndrome and intellectual disability in females and early demise or marked neurodevelopmental handicap in males. The consequences of a hypomorphic Mecp2 allele were recently investigated in a mouse model, which developed obesity, motor, social, learning, and behavioral deficits, predicting a human neurobehavioral syndrome. Here, we describe mutation analysis of a nondysmorphic female proband and her father who presented with primarily neuropsychiatric manifestations and obesity with relative sparing of intelligence, language, growth, and gross motor skills. We identified and characterized a novel missense mutation (c.454C>G; p.P152A) in the critical methyl-binding domain of MeCP2 that disrupts MeCP2 functional activity. We show that a gradient of impairment is present when the p.P152A mutation is compared with an allelic p.P152R mutation, which causes classic Rett syndrome and another Rett syndrome-causing mutation, such that protein-heterochromatin binding observed by immunofluorescence and immunoblotting is wild-type > P152A > P152R > T158 M, consistent with the severity of the observed phenotype. Our findings provide evidence for very mild phenotypes in humans associated with partial reduction of MeCP2 function arising from subtle variation in MECP2.

The variable transmembrane domain of Drosophila N-cadherin regulates adhesive activity.

Drosophila N-cadherin (CadN) is an evolutionarily conserved classic cadherin which has a large, complex extracellular domain and a catenin-binding cytoplasmic domain. The CadN locus contains three modules of alternative exons (7a/b, 13a/b, and 18a/b) and undergoes alternative splicing to generate multiple isoforms. Using quantitative transcript analyses and green fluorescent protein-based cell sorting, we found that during development CadN alternative splicing is regulated in a temporal but not cell-type-specific fashion. In particular, exon 18b is predominantly expressed during early developmental stages, while exon 18a is prevalent at the late developmental and adult stages. All CadN isoforms share the same molecular architecture but have different sequences in their extracellular and transmembrane domains, suggesting functional diversity. In vitro quantitative cell aggregation assays revealed that all CadN isoforms mediate homophilic interactions, but the isoforms encoded by exon 18b have a higher adhesive activity than those by its alternative, 18a. Domain-swapping experiments further revealed that the different sequences in the transmembrane domains of isoforms are responsible for their differential adhesive activities. CadN alternative splicing might provide a novel mechanism to fine-tune its adhesive activity at different developmental stages or to restrict the use of high-affinity 18b-type isoforms at the adult stage.