Identification of genetic variants that influence atherosclerosis in young persons

Atherosclerosis is widely accepted as a complex genetic phenotype and is the usual cause of cardiovascular disease, the world’s leading killer. Genetic factors have been proven to be important risk contributors for atherosclerosis and much work has been done to identify promising candidates that might play a role in the development of atherosclerosis. It is well known that many independent replications are needed to unequivocally establish a valid genotype-phenotype association across different populations before the findings are extended to clinical settings and to the expensive follow-up studies designed to identify causal genetic variants. Aiming to replicate the association with atherosclerosis in the Pathobiological Determinants of Atherosclerosis in Youth (PDAY) study, we assessed the relationship of 32 atherosclerosis candidate SNPs to atherosclerosis in the PDAY cohort, consisting of AA and EA young people aged 15-34 years who died of non-medical causes. Two association studies, a whole sample study and a 1:1 matched case control study were performed by use of multiple linear regression and logistic regression analyses, respectively. For the whole sample association study, 32 SNPs among 2,650 individuals (1,369 AA and 1,281 EA) were tested for the association with six early atherosclerosis phenotypes: abdominal aorta fatty streaks, abdominal aorta raised lesions, right coronary artery fatty streaks, right coronary artery raised lesions, thoracic aorta fatty streaks, and thoracic aorta raised lesions. For the matched case-control association study, 337 case-control paired samples were included; cases were chosen with the highest total raised lesion scores from the studied population, while controls were randomly selected from individuals that had no raised lesions and matched to cases by age, gender and race. Sixteen SNPs in 13 genes were found to be significantly associated with atherosclerosis in at least one of the PDAY association studies. Among these 16 findings: eight SNPs (rs9579646, rs6053733, rs3849150, rs10499903, rs2148079, rs5073691, rs10116277, and rs17228212) successfully replicated previous results, six SNPs (rs17222814, rs10811661, rs7028570, rs7291467, rs16996148 and rs10401969) were reported as new findings exclusive to our study, the last two of the 16 SNPs, rs501120 and rs6922269, showed either intriguing or conflicting result. SNP rs17222814 in ALOX5AP and SNP rs3849150 in LRRC18 were consistently associated with atherosclerosis in both prior and the two PDAY association studies. SNP rs3849150 was also identified to be highly correlated with a non-synonymous coding SNP, rs17772611, which may damage the protein (polyphen score = 0.996), suggesting that SNP rs17772611 may be the causal functional variant.^ In conclusion, our study added more support for the association of these candidate genes with atherosclerosis. SNPs rs3849150 and rs17772611 of LRRC18, as well as SNP rs17222814 of ALOX5AP, were the most significant findings from our study, and may be ranked among the best for further study.^

Comparison of two statistical methods for rare variant association analysis

Studies have shown that rare genetic variants have stronger effects in predisposing common diseases, and several statistical methods have been developed for association studies involving rare variants. In order to better understand how these statistical methods perform, we seek to compare two recently developed rare variant statistical methods (VT and C-alpha) on 10,000 simulated re-sequencing data sets with disease status and the corresponding 10,000 simulated null data sets. The SLC1A1 gene has been suggested to be associated with diastolic blood pressure (DBP) in previous studies. In the current study, we applied VT and C-alpha methods to the empirical re-sequencing data for the SLC1A1 gene from 300 whites and 200 blacks. We found that VT method obtains higher power and performs better than C-alpha method with the simulated data we used. The type I errors were well-controlled for both methods. In addition, both VT and C-alpha methods suggested no statistical evidence for the association between the SLC1A1 gene and DBP. Overall, our findings provided an important comparison of the two statistical methods for future reference and provided preliminary and pioneer findings on the association between the SLC1A1 gene and blood pressure.^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.

Genetic variants within the Wnt/ B-catenin signaling pathway as indicators of bladder cancer clinical outcomes

