Genetic variants in nuclear-encoded mitochondrial genes influence AIDS progression

Maternal and fetal mitochondrial gene dysregulation in hypertensive disorders of pregnancy

Mitochondrial dysfunction has been implicated in pregnancy-induced hypertension (PIH). The role of mitochondrial gene dysregulation in PIH, and consequences for maternal-fetal interactions, remain elusive. Here, we investigated mitochondrial gene expression and dysregulation in maternal and placental tissues from pregnancies with and without PIH; further, we measured circulating mitochondrial DNA (mtDNA) mutational load, an index of mtDNA integrity. Differential gene expression analysis followed by Time Course Gene Set Analysis (TcGSA) was conducted on publicly available high throughput sequencing transcriptomic data sets. Mutational load analysis was carried out on peripheral mononuclear blood cells from healthy pregnant individuals and individuals with preeclampsia. Thirty mitochondrial differentially expressed genes (mtDEGs) were detected in the maternal cell-free circulating transcriptome, whereas nine were detected in placental transcriptome from pregnancies with PIH. In PIH pregnancies, maternal mitochondrial dysregulation was associated with pathways involved in inflammation, cell death/survival, and placental development, whereas fetal mitochondrial dysregulation was associated with increased production of extracellular vesicles (EVs) at term. Mothers with preeclampsia did not exhibit a significantly different degree of mtDNA mutational load. Our findings support the involvement of maternal mitochondrial dysregulation in the pathophysiology of PIH and suggest that mitochondria may mediate maternal-fetal interactions during healthy pregnancy.NEW & NOTEWORTHY This study identifies aberrant maternal and fetal expression of mitochondrial genes in pregnancies with gestational hypertension and preeclampsia. Mitochondrial gene dysregulation may be a common etiological factor contributing to the development of de novo hypertension in pregnancy-associated hypertensive disorders.

Identification of Differential Expression Genes between Volume and Pressure Overloaded Hearts Based on Bioinformatics Analysis

Volume overload (VO) and pressure overload (PO) are two common pathophysiological conditions associated with cardiac disease. VO, in particular, often occurs in a number of diseases, and no clinically meaningful molecular marker has yet been established. We intend to find the main differential gene expression using bioinformatics analysis. GSE97363 and GSE52796 are the two gene expression array datasets related with VO and PO, respectively. The LIMMA algorithm was used to identify differentially expressed genes (DEGs) of VO and PO. The DEGs were divided into three groups and subjected to functional enrichment analysis, which comprised GO analysis, KEGG analysis, and the protein–protein interaction (PPI) network. To validate the sequencing data, cardiomyocytes from AR and TAC mouse models were used to extract RNA for qRT-PCR. The three genes with random absolute values of LogFC and indicators of heart failure (natriuretic peptide B, NPPB) were detected: carboxylesterase 1D (CES1D), whirlin (WHRN), and WNK lysine deficient protein kinase 2 (WNK2). The DEGs in VO and PO were determined to be 2761 and 1093, respectively, in this study. Following the intersection, 305 genes were obtained, 255 of which expressed the opposing regulation and 50 of which expressed the same regulation. According to the GO and pathway enrichment studies, DEGs with opposing regulation are mostly common in fatty acid degradation, propanoate metabolism, and other signaling pathways. Finally, we used Cytoscape’s three techniques to identify six hub genes by intersecting 255 with the opposite expression and constructing a PPI network. Peroxisome proliferator-activated receptor (PPARα), acyl-CoA dehydrogenase medium chain (ACADM), patatin-like phospholipase domain containing 2 (PNPLA2), isocitrate dehydrogenase 3 (IDH3), heat shock protein family D member 1 (HSPD1), and dihydrolipoamide S-acetyltransferase (DLAT) were identified as six potential genes. Furthermore, we predict that the hub genes PPARα, ACADM, and PNPLA2 regulate VO myocardial changes via fatty acid metabolism and acyl-Coa dehydrogenase activity, and that these genes could be employed as basic biomarkers for VO diagnosis and treatment.

A novel prognostic two-gene signature for triple negative breast cancer

The absence of a robust risk stratification tool for triple negative breast cancer (TNBC) underlies imprecise and nonselective treatment of these patients with cytotoxic chemotherapy. This study aimed to interrogate transcriptomes of TNBC resected samples using next generation sequencing to identify novel biomarkers associated with disease outcomes. A subset of cases (n = 112) from a large, well-characterized cohort of primary TNBC (n = 333) were subjected to RNA-sequencing. Reads were aligned to the human reference genome (GRCH38.83) using the STAR aligner and gene expression quantified using HTSEQ. We identified genes associated with distant metastasis-free survival and breast cancer-specific survival by applying supervised artificial neural network analysis with gene selection to the RNA-sequencing data. The prognostic ability of these genes was validated using the Breast Cancer Gene-Expression Miner v4. 0 and Genotype 2 outcome datasets. Multivariate Cox regression analysis identified a prognostic gene signature that was independently associated with poor prognosis. Finally, we corroborated our results from the two-gene prognostic signature by their protein expression using immunohistochemistry. Artificial neural network identified two gene panels that strongly predicted distant metastasis-free survival and breast cancer-specific survival. Univariate Cox regression analysis of 21 genes common to both panels revealed that the expression level of eight genes was independently associated with poor prognosis (p < 0.05). Adjusting for clinicopathological factors including patient's age, grade, nodal stage, tumor size, and lymphovascular invasion using multivariate Cox regression analysis yielded a two-gene prognostic signature (ACSM4 and SPDYC), which was associated with poor prognosis (p < 0.05) independent of other prognostic variables. We validated the protein expression of these two genes, and it was significantly associated with patient outcome in both independent and combined manner (p < 0.05). Our study identifies a prognostic gene signature that can predict prognosis in TNBC patients and could potentially be used to guide the clinical management of TNBC patients.

Rat liver epigenome programing by perinatal exposure to 2,2',4'4'-tetrabromodiphenyl ether

