Transcriptome-wide survey of mouse CNS-derived cells reveals monoallelic expression within novel gene families

The counterpart congenital overgrowth syndromes Beckwith-Wiedemann Syndrome in human and large offspring syndrome in bovine involve alterations in DNA methylation, transcription, and chromatin configuration

Beckwith-Wiedemann Syndrome (BWS, OMIM #130650) is a congenital epigenetic disorder in humans which affects approximately 1 in 10,340 children. The incidence is likely an underestimation as the condition is usually recognized based on observable phenotypes at birth. BWS children have up to a 28% risk of developing tumors and currently, only 80% of patients can be corroborated molecularly (epimutations/variants). It is unknown how the subtypes of this condition are molecularly similar/dissimilar globally, therefore there is a need to deeply characterize the syndrome at the molecular level. Here we characterize the methylome, transcriptome and chromatin configuration of 18 BWS individuals together with the animal model of the condition, the bovine large offspring syndrome (LOS). Sex specific comparisons are performed for a subset of the BWS patients and LOS. Given that this epigenetic overgrowth syndrome has been characterized as a loss-of-imprinting condition, parental allele-specific comparisons were performed using the bovine animal model. In general, the differentially methylated regions (DMRs) detected in BWS and LOS showed significant enrichment for CTCF binding sites. Altered chromosome compartments in BWS and LOS were positively correlated with gene expression changes, and the promoters of differentially expressed genes showed significant enrichment for DMRs, differential topologically associating domains, and differential A/B compartments in some comparisons of BWS subtypes and LOS. We show shared regions of dysregulation between BWS and LOS, including several HOX gene clusters, and also demonstrate that altered DNA methylation differs between the clinically epigenetically identified BWS patients and those identified as having DNA variants (i.e. CDKN1C microdeletion). Lastly, we highlight additional genes and genomic regions that have the potential to serve as targets for biomarker development to improve current molecular methodologies. In summary, our results suggest that genome-wide alternation of chromosome architecture, which is partially caused by DNA methylation changes, also contribute to the development of BWS and LOS.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

In Vivo Clonal Analysis Reveals Random Monoallelic Expression in Lymphocytes That Traces Back to Hematopoietic Stem Cells

Evaluating the epigenetic landscape in the stem cell compartment at the single-cell level is essential to assess the cells’ heterogeneity and predict their fate. Here, using a genome-wide transcriptomics approach in vivo, we evaluated the allelic expression imbalance in the progeny of single hematopoietic cells (HSCs) as a read-out of epigenetic marking. After 4 months of extensive proliferation and differentiation, we found that X-chromosome inactivation (XCI) is tightly maintained in all single-HSC derived hematopoietic cells. In contrast, the vast majority of the autosomal genes did not show clonal patterns of random monoallelic expression (RME). However, a persistent allele-specific autosomal transcription in HSCs and their progeny was found in a rare number of cases, none of which has been previously reported. These data show that: 1) XCI and RME in the autosomal chromosomes are driven by different mechanisms; 2) the previously reported high frequency of genes under RME in clones expanded in vitro (up to 15%) is not found in clones undergoing multiple differentiation steps in vivo; 3) prior to differentiation, HSCs have stable patterns of autosomal RME. We propose that most RME patterns in autosomal chromosomes are erased and established de novo during cell lineage differentiation.

RNA sequencing-based screen for reactivation of silenced alleles of autosomal genes

In mammalian cells, maternal and paternal alleles usually have similar transcriptional activity. Epigenetic mechanisms such as X-chromosome inactivation (XCI) and imprinting were historically viewed as rare exceptions to this rule. Discovery of autosomal monoallelic autosomal expression (MAE) a decade ago revealed an additional allele-specific mode regulating thousands of mammalian genes. Despite MAE prevalence, its mechanistic basis remains unknown. Using an RNA sequencing-based screen for reactivation of silenced alleles, we identified DNA methylation as key mechanism of MAE mitotic maintenance. In contrast with the all-or-nothing allelic choice in XCI, allele-specific expression in MAE loci is tunable, with exact allelic imbalance dependent on the extent of DNA methylation. In a subset of MAE genes, allelic imbalance was insensitive to DNA demethylation, implicating additional mechanisms in MAE maintenance in these loci. Our findings identify a key mechanism of MAE maintenance and provide basis for understanding the biological role of MAE.

Introns control stochastic allele expression bias

Monoallelic expression (MAE) or extreme allele bias can account for incomplete penetrance, missing heritability and non-Mendelian diseases. In cancer, MAE is associated with shorter patient survival times and higher tumor grade. Prior studies showed that stochastic MAE is caused by stochastic epigenetic silencing, in a gene and tissue-specific manner. Here, we used C. elegans to study stochastic MAE in vivo. We found allele bias/MAE to be widespread within C. elegans tissues, presenting as a continuum from fully biallelic to MAE. We discovered that the presence of introns within alleles robustly decreases MAE. We determined that introns control MAE at distinct loci, in distinct cell types, with distinct promoters, and within distinct coding sequences, using a 5'-intron position-dependent mechanism. Bioinformatic analysis showed human intronless genes are significantly enriched for MAE. Our experimental evidence demonstrates a role for introns in regulating MAE, possibly explaining why some mutations within introns result in disease.

Differential Allelic Expression among Long Non-Coding RNAs

Long non-coding RNAs (lncRNA) comprise a diverse group of non-protein-coding RNAs >200 bp in length that are involved in various normal cellular processes and disease states, and can affect coding gene expression through mechanisms in cis or in trans. Since the discovery of the first functional lncRNAs transcribed by RNA Polymerase II, H19 and Xist, many others have been identified and noted for their unusual transcriptional pattern, whereby expression from one chromosome homolog is strongly favored over the other, also known as mono-allelic or differential allelic expression. lncRNAs with differential allelic expression have been observed to play critical roles in developmental gene regulation, chromosome structure, and disease. Here, we will focus on known examples of differential allelic expression of lncRNAs and highlight recent research describing functional lncRNAs expressed from both imprinted and random mono-allelic expression domains.

