Allele-specific Expression Analysis Research Articles

Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes.

A major goal in evolutionary biology and biomedicine is to understand the complex interactions between genetic variants, the epigenome, and gene expression. However, the causal relationships between these factors remain poorly understood. mSTARR-seq, a methylation-sensitive massively parallel reporter assay, is capable of identifying methylation-dependent regulatory activity at many thousands of genomic regions simultaneously, and allows for the testing of causal relationships between DNA methylation and gene expression on a region-by-region basis. Here, we developed a multiplexed mSTARR-seq protocol to assay naturally occurring human genetic variation from 25 individuals sampled from 10 localities in Europe and Africa. We identified 6,957 regulatory elements in either the unmethylated or methylated state, and this set was enriched for enhancer and promoter annotations, as expected. The expression of 58% of these regulatory elements was modulated by methylation, which was generally associated with decreased RNA expression. Within our set of regulatory elements, we used allele-specific expression analyses to identify 8,020 sites with genetic effects on gene regulation; further, we found that 42.3% of these genetic effects varied between methylated and unmethylated states. Sites exhibiting methylation-dependent genetic effects were enriched for GWAS and EWAS annotations, implicating them in human disease. Compared to datasets that assay DNA from a single European individual, our multiplexed assay uncovers dramatically more genetic effects and methylation-dependent genetic effects, highlighting the importance of including diverse individuals in assays which aim to understand gene regulatory processes.

Bayesian Estimation of Allele-Specific Expression in the Presence of Phasing Uncertainty.

Allele-specific expression (ASE) analyses aim to detect imbalanced expression of maternal versus paternal copies of an autosomal gene. Such allelic imbalance can result from a variety of cis-acting causes, including disruptive mutations within one copy of a gene that impact the stability of transcripts, as well as regulatory variants outside the gene that impact transcription initiation. Current methods for ASE estimation suffer from a number of shortcomings, such as relying on only one variant within a gene, assuming perfect phasing information across multiple variants within a gene, or failing to account for alignment biases and possible genotyping errors. We developed BEASTIE, a Bayesian hierarchical model designed for precise ASE quantification at the gene level, based on given genotypes and RNA-Seq data. BEASTIE addresses the complexities of allelic mapping bias, genotyping error, and phasing errors by incorporating empirical phasing error rates derived from Genome-in-a-Bottle individual NA12878. BEASTIE surpasses existing methods in accuracy, especially in scenarios with high phasing errors. This improvement is critical for identifying rare genetic variants often obscured by such errors. Through rigorous validation on simulated data and application to real data from the 1000 Genomes Project, we establish the robustness of BEASTIE. These findings underscore the value of BEASTIE in revealing patterns of ASE across gene sets and pathways. The software is freely available from https://github.com/x811zou/BEASTIE . BEASTIE is available as Python source code and as a Docker image. Additional information is available online.

Bi-parental 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 two hybrids, an intraspecific hybrid between two maize (Zea may ssp. mays) inbred lines and an interspecific hybrid between maize and its wild relative teosinte (Zea may ssp. parviglumis), utilizing a combination of PacBio High Fidelity (HiFi) 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 bi-parental genome graph, the haplotypic assemblies can facilitate downstream short-reads-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.

The structure of the TH/INS locus and the parental allele expressed are not conserved between mammals

Parent-of-origin-specific expression of imprinted genes is critical for successful mammalian growth and development. Insulin, coded by the INS gene, is an important growth factor expressed from the paternal allele in the yolk sac placenta of therian mammals. The tyrosine hydroxylase gene TH encodes an enzyme involved in dopamine synthesis. TH and INS are closely associated in most vertebrates, but the mouse orthologues, Th and Ins2, are separated by repeated DNA. In mice, Th is expressed from the maternal allele, but the parental origin of expression is not known for any other mammal so it is unclear whether the maternal expression observed in the mouse represents an evolutionary divergence or an ancestral condition. We compared the length of the DNA segment between TH and INS across species and show that separation of these genes occurred in the rodent lineage with an accumulation of repeated DNA. We found that the region containing TH and INS in the tammar wallaby produces at least five distinct RNA transcripts: TH, TH-INS1, TH-INS2, lncINS and INS. Using allele-specific expression analysis, we show that the TH/INS locus is expressed from the paternal allele in pre- and postnatal tammar wallaby tissues. Determining the imprinting pattern of TH/INS in other mammals might clarify if paternal expression is the ancestral condition which has been flipped to maternal expression in rodents by the accumulation of repeat sequences.

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.

