Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas

To identify genetic susceptibility loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas completed detailed physical and ophthalmologic examinations including fundus photography for diabetic retinopathy grading. 103 individuals with moderate-to-severe non-proliferative diabetic retinopathy or proliferative diabetic retinopathy were defined as cases for this study. DNA samples extracted from study subjects were genotyped using the Affymetrix GeneChip® Human Mapping 100K Set, which includes 116,204 single nucleotide polymorphisms (SNPs) across the whole genome. Single-marker allelic tests and 2- to 8-SNP sliding-window Haplotype Trend Regression implemented in HelixTreeTM were first performed with these direct genotypes to identify genes/regions contributing to the risk of severe diabetic retinopathy. An additional 1,885,781 HapMap Phase II SNPs were imputed from the direct genotypes to expand the genomic coverage for a more detailed exploration of genetic susceptibility to diabetic retinopathy. The average estimated allelic dosage and imputed genotypes with the highest posterior probabilities were subsequently analyzed for associations using logistic regression and Fisher's Exact allelic tests, respectively. To move beyond these SNP-based approaches, 104,572 directly genotyped and 333,375 well-imputed SNPs were used to construct genetic distance matrices based on 262 retinopathy candidate genes and their 112 related biological pathways. Multivariate distance matrix regression was then used to test hypotheses with genes and pathways as the units of inference in the context of susceptibility to diabetic retinopathy. This study provides a framework for genome-wide association analyses, and implicated several genes involved in the regulation of oxidative stress, inflammatory processes, histidine metabolism, and pancreatic cancer pathways associated with severe diabetic retinopathy. Many of these loci have not previously been implicated in either diabetic retinopathy or diabetes. In summary, CDC73, IL12RB2, and SULF1 had the best evidence as candidates to influence diabetic retinopathy, possibly through novel biological mechanisms related to VEGF-mediated signaling pathway or inflammatory processes. While this study uncovered some genes for diabetic retinopathy, a comprehensive picture of the genetic architecture of diabetic retinopathy has not yet been achieved. Once fully understood, the genetics and biology of diabetic retinopathy will contribute to better strategies for diagnosis, treatment and prevention of this disease.^

Expression study of the atherosclerosis mouse aorta reveals significant disturbance of calcium signaling pathway

Atherosclerosis is a complex disease resulting from interactions of genetic and environmental risk factors leading to heart failure and stroke. Using an atherosclerotic mouse model (ldlr-/-, apobec1-/- designated as LDb), we performed microarray analysis to identify candidate genes and pathways, which are most perturbed in changes in the following risk factors: genetics (control C57BL/6 vs. LDb mice), shearstress (lesion-prone vs. lesion-resistant regions in LDb mice), diet (chow vs. high fat fed LDb mice) and age (2-month-old vs. 8-month old LDb mice). ^ Atherosclerotic lesion quantification and lipid profile studies were performed to assess the disease phenotype. A microarray study was performed on lesion-prone and lesion-resistant regions of each aorta. Briefly, 32 male C57BL/6 and LDb mice (n =16/each) were fed on either chow or high fat diet, sacrificed at 2- and 8-months old, and RNA isolated from the aortic lesion-prone and aortic lesion-resistant segments. Using 64 Affymetrix Murine 430 2.0 chips, we profiled differentially expressed genes with the cut off value of FDR ≤ 0.15 for t-test, and q <0.0001 for the ANOVA. The data were normalized using two normalization methods---invariant probe sets (Loess) and Quantile normalization, the statistical analysis was performed using t-tests and ANOVA, and pathway characterization was done using Pathway Express (Wayne State). The result identified the calcium signaling pathway as the most significant overrepresented pathway, followed by focal adhesion. In the calcium signaling pathway, 56 genes were found to be significantly differentially expressed out of 180 genes listed in the KEGG calcium signaling pathway. Nineteen of these genes were consistently identified by both statistical tests, 11 of which were unique to the test, and 26 were unique to the ANOVA test, using the cutoffs noted above. ^ In conclusion, this finding suggested that hypercholesterolemia drives the disease progression by altering the expression of calcium channels and regulators which subsequently results in cell differentiation, growth, adhesion, cytoskeletal change and death. Clinically, this pathway may serve as an important target for future therapeutic intervention, and thus the calcium signaling pathway may serve as an important target for future diagnostic and therapeutic intervention. ^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

GENETIC ANALYSIS OF THE HIPPO PATHWAY IN MOUSE LIVER

