Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

ROLE OF SOX9 IN UTERINE GLAND DEVELOPMENT AND DISEASE INITIATION

The female reproductive tract (FRT) develops midway through embryogenesis, and consists of oviducts, uterine horns, cervix and upper part of the vagina. The uterine horns are composed of an epithelial layer, luminal (LE) and glandular epithelium (GE), surrounded by a mesenchymal layer, the stroma and myometrium. Interestingly, in most mammals the GE forms after birth and it only becomes fully differentiated as the female reaches sexual maturity. Uterine glands (UG) are made up of GE and are present in all mammals. They secrete nutrients, cytokines and several other proteins, termed histotroph, that are necessary for embryo implantation and development. Experiments in ewes and mice have revealed that females who lack UGs are infertile mainly due to impaired implantation and early pregnancy loss, suggesting that UGs are essential for fertility. Fortunately for us, UGs develop after birth allowing us to peer into the genetic mechanism of tubulogenesis and branching morphogenesis; two processes that are disrupted in various adenocarcinomas (cancer derived from glands). We created 3D replicas of the epithelium lining the FRT using optical projection tomography and characterized UG development in mice using lineagetracing experiments. Our findings indicate that mouse UGs develop as simple tubular structures and later grow multiple secretory units that stem from the main duct. The main aim of this project was to study the role of SOX9 in the UGs. Preliminary studies revealed that Sox9 is mostly found in the nucleus of the GE. vii This observation led to the hypothesis that Sox9 plays a role in the formation and/or differentiation of the GE. To study the role of Sox9 in UGs differentiation, we conditionally knocked out and overexpressed Sox9 in both the LE and GE using the progesterone receptor (Pgr) promoter. Overexpressing Sox9 in the uterine epithelium, parts of the stroma, and myometrium led to formation of multiple cystic structures inside the endometrium. Histological analysis revealed that these structures appeared morphologically similar to structures present in histological tissue sections obtained from patients with endometrial polyps. We have accounted for the presence of simple and complex hyperplasia with atypia, metaplasia, thick-walled blood vessels, and stromal fibrosis; all “hallmarks” that indicate overexpressing Sox9 leads to development of a polyp-like morphology. Therefore, we can propose the use of Sox9-cOE mice to study development of endometrial cystic lesions and disease progression into hyperplastic lesions.

Retrospective assesment of the levonorgestrel intrauterine device to treat complex atypical hyperplasia and early endometrioid endometrial cancer

Objective: The primary objective of this project was to describe the efficacy of the Levonorgestrel Intrauterine Device (LIUD) for treatment of Complex Endometrial Cancer (CAH) and Grade 1 Endometrial Cancer (G1EEC) in terms of rate of Complete Response (CR) and Partial Response (PR) after 6 months of therapy. Finally, we assessed if any clinical or pathologic features were associated with response to the LIUD. ^ Methods: This study was a retrospective case series designed to report the response rate of patients with CAH or G1EEC treated with LIUD therapy. In addition, this study has a laboratory component to assess molecular predictors of response to LIUD therapy. Retrospective data already collected from patients diagnosed with CAH or EEC grade 1 and treated with LIUD therapy at MD Anderson Cancer Center (MDACC) were used for this study. Patients from all ethnic and race groups were included. A Complete Response (CR) was defined in patients diagnosed with CAH if pathologic report at 6 months demonstrated either no evidence of hyperplasia or no atypia in the setting of simple or complex hyperplasia. Partial Response (PR) was recorded if disease downgraded to only CAH from G1EEC. No Response (NR) was recorded if pathologic report demonstrates no change (Stable Disease, SD) or progression to cancer (Progressive Disease, PD). We calculated the proportion of patients with complete response to LIUD therapy with 95% confidence interval. We compared the response rates (CR/PR vs NR) by obesity status (Obese if BMI > 40 kg/m2 vs non-obese if BMI <= 40 kg/m2) as well as other clinical and pathologic factors, such as age, uterine size (median size), and presence of exogenous progesterone effect. ^ Results: There were 39 patients diagnosed with either CAH or G1EEC treated with the LIUD. Of 39 patients, 12 did not have pathological results of biopsy at 6months time period. Of 27 evaluable patients, 17 were diagnosed with CAH and 10 with G1EEC. Overall response rate (RR) was 78% (95% CI = 62-94%) at 6 months, 18 patients had CR (4 in G1EEC; 14 in CAH), 3 patients had PR (3 in G1EEC), 3 had SD (1 in CAH; 2 in G1EEC), 3 had PD (2 in CAH; 1 in G1EEC). After histology stratification, RR at 6 months was 82.35% (14/17; 95%CI = 67.4-97.3%) in CAH and 70% (7/10; 95% CI = 41-98.4%) in G1EEC. ^ There was no difference in response (R) and no response (NR) based on BMI (p=0.56). He observed a trend showing association between age with response (p=0.1). There was no association between uterine size and response to therapy (p=0.17). We recorded strong association between exogenous progesterone effect and response. ^ Conclusion: LIUD therapy for the treatment of CAH and G1EEC may be effective and safe. Presence of exogenous progesterone effect may predict the response to LIUD therapy at earlier time points. There is need of further studies with larger sample size to explore the relationship of response with other clinical and pathologic factors^

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

