Widespread monoallelic expression on human autosomes

Allelic expression patterns of imprinted and non-imprinted genes in cancer cell lines from multiple histologies

BACKGROUND

Imprinted genes are epigenetically regulated in normal tissues to follow monoallelic expression according to the parent of origin of each allele. Some of these patterns are dysregulated in cancer.

RESULTS

We developed a novel computational multi-omic pipeline to evaluate monoallelic and biallelic expression patterns based on matched RNA-seq expression data, whole-exome sequencing information, and copy number data. We analyzed allelic expression of the entire genes, individual isoforms, and each exon of 59,283 autosomal protein-coding and ncRNA genes, with a focus on 94 genes previously reported to be imprinted. We analyzed 108 cell lines from 9 different tumor histologies using molecular data from the DepMap Portal for the Cancer Cell Line Encyclopedia. Allelic expression patterns of imprinted genes and isoforms in tumor cells were variable. We also identified additional genes and isoforms with predominantly monoallelic expression due to a variety of potential mechanisms. We provide a novel public dataset of transcriptome-wide allelic expression patterns in cell lines from diverse tumor categories, which can serve as a resource for future cancer studies. We examined associations of in vitro cell line response to antitumor agents and repurposed drugs with allelic patterns and overall levels of isoform expression of imprinted genes and of additional genes with predominantly monoallelic expression. Drug response was associated with isoform expression patterns of multiple imprinted genes including CPA4, DGCR6, DNMT1, GNAS, GRB10, H19, NAA60, OSBPL5, PHACTR2, and ZFAT, predominantly monoallelically expressed MAP2K5 and BCLAF1, and additional predominantly monoallelically expressed genes. Multiple associations may be related to mechanisms of drug activity, including associations between the response to the DNA damaging agents and allelic expression of ZFAT, CDC27, and BCLAF1 isoforms, and the response to inhibitors of multiple signaling pathways with expression patterns of GNAS isoforms.

CONCLUSIONS

Tumor cells have a range of monoallelic and biallelic expression patterns in both imprinted and non-imprinted genes and are likely affected by the complex interplay among changes in allelic expression, sequence variants, copy number changes, and expression changes of biologically important genes. Multiple isoform-specific patterns of allelic expression were associated with drug response, indicating complex mechanisms of cancer chemoresistance.

Beyond genes and environment: mapping biological stochasticity in aging

Aging is characterized by extensive variability in the onset of morbidity and mortality, even in genetically identical populations with carefully controlled environments. This points to the important role stochasticity plays in shaping the divergent aging process between individual organisms. Here, we survey how stochastic factors at the level of molecules, cells, tissues, and organisms manifest in and impact the aging process, with a focus on the nematode Caenorhabditis elegans. Findings of stochasticity in C. elegans give additional insights for aspects of aging in the more complex settings of mammals with parallels drawn between organisms when appropriate. The emerging understanding of the stochastic contributors to longevity will enhance research strategies and medical interventions for personalized medicine.

Genome wide analysis of allele-specific circular RNAs in mammals and their role in cell proliferation

Circular RNAs (circRNAs) are a large class of widely expressed RNAs with covalently closed continuous structures. However, it is currently unknown if circRNAs shows allele-specific expression, as are the consequences of genetic variation on their circularization efficiency and subsequent biological function. Here, we propose a novel pipeline, ASE-circRNA, to accurately quantify both circRNA and their related linear RNA for each allele, and then assess the allele-specificity of the expression of a circular RNA. We identified and analyzed allele-specific circRNAs from human tissue, as well as brains from reciprocal crosses between pairs of highly divergent strains of both mice and pigs by next generation sequencing. Droplet digital PCR (ddPCR) was used to confirm the circularization efficiency measured by next generation sequencing. We found that variation in intron sequences affect the circularization efficiency of circRNAs. Furthermore, we demonstrate that a circRNA, circHK1, regulates the expression of POLR2A to influence the rate of cell proliferation. Our study provides new insight into the molecular mechanisms impacted by variation in genome sequence in the origin of human disease and phenotype.

A potential role for monoallelic expression in penetrance of autosomal dominant inborn errors of immunity

In this article, we discuss a recent study, where autosomal monoallelic expression of genes underlying Inborn Errors of Immunity were investigated. About 2-10% of genes are predominantly transcribed from a single allele leading to autosomal random monoallelic expression (I). If this is skewed in a cell population from an individual with an autosomal dominant inborn error of immunity, this can lead to a mild to no phenotype (incomplete penetrance) if the wildtype allele is favored (II), or to more severe disease presentation if the variant allele is favored (III).

Monoallelic expression can govern penetrance of inborn errors of immunity

Inborn errors of immunity (IEIs) are genetic disorders that underlie susceptibility to infection, autoimmunity, autoinflammation, allergy and/or malignancy1. Incomplete penetrance is common among IEIs despite their monogenic basis2. Here we investigate the contribution of autosomal random monoallelic expression (aRMAE), a somatic commitment to the expression of one allele3,4, to phenotypic variability observed in families with IEIs. Using a clonal primary T cell system to assess aRMAE status of genes in healthy individuals, we find that 4.30% of IEI genes and 5.20% of all genes undergo aRMAE. Perturbing H3K27me3 and DNA methylation alters allele expression commitment, in support of two proposed mechanisms5,6 for the regulation of aRMAE. We tested peripheral blood mononuclear cells from individuals with IEIs with shared genetic lesions but discordant clinical phenotypes for aRMAE. Among two relatives who were heterozygous for a mutation in PLCG2 (delEx19), an antibody deficiency phenotype corresponds to selective mutant allele expression in B cells. By contrast, among relatives who were heterozygous for a mutation in JAK1 (c.2099G>A; p.S700N), the unaffected carrier T cells predominantly expressed the wild-type JAK1 allele, whereas the affected carrier T cells exhibited biallelic expression. Allelic expression bias was also documented in phenotypically discordant family members with mutations in STAT1 and CARD11. This study highlights the importance of considering both the genotype and the 'transcriptotype' in analyses of the penetrance and expressivity of monogenic disorders.