Perinatal exposures to polybrominated diphenyl ethers permanently reprogram liver metabolism and induce a nonalcoholic fatty liver disease-like phenotype and insulin resistance in rodents. Aim: To test if these changes are associated with altered liver epigenome. Materials & methods: Expression of small RNA and changes in DNA methylation in livers of adult rats were analyzed following perinatal exposure to 2,2',4,4'-tetrabromodiphenyl ether, the polybrominated diphenyl ether congener most prevalent in human tissues. Results: We identified 33 differentially methylated DNA regions and 15 differentially expressed miRNAs. These changes were enriched for terms related to lipid and carbohydrate metabolism, insulin signaling, Type-2 diabetes and nonalcoholic fatty liver disease. Conclusion: Changes in the liver epigenome are a likely candidate mechanism of long-term maintenance of an aberrant metabolic phenotype.

Identification of Noncanonical Transcripts Produced by Systematic Nucleotide Exchanges in HIV-Associated Centroblastic Lymphoma

Noncanonical transcriptions include transcriptions that systematically exchange nucleotides, also called bijective transformations or swinger transformations. Swinger transformation A↔T+C↔G recovers identities of 8 among 9 unknown RNAs differentially expressed in centroblastic lymphoma, a human immunodeficiency virus (HIV)-associated non-Hodgkin's lymphoma. The identified RNAs align with human genes with known anti-HIV1 or oncogenic activities. Function disruption through swinger-transformed transcription potentially enables avoiding antiviral responses and contributes to cancer induction.

A new risk score based on twelve hepatocellular carcinoma-specific gene expression can predict the patients' prognosis

A large panel of molecular biomarkers have been identified to predict the prognosis of hepatocellular carcinoma (HCC), yet with limited clinical application due to difficult extrapolation. We here generated a genetic risk score system comprised of 12 HCC-specific genes to better predict the prognosis of HCC patients. Four genomics profiling datasets (GSE5851, GSE28691, GSE15765 and GSE14323) were searched to seek HCC-specific genes by comparisons between cancer samples and normal liver tissues and between different subtypes of hepatic neoplasms. Univariate survival analysis screened HCC-specific genes associated with overall survival (OS) in the training dataset for next-step risk model construction. The prognostic value of the constructed HCC risk score system was then validated in the TCGA dataset. Stratified analysis indicated this scoring system showed better performance in elderly male patients with HBV infection and preoperative lower levels of creatinine, alpha-fetoprotein and platelet and higher level of albumin. Functional annotation of this risk model in high-risk patients revealed that pathways associated with cell cycle, cell migration and inflammation were significantly enriched. In summary, our constructed HCC-specific gene risk model demonstrated robustness and potentiality in predicting the prognosis of HCC patients, especially among elderly male patients with HBV infection and relatively better general conditions.

Characterization and Analysis of Mammalian AKR7A Gene Promoters: Implications for Transcriptional Regulation

Aldo-keto reductase (AKR) superfamily is responsible for preventing mammalian cells from the toxic and carcinogenic effect of different genotoxic and non-genotoxic chemicals by reducing them, though the inducibility of these genes are different in different species. The aim of this paper is to compare the gene regulation mechanisms of AKR superfamily genes in different species and to identify the conserved areas, which are responsible for gene regulations in the presence of antioxidant, toxicants, and non-genotoxic carcinogens. At the beginning of the analysis AKR genes found in different species were divided into two groups based on their amino acid sequence similarities. Comparison of AKR7A gene clusters between different species revealed that Human AKR7A2 has orthologues in mammalians like rat, mouse, pigs, and other primates. On the other hand, AKR7A3 has orthologues only in rat and AKR7L is present only in primates. All the genes of AKR superfamily have a trend to stay in clusters in mammalian chromosomes having repeated sequences in between them. Transcription start site analysis revealed that genes like human AKR7A2 and rat Akr7a4 do not have conventional promoter regions such as TATA box, CAAT box and have several GC-rich regions, whereas gene like Akr7a1 contains a TATA box 25 bp upstream of transcription start site instead of having CpG islands. Putative orthologous genes i.e., rat AKR7A4, human AKR7A2, and mouse AKR7A5 share more common features such as common transcription factor binding site for specificity protein 1 (SP1), GATA binding factor family, Selenocysteine tRNA gene transcription activating factor (STAF) zinc finger protein, Krüppel-like C2H2 zinc finger (HICF) protein, negative glucocorticoid response element (NGRE) etc. Similarly, genes like rat AKR7A1, human AKR7A3, and human AKR7L share common sequence and transcription factor binding sites. Among those, Nuclear factor erythroid 2-related factor 2 (Nrf2) is thought to be responsible for the inducibility of these genes in the presence of antioxidants. Our analysis revealed that AKR7A gene family consists of genes having a large number of variations in them. Some of these, such as AKR7A2 are housekeeping genes, on the other hand, genes like AKR7A3 are highly inducible in the presence of antioxidants because of the presence of Nrf2 binding site in their promoter. AKR7A1 has a different promoter than others and function of AKR7L gene is still unknown.

Genetic and clinical predictors of CD4 lymphocyte recovery during suppressive antiretroviral therapy: Whole exome sequencing and antiretroviral therapy response phenotypes

