Polymorphisms affecting gene transcription and mRNA processing in pharmacogenetic candidate genes: detection through allelic expression imbalance in human target tissues

The role of inflammatory gene polymorphisms in severe COVID-19: a review

The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, has profoundly impacted global healthcare systems and spurred extensive research efforts over the past three years. One critical aspect of the disease is the intricate interplay between the virus and the host immune response, particularly the role of inflammatory gene expression in severe COVID-19. While numerous previous studies have explored the role of genetic polymorphisms in COVID-19, research specifically focusing on inflammatory genes and their associations with disease severity remains limited. This review explores the relationship between severe COVID-19 outcomes and genetic polymorphisms within key inflammatory genes. By investigating the impact of genetic variations on immune responses, which include cytokine production and downstream signalling pathways, we aim to provide a comprehensive overview of how genetic polymorphisms contribute to the variability in disease presentation. Through an in-depth analysis of existing literature, we shed light on potential therapeutic targets and personalized approaches that may enhance our understanding of disease pathogenesis and treatment strategies.

BER genes expression in oral and pre-oral cancer: Combinatorial approach to propose potential biomarker

OBJECTIVE

DNA repair genes and their variants have been found to alter the risk of oral cancer.

METHOD

The level of expression of XRCC3, NBS1, and OGG1 genes among 20 cases of oral cancer, 6 pre-oral cancer, and 50 healthy control subjects was measured with RT-PCR. All the subjects were also genotyped for XRCC3 rs861539 C>T, NBS1 rs1805794 C>G, and OGG1 rs1052133 C>G polymorphisms by the PCR-RFLP method; their genotypes were correlated with their level of expression. Further, a localized fold structure analysis of the mRNA sequence surrounding the studied SNPs was performed with RNAfold.

RESULTS

Results showed increased expression of XRCC3, NBS1, and OGG1 transcripts among oral cancer (4.49 fold, 3.45 fold, and 3.27 fold) as well as pre-oral cancer (3.04 fold, 5.32 fold, and 1.74 fold) as compared to control subjects. The transcript level of OGG1 was found to be significantly increased (6.68 fold, p-value 0.009) with the GG genotype compared to the CC genotype. The C>T polymorphism of XRCC3 and the C>G polymorphism of OGG1 result in an apparent change in its mRNA secondary structure. Folding energy with the C allele for XRCC3 C>T polymorphism was lower than that of the T allele (MFE C vs T: -50.20 kcal/mol vs -48.70 kcal/mol). In the case of OGG1 C>G polymorphism MFE for the C allele was higher (-23.30 kcal/mole) than with the G allele (-24.80 kcal/mol).

CONCLUSION

Our results showed elevated levels of XRCC3, NBS1, and OGG1 both in oral cancer and pre-oral cancer conditions, which indicates their role as prospective biomarkers of oral cancer and pre-cancerous lesions. SNPs in these genes alter their level of expression, possibly by altering the secondary structure of their transcript. However, due to the small sample size our study can only provide a suggestive conclusion and warned future study with large sample size to verify our findings.

Gene expression profiling of vitamin D metabolism enzymes in leukemia and lymphoma patients: molecular aspect interplay of VDR, CYP2R1, and CYP24A1

BACKGROUND

Vitamin D deficiency is prevalent among the Indonesian population, particularly in individuals diagnosed with leukemia-lymphoma. The regulation of vitamin D metabolism is influenced by the expression of several enzymes, such as CYP2R1, CYP24A1, and the vitamin D receptor (VDR). This study aimed to scrutinize the gene expression profiles in both mRNA and protein levels of VDR, CYP2R1, and CYP24A1 in leukemia and lymphoma patients.

METHOD

The research was a cross-sectional study conducted at Cipto Mangunkusumo Hospital (RSCM) in Jakarta, Indonesia. The study included a total of 45 patients aged over 18 years old who have received a diagnosis of lymphoma or leukemia. Vitamin D status was measured by examining serum 25 (OH) D levels. The analysis of VDR, CYP2R1, and CYP24A1 mRNA expression utilized the qRT-PCR method, while protein levels were measured through the ELISA method.

CONCLUSION

The study revealed a noteworthy difference in VDR protein levels between men and women. The highest mean CYP24A1 protein levels were observed in the age group > 60 years. This study found a significant, moderately positive correlation between VDR protein levels and CYP24A1 protein levels in the male and vitamin D sufficiency groups. In addition, a significant positive correlation was found between VDR mRNA levels and CYP2R1 mRNA levels, VDR mRNA levels and CYP2R1 mRNA levels, and CYP2R1 mRNA levels and CYP24A1 mRNA levels. However, the expression of these genes does not correlate with the protein levels of its mRNA translation products in blood circulation.

Modulation of gene expression and influence of gene polymorphisms related to genotoxicity and redox status on occupational exposure to inhaled anesthetics

The extensive use of inhalational anesthetics contributes to both indoor and outdoor (environmental) pollution. The influence of genetic susceptibility on DNA damage and oxidative stress and the possible modulation of gene expression have not yet been investigated upon occupational exposure to waste anesthetic gases (WAGs). This study assessed 8-oxoguanine DNA glycosylase 1 (OGG1) and superoxide dismutase 2 (SOD2) gene expression, which are related to oxidized DNA repair and antioxidant capacity, respectively, and the influence of their polymorphisms (OGG1 rs1052133 and SOD2 rs4880) in 100 professionals highly exposed to WAGs and 93 unexposed volunteers (control group). Additionally, X-ray repair cross complementing 1 (XRCC1 rs25487 and rs1799782) and ataxia telangiectasia mutated (ATM rs600931) gene polymorphisms as well as genetic instability (micronucleus-MN and nuclear bud-NBUD) and oxidative stress (malondialdehyde-MDA and ferric reducing antioxidant power-FRAP) biomarkers were assessed in the groups (control and exposed) and in the subgroups of the exposed group according to job occupation (anesthesiologists versus surgeons/technicians). Except for the ATM TT controls (associated with increased FRAP), there were no influences of OGG1, XRCC1, ATM, and SOD2 polymorphisms on MN, NBUD, MDA, and FRAP values in exposed or control subjects. No significant difference in the expression of either gene evaluated (OGG1 and SOD2) was found between the exposed and control groups. Increased OGG1 expression was observed among OGG1 -/Cys individuals only in the control group. Among the exposed group, anesthesiologists had a greater duration of WAG exposure (both h/week and years) than surgeons/technicians, which was associated with increased MDA and decreased antioxidant capacity (FRAP) and SOD2 expression (redox status). Higher expression of OGG1 was found in -/Cys surgeons/technicians than in anesthesiologists with the same genotype. Increased antioxidant capacity was noted in the surgeons/technicians carrying the ATM T allele and in those carrying XRCC1 -/Gln. Increased MN was influenced by OGG1 -/Cys in surgeons/technicians. Anesthesiologists with ATM CC exhibited increased MN, and those carrying the C allele (CC/CT genotype) exhibited increased NBUD. SOD2 polymorphism did not seem to be relevant for WAG exposure. These findings contribute to advancing the knowledge on genetic susceptibility/gene expression/genetic instability/oxidative stress, including differences in job occupation considering the workload, in response to occupational exposure to WAGs.

Genetic polymorphisms as potential pharmacogenetic biomarkers for platinum-based chemotherapy in non-small cell lung cancer