X-Chromosome Inactivation and Autosomal Random Monoallelic Expression as "Faux Amis"

X-chromosome inactivation (XCI) and random monoallelic expression of autosomal genes (RMAE) are two paradigms of gene expression regulation where, at the single cell level, genes can be expressed from either the maternal or paternal alleles. X-chromosome inactivation takes place in female marsupial and placental mammals, while RMAE has been described in mammals and also other species. Although the outcome of both processes results in random monoallelic expression and mosaicism at the cellular level, there are many important differences. We provide here a brief sketch of the history behind the discovery of XCI and RMAE. Moreover, we review some of the distinctive features of these two phenomena, with respect to when in development they are established, their roles in dosage compensation and cellular phenotypic diversity, and the molecular mechanisms underlying their initiation and stability.

Locus specific epigenetic modalities of random allelic expression imbalance

Most autosomal genes are thought to be expressed from both alleles, with some notable exceptions, including imprinted genes and genes showing random monoallelic expression (RME). The extent and nature of RME has been the subject of debate. Here we investigate the expression of several candidate RME genes in F1 hybrid mouse cells before and after differentiation, to define how they become persistently, monoallelically expressed. Clonal monoallelic expression is not present in embryonic stem cells, but we observe high frequencies of monoallelism in neuronal progenitor cells by assessing expression status in more than 200 clones. We uncover unforeseen modes of allelic expression that appear to be gene-specific and epigenetically regulated. This non-canonical allelic regulation has important implications for development and disease, including autosomal dominant disorders and opens up therapeutic perspectives.

Parental bias in expression and interaction of genes in the equine placenta

Most autosomal genes in the placenta show a biallelic expression pattern. However, some genes exhibit allele-specific transcription depending on the parental origin of the chromosomes on which the copy of the gene resides. Parentally expressed genes are involved in the reciprocal interaction between maternal and paternal genes, coordinating the allocation of resources between fetus and mother. One of the main challenges of studying parental-specific allelic expression (allele-specific expression [ASE]) in the placenta is the maternal cellular remnant at the fetomaternal interface. Horses (Equus caballus) have an epitheliochorial placenta in which both the endometrial epithelium and the epithelium of the chorionic villi are juxtaposed with minimal extension into the uterine mucosa, yet there is no information available on the allelic gene expression of equine chorioallantois (CA). In the current study, we present a dataset of 1,336 genes showing ASE in the equine CA (https://pouya-dini.github.io/equine-gene-db/) along with a workflow for analyzing ASE genes. We further identified 254 potentially imprinted genes among the parentally expressed genes in the equine CA and evaluated the expression pattern of these genes throughout gestation. Our gene ontology analysis implies that maternally expressed genes tend to decrease the length of gestation, while paternally expressed genes extend the length of gestation. This study provides fundamental information regarding parental gene expression during equine pregnancy, a species with a negligible amount of maternal cellular remnant in its placenta. This information will provide the basis for a better understanding of the role of parental gene expression in the placenta during gestation.

Reciprocal monoallelic expression of ASAR lncRNA genes controls replication timing of human chromosome 6

DNA replication occurs on mammalian chromosomes in a cell-type distinctive temporal order known as the replication timing program. We previously found that disruption of the noncanonical lncRNA genes ASAR6 and ASAR15 results in delayed replication timing and delayed mitotic chromosome condensation of human chromosomes 6 and 15, respectively. ASAR6 and ASAR15 display random monoallelic expression and display asynchronous replication between alleles that is coordinated with other random monoallelic genes on their respective chromosomes. Disruption of the expressed allele, but not the silent allele, of ASAR6 leads to delayed replication, activation of the previously silent alleles of linked monoallelic genes, and structural instability of human chromosome 6. In this report, we describe a second lncRNA gene (ASAR6-141) on human chromosome 6 that when disrupted results in delayed replication timing in cisASAR6-141 is subject to random monoallelic expression and asynchronous replication and is expressed from the opposite chromosome 6 homolog as ASAR6 ASAR6-141 RNA, like ASAR6 and ASAR15 RNAs, contains a high L1 content and remains associated with the chromosome territory where it is transcribed. Three classes of cis-acting elements control proper chromosome function in mammals: origins of replication, centromeres, and telomeres, which are responsible for replication, segregation, and stability of all chromosomes. Our work supports a fourth type of essential chromosomal element, the "Inactivation/Stability Center," which expresses ASAR lncRNAs responsible for proper replication timing, monoallelic expression, and structural stability of each chromosome.

Role of the Chromosome Architectural Factor SMCHD1 in X-Chromosome Inactivation, Gene Regulation, and Disease in Humans

Structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1) is an architectural factor critical for X-chromosome inactivation (XCI) and the repression of select autosomal gene clusters. In mice, homozygous nonsense mutations in Smchd1 cause female-specific embryonic lethality due to an XCI defect. However, although human mutations in SMCHD1 are associated with congenital arhinia and facioscapulohumeral muscular dystrophy type 2 (FSHD2), the diseases do not show a sex-specific bias, despite the essential nature of XCI in humans. To investigate whether there is a dosage imbalance for the sex chromosomes, we here analyze transcriptomic data from arhinia and FSHD2 patient blood and muscle cells. We find that X-linked dosage compensation is maintained in these patients. In mice, SMCHD1 controls not only protocadherin (Pcdh) gene clusters, but also Hox genes critical for craniofacial development. Ablating Smchd1 results in aberrant expression of these genes, coinciding with altered chromatin states and three-dimensional (3D) topological organization. In a subset of FSHD2 and arhinia patients, we also found dysregulation of clustered PCDH, but not HOX genes. Overall, our study demonstrates preservation of XCI in arhinia and FSHD2, and implicates SMCHD1 in the regulation of the 3D organization of select autosomal gene clusters.

Allele-specific RNA imaging shows that allelic imbalances can arise in tissues through transcriptional bursting

