Allele-specific Expression Studies Research Articles

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans.

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

GPNMB confers risk for Parkinson's disease through interaction with α-synuclein.

Many risk loci for Parkinson's disease (PD) have been identified by genome-wide association studies (GWASs), but target genes and mechanisms remain largely unknown. We linked the GWAS-derived chromosome 7 locus (sentinel single-nucleotide polymorphism rs199347) to GPNMB through colocalization analyses of expression quantitative trait locus and PD risk signals, confirmed by allele-specific expression studies in the human brain. In cells, glycoprotein nonmetastatic melanoma protein B (GPNMB) coimmunoprecipitated and colocalized with α-synuclein (aSyn). In induced pluripotent stem cell-derived neurons, loss of GPNMB resulted in loss of ability to internalize aSyn fibrils and develop aSyn pathology. In 731 PD and 59 control biosamples, GPNMB was elevated in PD plasma, associating with disease severity. Thus, GPNMB represents a PD risk gene with potential for biomarker development and therapeutic targeting.

AScan: A Novel Method for the Study of Allele Specific Expression in Single Individuals

In diploid organisms, two copies of each allele are normally inherited from parents. Paternal and maternal alleles can be regulated and expressed unequally, which is referred to as allele-specific expression (ASE). In this work, we present aScan, a novel method for the identification of ASE from the analysis of matched individual genomic and RNA sequencing data. By performing extensive analyses of both real and simulated data, we demonstrate that aScan can correctly identify ASE with high accuracy and sensitivity in different experimental settings. Additionally, by applying our method to a small cohort of individuals that are not included in publicly available databases of human genetic variation, we outline the value of possible applications of ASE analysis in single individuals for deriving a more accurate annotation of “private” low-frequency genetic variants associated with regulatory effects on transcription. All in all, we believe that aScan will represent a beneficial addition to the set of bioinformatics tools for the analysis of ASE. Finally, while our method was initially conceived for the analysis of RNA-seq data, it can in principle be applied to any quantitative NGS assay for which matched genotypic and expression data are available. AvailabilityaScan is currently available in the form of an open source standalone software package at: https://github.com/Federico77z/aScan/. aScan version 1.0.3, available at https://github.com/Federico77z/aScan/releases/tag/1.0.3, has been used for all the analyses included in this manuscript. A Docker image of the tool has also been made available at https://github.com/pmandreoli/aScanDocker.

The contribution of cis- and trans-acting variants to gene regulation in wild and domesticated barley under cold stress and control conditions.

Barley, like other crops, has experienced a series of genetic changes that have impacted its architecture and growth habit to suit the needs of humans, termed the domestication syndrome. Domestication also resulted in a concomitant bottleneck that reduced sequence diversity in genes and regulatory regions. Little is known about regulatory changes resulting from domestication in barley. We used RNA sequencing to examine allele-specific expression in hybrids between wild and domesticated barley. Our results show that most genes have conserved regulation. In contrast to studies of allele-specific expression in interspecific hybrids, we find almost a complete absence of trans effects. We also find that cis regulation is largely stable in response to short-term cold stress. Our study has practical implications for crop improvement using wild relatives. Genes regulated in cis are more likely to be expressed in a new genetic background at the same level as in their native background.

A Genome-Wide Study of Allele-Specific Expression in Colorectal Cancer.

Accumulating evidence from small-scale studies has suggested that allele-specific expression (ASE) plays an important role in tumor initiation and progression. However, little is known about genome-wide ASE in tumors. In this study, we conducted a comprehensive analysis of ASE in individuals with colorectal cancer (CRC) on a genome-wide scale. We identified 5.4 thousand genome-wide ASEs of single nucleotide variations (SNVs) from tumor and normal tissues of 59 individuals with CRC. We observed an increased ASE level in tumor samples and the ASEs enriched as hotspots on the genome. Around 63% of the genes located there were previously reported to contain complex regulatory elements, e.g., human leukocyte antigen (HLA), or were implicated in tumor progression. Focussing on the allelic expression of somatic mutations, we found that 37.5% of them exhibited ASE, and genes harboring such somatic mutations, were enriched in important pathways implicated in cancers. In addition, by comparing the expected and observed ASE events in tumor samples, we identified 50 tumor specific ASEs which possibly contributed to the somatic events in the regulatory regions of the genes and significantly enriched known cancer driver genes. By analyzing CRC ASEs from several perspectives, we provided a systematic understanding of how ASE is implicated in both tumor and normal tissues and will be of critical value in guiding ASE studies in cancer.

Cis-regulatory evolution in prokaryotes revealed by interspecific archaeal hybrids