Comprehensive body fluid identification and contributor assignment by combining targeted sequencing of mRNA and coding region SNPs

Forensic genetic analyses aim to retrieve as much information as possible from biological trace material recovered from crime scenes. While standard short tandem repeat (STR) profiling is essential to individualize biological traces, its significance is diminished in crime scenarios where the presence of a suspect's DNA is acknowledged by all parties. In such cases, forensic (m)RNA analysis can provide crucial contextualizing information on the source level about a trace's composition, i.e., body fluids/tissues, and has therefore emerged as a powerful tool for modern forensic investigations. However, the question which of several suspects contributed a specific component (body fluid) to a mixed trace cannot be answered by RNA analysis using conventional methods. This individualizing information is stored within the sequence of the mRNA transcripts. Massively parallel sequencing (MPS) represents a promising alternative, offering not only higher multiplex capacity, but also the typing of individual coding region SNPs (cSNPs) to enable the assignment of contributors to mixture components, thereby reducing the risk of association fallacies. Herein, we describe the development of an extensive mRNA/cSNP panel for targeted sequencing on the IonTorrent S5 platform. Our panel comprises 30 markers for the detection of six body fluids/tissues (blood, saliva, semen, skin, vaginal and menstrual secretion), along with 70 linkage-controlled cSNPs for contributor assignment. It exhibited high reliable detection sensitivity with RNA inputs down to 0.75 ng and a conservatively calculated probability of identity of 0.03 - 6 % for individual body fluid-specific cSNP profiles. Limitations and areas for future work include RNA-related allele imbalances, inclusion of markers to correctly identify rectal mucosa and the optimization of specific markers. In summary, our new panel is intended to be a major step forward to interpret biological evidence at sub-source and source level based on cSNP attribution of a body fluid component to a suspect and victim, respectively.

Identification of candidate causal variants and target genes at 41 breast cancer risk loci through differential allelic expression analysis

Understanding breast cancer genetic risk relies on identifying causal variants and candidate target genes in risk loci identified by genome-wide association studies (GWAS), which remains challenging. Since most loci fall in active gene regulatory regions, we developed a novel approach facilitated by pinpointing the variants with greater regulatory potential in the disease's tissue of origin. Through genome-wide differential allelic expression (DAE) analysis, using microarray data from 64 normal breast tissue samples, we mapped the variants associated with DAE (daeQTLs). Then, we intersected these with GWAS data to reveal candidate risk regulatory variants and analysed their cis-acting regulatory potential. Finally, we validated our approach by extensive functional analysis of the 5q14.1 breast cancer risk locus. We observed widespread gene expression regulation by cis-acting variants in breast tissue, with 65% of coding and noncoding expressed genes displaying DAE (daeGenes). We identified over 54 K daeQTLs for 6761 (26%) daeGenes, including 385 daeGenes harbouring variants previously associated with BC risk. We found 1431 daeQTLs mapped to 93 different loci in strong linkage disequilibrium with risk-associated variants (risk-daeQTLs), suggesting a link between risk-causing variants and cis-regulation. There were 122 risk-daeQTL with stronger cis-acting potential in active regulatory regions with protein binding evidence. These variants mapped to 41 risk loci, of which 29 had no previous report of target genes and were candidates for regulating the expression levels of 65 genes. As validation, we identified and functionally characterised five candidate causal variants at the 5q14.1 risk locus targeting the ATG10 and ATP6AP1L genes, likely acting via modulation of alternative transcription and transcription factor binding. Our study demonstrates the power of DAE analysis and daeQTL mapping to identify causal regulatory variants and target genes at breast cancer risk loci, including those with complex regulatory landscapes. It additionally provides a genome-wide resource of variants associated with DAE for future functional studies.

Detection of allele-specific expression in spatial transcriptomics with spASE

Spatial transcriptomics technologies permit the study of the spatial distribution of RNA at near-single-cell resolution genome-wide. However, the feasibility of studying spatial allele-specific expression (ASE) from these data remains uncharacterized. Here, we introduce spASE, a computational framework for detecting and estimating spatial ASE. To tackle the challenges presented by cell type mixtures and a low signal to noise ratio, we implement a hierarchical model involving additive mixtures of spatial smoothing splines. We apply our method to allele-resolved Visium and Slide-seq from the mouse cerebellum and hippocampus and report new insight into the landscape of spatial and cell type-specific ASE therein.

Gene-specific somatic epigenetic mosaicism of FDFT1 underlies a non-hereditary localized form of porokeratosis

