Genome sequence, comparative analysis, and population genetics of the domestic horse

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

Regulation of survivin gene expression in the human endometrium and endometrial cancer

In the United States, endometrial cancer is the leading cancer of the female reproductive tract. There are 40,100 new cases and 7,470 deaths from endometrial cancer estimated for 2008 (47). The average five year survival rate for endometrial cancer is 84% however, this figure is substantially lower in patients diagnosed with late stage, advanced disease and much higher for patients diagnosed in early stage disease (47). Endometrial cancer (EC) has been associated with several risk factors including obesity, diabetes, hypertension, previously documented occurrence of hereditary non-polyposis colorectal cancer (HNPCC), and heightened exposure to estrogen (25). As of yet, there has not been a dependable molecular predictor of endometrial cancer occurrence in women with these predisposing factors. The goal of our lab is to identify genes that are aberrantly expressed in EC and may serve as molecular biomarkers of EC progression. One candidate protein that we are exploring as a biomarker of EC progression is the cell survival protein survivin.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

THE ROLE OF MAP KINASES ON THE FUNCTIONAL HETEROGENEITY OF HUMAN CD8+ T CELL MATURATION SUBSETS

While prior studies have focused on naïve (CD45RA+CD27+) and early stage memory (CD45RA-CD27+) CD8+ T cells, late memory CD8+ T cells (CD45RA+CD27) have received less interest because this subset of T cells is generally recognized as effectors, which produce IFNγ (but no IL-2) and perforin. However, multiple studies suggest that late memory CD8+ T cells may provide inadequate protection in infectious diseases and cancer models. To better understand the unique function of late memory CD8+ T cells, I optimized multi-color flow cytometry techniques to assess the cytokine production of each human CD8+ T cell maturation subset. I demonstrated that late memory CD8+ T cells are the predominant producer of CC chemokines (e.g. MIP-1β), but rarely produce IL-2; therefore they do not co-produce IL-2/IFNγ (polyfunctionality), which has been shown to be critical for protective immunity against chronic viral infection. These data suggest that late memory CD8+ T cells are not just cytotoxic effectors, but may have unique functional properties. Determining the molecular signature of each CD8+ T cell maturation subset will help characterize the role of late memory CD8+ T cells. Prior studies suggest that ERK1 and ERK2 play a role in cytokine production including IL-2 in T cells. Therefore, I tested whether differential expression of ERK1 and ERK2 in CD8+ T cell maturation subsets contributes to their functional signature by a novel flow cytometry technique. I found that the expression of total ERK1, but not ERK2, is significantly diminished in late memory CD8+ T cells and that ERK1 expression is strongly associated with IL-2 production and CD28 expression. I also found that IL-2 production is increased in late memory CD8+ T cells by over-expressing ERK1. Collectively, these data suggest that ERK1 is required for IL-2 production in human CD8+ T cells. In summary, this dissertation demonstrated that ERK1 is down-regulated in human late memory CD8+ T cells, leading to decreased production of IL-2. The data in this dissertation also suggested that the functional heterogeneity in human CD8+ T cell maturation subsets results from their differential ERK1 expression.

ROLE OF VITAMIN-D3 AND RETINOIC ACID IN A HUMAN THP-1 MACROPHAGE MODEL OF MYCOBACTERIUM TUBERCULOSIS INFECTION

Mycobacterium tuberculosis (Mtb) replicates within the human macrophages and we investigated the activating effects of retinoic acid (RA) and vitamin D3 (VD) on macrophages in relation to the viability of Mtb. A combination of these vitamins (RAVD) enhanced the receptors on THP-1 macrophage (Mannose receptor and DC-SIGN) that increased mycobacterial uptake but inhibited thesubsequent intracellular growth of Mtb by inducing reactive oxygen species (ROS) and autophagy. RAVD also enhanced antigen presenting and homing receptors in THPs that suggested an activated phenotype for THPs following RAVD treatment. RAVD mediated activation was also associated with a marked phenotypic change in Mtb infected THPs that fused with adjacent cells to formmultinucleate giant cells (MNGCs). Typically MNGCs occurred over 30 days of in vitro culture and contained non-replicating persisting Mtb for as long as 60 days in culture. We propose that the RAVD mediated inhibition of replicating Mtb leading to persistence of non-replicating Mtb within THPs may provide a novel human macrophage model simulating formation of MNGCs in humanlungs.

A NOVEL FUNCTION FOR AURORA B KINASE IN THE REGULATION OF P53 BY PHOSPHORYLATION

The mitotic kinase Aurora B plays a pivotal role in mitosis and cytokinesis and governs the spindle assembly checkpoint which ensures correct chromosome segregation and normal progression through mitosis. Aurora B is overexpressed in breast and other cancers and may be an important molecular target for chemotherapy. Tumor suppressor p53 is the guardian of the genome and an important negative regulator of the cell cycle. Previously, it was unknown whether Aurora B and p53 had mutual regulation during the cell cycle. A small molecule specific inhibitor of Aurora B, AZD1152, gave us an indication that Aurora B negatively impacted p53 during interphase and mitosis. Here, we show the antineoplastic activity of AZD1152 in six human breast cancer cell lines, three of which overexpress HER2. AZD1152 specifically inhibited Aurora B kinase activity, thereby causing mitotic catastrophe, polyploidy and apoptosis, which in turn led to apoptotic death. Further, AZD1152 administration efficiently suppressed tumor growth in orthotopic and metastatic breast cancer cell xenograft models. Notably, it was found that the protein level of Aurora B kinase declined after inhibition of Aurora B kinase activity. Investigation of the underlying mechanism suggested that AZD1152 accelerated the protein turnover of Aurora B by enhancing its ubiquitination. As a consequence of inhibition of Aurora B, p53 levels were increased in tissue culture and murine models. This hinted at a possible direct interaction between p53 and Aurora B. Indeed, it was found that p53 and Aurora B exist in complex and interact directly during interphase and at the centromere in mitosis. Further, Aurora B was shown to phosphorylate p53 at several serine/threonine residues in the DNA binding domain and these events caused downregulation of p53 levels via ubiquitination mediated by Mdm2. Importantly, phosphorylation of threonine 211 was shown to reduce p53’s transcriptional activity while other phosphorylation sites did not. On a functional level, Aurora B was shown to reduce p53’s capacity to mediate apoptosis in response to the DNA damaging agent, cisplatin. These results define a novel mechanism for p53 inactivation by Aurora B and imply that oncogenic hyperactivation or overexpression of Aurora B may compromise p53’s tumor suppressor function.

Genome sequence of Fusobacterium nucleatum subspecies polymorphum - a genetically tractable fusobacterium.

Fusobacterium nucleatum is a prominent member of the oral microbiota and is a common cause of human infection. F. nucleatum includes five subspecies: polymorphum, nucleatum, vincentii, fusiforme, and animalis. F. nucleatum subsp. polymorphum ATCC 10953 has been well characterized phenotypically and, in contrast to previously sequenced strains, is amenable to gene transfer. We sequenced and annotated the 2,429,698 bp genome of F. nucleatum subsp. polymorphum ATCC 10953. Plasmid pFN3 from the strain was also sequenced and analyzed. When compared to the other two available fusobacterial genomes (F. nucleatum subsp. nucleatum, and F. nucleatum subsp. vincentii) 627 open reading frames unique to F. nucleatum subsp. polymorphum ATCC 10953 were identified. A large percentage of these mapped within one of 28 regions or islands containing five or more genes. Seventeen percent of the clustered proteins that demonstrated similarity were most similar to proteins from the clostridia, with others being most similar to proteins from other gram-positive organisms such as Bacillus and Streptococcus. A ten kilobase region homologous to the Salmonella typhimurium propanediol utilization locus was identified, as was a prophage and integrated conjugal plasmid. The genome contains five composite ribozyme/transposons, similar to the CdISt IStrons described in Clostridium difficile. IStrons are not present in the other fusobacterial genomes. These findings indicate that F. nucleatum subsp. polymorphum is proficient at horizontal gene transfer and that exchange with the Firmicutes, particularly the Clostridia, is common.

HLA-DPB1 and DPB2 are genetic loci for systemic sclerosis: a genome-wide association study in Koreans with replication in North Americans.

OBJECTIVE: To identify systemic sclerosis (SSc) susceptibility loci via a genome-wide association study. METHODS: A genome-wide association study was performed in 137 patients with SSc and 564 controls from Korea using the Affymetrix Human SNP Array 5.0. After fine-mapping studies, the results were replicated in 1,107 SSc patients and 2,747 controls from a US Caucasian population. RESULTS: The single-nucleotide polymorphisms (SNPs) (rs3128930, rs7763822, rs7764491, rs3117230, and rs3128965) of HLA-DPB1 and DPB2 on chromosome 6 formed a distinctive peak with log P values for association with SSc susceptibility (P=8.16x10(-13)). Subtyping analysis of HLA-DPB1 showed that DPB1*1301 (P=7.61x10(-8)) and DPB1*0901 (P=2.55x10(-5)) were the subtypes most susceptible to SSc in Korean subjects. In US Caucasians, 2 pairs of SNPs, rs7763822/rs7764491 and rs3117230/rs3128965, showed strong association with SSc patients who had either circulating anti-DNA topoisomerase I (P=7.58x10(-17)/4.84x10(-16)) or anticentromere autoantibodies (P=1.12x10(-3)/3.2x10(-5)), respectively. CONCLUSION: The results of our genome-wide association study in Korean subjects indicate that the region of HLA-DPB1 and DPB2 contains the loci most susceptible to SSc in a Korean population. The confirmatory studies in US Caucasians indicate that specific SNPs of HLA-DPB1 and/or DPB2 are strongly associated with US Caucasian patients with SSc who are positive for anti-DNA topoisomerase I or anticentromere autoantibodies.

Expression and Regulation of Human Cytochrome P450 4F Isoforms in Tissue Samples and Under TNF-Alpha Challenges

CYP4F (Cytochrome P4504F) enzymes metabolize endogenous molecules including leukotrienes, prostaglandins and arachidonic acid. The involvement of these endogenous compounds in inflammation has led to the hypothesis that changes in the inflamed tissue environment may affect the expression of CYP4Fs during the pro-inflammatory state, which in turn may modulate inflammatory conditions during the anti-inflammatory state. We demonstrated that inflamed tissues have different levels of CYP4F isoform expression profiles in a number of human samples when compared to the average population. The CYP4F isoform expression levels change with the degree of inflammation present in tissue. Further investigation in cell culture studies revealed that inflammatory cytokines, in particular TNF-α, play a role in regulating the expression of the CYP4F family. One of the isoforms, CYP4F11, had different characteristics than that of the other five CYP4F family members. CYP4F11 metabolizes xenobiotics while the other isoforms metabolize endogenous compounds with higher affinity. CYP4F11 also was expressed at high quantities in the brain, and was up-regulated by TNF-α, while the other isoforms were not expressed at high quantities in the brain and were down-regulated by TNF-α. We identified the AP-1 protein of the JNK pathway as the signaling protein that causes significant increase in CYP4F11 expression. Since TNF-α stimulation causes a simultaneous activation of both JNK pathway and NF-κB signaling, we investigated further the role that NF-κB plays on expression of the CYP4F11 gene. We concluded that although there is a significant increase in CYP4F11 expression in the presence of TNF-α, the activation of NF-κB signaling inhibits CYP4F11 expression in a time dependent manner. The expression of CYP4F11 is only significantly increased after 24 hours of treatment with TNF-α; at shorter time points NF-κB signaling overpowers the JNK pathway activation. We believe that these findings may in the future lead to improved drug design for modulating inflammation.

A Cell Biological Determination of Integrator Subunit Localization

Uridine-rich small nuclear (U snRNAs), with the exception of the U6 snRNA, are RNA polymerase II (RNAPII) transcripts. The mechanism of 3’ cleavage of snRNAs has been unknown until recently. This area was greatly advanced when 12 of the Integrator complex subunits (IntS) were purified in 2005 through their interaction with the C-terminal domain (CTD) of the large subunit (RpbI) of RNAPII. Subsequently, our lab performed a genome-wide RNAi screen that identified two more members of the complex that we have termed IntS13 and IntS14. We have determined that IntS9 and 11 mediate the 3’ cleavage of snRNAs, but the exact function of the other subunits remains unknown. However, through the use of a U7 snRNA-GFP reporter and RNAi knockdown of the Integrator subunits in Drosophila S2 cells, we have shown that all subunits are required for the proper processing of snRNAs, albeit to differing degrees. Because snRNA transcription takes place in the nucleus of the cell, it is expected that all of the Integrator subunits would exhibit nuclear localization, but the knowledge of discrete subnuclear localization (i.e. to Cajal bodies) of any of the subunits could provide important clues to the function of that subunit. In this study, we used a cell biological approach to determine the localization of the 14 Integrator subunits. We hypothesized that the majority of the subunits would be nuclear, however, a few would display distinct localization to the Cajal bodies, as this is where snRNA genes are localized and transcribed. The specific aims and results are: 1. To determine the subcellular localization of the 14 Integrator subunits. To accomplish this, mCherry and GFP tagged clones were generated for each of the 14 Drosophila and human Integrator subunits. Confocal microscopy studies revealed that the majority of the subunits were diffuse in the nucleus, however, IntS3 formed discrete subnuclear foci. Surprisingly, two of the subunits, IntS2 and 7 were observed in cytoplasmic foci. 2. To further characterize Integrator subunits with unique subcellular localizations. Colocalization studies with endogenous IntS3 and Cajal body marker, coilin, showed that these two proteins overlap, and from this we concluded that IntS3 localized to Cajal bodies. Additionally, colocalization studies with mCherry-tagged IntS2 and 7 and the P body marker, Dcp1, revealed that these proteins colocalize as well. IntS7, however, is more stable in cytoplasmic foci than Dcp1. It was also shown through RNAi knockdown of Integrator subunits, that the cytoplasmic localization of IntS2 and 7 is dependent on the expression of IntS1 and 11 in S2 cells.

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

EXPRESSION AND IMMUNOLOGICAL CHARACTERIZATION OF HERV-K TRANSMEMBRANE ENVELOPE PROTEIN IN HUMAN BREAST CANCER

The human endogenous retrovirus K (HERV-K) env gene encodes envelope protein comprising surface (SU) and transmembrane (TM) domains. Having shown the exclusive expression of SU in human breast cancer and the stimulation of SU-specific immune responses in patients with breast cancer, our research here confirmed and extended the data by investigating the expression of HERV-K TM envelope domain and the induction of specific immune responses against TM in breast cancer patients. We found HERV-K TM mRNA and protein expression only in human breast cancer cells but not in normal controls. The specific immune responses against TM domain were induced in mice determined by enzyme-linked immunosorbent assay (ELISA) and IFN-γ enzyme-linked immunosorbent spot (ELISPOT) assay. Furthermore, ELISA detected higher titers of anti-HERV-K TM Env IgG antibodies in sera of breast cancer patients. In addition, the magnitude of the anti-HERV TM B cell response was correlated with the disease stage. Peripheral blood mononuclear cells (PBMCs) before and after in vitro stimulation (IVS) with HERV-K TM from patients with breast cancer as well as healthy controls were tested for T cell responses against HERV-K TM domain by ELISPOT assay. Breast cancer patients (n=21) had stronger HERV-K TM-specific cellular responses than healthy controls (n=12) (P < 0.05). These findings suggest, for the first time, that HERV-K TM expression was enhanced in human breast cancer cells and was able to induce specific B cell and T cell immune responses in breast cancer patients. This study provides support for HERV-K TM as a promising source of antigen for anti-tumor immunotherapy, prevention, diagnosis, and prognosis.