Diffusion Tensor Imaging in Alzheimer's disease

Attentional control and Information processing speed are central concepts in cognitive psychology and neuropsychology. Functional neuroimaging and neuropsychological assessment have depicted theoretical models considering attention as a complex and non-unitary process. One of its component processes, Attentional set-shifting ability, is commonly assessed using the Trail Making Test (TMT). Performance in the TMT decreases with increasing age in adults, Mild Cognitive Impairment (MCI) and Alzheimer’s Disease (AD). Besides, speed of information processing (SIP) seems to modulate attentional performance. While neural correlates of attentional control have been widely studied, there are few evidences about the neural substrates of SIP in these groups of patients. Different authors have suggested that it could be a property of cerebral white matter, thus, deterioration of the white matter tracts that connect brain regions related to set-shifting may underlie the age-related, MCI and AD decrease in performance. The aim of this study was to study the anatomical dissociation of attentional and speed mechanisms. Diffusion tensor imaging (DTI) provides a unique insight into the cellular integrity of the brain, offering an in vivo view into the microarchitecture of cerebral white matter. At the same time, the study of ageing, characterized by white matter decline, provides the opportunity to study the anatomical substrates speeded or slowed information processing. We hypothesized that FA values would be inversely correlated with time to completion on Parts A and B of the TMT, but not the derived scores B/A and B-A.

Preprocessing and analyzing genetic data with complex networks: An application to Obstructive Nephropathy

Many diseases have a genetic origin, and a great effort is being made to detect the genes that are responsible for their insurgence. One of the most promising techniques is the analysis of genetic information through the use of complex networks theory. Yet, a practical problem of this approach is its computational cost, which scales as the square of the number of features included in the initial dataset. In this paper, we propose the use of an iterative feature selection strategy to identify reduced subsets of relevant features, and show an application to the analysis of congenital Obstructive Nephropathy. Results demonstrate that, besides achieving a drastic reduction of the computational cost, the topologies of the obtained networks still hold all the relevant information, and are thus able to fully characterize the severity of the disease.

Knowledge Discovery in Spectral Data by Means of Complex Networks

In the last decade, complex networks have widely been applied to the study of many natural and man-made systems, and to the extraction of meaningful information from the interaction structures created by genes and proteins. Nevertheless, less attention has been devoted to metabonomics, due to the lack of a natural network representation of spectral data. Here we define a technique for reconstructing networks from spectral data sets, where nodes represent spectral bins, and pairs of them are connected when their intensities follow a pattern associated with a disease. The structural analysis of the resulting network can then be used to feed standard data-mining algorithms, for instance for the classification of new (unlabeled) subjects. Furthermore, we show how the structure of the network is resilient to the presence of external additive noise, and how it can be used to extract relevant knowledge about the development of the disease.

Network-based biomarkers in Alzheimer's disease: review and future directions

By 2050 it is estimated that the number of worldwide Alzheimer?s disease (AD) patients will quadruple from the current number of 36 million people. To date, no single test, prior to postmortem examination, can confirm that a person suffers from AD. Therefore, there is a strong need for accurate and sensitive tools for the early diagnoses of AD. The complex etiology and multiple pathogenesis of AD call for a system-level understanding of the currently available biomarkers and the study of new biomarkers via network-based modeling of heterogeneous data types. In this review, we summarize recent research on the study of AD as a connectivity syndrome. We argue that a network-based approach in biomarker discovery will provide key insights to fully understand the network degeneration hypothesis (disease starts in specific network areas and progressively spreads to connected areas of the initial loci-networks) with a potential impact for early diagnosis and disease-modifying treatments. We introduce a new framework for the quantitative study of biomarkers that can help shorten the transition between academic research and clinical diagnosis in AD.

Translocation and Functional Analysis of Pseudomonas savastanoi pv. savastanoi NCPPB 3335 Type III Secretion System Effectors Reveals Two Novel Effector Families of the Pseudomonas syringae Complex

