Allele-specific Gene Expression Research Articles

Characterizing the allele-specific gene expression landscape in high hyperdiploid acute lymphoblastic leukemia with BASE.

Somatic copy number variations (CNVs), including abnormal chromosome numbers and structural changes leading to gain or loss of genetic material, play a crucial role in initiation and progression of cancer. CNVs are believed to cause gene dosage imbalances and modify cis-regulatory elements, leading to allelic expression imbalances in genes that influence cell division and thereby contribute to cancer development. However, the impact of CNVs on allelic gene expression in cancer remains unclear. Allele-specific expression (ASE) analysis, a potent method for investigating genome-wide allelic imbalance profiles in tumors, assesses the relative expression of two alleles using high-throughput sequencing data. However, many existing methods for gene-level ASE detection rely on only RNA sequencing data, which present challenges in interpreting the genetic mechanisms underlying ASE in cancer. To address this issue, we developed a robust framework that integrates allele-specific copy number calls into ASE calling algorithms by leveraging paired genome and transcriptome data from the same sample. This integration enhances the interpretability of the genetic mechanisms driving ASE, thereby facilitating the identification of driver events triggered by CNVs in cancer. In this study, we utilized BASE to conduct a comprehensive analysis of ASE in high hyperdiploid acute lymphoblastic leukemia (HeH ALL), a prevalent childhood malignancy characterized by gains of chromosomes X, 4, 6, 10, 14, 17, 18, and 21. Our analysis unveiled the comprehensive ASE landscape in HeH ALL. Through a multi-perspective examination of HeH ASEs, we offer a systematic understanding of how CNVs impact ASE in HeH, providing valuable insights to guide ASE studies in cancer.

Maintenance of X chromosome inactivation after T cell activation requires NF-κB signaling.

X chromosome inactivation (XCI) balances X-linked gene dosage between sexes. Unstimulated T cells lack cytological enrichment of X-inactive specific transcript (Xist) RNA and heterochromatic modifications on the inactive X chromosome (Xi), which are involved in maintenance of XCI, and these modifications return to the Xi after stimulation. Here, we examined allele-specific gene expression and epigenomic profiles of the Xi in T cells. We found that the Xi in unstimulated T cells is largely dosage compensated and enriched with the repressive H3K27me3 modification but not the H2AK119-ubiquitin (Ub) mark. Upon T cell stimulation mediated by both CD3 and CD28, the Xi accumulated H2AK119-Ub at gene regions of previous H3K27me3 enrichment. T cell receptor (TCR) engagement, specifically NF-κB signaling downstream of the TCR, was required for Xist RNA localization to the Xi. Disruption of NF-κB signaling in mouse and human T cells using genetic deletion, chemical inhibitors, and patients with immunodeficiencies prevented Xist/XIST RNA accumulation at the Xi and altered X-linked gene expression. Our findings reveal a previously undescribed connection between NF-κB signaling pathways, which affects XCI maintenance in T cells in females.

Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs.

Transcription factors (TFs) regulate gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because TF occupancy is driven in part by recognition of DNA sequence, genetic variation can influence TF-DNA associations and gene regulation. To identify variants that impact TF binding in human brain tissues, we assessed allele-specific binding (ASB) at heterozygous variants for 94 TFs in nine brain regions from two donors. Leveraging graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signals between alleles at heterozygous variants within each brain region and identified thousands of variants exhibiting ASB for at least one TF. ASB reproducibility was measured by comparisons between independent experiments both within and between donors. We found that rare alleles in the general population more frequently led to reduced TF binding, whereas common alleles had an equal likelihood of increasing or decreasing binding. Further, for ASB variants in predicted binding motifs, the favored allele tended to be the one with the stronger expected motif match, but this concordance was not observed within highly occupied sites. We also found that neuron-specific cis-regulatory elements (cCREs), in contrast with oligodendrocyte-specific cCREs, showed depletion of ASB variants. We identified 2670 ASB variants associated with evidence for allele-specific gene expression in the brain from GTEx data and observed increasing eQTL effect direction concordance as ASB significance increases. These results provide a valuable and unique resource for mechanistic analysis of cis-regulatory variation in human brain tissue.

PSVII-1 Widespread allele-specific expression across tissues and cells of the chicken genome