Platinum-based chemotherapy (PBC) is a widely used treatment for various solid tumors, including non-small cell lung cancer (NSCLC). However, its efficacy is often compromised by the emergence of drug resistance in patients. There is growing evidence that genetic variations may influence the susceptibility of NSCLC patients to develop resistance to PBC. Here, we provide a comprehensive overview of the mechanisms underlying platinum drug resistance and highlight the important role that genetic polymorphisms play in this process. This paper discussed the genetic variants that regulate DNA repair, cellular movement, drug transport, metabolic processing, and immune response, with a focus on their effects on response to PBC. The potential applications of these genetic polymorphisms as predictive indicators in clinical practice are explored, as are the challenges associated with their implementation.

The role of SOD2 and NOS2 genes in the molecular aspect of bladder cancer pathophysiology

Bladder cancer (BC) is a severe health problem of the genitourinary system and is characterised by a high risk of recurrence. According to the recent GLOBOCAN report, bladder cancer accounts for 3% of diagnosed cancers in the world, taking 10th place on the list of the most common cancers. Despite numerous studies, the full mechanism of BC development remains unknown. Nevertheless, precious results suggest a crucial role of oxidative stress in the development of BC. Therefore, this study explores whether the c. 47 C > T (rs4880)-SOD2, (c. 1823 C > T (rs2297518) and g.-1026 C > A (rs2779249)-NOS2(iNOS) polymorphisms are associated with BC occurrence and whether the bladder carcinogenesis induces changes in SOD2 and NOS2 expression and methylation status in peripheral blood mononuclear cells (PBMCs). In this aim, the TaqMan SNP genotyping assay, TaqMan Gene Expression Assay, and methylation-sensitive high-resolution melting techniques were used to genotype profiling and evaluate the expression of the genes and the methylation status of their promoters, respectively. Our findings confirm that heterozygote of the g.-1026 C > A SNP was associated with a decreased risk of BC. Moreover, we detected that BC development influenced the expression level and methylation status of the promoter region of investigated genes in PBMCs. Concluding, our results confirmed that oxidative stress, especially NOS2 polymorphisms and changes in the expression and methylation of the promoters of SOD2 and NOS2 are involved in the cancer transformation initiation of the cell urinary bladder.

Pharmacogenomics: Driving Personalized Medicine

Personalized medicine tailors therapies, disease prevention, and health maintenance to the individual, with pharmacogenomics serving as a key tool to improve outcomes and prevent adverse effects. Advances in genomics have transformed pharmacogenetics, traditionally focused on single gene-drug pairs, into pharmacogenomics, encompassing all "-omics" fields (e.g., proteomics, transcriptomics, metabolomics, and metagenomics). This review summarizes basic genomics principles relevant to translation into therapies, assessing pharmacogenomics' central role in converging diverse elements of personalized medicine. We discuss genetic variations in pharmacogenes (drug-metabolizing enzymes, drug transporters, and receptors), their clinical relevance as biomarkers, and the legacy of decades of research in pharmacogenetics. All types of therapies, including proteins, nucleic acids, viruses, cells, genes, and irradiation, can benefit from genomics, expanding the role of pharmacogenomics across medicine. Food and Drug Administration approvals of personalized therapeutics involving biomarkers increase rapidly, demonstrating the growing impact of pharmacogenomics. A beacon for all therapeutic approaches, molecularly targeted cancer therapies highlight trends in drug discovery and clinical applications. To account for human complexity, multicomponent biomarker panels encompassing genetic, personal, and environmental factors can guide diagnosis and therapies, increasingly involving artificial intelligence to cope with extreme data complexities. However, clinical application encounters substantial hurdles, such as unknown validity across ethnic groups, underlying bias in health care, and real-world validation. This review address the underlying science and technologies germane to pharmacogenomics and personalized medicine, integrated with economic, ethical, and regulatory issues, providing insights into the current status and future direction of health care. SIGNIFICANCE STATEMENT: Personalized medicine aims to optimize health care for the individual patients with use of predictive biomarkers to improve outcomes and prevent adverse effects. Pharmacogenomics drives biomarker discovery and guides the development of targeted therapeutics. This review addresses basic principles and current trends in pharmacogenomics, with large-scale data repositories accelerating medical advances. The impact of pharmacogenomics is discussed, along with hurdles impeding broad clinical implementation, in the context of clinical care, ethics, economics, and regulatory affairs.

Functional single nucleotide polymorphisms in dopaminergic receptors D2 predict clinical response to Cariprazine

Cariprazine (CAR) is an antipsychotic drug for the treatment of schizophrenia (SCZ) and bipolar disorder (BD), and it acts as a partial agonist on the dopamine receptors (DR), D2, and D3. Although many single nucleotide polymorphisms (SNPs) in genes coding for these receptors are known to influence response to antipsychotics, to date, no study on CAR pharmacogenetics exists. In this pilot study, we investigated the relationship between SNPs in DRD2 (rs1800497 and rs6277) and DRD3 (rs6280), and response to CAR treatment, evaluated by the psychometric Brief Psychiatric Rating Scale (BPRS), in a cohort of Caucasian patients. We found a significant association between DRD2 rs1800497 and rs6277 and response to CAR treatment. When genotypes were combined into an arbitrary score, the receiver operating characteristic curve analysis showed that using a cut-off value of -2.5 the response to CAR treatment could be predicted with a positive likelihood ratio of 8.0. Our study report, for the first time, a correlation between SNPs in DRD2 and response to CAR treatment. After confirmation in a larger cohort of patients, our results could open the way for the identification of new tools for the provision of response to CAR treatment.

Integrative analysis of omics summary data reveals putative mechanisms linked to different cell populations in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex disease involving many interactions at the molecular level, the details of which remain unclear. Here, we demonstrated an analytical paradigm of prioritizing genes and regulatory elements based on GWAS loci at the single-cell levels. Our initial step was to apply TWMR to identify causal genes and causal methylation sites in SLE. Based on the eQTL, LD and mQTL, we calculated the correlation between these genes and methylation sites. Next, we separately used gene expression and DNAm as exposure variables and outcome variables to analyze the regulatory mechanisms. We identified two mediating modes for SLE: 1) transcription mediation model and 2) epigenetic mediation model. Further, using single-cell RNA sequencing data, we revealed the cell subclusters associated with these mechanisms. Our identification of the mechanisms of SLE in different cell populations is of great significance for understanding the heterogeneity of disease in different cell populations.

Role of metabolizing MTHFR gene polymorphism (rs1801133) and its mRNA expression among Type 2 Diabetes

OBJECTIVE

Type 2 Diabetes is a glucose metabolic disorder occurred by insulin insensitivity in which folate metabolism plays an important role. it is believed that polymorphism of Methylenetetrahydrofolate reductase (MTHFR) C677T linked with type 2 diabetes mellitus. However, results are conflicted. therefore, in this study we re-examine the relationship between MTHFR C677T in type 2 diabetes mellitus patients.

METHODS

Present research work included 100 newly diagnosed type 2 diabetic mellitus (T2DM) cases and 100 healthy individuals. After the blood sample collection all the biochemical parameters were evaluated among the T2DM cases and healthy individuals. DNA and RNA extraction from whole blood was done to study the MTHFR gene polymorphism by allele specific polymerase chain reaction method and its expression analysis was done by quantitative real time polymerase chain reaction method.

RESULTS

The significant difference was observed in genotype distribution among case and control group (p=0.0002). Compared with wildtype CC genotype, CT heterozygous (OR=2.95, 95% Cl=1.62-5.38) and TT homozygous (OR=3.20, CI=1.79-5.73) suggest to have effect of MTHFR polymorphism on type 2 mellitus risk. Moreover, relative MTHFR mRNA expression was found for wild type CC genotype 3.02-fold, CT heterozygous genotype 2.57 fold and mutant TT homozygous genotype 0.50-fold which is down regulated (p<0.0001).

CONCLUSION