Cancer therapy and tumor treatment remain unsolved puzzles. Genetic screening for tumor suppressor genes in Drosophila revealed the Hippo-signaling pathway as a kinase cascade consisting of five core components. Disrupting the pathway by deleting the main component genes breaks the balance of cell proliferation and apoptosis and results in epithelial tissue tumorigenesis. The pathway is therefore believed to be a tumor suppressor pathway. However, a corresponding role in mammals is yet to be determined. Our lab began to investigate the tumor suppression function of the potent mammalian Hippo pathway by putting floxed alleles into the mouse genome flanking the functional-domain-expressing exons in each component (Mst1, Mst2, Sav1, Lats1 and Lats2). These mice were then crossed with different cre-mouse lines to generate conditional knockout mice. Results indicate a ubiquitous tumor suppression function of these components, predominantly in the liver. A further liver specific analysis of the deletion mutation of these components, as well as the Yap/Taz double deletion mutation, reveals essential roles of the Hippo pathway in regulating hepatic quiescence and embryonic liver development. One of the key cellular mechanisms for the Hippo pathway’s involvement in these liver biological events is likely its cell cycle regulation function. Our work will help to develop potential therapeutic approaches for liver cancer.

Role of obesity, weight gain, physical activity, and polymorphisms in the mTOR pathway and the risk of renal cell carcinoma

The interplay between obesity, physical activity, weight gain and genetic variants in mTOR pathway have not been studied in renal cell carcinoma (RCC). We examined the associations between obesity, weight gain, physical activity and RCC risk. We also analyzed whether genetic variants in the mTOR pathway could modify the association. Incident renal cell carcinoma cases and healthy controls were recruited from the University of Texas MD Anderson Cancer Center in Houston, Texas. Cases and controls were frequency-matched by age (±5 years), ethnicity, sex, and county of residence. Epidemiologic data were collected via in-person interview. A total of 577 cases and 593 healthy controls (all white) were included. One hundred ninety-two (192) SNPs from 22 genes were available and their genotyping data were extracted from previous genome-wide association studies. Logistic regression and regression spline were performed to obtain odds ratios. Obesity at age 20, 40, and 3 years prior to diagnosis/recruitment, and moderate and large weight gain from age 20 to 40 were each significantly associated with increased RCC risk. Low physical activity was associated with a 4.08-fold (95% CI: 2.92-5.70) increased risk. Five single nucleotide polymorphisms (SNPs) were significantly associated with RCC risk and their cumulative effect increased the risk by up to 72% (95% CI: 1.20-2.46). Strata specific effects for weight change and genotyping cumulative groups were observed. However, no interaction was suggested by our study. In conclusion, energy balance related risk factors and genetic variants in the mTOR pathway may jointly influence susceptibility to RCC. ^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

INVESTIGATING THE ROLES OF THE p63 ISOFORMS IN THE MICRORNA BIOGENESIS PATHWAY

MicroRNAs play roles in various biological processes like development, tumorigenesis, metastasis and pluripotency. My thesis work has demonstrated roles for p63, a p53 family member, in the upstream regulation of microRNA biogenesis. The p63 gene has a complex gene structure and has multiple isoforms. The TAp63 isoforms contain an acidic transcription activation domain. The ΔNp63 isoforms, lack the TA domain, but have a proline rich region critical for gene transactivation. To understand the functions of these isoforms, the Flores lab generated TAp63 and ΔNp63 conditional knock out mice. Using these mice and tissues and cells from these mice we have found that TAp63 transcriptionally regulates Dicer while ΔNp63 transcriptionally regulates DGCR8. TAp63 -/- mice are highly tumor prone. These mice develop metastatic mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas to distant sites including the liver, lungs, and brain. I found that TAp63 suppresses metastasis by transcriptionally activating Dicer. TAp63 and Dicer levels were very low or lost in high grade human tumors like mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas. Expression of Dicer in these tumor cell lines reduced their invasiveness. Using ΔNp63 -/- mice, I found that ΔNp63 transcriptionally activates DGCR8, resulting in a miRNA profile that is critical to reprogram cells to pluripotency. Analysis of epidermal cells derived from ΔNp63 -/- mice revealed that these cells expressed markers of pluripotency, including Sox2, Oct 4 and Nanog; however, genome-wide analysis revealed a novel profile of genes that are common between ΔNp63 -/- epidermal cells and embryonic stem cells. I also found that mouse cells depleted of ΔNp63 form chimeric mice and teratomas in SCID mice, demonstrating that ΔNp63 deficient cells are pluripotent. Further, I found that restoration of DGCR8 in ΔNp63 -/- epidermal cells reduces their pluripotency and induces terminal differentiation. I also demonstrated that iMS (induced multipotent stem) cells could be generated using human keratinocytes by knockdown of ∆Np63 or DGCR8. Taken together, my work has placed p63 and its isoforms at a critical node in controlling miRNA biogenesis.