Abstract C050: Allele-specific expression analysis identifies novel dysregulated genes in high-risk prostate tumors from African-American men

Abstract African-American (AA) patients are twice as likely to die from prostate cancer (PCa) compared to European American (EA) patients. Racism, socioeconomic status, access to healthcare, and tumor biology have been attributed to this mortality gap. Tumor biology has not been fully characterized in AA patients, particularly the landscape of gene-regulatory alterations and how they may contribute to tumor aggressiveness. We use allele specific expression (ASE) analysis, which detects genes showing imbalanced expression of alleles often due to altered gene regulation on one allele, to comprehensively identify dysregulated genes in high-risk PCa tumors from AA patients. We develop a framework for identifying candidate cancer genes by integrating ASE data, clinical data, and epigenetics. Our approach identifies genes that recurrently undergo ASE across AA patient tumors. Further, we describe genes that more frequently undergo ASE in AA patient tumors compared to EA patient tumors, some of which are implicated in hormone therapy resistance, aggressive disease, and worse disease-free survival. Our work demonstrates the potential for ASE to uncover dysregulated genes that may not have been identified through RNAseq alone. Ultimately, this work increases our understanding of the effects of regulatory alterations in tumors and how they may contribute to disparities in PCa. Citation Format: Margaret Tsui, Kevin Hu, Hanbing Song, Franklin W. Huang. Allele-specific expression analysis identifies novel dysregulated genes in high-risk prostate tumors from African-American men [abstract]. In: Proceedings of the 16th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2023 Sep 29-Oct 2;Orlando, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2023;32(12 Suppl):Abstract nr C050.

Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects

Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.

Transcriptome-wide identification and characterization of genes exhibit allele-specific imprinting in maize embryo and endosperm

BackgroundGenomic imprinting refers to a subset of genes that are expressed from only one parental allele during seed development in plants. Studies on genomic imprinting have revealed that intraspecific variations in genomic imprinting expression exist in naturally genetic varieties. However, there have been few studies on the functional analysis of allele-specific imprinted genes.ResultsHere, we generated three reciprocal crosses among the B73, Mo17 and CAU5 inbred lines. Based on the transcriptome-wide analysis of allele-specific expression using RNA sequencing technology, 305 allele-specific imprinting genes (ASIGs) were identified in embryos, and 655 ASIGs were identified in endosperms from three maize F1 hybrids. Of these ASIGs, most did not show consistent maternal or paternal bias between the same tissue from different hybrids or different tissues from one hybrid cross. By gene ontology (GO) analysis, five and eight categories of GO exhibited significantly higher functional enrichments for ASIGs identified in embryo and endosperm, respectively. These functional categories indicated that ASIGs are involved in intercellular nutrient transport, signaling pathways, and transcriptional regulation of kernel development. Finally, the mutation and overexpression of one ASIG (Zm305) affected the length and width of the kernel.ConclusionIn this study, our data will be helpful in gaining further knowledge of genes exhibiting allele-specific imprinting patterns in seeds. The gain- and loss-of-function phenotypes of ASIGs associated with agronomically important seed traits provide compelling evidence for ASIGs as crucial targets to optimize seed traits in crop plants.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

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

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. 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%. Analysis pipeline for this approach is available at GitHub as R package controlFreq (github.com/gimelbrantlab/controlFreq).

Dynamic Changes in the Global Transcriptome of Postnatal Skeletal Muscle in Different Sheep.

Sheep growth performance, mainly skeletal muscle growth, provides direct economic benefits to the animal husbandry industry. However, the underlying genetic mechanisms of different breeds remain unclear. We found that the cross-sectional area (CSA) of skeletal muscle in Dorper (D) and binary cross-breeding (HD) was higher than that in Hu sheep (H) from 3 months to 12 months after birth. The transcriptomic analysis of 42 quadriceps femoris samples showed that a total of 5053 differential expression genes (DEGs) were identified. The differences in the global gene expression patterns, the dynamic transcriptome of skeletal muscle development, and the transcriptome of the transformation of fast and slow muscles were explored using weighted correlation network analysis (WGCNA) and allele-specific expression analysis. Moreover, the gene expression patterns of HD were more similar to D rather than H from 3 months to 12 months, which might be the reason for the difference in muscle growth in the three breeds. Additionally, several genes (GNB2L1, RPL15, DVL1, FBXO31, etc.) were identified as candidates related to skeletal muscle growth. These results should serve as an important resource revealing the molecular basis of muscle growth and development in sheep.