Our results indicates that the polymorphism in MTHFR C677T plays significant role in type II diabetes risk. MTHFR CT heterozygous and mutant TT genotype showed reduced mRNA expression among the T2DM patients. However, large scale case-control studies are needed to strengthen such conclusion in the future.

Interpreting coronary artery disease GWAS results: A functional genomics approach assessing biological significance

Genome-wide association studies (GWAS) have implicated 58 loci in coronary artery disease (CAD). However, the biological basis for these associations, the relevant genes, and causative variants often remain uncertain. Since the vast majority of GWAS loci reside outside coding regions, most exert regulatory functions. Here we explore the complexity of each of these loci, using tissue specific RNA sequencing data from GTEx to identify genes that exhibit altered expression patterns in the context of GWAS-significant loci, expanding the list of candidate genes from the 75 currently annotated by GWAS to 245, with almost half of these transcripts being non-coding. Tissue specific allelic expression imbalance data, also from GTEx, allows us to uncover GWAS variants that mark functional variation in a locus, e.g., rs7528419 residing in the SORT1 locus, in liver specifically, and rs72689147 in the GUYC1A1 locus, across a variety of tissues. We consider the GWAS variant rs1412444 in the LIPA locus in more detail as an example, probing tissue and transcript specific effects of genetic variation in the region. By evaluating linkage disequilibrium (LD) between tissue specific eQTLs, we reveal evidence for multiple functional variants within loci. We identify 3 variants (rs1412444, rs1051338, rs2250781) that when considered together, each improve the ability to account for LIPA gene expression, suggesting multiple interacting factors. These results refine the assignment of 58 GWAS loci to likely causative variants in a handful of cases and for the remainder help to re-prioritize associated genes and RNA isoforms, suggesting that ncRNAs maybe a relevant transcript in almost half of CAD GWAS results. Our findings support a multi-factorial system where a single variant can influence multiple genes and each genes is regulated by multiple variants.

Correlation between adiponectin rs2241766 and rs266729 polymorphisms and risk of papillary thyroid cancer

About 60-80% of thyroid cancer (TC) cases are papillary thyroid cancer (PTC). Studies have shown that serum adiponectin levels are inversely related to the risk of TC and PTC. Aim of the present study was to evaluate the association between adiponectin rs2241766 and rs266729 polymorphisms and risk of PTC. 122 PTC patients and 128 healthy subjects were enrolled in the study. PCR-RFLP and ARMS-PCR methods were used for genotype analysis. The rs266729 polymorphism did not correlate with risk of PTC. As regard rs2241766 polymorphism, the frequency of the GG genotype did not have a significant difference between the two groups, although, PTC cases showed higher frequency of GT genotype compared to controls (OR=2.87, 95% CI=1.56-5.28, P=0.001). We observed a significant association between adiponectin rs2241766 polymorphism and PTC, however, our result showed no significant relationship between adiponectin rs266729 polymorphism and risk of PTC.

Possible interplay between estrogen and the BAFF may modify thyroid activity in Graves' disease

A link between sex hormones and B-cell activating factor (BAFF), a crucial immunoregulator of autoimmune thyroid disease (AITD), may exist. The study aimed to elucidate the role of estrogen (E2) in regulating BAFF in Graves' disease (GD). In clinical samples, serum BAFF levels were higher in women than in men in both the GD and control groups. serum BAFF levels were associated with thyroid-stimulating hormone receptor antibody levels and thyroid function only in women and not in men. BAFF transcripts in peripheral blood mononuclear cells were higher in women with GD than those in the control group. Among GD patients with the AA genotype of rs2893321, women had higher BAFF transcripts and protein levels than men. In the progression of a spontaneous autoimmune thyroiditis (SAT) murine model, NOD.H-2h4, serum free thyroxine and BAFF levels were higher in female than in male mice. Moreover, exogenous E2 treatment increased serum BAFF levels in male SAT mice. Meanwhile, female SAT mice exhibited higher thyroid BAFF transcripts levels than either the E2-treated or untreated male SAT mouse groups. Our results showed that E2 might be implicated in modulating BAFF expression, and support a possible mechanism for the higher incidence of AITD in women.

Association of Vitamin D Receptor Polymorphism (rs2228570, rs1544410, rs7975232, and rs731236) and Macrophage Migration Inhibitory Factor -173 G/C (rs755622) with the Susceptibility of Active Pulmonary Tuberculosis in Makassar, Indonesia

BACKGROUND: The study of Vitamin D Receptor (VDR) and Macrophage Migration Inhibitory Factor (MIF) polymorphisms, associated with active pulmonary tuberculosis (ATB) presents varying results. AIMS: This study aimed to investigate the association between VDR rs2228570, rs1544410, rs7975232, rs731236 and MIF -173 G/C (rs755622) single nucleotide polymorphism (SNP), with susceptibility of developing ATB, and positivity of Interferon Gamma Release Assay (IGRA) results (in household contact). METHODS AND MATERIAL: This study involved 83 ATB and 73 household contacts in Makassar. We checked IGRA based on ELISA in household contacts by using QuantiFERON TB Gold Plus test, and we found that 61.64% (n = 45) of household contacts had positive IGRA. Polymorphism examination was carried out by Sanger sequencing. RESULTS: VDR rs2228570 T/T and T/C-T/T were significantly associated with higher risk of active tuberculosis. VDR rs7975232 G/G genotype was associated with an increased risk of developing active TB compared to T/T-T/G. Haplotype analysis of VDR rs2228570, rs1544410, rs7975232, rs731236 and combination with MIF rs755622 demonstrated that TGGTG was observed to have a higher risk of tuberculosis. CONCLUSIONS: The combination of VDR and MIF variants may contribute to the susceptibility of active tuberculosis disease.

A Noncoding Variant Near PPP1R3B Promotes Liver Glycogen Storage and MetS, but Protects Against Myocardial Infarction

CONTEXT

Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation.

OBJECTIVE

Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease.

DESIGN

Genetics of Obesity-associated Liver Disease Consortium.

SETTING

Population-based.

MAIN OUTCOME

Computed tomography measured liver attenuation.

RESULTS

Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate.

CONCLUSIONS

These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.

Association study of rs7799039, rs1137101 and rs8192678 gene variants with disease susceptibility/severity and corresponding LEP, LEPR and PGC1A gene expression in multiple sclerosis

BACKGROUND

Leptin (LEP), leptin receptor (LEPR) and peroxisome proliferator-activated receptor gamma co-activator 1-alpha (PGC1A) are involved in the pathogenesis of multiple sclerosis (MS) by affecting the inflammatory response and reactive oxygen species production. LEP rs7799039 and LEPR rs1137101 genetic variants modify the serum LEP levels and PGC1A rs8192678 alters the PGC1A activity. The study objective was to explore the associations of these variants with susceptibility to MS, disease course/clinical parameters and also with peripheral blood mononuclear cell expression of the target genes and plasma LEP concentrations, in the study subjects.

METHODS

The study groups included 528 patients with MS and 429 controls. TaqMan® assays were used for genotyping and gene expression quantification. The Chi-square, parametric and nonparametric tests and simple/multiple logistic regression were performed for the statistical analysis of data.

RESULTS

A multiple logistic regression model including all three investigated variants, applied to patients (RRMS + SPMS) and controls, showed that PGC1A rs8192678 minor allele had an increased risk for the occurrence of disease, with OR (95%CI) = 1,32 (1,01-1,73), P = 0,04. Between-effect of gender and LEPR variant on the multiple sclerosis severity score (MSSS) was identified (P = 0,005). In male patients (relapsing-remitting and secondary progressive), LEPR minor allele carriers had increased MSSS (GG + AG vs AA, median (minimum-maximum) = 5,38 (0,64-9,88) vs 4,27 (0,78-9,63); P = 0,01, P_adj = 0,03). In relapsing-remitting patients LEP rs7799039 affected the LEP gene expression (P = 0,006; P_adj = 0,04).

