Analysis of TAP2 polymorphisms in Finnish individuals with type I diabetes

Identification of biomarkers, pathways, and potential therapeutic targets for heart failure using next-generation sequencing data and bioinformatics analysis

BACKGROUND

Heart failure (HF) is the most common cardiovascular diseases and the leading cause of cardiovascular diseases related deaths. Increasing molecular targets have been discovered for HF prognosis and therapy. However, there is still an urgent need to identify novel biomarkers. Therefore, we evaluated biomarkers that might aid the diagnosis and treatment of HF.

METHODS

We searched next-generation sequencing (NGS) dataset (GSE161472) and identified differentially expressed genes (DEGs) by comparing 47 HF samples and 37 normal control samples using limma in R package. Gene ontology (GO) and pathway enrichment analyses of the DEGs were performed using the g: Profiler database. The protein-protein interaction (PPI) network was plotted with Human Integrated Protein-Protein Interaction rEference (HiPPIE) and visualized using Cytoscape. Module analysis of the PPI network was done using PEWCC1. Then, miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed by Cytoscape software. Finally, we performed receiver operating characteristic (ROC) curve analysis to predict the diagnostic effectiveness of the hub genes.

RESULTS

A total of 930 DEGs, 464 upregulated genes and 466 downregulated genes, were identified in HF. GO and REACTOME pathway enrichment results showed that DEGs mainly enriched in localization, small molecule metabolic process, SARS-CoV infections, and the citric acid tricarboxylic acid (TCA) cycle and respiratory electron transport. After combining the results of the PPI network miRNA-hub gene regulatory network and TF-hub gene regulatory network, 10 hub genes were selected, including heat shock protein 90 alpha family class A member 1 (HSP90AA1), arrestin beta 2 (ARRB2), myosin heavy chain 9 (MYH9), heat shock protein 90 alpha family class B member 1 (HSP90AB1), filamin A (FLNA), epidermal growth factor receptor (EGFR), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), cullin 4A (CUL4A), YEATS domain containing 4 (YEATS4), and lysine acetyltransferase 2B (KAT2B).

CONCLUSIONS

This discovery-driven study might be useful to provide a novel insight into the diagnosis and treatment of HF. However, more experiments are needed in the future to investigate the functional roles of these genes in HF.

Identification and Characterization of Renal Cell Carcinoma Gene Markers

Microarray gene expression profiling has been used to distinguish histological subtypes of renal cell carcinoma (RCC), and consequently to identify specific tumor markers. The analytical procedures currently in use find sets of genes whose average differential expression across the two categories differ significantly. In general each of the markers thus identified does not distinguish tumor from normal with 100% accuracy, although the group as a whole might be able to do so. For the purpose of developing a widely used economically viable diagnostic signature, however, large groups of genes are not likely to be useful. Here we use two different methods, one a support vector machine variant, and the other an exhaustive search, to reanalyze data previously generated in our Lab (Lenburg et al. 2003). We identify 158 genes, each having an expression level that is higher (lower) in every tumor sample than in any normal sample, and each having a minimum differential expression across the two categories at a significance of 0.01. The set is highly enriched in cancer related genes (p = 1.6 × 10–12), containing 43 genes previously associated with either RCC or other types of cancer. Many of the biomarkers appear to be associated with the central alterations known to be required for cancer transformation. These include the oncogenes JAZF1, AXL, ABL2; tumor suppressors RASD1, PTPRO, TFAP2A, CDKN1C; and genes involved in proteolysis or cell-adhesion such as WASF2, and PAPPA.

A genetic risk score for thyroid peroxidase antibodies associates with clinical thyroid disease in community-based populations

CONTEXT

Antibodies against thyroid peroxidase (TPOAbs) are detected in 90% of all patients with Hashimoto thyroiditis, the most common cause of hypothyroidism. Hypothyroidism is associated with a range of adverse outcomes. The current knowledge of its genetic underpinnings is limited.

OBJECTIVE

The purpose of this study was to identify novel genetic variants associated with TPOAb concentrations and positivity using genome-wide association data and to characterize their association with thyroid function and disease.

DESIGN, SETTING, AND PARTICIPANTS