Exploration of genotype-by-environment interactions affecting gene expression responses in porcine immune cells.

As one of the keys to healthy performance, robustness of farm animals is gaining importance, and with this comes increasing interest in genetic dissection of genotype-by-environment interactions (G×E). Changes in gene expression are among the most sensitive responses conveying adaptation to environmental stimuli. Environmentally responsive regulatory variation thus likely plays a central role in G×E. In the present study, we set out to detect action of environmentally responsive cis-regulatory variation by the analysis of condition-dependent allele specific expression (cd-ASE) in porcine immune cells. For this, we harnessed mRNA-sequencing data of peripheral blood mononuclear cells (PBMCs) stimulated in vitro with lipopolysaccharide, dexamethasone, or their combination. These treatments mimic common challenges such as bacterial infection or stress, and induce vast transcriptome changes. About two thirds of the examined loci showed significant ASE in at least one treatment, and out of those about ten percent exhibited cd-ASE. Most of the ASE variants were not yet reported in the PigGTEx Atlas. Genes showing cd-ASE were enriched in cytokine signaling in immune system and include several key candidates for animal health. In contrast, genes showing no ASE featured cell-cycle related functions. We confirmed LPS-dependent ASE for one of the top candidates, SOD2, which ranks among the major response genes in LPS-stimulated monocytes. The results of the present study demonstrate the potential of in vitro cell models coupled with cd-ASE analysis for the investigation of G×E in farm animals. The identified loci may benefit efforts to unravel the genetic basis of robustness and improvement of health and welfare in pigs.

Characterization of Genes That Exhibit Genotype-Dependent Allele-Specific Expression and Its Implications for the Development of Maize Kernel

Heterosis or hybrid vigor refers to the superior phenotypic traits of hybrids relative to their parental inbred lines. An imbalance between the expression levels of two parental alleles in the F1 hybrid has been suggested as a mechanism of heterosis. Here, based on genome-wide allele-specific expression analysis using RNA sequencing technology, 1689 genes exhibiting genotype-dependent allele-specific expression (genotype-dependent ASEGs) were identified in the embryos, and 1390 genotype-dependent ASEGs in the endosperm, of three maize F1 hybrids. Of these ASEGs, most were consistent in different tissues from one hybrid cross, but nearly 50% showed allele-specific expression from some genotypes but not others. These genotype-dependent ASEGs were mostly enriched in metabolic pathways of substances and energy, including the tricarboxylic acid cycle, aerobic respiration, and energy derivation by oxidation of organic compounds and ADP binding. Mutation and overexpression of one ASEG affected kernel size, which indicates that these genotype-dependent ASEGs may make important contributions to kernel development. Finally, the allele-specific methylation pattern on genotype-dependent ASEGs indicated that DNA methylation plays a potential role in the regulation of allelic expression for some ASEGs. In this study, a detailed analysis of genotype-dependent ASEGs in the embryo and endosperm of three different maize F1 hybrids will provide an index of genes for future research on the genetic and molecular mechanism of heterosis.

Red-TE Homozygous Alleles of MdMYB10 Confer Full-Red Apple Fruit Skin in a High-Temperature Region

Apple is a major fruit crop grown worldwide and provides humans with an essential diet and health benefits. One of the health benefits is related to the accumulation of fruit anthocyanin, which also provides fruit with an attractive red colour. It is known that an MdMYB10 allele containing a transposable element (TE) insertion in its promoter (termed Red-TE allele) underlies anthocyanin accumulation and red colouration in the fruit skin of cultivated apples. However, the distribution of this Red-TE allele in wider Malus germplasm accessions is not clear. In this study, we showed that MdMYB10 RNA in fruit skin was specifically expressed from the Red-TE allele by using allele-specific expression analysis of transcriptome data. Apple cultivars and hybrids with homozygous Red-TE alleles showed stronger red colour than those with heterozygous alleles after analysing 65 cultivars and 337 hybrids. Furthermore, both hetero- and homozygous plants growing in the same high-temperature conditions had different colourations. However, the Red-TE allele was not detected in 16 wild apple accessions showing red skin, indicating that the red skin colour of these wild apples was not conferred by the Red-TE allele. These findings provide guidance for selecting cultivars able to develop consistent red colouration under high growth temperature conditions and open the opportunity for identifying novel genetic variants underpinning fruit red colouration in wild apple species.

Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines.

Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs. We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments. We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes. Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

