Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

Extending in vitro conditioning in Aplysia to analyze operant and classical processes in the same preparation.

Operant and classical conditioning are major processes shaping behavioral responses in all animals. Although the understanding of the mechanisms of classical conditioning has expanded significantly, the understanding of the mechanisms of operant conditioning is more limited. Recent developments in Aplysia are helping to narrow the gap in the level of understanding between operant and classical conditioning, and have raised the possibility of studying the neuronal processes underlying the interaction of operant and classical components in a relatively complex learning task. In the present study, we describe a first step toward realizing this goal, by developing a single in vitro preparation in which both operant and classical conditioning can be studied concurrently. The new paradigm reproduced previously published results, even under more conservative and homogenous selection criteria and tonic stimulation regime. Moreover, the observed learning was resistant to delay, shortening, and signaling of reinforcement.

In vitro analog of classical conditioning of feeding behavior in aplysia.

The feeding behavior of Aplysia californica can be classically conditioned using tactile stimulation of the lips as a conditioned stimulus (CS) and food as an unconditioned stimulus (US). Moreover, several neural correlates of classical conditioning have been identified. The present study extended previous work by developing an in vitro analog of classical conditioning and by investigating pairing-specific changes in neuronal and synaptic properties. The preparation consisted of the isolated cerebral and buccal ganglia. Electrical stimulation of a lip nerve (AT4) and a branch of the esophageal nerve (En2) served as the CS and US, respectively. Three protocols were used: paired, unpaired, and US alone. Only the paired protocol produced a significant increase in CS-evoked fictive feeding. At the cellular level, classical conditioning enhanced the magnitude of the CS-evoked synaptic input to pattern-initiating neuron B31/32. In addition, paired training enhanced both the magnitude of the CS-evoked synaptic input and the CS-evoked spike activity in command-like neuron CBI-2. The in vitro analog of classical conditioning reproduced all of the cellular changes that previously were identified following behavioral conditioning and has led to the identification of several new learning-related neural changes. In addition, the pairing-specific enhancement of the CS response in CBI-2 indicates that some aspects of associative plasticity may occur at the level of the cerebral sensory neurons.

Genetics and epidemiology of gallbladder disease in New World native peoples.

Native peoples of the New World, including Amerindians and admixed Latin Americans such as Mexican-Americans, are highly susceptible to diseases of the gallbladder. These include cholesterol cholelithiasis (gallstones) and its complications, as well as cancer of the gallbladder. Although there is clearly some necessary dietary or other environmental risk factor involved, the pattern of disease prevalence is geographically associated with the distribution of genes of aboriginal Amerindian origin, and levels of risk generally correspond to the degree of Amerindian admixture. This pattern differs from that generally associated with Westernization, which suggests a gene-environment interaction, and that within an admixed population there is a subset whose risk is underestimated when admixture is ignored. The risk that an individual of a susceptible New World genotype will undergo a cholecystectomy by age 85 can approach 40% in Mexican-American females, and their risk of gallbladder cancer can reach several percent. These are heretofore unrecognized levels of risk, especially of the latter, because previous studies have not accounted for admixture or for the loss of at-risk individuals due to cholecystectomy. A genetic susceptibility may, thus, be as "carcinogenic" in New World peoples as any known major environmental exposure; yet, while the risk has a genetic basis, its expression as gallbladder cancer is so delayed as to lead only very rarely to multiply-affected families. Estimates in this paper are derived in part from two studies of Mexican-Americans in Starr County and Laredo, Texas.

Contribution of Ectodomain Mutations in Epidermal Growth Factor Receptor to Signaling in Glioblastoma Multiforme