CONCLUSION

The current findings implicate an impact of investigated genetic variants on the pathogenesis of MS.

Epigenetics of Primary Biliary Cholangitis

Primary biliary cholangitis (PBC) is a chronic cholestatic liver disease with non-suppurative destruction of the intrahepatic bile ducts. The interplay of genetics and environmental triggers contributes to the onset of the disease and subsequently results in cholestasis and progressive fibrosis. Recently, genome-wide association studies (GWAS) have identified multiple genes influencing the susceptibility to PBC in HLA and non-HLA loci. However, it is estimated that the known risk variants merely account for no more than 20% of the heritability of PBC and causes of the remaining heritability remain uncertain. Increasing evidence suggests that the presence of epigenetic abnormalities may explain the "missing heritability" that cannot be captured by GWAS. Among these epigenetic mechanisms, DNA methylation, histone modification, and noncoding RNAs (i.e. miRNA and lncRNA) are involved in the pathogenesis of PBC. Additionally, telomere dysregulation in biliary epithelial cells (BECs) may play a role in disease onset, whereas a deficiency in sex chromosome and skewed gene expression in the X chromosome may to some extent explain the female dominance in PBC.

Analysis of Four Polymorphisms Located at the Promoter of the Estrogen Receptor Alpha ESR1 Gene in a Population With Gender Incongruence

INTRODUCTION

Gender incongruence defines a state in which individuals feel discrepancy between the sex assigned at birth and their gender. Some of these people make a social transition from male to female (trans women) or from female to male (trans men). By contrast, the word cisgender describes a person whose gender identity is consistent with their sex assigned at birth.

AIM

To analyze the implication of the estrogen receptor α gene (ESR1) in the genetic basis of gender incongruence.

MAIN OUTCOME MEASURES

Polymorphisms rs9478245, rs3138774, rs2234693, rs9340799.

METHOD

We carried out the analysis of 4 polymorphisms located at the promoter of the ESR1 gene (C1 = rs9478245, C2 = rs3138774, C3 = rs2234693, and C4 = rs9340799) in a population of 273 trans women, 226 trans men, and 537 cis gender controls. For SNP polymorphisms, the allele and genotype frequencies were analyzed by χ² test. The strength of the SNP associations with gender incongruence was measured by binary logistic regression. For the STR polymorphism, the mean number of repeats were analyzed by the Mann-Whitney U test. Measurement of linkage disequilibrium and haplotype frequencies were also performed.

RESULTS

The C2 median repeats were shorter in the trans men population. Genotypes S/S and S/L for the C2 polymorphism were overrepresented in the trans men group (P = .012 and P = .003 respectively). We also found overtransmission of the A/A genotype (C4) in the trans men population (P = .017), while the A/G genotype (C4) was subrepresented (P = .009]. The analyzed polymorphisms were in linkage disequilibrium. In the trans men population, the T(C1)-L(C2)-C(C3)-A(C4) haplotype was overrepresented (P = .019) while the T(C1)-L(C2)-C(C3)-G(C4) was subrepresented (P = .005).

CONCLUSION

The ESR1 is associated with gender incongruence in the trans men population. Fernández R, Delgado-Zayas E,RamírezK, et al. Analysis of Four Polymorphisms Located at the Promoter of the Estrogen Receptor Alpha ESR1 Gene in a Population With Gender Incongruence. Sex Med 2020;8:490-500.

Hypertrophic cardiomyopathy MYH7 mutation R723G alters mRNA secondary structure

Point mutation R723G in the MYH7 gene causes hypertrophic cardiomyopathy (HCM). Heterozygous patients with this mutation exhibit a comparable allelic imbalance of the MYH7 gene. On average 67% of the total MYH7 mRNA are derived from the MYH7^R723G-allele and 33% from the MYH7^WT allele. Mechanisms underlying mRNA allelic imbalance are largely unknown. We suggest that a different mRNA lifetime of the alleles may cause the allelic drift in R723G patients. A potent regulator of mRNA lifetime is its secondary structure. To test for alterations in the MYH7^R723G mRNA structure we used selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis. We show significantly different SHAPE reactivity of wild-type and MYH7^R723G RNA, which is in accordance with bioinformatically predicted structures. Thus, we provide the first experimental evidence for mRNA secondary structure alterations by the HCM point mutation. We assume that this may result in a prolonged lifetime of MYH7^R723G mRNA in vivo and subsequently in the determined allelic imbalance.

Haplotype-Based Association and In Silico Studies of OPRM1 Gene Variants with Susceptibility to Opioid Dependence Among Addicted Iranians Undergoing Methadone Treatment

The associations of OPRM1 gene variants with opioid dependence have been demonstrated. This study investigated the association of rs495491, rs1799971 (A118G), rs589046, and rs10457090 variants of OPRM1 gene with opium dependence and their haplotypes among addicted individuals undergoing methadone treatment. Moreover, we investigated whether any of these variants were associated with libido dysfunction or insomnia among addicted people. A total of 404 individuals were genotyped by amplification refractory mutation system (ARMS) PCR. In silico studies were designed through homology modeling of A118G structures (N40 and D40) and docked with 41 FDA-approved drugs of OPRM1 protein by SWISS-MODEL, COACH, MolProbity, ProSA, Errat, Glide XP, and Autodock 4. Results revealed that rs495491, A118G, rs589046, and rs10457090 were significantly associated with opium dependence under recessive (P = 6.66E-10), dominant (P = 0.017), co-dominant (P = 0.001), and recessive (P = 9.28E-6) models of inheritance, respectively. Further analyses indicated three significant haplotypes including A-A-A-C (P-permutation < 1E-9), G-G-A-C (P-permutation = 0.04), and G-A-G-C (P-permutation = 8.69E-4). Genotype-phenotype associations of OPRM1 variants with insomnia and libido dysfunction showed no significant association. Docking showed the higher binding affinity of N40 rather than D40 model; however, methadone and morphine were bonded with D40 structure more powerful. Consequently, rs495491, A118G, rs589046, and rs10457090 were associated with opioid dependence among Iranians; also, A118G might be the most remarkable marker of OPRM1 owing to its vital structural roles.

Allele-specific expression identified rs2509956 as a novel long-distance cis-regulatory SNP for SCGB1A1, an important gene for multiple pulmonary diseases

SCGB1A1 (secretoglobin family 1A member 1) is an important protein for multiple pulmonary diseases, especially asthma, chronic obstructive pulmonary disease, and lung cancer. One single-nucleotide polymorphism (SNP) at 5'-untranslated region of SCGB1A1, rs3741240, has been suggested to be associated with reduced protein expression and further asthma susceptibility. However, it was still unclear whether there were other cis-regulatory elements for SCGB1A1 that might further contribute to pulmonary diseases. Allele-specific expression (ASE) is a novel approach to identify the functional region in human genome. In the present study, we measured ASE on rs3741240 in lung tissues and observed a consistent excess of G allele over A (P < 10^-6), which indicated that this SNP or the one(s) in linkage disequilibrium (LD) could regulate SCGB1A1 expression. By analyzing 1000 Genomes Project data for Chinese, one SNP locating ~10.2 kb away and downstream of SCGB1A1, rs2509956, was identified to be in strong LD with rs3741240. Reporter gene assay confirmed that both SNPs could regulate gene expression in the lung cell. By chromosome conformation capture, it was verified that the region surrounding rs2509956 could interact with SCGB1A1 promoter region and act as an enhancer. Through chromatin immunoprecipitation and overexpression assay, the related transcription factor RELA (RELA proto-oncogene, NF-kB subunit) was recognized to bind the region spanning rs2509956. Our work identified a novel long-distance cis-regulatory SNP for SCGB1A1, which might contribute to multiple pulmonary diseases.