Allele-specific expression analysis for complex genetic phenotypes applied to a unique dilated cardiomyopathy cohort

Allele-specific expression (ASE) analysis detects the relative abundance of alleles at heterozygous loci as a proxy for cis-regulatory variation, which affects the personal transcriptome and proteome. This study describes the development and application of an ASE analysis pipeline on a unique cohort of 87 well phenotyped and RNA sequenced patients from the Maastricht Cardiomyopathy Registry with dilated cardiomyopathy (DCM), a complex genetic disorder with a remaining gap in explained heritability. Regulatory processes for which ASE is a proxy might explain this gap. We found an overrepresentation of known DCM-associated genes among the significant results across the cohort. In addition, we were able to find genes of interest that have not been associated with DCM through conventional methods such as genome-wide association or differential gene expression studies. The pipeline offers RNA sequencing data processing, individual and population level ASE analyses as well as group comparisons and several intuitive visualizations such as Manhattan plots and protein–protein interaction networks. With this pipeline, we found evidence supporting the case that cis-regulatory variation contributes to the phenotypic heterogeneity of DCM. Additionally, our results highlight that ASE analysis offers an additional layer to conventional genomic and transcriptomic analyses for candidate gene identification and biological insight.

Testing for Allele-specific Expression from Human Brain Samples.

Many single nucleotide polymorphisms (SNPs) identified by genome-wide association studies exert their effects on disease risk as expression quantitative trait loci (eQTL) via allele-specific expression (ASE). While databases for probing eQTLs in tissues from normal individuals exist, one may wish to ascertain eQTLs or ASE in specific tissues or disease-states not characterized in these databases. Here, we present a protocol to assess ASE of two possible target genes (GPNMB and KLHL7) of a known genome-wide association study (GWAS) Parkinson's disease (PD) risk locus in postmortem human brain tissue from PD and neurologically normal individuals. This was done using a sequence of RNA isolation, cDNA library generation, enrichment for transcripts of interest using customizable cDNA capture probes, paired-end RNA sequencing, and subsequent analysis. This method provides increased sensitivity relative to traditional bulk RNAseq-based and a blueprint that can be extended to the study of other genes, tissues, and disease states. Key features • Analysis of GPNMB allele-specific expression (ASE) in brain lysates from cognitively normal controls (NC) and Parkinson's disease (PD) individuals. • Builds on the ASE protocol of Mayba et al. (2014) and extends application from cells to human tissue. • Increased sensitivity by enrichment for desired transcript via RNA CaptureSeq (Mercer et al., 2014). • Optimized for human brain lysates from cingulate gyrus, caudate nucleus, and cerebellum.

Insertion of a TRIM-like sequence in MdFLS2-1 promoter is associated with its allele-specific expression in response to Alternaria alternata in apple.

Alternaria blotch disease, caused by Alternaria alternata apple pathotype (AAAP), is one of the major fungal diseases in apple. Early field observations revealed, the anther-derived homozygote Hanfu line (HFTH1) was highly susceptible to AAAP, whereas Hanfu (HF) exhibited resistance to AAAP. To understand the molecular mechanisms underlying the difference in sensitivity of HF and HFTH1 to AAAP, we performed allele-specific expression (ASE) analysis and comparative transcriptomic analysis before and after AAAP inoculation. We reported an important immune gene, namely, MdFLS2, which displayed strong ASE in HF with much lower expression levels of HFTH1-derived alleles. Transient overexpression of the dominant allele of MdFLS2-1 from HF in GL-3 apple leaves could enhance resistance to AAAP and induce expression of genes related to salicylic acid pathway. In addition, MdFLS2-1 was identified with an insertion of an 85-bp terminal-repeat retrotransposon in miniature (TRIM) element-like sequence in the upstream region of the nonreference allele. In contrast, only one terminal direct repeat (TDR) from TRIM-like sequence was present in the upstream region of the HFTH1-derived allele MdFLS2-2. Furthermore, the results of luciferase and β-glucuronidase reporter assays demonstrated that the intact TRIM-like sequence has enhancer activity. This suggested that insertion of the TRIM-like sequence regulates the expression level of the allele of MdFLS2, in turn, affecting the sensitivity of HF and HFTH1 to AAAP.

Pklr Variants Associated with Acute Pain in Sickle Cell Disease Influence ATP Concentrations in Red Blood Cells