CONTRIBUTION OF ECTODOMAIN MUTATIONS IN EPIDERMAL GROWTH FACTOR RECEPTOR TO SIGNALING IN GLIOBLASTOMA MULTIFORME Publication No._________ Marta Rojas, M.S. Supervisory Professor: Oliver Bögler, Ph.D. The Cancer Genome Atlas (TCGA) has conducted a comprehensive analysis of a large tumor cohort and has cataloged genetic alterations involving primary sequence variations and copy number aberrations of genes involved in key signaling pathways in glioblastoma (GBM). This dataset revealed missense ectodomain point mutations in epidermal growth factor receptor (EGFR), but the biological and clinical significance of these mutations is not well defined in the context of gliomas. In our study, we focused on understanding and defining the molecular mechanisms underlying the functions of EGFR ectodomain mutants. Using proteomic approaches to broadly analyze cell signaling, including antibody array and mass spectrometry-based methods, we found a differential spectrum of tyrosine phosphorylation across the EGFR ectodomain mutations that enabled us to stratify them into three main groups that correlate with either wild type EGFR (EGFR) or the long-studied mutant, EGFRvIII. Interestingly, one mutant shared characteristics of both groups suggesting a continuum of behaviors along which different mutants fall. Surprisingly, no substantial differences were seen in activation of classical downstream signaling pathways such as Akt and S6 pathways between these classes of mutants. Importantly, we demonstrated that ectodomain mutations lead to differential tumor growth capabilities in both in vitro (anchorage independent colony formation) and in vivo conditions (xenografts). Our data from the biological characterization allowed us to categorize the mutants into three main groups: the first group typified by EGFRvIII are mutations with a more aggressive phenotype including R108K and A289T; a second group characterized by a less aggressive phenotype exemplified by EGFR and the T263P mutation; and a third group which shared characteristics from both groups and is exemplified by the mutation A289D. In addition, we treated cells overexpressing the mutants with various agents employed in the clinic including temozolomide, cisplatin and tarceva. We found that cells overexpressing the mutants in general displayed resistance to the treatments. Our findings yield insights that help with the molecular characterization of these mutants. In addition, our results from the drug studies might be valuable in explaining differential responses to specific treatments in GBM patients.

Cerebellar cortex involvement in classical eyelid conditioning

A model for cerebellar involvement in motor learning was tested using classical eyelid conditioning in the rabbit. Briefly, we assume that modifications of the strength of granule cell synapses at Purkinje cells in the cerebellar cortex and mossy fiber (MF) synapses at cerebellar interpositus nuclei are responsible for the acquisition, adaptively-timed expression, and extinction of conditioned eyelid responses (CRs). A corollary of these assumptions is that the cerebellar cortex is necessary for acquisition and extinction. This model also suggests a mechanism whereby the cerebellar cortex can discriminate different times during a conditioned stimulus (CS) and thus mediate the learned timing of CRs. Therefore, experiments were done to determine the role of the cerebellar cortex in the timing, extinction, and acquisition of CRs. Lesions of the cerebellar cortex that included the anterior lobe disrupted the learned timing of CRs such that they occurred at extremely short latencies. Stimulation of MFs in the middle cerebellar peduncle as the CS could support differently timed CRs in the same animal. These data indicate that synaptic plasticity in the cerebellar cortex mediates the learned timing of CRs. These short-latency CRs which resulted from anterior lobe damage did not extinguish, while CRs in animals receiving lesions which did not include the anterior lobe extinguished normally. Preliminary data suggests that lesions of the anterior lobe which produce short-latency responses prevent the acquisition of CRs to a novel CS. These findings indicate that the anterior lobe of cerebellar cortex is necessary for eyelid conditioning. The involvement of the anterior lobe in eyelid conditioning has not been previously reported, however, the anterior lobe has generally been spared in lesion studies examining cerebellar cortex involvement in eyelid conditioning due to its relatively inaccessible location. The observation that the anterior lobe of the cerebellar cortex is not always required for the basic expression of CRs, but is necessary for response timing, extinction, and acquisition, is consistent with the hypothesis that eyelid conditioning can involve plasticity in both the cerebellar cortex and interpositus nucleus and that plasticity in the nucleus is controlled by Purkinje cell activity. ^

Behavioral and cellular analysis of classical conditioning of feeding behavior in Aplysia