Allelic Expression Imbalance Promoting a Mutant PEX6 Allele Causes Zellweger Spectrum Disorder

Zellweger spectrum disorders (ZSDs) are autosomal-recessive disorders that are caused by defects in peroxisome biogenesis due to bi-allelic mutations in any of 13 different PEX genes. Here, we identified seven unrelated individuals affected with an apparent dominant ZSD in whom a heterozygous mutant PEX6 allele (c.2578C>T [p.Arg860Trp]) was overrepresented due to allelic expression imbalance (AEI). We demonstrated that AEI of PEX6 is a common phenomenon and is correlated with heterozygosity for a frequent variant in the 3' untranslated region (UTR) of the mutant allele, which disrupts the most distal of two polyadenylation sites. Asymptomatic parents, who were heterozygous for PEX c.2578C>T, did not show AEI and were homozygous for the 3' UTR variant. Overexpression models confirmed that the overrepresentation of the pathogenic PEX6 c.2578T variant compared to wild-type PEX6 c.2578C results in a peroxisome biogenesis defect and thus constitutes the cause of disease in the affected individuals. AEI promoting the overrepresentation of a mutant allele might also play a role in other autosomal-recessive disorders, in which only one heterozygous pathogenic variant is identified.

An xQTL map integrates the genetic architecture of the human brain's transcriptome and epigenome

We report a multi-omic resource generated by applying quantitative trait locus (xQTL) analyses to RNA sequence, DNA methylation and histone acetylation data from the dorsolateral prefrontal cortex of 411 older adults who have all three data types. We identify SNPs significantly associated with gene expression, DNA methylation and histone modification levels. Many of these SNPs influence multiple molecular features, and we demonstrate that SNP effects on RNA expression are fully mediated by epigenetic features in 9% of these loci. Further, we illustrate the utility of our new resource, xQTL Serve, by using it to prioritize the cell type(s) most affected by an xQTL. We also reanalyze published genome wide association studies using an xQTL-weighted analysis approach and identify 18 new schizophrenia and 2 new bipolar susceptibility variants, which is more than double the number of loci that can be discovered with a larger blood-based expression eQTL resource.

Dopamine transporter (DAT1/SLC6A3) polymorphism and the association between being born small for gestational age and symptoms of ADHD

Being small for gestational age (SGA) has been established as a risk factor for Attention Deficit Hyperactivity Disorder (ADHD). Likewise, several molecular genetic studies have found a link between DAT1 and ADHD. This study investigated whether SGA moderates the effect of dopamine transporter gene variants on the risk of ADHD. A total of 546 children of European descent were genotyped at age 11 for seven DAT1 SNPs (rs6347, rs11564774, rs40184, rs1042098, rs2702, rs8179029 and rs3863145). The Strengths and Difficulties Questionnaire was used to measure symptoms of ADHD at ages 3.5, 7 and 11. We found significant gene-environment interactions between birth weight and DAT1 SNPs (rs6347, rs40184, rs1042098, rs3863145) on ADHD symptoms at 3.5 years only. Results suggest that genotypic variation of DAT1 may confer a relative protective effect against ADHD in SGA individuals. This study supports the idea that being born SGA moderates the effect of the DAT1 gene on ADHD symptoms in the preschool years and may help to explain some of the heterogeneity in ADHD outcomes.

Intronic SNP in ESR1 encoding human estrogen receptor alpha is associated with brain ESR1 mRNA isoform expression and behavioral traits

Genetic variants of ESR1 have been implicated in multiple diseases, including behavioral disorders, but causative variants remain uncertain. We have searched for regulatory variants affecting ESR1 expression in human brain, measuring allelic ESR1 mRNA expression in human brain tissues with marker SNPs in exon4 representing ESR1-008 (or ESRα-36), and in the 3'UTR of ESR1-203, two main ESR1 isoforms in brain. In prefrontal cortex from subjects with bipolar disorder, schizophrenia, and controls (n = 35 each; Stanley Foundation brain bank), allelic ESR1 mRNA ratios deviated from unity up to tenfold at the exon4 marker SNP, with large allelic ratios observed primarily in bipolar and schizophrenic subjects. SNP scanning and targeted sequencing identified rs2144025, associated with large allelic mRNA ratios (p = 1.6E10-6). Moreover, rs2144025 was significantly associated with ESR1 mRNA levels in the Brain eQTL Almanac and in brain regions in the Genotype-Tissue Expression project. In four GWAS cohorts, rs2104425 was significantly associated with behavioral traits, including: hypomanic episodes in female bipolar disorder subjects (GAIN bipolar disorder study; p = 0.0004), comorbid psychological symptoms in both males and females with attention deficit hyperactivity disorder (GAIN ADHD, p = 0.00002), psychological diagnoses in female children (eMERGE study of childhood health, subject age ≥9, p = 0.0009), and traits in schizophrenia (e.g., grandiose delusions, GAIN schizophrenia, p = 0.0004). The first common ESR1 variant (MAF 12-33% across races) linked to regulatory functions, rs2144025 appears conditionally to affect ESR1 mRNA expression in the brain and modulate traits in behavioral disorders.

Effect of Sodium Arsenite on the Expression of Antioxidant Genes (SOD2 and CAT) in MCF-7 and Jurkat Cell Lines

BACKGROUND

Sodium arsenite (NaAsO2) has potent cytotoxic activity in human cancer cells. Oxidative stress has been suggested as a mechanism for arsenic-induced carcinogenesis. The purpose of the present study was to evaluate the alteration of mRNA levels of catalase (CAT) and superoxide dismutase 2 (SOD2) in MCF-7 and Jurkat cells after exposure to NaAsO2.

METHODS

Methylthiazol tetrazolium (MTT) viability assay was performed to evaluate cytotoxicity of NaAsO2 in MCF-7 and Jurkat cells. For evaluating the expression levels of the CAT and SOD2, we used two concentrations of NaAsO2 (5 and 15 μM), lower than the concentrations at which 50% of cell viability were lost. The cells were treated with co-treatment of NaAsO2 (15 μM) and N-acetyl-cysteine (NAC; 5 μM) in the media for 24 h. The control cells were maintained in sodium arsenite free growth medium. The experiments were done triplicate. Using quantitative real-time PCR, the expression levels of CAT and SOD2 were quantified. One sample student's t test was performed for comparisons of mRNA levels between treatment groups and their corresponding untreated control cells.

RESULTS

CAT mRNA level decreased significantly in both cell lines following exposure to NaAsO2 (P<0.05). Expression levels of SOD2 decreased in Jurkat cells and increased in MCF-7 cells after treatment with NaAsO2 (P<0.05).

CONCLUSION

After cells exposure to NaAsO2, CAT mRNA level decreased in both examined cell lines but the alterations of SOD2 mRNA level is cell specific. The NAC modulated the NaAsO2 associated alterations of CAT and SOD2 mRNA levels, therefore, the NaAsO2 might act through inducing reactive oxygen species.

Polymorphisms of Dopamine Receptor Genes and Risk of L-Dopa-Induced Dyskinesia in Parkinson's Disease