Extensive cell-to-cell variation exists even among putatively identical cells, and there is great interest in understanding how the properties of transcription relate to this heterogeneity. Differential expression from the two gene copies in diploid cells could potentially contribute, yet our ability to measure from which gene copy individual RNAs originated remains limited, particularly in the context of tissues. Here, we demonstrate quantitative, single molecule allele-specific RNA FISH adapted for use on tissue sections, allowing us to determine the chromosome of origin of individual RNA molecules in formaldehyde-fixed tissues. We used this method to visualize the allele-specific expression of Xist and multiple autosomal genes in mouse kidney. By combining these data with mathematical modeling, we evaluated models for allele-specific heterogeneity, in particular demonstrating that apparent expression from only one of the alleles in single cells can arise as a consequence of low-level mRNA abundance and transcriptional bursting.

In mammals, most cells of the body contain two genetic datasets: one from the mother and one from the father, and—in theory—these two sets of information could contribute equally to produce the molecules in a given cell. In practice, however, this is not always the case, which can have dramatic implications for many traits, including visible features (such as fur color) and even disease outcomes. However, it remains difficult to measure the parental origin of individual molecules in a given cell and thus to assess what processes lead to an imbalance of the maternal and paternal contribution. We adapted a microscopy technique—called allele-specific single molecule RNA FISH—that uses a combination of fluorescent tags to specifically label one type of molecule, RNA, depending on its origin, for use in mouse kidney sections. Focusing on RNAs that were previously reported to show imbalance, we performed measurements and combined these with mathematical modeling to quantify imbalance in tissues and explain how these can arise. We found that we could recapitulate the observed imbalances using models that only take into account the random processes that produce RNA, without the need to invoke special regulatory mechanisms to create unequal contributions.

Frequent monoallelic or skewed expression for developmental genes in CNS-derived cells and evidence for balancing selection

Cellular mosaicism due to monoallelic autosomal expression (MAE), with cell selection during development, is becoming increasingly recognized as prevalent in mammals, leading to interest in understanding its extent and mechanism(s). We report here use of clonal cell lines derived from the CNS of adult female [Formula: see text] hybrid (C57BL/6 X JF1) mice to characterize MAE as neural stem cells (nscs) differentiate to astrocyte-like cells (asls). We found that different subsets of genes show MAE in the two populations of cells; in each case, there is strong enrichment for genes specific to the respective developmental state. Genes that exhibit MAE are 22% of nsc-specific genes and 26% of asl-specific genes. Moreover, the promoters of genes with MAE have reduced CpG dinucleotides but increased CpG differences between the two parental mouse strains. Extending the study of variability to wild populations of mice, we found evidence for balancing selection as a contributing force in evolution of those genes showing developmental specificity (i.e., expressed in either nsc or asl), not just for genes showing MAE. Furthermore, we found that genes showing skewed allelic expression (SKE) were similarly enriched among cell type-specific genes and also showed a heightened probability of balancing selection. Thus, developmental stage-specific genes and genes with MAE or SKE seem to make up overlapping classes subject to selection for increased diversity. The implications of these results for development and evolution are discussed in the context of a model with stochastic epigenetic modifications taking place only during a relatively brief developmental window.

Dosage compensation and DNA methylation landscape of the X chromosome in mouse liver

DNA methylation plays a key role in X-chromosome inactivation (XCI), a process that achieves dosage compensation for X-encoded gene products between mammalian female and male cells. However, differential sex chromosome dosage complicates genome-wide epigenomic assessments, and the X chromosome is frequently excluded from female-to-male comparative analyses. Using the X chromosome in the sexually dimorphic mouse liver as a model, we provide a general framework for comparing base-resolution DNA methylation patterns across samples that have different chromosome numbers and ask at a systematic level if predictions by historical analyses of X-linked DNA methylation hold true at a base-resolution chromosome-wide level. We demonstrate that sex-specific methylation patterns on the X chromosome largely reflect the effects of XCI. While our observations concur with longstanding observations of XCI at promoter-proximal CpG islands, we provide evidence that sex-specific DNA methylation differences are not limited to CpG island boundaries. Moreover, these data support a model in which maintenance of CpG islands in the inactive state does not require complete regional methylation. Further, we validate an intragenic non-CpG methylation signature in genes escaping XCI in mouse liver. Our analyses provide insight into underlying methylation patterns that should be considered when assessing sex differences in genome-wide methylation analyses.

Allele-Specific Biased Expression of the CNTN6 Gene in iPS Cell-Derived Neurons from a Patient with Intellectual Disability and 3p26.3 Microduplication Involving the CNTN6 Gene

Copy number variations (CNVs) of the human CNTN6 gene caused by megabase-scale microdeletions or microduplications in the 3p26.3 region are often the cause of neurodevelopmental disorders, including intellectual disability and developmental delay. Surprisingly, patients with different copy numbers of this gene display notable overlapping of neuropsychiatric symptoms. The complexity of the study of human neuropathologies is associated with the inaccessibility of brain material. This problem can be overcome through the use of reprogramming technologies that permit the generation of induced pluripotent stem (iPS) cells from fibroblasts and their subsequent in vitro differentiation into neurons. We obtained a set of iPS cell lines derived from a patient carrier of the CNTN6 gene duplication and from two healthy donors. All iPS cell lines displayed the characteristics of pluripotent cells. Some iPS cell lines derived from the patient and from healthy donors were differentiated in vitro by exogenous expression of the Ngn2 transcription factor or by spontaneous neural differentiation of iPS cells through the neural rosette stage. The obtained neurons showed the characteristics of mature neurons as judged by the presence of neuronal markers and by their electrophysiological characteristics. Analysis of allele-specific expression of the CNTN6 gene in these neuronal cells by droplet digital PCR demonstrated that the level of expression of the duplicated allele was significantly reduced compared to that of the wild-type allele. Importantly, according to the sequencing data, both copies of the CNTN6 gene, which were approximately 1 Mb in size, showed no any additional structural rearrangements.

L1 retrotransposon antisense RNA within ASAR lncRNAs controls chromosome-wide replication timing

Proper chromosome duplication is critical for genome integrity and normal cellular function. Platt et al. show that the lncRNA genes ASAR6 and ASAR15 control chromosome-wide replication timing via the antisense strand of L1 retrotransposons located within ASAR6 and ASAR15 RNAs.