In classical conditioning, an associative form of learning, animals learn to associate two stimuli. Cellular and molecular mechanisms for the induction and consolidation of associative learning and memory at the level of single cells and synaptic connections have been studied in both vertebrate and invertebrate animals. The majority of studies, however, relied on aversive stimuli to induce learning. This bias may limit the extent to which identified mechanisms generalize to other forms of associative learning and memory, such as appetitive forms. The goal of the present study was to develop a classical conditioning procedure for the marine mollusk Aplysia californica using appetitive reinforcement, and to analyze associative learning using behavioral and electrophysiological techniques. ^ Using tactile stimulation of the lips as the conditional stimulus (CS) and food as the unconditional stimulus (US) a training protocol was developed that reliably induced classical conditioning of feeding behavior. Memory persisted for at least 24 hours. The gross organization of reinforcement-mediating pathways was analyzed in additional behavioral experiments. Moreover, neurophysiological correlates of classical conditioning were identified and characterized in an in vitro preparation containing the circuitry for feeding behavior. In vitro stimulation of a nerve (AT4) that may mediate the CS during training, resulted in a greater number of buccal motor patterns (BMPs) in brains from conditioned animals, as compared to control animals. The majority of these BMPs were ingestion-like, consistent with the increased number of bites in response to the CS after classical conditioning. Moreover, classical conditioning correlated with increased excitatory synaptic input to BMP-initiating neuron B31/32, in response to stimulation of AT 4, as compared to controls. The expression of the correlates of classical conditioning identified in this study was specific to stimulation of AT 4, which is consistent the stimulus specificity that is characteristic for classical conditioning. ^ The identification of cellular correlates of classical conditioning documented here provides the basis for future, more detailed analyses of an appetitive form of associative learning and memory, that may extend the working knowledge of the cellular and molecular mechanisms for associative plasticity in general. ^

A novel myeloid leukemia suppressor locus at human chromosome 5q13.3

Molecular mechanisms that underlie preleukemic myelodysplasia (MDS) and acute myelogenous leukemia (AML) are poorly understood. In MDS or AML with a refractory clinical course, more than 30% of patients have acquired interstitial or complete deletions of chromosome 5. The 5q13.3 chromosomal segment is commonly lost as the result of 5q deletion. Reciprocal and unbalanced translocations of 5q13.3 can also occur as sole anomalies associated with refractory AML or MDS. This study addresses the hypothesis that a critical gene at 5q13.3 functions either as a classical tumor suppressor or as a chromosomal translocation partner and contributes to leukemogenesis. ^ Previous studies from our laboratory delineated a critical region of loss to a 2.5–3.0Mb interval at 5q13.3 between microsatellite markers D5S672 and GATA-P18104. The critical region of loss was later resolved to an interval of approximately 2Mb between the markers D5S672 and D5S2029. I, then generated a long range physical map of yeast artificial chromosomes (YACs) and developed novel sequence tagged sites (STS). To enhance the resolution of this map, bacterial artificial chromosomes (BACs) were used to construct a triply linked contig across a 1 Mb interval. These BACs were used as probes for fluorescent in situ hybridization (FISH) on an AML cell line to define the 5q13.3 critical region. A 200kb BAC, 484a9, spans the translocation breakpoint in this cell line. A novel gene, SSDP2 (single stranded DNA binding protein), is disrupted at the breakpoint because its first four exons are encoded within 140kb of BAC 484a9. This finding suggests that SSDP2 is the critical gene at 5q13.3. ^ In addition, I made an observation that deletions of chromosome 5q13 co-segregate with loss of the chromosome 17p. In some cases the deletions result from unbalanced translocations between 5q13 and 17p13. It was confirmed that the TP53 gene is deleted in patients with 17p loss, and the remaining allele harbors somatic mutation. Thus, the genetic basis for the aggressive clinical course in AML and MDS may be caused by functional cooperation between deletion or disruption of the 5q13.3 critical gene and inactivation of TP53. ^

Familial aggregation, candidate genes and genome scans: Analyzing the role of genetics in stroke

The genetic etiology of stroke likely reflects the influence of multiple loci with small effects, each modulating different pathophysiological processes. This research project utilized three analytical strategies to address the paucity of information related to the identification and characterization of genetic variation associated with stroke in the general population. ^ First, the general contribution of familial factors to stroke susceptibility was evaluated in a population-based sample of unrelated individuals. Increased risk of subclinical cerebral infarction was observed among individuals with a positive parental history of stroke. This association did not appear to be mediated by established stroke risk factors, specifically blood pressure levels or hypertension status. ^ The need to identify specific gene variation associated with stroke in the general population was addressed by evaluating seven candidate gene polymorphisms in a population-based sample of unrelated individuals. Three polymorphisms were significantly associated with increased subclinical cerebral infarction or incident clinical ischemic stroke risk. These relationships include the G-protein β3 subunit 825C/T polymorphism and clinical stroke in Whites, the lipoprotein lipase S/X447 polymorphism and subclinical and clinical stroke in men, and the angiotensin I-converting enzyme Ins/Del polymorphism and subclinical stroke in White men. These associations did not appear to be obfuscated by the stroke risk factors adjusted for in the analysis models specifically blood pressure levels or anti-hypertensive medication use. ^ The final research strategy considered, on a genome-wide scale, the idea that genetic variation may contribute to the occurrence of hypertension or stroke through a common etiologic pathway. Genomic regions were identified for which significant evidence of heterogeneity was observed among hypertensive sibpairs stratified by family history of stroke information. Regions identified on chromosome 15 in African Americans, and chromosome 13 in Whites and African Americans, suggest the presence of genes influencing hypertension and stroke susceptibility. ^ Insight into the role of genetics in stroke is useful for the potential early identification of individuals at increased risk for stroke and improved understanding of the etiology of the disease. The ultimate goal of these endeavors is to guide the development of therapeutic intervention and informed prevention to provide a lasting and positive impact on public health. ^