Increase of peripheral blood CD4 lymphocyte counts is a key goal of combined antiretroviral therapy (cART); most, but not all, recipients respond adequately and promptly. A small number of studies have examined specific genetic factors associated with the extent of CD4 recovery. We report a genome-wide examination of factors that predict CD4 recovery in HIV-infected women. We identified women in in a cohort study who were on cART with viral load below 400 copies, and drew racially and ethnically matched samples of those with good CD4 response over 2 years or poor response. We analyzed the exomes of those women employing next generation sequencing for genes associated with CD4 recovery after controlling for non-genetic factors identified through forward stepwise selection as important. We studied 48 women with good CD4 recovery and 42 with poor CD4 recovery during virologically-suppressive cART. Stepwise logistic regression selected only age as a statistically significant (p<0.05) non-genetic predictor of response type (each additional year of age reduced the odds of good recovery by 11% (OR = 0.89, CI = 0.84–0.96, p = 0.0009). After adjustment for age and genomic estimates of race and ethnicity, 41 genes harbored variations associated with CD4 recovery group (p≤0.001); 5 of these have been previously reported to be associated with HIV infection, 4 genes would likely influence CD4 homeostasis, and 13 genes either had known functions or were members of product families that had functions for which interactions with HIV or effects on lymphocyte homeostasis were biologically plausible. Greater age was the strongest acquired factor that predicted poor CD4 cell recovery. Sequence variations spanning 41 genes were independently predictive of CD4 recovery. Many of these genes have functions that impact the cell cycle, apoptosis, lymphocyte migration, or have known interactions with HIV. These findings may help inform new hypotheses related to responses to HIV therapy and CD4 lymphocyte homeostasis.

Overview of mitochondrial germline variants and mutations in human disease: Focus on breast cancer (Review)

High lactate production in cells during growth under oxygen-rich conditions (aerobic glycolysis) is a hallmark of tumor cells, indicating the role of mitochondrial function in tumorigenesis. In fact, enhanced mitochondrial biogenesis and impaired quality control are frequently observed in cancer cells. Mitochondrial DNA (mtDNA) encodes 13 subunits of oxidative phosphorylation (OXPHOS), is present in thousands of copies per cell, and has a very high mutation rate. Mutations in mtDNA and nuclear DNA (nDNA) genes encoding proteins that are important players in mitochondrial biogenesis and function are involved in oncogenic processes. A wide range of germline mtDNA polymorphisms, as well as tumor mtDNA somatic mutations have been identified in diverse cancer types. Approximately 72% of supposed tumor-specific somatic mtDNA mutations reported, have also been found as polymorphisms in the general population. The ATPase 6 and NADH dehydrogenase subunit genes of mtDNA are the most commonly mutated genes in breast cancer (BC). Furthermore, nuclear genes playing a role in mitochondrial biogenesis and function, such as peroxisome proliferators-activated receptor gamma coactivator-1 (PGC-1), fumarate hydratase (FH) and succinate dehydrogenase (SDH) are frequently mutated in cancer. In this review, we provide an overview of the mitochondrial germline variants and mutations in cancer, with particular focus on those found in BC.

Analyses of the genetic diversity and protein expression variation of the acyl: CoA medium-chain ligases, ACSM2A and ACSM2B

Benzoate (found in milk and widely used as preservative), salicylate (present in fruits and the active component of aspirin), dietary polyphenols produced by gut microbiota, metabolites from organic acidemias, and medium-chain fatty acids (MCFAs) are all metabolised/detoxified by the glycine conjugation pathway. Xenobiotics are first activated to an acyl-CoA by the mitochondrial xenobiotic/medium-chain fatty acid: CoA ligases (ACSMs) and subsequently conjugated to glycine by glycine N-acyltransferase (GLYAT). The MCFAs are activated to acyl-CoA by the ACSMs before entering mitochondrial β-oxidation. This two-step enzymatic pathway has, however, not been thoroughly investigated and the biggest gap in the literature remains the fact that studies continuously characterise the pathway as a one-step reaction. There are no studies available on the interaction/competition of the various substrates involved in the pathway, whilst very little research has been done on the ACSM ligases. To identify variants/haplotypes that should be characterised in future detoxification association studies, this study assessed the naturally observed sequence diversity and protein expression variation of ACSM2A and ACSM2B. The allelic variation, haplotype diversity, Tajima's D values, and phylogenetic analyses indicated that ACSM2A and ACSM2B are highly conserved. This confirmed an earlier hypothesis that the glycine conjugation pathway is highly conserved and essential for life as it maintains the CoA and glycine homeostasis in the liver mitochondria. The protein expression analyses showed that ACSM2A is the predominant transcript in liver. Future studies should investigate the effect of the variants identified in this study on the substrate specificity of these proteins.

Ancestry and different rates of suicide and homicide in European countries: A study with population-level data

INTRODUCTION

There are large differences in suicide rates across Europe. The current study investigated the relationship of suicide and homicide rates in different countries of Europe with ancestry as it is defined with the haplotype frequencies of Y-DNA and mtDNA.

MATERIAL AND METHODS

The mortality data were retrieved from the WHO online database. The genetic data were retrieved from http://www.eupedia.com. The statistical analysis included Forward Stepwise Multiple Linear Regression analysis and Pearson Correlation Coefficient (R).

RESULTS

In males, N and R1a Y-DNA haplotypes were positively related to both homicidal and suicidal behaviors while I1 was negatively related. The Q was positively related to the homicidal rate. Overall, 60-75% of the observed variance was explained. L, J and X mtDNA haplogroups were negatively related with suicide in females alone, with 82-85% of the observed variance described.

DISCUSSION

The current study should not be considered as a study of genetic markers but rather a study of human ancestry. Its results could mean that research on suicidality has a strong biological but locally restricted component and could be limited by the study population; generalizability of the results at an international level might not be possible. Further research with patient-level data are needed to verify whether these haplotypes could serve as biological markers to identify persons at risk to commit suicide or homicide and whether biologically-determined ancestry could serve as an intermediate grouping method or even as an endophenotype in suicide research.

Ornithine Aminotransferase, an Important Glutamate-Metabolizing Enzyme at the Crossroads of Multiple Metabolic Pathways

Ornithine δ-aminotransferase (OAT, E.C. 2.6.1.13) catalyzes the transfer of the δ-amino group from ornithine (Orn) to α-ketoglutarate (aKG), yielding glutamate-5-semialdehyde and glutamate (Glu), and vice versa. In mammals, OAT is a mitochondrial enzyme, mainly located in the liver, intestine, brain, and kidney. In general, OAT serves to form glutamate from ornithine, with the notable exception of the intestine, where citrulline (Cit) or arginine (Arg) are end products. Its main function is to control the production of signaling molecules and mediators, such as Glu itself, Cit, GABA, and aliphatic polyamines. It is also involved in proline (Pro) synthesis. Deficiency in OAT causes gyrate atrophy, a rare but serious inherited disease, a further measure of the importance of this enzyme.

Gene-centric analysis implicates nuclear encoded mitochondrial protein gene variants in migraine susceptibility

Background

Migraine is a common neurological disorder which affects a large proportion of the population. The Norfolk Island population is a genetically isolated population and is an ideal discovery cohort for genetic variants involved in complex disease susceptibility given the reduced genetic and environmental heterogeneity. Given that the majority of proteins responsible for mitochondrial function are nuclear encoded, this study aimed to investigate the role of Nuclear Encoded Mitochondrial Protein (NEMP) genes in relation to migraine susceptibility.

Methods

A gene‐centric association analysis of NEMP genes was undertaken in the most related individuals (n = 315) within the genetically isolated Norfolk Island population. The discovery phase included genes with three or more SNP associations (P < 0.005), which were investigated further in a replication phase using an unrelated migraine case–control cohort (544 patients and 584 controls).

Results

The discovery phase of the study implicated SNPs in 5 NEMP genes to be associated with migraine susceptibility (P < 0.005). Replication analysis validated some of these implicated genes with SNPs in three NEMP genes shown to be associated with migraine in the replication cohort. These were CSNK1G3 (P = 0.00037), ELOVL6 (P = 0.00035) and SARDH (P = 0.00081), which are involved in phosphorylation, fatty acid metabolism, and oxidative demethylation, respectively.

Conclusion

Here we provide evidence that variation in NEMP genes is associated with migraine susceptibility. This study provides evidence for a link between mitochondrial function and migraine susceptibility.

Mitochondrial fatty acid synthesis, fatty acids and mitochondrial physiology

Mitochondria and fatty acids are tightly connected to a multiplicity of cellular processes that go far beyond mitochondrial fatty acid metabolism. In line with this view, there is hardly any common metabolic disorder that is not associated with disturbed mitochondrial lipid handling. Among other aspects of mitochondrial lipid metabolism, apparently all eukaryotes are capable of carrying out de novo fatty acid synthesis (FAS) in this cellular compartment in an acyl carrier protein (ACP)-dependent manner. The dual localization of FAS in eukaryotic cells raises the questions why eukaryotes have maintained the FAS in mitochondria in addition to the "classic" cytoplasmic FAS and what the products are that cannot be substituted by delivery of fatty acids of extramitochondrial origin. The current evidence indicates that mitochondrial FAS is essential for cellular respiration and mitochondrial biogenesis. Although both β-oxidation and FAS utilize thioester chemistry, CoA acts as acyl-group carrier in the breakdown pathway whereas ACP assumes this role in the synthetic direction. This arrangement metabolically separates these two pathways running towards opposite directions and prevents futile cycling. A role of this pathway in mitochondrial metabolic sensing has recently been proposed. This article is part of a Special Issue entitled: Lipids of Mitochondria edited by Guenther Daum.

Xenobiotic/medium chain fatty acid: CoA ligase - a critical review on its role in fatty acid metabolism and the detoxification of benzoic acid and aspirin

INTRODUCTION

Activation of fatty acids by the acyl-CoA synthetases (ACSs) is the vital first step in fatty acid metabolism. The enzymatic and physiological characterization of the human xenobiotic/medium chain fatty acid: CoA ligases (ACSMs) has been severely neglected even though xenobiotics, such as benzoate and salicylate, are detoxified through this pathway.

AREAS COVERED

This review will focus on the nomenclature and substrate specificity of the human ACSM ligases; the biochemical and enzymatic characterization of ACSM1 and ACSM2B; the high sequence homology of the ACSM2 genes (ACSM2A and ACSM2B) as well as what is currently known regarding disease association studies.

EXPERT OPINION

Several discrepancies exist in the current literature that should be taken note of. For example, the single nucleotide polymorphisms (SNPs) reported to be associated with aspirin metabolism and multiple risk factors of metabolic syndrome are incorrect. Kinetic data on the substrate specificity of the human ACSM ligases are non-existent and currently no data exist on the influence of SNPs on the enzyme activity of these ligases. One of the biggest obstacles currently in the field is that glycine conjugation is continuously studied as a one-step process, which means that key regulatory factors of the two individual steps remain unknown.

Molecular Characterization and Transcriptional Regulation Analysis of the Bovine PDHB Gene

The pyruvate dehydrogenase beta subunit (PDHB) is a subunit of pyruvate dehydrogenase (E1), which catalyzes pyruvate into acetyl-CoA and provides a linkage between the tricarboxylic acid cycle (TCA) and the glycolysis pathway. Previous studies demonstrated PDHB to be positively related to the intramuscular fat (IMF) content. However, the transcriptional regulation of PDHB remains unclear. In our present study, the cDNA of bovine PDHB was cloned and the genomic structure was analyzed. The phylogenetic tree showed bovine PDHB to be closely related to goat and sheep, and least related to chicken. Spatial expression pattern analysis revealed the products of bovine PDHB to be widely expressed with the highest level in the fat of testis. To understand the transcriptional regulation of bovine PDHB, 1899 base pairs (bp) of the 5’-regulatory region was cloned. Sequence analysis neither found consensus TATA-box nor CCAAT-box in the 5’-flanking region of bovine PDHB. However, a CpG island was predicted from nucleotides -284 to +117. Serial deletion constructs of the 5’-flanking region, evaluated in dual-luciferase reporter assay, revealed the core promoter to be located 490bp upstream from the transcription initiation site (+1). Electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation assay (ChIP) in combination with asite-directed mutation experiment indicated both myogenin (MYOG) and the CCAAT/enhancer-binding protein beta (C/EBPß) to be important transcription factors for bovine PDHB in skeletal muscle cells and adipocytes. Our results provide an important basis for further investigation of the bovine PDHB function and regulation in cattle.

HIV-1 transgenic rats display mitochondrial abnormalities consistent with abnormal energy generation and distribution

With the advent of the combination antiretroviral therapy era (cART), the development of AIDS has been largely limited in the USA. Unfortunately, despite the development of efficacious treatments, HIV-1-associated neurocognitive disorders (HAND) can still develop, and as many HIV-1 positive individuals age, the prevalence of HAND is likely to rise because HAND manifests in the brain with very low levels of virus. However, the mechanism producing this viral disorder is still debated. Interestingly, HIV-1 infection exposes neurons to proteins including Tat, Nef, and Vpr which can drastically alter mitochondrial properties. Mitochondrial dysfunction has been posited to be a cornerstone of the development of numerous neurodegenerative diseases. Therefore, we investigated mitochondria in an animal model of HAND. Using an HIV-1 transgenic rat model expressing seven of the nine HIV-1 viral proteins, mitochondrial functional and proteomic analysis were performed on a subset of mitochondria that are particularly sensitive to cellular changes, the neuronal synaptic mitochondria. Quantitative mass spectroscopic studies followed by statistical analysis revealed extensive proteome alteration in this model paralleling mitochondrial abnormalities identified in HIV-1 animal models and HIV-1-infected humans. Novel mitochondrial protein changes were discovered in the electron transport chain (ETC), the glycolytic pathways, mitochondrial trafficking proteins, and proteins involved in various energy pathways, and these findings correlated well with the function of the mitochondria as assessed by a mitochondrial coupling and flux assay. By targeting these proteins and proteins upstream in the same pathway, we may be able to limit the development of HAND.

Mitochondrial and apoptotic in vitro modelling of differential HIV-1 progression and antiretroviral toxicity

OBJECTIVES

Ex vivo analysis of mitochondrial function may reveal HIV progression and the impact of ART. We propose a mitochondrial and apoptotic in vitro model using Jurkat T cells incubated with plasma. The objectives of this study were to evaluate mitochondrial and apoptotic lesions in this model in relation to HIV progression, and to assess the effect of >1 year of standard non-thymidine-containing therapy.

METHODS

This was a cross-sectional comparison among three age- and gender-matched groups (n = 19 × 3): healthy non-HIV-infected participants, HIV-infected long-term non-progressors (LTNPs) and standard antiretroviral-naive chronically infected patients [standard progressors (Sps)], longitudinally evaluated before (Sp1) and after (Sp2) >1 year of efavirenz + tenofovir + emtricitabine therapy. We analysed mitochondrial DNA content by RT-PCR, mitochondrial function by spectrophotometry, mitochondrial protein synthesis by western blot analysis, mitochondrial dynamics by western blot analysis (MFN2), apoptotic transition pore formation by western blot analysis (VDAC-1) and mitochondrial membrane potential and annexin V/propidium iodide fluorescence by flow cytometry.

RESULTS

There was a decreasing non-significant trend towards lower mitochondrial parameters for HIV-infected values with respect to uninfected control reference values. HIV progression (LTNP versus Sp1) was associated with decreased mitochondrial genetic, functional and translational parameters, which partially recovered after treatment intervention (Sp2). Mitochondrial fusion showed a trend to decrease non-significantly in Sp patients compared with LTNP patients, especially after therapy. All apoptotic parameters showed a trend to increase in Sp1 with respect to LTNP, followed by recovery in Sp2.

CONCLUSIONS

We proposed an in vitro model for mitochondrial and apoptotic assessment to test the effects of HIV infection and its therapy, resembling in vivo conditions. This model could be useful for clinical research purposes.

GWATCH: a web platform for automated gene association discovery analysis

Background

As genome-wide sequence analyses for complex human disease determinants are expanding, it is increasingly necessary to develop strategies to promote discovery and validation of potential disease-gene associations.

Findings

Here we present a dynamic web-based platform – GWATCH – that automates and facilitates four steps in genetic epidemiological discovery: 1) Rapid gene association search and discovery analysis of large genome-wide datasets; 2) Expanded visual display of gene associations for genome-wide variants (SNPs, indels, CNVs), including Manhattan plots, 2D and 3D snapshots of any gene region, and a dynamic genome browser illustrating gene association chromosomal regions; 3) Real-time validation/replication of candidate or putative genes suggested from other sources, limiting Bonferroni genome-wide association study (GWAS) penalties; 4) Open data release and sharing by eliminating privacy constraints (The National Human Genome Research Institute (NHGRI) Institutional Review Board (IRB), informed consent, The Health Insurance Portability and Accountability Act (HIPAA) of 1996 etc.) on unabridged results, which allows for open access comparative and meta-analysis.