Mammalian cells replicate their chromosomes via a temporal replication program. The ASAR6 and ASAR15 genes were identified as loci that when disrupted result in delayed replication and condensation of entire human chromosomes. ASAR6 and ASAR15 are monoallelically expressed long noncoding RNAs that remain associated with the chromosome from which they are transcribed. The chromosome-wide effects of ASAR6 map to the antisense strand of an L1 retrotransposon within ASAR6 RNA, deletion or inversion of which delayed replication of human chromosome 6. Furthermore, ectopic integration of ASAR6 or ASAR15 transgenes into mouse chromosomes resulted in delayed replication and condensation, an increase in H3K27me3, coating of the mouse chromosome with ASAR RNA, and a loss of mouse Cot-1 RNA expression in cis. Targeting the antisense strand of the L1 within ectopically expressed ASAR6 RNA restored normal replication timing. Our results provide direct evidence that L1 antisense RNA plays a functional role in chromosome-wide replication timing of mammalian chromosomes.

CD-tagging-MS2: detecting allelic expression of endogenous mRNAs and their protein products in single cells

Discriminating between the mRNA and protein outputs of each of the alleles of an endogenous gene in intact cells, is a difficult task. To examine endogenous transcripts originating from a specific allele, we applied Central Dogma tagging (CD-tagging), which is based on a tag insertion into an endogenous gene by creation of a new exon. Previously, CD-tagging was used to tag endogenous proteins. Here we developed a CD-tagging-MS2 approach in which two tags were inserted in tandem; a fluorescent protein tag in conjunction with the mRNA MS2 tag used for tagging mRNAs in cells. A cell clone library of CD-tagged-MS2 genes was generated, and protein and mRNA distributions were examined and characterized in single cells. Taking advantage of having one allele tagged, we demonstrate how the transcriptional activity of all alleles, tagged and untagged, can be identified using single molecule RNA fluorescence in situ hybridization (smFISH). Allele-specific mRNA expression and localization were quantified under normal and stress conditions. The latter generate cytoplasmic stress granules (SGs) that can store mRNAs, and the distribution of the mRNAs within and outside of the SGs was measured. Altogether, CD-tagging-MS2 is a robust and inexpensive approach for direct simultaneous detection of an endogenous mRNA and its translated protein product in the same cell.

The X chromosome in space

Extensive 3D folding is required to package a genome into the tiny nuclear space, and this packaging must be compatible with proper gene expression. Thus, in the well-hierarchized nucleus, chromosomes occupy discrete territories and adopt specific 3D organizational structures that facilitate interactions between regulatory elements for gene expression. The mammalian X chromosome exemplifies this structure-function relationship. Recent studies have shown that, upon X-chromosome inactivation, active and inactive X chromosomes localize to different subnuclear positions and adopt distinct chromosomal architectures that reflect their activity states. Here, we review the roles of long non-coding RNAs, chromosomal organizational structures and the subnuclear localization of chromosomes as they relate to X-linked gene expression.

Erasure and reestablishment of random allelic expression imbalance after epigenetic reprogramming

Clonal level random allelic expression imbalance and random monoallelic expression provides cellular heterogeneity within tissues by modulating allelic dosage. Although such expression patterns have been observed in multiple cell types, little is known about when in development these stochastic allelic choices are made. We examine allelic expression patterns in human neural progenitor cells before and after epigenetic reprogramming to induced pluripotency, observing that loci previously characterized by random allelic expression imbalance (0.63% of expressed genes) are generally reset to a biallelic state in induced pluripotent stem cells (iPSCs). We subsequently neuralized the iPSCs and profiled isolated clonal neural stem cells, observing that significant random allelic expression imbalance is reestablished at 0.65% of expressed genes, including novel loci not found to show allelic expression imbalance in the original parental neural progenitor cells. Allelic expression imbalance was associated with altered DNA methylation across promoter regulatory regions, with clones characterized by skewed allelic expression being hypermethylated compared to their biallelic sister clones. Our results suggest that random allelic expression imbalance is established during lineage commitment and is associated with increased DNA methylation at the gene promoter.

Analysis of allelic expression patterns in clonal somatic cells by single-cell RNA-seq

Cellular heterogeneity can emerge from the expression of only one parental allele. However, it has remained controversial whether, or to what degree, random monoallelic expression of autosomal genes (aRME) is mitotically inherited (clonal) or stochastic (dynamic) in somatic cells, particularly in vivo. Here we used allele-sensitive single-cell RNA-seq on clonal primary mouse fibroblasts and freshly isolated human CD8⁺ T cells to dissect clonal and dynamic monoallelic expression patterns. Dynamic aRME affected a considerable portion of the cells' transcriptomes, with levels dependent on the cells' transcriptional activity. Notably, clonal aRME was detected, but it was surprisingly scarce (<1% of genes) and mainly affected the most weakly expressed genes. Consequently, the overwhelming majority of aRME occurs transiently within individual cells, and patterns of aRME are thus primarily scattered throughout somatic cell populations rather than, as previously hypothesized, confined to patches of clonally related cells.

Random monoallelic expression of genes on autosomes: Parallels with X-chromosome inactivation

Genes are generally expressed from their two alleles, except in some particular cases such as random inactivation of one of the two X chromosomes in female mammals or imprinted genes which are expressed only from the maternal or the paternal allele. A lesser-known phenomenon is random monoallelic expression (RME) of autosomal genes, where genes can be stably expressed in a monoallelic manner, from either one of the parental alleles. Studies on autosomal RME face several challenges. First, RME that is based on epigenetic mechanisms has to be distinguished from biased expression of one allele caused by a DNA sequence polymorphism in a regulatory element. Second, RME should not be confused with transient monoallelic expression often observed in single cell analyses, and that often corresponds to dynamic bursting of expression. Thanks to analyses on clonal cell populations, the existence of RME in cultured cells is now well established. Future studies of RME in vivo will have to overcome tissue heterogeneity and certain technical limitations. Here, we discuss current knowledge on autosomal RME, as well as possible mechanisms controlling these expression patterns and potential implications for development and disease, drawing parallels with what is known for X-chromosome inactivation, a paradigm of random monoallelic expression.