We studied European ancestry participants of 3 independent prospective population-based studies: Atherosclerosis Risk In Communities study (n = 7524), Study of Health in Pomerania (n = 3803), and Study of Health in Pomerania-TREND (n = 887).

EXPOSURE

Single nucleotide polymorphisms (SNPs), individually and combined into a genetic risk score (GRS), were examined.

MAIN OUTCOMES

The main outcomes were TPOAb concentrations and positivity, thyroid hormone concentrations (TSH, free T4), and clinical thyroid diseases (subclinical and overt hypothyroidism and goiter).

RESULTS

Significantly associated single nucleotide polymorphisms (P < 5 · 10(-8)) mapped into 4 genomic regions not previously implicated for TPOAbs (RERE, extended HLA region) and into 5 previously described loci. A higher Genetic Risk Score (GRS) based on these 9 SNPs showed strong and graded associations with higher TPOAb, TSH, and lower free T4 concentrations (P < .001). Compared with individuals in the lowest GRS quartile, those in the highest quartile had 1.80-fold higher odds of subclinical hypothyroidism (95% confidence interval, 1.27-2.55) and 1.89-fold higher odds of overt hypothyroidism (95% confidence interval, 1.24-2.87).

CONCLUSION

The identification of 4 novel genetic loci associated with TPOAb concentrations and positivity gives further insight into the genetic underpinnings of hypothyroidism. A GRS showed strong and graded associations with markers of thyroid function and disease in independent population-based studies.

Association of polymorphisms in HLA antigen presentation-related genes with the outcomes of HCV infection

Antigen-presentation genes play a vital role in the pathogenesis of HCV infection. However, the relationship of variants of these genes with spontaneous outcomes of HCV infection has not been fully investigated. To explore novel loci in the Chinese population, 34 tagging-SNPs in 9 candidate genes were genotyped for their associations with the outcomes of HCV infection. The distributions of different genotypes and haplotypes were compared among 773 HCV-negative controls, 246 subjects with HCV natural clearance, and 218 HCV persistent carriers recruited from hemodialysis patients and intravenous drug users. Our study implicated that TAP2, HLA-DOA, HLA-DOB, and tapasin loci were novel candidate regions for susceptibility to HCV infection and viral clearance in the Chinese population. Logistic regression analyses showed that TAP2 rs1800454 A (OR = 1.48, P = 0.002) and HLA-DOB rs2071469 G (OR = 1.23, P = 0.048) were significantly associated with increased susceptibility to establishment of HCV infection. However, high-risk behavior exposure and age were stronger predictors of HCV infection. Mutation of tapasin rs9277972 T (OR = 1.57, P =0.043) increased the risk of HCV chronicity, and HLA-DOA rs3128935 C (OR = 0.62, P = 0.019) increased the chance of viral resolution. With regards to the effect of rs3128925, interactions were found with high-risk behavior (P = 0.013) and age (P = 0.035). The risk effect of rs3128925 T for persistent HCV infection was higher in injecting drug users (vs. dialysis patients) and in subjects ≥ 40 years old (vs. < 40 years old).

Innate immune activity is detected prior to seroconversion in children with HLA-conferred type 1 diabetes susceptibility

The insult leading to autoantibody development in children who will progress to develop type 1 diabetes (T1D) has remained elusive. To investigate the genes and molecular pathways in the pathogenesis of this disease, we performed genome-wide transcriptomics analysis on a unique series of prospective whole-blood RNA samples from at-risk children collected in the Finnish Type 1 Diabetes Prediction and Prevention study. We studied 28 autoantibody-positive children, out of which 22 progressed to clinical disease. Collectively, the samples covered the time span from before the development of autoantibodies (seroconversion) through the diagnosis of diabetes. Healthy control subjects matched for date and place of birth, sex, and HLA-DQB1 susceptibility were selected for each case. Additionally, we genotyped the study subjects with Immunochip to identify potential genetic variants associated with the observed transcriptional signatures. Genes and pathways related to innate immunity functions, such as the type 1 interferon (IFN) response, were active, and IFN response factors were identified as central mediators of the IFN-related transcriptional changes. Importantly, this signature was detected already before the T1D-associated autoantibodies were detected. Together, these data provide a unique resource for new hypotheses explaining T1D biology.

Autoimmune disease classification by inverse association with SNP alleles