Conclusions

GWATCH is suitable for both GWAS and whole genome sequence association datasets. We illustrate the utility of GWATCH with three large genome-wide association studies for HIV-AIDS resistance genes screened in large multicenter cohorts; however, association datasets from any study can be uploaded and analyzed by GWATCH.

Mitochondrial DNA variation and virologic and immunological HIV outcomes in African Americans

OBJECTIVE

To evaluate the impact of mitochondrial DNA (mtDNA) haplogroups on virologic and immunological outcomes of HIV infection.

DESIGN

HAART-naive African American adolescent participants to the Reaching for Excellence in Adolescent Care and Health study.

METHODS

The mtDNA haplogroups were inferred from sequenced mtDNA hypervariable regions HV1 and HV2 and their predictive value on HIV outcomes were evaluated in linear mixed models, controlled for human leukocyte antigen (HLA)-B27, HLA-B57 and HLA-B35-Px alleles and other covariates.

RESULTS

We report data showing that the mtDNA L2 lineage, a group composed of L2a, L2b and L2e mtDNA haplogroups in the studied population, is significantly associated (beta  = -0.08; Bonferroni-adjusted P = 0.004) with decline of CD4 T cells (median loss of 8 ± 1 cells per month) in HAART-naive HIV-infected individuals of African American descent (n = 133). No significant association (P < 0.05) with set-point viral load was observed with any of the tested mtDNA haplogroups. The present data concur with previous findings in the AIDS Clinical Trials Group study 384, implicating the L2 lineage with slower CD4 T-cell recovery after antiretroviral therapy in African Americans.