Genes with monoallelic expression contribute disproportionately to genetic diversity in humans

An unexpectedly large number of human autosomal genes are subject to monoallelic expression (MAE). Our analysis of 4,227 such genes uncovers surprisingly high genetic variation across human populations. This increased diversity is unlikely to reflect relaxed purifying selection. Remarkably, MAE genes exhibit an elevated recombination rate and an increased density of hypermutable sequence contexts. However, these factors do not fully account for the increased diversity. We find that the elevated nucleotide diversity of MAE genes is also associated with greater allelic age: variants in these genes tend to be older and are enriched in polymorphisms shared by Neanderthals and chimpanzees. Both synonymous and nonsynonymous alleles of MAE genes have elevated average population frequencies. We also observed strong enrichment of the MAE signature among genes reported to evolve under balancing selection. We propose that an important biological function of widespread MAE might be the generation of cell-to-cell heterogeneity; the increased genetic variation contributes to this heterogeneity.

Random monoallelic expression of autosomal genes: stochastic transcription and allele-level regulation

Random monoallelic expression (RME) of genes represents a striking example of how stochastic molecular processes can result in cellular heterogeneity. Recent transcriptome-wide studies have revealed both mitotically stable and cell-to-cell dynamic forms of autosomal RME, with the latter presumably resulting from burst-like stochastic transcription. Here, we discuss the distinguishing features of these two forms of RME and revisit literature on their nature, pervasiveness and regulation. Finally, we explore how RME may contribute to phenotypic variation, including the incomplete penetrance and variable expressivity often seen in genetic disease.

Dynamics of gene silencing during X inactivation using allele-specific RNA-seq

BACKGROUND

During early embryonic development, one of the two X chromosomes in mammalian female cells is inactivated to compensate for a potential imbalance in transcript levels with male cells, which contain a single X chromosome. Here, we use mouse female embryonic stem cells (ESCs) with non-random X chromosome inactivation (XCI) and polymorphic X chromosomes to study the dynamics of gene silencing over the inactive X chromosome by high-resolution allele-specific RNA-seq.

RESULTS

Induction of XCI by differentiation of female ESCs shows that genes proximal to the X-inactivation center are silenced earlier than distal genes, while lowly expressed genes show faster XCI dynamics than highly expressed genes. The active X chromosome shows a minor but significant increase in gene activity during differentiation, resulting in complete dosage compensation in differentiated cell types. Genes escaping XCI show little or no silencing during early propagation of XCI. Allele-specific RNA-seq of neural progenitor cells generated from the female ESCs identifies three regions distal to the X-inactivation center that escape XCI. These regions, which stably escape during propagation and maintenance of XCI, coincide with topologically associating domains (TADs) as present in the female ESCs. Also, the previously characterized gene clusters escaping XCI in human fibroblasts correlate with TADs.

CONCLUSIONS

The gene silencing observed during XCI provides further insight in the establishment of the repressive complex formed by the inactive X chromosome. The association of escape regions with TADs, in mouse and human, suggests that TADs are the primary targets during propagation of XCI over the X chromosome.

Epigenetic events regulating monoallelic gene expression

In mammals, generally it is assumed that the genes inherited from each parent are expressed to similar levels. However, it is now apparent that in non-sex chromosomes, 6-10% of genes are selected for monoallelic expression. Monoallelic expression or allelic exclusion is established either in an imprinted (parent-of-origin) or a stochastic manner. The stochastic model explains random selection while the imprinted model describes parent-of-origin specific selection of alleles for expression. Allelic exclusion occurs during X chromosome inactivation, parent-of-origin expression of imprinted genes and stochastic monoallelic expression of cell surface molecules, clustered protocadherin (PCDH) genes. Mis-regulation or loss of allelic exclusion contributes to developmental diseases. Epigenetic mechanisms are fundamental players that determine this type of expression despite a homogenous genetic background. DNA methylation and histone modifications are two mediators of the epigenetic phenomena. The majority of DNA methylation is found on cytosines of the CpG dinucleotide in mammals. Several covalent modifications of histones change the electrostatic forces between DNA and histones modifying gene expression. Long-range chromatin interactions organize chromatin into transcriptionally permissive and prohibitive regions leading to simultaneous regulation of gene expression and repression. Non-coding RNAs (ncRNAs) are also players in regulating gene expression. Together, these epigenetic mechanisms fine-tune gene expression levels essential for normal development and survival. In this review, first we discuss what is known about monoallelic gene expression. Then, we focus on the molecular mechanisms that regulate expression of three monoallelically expressed gene classes: the X-linked genes, selected imprinted genes and PCDH genes.

Single-cell transcriptogenomics reveals transcriptional exclusion of ENU-mutated alleles

Recently, great progress has been made in single cell genomics and transcriptomics. Here, we present an integrative method, termed single-cell transcriptogenomics (SCTG), in which whole exome sequencing and RNA-seq is performed concurrently on single cells. This methodology enables one to track germline and somatic variants directly from the genome to the transcriptome in individual cells. Mouse embryonic fibroblasts were treated with the powerful mutagen ethylnitrosourea (ENU) and subjected to SCTG. Interestingly, while germline variants were found to be transcribed in an allelically balanced fashion, a significantly different pattern of allelic exclusion was observed for ENU-mutant variants. These results suggest that the adverse effects of induced mutations, in contrast to germline variants, may be mitigated by allelically biased transcription. They also illustrate how SCTG can be instrumental in the direct assessment of phenotypic consequences of genomic variants.

ASAR15, A cis-acting locus that controls chromosome-wide replication timing and stability of human chromosome 15