With multiple genome-wide association studies (GWAS) performed across autoimmune diseases, there is a great opportunity to study the homogeneity of genetic architectures across autoimmune disease. Previous approaches have been limited in the scope of their analysis and have failed to properly incorporate the direction of allele-specific disease associations for SNPs. In this work, we refine the notion of a genetic variation profile for a given disease to capture strength of association with multiple SNPs in an allele-specific fashion. We apply this method to compare genetic variation profiles of six autoimmune diseases: multiple sclerosis (MS), ankylosing spondylitis (AS), autoimmune thyroid disease (ATD), rheumatoid arthritis (RA), Crohn's disease (CD), and type 1 diabetes (T1D), as well as five non-autoimmune diseases. We quantify pair-wise relationships between these diseases and find two broad clusters of autoimmune disease where SNPs that make an individual susceptible to one class of autoimmune disease also protect from diseases in the other autoimmune class. We find that RA and AS form one such class, and MS and ATD another. We identify specific SNPs and genes with opposite risk profiles for these two classes. We furthermore explore individual SNPs that play an important role in defining similarities and differences between disease pairs. We present a novel, systematic, cross-platform approach to identify allele-specific relationships between disease pairs based on genetic variation as well as the individual SNPs which drive the relationships. While recognizing similarities between diseases might lead to identifying novel treatment options, detecting differences between diseases previously thought to be similar may point to key novel disease-specific genes and pathways.

Genetic control of alternative splicing in the TAP2 gene: possible implication in the genetics of type 1 diabetes

The transporter 2, ATP-binding cassette, subfamily B (TAP2) is involved in the transport of antigenic peptides to HLA molecules. Coding TAP2 polymorphisms shows a strong association with type 1 diabetes, but it is not clear whether this association may be entirely due to linkage disequilibrium with HLA DR and DQ. Functionally, rat Tap2 nonsynonymous single-nucleotide polymorphisms (nsSNPs) confer differential selectivity for antigenic peptides, but this was not shown to be the case for human TAP2 nsSNPs. In the human, differential peptide selectivity is rather conferred by two splicing isoforms with alternative carboxy terminals. Here, we tested the hypothesis that alleles at the coding SNPs favor different splicing isoforms, thus determining peptide selectivity indirectly. This may be the basis for independent contribution to the type 1 diabetes association. In RNA from heterozygous lymphoblastoid lines, we measured the relative abundance of each SNP haplotype in each isoform. In isoform NM_000544, the G (Ala) allele at 665 Thr>Ala (rs241447) is more than twice as abundant as A (Thr) (GA = 2.2 +/- 0.4, P = 1.5 x 10(-4)), while isoform NM_018833 is derived almost exclusively from chromosomes carrying A (AG = 18.1 +/- 5.6, P = 2.04 x 10(-7)). In 889 Canadian children with type 1 diabetes, differential transmission of parental TAP2 alleles persisted (P = 0.011) when analysis was confined to chromosomes carrying only DQ*02 alleles, which mark a conserved DR-DQ haplotype, thus eliminating most of the variation at DR-DQ. Thus, we present evidence of TAP2 association with type 1 diabetes that is independent of HLA DR-DQ and describe a plausible functional mechanism based on allele dependence of splicing into isoforms known to have differential peptide selectivities.

Genetic susceptibility to type 1 diabetes in the intracellular pathway of antigen processing - a subject review and cross-study comparison

Ligand binding grooves of MHC class I molecules are able to load a panel of endogenous peptides of varying length and sequence derived from self or foreign origin to activate or deactivate cytotoxic CD8(+) T cells. Peptides are assembled with class I molecules by pathways that are either dependent or independent of transport by ABC proteins (TAP) and degradation in the immunoproteasome by its subunits LMP2 and LMP7. Those peptides that require TAP and LMP treatment appear to be subject to control and optimization by TAP for proper customizing and efficient presentation. Therefore, allelic variations in the coding sequences of TAP and LMP were suspected for a long time to be responsible for improper antigen processing, interruption of self-peptide presentation and reduced cell surface expression of MHC class I molecules resulting in the activation of autoreactive CD8(+) T cells. In this article we reviewed the controversial findings regarding the role of TAP and LMP genes in autoimmune diabetes and reevaluated data of eleven separate studies in a cross-study analysis by genotype and HLA haplotype matching. We could confirm previous results by showing that TAP2*651-A/F and TAP2*687-A/A are significantly associated with disease, independently of linkage disequilibrium (LD). LMP2-R/H surprisingly seems to be primarily disease-conferring although a weak association with DR4 serotypes can be observed. Our analysis also suggests that LMP7-B/B, TAP1-A/A and TAP2*687-A/B are the protective genotypes and that these associations are not secondary to LD with DRB1. Consequently, intracellular antigen processing associated with TAP- and proteasome-dependent pathways seems to be a critical element in T cell selection for the retention of a balanced immunity.