Porokeratosis is a clonal keratinization disorder characterized by solitary, linearly arranged, or generally distributed multiple skin lesions. Previous studies showed that genetic alterations in MVK, PMVK, MVD, or FDPS-genes in the mevalonate pathway-cause hereditary porokeratosis, with skin lesions harboring germline and lesion-specific somatic variants on opposite alleles. Here, we identified non-hereditary porokeratosis associated with epigenetic silencing of FDFT1, another gene in the mevalonate pathway. Skin lesions of the generalized form had germline and lesion-specific somatic variants on opposite alleles in FDFT1, representing FDFT1-associated hereditary porokeratosis identified in this study. Conversely, lesions of the solitary or linearly arranged localized form had somatic bi-allelic promoter hypermethylation or mono-allelic promoter hypermethylation with somatic genetic alterations on opposite alleles in FDFT1, indicating non-hereditary porokeratosis. FDFT1 localization was uniformly diminished within the lesions, and lesion-derived keratinocytes showed cholesterol dependence for cell growth and altered expression of genes related to cell-cycle and epidermal development, confirming that lesions form by clonal expansion of FDFT1-deficient keratinocytes. In some individuals with the localized form, gene-specific promoter hypermethylation of FDFT1 was detected in morphologically normal epidermis adjacent to methylation-related lesions but not distal to these lesions, suggesting that asymptomatic somatic epigenetic mosaicism of FDFT1 predisposes certain skin areas to the disease. Finally, consistent with its genetic etiology, topical statin treatment ameliorated lesions in FDFT1-deficient porokeratosis. In conclusion, we identified bi-allelic genetic and/or epigenetic alterations of FDFT1 as a cause of porokeratosis and shed light on the pathogenesis of skin mosaicism involving clonal expansion of epigenetically altered cells.

Extrinsic and intrinsic drivers of natural killer cell clonality

Clonal expansion of antigen-specific lymphocytes is the fundamental mechanism enabling potent adaptive immune responses and the generation of immune memory. Accompanied by pronounced epigenetic remodeling, the massive proliferation of individual cells generates a critical mass of effectors for the control of acute infections, as well as a pool of memory cells protecting against future pathogen encounters. Classically associated with the adaptive immune system, recent work has demonstrated that innate immune memory to human cytomegalovirus (CMV) infection is stably maintained as large clonal expansions of natural killer (NK) cells, raising questions on the mechanisms for clonal selection and expansion in the absence of re-arranged antigen receptors. Here, we discuss clonal NK cell memory in the context of the mechanisms underlying clonal competition of adaptive lymphocytes and propose alternative selection mechanisms that might decide on the clonal success of their innate counterparts. We propose that the integration of external cues with cell-intrinsic sources of heterogeneity, such as variegated receptor expression, transcriptional states, and somatic variants, compose a bottleneck for clonal selection, contributing to the large size of memory NK cell clones.

Inborn errors of immunity: an expanding universe of disease and genetic architecture

Inborn errors of immunity (IEIs) are generally considered to be rare monogenic disorders of the immune system that cause immunodeficiency, autoinflammation, autoimmunity, allergy and/or cancer. Here, we discuss evidence that IEIs need not be rare disorders or exclusively affect the immune system. Namely, an increasing number of patients with IEIs present with severe dysregulations of the central nervous, digestive, renal or pulmonary systems. Current challenges in the diagnosis of IEIs that result from the segregated practice of specialized medicine could thus be mitigated, in part, by immunogenetic approaches. Starting with a brief historical overview of IEIs, we then discuss the technological advances that are facilitating the immunogenetic study of IEIs, progress in understanding disease penetrance in IEIs, the expanding universe of IEIs affecting distal organ systems and the future of genetic, biochemical and medical discoveries in this field.

Maternal histone methyltransferases antagonistically regulate monoallelic expression in C. elegans

Undefined epigenetic programs act to probabilistically silence individual autosomal alleles, generating unique individuals, even from genetic clones. This sort of random monoallelic expression can explain variation in traits and diseases that differences in genes and environments cannot. Here, we developed the nematode Caenorhabditis elegans to study monoallelic expression in whole tissues, and defined a developmental genetic regulation pathway. We found maternal H3K9 histone methyltransferase (HMT) SET-25/SUV39/G9a works with HPL-2/HP1 and LIN-61/L3MBTL2 to randomly silence alleles in the intestinal progenitor E-cell of 8-cell embryos to cause monoallelic expression. SET-25 was antagonized by another maternal H3K9 HMT, MET-2/SETDB1, which works with LIN-65/ATF7ZIP and ARLE-14/ARL14EP to prevent monoallelic expression. The HMT-catalytic SET domains of both MET-2 and SET-25 were required for regulating monoallelic expression. Our data support a model wherein SET-25 and MET-2 regulate histones during development to generate patterns of somatic monoallelic expression that are persistent but not heritable.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Genomic Imprinting and Random Monoallelic Expression

The review discusses the mechanisms of monoallelic expression, such as genomic imprinting, in which gene transcription depends on the parental origin of the allele, and random monoallelic transcription. Data on the regulation of gene activity in the imprinted regions are summarized with a particular focus on the areas controlling imprinting and factors influencing the variability of the imprintome. The prospects of studies of the monoallelic expression are discussed.

Random Allelic Expression in Inherited Retinal Disease Genes