The study of allele-specific expression (ASE) in interspecific hybrids has played a central role in our understanding of a wide range of phenomena, including genomic imprinting, X-chromosome inactivation, and cis-regulatory evolution. However across the hundreds of studies of hybrid ASE, all have been restricted to sexually reproducing eukaryotes, leaving a major gap in our understanding of the genomic patterns of cis-regulatory evolution in prokaryotes. Here we introduce a method to generate stable hybrids between two species of halophilic archaea, and measure genome-wide ASE in these hybrids with RNA-seq. We found that over half of all genes have significant ASE, and that genes encoding kinases show evidence of lineage-specific selection on their cis-regulation. This pattern of polygenic selection suggested species-specific adaptation to low phosphate conditions, which we confirmed with growth experiments. Altogether, our work extends the study of ASE to archaea, and suggests that cis-regulation can evolve under polygenic lineage-specific selection in prokaryotes.

RNA-Sequencing: A tool to explore new frontiers in animal genetics

Global gene expression profiling by RNA-Sequencing (RNA-Seq) was introduced in 2004 with the advent of next generation sequencing technologies and has made a revolutionary impact on transcriptome analysis. When compared to the traditional microarray technologies, RNA-Seq has major advantages such as a higher dynamic range, a precise unbiased quantification of transcripts and a comprehensive coverage of all expressed sequences in a given tissue. Besides transcript quantification, RNA-Seq is utilized for comparison of transcriptomes, mapping the 5′ and 3′ boundaries of mRNA and introns, finding of novel exons, study of allele-specific expression, SNP discovery, study of non-coding RNA and de novo gene model building. Two major challenges faced with the RNA-Seq are the differences in sequencing technologies and bioinformatics analysis. Recently, RNA-Seq technology has increasingly been used in livestock, wild and aquatic animal studies to measure gene expression analysis and SNP discovery. This review summarizes the capabilities and challenges of RNA-Seq and the impact of this technology on animal genetics research.

Quantifying RNA allelic ratios by microfluidic multiplex PCR and sequencing

We developed a targeted RNA sequencing method that couples microfluidics-based multiplex PCR and deep sequencing (mmPCR-seq) to uniformly and simultaneously amplify up to 960 loci in 48 samples independently of their gene expression levels, and accurately and cost-effectively measure allelic ratios even for low-quantity or low-quality RNA samples. We applied mmPCR-seq to RNA editing and allele-specific expression studies. mmPCR-seq complements RNA-seq and provides a highly desirable solution for future applications.

The 14q22.2 colorectal cancer variant rs4444235 shows cis-acting regulation of BMP4

Common genetic variation at human 14q22.2 tagged by rs4444235 is significantly associated with colorectal cancer (CRC) risk. Re-sequencing was used to comprehensively annotate the 17kb region of strong linkage disequilibrium encompassing rs4444235. Through bioinformatic analyses using H3K4Me1, H3K4Me3, and DNase-I hypersensitivity chromatin signatures and evolutionary conservation we identified seven candidate disease-causing single-nucleotide polymorphisms mapping to six regions within the 17-kb region predicted to have regulatory potential. Reporter gene studies of these regions demonstrated that the element to which rs4444235 maps acts as an allele-specific transcriptional enhancer. Allele-specific expression studies in CRC cell lines heterozygous for rs4444235 showed significantly increased expression of bone morphogenetic protein-4 (BMP4) associated with the risk allele (P<0.001). These data provide evidence for a functional basis for the non-coding risk variant rs4444235 at 14q22.2 and emphasizes the importance of genetic variation in the BMP pathway genes as determinants of CRC risk.

The genetic basis of evolutionary change in gene expression levels

The regulation of gene expression is an important determinant of organismal phenotype and evolution. However, the widespread recognition of this fact occurred long after the synthesis of evolution and genetics. Here, we give a brief sketch of thoughts regarding gene regulation in the history of evolution and genetics. We then review the development of genome-wide studies of gene regulatory variation in the context of the location and mode of action of the causative genetic changes. In particular, we review mapping of the genetic basis of expression variation through expression quantitative trait locus studies and measuring the cis/trans component of expression variation in allele-specific expression studies. We conclude by proposing a systematic integration of ideas that combines global mapping studies, cis/trans tests and modern population genetics methodologies, in order to directly estimate the forces acting on regulatory variation within and between species.

Influence of human genome polymorphism on gene expression

Genetic variation, through its effects on gene expression, plays a crucial role in phenotypic variation and disease susceptibility. Recent studies from our group and others have integrated a number of resources and technologies to assess several aspects of genome variation affecting gene expression. Some of these large-scale mapping studies involving expression quantitative traits have recently been reviewed [Gibson, G. and Weir, B. (2005) The quantitative genetics of transcription. Trends Genet., 21, 616-623; de Koning, D.J. and Haley, C.S. (2005) Genetical genomics in humans and model organisms. Trends Genet., 21, 377-381], with particular attention to the statistical issues. In this review, we compare allele-specific expression studies in human samples (primarily lymphoblastoid cell lines from the CEPH HapMap panel), as a prelude to a discussion on study design issues and sources of variation, in order to propose the steps required to build a detailed map of cis-acting regulatory variation in the human genome. Obtaining panels of tissues from large numbers of individuals remains an important limitation. We also conclude that there is insufficient knowledge as to the feasibility of comprehensive studies of trans-acting variation in the human genome.