Abstract Chickens are a valuable livestock species, providing an affordable and nutritious food source worldwide through eggs and meat. Additionally, chickens are extensively used in scientific research as model organisms in virology, immunology, epigenetics, development, and conservation biology. Considering this, there are increasing efforts to deepen our understanding of their genome complexity, such as efforts by the Functional Annotation of Animal Genomes (FAANG) Consortium. Toward this goal, we are working to contribute to the annotation of the chicken transcriptome. Allele-specific expression (ASE) is an imbalance of allelic gene expression often driven by genetic and epigenetic changes in cis-regulatory regions. We, therefore, aimed to determine allelic imbalances in gene expression across 20 tissues and cells to shed light on mechanisms regulating gene expression. Using replicate paired-end (PE) RNA-seq samples, we identified around 7 million variants across 20 tissues and cells (e.g., reproductive tissues, intestine tissues, muscle, and immune tissues and cells). We applied quality control measures and filtered out monoallelic and homozygous variants. After that, we calculated read counts per allele to determine allelic expression. We found 365,894 significant ASE SNPs and 11,530 ASE genes in at least one tissue or cell type (FDR ≤ 0.05). We discovered that ASE SNPs are widely distributed throughout the chicken genome, and ASE SNPs and ASE genes showed a varied distribution across tissues and cells. Primary macrophage exhibited the most abundant, while the pectoralis major (breast muscle) had the lowest ASE genes and ASE SNPs, ranging from 7,592 to 32,505 ASE SNPs and 773 to 2,387 ASE genes. Our findings reveal several important pathways affected by allelic imbalance of expression. These pathways included fatty acid biosynthesis and regulation, cell proliferation and differentiation, cell metabolism, adipocytokine signaling, proteolysis and autophagy, and focal adhesion. In summary, our study provides insight into relevant biological processes critical to chickens and can contribute to improving genome annotation, helping to expand the transcriptomic reference source.

Allele-specific DNA demethylation editing leads to stable upregulation of allele-specific gene expression

Epigenome editing is an emerging technology that allows to rewrite epigenome states and reprogram gene expression. Here, we have developed allele-specific DNA demethylation editing at gene promoters containing an SNP by sgRNA/dCas9 mediated recuitment of TET1. Maximal DNA demethylation (up to 90%) was observed 6days after transient transfection of the epigenome editors and it was almost stable for 15days. After allele-specific targeting, DNA demethylation was up to 15-fold more efficient at the targeted allele. Our data show that locus-specific and allele-specific DNA demethylation can trigger the expression of previously silenced genes. Allele-specific DNA demethylation shifted allelic expression ratios about 4-fold. Allele-specific DNA demethylation could be used to correct aberrant imprinting in patients suffering from imprinting disorders and to study the roles of individual alleles of a gene in a given cellular context.

Biparental graph strategy to represent and analyze hybrid plant genomes.

Hybrid plants are found extensively in the wild, and they often demonstrate superior performance of complex traits over their parents and other selfing plants. This phenomenon, known as heterosis, has been extensively applied in plant breeding for decades. However, the process of decoding hybrid plant genomes has seriously lagged due to the challenges associated with genome assembly and the lack of appropriate methodologies for their subsequent representation and analysis. Here, we present the assembly and analysis of 2 hybrids, an intraspecific hybrid between 2 maize (Zea mays ssp. mays) inbred lines and an interspecific hybrid between maize and its wild relative teosinte (Z. mays ssp. parviglumis), utilizing a combination of PacBio High Fidelity sequencing and chromatin conformation capture sequencing data. The haplotypic assemblies are well phased at chromosomal scale, successfully resolving the complex loci with extensive parental structural variations (SVs). By integrating into a biparental genome graph, the haplotypic assemblies can facilitate downstream short-read-based SV calling and allele-specific gene expression analysis, demonstrating outstanding advantages over a single linear genome. Our work offers a comprehensive workflow that aims to facilitate the decoding of numerous hybrid plant genomes, particularly those with unknown or inaccessible parentage, thereby enhancing our understanding of genome evolution and heterosis.

Molecular mechanisms underlying gene regulatory variation of maize metabolic traits.