CONCLUSIONS

Whereas the L2 lineage showed an association with unfavorable immunological outcomes of HIV infection, its phylogenetic divergence from J and U5a, two lineages associated with accelerated HIV progression in European Americans, raises the possibility that interactions with common nucleus-encoded variants drive HIV progression. Disentangling the effects of mitochondrial and nuclear gene variants on the outcomes of HIV infection is an important step to be taken toward a better understanding of HIV/AIDS pathogenesis and pharmacogenomics.

ACSM4 polymorphisms are associated with rapid AIDS progression in HIV-infected patients

: Our aim was to explore the association among ACSM4 and PECI polymorphisms and AIDS progression in 454 HIV-infected patients never treated with antiretroviral drugs (146 long-term nonprogressors, 228 moderate progressors, and 80 rapid progressors). For ACSM4 polymorphisms, rs7137120 AA/AG and rs7961991 CC/CT genotypes had higher odds of having a rapid AIDS progression [odds ratio (OR) = 3.21; 95% of confidence interval (95% CI) = 1.26 to 8.16; P = 0.014 and OR = 3.60; 95% CI = 1.38 to 9.36; P = 0.009, respectively]. Additionally, the ACSM4 haplotype integrated for both rs7961991 A and rs7137120 C alleles had higher odds of having a rapid AIDS progression (OR = 2.85; 95% CI = 1.28 to 6.25; P = 0.010). For PECI polymorphisms, no significant associations were found. In conclusion, ACSM4 polymorphisms might play a significant role in AIDS progression.