The genetic factors that influence bladder cancer clinical outcomes are largely unknown. In this clinical outcomes study, I assessed genetic variations in the Wnt/β-catenin stem-cell pathway genes for association with recurrence and progression. A total of 230 SNPS in 40 genes from the Wnt/β-catenin pathway were genotyped in 419 histologically confirmed non-muscle invasive bladder cancer cases. Several significant associations were observed in the clinical outcomes analysis. Under the dominant model WNT8B: rs4919464 (HR: 1.55, 95% CI: 1.17-2.06, P=2.2x10-3) and WNT8B: rs3793771 (HR: 1.54, 95% CI: 1.09-1.62, P=4.6x10-3 ) were statistically significantly associated with an increase risk of recurrence while two other variants, APC2: rs11668593 (HR: 2.50, 95% CI: 1.43-4.35, P=1.2x10-3) and LRP5 : rs312778 (HR: 1.81, 95% CI: 1.23-2.65, P=2.7x10-3), were significantly associated with recurrence risk under the recessive model of inheritance. Four SNPs in the recessive model were associated with an increased risk of progression (AXIN2: rs1544427, LRP5: rs312778, AXIN1: rs370681, AXIN1: rs2301522). LRP5: rs312778 had the most significant increased risk of progression with a 2.68 (95% CI: 1.52-4.72, P=6.4x10-4)-fold increased risk. Stratification analysis based on treatment regimen (transurethral resection (TUR) and Bacillus Calmette-Guérin (BCG)) was also performed. Individuals with at least one variant in AXIN2: rs2007085 were found to have a 2.09 (95% CI: 1.24-3.52, P=5.4x10-3) -fold increased risk of recurrence in those that received TUR only, and no statistically significant effect was seen in those that received BCG. Individuals who received TUR with at least one variant in LEF1: rs10516550 were found to have a 2.26 (95% CI: 1.22-4.18, P=9.7x10-3)-fold increase risk of recurrence and no statistically significant effect was found in individuals who received BCG. Also, the recessive model of LRP6: rs2302684 in TUR only treatment was shown to have a 1.95 (95%CI: 1.18-3.21, P=8.8x10 -3)-fold increased risk of recurrence, and a suggested protective effect associated with a (HR: 0.83, 95% CI: 0.51-1.37, P=0.468) decreased risk of recurrence. Together, these findings implicate the Wnt/β-catenin stem-cell pathway as playing a role in bladder cancer clinical outcomes and have important implications for personalization of future treatment regimens. ^

Gilbert syndrome and glucose-6-phosphate dehydrogenase deficiency: A dose-dependent genetic interaction crucial to neonatal hyperbilirubinemia

Severe jaundice leading to kernicterus or death in the newborn is the most devastating consequence of glucose-6-phosphate dehydrogenase (EC 1.1.1.49; G-6-PD) deficiency. We asked whether the TA repeat promoter polymorphism in the gene for uridinediphosphoglucuronate glucuronosyltransferase 1 (EC 2.4.1.17; UDPGT1), associated with benign jaundice in adults (Gilbert syndrome), increases the incidence of neonatal hyperbilirubinemia in G-6-PD deficiency. DNA from term neonates was analyzed for UDPGT1 polymorphism (normal homozygotes, heterozygotes, variant homozygotes), and for G-6-PD Mediterranean deficiency. The variant UDPGT1 promoter allele frequency was similar in G-6-PD-deficient and normal neonates. Thirty (22.9%) G-6-PD deficient neonates developed serum total bilirubin ≥ 257 μmol/liter, vs. 22 (9.2%) normals (P = 0.0005). Of those with the normal homozygous UDPGT1 genotype, the incidence of hyperbilirubinemia was similar in G-6-PD-deficients and controls (9.7% and 9.9%). In contrast, in the G-6-PD-deficient neonates, those with the heterozygous or homozygous variant UDPGT1 genotype had a higher incidence of hyperbilirubinemia than corresponding controls (heterozygotes: 31.6% vs. 6.7%, P < 0.0001; variant homozygotes: 50% vs. 14.7%, P = 0.02). Among G-6-PD-deficient infants the incidence of hyperbilirubinemia was greater in those with the heterozygous (31.6%, P = 0.006) or variant homozygous (50%, P = 0.003) UDPGT1 genotype than in normal homozygotes. In contrast, among those normal for G-6-PD, the UDPGT1 polymorphism had no significant effect (heterozygotes: 6.7%; variant homozygotes: 14.7%). Thus, neither G-6-PD deficiency nor the variant UDPGT1 promoter, alone, increased the incidence of hyperbilirubinemia, but both in combination did. This gene interaction may serve as a paradigm of the interaction of benign genetic polymorphisms in the causation of disease.

Naturally variant autosomal and sex-linked loci determine the severity of iron overload in β2-microglobulin-deficient mice

Hereditary hemochromatosis (HH) is a common chronic human genetic disorder whose hallmark is systemic iron overload. Homozygosity for a mutation in the MHC class I heavy chain paralogue gene HFE has been found to be a primary cause of HH. However, many individuals homozygous for the defective allele of HFE do not develop iron overload, raising the possibility that genetic variation in modifier loci contributes to the HH phenotype. Mice deficient in the product of the β2-microglobulin (β2M) class I light chain fail to express HFE and other MHC class I family proteins, and they have been found to manifest many characteristics of the HH phenotype. To determine whether natural genetic variation plays a role in controlling iron overload, we performed classical genetic analysis of the iron-loading phenotype in β2M-deficient mice in the context of different genetic backgrounds. Strain background was found to be a major determinant in iron loading. Sex played a role that was less than that of strain background but still significant. Resistance and susceptibility to iron overload segregated as complex genetic traits in F1 and back-cross progeny. These results suggest the existence of naturally variant autosomal and Y chromosome-linked modifier loci that, in the context of mice genetically predisposed by virtue of a β2M deficiency, can profoundly influence the severity of iron loading. These results thus provide a genetic explanation for some of the variability of the HH phenotype.