L-dopa–induced dyskinesia (LID) is a frequent motor complication of Parkinson’s disease (PD), associated with a negative prognosis. Previous studies showed an association between dopamine receptor (DR) gene (DR) variants and LID, the results of which have not been confirmed. The present study is aimed to determine whether genetic differences of DR are associated with LID in a small but well-characterized cohort of PD patients. To this end we enrolled 100 PD subjects, 50 with and 50 without LID, matched for age, gender, disease duration and dopaminergic medication in a case-control study. We conducted polymerase chain reaction for single nucleotide polymorphisms (SNP) in both D1-like (DRD1A48G; DRD1C62T and DRD5T798C) and D2-like DR (DRD2G2137A, DRD2C957T, DRD3G25A, DRD3G712C, DRD4C616G and DRD4nR VNTR 48bp) analyzed genomic DNA. Our results showed that PD patients carrying allele A at DRD3G3127A had an increased risk of LID (OR 4.9; 95% CI 1.7–13.9; p = 0.004). The present findings may provide valuable information for personalizing pharmacological therapy in PD patients.

Rapporteur summaries of plenary, symposia, and oral sessions from the XXIIIrd World Congress of Psychiatric Genetics Meeting in Toronto, Canada, 16-20 October 2015

The XXIIIrd World Congress of Psychiatric Genetics meeting, sponsored by the International Society of Psychiatric Genetics, was held in Toronto, ON, Canada, on 16-20 October 2015. Approximately 700 participants attended to discuss the latest state-of-the-art findings in this rapidly advancing and evolving field. The following report was written by trainee travel awardees. Each was assigned one session as a rapporteur. This manuscript represents the highlights and topics that were covered in the plenary sessions, symposia, and oral sessions during the conference, and contains major notable and new findings.

The epigenetics of PBC: The link between genetic susceptibility and environment

Primary biliary cholangitis (PBC) previously known as primary biliary cirrhosis is an autoimmune disease-associated with progressive cholestasis, the presence of autoreactive T cell and characteristic serological autoantibodies. Genetic and genome-wide association studies (GWAS) have recently shed light on the genetic background of PBC. Besides that some causal nucleotide changes and mechanisms remain largely unknown as suggested for example, by the observation that monozygotic twins have an identical DNA sequence even if presents some phenotypic differences that may be consequences of different exposures to environmental stressors. For this reason, it is believed that epigenetic mechanisms may be involved in PBC pathogenesis, as already demonstrated in many autoimmune diseases and can eventually provide an understanding that has been missed from genetics alone. This review will focus on the most commonly studied epigenetic modifications already demonstrated in PBC; special attention will be paid also to other epigenetic mechanisms so far not demonstrated in PBC patients, but that could increase our understanding in PBC pathogenesis.

The genetic architecture of autism spectrum disorders (ASDs) and the potential importance of common regulatory genetic variants

Currently, there is great interest in identifying genetic variants that contribute to the risk of developing autism spectrum disorders (ASDs), due in part to recent increases in the frequency of diagnosis of these disorders worldwide. While there is nearly universal agreement that ASDs are complex diseases, with multiple genetic and environmental contributing factors, there is less agreement concerning the relative importance of common vs rare genetic variants in ASD liability. Recent observations that rare mutations and copy number variants (CNVs) are frequently associated with ASDs, combined with reduced fecundity of individuals with these disorders, has led to the hypothesis that ASDs are caused primarily by de novo or rare genetic mutations. Based on this model, large-scale whole-genome DNA sequencing has been proposed as the most appropriate method for discovering ASD liability genes. While this approach will undoubtedly identify many novel candidate genes and produce important new insights concerning the genetic causes of these disorders, a full accounting of the genetics of ASDs will be incomplete absent an understanding of the contributions of common regulatory variants, which are likely to influence ASD liability by modifying the effects of rare variants or, by assuming unfavorable combinations, directly produce these disorders. Because it is not yet possible to identify regulatory genetic variants by examination of DNA sequences alone, their identification will require experimentation. In this essay, I discuss these issues and describe the advantages of measurements of allelic expression imbalance (AEI) of mRNA expression for identifying cis-acting regulatory variants that contribute to ASDs.

Polymorphisms of dopamine receptor genes and risk of visual hallucinations in Parkinson's patients

BACKGROUND

Visual hallucinations (VHs) are frequent non-motor complication of Parkinson's disease (PD), associated to a negative prognosis. Previous studies showed an association between dopamine receptor (DR) gene (DR) variants and psychosis in Alzheimer's disease, addictions, schizophrenia, and bipolar disorder. However, there are only a few studies on DR variants and VHs in PD, which did not provide conclusive results.

OBJECTIVES

The present study aimed to determine whether genetic differences of DR are associated with visual hallucinations (VHs) in a cohort of Parkinson's disease (PD) patients.

METHODS

A case-control study of 84 PD subjects, 42 with and 42 without VHs,that were matched for age, gender, disease duration, and dopaminergic medication was conducted. Polymerase chain reaction for SNPs in both D1-like (DRD1A-48G [rs4532] and C62T [rs686], DRD5T798C [rs6283]) and D2-like DR (DRD2G2137A [rs1800497] and C957T [rs6277], DRD3G25A [rs6280] and G712C [rs1800828], DRD4C616G [rs747302] and nR VNTR 48bp) analyzed genomic DNA.

RESULTS

Patients carrying allele T at DRD1C62T had an increased risk of VHs, expressed as OR (95 % CI, p value), of 10.7 (2.9-40, p = 0.0001). Moreover, patients with DRD1-48 GG and 62TT genotype displayed shorter time to VHs, whereas a longer time to VHs was found in subjects carrying the DRD4 CG alleles.

CONCLUSIONS

PD patients with VHs display higher frequency of DR SNPs associated with increased D1-like activity and decreased D2-like activity. Our data are in line with associations reported in other neurodegenerative and psychiatric conditions. Results likely provide valuable information for personalizing pharmacological therapy in PD patients.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Uncovering drug-responsive regulatory elements

Nucleotide changes in gene regulatory elements can have a major effect on interindividual differences in drug response. For example, by reviewing all published pharmacogenomic genome-wide association studies, we show here that 96.4% of the associated single nucleotide polymorphisms reside in noncoding regions. We discuss how sequencing technologies are improving our ability to identify drug response-associated regulatory elements genome-wide and to annotate nucleotide variants within them. We highlight specific examples of how nucleotide changes in these elements can affect drug response and illustrate the techniques used to find them and functionally characterize them. Finally, we also discuss challenges in the field of drug-responsive regulatory elements that need to be considered in order to translate these findings into the clinic.

Allele-Selective Transcriptome Recruitment to Polysomes Primed for Translation: Protein-Coding and Noncoding RNAs, and RNA Isoforms

mRNA translation into proteins is highly regulated, but the role of mRNA isoforms, noncoding RNAs (ncRNAs), and genetic variants remains poorly understood. mRNA levels on polysomes have been shown to correlate well with expressed protein levels, pointing to polysomal loading as a critical factor. To study regulation and genetic factors of protein translation we measured levels and allelic ratios of mRNAs and ncRNAs (including microRNAs) in lymphoblast cell lines (LCL) and in polysomal fractions. We first used targeted assays to measure polysomal loading of mRNA alleles, confirming reported genetic effects on translation of OPRM1 and NAT1, and detecting no effect of rs1045642 (3435C>T) in ABCB1 (MDR1) on polysomal loading while supporting previous results showing increased mRNA turnover of the 3435T allele. Use of high-throughput sequencing of complete transcript profiles (RNA-Seq) in three LCLs revealed significant differences in polysomal loading of individual RNA classes and isoforms. Correlated polysomal distribution between protein-coding and non-coding RNAs suggests interactions between them. Allele-selective polysome recruitment revealed strong genetic influence for multiple RNAs, attributable either to differential expression of RNA isoforms or to differential loading onto polysomes, the latter defining a direct genetic effect on translation. Genes identified by different allelic RNA ratios between cytosol and polysomes were enriched with published expression quantitative trait loci (eQTLs) affecting RNA functions, and associations with clinical phenotypes. Polysomal RNA-Seq combined with allelic ratio analysis provides a powerful approach to study polysomal RNA recruitment and regulatory variants affecting protein translation.