Inherited retinal diseases (IRDs) are a significant contributor to visual loss in children and young adults, falling second only to diabetic retinopathy. Understanding the pathogenic mechanisms of IRDs remains paramount. Some autosomal genes exhibit random allelic expression (RAE), similar to X-chromosome inactivation. This study identifies RAE genes in IRDs. Genes in the Retinal Information Network were cross-referenced with the recent literature to identify expression profiles, RAE, or biallelic expression (BAE). Loss-of-function intolerance (LOFI) was determined by cross-referencing the existing literature. Molecular and biological pathways that are significantly enriched were evaluated using gene ontology. A total of 184 IRD-causing genes were evaluated. Of these, 31 (16.8%) genes exhibited RAE. LOFI was exhibited in 6/31 (19.4%) of the RAE genes and 18/153 (11.8%) of the BAE genes. Brain tissue exhibited BAE in 107/128 (83.6%) genes for both sexes. The molecular pathways significantly enriched among BAE genes were photoreceptor activity, tubulin binding, and nucleotide/ribonucleotide binding. The biologic pathways significantly enriched for RAE genes were equilibrioception, parallel actin filament bundle assembly, photoreceptor cell outer segment organization, and protein depalmitoylation. Allele-specific expression may be a mechanism underlying IRD phenotypic variability, with clonal populations of embryologic precursor cells exhibiting RAE. Brain tissue preferentially exhibited BAE, possibly due to selective pressures against RAE. Pathways critical for cellular and visual function were enriched in BAE, which may offer a survival benefit.

The transcriptional legacy of developmental stochasticity

Genetic and environmental variation are key contributors during organism development, but the influence of minor perturbations or noise is difficult to assess. This study focuses on the stochastic variation in allele-specific expression that persists through cell divisions in the nine-banded armadillo (Dasypus novemcinctus). We investigated the blood transcriptome of five wild monozygotic quadruplets over time to explore the influence of developmental stochasticity on gene expression. We identify an enduring signal of autosomal allelic variability that distinguishes individuals within a quadruplet despite their genetic similarity. This stochastic allelic variation, akin to X-inactivation but broader, provides insight into non-genetic influences on phenotype. The presence of stochastically canalized allelic signatures represents a novel axis for characterizing organismal variability, complementing traditional approaches based on genetic and environmental factors. We also developed a model to explain the inconsistent penetrance associated with these stochastically canalized allelic expressions. By elucidating mechanisms underlying the persistence of allele-specific expression, we enhance understanding of development's role in shaping organismal diversity.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Nature and nurture: understanding phenotypic variation in inborn errors of immunity

The overall disease burden of pediatric infection is high, with widely varying clinical outcomes including death. Among the most vulnerable children, those with inborn errors of immunity, reduced penetrance and variable expressivity are common but poorly understood. There are several genetic mechanisms that influence phenotypic variation in inborn errors of immunity, as well as a body of knowledge on environmental influences and specific pathogen triggers. Critically, recent advances are illuminating novel nuances for fundamental concepts on disease penetrance, as well as raising new areas of inquiry. The last few decades have seen the identification of almost 500 causes of inborn errors of immunity, as well as major advancements in our ability to characterize somatic events, the microbiome, and genotypes across large populations. The progress has not been linear, and yet, these developments have accumulated into an enhanced ability to diagnose and treat inborn errors of immunity, in some cases with precision therapy. Nonetheless, many questions remain regarding the genetic and environmental contributions to phenotypic variation both within and among families. The purpose of this review is to provide an updated summary of key concepts in genetic and environmental contributions to phenotypic variation within inborn errors of immunity, conceptualized as including dynamic, reciprocal interplay among factors unfolding across the key dimension of time. The associated findings, potential gaps, and implications for research are discussed in turn for each major influencing factor. The substantial challenge ahead will be to organize and integrate information in such a way that accommodates the heterogeneity within inborn errors of immunity to arrive at a more comprehensive and accurate understanding of how the immune system operates in health and disease. And, crucially, to translate this understanding into improved patient care for the millions at risk for serious infection and other immune-related morbidity.

Foreign RNA spike-ins enable accurate allele-specific expression analysis at scale

MOTIVATION

Analysis of allele-specific expression is strongly affected by the technical noise present in RNA-seq experiments. Previously, we showed that technical replicates can be used for precise estimates of this noise, and we provided a tool for correction of technical noise in allele-specific expression analysis. This approach is very accurate but costly due to the need for two or more replicates of each library. Here, we develop a spike-in approach which is highly accurate at only a small fraction of the cost.

RESULTS

We show that a distinct RNA added as a spike-in before library preparation reflects technical noise of the whole library and can be used in large batches of samples. We experimentally demonstrate the effectiveness of this approach using combinations of RNA from species distinguishable by alignment, namely, mouse, human, and Caenorhabditis elegans. Our new approach, controlFreq, enables highly accurate and computationally efficient analysis of allele-specific expression in (and between) arbitrarily large studies at an overall cost increase of ∼5%.

AVAILABILITY AND IMPLEMENTATION

Analysis pipeline for this approach is available at GitHub as R package controlFreq (github.com/gimelbrantlab/controlFreq).

Examining parent-of-origin effects on transcription and RNA methylation in mediating aggressive behavior in honey bees (Apis mellifera)

Conflict between genes inherited from the mother (matrigenes) and the father (patrigenes) is predicted to arise during social interactions among offspring if these genes are not evenly distributed among offspring genotypes. This intragenomic conflict drives parent-specific transcription patterns in offspring resulting from parent-specific epigenetic modifications. Previous tests of the kinship theory of intragenomic conflict in honey bees (Apis mellifera) provided evidence in support of theoretical predictions for variation in worker reproduction, which is associated with extreme variation in morphology and behavior. However, more subtle behaviors - such as aggression - have not been extensively studied. Additionally, the canonical epigenetic mark (DNA methylation) associated with parent-specific transcription in plant and mammalian model species does not appear to play the same role as in honey bees, and thus the molecular mechanisms underlying intragenomic conflict in this species is an open area of investigation. Here, we examined the role of intragenomic conflict in shaping aggression in honey bee workers through a reciprocal cross design and Oxford Nanopore direct RNA sequencing. We attempted to probe the underlying regulatory basis of this conflict through analyses of parent-specific RNA m6A and alternative splicing patterns. We report evidence that intragenomic conflict occurs in the context of honey bee aggression, with increased paternal and maternal allele-biased transcription in aggressive compared to non-aggressive bees, and higher paternal allele-biased transcription overall. However, we found no evidence to suggest that RNA m6A or alternative splicing mediate intragenomic conflict in this species.