New transporter associated with antigen processing (TAP-2) polymorphisms in the Shona people of Zimbabwe

Most studies, to date, on transporter associated with antigen processing (TAP2) polymorphism have been conducted in Caucasians or Asians from industrialized countries. Because of the essential role of this molecule in antigen processing, the implication that polymorphism could be a major factor in human disease and the possible genetic variation at this locus among ethnically diverse populations, we undertook a study to analyze the full extent of TAP2 polymorphism in an indigenous Zimbabwean population (Shona ethnic group). Using single-stranded conformation polymorphism and DNA direct sequencing procedures, we detected the presence of 17 nucleotide sequence variations in the entire coding region of TAP2. Of these variants, 11 are nonconservative substitutions with respect to amino acid composition and are located in a region of the protein that could modulate its function. Six new polymorphic sites were identified in exon 1 (codons 15 Val-->Ala, 53 Leu-->Val), exon 3 (codon 220 Arg-->Arg), exon 4 (codons 257 Thr-->Ile, 313 Arg-->His), and exon10 (codon 609 Ala-->Val). Significant differences were seen in the distribution of the known 374Thr, 565Thr and 651Cys variants between African and non-African populations. These differences may reflect evolutionary pressures generated by environmental factors, such as prevalent pathogens in these geographically distinct regions. Further studies are needed to elucidate the net impact of TAP2 polymorphism on the protein's function and it's role in disease pathogenesis.

Polymorphisms of the TAP1 and TAP2 genes in human alveolar echinococcosis

We postulated that TAP genes may influence the susceptibility of some individuals to Echinococcus multilocularis infection. Six coding region variants (codons 333 and 637 in TAP1, and 379, 565, 651 and 665 in TAP2) were typed in 94 patients and 100 controls. Thr/Thr homozygosity at TAP2/665 was more prevalent in patients than in controls [64% vs. 45%, respectively; odds ratio (OR) = 2.1 (95% confidence interval (CI) 1.1; 2.7)] and Thr/Ala heterozygozity was less prevalent (32% vs. 50%, respectively) (P = 0.014). Of the 38 patients with progressive lesions, 76% were Thr/Thr, as compared with 55% of patients without progressive lesions and 45% of controls (P = 0.058 and 0.02, respectively), independent of HLA status. To determine whether this association is functionally relevant, functional analyses and/or confirmation in distinct populations of patients with alveolar echinococcosis would be required.

Polymorphisms of human TAP2 detected by denaturing gradient gel electrophoresis

The human transporter associated with antigen processing (TAP1 and TAP2) genes are located in the human leukocyte antigen (HLA) class II region of the genome and encode proteins that form a heterodimer essential for the transport of endogenous peptides into the endoplasmic reticulum for assembly with HLA class I molecules. Type 1 diabetes is an autoimmune disease that is associated with the HLA region of the genome, with HLA class II genes conferring the greatest statistical risk. The presentation of self-peptides by HLA class I molecules is defective in individuals with this disease, and both TAP1 and TAP2 are potential contributors to this defect. Denaturing gradient gel electrophoresis (DGGE) was applied to screen all 11 exons and the 3' flanking region of TAP2 for polymorphisms in individuals with type 1 diabetes patients and controls. Seventy polymorphisms, including 51 in introns, 4 in the 3' flanking region, and 15 in exons, were identified. Sequencing of polymorphic DNA fragments revealed several new polymorphisms, including a Gln --> Arg substitution at codon 611 and a GT --> GC polymorphism affecting the donor splice site of intron 4, that might be of functional significance. None of the polymorphisms examined differed in frequency between individuals with type 1 diabetes and controls.