Regulation of cholesteryl ester transfer protein expression by upstream polymorphisms: reduced expression associated with rs247616

BACKGROUND

Cholesteryl ester transfer protein (CETP) is involved in reverse cholesterol transport by exchanging cholesteryl esters for triglycerides between high-density lipoprotein and low-density lipoprotein particles, effectively decreasing high-density lipoprotein cholesterol levels. Variants within a large haplotype block upstream of CETP (rs247616, rs173539) have been shown to be significantly associated with reduced expression; however, the underlying mechanism has not been identified.

METHODS

We analyzed the linkage structure of our top candidate single-nucleotide polymorphism (SNP), rs247616, and assessed each SNP of the haplotype block for potential interactions with transcription factor binding sites. We then used a reporter gene assay to assess the effect of three SNPs (rs247616, rs173539, and rs1723150) on expression in vitro.

RESULTS

Several variants in the upstream haplotype, including rs247616, rs173539, and rs1723150, disrupt or generate transcription factor binding sites. In reporter gene assays, rs247616 and rs173539 were found to significantly affected expression in HepG2 cells, whereas rs17231506 had no effect. rs247616 decreased expression by 1.7-fold (P<0.0001), whereas rs173539 increased expression by 2.2-fold (P=0.0006).

CONCLUSION