Background: We previously demonstrated association of PKLR intron 2 variants with hospitalization rate for acute pain in people with sickle cell disease (SCD) and utilized allele-specific expression analysis to show that these variants influence PKLR activity (Wang et al, 2022). Pyruvate kinase (gene PKLR, protein PKR) is a rate-limiting enzyme in glycolysis. Reduced PKR activity leads to an increase in intracellular 2,3-diphosphoglycerate (2,3-DPG), and a concomitant decrease of adenosine triphosphate (ATP), factors that promote sickle hemoglobin (HbS) polymerization and red blood cell (RBC) sickling. On this basis, we explored how the "high risk" PKLR intron 2 variants affect levels of PKR protein expression, 2,3-DPG, ATP, oxygen dissociation (p50) and sickling kinetics (t50). Methods: We studied a total of 416 subjects which included 125 HbSS, 89 HbAS, and 202 ethnic-matched healthy controls (HbAA) enrolled under protocol NCT03685721 approved by the NHLBI Institutional Review Board. DNA was genome scanned using Ilumina's Infinium "MEGA" chip (1.7m markers). The results were quality controlled followed by genotype imputation based on the 1000 Genomes Project phase 3 data using the Michigan Imputation Server. Whole blood levels of ATP and 2,3-DPG were measured using LC-MS/MS with LLOQ at 50.0 μg/mL and converted to intracellular concentrations by dividing by the hematocrit (as a fraction). PKR was detected by antibody-based capture and detection using an electrochemiluminescence immunoassay (Meso Scale Discovery, Rockville, MD, USA) and the signal normalized to a control sample. A Hemox Analyzer [TCS, Scientific Corp, New Hope, PA] was used to obtain p50 values (defined as the partial pressure of oxygen at which Hb is 50% saturated with oxygen). t50 is single parameter measure of sickling kinetics (Dunkelberger et al, 2018). Single variant analyses (under additive genetic model adjusted for age and sex) were conducted in each of the Hb sub-populations for associations between the 7 PKLR variants (intron 4 - rs2071053, and intron 2 -rs8177970, rs116244351, rs114455416, rs12741350, rs3020781, and rs8177964) and the 5 parameters (PKR, 2,3-DPG, ATP, p50 and t50). False-discovery rate (FDR) adjusted p-values were calculated to account for the multiple testing issue. Results: Table 1A shows levels of the 5 parameters (ATP/HCT, 2,3-DPG/HCT, PKR, p50 and t50) in the 3 Hb genotypic groups (HbAA, HbAS and HbSS). Because patients with SCD have variable anemia, ATP and 2,3-DPG were normalized against HCT. The observed trending higher ATP / HCT and significantly elevated 2,3-DPG /HCT in HbSS compared to HbAS and to HbAA are consistent with previous reports (Poillon et al, 1990). HbSS subjects also had higher PKR protein expression, as well as right-shifted p50 values, compared to HbAS and HbAA (Table 1A), which suggest HbSS RBCs have lower specific activity of PKR. A one-way ANOVA showed that the differences in these mean values were highly significant (Table 1A). All 7 PKLR variants were associated with ATP levels in the HbAS group, of which 4 (rs8177970, rs116244351, rs114455416 and rs8177964) in intron 2 remained statistically significant after multiple testing adjustment (Table 1B). We also identified association between 3 variants (intron 2 - rs12741350, rs3020781, intron 4 - rs071053) and 2,3-DPG levels in the HbSS, as well as t50 in the HbAS group. Conclusion: ATP is essential for maintaining RBC membrane integrity and shape; reduced ATP dehydrates RBCs and promotes sickling. In a previous report, a PKLR intron 2 'high-risk" haplotype was implicated in acute sickle pain by affecting PKLR expression (Wang et al, 2022). In this follow-up study, we show that the same intron 2 "high-risk" haplotype negatively affects ATP levels when accounting for PKR expression, thus providing a biological basis for PKLR as a genetic modifier of SCD pain presentation. 3 of the 7 variants are also associated with elevated 2,3-DPG levels in HbSS as well as t50 (a single parameter that depends primarily on sickling kinetics) in HbAS. The variant associations are observed more clearly in HbAS but not in HbSS, possibly because of the much more homogeneous cell population in HbAS. We conclude that PKLR intron variants may directly contribute to the severity and frequency of acute pain episodes in SCD, thus additional studies with larger sample sizes are warranted. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal