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.^

Association of TNF-alpha genetic variants with HPV-associated oral cancer risk---a case-control study

The incidence of OSCC in younger population and in those who never smoked or drank has increased since the last decade. This increase may be attributable to increase of infection with HPV. The pro-inflammatory cytokine TNF-&agr; has the role in the pathogenesis of chronic inflammatory diseases and was found to control HPV infection in cervical cancer studies. Our study aimed to investigate the association between the four polymorphisms located in TNF-&agr; promoter region, -308(rs1800629), -857(rs1799724), -863(rs1800630) and -1031(rs1799964), and the risk of HPV-related OSCC. In this hospital-based case-control study, 325 cases and 335 controls were included. We found that HPV 16 seropositivity was associated with an increased risk of oral cancer (OR = 3.1, 95% CI, 2.1–4.6). Each of the polymorphism showed to increase the risk of HPV-related OSCC. And after combining the risk genotypes and using the low-risk group (0–1 combined risk genotypes) and HPV16 seronegativity as the reference group, only the high-risk groups (3–4 combined risk genotypes) and HPV16 seronegativity were associated with a low OR of 1.8 (95% CI, 1.1–2.8), while the low-risk and high-risk groups and HPV16 seropositivity were significantly associated with a higher OR of 2.7 (95% CI, 1.3–5.8) and 8.5 (95% CI, 3.7–19.4), respectively. In addition, the joint effects were greater among the young subjects (aged<50), males, never smokers or never drinkers, and patients with oropharyngeal cancer. Overall, the four TNF-&agr; polymorphisms, individually or collectively, would result in a significantly increased risk for HPV16-associated oral cancer in a non-Hispanic white population. More large sized studies are needed for future investigation.^

Genetic variants associated with metastatic prostate cancer

Introduction. Distant metastasis remains the leading cause of death among prostate cancer patients. Several genetic susceptibility loci associated with Prostate cancer have been identified by the Genome Wide Association Studies (GWAS). To date, few studies have explored the ability of these SNPs to identify metastatic prostate cancer. Based on the identification of genetic variants as predictors of aggressive disease, a case comparison study of prostate cancer patients was designed to explore the association of 96 GWAS single nucleotide polymorphisms (SNPs) with metastatic disease. ^ Method. 1242 histologically confirmed prostate cancer patients, with and without metastatic disease, were enrolled into the study. Data were collected from personal interviews, hospital database and abstraction of medical records. Ninety six SNPs identified from GWAS studies based on their associations with prostate cancer risk were genotyped in the study population. Univariate and multivariate logistic regression analyses were used to explore the relationships of the variants with metastatic prostate cancer in Whites and African American men. ^ Results. Four SNPs showed independent associations with metastatic prostate cancer (rs721048 in EHBP1 (2p15), rs3025039 in VEGF (6p12), rs11228565 in Intergenic(11q13.2) and rs2735839 in KLK3(19q13.33)) in the White population. For SNP rs2735839 in KLK3, genotype GA was 1.71 times as likely to be associated with metastatic prostate cancer diagnosis as genotype AA after adjusting for other significant SNPs and covariates (95% CI, 1.12-2.60; p=0.012). In men of African descent, three SNPs: rs1512268 in NKX3-1(8p21.2), rs12155172 in intergenic (7p15.3) & rs10486567 in JAZF1 (7p15.2) were positively associated with metastatic disease in the multivariate analysis. The strongest SNP was rs1512268 heterozygous genotype AG in NKX3-1(8p21.2) which was associated with 3.97-fold increased risk of metastatic prostate cancer diagnosis (95% CI, 1.69-9.34; p =0.002). ^ Conclusion. Genetic variants associated with metastatic prostate cancer were different in Whites and African American men. Given the high mortality rate recorded in men diagnosed with metastatic prostate tumor, further studies are needed to validate associations and establish their clinical application.^

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. ^

GENETIC INTERACTIONS OF FACTORS THAT REGULATE ALTERNATIVE RNA SPLICING IN THE MALE GERM LINE OF DROSOPHILA

Alternative RNA splicing is a critical process that contributes variety to protein functions, and further controls cell differentiation and normal development. Although it is known that most eukaryotic genes produce multiple transcripts in which splice site selection is regulated, how RNA binding proteins cooperate to activate and repress specific splice sites is still poorly understood. In addition how the regulation of alternative splicing affects germ cell development is also not well known. In this study, Drosophila Transformer 2 (Tra2) was used as a model to explore both the mechanism of its repressive function on its own pre-mRNA splicing, and the effect of the splicing regulation on spermatogenesis in testis. Half-pint (Hfp), a protein known as splicing activator, was identified in an S2 cell-based RNAi screen as a co-repressor that functions in combination with Tra2 in the splicing repression of the M1 intron. Its repressive splicing function is found to be sequence specific and is dependent on both the weak 3’ splice site and an intronic splicing silencer within the M1 intron. In addition we found that in vivo, two forms of Hfp are expressed in a cell type specific manner. These alternative forms differ at their amino terminus affecting the presence of a region with four RS dipeptides. Using assays in Drosophila S2 cells, we determined that the alternative N terminal domain is necessary in repression. This difference is probably due to differential localization of the two isoforms in the nucleus and cytoplasm. Our in vivo studies show that both Hfp and Tra2 are required for normal spermatogenesis and cooperate in repression of M1 splicing in spermatocytes. But interestingly, Tra2 and Hfp antagonize each other’s function in regulating germline specific alternative splicing of Taf1 (TBP associated factor 1). Genetic and cytological studies showed that mutants of Hfp and Taf1 both cause similar defects in meiosis and spermatogenesis. These results suggest Hfp regulates normal spermatogenesis partially through the regulation of taf1 splicing. These observations indicate that Hfp regulates tra2 and taf1 activity and play an important role in germ cell differentiation of male flies.

IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.

Molecular Mechanisms of Vascular Disease in Patients with Rare Variants in MYH11

Thoracic aortic aneurysms and dissections (TAAD) are the primary disease affecting the thoracic ascending aorta, with an incidence rate of 10.4/100,000. Although about 20% of patients carry a mutation in a single gene that causes their disease, the remaining 80% of patients may also have genetic factors that increase their risk for developing TAAD. Many of the genes that predispose to TAAD encode proteins involved in smooth muscle cell (SMC) contraction and the disease-causing mutations are predicted to disrupt contractile function. SMCs are the predominant cell type in the ascending aortic wall. Mutations in MYH11, encoding the smooth muscle specific myosin heavy chain, are a rare cause of inherited TAAD. However, rare but recurrent non-synonymous variants in MYH11 are present in the general population but do not cause inherited TAAD. The goal of this study was to assess the potential role of these rare variants in vascular diseases. Two distinct variants were selected: the most commonly seen rare variant, MYH11 R247C, and a duplication of the chromosomal region spanning the MYH11 locus at 16p13.1. Genetic analyses indicated that both of these variants were significantly enriched in patients with TAAD compared with controls. A knock-in mouse model of the Myh11 R247C rare variant was generated, and these mice survive and reproduce normally. They have no structural abnormalities of the aorta or signs of aortic disease, but do have decreased aortic contractility. Myh11R247C/R247C mice also have increased proliferative response to vascular injury in vivo and increased proliferation of SMCs in vitro. Myh11R247C/R247C SMCs have decreased contractile gene and protein expression and are dedifferentiated. In fibroblasts, myosin force generation is required for maturation of focal adhesions, and enhancers of RhoA activity replace enhancers of Rac1 activity as maturation occurs. Consistent with these previous findings, focal adhesions are smaller in Myh11R247C/R247C SMCs, and there is decreased RhoA activation. A RhoA activator (CN03) rescues the dedifferentiated phenotype of Myh11R247C/R247C SMCs. Myh11R247C/R247C mice were bred with an existing murine model of aneurysm formation, the Acta2-/- mouse. Over time, mice carrying the R247C allele in conjunction with heterozygous or homozygous loss of Acta2 had significantly increased aortic diameter, and a more rapid accumulation of pathologic markers. These results suggest that the Myh11 R247C rare variant acts as a modifier gene increasing the risk for and severity of TAAD in mice. In patients with 16p13.1 duplications, aortic MYH11 expression is increased, but there is no corresponding increase in smooth muscle myosin heavy chain protein. Using SMCs that overexpress Myh11, we identified alterations in SMC phenotype leading to excessive protein turnover. All contractile proteins, not just myosin, are affected, and the proteins are turned over by autophagic degradation. Surprisingly, these cells are also more contractile compared with wild-type SMCs. The results described in this dissertation firmly establish that rare variants in MYH11 significantly affect the phenotype of SMCs. Further, the data suggests that these rare variants do increase the risk of TAAD via pathways involving altered SMC phenotype and contraction. Therefore, this study validates that these rare genetic variants alter vascular SMCs and provides model systems to explore the contribution of rare variants to disease.