Allele-specific gene expression differences in humans.

In the last decade, the search for the genetic origins of phenotypic variation has expanded beyond the non-synonymous variants which alter the amino acid sequence of the encoded protein, and many examples of sequence variants which alter gene expression have been found. Recently, using both traditional and novel technologies, a number of surveys have been carried out to examine the frequency with which cis-acting sequence variants or other cis-acting effects, alter gene expression either in vitro or in vivo. Microarray data have shown that the expression of many genes varies markedly between individuals and allele-specific expression studies have shown that the source of much of this variation appears to be cis-acting effects. A significant proportion of the variation may originate in gene promoter regions and a large number of sequence variants which have functional effect in vitro have been found. The evidence suggests that given a large enough population, most, if not all genes may have allele-specific expression differences in at least some individuals and finding the genetic origins of each of these and linking the former to a possible phenotype must be a major long term goal of the biomedical community.

Allele-specific quantification of HLA-DQB1 gene expression by real-time reverse transcriptase-polymerase chain reaction.

In addition to coding region polymorphism, allele-specific variation in the upstream regulatory region of the HLA-DQB1 gene has been detected. Reporter gene assays and transfection studies have indicated that HLA-DQB1 promoter polymorphism may be of functional significance. The aim of this study was to utilize real-time reverse transcriptase-polymerase chain reaction (RT-PCR) for allele-specific quantification of HLA-DQB1 expression and to analyze cell-specific HLA-DQB1 expression in vivo. For the allele-specific quantification of DQB1 gene products, a real-time RT-PCR set of primer pairs (n=27) and probes (n=5) targeting exon 2 variability was established. The robustness and integrity of the assay system were confirmed by using recombinant DQB1 exon 2 plasmid clones as active exogenous controls. Sensitivity and reproducibility were assessed by serial dilution and allelic mixing analyses. In application to the study of allele-specific expression of DQB1 gene products during cytokine-driven maturation of monocyte-derived dendritic cells, differential patterns of allelic expression in heterozygous individuals were observed for DQB1*0301, compared to DQB1*0501 and DQB1*0602. At maximum, 1.9-fold (*0301/*0501) and 2.5-fold (*0301/*0602) higher induction was seen for DQB*0301. In conclusion, HLA-DQB1 expression can be analyzed by real-time RT-PCR suitable for cell- and allele-specific detection of HLA-DQB1 transcripts in homo- and heterozygous combinations.

Allele-specific gene expression uncovered

Genetic variation in populations can result in variation in levels of gene expression but the extent to which this occurs has been unclear. In this article, recent studies of allele-specific expression among autosomal non-imprinted genes are reviewed. These new data provide evidence that differential expression is relatively common and that allelic differences are heritable and can be highly context specific.

Uncoupling protein-1 mRNA expression in obese human subjects: the role of sequence variations at the uncoupling protein-1 gene locus

Uncoupling protein-1 (UCP-1) activity in brown adipose tissue increases thermogenesis, contributes to facultative energy expenditure in humans, and has been implicated in the pathogenesis of rodent obesity. To determine genetic factors controlling UCP-1 expression in humans, we measured intra- and extraperitoneal UCP-1 mRNA abundance levels by a competitive RT-PCR method and compared expression levels with common sequence variations in the β3-adrenergic receptor gene and the distal UCP-1 gene promoter in obese human subjects. While median and average UCP-1 mRNA levels in both the intra- and extraperitoneal tissue were lower in subjects heterozygous for the Trp64Arg mutation in the β3-adrenergic receptor gene, this difference was not statistically significant. However, a strong association of intraperitoneal UCP-1 mRNA abundance with the UCP-1 gene polymorphism at –3826 relative to the transcription start site was observed that explained 19.3% of the interindividual variability. The minor allele imparted a dose-dependent reduction on UCP gene expression. The importance of sequence variations at the UCP-1 gene locus as a common source of UCP-1 mRNA abundance variability was supported by allele-specific expression studies utilizing a newly identified polymorphism in exon 2 of the UCP-1 gene that predicts a substitution of alanine by threonine. In four subjects heterozygous for the –3826 polymorphism, the mRNA species transcribed from the wild-type allele accounted for 63 ± 6% percent of total intraperitoneal mRNA abundance. In one subject homozygous for the minor promoter allele, wild-type mRNA was also more abundant than variant mRNA. Thus, the UCP-1 polymorphism at –3826 is probably only a marker for a frequent mutation causing reduced mRNA expression.—Esterbauer, H., H. Oberkofler, Y-M. Liu, D. Breban, E. Hell, F. Krempler, and W. Patsch. Uncoupling protein-1 mRNA expression in obese human subjects: the role of sequence variations at the uncoupling protein-1 gene locus. J. Lipid Res. 1998. 39: 834–844.