Pseudomonas savastanoi pv. savastanoi NCPPB 3335 causes olive knot disease and is a model pathogen for exploring bacterial infection of woody hosts. The type III secretion system (T3SS) effector repertoire of this strain includes 31 effector candidates plus two novel candidates identified in this study which have not been reported to translocate into plant cells. In this work, we demonstrate the delivery of seven NCPPB 3335 effectors into Nicotiana tabacum leaves, including three proteins from two novel families of the P. syringae complex effector super-repertoire (HopBK and HopBL), one of which comprises two proteins (HopBL1 and HopBL2) that harbor a SUMO protease domain. When delivered by P. fluorescens heterologously expressing a P. syringae T3SS, all seven effectors were found to suppress the production of defense-associated reactive oxygen species. Moreover, six of these effectors, including the truncated versions of HopAA1 and HopAZ1 encoded by NCPPB 3335, suppressed callose deposition. The expression of HopAZ1 and HopBL1 by functionally effectorless P. syringae pv. tomato DC3000D28E inhibited the hypersensitive response in tobacco and, additionally, expression of HopBL2 by this strain significantly increased its competitiveness in N. benthamiana. DNA sequences encoding HopBL1 and HopBL2 were uniquely detected in a collection of 31 P. savastanoi pv. savastanoi strains and other P. syringae strains isolated from woody hosts, suggesting a relevant role of these two effectors in bacterial interactions with olive and other woody plants.

Identification of novel susceptibility loci for inflammatory bowel disease on chromosomes 1p, 3q, and 4q: Evidence for epistasis between 1p and IBD1

The idiopathic inflammatory bowel diseases, Crohn’s disease (CD) and ulcerative colitis (UC), are chronic, frequently disabling diseases of the intestines. Segregation analyses, twin concordance, and ethnic differences in familial risks have established that CD and UC are complex, non-Mendelian, related genetic disorders. We performed a genome-wide screen using 377 autosomal markers, on 297 CD, UC, or mixed relative pairs from 174 families, 37% Ashkenazim. We observed evidence for linkage at 3q for all families (multipoint logarithm of the odds score (MLod) = 2.29, P = 5.7 × 10−4), with greatest significance for non-Ashkenazim Caucasians (MLod = 3.39, P = 3.92 × 10−5), and at chromosome 1p (MLod = 2.65, P = 2.4 × 10−4) for all families. In a limited subset of mixed families (containing one member with CD and another with UC), evidence for linkage was observed on chromosome 4q (MLod = 2.76, P = 1.9 × 10−4), especially among Ashkenazim. There was confirmatory evidence for a CD locus, overlapping IBD1, in the pericentromeric region of chromosome 16 (MLod = 1.69, P = 2.6 × 10−3), particularly among Ashkenazim (MLod = 1.51, P = 7.8 × 10−3); however, positive MLod scores were observed over a very broad region of chromosome 16. Furthermore, evidence for epistasis between IBD1 and chromosome 1p was observed. Thirteen additional loci demonstrated nominal (MLod > 1.0, P < 0.016) evidence for linkage. This screen provides strong evidence that there are several major susceptibility loci contributing to the genetic risk for CD and UC.

Differential expression of major histocompatibility complex class II genes on murine macrophages associated with T cell cytokine profile and protective/suppressive effects

Protective/suppressive major histocompatibility complex (MHC) class II alleles have been identified in humans and mice where they exert a disease-protective and immunosuppressive effect. Various modes of action have been proposed, among them differential expression of MHC class II genes in different types of antigen-presenting cells impacting on the T helper type 1 (Th1)–Th2 balance. To test this possibility, the expression of H-2 molecules from the four haplotypes H-2b, H-2d, H-2k, and H-2q was determined on bone marrow-derived macrophages (BMDMs) and splenic B cells. The I-Ab and I-Ek molecules, both well characterized as protective/suppressive, are expressed at a high level on almost all CD11b+ BMDMs for 5–8 days, after which expression slowly declines. In contrast, I-Ad, I-Ak, and I-Aq expression is lower, peaks over a shorter period, and declines more rapidly. No differential expression could be detected on B cells. In addition, the differential MHC class II expression found on macrophages skews the cytokine response of T cells as shown by an in vitro restimulation assay with BMDMs as antigen-presenting cells. The results indicate that macrophages of the protective/suppressive haplotypes express MHC class II molecules at a high level and exert Th1 bias, whereas low-level expression favors a Th2 response. We suggest that the extent of expression of the class II gene gates the back signal from T cells and in this way controls the activity of macrophages. This effect mediated by polymorphic nonexon segments of MHC class II genes may play a role in determining disease susceptibility in humans and mice.

Identification of the von Hippel–Lindau tumor-suppressor protein as part of an active E3 ubiquitin ligase complex