A Genomic Approach to Identify The Notch Pathway as a Putative Tumor Suppressor in Endometrial Cancer

Endometrial cancer is the most common gynecological malignancy and the fourth most frequently diagnosed cancer among women. The molecular changes that distinguish normal endometrium from endometrial carcinoma are not thoroughly understood. Identification of these changes could potentially aid in identifying at-risk women who are especially prone to develop endometrial cancer, such as obese women and women with Lynch Syndrome. A microarray analysis was performed using normal endometrium from thin and obese women and cancerous endometrium from obese women. We validated the differential expression of ten genes whose expression was significantly up-regulated or down-regulated using qRT-PCR. All of the genes had distinct expression levels depending on the endometrial carcinoma histotype. As a result, they could serve as molecular markers to distinguish between normal endometrium and endometrial cancer, as well as between low grade endometrial carcinomas and high grade endometrial carcinomas. Two of the ten genes validated, HEYL and HES1, are down-stream targets of the Notch signaling pathway. HEYL and HES1 were identified by microarray and qRT-PCR to have a significant decrease in expression in endometrial carcinomas compared to normal endometrium. We further analyzed the differential expression of other components of the Notch signaling pathway, Notch4 and Jagged1. They were also identified by qRT-PCR to be significantly down-regulated in endometrial carcinomas compared to normal endometrium. Therefore, we believe the Notch signaling pathway to act as a tumor suppressor in endometrial carcinomas.

The p53 tumor suppressor pathway in culture and {\it in vivo\/}

p53 functions as a tumor suppressor through its ability to initiate either growth arrest or apoptosis in cells which have sustained DNA damage. p53 elicits these cellular phenotypes through its biochemical function as a transcriptional activator. By inducing the expression of a battery of target genes, p53 is able to prevent the propagation of cells with damaged DNA. However, the genes transcriptionally induced by p53 which have been identified to date do not fully explain p53 function. p53 has been demonstrated to activate genes involved in cell cycle inhibition, apoptosis and cell proliferation. The reasons for simultaneous activation of p53 targets with disparate, opposing functions are not clear, but may be due to the use of transformed cell lines in previous experiments. In the studies presented in this thesis, the pathway of p53 tumor suppression has been studied in detail in two systems chosen for their relevance to the natural cell environment. One utilizes a normal, unaltered cultured cell system; the other the whole mouse. In order to better understand the role of the known p53 targets in effecting p53 function in normal cells, early rat embryo fibroblasts were irradiated with ultraviolet light to induce DNA damage. It was discovered that p53 protein levels increased in response to irradiation. The known targets of p53, namely, $p21\sp{WAF1/CIP1},\ mdm2,\ cyclin\ G,$ and bax, were shown for the first time to have a differential temporal induction. The growth suppressor $p21\sp{WAF1/CIP1}$ was induced first, followed by cyclin G then mdm2, which is involved in proliferation through its inactivation of p53, and finally, the apoptosis promoter, bax. These findings indicated that p53 activates its target genes in a manner to allow maximum effectiveness of target function. The rat embryo fibroblasts were shown to undergo apoptosis 24 h after irradiation. Additionally, investigation of these cells for cell cycle alterations demonstrated a brief arrest in G1. In the second study, thymocytes from mice with wild type p53 were shown to undergo apoptosis and activate p53 target genes upon ionizing radiation treatment, while thymocytes from mice deficient in p53 could not. The p53 target genes mdm2 and fas were tested in vivo for their ability to mediate p53-regulated apoptosis, and were found dispensible for that cellular function. Therefore, the p53 targets identified to date do not fully explain the ability of p53 to function as a tumor suppressor. Potentially, functional redundancy between the known targets would account for the data seen in these experiments. Additionally, identification of additional target genes should add further understanding of the p53 pathway of tumor suppression. ^

Identification and characterization of cellular genes differentially expressed in a subset of tumors with non -functional WT1