JAKs and STATs from a Clinical Perspective: Loss-of-Function Mutations, Gain-of-Function Mutations, and Their Multidimensional Consequences

The JAK/STAT signaling pathway plays a key role in cytokine signaling and is involved in development, immunity, and tumorigenesis for nearly any cell. At first glance, the JAK/STAT signaling pathway appears to be straightforward. However, on closer examination, the factors influencing the JAK/STAT signaling activity, such as cytokine diversity, receptor profile, overlapping JAK and STAT specificity among non-redundant functions of the JAK/STAT complexes, positive regulators (e.g., cooperating transcription factors), and negative regulators (e.g., SOCS, PIAS, PTP), demonstrate the complexity of the pathway's architecture, which can be quickly disturbed by mutations. The JAK/STAT signaling pathway has been, and still is, subject of basic research and offers an enormous potential for the development of new methods of personalized medicine and thus the translation of basic molecular research into clinical practice beyond the use of JAK inhibitors. Gain-of-function and loss-of-function mutations in the three immunologically particularly relevant signal transducers STAT1, STAT3, and STAT6 as well as JAK1 and JAK3 present themselves through individual phenotypic clinical pictures. The established, traditional paradigm of loss-of-function mutations leading to immunodeficiency and gain-of-function mutation leading to autoimmunity breaks down and a more differentiated picture of disease patterns evolve. This review is intended to provide an overview of these specific syndromes from a clinical perspective and to summarize current findings on pathomechanism, symptoms, immunological features, and therapeutic options of STAT1, STAT3, STAT6, JAK1, and JAK3 loss-of-function and gain-of-function diseases.

Hierarchical DNA branch assembly-encoded fluorescent nanoladders for single-cell transcripts imaging

Spatial visualization of single-cell transcripts is limited by signal specificity and multiplexing. Here, we report hierarchical DNA branch assembly-encoded fluorescent nanoladders, which achieve denoised and highly multiplexed signal amplification for single-molecule transcript imaging. This method first offers independent RNA-primed rolling circle amplification without nonspecific amplification based on circular DNAzyme. It then executes programmable DNA branch assembly on these amplicons to encode virtual signals for visualizing numbers of targets by FISH. In theory, more virtual signals can be encoded via the increase of detection spectral channels and repeats of the same sequences on barcode. Our method almost eliminates nonspecific amplification in fixed cells (reducing nonspecific spots of single cells from 16 to nearly zero), and achieves simultaneous quantitation of nine transcripts by using only two detection spectral channels. We demonstrate accurate RNA profiling in different cancer cells, and reveal diverse localization patterns for spatial regulation of transcripts.

Foreign RNA spike-ins enable accurate allele-specific expression analysis at scale

Motivation

Analysis of allele-specific expression is strongly affected by the technical noise present in RNA-seq experiments. Previously, we showed that technical replicates can be used for precise estimates of this noise, and we provided a tool for correction of technical noise in allele-specific expression analysis. This approach is very accurate but costly due to the need for two or more replicates of each library. Here, we develop a spike-in approach that is highly accurate at only a small fraction of the cost.

Results

We show that a distinct RNA added as a spike-in before library preparation reflects technical noise of the whole library and can be used in large batches of samples. We experimentally demonstrate the effectiveness of this approach using combinations of RNA from species distinguishable by alignment, namely, mouse, human, and C.elegans . Our new approach, controlFreq , enables highly accurate and computationally efficient analysis of allele-specific expression in (and between) arbitrarily large studies at an overall cost increase of ~ 5%.

Availability

Analysis pipeline for this approach is available at GitHub as R package controlFreq ( github.com/gimelbrantlab/controlFreq ).

Contact

agimelbrant@altius.org.

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.

Targeted S5 RNA sequencing assay for the identification and direct association of common body fluids with DNA donors in mixtures

The evidentiary value of DNA profiles varies depending upon the context in which the DNA was found. Linking a DNA profile to a particular cellular phenotype in mixtures may aid in assessing its evidentiary relevance and value. We report the development of two dual-function high-resolution messenger RNA (mRNA) sequencing assays that can each identify the presence of 6 body fluids/tissues (blood, semen, saliva, vaginal secretions, menstrual blood, skin) and, via coding region SNPs (cSNPs) present in the body fluid-specific mRNA transcripts, directly associate particular body fluids with their specific DNA donors in mixtures. The original blood, semen, and saliva (BSS) assay contains 23 cSNPs for blood, semen, and saliva, while the expanded 6F (all 6 fluids/tissues) assay encompasses the BSS assay and also contains 23 additional cSNPs for vaginal secretions, menstrual blood, and skin. Software tools were developed to infer the identity of the body fluids present as well as providing the corresponding cSNP genotypes. Concomitant genomic DNA assays (BSS-d and 6F-d), required to genotype the same cSNPs from persons of interest/inferred contributors to the body fluid mixture, were also developed. Body fluid specificity was demonstrated by the ability to identify the body fluid origin of single-source and two-fluid admixtures. The discriminatory power (European Caucasians) for all body fluids is 0.957-0.997, with linkage disequilibrium considered. Reciprocal body fluid admixtures (mixture pairs with the same two donors but reversed body fluid types) were used to demonstrate the ability to identify the body fluid source of origin as well as associate the donor of each of the two fluids.

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.

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.