Evaluation of the pharmacogenetics of immune recovery in treated HIV-infected patients

INTRODUCTION

Combination antiretroviral therapy has markedly improved the survival rate and quality of life in patients infected with HIV due to the powerful suppressor effect that current antiretroviral drugs have on the viral load. Consequently, the immune system undergoes a substantial qualitative and quantitative improvement; and this leads to an increase in the absolute CD4(+) T-lymphocyte count and the restoration of lost T-cell responses against certain opportunistic pathogens. Unfortunately, not all patients who successfully suppress plasma viremia experience sufficient CD4(+) T-cell gain and these patients, in turn, are associated with worse outcomes. Pharmacogenetic studies have been used to investigate how a patient's genetic predisposition may affect their response to antiretroviral drugs.

AREAS COVERED

This article reviews the investigations that have been published on the association between host genetic determinants of CD4(+) T-cell gain in treated HIV-infected patients. Studies were identified through a PubMed database search. Longitudinal studies into pharmacogenetic association were specifically selected.

EXPERT OPINION

While the possibility of genetic predisposition to HIV therapeutics has potential, most studies provide inconsistent data. Inconsistency is often due to partial genetic evaluation, different categorization of poor immune recovery or due to small numbers of patients evaluated. Currently, studies still belong to the research laboratory stage and more studies are required to improve our understanding.

HIV-1 infection is blocked at an early stage in cells devoid of mitochondrial DNA

Human immunodeficiency virus type I (HIV-1) exploits various host cellular pathways for efficient infection. Here we report that the absence of mitochondrial DNA (mtDNA) in ρ(0) cells markedly attenuates HIV-1 infection. Importantly, reduced infection efficiency in ρ(0) cells is not simply the result of impaired oxidative phosphorylation (OXPHOS) because pharmacological OXPHOS inhibition did not inhibit HIV-1 infection. Analysis of the early steps of virus infection by real-time PCR quantification of stage-specific HIV-1 DNA products in the infected ρ(0) and parental cell line have allowed us to conclude that HIV-1 infection in ρ(0) cells is blocked at the steps that occur after reverse transcription and prior to nuclear import. Additionally, confocal fluorescence microscope analysis showed that the majority of viral complexes containing HIV-1 p24 co-localize with mitochondria in target cells, suggesting an interaction between the two. Collectively, our data strongly indicate that mitochondria play an important role during early stages of HIV-1 infection, probably through direct association with HIV-1 intracellular complexes.

Genome-wide and candidate gene association studies of placental abruption

Placental abruption (PA), a pregnancy-related vascular disorder, is a leading cause of maternal and perinatal morbidity and mortality. The success of identifying genetic susceptibility loci for PA, a multi-factorial heritable disorder, has been limited. We conducted a genome-wide association study (GWAS) and candidate gene association study using 470 PA cases and 473 controls from Lima, Peru. Genotyping for common genetic variations (single nucleotide polymorphisms, SNPs) was conducted using the Illumina Cardio-Metabo Chip platform. Common variations in 35 genes that participate in mitochondrial biogenesis (MB) and oxidative phosphorylation (OS) were selected for the candidate gene study. Regression models were fit to examine associations of each SNP with risk of PA. In pathway analyses, we examined functions and functional relationships of genes represented by the top GWAS hits. Genetic risk scores (GRS), based on top hits of the GWAS and candidate gene analyses, respectively, were computed using the risk allele counting method. The top hit in the GWAS analyses was rs1238566 (empirical P-value=1.04e-4 and FDR-adjusted P-value=5.65E-04) in FLI-1 gene, a megakaryocyte-specific transcription factor. Networks of genes involved in lipid metabolism and cell signaling were significantly enriched by the 51 genes whose SNPs were among the top 200 GWAS hits (P-value <2.1e-3). SNPs known to regulate MB (e.g. CAMK2B, NR1H3, PPARG, PRKCA, and THRB) and OP (e.g., COX5A, and NDUF family of genes) were associated with PA risk (P-value <0.05). GRS was significantly associated with PA risk (trend P-value <0.001 and 0.01 for GWAS and candidate gene based GRS, respectively). Our study suggests that integrating multiple analytical strategies in genetic association studies can provide opportunities for identifying genetic risk factors and novel molecular mechanisms that underlie PA.

IL7RA haplotype-associated alterations in cellular immune function and gene expression patterns in multiple sclerosis

Interleukin-7 receptor alpha (IL7RA) is among the top listed candidate genes influencing the risk to develop multiple sclerosis (MS), an inflammatory demyelinating disease of the central nervous system. Soluble IL-7RA (sIL-7RA) protein and mRNA levels vary among the four common IL7RA haplotypes. Here we show and confirm that protective haplotype carriers have three times lower sIL-7RA serum levels than the other three haplotypes. High sIL-7RA concentrations significantly decrease IL-7-mediated STAT5 phosphorylation in CD4(+) T cells. Transcriptome analysis of unstimulated and stimulated CD4(+) T cells of MS patients carrying the different IL7RA haplotypes revealed complex and overlapping patterns in genes participating in cytokine signaling networks, apoptosis, cell cycle progression and cell differentiation. Our findings indicate that genetic variants of IL7RA result in haplotype-associated differential responsiveness to immunological stimuli that influence MS susceptibility not exclusively by varying levels of sIL-7RA.

Bridging HIV-1 cellular latency and clinical long-term non-progressor: an interactomic view