Genetic and environmental factors modify bovine spongiform encephalopathy incubation period in mice

The incubation period (IP) and the neuropathology of transmissible spongiform encephalopathies (TSEs) have been extensively used to distinguish prion isolates (or strains) inoculated into panels of inbred mouse strains. Such studies have shown that the bovine spongiform encephalopathy (BSE) agent is indistinguishable from the agent causing variant Creutzfeldt–Jakob disease (vCJD), but differs from isolates of sporadic CJD, reinforcing the idea that the vCJD epidemic in Britain results from consumption of contaminated beef products. We present a mouse model for genetic and environmental factors that modify the incubation period of BSE cross-species transmission. We have used two mouse strains that carry the same prion protein (PrP) allele, but display a 100-day difference in their mean IP following intracerebral inoculation with primary BSE isolate. We report genetic effects on IP that map to four chromosomal regions, and in addition we find significant factors of host environment, namely the age of the host's mother, the age of the host at infection, and an X-cytoplasm interaction in the host.

Enhanced antitumor efficacy of a herpes simplex virus mutant isolated by genetic selection in cancer cells

Replication-competent, attenuated herpes simplex virus-1 (HSV-1) derivatives that contain engineered mutations into the viral γ34.5 virulence gene have been used as oncolytic agents. However, as attenuated mutants often grow poorly, they may not completely destroy some tumors and surviving cancer cells simply regrow. Thus, although HSV-1 γ34.5 mutants can reduce the growth of human tumor xenografts in mice and have passed phase I safety studies, their efficacy is limited because they replicate poorly in many human tumor cells. Previously, we selected for a γ34.5 deletion mutant variant that regained the ability to replicate efficiently in tumor cells. Although this virus contains an extragenic suppressor mutation that confers enhanced growth in tumor cells, it remains attenuated. Here, we demonstrate that the suppressor virus replicates to greater levels in prostate carcinoma cells and, importantly, is a more potent inhibitor of tumor growth in an animal model of human prostate cancer than the γ34.5 parent virus. Thus, genetic selection in cancer cells can be used as a tool to enhance the antitumor activity of a replication-competent virus. The increased therapeutic potency of this oncolytic virus may be useful in the treatment of a wide variety of cancers.

Mutagenesis by third-strand-directed psoralen adducts in repair-deficient human cells: high frequency and altered spectrum in a xeroderma pigmentosum variant.

Psoralen-conjugated triple-helix-forming oligonucleotides have been used to generate site-specific mutations within mammalian cells. To investigate factors influencing the efficiency of oligonucleotide-mediated gene targeting, the processing of third-strand-directed psoralen adducts was compared in normal and repair-deficient human cells. An unusually high mutation frequency and an altered mutation pattern were seen in xeroderma pigmentosum variant (XPV) cells compared with normal, xeroderma pigmentosum group A (XPA), and Fanconi anemia cells. In XPV, targeted mutations were produced in the supF reporter gene carried in a simian virus 40 vector at a frequency of 30%, 3-fold above that in normal or Fanconi anemia cells and 6-fold above that in XPA. The mutations generated by targeted psoralen crosslinks and monoadducts in the XPV cells formed a pattern distinct from that in the other three cell lines, with mutations occurring not just at the damaged site but also at adjacent base pairs. Hence, the XPV cells may have an abnormality in trans-lesion bypass synthesis during repair and/or replication, implicating a DNA polymerase or an accessory factor as a basis of the defect in XPV. These results may help to elucidate the repair deficiency in XPV, and they raise the possibility that genetic manipulation via triplex-targeted mutagenesis may be enhanced by modulation of the XPV-associated activity in normal cells.

Genetic heterogeneity in isogenic homozygous clonal zebrafish.

The C32 isogenic homozygous diploid (IHD) strain of the zebrafish (Danio rerio) was found to be polyallelic at a malate dehydrogenase locus (sMdh-A). A variant allele is thought to have arisen via mutation within the past 10 bisexual generations that have maintained the strain since its last gynogenetic cloning event; this unique allele now predominates at the sMdh-A locus. The estimated mutation rate in this species is sufficiently high that long-term genetic homogeneity of its IHD clones cannot be assumed. Researchers using such bisexually maintained clones should be aware that they are not necessarily using genetically uniform subjects. Genetic uniformity of cloned IHD zebrafish will be maximized if experimental subjects are obtained soon after a cloning event.