DNA replication initiates at multiple sites along each mammalian chromosome at different times during each S phase, following a temporal replication program. We have used a Cre/loxP-based strategy to identify cis-acting elements that control this replication-timing program on individual human chromosomes. In this report, we show that rearrangements at a complex locus at chromosome 15q24.3 result in delayed replication and structural instability of human chromosome 15. Characterization of this locus identified long, RNA transcripts that are retained in the nucleus and form a “cloud” on one homolog of chromosome 15. We also found that this locus displays asynchronous replication that is coordinated with other random monoallelic genes on chromosome 15. We have named this locus ASynchronous replication and Autosomal RNA on chromosome 15, or ASAR15. Previously, we found that disruption of the ASAR6 lincRNA gene results in delayed replication, delayed mitotic condensation and structural instability of human chromosome 6. Previous studies in the mouse found that deletion of the Xist gene, from the X chromosome in adult somatic cells, results in a delayed replication and instability phenotype that is indistinguishable from the phenotype caused by disruption of either ASAR6 or ASAR15. In addition, delayed replication and chromosome instability were detected following structural rearrangement of many different human or mouse chromosomes. These observations suggest that all mammalian chromosomes contain similar cis-acting loci. Thus, under this scenario, all mammalian chromosomes contain four distinct types of essential cis-acting elements: origins, telomeres, centromeres and “inactivation/stability centers”, all functioning to promote proper replication, segregation and structural stability of each chromosome.

Mammalian cells replicate their DNA along each chromosome during a precise temporal replication program. In this report, we used a novel “chromosome-engineering” strategy to identify a DNA element that controls this replication-timing program of human chromosome 15. Characterization of this element indicated that it encodes large non-protein-coding RNAs that are retained in the nucleus and form a “cloud” on one copy of chromosome 15. Previously, we found that structural rearrangements of a similar element on human chromosome 6 causes delayed replication and structural instability of chromosome 6. Mammalian chromosomes are known to contain three distinct types of essential DNA elements that promote proper chromosome function. Thus, every chromosome contains: 1) origins of replication, which are responsible for proper initiation of DNA synthesis; 2) centromeres, which are responsible for proper chromosome separation during cell division; and 3) telomeres, which are responsible for replication and protection of the ends of linear chromosomes. Our work supports a model in which all mammalian chromosomes contain a fourth type of essential DNA element, the “inactivation/stability center”, which is responsible for proper DNA replication timing and structural stability of each chromosome.

Known unknowns for allele-specific expression and genomic imprinting effects

Recent studies have provided evidence for non-canonical imprinting effects that are associated with allele-specific expression biases at the tissue level in mice. These imprinting effects have features that are distinct from canonical imprinting effects that involve allele silencing. Here, I discuss some of the evidence for non-canonical imprinting effects in the context of random X-inactivation and epigenetic allele-specific expression effects on the autosomes. I propose several mechanisms that may underlie non-canonical imprinting effects and outline future directions and approaches to study these effects at the cellular level in vivo. The growing evidence for complex allele-specific expression effects that are cell- and developmental stage-specific has opened a new frontier for study. Currently, the function of these effects and the underlying regulatory mechanisms are largely unknown.

Variable escape from X-chromosome inactivation: identifying factors that tip the scales towards expression

In humans over 15% of X-linked genes have been shown to ‘escape’ from X-chromosome inactivation (XCI): they continue to be expressed to some extent from the inactive X chromosome. Mono-allelic expression is anticipated within a cell for genes subject to XCI, but random XCI usually results in expression of both alleles in a cell population. Using a study of allelic expression from cultured lymphoblasts and fibroblasts, many of which showed substantial skewing of XCI, we recently reported that the expression of genes lies on a contiunuum between those that are subject to inactivation, and those that escape. We now review allelic expression studies from mouse, and discuss the variability in escape seen in both humans and mice in genic expression levels, between X chromosomes and between tissues. We also discuss current knowledge of the heterochromatic features, DNA elements and three-dimensional topology of the inactive X that contribute to the balance of expression from the otherwise inactive X chromosome.

Random monoallelic expression: regulating gene expression one allele at a time

Monoallelic gene expression is a remarkable process in which transcription occurs from only one of two homologous alleles in a diploid cell. Interestingly, between 0.5% and 15% of autosomal genes exhibit random monoallelic gene expression, in which different cells express only one allele independently of the underlying genomic sequence, in a cell type-specific manner. Recently, genome-wide studies have increased our understanding of the cell type-specific incidence of random monoallelic gene expression, and how the imbalance in allelic expression is distinguished within the cell and potentially maintained across cell generations. Monoallelic gene expression is likely generated through stochastic independent regulation of the two alleles upon differentiation, and has varied implications for the cell and organism, in particular with respect to disease.

Developmental dynamics and disease potential of random monoallelic gene expression

X chromosome inactivation (XCI) and allelic exclusion of olfactory receptors or immunoglobulin loci represent classic examples of random monoallelic expression (RME). RME of some single copy genes has also been reported, but the in vivo relevance of this remains unclear. Here we identify several hundred RME genes in clonal neural progenitor cell lines derived from embryonic stem cells. RME occurs during differentiation, and, once established, the monoallelic state can be highly stable. We show that monoallelic expression also occurs in vivo, in the absence of DNA sequence polymorphism. Several of the RME genes identified play important roles in development and have been implicated in human autosomal-dominant disorders. We propose that monoallelic expression of such genes contributes to the fine-tuning of the developmental regulatory pathways they control, and, in the context of a mutation, RME can predispose to loss of function in a proportion of cells and thus contribute to disease.

Cellular resolution maps of X chromosome inactivation: implications for neural development, function, and disease

Female eutherian mammals use X chromosome inactivation (XCI) to epigenetically regulate gene expression from ∼4% of the genome. To quantitatively map the topography of XCI for defined cell types at single cell resolution, we have generated female mice that carry X-linked, Cre-activated, and nuclear-localized fluorescent reporters--GFP on one X chromosome and tdTomato on the other. Using these reporters in combination with different Cre drivers, we have defined the topographies of XCI mosaicism for multiple CNS cell types and of retinal vascular dysfunction in a model of Norrie disease. Depending on cell type, fluctuations in the XCI mosaic are observed over a wide range of spatial scales, from neighboring cells to left versus right sides of the body. These data imply a major role for XCI in generating female-specific, genetically directed, stochastic diversity in eutherian mammals on spatial scales that would be predicted to affect CNS function within and between individuals.

Random monoallelic gene expression increases upon embryonic stem cell differentiation