Variation in gene expression levels is pervasive among individuals and races or varieties, and has substantial agronomic consequences, for example, by contributing to hybrid vigor. Gene expression level variation results from mutations in regulatory sequences (cis) and/or transcription factor (TF) activity (trans), but the mechanisms underlying cis- and/or trans-regulatory variation of complex phenotypes remain largely unknown. Here, we investigated gene expression variation mechanisms underlying the differential accumulation of the insecticidal compounds maysin and chlorogenic acid in silks of widely used maize (Zea mays) inbreds, B73 and A632. By combining transcriptomics and cistromics, we identified 1,338 silk direct targets of the maize R2R3-MYB TF Pericarp color1 (P1), consistent with it being a regulator of maysin and chlorogenic acid biosynthesis. Among these P1 targets, 464 showed allele-specific expression (ASE) between B73 and A632 silks. Allelic DNA-affinity purification sequencing identified 34 examples in which P1 allelic specific binding (ASB) correlated with cis-expression variation. From previous yeast one-hybrid studies, we identified 9 TFs potentially implicated in the control of P1 targets, with ASB to 83 out of 464 ASE genes (cis) and differential expression of 4 out of 9 TFs between B73 and A632 silks (trans). These results provide a molecular framework for understanding universal mechanisms underlying natural variation of gene expression levels, and how the regulation of metabolic diversity is established.

Haplotype-resolved T2T genome assemblies and pangenome graph of pear reveal diverse patterns of allele-specific expression and the genomic basis of fruit quality traits

Hybrid crops often exhibit increased yield and greater resilience, yet the genomic mechanism(s) underlying hybrid vigor or heterosis remain unclear, hindering our ability to predict the expression of phenotypic traits in hybrid breeding. Here, we generated haplotype-resolved T2T genome assemblies of two pear hybrid varieties, ‘Yuluxiang’ (YLX) and ‘Hongxiangsu’ (HXS), which share the same maternal parent but differ in their paternal parents. We then used these assemblies to explore the genome-scale landscape of allele-specific expression (ASE) and create a pangenome graph for pear. ASE was observed for close to 6000 genes in both hybrid cultivars. A subset of ASE genes related to aspects of fruit quality such as sugars, organic acids, and cuticular wax were identified, suggesting their important contributions to heterosis. Specifically, Ma1, a gene regulating fruit acidity, is absent in the paternal haplotypes of HXS and YLX. A pangenome graph was built based on our assemblies and seven published pear genomes. Resequencing data for 139 cultivated pear genotypes (including 97 genotypes sequenced here) were subsequently aligned to the pangenome graph, revealing numerous structural variant hotspots and selective sweeps during pear diversification. As predicted, the Ma1 allele was found to be absent in varieties with low organic acid content, and this association was functionally validated by Ma1 overexpression in pear fruit and calli. Overall, these results reveal the contributions of ASE to fruit-quality heterosis and provide a robust pangenome reference for high-resolution allele discovery and association mapping.

Differential 3D genome architecture and imprinted gene expression: cause or consequence?

Imprinted genes provide an attractive paradigm to unravel links between transcription and genome architecture. The parental allele-specific expression of these essential genes - which are clustered in chromosomal domains - is mediated by parental methylation imprints at key regulatory DNA sequences. Recent chromatin conformation capture (3C)-based studies show differential organization of topologically associating domains between the parental chromosomes at imprinted domains, in embryonic stem and differentiated cells. At several imprinted domains, differentially methylated regions show allelic binding of the insulator protein CTCF, and linked focal retention of cohesin, at the non-methylated allele only. This generates differential patterns of chromatin looping between the parental chromosomes, already in the early embryo, and thereby facilitates the allelic gene expression. Recent research evokes also the opposite scenario, in which allelic transcription contributes to the differential genome organization, similarly as reported for imprinted X chromosome inactivation. This may occur through epigenetic effects on CTCF binding, through structural effects of RNA Polymerase II, or through imprinted long non-coding RNAs that have chromatin repressive functions. The emerging picture is that epigenetically-controlled differential genome architecture precedes and facilitates imprinted gene expression during development, and that at some domains, conversely, the mono-allelic gene expression also influences genome architecture.

Natural variations of heterosis-related allele-specific expression genes in promoter regions lead to allele-specific expression in maize