Wilms tumor (WT) or nephroblastoma is a genetically heterogeneous pediatric renal tumor that accounts for 6–7% of all childhood cancers in the U.S. WT1, located at 11p13, is the sole WT gene cloned to date. Additional genomic regions containing genes that play a role in the development of Wilms tumor include 11p15, 7p, 16q, 1p, 17q and 19q. This heterogeneity has made it extremely difficult to develop an understanding of the pathways involved in the development of WT, even in the 5–20% of tumors that show mutations at the WT1 locus. My research addresses this gap in our current comprehension of the development of WT. ^ I have used two complementary approaches to extend the current understanding of molecular changes involved in the development of WT. In order to minimize complexities due to genetic heterogeneity, I confined my analysis to the WT1 pathway by assessing those genetically defined tumors that carry WT1 mutations. WT1 encodes a zinc finger transcription factor, and in vitro studies have identified many genes that are potentially regulated in vivo by WT1. However, there is very little in vivo data that suggests that they are transcriptionally regulated endogenously by WT1. In one approach I assessed the role of WT1 in the in vivo regulation of PDGFA and IGF2, two genes that are strong contenders for endogenous regulation by WT1. Using primary tissue samples, I found no correlation between the level of RNA expression of WT1 with either PDGFA or IGF2, suggesting that WT1 does not play a critical role in their expression in either normal kidney or WT. ^ In a parallel strategy, using differential display analysis I compared global gene expression in a subset of tumors with known homozygous inactivating WT1 mutations (WT1-tumors) to the gene expression in a panel of appropriate control tissues (fetal kidney, normal kidney, rhabdoid tumor and pediatric renal cell carcinoma). Transcripts that are aberrantly expressed in this subset of Wilms tumors are candidates for endogenous transcriptional regulation by WT1 as well as for potentially functioning in the development of WT. By this approach I identified several differentially expressed transcripts. I further characterized two of these transcripts, identifying a candidate WT gene in the process. I then performed a detailed analysis of this WT candidate gene, which maps to 7p. Future studies will shed more light on the role of these differentially expressed genes in WT. ^

Rescue of embryonic lethality in mdmx null mice by loss of p53 suggests a non-overlapping pathway with MDM2 to regulate p53

The tumor suppressor p53 is mutated in over 50% of human sporadic tumors originating from diverse tissues. p53 responds to DNA damage and cell stress by activating the transcription of a variety of target genes, the protein products of which then initiate either growth arrest or apoptosis. ^ A p53 target with a particularly intriguing function is the oncogene MDM2. MDM2 functions, in part, by binding to and inhibiting p53's activity. Overexpression of MDM2, by gene amplification, has been found in 30% of human sarcomas harboring a wild type p53, indicating that an increase in MDM2 levels is sufficient for p53 inactivation. Mice carrying a homozygous null allele for mdm2 exhibit an early embryonic lethality that is completely rescued in a p53-null background. These data indicate that MDM2's only critical function in early mouse embryogenesis is the negative regulation of p53. ^ The mdmx gene is the first additional member of the mdm2 gene family to be isolated. MDMX, like MDM2, contains a RING-finger domain, ATP binding domain and a p53 binding domain, which retains the ability to bind and inhibit p53 transactivation in vitro. However, mdmx does not appear to be transcriptionally regulated by p53. We have cloned and characterized the murine mdmx genomic locus from a mouse 129 genomic library. The mdmx gene contains 11 exons, spans approximately 37 Kb of DNA, and is located on mouse chromosome 1. The genomic organization of the mdmx gene is identical to that of mdm2 except at the 5′ end of the gene near the p53 responsive element. Northern expression analysis of mdmx transcripts during mouse embryogenesis and in adult tissues revealed constitutive and ubiquitous expression throughout adult tissues and embryonic development. To determine the in vivo function of MDMX, mice carrying a null allele of mdmx have been generated. Mdmx homozygous null mice are early embryonic lethal. Mdmx null mice do not develop beyond 9.5 dpc and can be discerned by gross dissection as early as 7.5 dpc. Utilizing TUNEL and BrdU assays on 7.5 dpc histological sections we have determined that the mutant embryos are dying due to increased levels of growth arrest, but not apoptosis. Surprisingly, Mdmx homozygous null mice are viable in a p53 null background, indicating that MDMX is also very important in the negative regulation of p53. ^

The role of Ski/Sno oncoproteins as negative regulators of the TGF-beta/SMAD signaling pathway in B -cell non -Hodgkin's lymphoma (NHL-B)