Random autosomal monoallelic gene expression refers to the transcription of a gene from one of two homologous alleles. We assessed the dynamics of monoallelic expression during development through an allele-specific RNA-sequencing screen in clonal populations of hybrid mouse embryonic stem cells (ESCs) and neural progenitor cells (NPCs). We identified 67 and 376 inheritable autosomal random monoallelically expressed genes in ESCs and NPCs, respectively, a 5.6-fold increase upon differentiation. Although DNA methylation and nuclear positioning did not distinguish the active and inactive alleles, specific histone modifications were differentially enriched between the two alleles. Interestingly, expression levels of 8% of the monoallelically expressed genes remained similar between monoallelic and biallelic clones. These results support a model in which random monoallelic expression occurs stochastically during differentiation and, for some genes, is compensated for by the cell to maintain the required transcriptional output of these genes.

Chromatin signature of widespread monoallelic expression

In mammals, numerous autosomal genes are subject to mitotically stable monoallelic expression (MAE), including genes that play critical roles in a variety of human diseases. Due to challenges posed by the clonal nature of MAE, very little is known about its regulation; in particular, no molecular features have been specifically linked to MAE. In this study, we report an approach that distinguishes MAE genes in human cells with great accuracy: a chromatin signature consisting of chromatin marks associated with active transcription (H3K36me3) and silencing (H3K27me3) simultaneously occurring in the gene body. The MAE signature is present in ∼20% of ubiquitously expressed genes and over 30% of tissue-specific genes across cell types. Notably, it is enriched among key developmental genes that have bivalent chromatin structure in pluripotent cells. Our results open a new approach to the study of MAE that is independent of polymorphisms, and suggest that MAE is linked to cell differentiation.

DOI:http://dx.doi.org/10.7554/eLife.01256.001

Understanding how genes are activated and silenced is one of the central challenges in modern biology. These processes underpin the development of a fertilized egg into a complex organism, and they can also lead to life-threatening diseases when they go wrong. There are two copies of each gene in a human cell, a maternal copy and a paternal copy, and it is thought that both copies are usually regulated together. However, there are exceptions to this rule: for certain genes only the maternal copy is expressed as a protein in some cells, whereas the paternal copy is expressed in other cells.

This form of gene regulation, which is called monoallelic expression, can result in neighboring cells heading down very different paths. In extreme cases, depending on the differences between the two copies of the gene, cells that express one copy may function normally, while cells where the other copy is activated will start forming tumors. However, despite these potentially grave consequences, and early results which suggested that monoallelic expression affected a large number of human and mouse genes, it has proved to be a major technical challenge to identify these genes in most cell types.

Now, Nag, Savova et al. have discovered a molecular signature that can be used to detect monoallelic expression. The signature was found in chromatin, the densely packed structure formed by DNA and proteins inside the cell nucleus. Nag, Savova et al. discovered that the genes that are subject to monoallelic expression are bound with proteins that are modified in two contrasting ways. One modification, which is usually a sign of gene silencing, is prevalent on the inactive copy of the gene, and the other, which often marks active genes, is chiefly present on the active copy.

Nag, Savova et al. report that these modifications are found in different sets of genes in different cell types, indicating distinct genome-wide patterns of monoallelic expression. The chromatin signature approach lets them estimate the fraction of human genes that are subject to monoallelic expression. This number is surprisingly high: about 20% of commonly expressed genes and more than one-third of tissue-specific genes. In a particularly intriguing finding, almost all bivalent genes—a subset of genes that are involved in determining the fate of cell during development—are estimated to become monoallelic when they are activated.

In addition to these unexpected findings, the chromatin signature approach opens the door to exploring monoallelic expression as a form of gene regulation in all types of cells and, ultimately, to understanding how it is involved in both normal development and in disease.

DOI:http://dx.doi.org/10.7554/eLife.01256.002

Random or stochastic monoallelic expressed genes are enriched for neurodevelopmental disorder candidate genes

Random or stochastic monoallelic expressed genes (StMA genes) represent a unique form of monoallelic expression where allelic choice is made at random early in development. The consequential clonal diversity provides opportunity for functional heterozygosity in tissues such as the brain, and can impact on both development and disease. We investigate the relationship of StMA expressed genes previously identified in clonal neural stem cells with the neurodevelopmental disorders autism and schizophrenia. We found that StMA genes show an overrepresentation of schizophrenia risk candidates identified by genome wide association studies from the genetic association database. Similar suggestive enrichment was also found for genes from the NHGRI genome-wide association study catalog and a psychiatric genetics consortium schizophrenia dataset although these latter more robust gene lists did not achieve statistical significance. We also examined multiple sources of copy number variation (CNV) datasets from autism and schizophrenia cohorts. After taking into account total gene numbers and CNV size, both autism and schizophrenia associated CNVs appeared to show an enrichment of StMA genes relative to the control CNV datasets. Since the StMA genes were originally identified in neural stem cells, bias due to the neural transcriptome is possible. To address this, we randomly sampled neural stem cell expressed genes and repeated the tests. After a significant number of iterations, neural stem cell expressed genes did not show an overrepresentation in autism or schizophrenia CNV datasets. Therefore, irrespective of the neural derived transcriptome, StMA genes originally identified in neural stem cells show an overrepresentation in CNVs associated with autism and schizophrenia. If this association is functional, then the regulation (or dysregulation) of this form of allelic expression status within tissues such as the brain may be a contributory risk factor for neurodevelopmental disorders and may also influence disease discordance sometimes observed in monozygotic twins.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data

MOTIVATION

High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts.

RESULTS

We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.

AVAILABILITY

The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter

Autosomal monoallelic expression: genetics of epigenetic diversity?

In mammals, relative expression of the two parental alleles of many genes is controlled by one of three major epigenetic phenomena: X chromosome inactivation, imprinting, and mitotically stable autosomal monoallelic expression (MAE). MAE affects a large fraction of human autosomal genes and introduces enormous epigenetic heterogeneity in otherwise similar cell populations. Despite its prevalence, many functional and mechanistic aspects of MAE biology remain unknown. Several lines of evidence imply that MAE establishment and maintenance are controlled by a variety of genetic elements. Based on known genomic features regulating X-inactivation and imprinting, we outline likely features of MAE-controlling elements. We also assess implications of MAE for genotype-phenotype relationship, with a focus on haploinsufficiency.