Analysis of RNA-binding proteins by in vitro genetic selection: identification of an amino acid residue important for locking U1A onto its RNA target.

An in vitro genetic system was developed as a rapid means for studying the specificity determinants of RNA-binding proteins. This system was used to investigate the origin of the RNA-binding specificity of the mammalian spliceosomal protein U1A. The U1A domain responsible for binding to U1 small nuclear RNA was locally mutagenized and displayed as a combinatorial library on filamentous bacteriophage. Affinity selection identified four U1A residues in the mutagenized region that are important for specific binding to U1 hairpin II. One of these residues (Leu-49) disproportionately affects the rates of binding and release and appears to play a critical role in locking the protein onto the RNA. Interestingly, a protein variant that binds more tightly than U1A emerged during the selection, showing that the affinity of U1A for U1 RNA has not been optimized during evolution.

Next-generation sequencing of iron-metabolism related genes in Portuguese patients with iron overload: novel pathogenic genetic variants

Objective: In Southern European countries up to one-third of the patients with hereditary hemochromatosis (HH) do not present the common HFE risk genotype. In order to investigate the molecular basis of these cases we have designed a gene panel for rapid and simultaneous analysis of 6 HH-related genes (HFE, TFR2, HJV, HAMP, SLC40A1 and FTL) by next-generation sequencing (NGS). Materials and Methods: Eighty-eight iron overload Portuguese patients, negative for the common HFE mutations, were analysed. A TruSeq Custom Amplicon kit (TSCA, by Illumina) was designed in order to generate 97 amplicons covering exons, intron/exon junctions and UTRs of the mentioned genes with a cumulative target sequence of 12115bp. Amplicons were sequenced in the MiSeq instrument (IIlumina) using 250bp paired-end reads. Sequences were aligned against human genome reference hg19 using alignment and variant caller algorithms in the MiSeq reporter software. Novel variants were validated by Sanger sequencing and their pathogenic significance were assessed by in silico studies. Results: We found a total of 55 different genetic variants. These include novel pathogenic missense and splicing variants (in HFE and TFR2), a very rare variant in IRE of FTL, a variant that originates a novel translation initiation codon in the HAMP gene, among others. Conclusion: The merging of TSCA methodology and NGS technology appears to be an appropriate tool for simultaneous and fast analysis of HH-related genes in a large number of samples. However, establishing the clinical relevance of NGS-detected variants for HH development remains a hard-working task, requiring further functional studies.

Construction of a New Familial Hypercholesterolemia Variant Data Base. A Systematic Review for a 2015 Update

Aims: Familial hypercholesterolemia (FH) is an autosomal dominant disorder with increased cardiovascular risk, caused by mutations in LDLR, APOB and PCSK 9 genes. Although it is described that over 1700 variants have been found, none of the existing databases are completely updated. The aim of this work is to construct a FH database in order to provide a unique source of verified information about variants associated with FH for a more accurate genetic diagnosis.

Hospitalization risk due to respiratory illness associated with genetic variation at IFITM3 in patients with influenza A(H1N1)pdm09 infection: a case-control study

Recent studies suggest an association between the Interferon Inducible Transmembrane 3 (IFITM3) rs12252 variant and the course of influenza infection. However, it is not clear whether the reported association relates to influenza infection severity. The aim of this study was to estimate the hospitalization risk associated with this variant in Influenza Like Illness (ILI) patients during the H1N1 pandemic influenza. A case-control genetic association study was performed, using nasopharyngeal/oropharyngeal swabs collected during the H1N1 pandemic influenza. Laboratory diagnosis of influenza infection was performed by RT-PCR, the IFITM3 rs12252 was genotyped by RFLP and tested for association with hospitalization. Conditional logistic regression was performed to calculate the confounder-adjusted odds ratio of hospitalization associated with IFITM3 rs12252. We selected 312 ILI cases and 624 matched non-hospitalized controls. Within ILI Influenza A(H1N1)pdm09 positive patients, no statistical significant association was found between the variant and the hospitalization risk (Adjusted OR: 0.73 (95%CI: 0.33–1.50)). Regarding ILI Influenza A(H1N1)pdm09 negative patients, CT/CC genotype carriers had a higher risk of being hospitalized than patients with TT genotype (Adjusted OR: 2.54 (95%CI: 1.54–4.19)). The IFITM3 rs12252 variant was associated with respiratory infection hospitalization but not specifically in patients infected with Influenza A(H1N1)pdm09.