Non-Hodgkin's lymphomas are common tumors of the human immune system, primarily of B cell lineage (NHL-B). Negative growth regulation in the B cell lineage is mediated primarily through the TGF-β/SMAD signaling pathway that regulates a variety of tumor suppressor genes. Ski was originally identified as a transforming oncoprotein, whereas SnoN is an isoform of the Sno protein that shares a large region of homology with Ski. In this study, we show that Ski/SnoN are endogenously over-expressed both in patients' lymphoma cells and NHL-B cell lines. Exogenous TGF-β1 treatment induces down-regulation of Ski and SnoN oncoprotein expression in an NHL-B cell line, implying that Ski and SnoN modulate the TGF-β signaling pathway and are involved in cell growth regulation. Furthermore, we have developed an NHL-B cell line (DB) that has a null mutation in TGF-β receptor type II. In this mutant cell line, Ski/SnoN proteins are not down-regulated in response to TGF-β1 treatment, suggesting that downregulation of Ski and SnoN proteins in NHL-B require an intact functional TGF-β signaling pathway Resting normal B cells do not express Ski until activated by antigens and exogenous cytokines, whereas a low level of SnoN is also present in peripheral blood Go B cells. In contrast, autonomously growing NHL-B cells over-express Ski and SnoN, implying that Ski and SnoN are important cell cycle regulators. To further investigate a possible link between reduction of the Ski protein level and growth inhibition, Ski antisense oligodeoxynucleotides were transfected into NHL-B cells. The Ski protein level was found to decrease to less than 40%, resulting in restoring the effect of TGF-β and leading to cell growth inhibition and G1 cell cycle arrest. Co-immunoprecipitation experiments demonstrated that Ski associates with Smad4 in the nucleus, strongly suggesting that over-expression of the nuclear protein Ski and/or SnoN negatively regulates the TGF-β pathway, possibly by modulating Smad-mediated tumor suppressor gene expression. Together, in NHL-B, the TGF-β/SMAD growth inhibitory pathway is usually intact, but over-expression of the Ski and/or SnoN, which binds to Smad4, abrogates the negative regulatory effects of TGF-β/SMAD in lymphoma cell growth and potentiates the growth potential of neoplastic B cells. ^

Developmental functions of XKaiso, a transcriptional repressor associating with the p120 catenin in Xenopus laevis

The Armadillo family catenin proteins function in multiple capacities including cadherin-mediated cell-cell adhesion and nuclear signaling. The newest catenin, p120 catenin, differs from the classical catenins and binds to the membrane-proximal domain of cadherins. Recently, a novel transcription factor Kaiso was found to interact with p120 catenin, suggesting that p120 catenin also possesses a nuclear function. We isolated the Xenopus homolog of Kaiso, XKaiso, from a Xenopus stage 17 cDNA library. XKaiso contains an amino-terminal BTB/POZ domain and three carboxyl-terminal zinc fingers. The XKaiso transcript was present maternally and expressed throughout early embryonic development. XKaiso's spatial expression was defined via in situ hybridization and was found localized to the brain, eye, ear, branchial arches, and spinal cord. Co-immunoprecipitation of Xenopus p120 catenin and XKaiso demonstrated their mutual association, while related experiments employing differentially epitope-tagged XKaiso constructs suggest that XKaiso also self-associates. On the functional level, reporter assays employing a chimera of XKaiso fused to the GAL4 DNA binding domain indicated that XKaiso is a transcriptional repressor. To better understand the significance of the Kaiso-p120 catenin complex in vertebrate development, Kaiso knock-down experiments were undertaken, and the modulatory role of p120 catenin in Kaiso function examined during Xenopus development. Using morpholino antisense oligonucleotides to block translation of XKaiso, XKaiso was found to be essential for Xenopus gastrulation, being required for correct morphogenetic movements in early embryogenesis. Molecular marker analyses indicated that one target gene of the Wnt/β-catenin pathway, Siamois, is significantly increased in embryos depleted for XKaiso, while other dorsal, ventral, and mesodermal cell fate markers were unaltered. In addition, the non-canonical Wnt-11, known to participate in planar cell polarity/convergent extension processes, was significantly upregulated following depletion of XKaiso. Such increased Wnt-11 expression likely contributed to the XKaiso depletion phenotype because a dominant negative form of Wnt-11 or of the downstream effector Dishevelled partially rescued the observed gastrulation defects. These results show that XKaiso is essential for proper gastrulation movements, resulting at least in part from its modulation of non-canonical Wnt signaling. The significance of the XKaiso-p120 catenin interaction has yet to be determined, but appears to include a role in modulating genes promoting canonical and non-canonical Wnt signals. ^