Nuclear-specific gene expression in heterokaryons of the filamentous ascomycete Neurospora tetrasperma

Heterokaryosis is a system in which genetically distinct nuclei coexist within the same cytoplasm. While heterokaryosis dominates the life cycle of many fungal species, the transcriptomic changes associated with the transition from homokaryosis to heterokaryosis is not well understood. Here, we analyse gene expression profiles of homokaryons and heterokaryons from three phylogenetically and reproductively isolated lineages of the filamentous ascomycete Neurospora tetrasperma. We show that heterokaryons are transcriptionally distinct from homokaryons in the sexual stage of development, but not in the vegetative stage, suggesting that the phenotypic switch to fertility in heterokaryons is associated with major changes in gene expression. Heterokaryon expression is predominantly defined by additive effects of its two nuclear components. Furthermore, allele-specific expression analysis of heterokaryons with varying nuclear ratios show patterns of expression ratios strongly dependent on nuclear ratios in the vegetative stage. By contrast, in the sexual stage, strong deviations of expression ratios indicate a co-regulation of nuclear gene expression in all three lineages. Taken together, our results show two levels of expression control: additive effects suggest a nuclear level of expression, whereas co-regulation of gene expression indicate a heterokaryon level of control.

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.

Binary outcomes of enhancer activity underlie stable random monoallelic expression

Mitotically stable random monoallelic gene expression (RME) is documented for a small percentage of autosomal genes. We developed an in vivo genetic model to study the role of enhancers in RME using high-resolution single-cell analysis of natural killer (NK) cell receptor gene expression and enhancer deletions in the mouse germline. Enhancers of the RME NK receptor genes were accessible and enriched in H3K27ac on silent and active alleles alike in cells sorted according to allelic expression status, suggesting enhancer activation and gene expression status can be decoupled. In genes with multiple enhancers, enhancer deletion reduced gene expression frequency, in one instance converting the universally expressed gene encoding NKG2D into an RME gene, recapitulating all aspects of natural RME including mitotic stability of both the active and silent states. The results support the binary model of enhancer action, and suggest that RME is a consequence of general properties of gene regulation by enhancers rather than an RME-specific epigenetic program. Therefore, many and perhaps all genes may be subject to some degree of RME. Surprisingly, this was borne out by analysis of several genes that define different major hematopoietic lineages, that were previously thought to be universally expressed within those lineages: the genes encoding NKG2D, CD45, CD8α, and Thy-1. We propose that intrinsically probabilistic gene allele regulation is a general property of enhancer-controlled gene expression, with previously documented RME representing an extreme on a broad continuum.

Genetic, epigenetic, and environmental mechanisms govern allele-specific gene expression

Allele-specific expression (ASE) is a phenomenon in which one allele is preferentially expressed over the other. Genetic and epigenetic factors cause ASE by altering the final composition of a gene's product, leading to expression imbalances that can have functional consequences on phenotypes. Environmental signals also impact allele-specific expression, but how they contribute to this cross talk remains understudied. Here, we explored how genotype, parent-of-origin, tissue, sex, and dietary fat simultaneously influence ASE biases. Male and female mice from a F₁ reciprocal cross of the LG/J and SM/J strains were fed a high or low fat diet. We harnessed strain-specific variants to distinguish between two ASE classes: parent-of-origin-dependent (unequal expression based on parental origin) and sequence-dependent (unequal expression based on nucleotide identity). We present a comprehensive map of ASE patterns in 2853 genes across three tissues and nine environmental contexts. We found that both ASE classes are highly dependent on tissue and environmental context. They vary across metabolically relevant tissues, between males and females, and in response to dietary fat. We also found 45 genes with inconsistent ASE biases that switched direction across tissues and/or environments. Finally, we integrated ASE and QTL data from published intercrosses of the LG/J and SM/J strains. Our ASE genes are often enriched in QTLs for metabolic and musculoskeletal traits, highlighting how this orthogonal approach can prioritize candidate genes. Together, our results provide novel insights into how genetic, epigenetic, and environmental mechanisms govern allele-specific expression, which is an essential step toward deciphering the genotype-to-phenotype map.

Biologic characterization of ABCA3 variants in lung tissue from infants and children with ABCA3 deficiency

ABCA3 is a phospholipid transporter protein required for surfactant assembly in lamellar bodies of alveolar type II cells. Biallelic pathogenic ABCA3 variants cause severe neonatal respiratory distress syndrome or childhood interstitial lung disease. However, ABCA3 genotype alone does not explain the diversity in disease presentation, severity, and progression. Additionally, monoallelic ABCA3 variants have been reported in infants and children with ABCA3-deficient phenotypes. The effects of most ABCA3 variants identified in patients have not been characterized at the RNA level. ABCA3 allele-specific expression occurs in some cell types due to epigenetic regulation. We obtained lung tissue at transplant or autopsy from 16 infants and children with ABCA3 deficiency due to compound heterozygous ABCA3 variants for biologic characterization of the predicted effects of ABCA3 variants at the RNA level and determination of ABCA3 allele expression. We extracted DNA and RNA from frozen lung tissue and reverse-transcribed cDNA from mRNA. We performed Sanger sequencing to assess allele-specific expression by comparing the heights of variant nucleotide peaks in amplicons from genomic DNA and cDNA. We found similar genomic and cDNA variant nucleotide peak heights and no evidence of allele-specific expression among explant or autopsy samples with biallelic missense ABCA3 variants (n = 6). We observed allele-specific expression of missense alleles in trans with frameshift (n = 4) or nonsense (n = 1) variants, attributable to nonsense-mediated decay. The missense variant c.53 A > G;p.Gln18Arg, located near an exon-intron junction, encoded abnormal splicing with skipping of exon 4. Biologic characterization of ABCA3 variants can inform discovery of variant-specific disease mechanisms.