Immunologic characteristics of immunogenic variants generated by {\it in vitro\/} treatment of a methylcholanthrene-induced murine fibrosarcoma MCA-F with 1 methyl-3-nitro-1-nitrosoguanidine, 5 -aza-2$\sp\prime$-deoxycytidine, or ultraviolet radiation

This study addresses the questions of whether the frequency of generation and in vivo cross-reactivity of highly immunogenic tumor clones induced in a single parental murine fibrosarcoma cell line MCA-F is more closely related to the agent used to induce the Imm$\sp{+}$ clone or whether these characteristics are independent of the agents used. These questions were addressed by treating the parental tumor cell line MCA-F with UV-B radiation (UV-B), 1-methyl-3-nitro-1-nitrosoguanidine (MNNG), or 5-aza-2$\sp\prime$-deoxycytidine (5-azaCdR). The frequency of Imm$\sp{+}$ variant generation was similarly high for the three different agents, suggesting that the frequency of Imm$\sp{+}$ generation was related more closely to the cell line than to the inducing agent used. Cross-reactivity was tested with two Imm$\sp{+}$ clones from each treatment group in a modified immunoprotection assay that selectively engendered antivariant, but not antiparental immunity. Under these conditions each clone, except one, immunized against itself. The MNNG-induced clones engendered stronger antivariant immunity but a weaker variant cross-reactive immunity could also be detected.^ This study also characterized the lymphocyte populations responsible for antivariant and antiparental immunity in vivo. Using the local adoptive transfer assay (LATA) and antibody plus complement depletion of T-cell subsets, we showed that immunity induced by the Imm$\sp{+}$ variants against the parent MCA-F was transferred by the Thy1.2$\sp{+}$, L3T4a$\sp{+}$, Lyt2.1$\sp{-}$ (CD4$\sp{+}$) population, without an apparent contribution by Thy1.2$\sp{+}$, L3T4a$\sp{-}$, Lyt2.1$\sp{+}$ (CD8$\sp{+}$) cells. A role for Lyt2.1$\sp{+}$T lymphocytes in antivariant, but not antiparent immunity was supported by the results of LATA and CTL assays. Immunization with low numbers of viable Imm$\sp{+}$ cells, or with high numbers of non viable Imm$\sp{+}$ cells engendered only antivariant immunity without parental cross-protection. The associative recognition of parental antigens and variant neoantigens resulting in strong antiparent immunity was investigated using somatic cells hybrids of Imm$\sp{+}$ variants of MCA-F and an antigenically distinct tumor MCA-D. An unexpected result of these latter experiments was the expression of a unique tumor-specific antigen by the hybrid cells. These studies demonstrate that the parental tumor-specific antigen and the variant neoantigen must be coexpressed on the cell surface to engender parental cross-protective immunity. (Abstract shortened with permission of author.) ^

{\it cis\/} -acting determinants in the control of MuSVts110 RNA splicing