BackgroundHeterosis has successfully enhanced maize productivity and quality. Although significant progress has been made in delineating the genetic basis of heterosis, the molecular mechanisms underlying its genetic components remain less explored. Allele-specific expression (ASE), the imbalanced expression between two parental alleles in hybrids, is increasingly being recognized as a factor contributing to heterosis. ASE is a complex process regulated by both epigenetic and genetic variations in response to developmental and environmental conditions.ResultsIn this study, we explored the differential characteristics of ASE by analyzing the transcriptome data of two maize hybrids and their parents under four light conditions. On the basis of allele expression patterns in different hybrids under various conditions, ASE genes were divided into three categories: bias-consistent genes involved in basal metabolic processes in a functionally complementary manner, bias-reversal genes adapting to the light environment, and bias-specific genes maintaining cell homeostasis. We observed that 758 ASE genes (ASEGs) were significantly overlapped with heterosis quantitative trait loci (QTLs), and high-frequency variations in the promoter regions of heterosis-related ASEGs were identified between parents. In addition, 10 heterosis-related ASEGs participating in yield heterosis were selected during domestication.ConclusionsThe comprehensive analysis of ASEGs offers a distinctive perspective on how light quality influences gene expression patterns and gene-environment interactions, with implications for the identification of heterosis-related ASEGs to enhance maize yield.

A phased genome of the highly heterozygous 'Texas' almond uncovers patterns of allele-specific expression linked to heterozygous structural variants.

The vast majority of traditional almond varieties are self-incompatible, and the level of variability of the species is very high, resulting in a high-heterozygosity genome. Therefore, information on the different haplotypes is particularly relevant to understand the genetic basis of trait variability in this species. However, although reference genomes for several almond varieties exist, none of them is phased and has genome information at the haplotype level. Here, we present a phased assembly of genome of the almond cv. Texas. This new assembly has 13% more assembled sequence than the previous version of the Texas genome and has an increased contiguity, in particular in repetitive regions such as the centromeres. Our analysis shows that the 'Texas' genome has a high degree of heterozygosity, both at SNPs, short indels, and structural variants level. Many of the SVs are the result of heterozygous transposable element insertions, and in many cases, they also contain genic sequences. In addition to the direct consequences of this genic variability on the presence/absence of genes, our results show that variants located close to genes are often associated with allele-specific gene expression, which highlights the importance of heterozygous SVs in almond.

The OsWRKY72-OsAAT30/OsGSTU26 module mediates reactive oxygen species scavenging to drive heterosis for salt tolerance in hybrid rice.

Hybrid rice (Oryza sativa) generally outperforms its inbred parents in yield and stress tolerance, a phenomenon termed heterosis, but the underlying mechanism is not completely understood. Here, we combined transcriptome, proteome, physiological, and heterosis analyses to examine the salt response of super hybrid rice Chaoyou1000 (CY1000). In addition to surpassing the mean values for its two parents (mid-parent heterosis), CY1000 exhibited a higher reactive oxygen species scavenging ability than both its parents (over-parent heterosis or heterobeltiosis). Nonadditive expression and allele-specific gene expression assays showed that the glutathione S-transferase gene OsGSTU26 and the amino acid transporter gene OsAAT30 may have major roles in heterosis for salt tolerance, acting in an overdominant fashion in CY1000. Furthermore, we identified OsWRKY72 as a common transcription factor that binds and regulates OsGSTU26 and OsAAT30. The salt-sensitive phenotypes were associated with the OsWRKY72paternal genotype or the OsAAT30maternal genotype in core rice germplasm varieties. OsWRKY72paternal specifically repressed the expression of OsGSTU26 under salt stress, leading to salinity sensitivity, while OsWRKY72maternal specifically repressed OsAAT30, resulting in salinity tolerance. These results suggest that the OsWRKY72-OsAAT30/OsGSTU26 module may play an important role in heterosis for salt tolerance in an overdominant fashion in CY1000 hybrid rice, providing valuable clues to elucidate the mechanism of heterosis for salinity tolerance in hybrid rice.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq

BackgroundParent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference.ResultsWe used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair.ConclusionsFor the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

AStruct: detection of allele-specific RNA secondary structure in structuromic probing data

BackgroundUncovering functional genetic variants from an allele-specific perspective is of paramount importance in advancing our understanding of gene regulation and genetic diseases. Recently, various allele-specific events, such as allele-specific gene expression, allele-specific methylation, and allele-specific binding, have been explored on a genome-wide scale due to the development of high-throughput sequencing methods. RNA secondary structure, which plays a crucial role in multiple RNA-associated processes like RNA modification, translation and splicing, has emerged as an essential focus of relevant research. However, tools to identify genetic variants associated with allele-specific RNA secondary structures are still lacking.ResultsHere, we develop a computational tool called ‘AStruct’ that enables us to detect allele-specific RNA secondary structure (ASRS) from RT-stop based structuromic probing data. AStruct shows robust performance in both simulated datasets and public icSHAPE datasets. We reveal that single nucleotide polymorphisms (SNPs) with higher AStruct scores are enriched in coding regions and tend to be functional. These SNPs are highly conservative, have the potential to disrupt sites involved in m6A modification or protein binding, and are frequently associated with disease.ConclusionsAStruct is a tool dedicated to invoke allele-specific RNA secondary structure events at heterozygous SNPs in RT-stop based structuromic probing data. It utilizes allelic variants, base pairing and RT-stop information under different cell conditions to detect dynamic and functional ASRS. Compared to sequence-based tools, AStruct considers dynamic cell conditions and outperforms in detecting functional variants. AStruct is implemented in JAVA and is freely accessible at: https://github.com/canceromics/AStruct.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome.

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