Cellular mechanisms of operant and classical conditioning

The ability to associate a predictive stimulus with a subsequent salient event (i.e., classical conditioning) and the ability to associate an expressed behavior with the consequences (i.e., operant conditioning) allow for a predictive understanding of a changing environment. Although they are operationally distinct, there has been considerable debate whether at some fundamental level classical and operant conditioning are mechanistically distinct or similar. Feeding behavior of Aplysia (i.e., biting) was chosen as the model system and was successfully conditioned with appetitive forms of both operant and classical conditioning. The neuronal circuitry responsible for feeding is well understood and is suitable for cellular analyses, thus providing for a mechanistic comparison between these two forms of associative learning. ^ Neuron B51 is part of the feeding circuitry of Aplysia and is critical for the expression of ingestive behaviors. B51 also is a locus of plasticity following both operant and classical conditioning. Both in vivo and in vitro operant conditioning increased the input resistance and the excitability of B51. No pairing-specific changes in the input resistance were observed following both in vivo and in vitro classical conditioning. However, classical conditioning decreased the excitability of B51. Thus, both operant and classical conditioning modified the threshold level for activation of neuron B51, but in opposite directions, revealing key differences in the cellular mechanisms underlying these two forms of associative learning. ^ Next, the cellular mechanisms underlying operant conditioning were investigated in more detail using a single-cell analogue. The single-cell analogue successfully recapitulated the previous in vivo and in vitro operant conditioning results by increasing the input resistance and the excitability of B51. Both PKA and PKC were necessary for operant conditioning. Dopamine appears to be the transmitter mediating the reinforcement signal in this form of conditioning. A D1 dopamine receptor antibody revealed that the D1receptor localizes to the axon hillock, which is also the region that gives the strongest response when iontophoresing dopamine. ^ The studies presented herein, thus, provide for a greater understanding of the mechanisms underlying both of these forms of associative learning and demonstrate that they likely operate through distinct cellular mechanisms. ^

Genotypic spectrum and genotype-phenotype correlation of trimethylaminuria

Trimethylaminuria (TMAU) or Fish odor syndrome is an autosomal recessive disease that is characterized by pungent body odor with subsequent psychosocial complications. There are limited studies of the sequence variants causing TMAU in the literature with most studies describing only one or two patients and lacking genotype-phenotype correlations. Also to date, there is no laboratory in the US or Europe that offers TMA genetic testing on a clinical basis. We have recently validated genetic testing in the University of Colorado DNA Diagnostic Laboratory. We have a database of a few dozen patients with a biochemical diagnosis of TMA at the University of Colorado at Denver Health Sciences Center (UCDHSC) which includes a few patients with the classical form of the disease. We have used the newly established clinical test in our institution to attempt to characterize the genotype (sequence variants including mutations and polymorphisms) of classical TMAU patients and to establish a genotype-phenotype (biochemical and clinical) association. The questionnaire results confirmed most of the previously reported epidemiological findings of TMAU and also indicated that TMAU patients use multiple intervention measures in attempt to control their symptoms with dietary control being most effective. Despite the complexity of intervention, most patients did not have any medical follow up and there was underutilization of specialist care. In a set of our patients, two deleterious mutations were identified in 4/12 patients including a novel T237P sequence variant, while the majority of our patients (8/12) did not reveal any mutations. Some of the latter were double heterozygous for the E158K and E308G polymorphisms which could explain a mild phenotype while others had only the E158K variant which raised the question of undetected mutations. These results indicate that further experiments are needed to further delineate the full mutational spectrum of the FMO3 gene. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^