Like other simple retroviruses the murine sarcoma virus ts110 (MuSVts110) displays an inefficient mode of genome splicing. But, unlike the splicing phenotypic of other retroviruses, the splicing event effected upon the transcript of MuSVts110 is temperature sensitive. Previous work in this laboratory has established that the conditionally defective nature of MuSVts110 RNA splicing is mediated in cis by features in the viral transcript. Here we show that the 5$\sp\prime$ splice site of the MuSVts110 transcript acts as a point of control of the overall splicing efficiency at both permissive and nonpermissive temperatures for splicing. We strengthened and simultaneously weakened the nucleotide structure of the 5$\sp\prime$ splice site in an attempt to elucidate the differential effects each of the two known critical splicing components which interact with the 5$\sp\prime$ splice site have on the overall efficiency of intron excision. We found that a transversion of the sixth nucleotide, resulting in the formation of a near-consensus 5$\sp\prime$ splice site, dramatically increased the overall efficiency of MuSVts110 RNA splicing and abrogated the thermosensitive nature of this splicing event. Various secondary mutations within this original transversion mutant, designed to selectively decrease specific splicing component interactions, lead to recovery of inefficient and thermosensitive splicing. We have further shown that a sequence of 415 nucleotides lying in the downstream exon of the viral RNA and hypothesized to act as an element in the temperature-dependent inhibition of splicing displays a functional redundancy throughout its length; loss and/or replacement of any one sequence of 100 nucleotides within this sequence does not, with one exception detailed below, diminish the degree to which MuSVts110 RNA is inhibited to splice at the restrictive temperature. One specific deletion, though, fortuitously juxtaposed and activated cryptic consensus splicing signals for the excision of a cryptic intron within the downstream exon and markedly potentiated--across a newly defined cryptic exon--the splicing event effected upon the upstream, native intron. We have exploited this mutant of MuSVts110 to further an understanding of the process of exon definition and intron definition and show that the polypyrimidine tract and consensus 3$\sp\prime$ splice site, as well as the 5$\sp\prime$ splice site, within the intron at the 3$\sp\prime$ flank of the defined exon are required for the exon's definition; implying that definition of the downstream intron is required for the in vivo definition of the proximal, upstream exon. Finally; we have shown, through the construction of heterologous mutants of MuSVts110 employing a foreign 3$\sp\prime$ end-forming sequence, that efficiency of transcript splicing can be increased--to a degree which abrogates its thermosensitive nature--in direct proportion to increasing proximity of the 3$\sp\prime$ end-forming signal to the terminal 3$\sp\prime$ splice site. ^

The role of the conserved N -terminal domain of STAT3 in STAT3 dephosphorylation and nuclear export

Rapid redistribution of STAT subcellular localization is an essential feature of cytokine signaling. To elucidate the molecular basis of STAT3 function, which plays a critical role in controlling innate immune responses in vivo, we initiated studies to determine the mechanisms controlling STAT3 nuclear trafficking. We found that STAT3 is transported to the nucleus in the absence of cytokine treatment, as judged by indirect immunofluorescence studies in the presence of leptomycin B, an inhibitor of CRM1-dependent nuclear export, suggesting that the non-phosphorylated STAT3 protein contains a functional nuclear import signal. An isoform lacking the STAT3 N-terminal domain (Δ133STAT3) retains the ability to undergo constitutive nuclear localization, indicating that this region is not essential for cytokine-independent nuclear import. Δ133STAT3 is also transported to the nucleus following stimulation with interleukin-6 (IL-6). Interestingly, IL-6-dependent tyrosine phosphorylation of Δ133STAT3 appears to be prolonged and the nuclear export of the protein delayed in cells expressing endogenous STAT3, consistent with defective Δ133STAT3 dephosphorylation. Endogenous STAT3 does not promote the nuclear export of Δ133STAT3, although dimerization between endogenous Stat3 and Δ133STAT3 is detected readily. Thus, the STAT3 N-terminal domain is not required for dimerization with full-length STAT3, yet appears to play a role in proper export of Stat3 from the nucleus following cytokine stimulation. STAT3-deficient cells reconstituted with Δ133STAT3 show enhanced and prolonged Stat1 signaling in response to IL-6, suggesting that induction of the STAT3-dependent negative regulator SOCS3 is impaired. In fact, Δ133STAT3 fails to induce SOCS3 mRNA efficiently. These studies collectively indicate that the STAT3 N-terminal region may be important for IL-6-dependent target gene activation and nuclear dephosphorylation, while dispensable for nuclear import. STAT3 is an oncogene. STAT3 is constitutively activated in primary tumors of many types. Thus far, research in the design of STAT3 protein inhibitors has focused on the SH2 and DNA-binding domains of STAT3. Interference with these domains eliminates all signaling through STAT3. If the N-terminal domain is involved in tetramerization on a subset of target genes, inhibition of this region may lead to a more selective inhibition of some STAT3 functions while leaving others intact. ^