Allele-specific DNA methylation and gene expression during shoot organogenesis in tissue culture of hybrid poplar.

Plant tissue regeneration is critical for genetic transformation and genome editing techniques. During the regeneration process, changes in epigenetic modifications accompany the cell fate transition. However, how allele-specific DNA methylation in two haplotypes contributes to the transcriptional dynamics during regeneration remains elusive. Here we applied an inter-species hybrid poplar (Populus alba × P. glandulosa cv. 84K) as a system to characterize the DNA methylation landscape during de novo shoot organogenesis at allele level. Both direct and indirect shoot organogenesis showed a reduction in genome-wide DNA methylation. At gene level, non-expressed genes were hypermethylated in comparison with expressed genes. Among the genes exhibiting significant correlations between levels of DNA methylation and gene expression, the expression patterns of 75% of genes were negatively correlated with DNA methylation in the CG context, whereas the correlation patterns in the CHH context were the reverse. The allele-biased DNA methylation was consistent during shoot organogenesis, with fewer than one-thousandth of allele-specific methylation regions shifted. Analysis of allele-specific expression revealed that there were only 1909 genes showing phase-dependent allele-biased expression in the regeneration process, among which the allele pairs with greater differences in transcription factor binding sites at promoter regions exhibited greater differences in allele expression. Our results indicated a relatively independent transcriptional regulation in two subgenomes during shoot organogenesis, which was contributed by cis-acting genomic and epigenomic variations.

Regulatory features aid interpretation of 3′UTR variants

Our ability to determine the clinical impact of variants in 3' untranslated regions (UTRs) of genes remains poor. We provide a thorough analysis of 3' UTR variants from several datasets. Variants in putative regulatory elements, including RNA-binding protein motifs, eCLIP peaks, and microRNA sites, are up to 16 times more likely than variants not in these elements to have gene expression and phenotype associations. Variants in regulatory motifs result in allele-specific protein binding in cell lines and allele-specific gene expression differences in population studies. In addition, variants in shared regions of alternatively polyadenylated isoforms and those proximal to polyA sites are more likely to affect gene expression and phenotype. Finally, pathogenic 3' UTR variants in ClinVar are up to 20 times more likely than benign variants to fall in a regulatory site. We incorporated these findings into RegVar, a software tool that interprets regulatory elements and annotations for any 3' UTR variant and predicts whether the variant is likely to affect gene expression or phenotype. This tool will help prioritize variants for experimental studies and identify pathogenic variants in individuals.

Inter-subspecies mouse F1 hybrid embryonic stem cell lines newly established for studies of allelic imbalance in gene expression.

Allele-specific monoallelic gene expression is a unique phenomenon and a great resource for analyzing gene regulation. To study this phenomenon, we established new embryonic stem (ES) cell lines derived from F1 hybrid blastocysts from crosses between four mouse subspecies (Mus musculus domesticus, C57BL/6; M. musculus molossinus, MSM/Ms; M. musculus musculus, PWK; M. musculus castaneus, HMI/Ms) and analyzed the expression levels of undifferentiated pluripotent stem cell markers and karyotypes of each line. To demonstrate the utility of our cell lines, we analyzed the allele-specific expression pattern of the Inpp5d gene as an example. The allelic expression depended on the parental alleles; this dependence could be a consequence of differences in compatibility between cis- and trans-elements of the Inpp5d gene from different subspecies. The use of parental mice from four subspecies greatly enhanced genetic polymorphism. The F1 hybrid ES cells retained this polymorphism not only in the Inpp5d gene, but also at a genome-wide level. As we demonstrated for the Inpp5d gene, the established cell lines can contribute to the analysis of allelic expression imbalance based on the incompatibility between cis- and trans-elements and of phenotypes related to this incompatibility.