Mutations of von Hippel–Lindau disease (VHL) tumor-suppressor gene product (pVHL) are found in patients with dominant inherited VHL syndrome and in the vast majority of sporadic clear cell renal carcinomas. The function of the pVHL protein has not been clarified. pVHL has been shown to form a complex with elongin B and elongin C (VBC) and with cullin (CUL)-2. In light of the structural analogy of VBC-CUL-2 to SKP1-CUL-1-F-box ubiquitin ligases, the ubiquitin ligase activity of VBC-CUL-2 was examined in this study. We show that VBC-CUL-2 exhibits ubiquitin ligase activity, and we identified UbcH5a, b, and c, but not CDC34, as the ubiquitin-conjugating enzymes of the VBC-CUL-2 ubiquitin ligase. The protein Rbx1/ROC1 enhances ligase activity of VBC-CUL-2 as it does in the SKP1-CUL-1-F-box protein ligase complex. We also found that pVHL associates with two proteins, p100 and p220, which migrate at a similar molecular weight as two major bands in the ubiquitination assay. Furthermore, naturally occurring pVHL missense mutations, including mutants capable of forming a complex with elongin B–elongin C-CUL-2, fail to associate with p100 and p220 and cannot exhibit the E3 ligase activity. These results suggest that pVHL might be the substrate recognition subunit of the VBC-CUL-2 E3 ligase. This is also, to our knowledge, the first example of a human tumor-suppressor protein being directly involved in the ubiquitin conjugation system which leads to the targeted degradation of substrate proteins.

Aggregation and disaggregation of senile plaques in Alzheimer disease

We quantitatively analyzed, using laser scanning confocal microscopy, the three-dimensional structure of individual senile plaques in Alzheimer disease. We carried out the quantitative analysis using statistical methods to gain insights about the processes that govern Aβ peptide deposition. Our results show that plaques are complex porous structures with characteristic pore sizes. We interpret plaque morphology in the context of a new dynamical model based on competing aggregation and disaggregation processes in kinetic steady-state equilibrium with an additional diffusion process allowing Aβ deposits to diffuse over the surface of plaques.

Functional analyses of troponin T mutations that cause hypertrophic cardiomyopathy: Insights into disease pathogenesis and troponin function

Mutations in a number of cardiac sarcomeric protein genes cause hypertrophic cardiomyopathy (HCM). Previous findings indicate that HCM-causing mutations associated with a truncated cardiac troponin T (TnT) and missense mutations in the β-myosin heavy chain share abnormalities in common, acting as dominant negative alleles that impair contractile performance. In contrast, Lin et al. [Lin, D., Bobkova, A., Homsher, E. & Tobacman, L. S. (1996) J. Clin. Invest. 97, 2842–2848] characterized a TnT point mutation (Ile79Asn) and concluded that it might lead to hypercontractility and, thus, potentially a different mechanism for HCM pathogenesis. In this study, three HCM-causing cardiac TnT mutations (Ile79Asn, Arg92Gln, and ΔGlu160) were studied in a myotube expression system. Functional studies of wild-type and mutant transfected myotubes revealed that all three mutants decreased the calcium sensitivity of force production and that the two missense mutations (Ile79Asn and Arg92Gln) increased the unloaded shortening velocity nearly 2-fold. The data demonstrate that TnT can alter the rate of myosin cross-bridge detachment, and thus the troponin complex plays a greater role in modulating muscle contractile performance than was recognized previously. Furthermore, these data suggest that these TnT mutations may cause disease via an increased energetic load on the heart. This would represent a second paradigm for HCM pathogenesis.

Transglutaminase-catalyzed inactivation of glyceraldehyde 3-phosphate dehydrogenase and α-ketoglutarate dehydrogenase complex by polyglutamine domains of pathological length

Several adult-onset neurodegenerative diseases are caused by genes with expanded CAG triplet repeats within their coding regions and extended polyglutamine (Qn) domains within the expressed proteins. Generally, in clinically affected individuals n ≥ 40. Glyceraldehyde 3-phosphate dehydrogenase binds tightly to four Qn disease proteins, but the significance of this interaction is unknown. We now report that purified glyceraldehyde 3-phosphate dehydrogenase is inactivated by tissue transglutaminase in the presence of glutathione S-transferase constructs containing a Qn domain of pathological length (n = 62 or 81). The dehydrogenase is less strongly inhibited by tissue transglutaminase in the presence of constructs containing shorter Qn domains (n = 0 or 10). Purified α-ketoglutarate dehydrogenase complex also is inactivated by tissue transglutaminase plus glutathione S-transferase constructs containing pathological-length Qn domains (n = 62 or 81). The results suggest that tissue transglutaminase-catalyzed covalent linkages involving the larger poly-Q domains may disrupt cerebral energy metabolism in CAG/Qn expansion diseases.