Cubilin, the intrinsic factor-vitamin B12 receptor

Cubilin (CUBN), the intrinsic factor-vitamin B12 receptor is a large endocytic protein involved in various physiological functions: vitamin B12 uptake in the gut; reabsorption of albumin and maturation of vitamin D in the kidney; nutrient delivery during embryonic development. Cubilin is an atypical receptor, peripherally associated to the plasma membrane. The transmembrane proteins amnionless (AMN) and Lrp2/Megalin are the currently known molecular partners contributing to plasma membrane transport and internalization of Cubilin. The role of Cubilin/Amn complex in the handling of vitamin B12 in health and disease has extensively been studied and so is the role of the Cubilin-Lrp2 tandem in renal pathophysiology. Accumulating evidence strongly supports a role of Cubilin in some developmental defects including impaired closure of the neural tube. Are these defects primarily caused by the dysfunction of a specific Cubilin ligand or are they secondary to impaired vitamin B12 or protein uptake? We will present the established Cubilin functions, discuss the developmental data and provide an overview of the emerging implications of Cubilin in the field of cardiovascular disease and cancer pathogenesis.

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.

Identification of Monoallelically Expressed Genes Associated with Economic Traits in Hanwoo (Korean Native Cattle)

Hanwoo, an indigenous Korean cattle breed, has been genetically improved by selecting superior sires called Korean-proven bulls. However, cows still contribute half of the genetic stock of their offspring, and allelic-specific expressed genes have potential, as selective targets of cows, to enhance genetic gain. The aim of this study is to identify genes that have MAEs based on both the genome and transcriptome and to estimate their effects on breeding values (BVs) for economically important traits in Hanwoo. We generated resequencing data for the parents and RNA-sequencing data for the muscle, fat, and brain tissues of the offspring. A total of 3801 heterozygous single nucleotide polymorphisms (SNPs) in offspring were identified and they were located in 1569 genes. Only 14 genes showed MAE (seven expressing maternal alleles and seven expressing paternal alleles). Tissue-specific MAE was observed, and LANCL1 showed maternal allele expression across all tissues. MAE genes were enriched for the biological process of cell death and angiogenesis, which included ACKR3 and PDCL3 genes, whose SNPs were significantly associated with BVs of lean meat production-related traits, such as weight at 12 months of age, carcass weight, and loin eye area. In the current study, monoallelically expressed genes were identified in various adult tissues and these genes were associated with genetic capacity in Hanwoo.

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.

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.

Genome-Wide Assessment Characteristics of Genes Overlapping Copy Number Variation Regions in Duroc Purebred Population

Copy number variations (CNVs) are important structural variations that can cause significant phenotypic diversity. Reliable CNVs mapping can be achieved by identification of CNVs from different genetic backgrounds. Investigations on the characteristics of overlapping between CNV regions (CNVRs) and protein-coding genes (CNV genes) or miRNAs (CNV-miRNAs) can reveal the potential mechanisms of their regulation. In this study, we used 50 K SNP arrays to detect CNVs in Duroc purebred pig. A total number of 211 CNVRs were detected with a total length of 118.48 Mb, accounting for 5.23% of the autosomal genome sequence. Of these CNVRs, 32 were gains, 175 losses, and four contained both types (loss and gain within the same region). The CNVRs we detected were non-randomly distributed in the swine genome and were significantly enriched in the segmental duplication and gene density region. Additionally, these CNVRs were overlapping with 1,096 protein-coding genes (CNV-genes), and 39 miRNAs (CNV-miRNAs), respectively. The CNV-genes were enriched in terms of dosage-sensitive gene list. The expression of the CNV genes was significantly higher than that of the non-CNV genes in the adult Duroc prostate. Of all detected CNV genes, 22.99% genes were tissue-specific (TSI > 0.9). Strong negative selection had been underway in the CNV-genes as the ones that were located entirely within the loss CNVRs appeared to be evolving rapidly as determined by the median dN plus dS values. Non-CNV genes tended to be miRNA target than CNV-genes. Furthermore, CNV-miRNAs tended to target more genes compared to non-CNV-miRNAs, and a combination of two CNV-miRNAs preferentially synergistically regulated the same target genes. We also focused our efforts on examining CNV genes and CNV-miRNAs functions, which were also involved in the lipid metabolism, including DGAT1, DGAT2, MOGAT2, miR143, miR335, and miRLET7. Further molecular experiments and independent large studies are needed to confirm our findings.

A Comprehensive Characterization of Monoallelic Expression During Hematopoiesis and Leukemogenesis via Single-Cell RNA-Sequencing