SNPs rs247616 and rs173539 are in high linkage disequilibrium (R=0.96, D'=1.00) and have the potential to regulate CETP expression. Although opposing effects suggest that regulation of CETP expression could vary between tissues, the minor allele of rs247616 and SNPs in high linkage with it were found to be associated with reduced expression across all tissues.

MDR1 gene polymorphisms and imatinib response in chronic myeloid leukemia: a meta-analysis

BACKGROUND

MDR1 gene polymorphisms were demonstrated to be associated with interindividual variability of imatinib response for chronic myeloid leukemia (CML) patients in several studies; however, the results have been inconclusive.

MATERIALS & METHODS

To clarify the effect of common MDR1 variants on clinical response to imatinib, we performed a meta-analysis to quantify the accumulated information from genetic association studies. After a thorough search of the published literature, we undertook a meta-analysis to evaluate the effect of MDR1 C1236T, G2677T and C3435T polymorphisms on imatinib response.

RESULTS

Our pooled data showed a significant association between MDR1 C1236T polymorphism and the increasing risk of imatinib resistance in Asian CML patients. However, no significant association was found for the MDR1 G2677T or C3435T polymorphisms in an Asian CML population as well as a Caucasian CML population.

CONCLUSION

The synonymous MDR1 C1236T polymorphism might be a risk factor for nonoptimal clinical response to imatinib in Asian CML patients.

Endothelial nitric oxide synthase gene expression is associated with hypertension in autosomal dominant polycystic kidney disease

BACKGROUND/AIMS

Early occurrence of hypertension is the prominent feature of autosomal dominant polycystic kidney disease (ADPKD). The role of angiotensin-converting enzyme (ACE) gene polymorphism and endothelial nitric oxide synthase (eNOS) gene polymorphism in the clinical course of ADPKD is not well understood. However, data about the expression of these genes are lacking. Thus, we aimed to investigate the polymorphisms and expressions of both the ACE and eNOS genes that affect hypertension in ADPKD.

METHODS

Whole blood samples were obtained from all participants. ACE and eNOS gene polymorphisms and their expressions were analyzed in 78 ADPKD patients and 30 controls. Gene expressions were assessed by quantitative real-time PCR. Twenty-four-hour blood pressure monitoring was performed for the diagnosis of hypertension in all study participants.

RESULTS

eNOS expression and the estimated glomerular filtration rate were found to be significantly higher in ADPKD patients without hypertension than in those with hypertension. Each unit of increase in eNOS expression led to a 0.88-fold decrease (95% CI: 0.80-0.96) in the risk of hypertension in multiple logistic regression analysis.

CONCLUSIONS

eNOS gene expression is independently predictive of hypertension in the ADPKD population. This study showed, for the first time, a novel link between eNOS gene expression and hypertension in ADPKD.

Missing heritability of common diseases and treatments outside the protein-coding exome

Genetic factors strongly influence risk of common human diseases and treatment outcomes but the causative variants remain largely unknown; this gap has been called the 'missing heritability'. We propose several hypotheses that in combination have the potential to narrow the gap. First, given a multi-stage path from wellness to disease, we propose that common variants under positive evolutionary selection represent normal variation and gate the transition between wellness and an 'off-well' state, revealing adaptations to changing environmental conditions. In contrast, genome-wide association studies (GWAS) focus on deleterious variants conveying disease risk, accelerating the path from off-well to illness and finally specific diseases, while common 'normal' variants remain hidden in the noise. Second, epistasis (dynamic gene-gene interactions) likely assumes a central role in adaptations and evolution; yet, GWAS analyses currently are poorly designed to reveal epistasis. As gene regulation is germane to adaptation, we propose that epistasis among common normal regulatory variants, or between common variants and less frequent deleterious variants, can have strong protective or deleterious phenotypic effects. These gene-gene interactions can be highly sensitive to environmental stimuli and could account for large differences in drug response between individuals. Residing largely outside the protein-coding exome, common regulatory variants affect either transcription of coding and non-coding RNAs (regulatory SNPs, or rSNPs) or RNA functions and processing (structural RNA SNPs, or srSNPs). Third, with the vast majority of causative variants yet to be discovered, GWAS rely on surrogate markers, a confounding factor aggravated by the presence of more than one causative variant per gene and by epistasis. We propose that the confluence of these factors may be responsible to large extent for the observed heritability gap.

Association of OPRM1 A118G variant with risk of morphine-induced respiratory depression following spine fusion in adolescents

The μ1 opioid receptor (OPRM1) genetic variant A118G results in decreased μ-receptor binding potential in the brain and increases morphine requirement. We hypothesized that OPRM1 A118G polymorphism will affect morphine-induced respiratory depression (MIRD) risk in children receiving morphine. A prospective genotype-blinded study was conducted in 88 healthy adolescents (11-18 years; 67% female, 85% Caucasian) who underwent spine fusion for scoliosis. They were followed for 48 h postoperatively for MIRD, pain scores, morphine consumption and use of analgesic adjuvants. Patients were genotyped for OPRM1 A118G variant-76% were wild type (AA) and 24% heterozygous/homozygous for variant (AG/GG). Multivariable logistic regression showed that the risk of MIRD in patients with AA genotype was significantly higher (odds ratio 5.6, 95% CI: 1.4-37.2, P=0.030). Presence of G allele was associated with higher pain scores (effect size 0.73, P=0.045). This novel association is an important step toward predicting MIRD susceptibility and personalizing morphine use.

Synthesis of 53 tissue and cell line expression QTL datasets reveals master eQTLs

BACKGROUND

Gene expression genetic studies in human tissues and cells identify cis- and trans-acting expression quantitative trait loci (eQTLs). These eQTLs provide insights into regulatory mechanisms underlying disease risk. However, few studies systematically characterized eQTL results across cell and tissues types. We synthesized eQTL results from >50 datasets, including new primary data from human brain, peripheral plaque and kidney samples, in order to discover features of human eQTLs.

RESULTS

We find a substantial number of robust cis-eQTLs and far fewer trans-eQTLs consistent across tissues. Analysis of 45 full human GWAS scans indicates eQTLs are enriched overall, and above nSNPs, among positive statistical signals in genetic mapping studies, and account for a significant fraction of the strongest human trait effects. Expression QTLs are enriched for gene centricity, higher population allele frequencies, in housekeeping genes, and for coincidence with regulatory features, though there is little evidence of 5' or 3' positional bias. Several regulatory categories are not enriched including microRNAs and their predicted binding sites and long, intergenic non-coding RNAs. Among the most tissue-ubiquitous cis-eQTLs, there is enrichment for genes involved in xenobiotic metabolism and mitochondrial function, suggesting these eQTLs may have adaptive origins. Several strong eQTLs (CDK5RAP2, NBPFs) coincide with regions of reported human lineage selection. The intersection of new kidney and plaque eQTLs with related GWAS suggest possible gene prioritization. For example, butyrophilins are now linked to arterial pathogenesis via multiple genetic and expression studies. Expression QTL and GWAS results are made available as a community resource through the NHLBI GRASP database [http://apps.nhlbi.nih.gov/grasp/].

CONCLUSIONS

Expression QTLs inform the interpretation of human trait variability, and may account for a greater fraction of phenotypic variability than protein-coding variants. The synthesis of available tissue eQTL data highlights many strong cis-eQTLs that may have important biologic roles and could serve as positive controls in future studies. Our results indicate some strong tissue-ubiquitous eQTLs may have adaptive origins in humans. Efforts to expand the genetic, splicing and tissue coverage of known eQTLs will provide further insights into human gene regulation.

Genetic associations with expression for genes implicated in GWAS studies for atherosclerotic cardiovascular disease and blood phenotypes

Genome-wide association studies (GWAS) have uncovered many genetic associations for cardiovascular disease (CVD). However, data are limited regarding causal genetic variants within implicated loci. We sought to identify regulatory variants (cis- and trans-eQTLs) affecting expression levels of 93 genes selected by their proximity to SNPs with significant associations in prior GWAS for CVD traits. Expression levels were measured by qRT-PCR in leukocytes from 1846 Framingham Heart Study participants. An additive genetic model was applied to 2.5 million imputed SNPs for each gene. Approximately 45% of genes (N = 38) harbored at least one cis-eSNP after a regional multiple-test adjustment. Applying a more rigorous significance threshold (P < 5 × 10(-8)), we found the expression level of 10 genes was significantly associated with more than one cis-eSNP. The top cis-eSNPs for 7 of these 10 genes exhibited moderate-to-strong association with ≥ 1 CVD clinical phenotypes. Several eSNPs or proxy SNPs (r(2) = 1) were replicated by other eQTL studies. After adjusting for the lead GWAS SNPs for the 10 genes, expression variances explained by top cis-eSNPs were attenuated markedly for LPL, FADS2 and C6orf184, suggesting a shared genetic basis for the GWAS and expression trait. A significant association between cis-eSNPs, gene expression and lipid levels was discovered for LPL and C6orf184. In conclusion, strong cis-acting variants are localized within nearly half of the GWAS loci studied, with particularly strong evidence for a regulatory role of the top GWAS SNP for expression of LPL, FADS2 and C6orf184.

Whole transcriptome RNA-Seq allelic expression in human brain

Background

Measuring allelic RNA expression ratios is a powerful approach for detecting cis-acting regulatory variants, RNA editing, loss of heterozygosity in cancer, copy number variation, and allele-specific epigenetic gene silencing. Whole transcriptome RNA sequencing (RNA-Seq) has emerged as a genome-wide tool for identifying allelic expression imbalance (AEI), but numerous factors bias allelic RNA ratio measurements. Here, we compare RNA-Seq allelic ratios measured in nine different human brain regions with a highly sensitive and accurate SNaPshot measure of allelic RNA ratios, identifying factors affecting reliable allelic ratio measurement. Accounting for these factors, we subsequently surveyed the variability of RNA editing across brain regions and across individuals.

Results

We find that RNA-Seq allelic ratios from standard alignment methods correlate poorly with SNaPshot, but applying alternative alignment strategies and correcting for observed biases significantly improves correlations. Deploying these methods on a transcriptome-wide basis in nine brain regions from a single individual, we identified genes with AEI across all regions (SLC1A3, NHP2L1) and many others with region-specific AEI. In dorsolateral prefrontal cortex (DLPFC) tissues from 14 individuals, we found evidence for frequent regulatory variants affecting RNA expression in tens to hundreds of genes, depending on stringency for assigning AEI. Further, we find that the extent and variability of RNA editing is similar across brain regions and across individuals.

Conclusions

These results identify critical factors affecting allelic ratios measured by RNA-Seq and provide a foundation for using this technology to screen allelic RNA expression on a transcriptome-wide basis. Using this technology as a screening tool reveals tens to hundreds of genes harboring frequent functional variants affecting RNA expression in the human brain. With respect to RNA editing, the similarities within and between individuals leads us to conclude that this post-transcriptional process is under heavy regulatory influence to maintain an optimal degree of editing for normal biological function.

Consequences of the 118A>G polymorphism in the OPRM1 gene: translation from bench to bedside?

The 118A>G single nucleotide polymorphism (SNP) in the μ-opioid receptor (OPRM1) gene has been the most described variant in pharmacogenetic studies regarding opioid drugs. Despite evidence for an altered biological function encoded by this variant, this knowledge is not yet utilized clinically. The aim of the present review was to collect and discuss the available information on the 118A>G SNP in the OPRM1 gene, at the molecular level and in its clinical manifestations. In vitro biochemical and molecular assays have shown that the variant receptor has higher binding affinity for β-endorphins, that it has altered signal transduction cascade, and that it has a lower expression compared with wild-type OPRM1. Studies using animal models for 118A>G have revealed a double effect of the variant receptor, with an apparent gain of function with respect to the response to endogenous opioids but a loss of function with exogenous administered opioid drugs. Although patients with this variant have shown a lower pain threshold and a higher drug consumption in order to achieve the analgesic effect, clinical experiences have demonstrated that patients carrying the variant allele are not affected by the increased opioid consumption in terms of side effects.

Allelic expression analysis of the osteoarthritis susceptibility gene COL11A1 in human joint tissues

Background

The single nucleotide polymorphism (SNP) rs2615977 is associated with osteoarthritis (OA) and is located in intron 31 of COL11A1, a strong candidate gene for this degenerative musculoskeletal disease. Furthermore, the common non-synonymous COL11A1 SNP rs1676486 is associated with another degenerative musculoskeletal disease, lumbar disc herniation (LDH). rs1676486 is a C-T transition mediating its affect on LDH susceptibility by modulating COL11A1 expression. The risk T-allele of rs1676486 leads to reduced expression of the COL11A1 transcript, a phenomenon known as allelic expression imbalance (AEI). We were keen therefore to assess whether the effect that rs1676486 has on COL11A1 expression in LDH is also observed in OA and whether the rs2615977 association to OA also marked AEI.

Methods

Using RNA from OA cartilage, we assessed whether either SNP correlated with COL11A1 AEI by 1) measuring COL11A1 expression and stratifying the data by genotype at each SNP; and 2) quantifying the mRNA transcribed from each allele of the two SNPs. We also assessed whether rs1676486 was associated with OA susceptibility using a case–control cohort of over 18,000 individuals.

Results

We observed significant AEI at rs1676486 (p < 0.0001) with the T-allele correlating with reduced COL11A1 expression. This corresponded with observations in LDH but the SNP was not associated with OA. We did not observe AEI at rs2615977.

Conclusions

COL11A1 is subject to AEI in OA cartilage. AEI at rs1676486 is a risk factor for LDH, but not for OA. These two diseases therefore share a common functional phenotype, namely AEI of COL11A1, but this appears to be a disease risk only in LDH. Other functional effects on COL11A1 presumably account for the OA susceptibility that maps to this gene.