Development of an effective HIV management is enticed by the fact that long-term non-progressors (LTNP) restrict viral replication spontaneously, but is hindered by HIV-1 latency. Given that the most overlapping characteristics found between HIV-1 LTNP and latency, detailed analysis of the difference would disclose the essentials of latency. In this study, microarray data from our previous study was combined with HIV-1 latency and LTNP data obtained from NCBI GEO database. Principal variance component analysis and hierarchical clustering verified the removal of batch effect across platform. The analysis revealed a total of 456 differential expressed genes with >2-fold change and B-statistic >0. Bayesian inference was used to reconstitute the transcriptional network of HIV-1 latency or LTNP, respectively. Gene regulation was reprogrammed under different disease condition. By network interference, KPNA2 and ATP5G3 were identified as the hubs in latency network which mediate nuclear export and RNA processing. These data offer comparative insights into HIV-1 latency, which will facilitate the understanding of the genetic basis of HIV-1 latency in vivo and serve as a clue for future treatment dealing with key targets in HIV-1 latency.

Host genomic influences on HIV/AIDS

The AIDS era has seen multiple advances in the power of genetics research; scores of host genetic protective factors have been nominated and several have translated to the bedside. We discuss how genomics may inform HIV/AIDS prevention, treatment and eradication.

Analysis of differentially expressed genes and microRNAs in alcoholic liver disease

The purpose of this study was to screen differentially expressed genes and microRNAs in order to find a new target for the accurate diagnosis and effective therapy of alcoholic liver disease (ALD) at the gene and microRNA levels. The total RNA of liver tissues was extracted from four groups of patients, ten subjects each. Microarrays were utilized to detect differentially expressed genes and microRNAs. According to gene values, significance levels and false discovery rate with a random variance model, gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, node genes and key microRNAs in networks were obtained and analyzed. A total of 878 differentially expressed genes and 26 microRNAs were found. In co-expression genetic networks, node genes modulating the network were Acyl-coenzyme A synthetase-3 (ACSF3), Frizzled-5 (FZD5), LOC727987 and C1orf222. In microRNA-gene networks, the key microRNAs were hsa-miR-570, hsa-miR-122, hsa-miR-34b, hsa-miR-29c, hsa-miR-922 and hsa-miR-185, which negatively regulated approximately 79 downstream target genes. In the course of ALD, we found 4 differentially expressed node genes and analyzed ACSF3 and FZD5. ACSF3 was significantly upregulated, and was involved in fatty acid and lipid metabolism and accelerated liver injury. These two genes were involved in fatty acids and lipid metabolism. FZD5 was downregulated and reduced the synthesis of membrane transport protein in the hepatic membrane and the membrane stability, and accelerated the liver cell apoptosis process. Six key microRNAs regulated numerous biological functions such as the immune response, the inflammatory response and glutathione metabolism. This finding provides valuable insight into the diagnosis and treatment of ALD.

Generation of a genome scale lentiviral vector library for EF1α promoter-driven expression of human ORFs and identification of human genes affecting viral titer

The bottleneck in elucidating gene function through high-throughput gain-of-function genome screening is the limited availability of comprehensive libraries for gene overexpression. Lentiviral vectors are the most versatile and widely used vehicles for gene expression in mammalian cells. Lentiviral supernatant libraries for genome screening are commonly generated in the HEK293T cell line, yet very little is known about the effect of introduced sequences on the produced viral titer, which we have shown to be gene dependent. We have generated an arrayed lentiviral vector library for the expression of 17,030 human proteins by using the GATEWAY® cloning system to transfer ORFs from the Mammalian Gene Collection into an EF1alpha promoter-dependent lentiviral expression vector. This promoter was chosen instead of the more potent and widely used CMV promoter, because it is less prone to silencing and provides more stable long term expression. The arrayed lentiviral clones were used to generate viral supernatant by packaging in the HEK293T cell line. The efficiency of transfection and virus production was estimated by measuring the fluorescence of IRES driven GFP, co-expressed with the ORFs. More than 90% of cloned ORFs produced sufficient virus for downstream screening applications. We identified genes which consistently produced very high or very low viral titer. Supernatants from select clones that were either high or low virus producers were tested on a range of cell lines. Some of the low virus producers, including two previously uncharacterized proteins were cytotoxic to HEK293T cells. The library we have constructed presents a powerful resource for high-throughput gain-of-function screening of the human genome and drug-target discovery. Identification of human genes that affect lentivirus production may lead to improved technology for gene expression using lentiviral vectors.

Identification of critical host mitochondrion-associated genes during Ehrlichia chaffeensis infections

Ehrlichia chaffeensis is an obligate intracellular bacterium that causes human monocytic ehrlichiosis (HME). To determine what host components are important for bacterial replication, we performed microarray analysis on Drosophila melanogaster S2 cells by comparing host gene transcript levels between permissive and nonpermissive conditions for E. chaffeensis growth. Five-hundred twenty-seven genes had increased transcript levels unique to permissive growth conditions 24 h postinfection. We screened adult flies that were mutants for several of the "permissive" genes for the ability to support Ehrlichia replication. Three additional D. melanogaster fly lines with putative mutations in pyrimidine metabolism were also tested. Ten fly lines carrying mutations in the genes CG6479, separation anxiety, chitinase 11, CG6364 (Uck2), CG6543 (Echs1), withered (whd), CG15881 (Ccdc58), CG14806 (Apop1), CG11875 (Nup37), and dumpy (dp) had increased resistance to infection with Ehrlichia. Analysis of RNA by quantitative real-time reverse transcription-PCR (qRT-PCR) confirmed that the bacterial load was decreased in these mutant flies compared to wild-type infected control flies. Seven of these genes (san, Cht11, Uck2, Echs1, whd, Ccdc58, and Apop1) encoded proteins that had mitochondrial functions or could be associated with proteins with mitochondrial functions. Treatment of THP-1 cells with double-stranded RNA to silence the human UCK2 gene indicates that the disruption of the uridine-cytidine kinase affects E. chaffeensis replication in human macrophages. Experiments with cyclopentenyl cytosine (CPEC), a CTP synthetase inhibitor and cytosine, suggest that the nucleotide salvage pathway is essential for E. chaffeensis replication and that it may be important for the provision of CTP, uridine, and cytidine nucleotides.

Association of genetic copy number variations at 11 q14.2 with brain regional volume differences in an alcohol use disorder population