Single-cell RNA-sequencing (scRNA-seq) is becoming a powerful tool to investigate monoallelic expression (MAE) in various developmental and pathological processes. However, our knowledge of MAE during hematopoiesis and leukemogenesis is limited. In this study, we conducted a systematic interrogation of MAEs in bone marrow mononuclear cells (BMMCs) at single-cell resolution to construct a MAE atlas of BMMCs. We identified 1,020 constitutive MAEs in BMMCs, which included imprinted genes such as MEG8, NAP1L5, and IRAIN. We classified the BMMCs into six cell types and identified 74 cell type specific MAEs including MTSS1, MOB1A, and TCF12. We further identified 114 random MAEs (rMAEs) at single-cell level, with 78.1% single-allele rMAE and 21.9% biallelic mosaic rMAE. Many MAEs identified in BMMCs have not been reported and are potentially hematopoietic specific, supporting MAEs are functional relevance. Comparison of BMMC samples from a leukemia patient with multiple clinical stages showed the fractions of constitutive MAE were correlated with fractions of leukemia cells in BMMCs. Further separation of the BMMCs into leukemia cells and normal cells showed that leukemia cells have much higher constitutive MAE and rMAEs than normal cells. We identified the leukemia cell-specific MAEs and relapsed leukemia cell-specific MAEs, which were enriched in immune-related functions. These results indicate MAE is prevalent and is an important gene regulation mechanism during hematopoiesis and leukemogenesis. As the first systematical interrogation of constitutive MAEs, cell type specific MAEs, and rMAEs during hematopoiesis and leukemogenesis, the study significantly increased our knowledge about the features and functions of MAEs.

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.

Semicoordinated allelic-bursting shape dynamic random monoallelic expression in pregastrulation embryos

Recently, allele-specific single-cell RNA-seq analysis has demonstrated widespread dynamic random monoallelic expression of autosomal genes (aRME) in different cell types. However, the prevalence of dynamic aRME during pregastrulation remains unknown. Here, we show that dynamic aRME is widespread in different lineages of pregastrulation embryos. Additionally, the origin of dynamic aRME remains elusive. It is believed that independent transcriptional bursting from each allele leads to dynamic aRME. Here, we show that allelic burst is not perfectly independent; instead it happens in a semicoordinated fashion. Importantly, we show that semicoordinated allelic bursting of genes, particularly with low burst frequency, leads to frequent asynchronous allelic bursting, thereby contributing to dynamic aRME. Furthermore, we found that coordination of allelic bursting is lineage specific and genes regulating the development have a higher degree of coordination. Altogether, our study provides significant insights into the prevalence and origin of dynamic aRME and their developmental relevance during early development.

•

Dynamic aRME is widespread in different lineages of pregastrulation embryos

•

Semicoordinated bursting of genes with low burst frequency leads to dynamic aRME

•

Degree of coordination of allelic bursting is lineage specific

•

Developmental genes have higher degree of coordination of allelic bursting

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

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.

Increased antitumor efficacy of PD-1-deficient melanoma-specific human lymphocytes

Background

Genome editing offers unique perspectives for optimizing the functional properties of T cells for adoptive cell transfer purposes. So far, PDCD1 editing has been successfully tested mainly in chimeric antigen receptor T (CAR-T) cells and human primary T cells. Nonetheless, for patients with solid tumors, the adoptive transfer of effector memory T cells specific for tumor antigens remains a relevant option, and the use of high avidity T cells deficient for programmed cell death-1 (PD-1) expression is susceptible to improve the therapeutic benefit of these treatments.

Methods

Here we used the transfection of CAS9/sgRNA ribonucleoproteic complexes to edit PDCD1 gene in human effector memory CD8⁺ T cells specific for the melanoma antigen Melan-A. We cloned edited T cell populations and validated PDCD1 editing through sequencing and cytometry in each T cell clone, together with T-cell receptor (TCR) chain’s sequencing. We also performed whole transcriptomic analyses on wild-type (WT) and edited T cell clones. Finally, we documented in vitro and in vivo through adoptive transfer in NOD scid gamma (NSG) mice, the antitumor properties of WT and PD-1KO T cell clones, expressing the same TCR.

Results

Here we demonstrated the feasibility to edit PDCD1 gene in human effector memory melanoma-specific T lymphocytes. We showed that PD-1 expression was dramatically reduced or totally absent on PDCD1-edited T cell clones. Extensive characterization of a panel of T cell clones expressing the same TCR and exhibiting similar functional avidity demonstrated superior antitumor reactivity against a PD-L1 expressing melanoma cell line. Transcriptomic analysis revealed a downregulation of genes involved in proliferation and DNA replication in PD-1-deficient T cell clones, whereas genes involved in metabolism and cell signaling were upregulated. Finally, we documented the superior ability of PD-1-deficient T cells to significantly delay the growth of a PD-L1 expressing human melanoma tumor in an NSG mouse model.

Conclusion

The use of such lymphocytes for adoptive cell transfer purposes, associated with other approaches modulating the tumor microenvironment, would be a promising alternative to improve immunotherapy efficacy in solid tumors.

Assessment of escape from X chromosome inactivation and gene expression in single human immune cells

X-chromosomal genes escaping from X chromosome inactivation (XCI) in immune cells can contribute to sex-specific differences in immune responses. This protocol describes the specific steps to determine escape from XCI and to simultaneously quantify mRNA expression of multiple genes at the single immune cell level using a single-nucleotide polymorphism approach. The protocol furthermore allows the analysis of allele-specific expression of X-chromosomal genes.

For complete details on the use and execution of this protocol, please refer to Hagen et al. (2020).

•

Approach to investigate escape from XCI and gene expression in single cells

•

Simultaneous gene expression measurement of over 100 genes in one cell

•

Assessment of allele-specific expression of genes with monoallelic expression pattern

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.