Asynchronous replication, mono-allelic expression, and long range Cis-effects of ASAR6

Mammalian chromosomes initiate DNA replication at multiple sites along their length during each S phase following a temporal replication program. The majority of genes on homologous chromosomes replicate synchronously. However, mono-allelically expressed genes such as imprinted genes, allelically excluded genes, and genes on female X chromosomes replicate asynchronously. We have identified a cis-acting locus on human chromosome 6 that controls this replication-timing program. This locus encodes a large intergenic non-coding RNA gene named Asynchronous replication and Autosomal RNA on chromosome 6, or ASAR6. Disruption of ASAR6 results in delayed replication, delayed mitotic chromosome condensation, and activation of the previously silent alleles of mono-allelic genes on chromosome 6. The ASAR6 gene resides within an ∼1.2 megabase domain of asynchronously replicating DNA that is coordinated with other random asynchronously replicating loci along chromosome 6. In contrast to other nearby mono-allelic genes, ASAR6 RNA is expressed from the later-replicating allele. ASAR6 RNA is synthesized by RNA Polymerase II, is not polyadenlyated, is restricted to the nucleus, and is subject to random mono-allelic expression. Disruption of ASAR6 leads to the formation of bridged chromosomes, micronuclei, and structural instability of chromosome 6. Finally, ectopic integration of cloned genomic DNA containing ASAR6 causes delayed replication of entire mouse chromosomes.

Mammalian chromosomes are duplicated every cell cycle during a precise temporal DNA replication program. Thus, every chromosome contains regions that are replicated early and other regions that are replicated late during each S phase. Most of the genes, present in two copies on homologous chromosomes, replicate synchronously during each S phase. Exceptions to this rule are genes located on X chromosomes, genetically imprinted genes, and genes subject to allelic exclusion. Thus, all mono-allelically expressed genes are subject to asynchronous replication, where one allele replicates before the other. Perhaps the best-studied example of asynchronous replication in mammals occurs during X inactivation in female cells. A large non-coding RNA gene called XIST, located within the X inactivation center, controls the transcriptional silencing and late replication of the inactive X chromosome. We have identified a locus on human chromosome 6 that shares many characteristics with XIST. This chromosome 6 locus encodes a large intergenic non-coding RNA gene, ASAR6, which displays random mono-allelic expression, asynchronous replication, and controls the mono-allelic expression of other genes on chromosome 6. Our work supports a model in which all mammalian chromosomes contain similar cis-acting loci that function to ensure proper chromosome replication, mitotic condensation, mono-allelic expression, and stability of individual chromosomes.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR–affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype–phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR–susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals, and currently the molecular causes of the condition are not known. Here we report the mapping, identification, and functional analysis of a susceptibility locus for IAR in stallions. The candidate region was mapped to horse chromosome 13 by SNP genotyping and GWAS of 7 IAR affected and 44 control Thoroughbred stallions. Re-sequencing and case-control analysis of functionally relevant candidate genes in the region identified FKBP6 gene as a significantly associated locus. The association was confirmed by genotyping 265 male horses of multiple breeds. FKBP6 belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. We showed that the equine FKBP6 is exclusively and monoallelically expressed in testis and sperm and has 5 different transcripts, of which 4 were novel. Overall, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in other mammals. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Candida albicans: A Model Organism for Studying Fungal Pathogens

Candida albicans is an opportunistic human fungal pathogen that causes candidiasis. As healthcare has been improved worldwide, the number of immunocompromised patients has been increased to a greater extent and they are highly susceptible to various pathogenic microbes and C. albicans has been prominent among the fungal pathogens. The complete genome sequence of this pathogen is now available and has been extremely useful for the identification of repertoire of genes present in this pathogen. The major challenge is now to assign the functions to these genes of which 13% are specific to C. albicans. Due to its close relationship with yeast Saccharomyces cerevisiae, an edge over other fungal pathogens because most of the technologies can be directly transferred to C. albicans from S. cerevisiae and it is amenable to mutation, gene disruption, and transformation. The last two decades have witnessed enormous amount of research activities on this pathogen that leads to the understanding of host-parasite interaction, infections, and disease propagation. Clearly, C. albicans has emerged as a model organism for studying fungal pathogens along with other two fungi Aspergillus fumigatus and Cryptococcus neoformans. Understanding its complete life style of C. albicans will undoubtedly be useful for developing potential antifungal drugs and tackling Candida infections. This will also shed light on the functioning of other fungal pathogens.

Allele-biased expression in differentiating human neurons: implications for neuropsychiatric disorders

Stochastic processes and imprinting, along with genetic factors, lead to monoallelic or allele-biased gene expression. Stochastic monoallelic expression fine-tunes information processing in immune cells and the olfactory system, and imprinting plays an important role in development. Recent studies suggest that both stochastic events and imprinting may be more widespread than previously considered. We are interested in allele-biased gene expression occurring in the brain because parent-of-origin effects suggestive of imprinting appear to play a role in the transmission of schizophrenia (SZ) and autism spectrum disorders (ASD) in some families. In addition, allele-biased expression could help explain monozygotic (MZ) twin discordance and reduced penetrance. The ability to study allele-biased expression in human neurons has been transformed with the advent of induced pluripotent stem cell (iPSC) technology and next generation sequencing. Using transcriptome sequencing (RNA-Seq) we identified 801 genes in differentiating neurons that were expressed in an allele-biased manner. These included a number of putative SZ and ASD candidates, such as A2BP1 (RBFOX1), ERBB4, NLGN4X, NRG1, NRG3, NRXN1, and NLGN1. Overall, there was a modest enrichment for SZ and ASD candidate genes among those that showed evidence for allele-biased expression (chi-square, p = 0.02). In addition to helping explain MZ twin discordance and reduced penetrance, the capacity to group many candidate genes affecting a variety of molecular and cellular pathways under a common regulatory process – allele-biased expression – could have therapeutic implications.