This study investigates the relationship between genetic copy number variations and brain volume differences in an alcohol use disorder (AUD) population. We hypothesized that copy number variations may influence subject's risk for alcohol use disorders through variations in regional gray and white matter brain volumes. Since genetic influences upon behavior are the result of many complicated interactions we focus on differences in brain volume as a putative intermediate phenotype between genetic variation and behavior. Copy number variation, alcohol use assessments and brain structural magnetic resonance images from 283 subjects, 199 male and 84 females who were enrolled in two AUD studies were obtained and analyzed using a combination of the Freesurfer image analysis suite and independent component analysis. Because brain volume varies by age we compared participant's volume variation with that derived from a control cohort of 75 subjects. In addition we also regressed out the possible brain volume changes induced by long term alcohol consumption. Small cerebral cortex, cerebellar and caudate along with large cerebral white matter and 5th ventricle volumes are shown to be significantly associated with increased AUD severity. When these volume variations are compared with control subject volumes; the variations seen in subjects with AUD are markedly different from normal aging effects. CNVs at 11 q14.2 are marginally (p < 0.05 uncorrected) correlated with such brain volume variations and the correlation holds true after controlling for long-term alcohol consumption; deletion carriers have smaller cerebral cortex, cerebellar, caudate and larger cerebral white matter and 5th ventricle volumes than insertion carriers or subjects with no variation in this region. Similarly, deletion carriers also demonstrate higher AUD severity scores than insertion carriers or subjects with no variation. The results presented here suggest that copy number variation and in particular the variation at chromosome 11 q14.2 may have an impact in brain volume variation, potentially influencing AUD behavior.

Mitochondrial proteomic analysis of cisplatin resistance in ovarian cancer

Epithelial ovarian cancer (EOC) is the leading cause of death among women with gynecologic malignancies and accounts for approximately 6% of cancer deaths among women. Cisplatin and its analogues form the backbone of the most active chemotherapy regimens in advanced EOC; however, development of platinum resistance is common and typically marks a transition in which curing the patient is no longer possible. An emerging theme in many cancers is that mitochondrial dysfunction contributes to an aggressive carcinogenic phenotype. We hypothesized that changes in the mitochondrial proteome are required to support development of cisplatin resistance in human EOC. To investigate this hypothesis, an organellar proteomics approach was utilized to quantify alterations in protein abundance in mitochondria enriched from isogenic cisplatin-sensitive (A2780) and -resistant (A2780-CP20) human EOC cells. Protein isolates from mitochondria-enriched fractions were analyzed by high resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS), and relative abundance of identified proteins was quantified by spectral counting. Pathway analyses revealed significant increases in notch signaling pathways, cell survival, and alternate apoptotic pathways in the A2780-CP20 subtype. Among the alterations identified in the mitochondrial proteomic composition in cisplatin-resistant EOC cells, activated leukocyte cell adhesion molecule (AKAP12) and A kinase anchoring protein 12 (AKAP12) were elevated, while nestin was diminished in the mitochondrial fraction of A2780-CP20 relative to A2780. This was verified by immunoblot analysis. These results confirm that important changes in the mitochondrial proteome, many of which promote evasion of apoptosis and tumor invasiveness and metastasis, are present in cisplatin-resistant EOC.

Family-based association studies for next-generation sequencing

An individual's disease risk is determined by the compounded action of both common variants, inherited from remote ancestors, that segregated within the population and rare variants, inherited from recent ancestors, that segregated mainly within pedigrees. Next-generation sequencing (NGS) technologies generate high-dimensional data that allow a nearly complete evaluation of genetic variation. Despite their promise, NGS technologies also suffer from remarkable limitations: high error rates, enrichment of rare variants, and a large proportion of missing values, as well as the fact that most current analytical methods are designed for population-based association studies. To meet the analytical challenges raised by NGS, we propose a general framework for sequence-based association studies that can use various types of family and unrelated-individual data sampled from any population structure and a universal procedure that can transform any population-based association test statistic for use in family-based association tests. We develop family-based functional principal-component analysis (FPCA) with or without smoothing, a generalized T(2), combined multivariate and collapsing (CMC) method, and single-marker association test statistics. Through intensive simulations, we demonstrate that the family-based smoothed FPCA (SFPCA) has the correct type I error rates and much more power to detect association of (1) common variants, (2) rare variants, (3) both common and rare variants, and (4) variants with opposite directions of effect from other population-based or family-based association analysis methods. The proposed statistics are applied to two data sets with pedigree structures. The results show that the smoothed FPCA has a much smaller p value than other statistics.

An automated bayesian framework for integrative gene expression analysis and predictive medicine

MOTIVATION

This work constructs a closed loop Bayesian Network framework for predictive medicine via integrative analysis of publicly available gene expression findings pertaining to various diseases.

RESULTS

An automated pipeline was successfully constructed. Integrative models were made based on gene expression data obtained from GEO experiments relating to four different diseases using Bayesian statistical methods. Many of these models demonstrated a high level of accuracy and predictive ability. The approach described in this paper can be applied to any complex disorder and can include any number and type of genome-scale studies.

Alterations in metabolism-related genes induced in SHSY5Y cells by okadaic acid exposure

Okadaic acid (OA) is a widely distributed marine toxin produced by several phytoplanktonic species and responsible for diarrheic shellfish poisoning in humans. At the molecular level OA is a specific inhibitor of several types of serine/threonine protein phosphatases. Due to this enzymic inhibition, OA was reported to induce numerous alterations in relevant cellular physiological processes, including several metabolic pathways such as glucose uptake, lipolysis and glycolysis, heme metabolism, and glycogen and protein synthesis. In order to further understand the underlying mechanisms involved in OA-induced effects on cellular metabolism, the expression levels of six genes related to different catabolic and anabolic metabolism-related processes were analyzed by real-time polymerase chain reaction. Specifically, the expression patterns of GAPDH, TOMM5, SLC25A4, COII, QARS, and RGS5 genes were determined in SHSY5Y human neuroblastoma cells exposed to OA for 3, 24, or 48 h. All these genes showed alterations in their expression levels after at least one of the OA treatments tested. These alterations provide a basis to understand the mechanisms underlying the previously described OA-induced effects on different metabolic processes, mainly regarding glucose and mitochondrial metabolism. However, other OA-induced affected genes can not be ruled out, and further studies are required to more comprehensively characterize in the mechanisms of OA-induced interaction on cell metabolism.