Browsing by Subject "Biology, Genetics"

A fundamental goal in the diverse field of evolutionary biology is reconstructing the historical processes that facilitated lineage diversification and the current geographic distribution of species diversity. Oceanic islands provide a view of evolutionary processes that may otherwise be obscured by the complex biogeographic histories of continental systems, and have thus provided evolutionary biology with some of its most lasting and significant theories. The Indian Ocean island of Madagascar is home to an extraordinarily diverse and endemic biota, and reconstructing the historical processes responsible for this diversity has consumed countless academic careers. While the flowering plant genus Coffea is but one lineage contributing to Madagascar's staggering floral diversity, it is representative of the common evolutionary theme of adaptive radiation and local endemism on the island. In this dissertation, I employ the genus Coffea as a model for understanding historical biogeographic processes in the Indian Ocean using methods of molecular phylogenetics and population genetics. In the molecular phylogenetic study of Coffea presented in chapter 2, I show that Madagascan Coffea diversity is likely the product of at least two independent colonization events from Africa, a result that contradicts current hypotheses for the single origin of this group.

Species of Coffea are known to exhibit self-incompatibly, which can have a dramatic affect on the geographic distribution of plant genetic diversity. In chapter 3, I identify the genetic mechanism of self-incompatibility in Coffea as homologous to the canonical eudicot S-RNase system. Baker's Rule suggests that self-incompatible lineages are very unlikely to colonize oceanic islands, and in chapter 4, I test this hypothesis by characterizing the strength of self-incompatibility and comparing S-RNase polymorphism in Coffea populations endemic to isolated Indian Ocean islands (Grande Comore and Mauritius) with that of Madagascan/African species. My findings suggest that while island populations show little evidence for genetic bottleneck in S-RNase allelic diversity, Mauritian endemic Coffea may have evolved a type of "leaky" self-incompatibility allowing self-fertilization at some unknown rate. Through the application of traditional phylogenetic methods and novel data from the self-incompatibly locus, my dissertation contributes a wealth of new information regarding the evolutionary and biogeographic history of Coffea in the Indian Ocean.

Pathologies of sexual development are common in humans and reflect the precarious processes of sex determination and sexual differentiation. The gonad forms as a bipotential organ, and recent results from the Capel lab revealed that it is initially balanced between testis and ovarian fates by opposing and antagonistic signaling networks. In XY embryos, this balance is disrupted by the transient expression of the Y-linked gene, Sry, which activates genes that promote the testis pathway and oppose the ovarian pathway. While the roles of a few genes have been defined by mutation, current evidence suggests that the interactions of many genes and signaling pathways are involved in the establishment of sexual fate. For example, most cases of disorders of sexual development (DSDs) are unexplained by mutations in known sex determination genes. In addition, recent microarray studies in the mouse revealed that nearly half the transcriptome is expressed in the gonad at the time of sex determination (Embryonic day 11.5, or E11.5), and as many as 1,500 genes are expressed in a sexually dimorphic pattern at this early stage. Thus the sexual fate decision in the developing gonad likely depends on a complex network of interacting factors that converge on a critical threshold.

To begin to elucidate the transcription network topology underlying sex determination, we exploited two inbred mouse strains with well-characterized differences in sex reversal. The common inbred strain C57BL/6J (B6) is uniquely sensitive to XY male-to-female sex reversal in response to a number of genetic perturbations, while other strains, including 129S1/SvImJ (129S1) and DBA/2J (D2) are resistant to sex reversal. We hypothesized that these strain differences in gonad phenotype likely result from underlying expression differences in the gonad at the critical timepoint of E11.5. Using microarrays, we identified significant, reproducible differences in the transcriptome of the E11.5 XY gonad between B6 and 129S1 indicating that the reported sensitivity of B6 to sex reversal is consistent with a higher expression of a female-like transcriptome in B6 XY gonads. Surprisingly, a well-characterized master regulator of the testis pathway, Sox9, was found to be upregulated in the sensitive B6 background, which may serve as a compensatory mechanism to counteract the female-leaning transcriptome and activate the testis pathway in wild type B6 XY gonads.

We extended our expression analysis to a large set of F2 XY gonads from B6 and 129S1 intercrosses. From each pair of gonads, we analyzed the expression of 56 sex-associated genes by nanoliter-scale quantitative RT-PCR (qRT-PCR). The expression levels of most genes were highly variable across the F2 population, yet strong correlations among genes emerged. We employed a First-Order Conditional Independence (FOCI) algorithm to estimate the F2 coexpression network. From this unbiased analysis of XY expression data, we uncovered two subnetworks consisting of primarily male and female genes. Furthermore, we predicted roles for genes of unknown function based on their connectivity and position within the network.

To identify the genes responsible for these strain expression differences, we genotyped each F2 embryo at 128 single nucleotide polymorphisms (SNPs) located evenly throughout the 19 autosomes and X chromosome. We then employed linkage analysis to detect autosomal regions that control the expression of one or more of the 56 genes in the F2 population. These regions are termed expression quantitative trait loci, or eQTLs. We identified eQTLs for many sex-related genes, including Sry and Sox9, the key regulators of male sex determination. In addition, we identified multiple prominent trans-band eQTLs that controlled the expression of many genes. My work represents the first eQTL analysis of a developing vertebrate organ, the mouse gonad. This systems-level approach revealed the complex transcription architecture underlying sex determination, and provides a mechanistic explanation for sensitivity to sex reversal seen in some individuals.

Gustation, also known as taste perception, is critical for the survival of most animal species. The fruit fly Drosophila melanogaster employs 68 different gustatory receptors (GRs) for the detection of sugars, bitter or toxic compounds, and pheromones. However, with a few notable exceptions, the functions of most GRs involved in feeding are unknown. Our research has focused on a cluster of highly-related Drosophila Grs, known as the Gr64 family, that have been shown to be critical for the perception of multiple sugars. Furthermore, we have demonstrated that another gene related to the Gr64 genes, Gr61a, is a sugar receptor that is narrowly tuned to a subset of pyranose sugars and may (along with the Gr64 genes) be indispensable for early fly development.

As a complementary approach to our behavioral analysis, we have examined the expression pattern of the Drosophila sugar receptors using knock-in driver alleles created by homologous recombination. As expected, most of these drivers have shown strong expression in various taste tissues. Intriguingly, some of these knock-in alleles also show expression in the maxillary palp and antenna, tissues previously thought to be involved only in olfaction. These expression patterns raise interesting questions about the true range of function of these chemosensory receptors and whether or not they might be involved in olfaction as well as gustation.

The Bayesian approach to model selection allows for uncertainty in both model specific parameters and in the models themselves. Much of the recent Bayesian model uncertainty literature has focused on defining these prior distributions in an objective manner, providing conditions under which Bayes factors lead to the correct model selection, particularly in the situation where the number of variables, p, increases with the sample size, n. This is certainly the case in our area of motivation; the biological application of genetic association studies involving single nucleotide polymorphisms. While the most common approach to this problem has been to apply a marginal test to all genetic markers, we employ analytical strategies that improve upon these marginal methods by modeling the outcome variable as a function of a multivariate genetic profile using Bayesian variable selection. In doing so, we perform variable selection on a large number of correlated covariates within studies involving modest sample sizes.

In particular, we present an efficient Bayesian model search strategy that searches over the space of genetic markers and their genetic parametrization. The resulting method for Multilevel Inference of SNP Associations MISA, allows computation of multilevel posterior probabilities and Bayes factors at the global, gene and SNP level. We use simulated data sets to characterize MISA's statistical power, and show that MISA has higher power to detect association than standard procedures. Using data from the North Carolina Ovarian Cancer Study (NCOCS), MISA identifies variants that were not identified by standard methods and have been externally 'validated' in independent studies.

In the context of Bayesian model uncertainty for problems involving a large number of correlated covariates we characterize commonly used prior distributions on the model space and investigate their implicit multiplicity correction properties first in the extreme case where the model includes an increasing number of redundant covariates and then under the case of full rank design matrices. We provide conditions on the asymptotic (in n and p) behavior of the model space prior

required to achieve consistent selection of the global hypothesis of at least one associated variable in the analysis using global posterior probabilities (i.e. under 0-1 loss). In particular, under the assumption that the null model is true, we show that the commonly used uniform prior on the model space leads to inconsistent selection of the global hypothesis via global posterior probabilities (the posterior probability of at least one association goes to 1) when the rank of the design matrix is finite. In the full rank case, we also show inconsistency when p goes to infinity faster than the square root of n. Alternatively, we show that any model space prior such that the global prior odds of association increases at a rate slower than the square root of n results in consistent selection of the global hypothesis in terms of posterior probabilities.

The physiology of the tumor microenvironment is characterized by lower oxygen (hypoxia), higher lactate, extracellular acidosis and glucose starvation. We examined the global, transcriptional cellular responses to each of these microenvironmental stresses in vitro, projected them onto clinical breast cancer patients' samples in vivo, and returned to perform further in vitro experiments to investigate the potential mechanisms involved in these stress responses. The reciprocal exchange of information was critical and advanced our understanding of the potential clinical relevance of cellular responses.

Our expression array result showed that lactic acidosis induces a strong response, distinct from that of hypoxia in human mammalian epithelial cells (HMECs), indicating lactic acidosis is not only a by-product of hypoxia but has a unique role as a stimulant to cells in the tumor microenvironment. Cellular responses to lactosis and acidosis further demonstrated that acidosis was the main driving force in the lactic acidosis response. These responding gene signatures were then statistically projected into clinical breast cancer patients' expression data sets. The hypoxia response, as reported previously, was associated with bad prognosis, where as the lactic acidosis and acidosis responses, were associated with good prognosis. Additionally, the acidosis response could be used to separate breast tumors with high versus low aggressiveness based on its inversed correlation with metastatic character. We further discovered that lactic acidosis, in contrast to hypoxia, abolished Akt signaling. Moreover, it downregulated glycolysis and shifted energy utilization towards aerobic respiration.

We continued to examine the cellular response to lactic acidosis temporally in MCF7 cells, a breast cancer cell line. The lactic acidosis response of MCF7 cells also showed the prognostic result of better clinical outcomes in datasets of breast cancer patients. The lactic acidosis responses of HMEC and MCF cells were highly correlated. Strikingly in MCF7 cells, lactic acidosis and glucose deprivation actually induced similar transcriptional profiles, with only a few genes being oppositely regulated. Furthermore, lactic acidosis, similar to glucose starvation, induced AMPK signaling and abolished mTOR. However, lactic acidosis and glucose deprivation induced opposite glucose uptake phenotypes. Lactic acidosis significantly repressed glucose uptake whereas glucose deprivation significantly induced it. Among the genes differentially regulated by these two stresses, thioredoxin-interacting protein (TXNIP) was among the most different. The negative regulatory role of TXNIP on glucose uptake has been demonstrated previously. In the cancer research field, TXNIP is recognized as a tumor suppressor gene. We observed that lactic acidosis induced TXNIP strongly and most importantly, TXNIP played a critical role in regulating glucose uptake in cells under lactic acidosis. Furthermore, MondoA, the transcription factor and glucose sensor previously reported to regulate TXNIP induction upon glucose exposure, was also responsible for regulating TXNIP under lactic acidosis. We demonstrated that TXNIP not only plays an important role in the lactic acidosis response but also has strong prognostic power to separate breast cancer patients based on survival.

Disease caused by the Gram-negative Haemophilus cryptic genospecies begins with colonization of the maternal genital or neonatal respiratory tract. The primary goal of this work was to identify and characterize the molecular determinant(s) of Haemophilus cryptic genospecies adherence as a means to better understand the specific adaptation of this species to the urogenital tract and neonatal respiratory tract. Using transposon mutagenesis of prototype strain 1595, we identified a locus that is essential for Haemophilus cryptic genospecies adherence to a variety of epithelial cell lines of both genital and respiratory origin. This locus encodes a protein called Cha that shares homology with trimeric autotransporters. Trimeric autotransporters are composed of an N-terminal signal peptide, an internal passenger domain that harbors adhesive activity, and a short C-terminal membrane anchor domain and are classically characterized by head-stalk-anchor domain architecture. By generating chimeric proteins, we demonstrated that the C-terminus of Cha trimerizes in the bacterial outer membrane and is capable presenting a heterologous passenger domain (Hia) in a functional form, thus confirming that Cha is a trimeric autotransporter. Southern analysis revealed that cha is unique to the Haemophilus cryptic genospecies and is ubiquitous among these strains.

Similar to a number of trimeric autotransporters, the passenger domain of Cha contains scattered clusters of YadA-like head domains associated with head-to-stalk neck adaptor motifs, predicted coiled-coil stalks and a series of identical tandem coding repeats which are not required for adherence. By evaluating the adherence capacity of H. influenzae expressing Cha deletion derivatives, we established that the N-terminal 473 residues of Cha harbor the binding domains responsible for Cha-mediated adherence to epithelial cells. In additional studies, we demonstrated that this same N-terminal region mediates bacterial aggregation through inter-bacterial Cha-Cha binding.

Further analysis revealed that variable Cha-mediated adherence is linked to spontaneous changes in the number of identical tandem repeats predicted to comprise a coiled-coil stalk domain. Variation in repeat copy number has a direct effect on Cha adhesive and aggregative activity, independent of an impact on transcription of the cha locus or surface localization of Cha protein. Moreover, length of Cha surface fibers correlates with repeat copy number expansion. We propose two hypotheses to explain how repeat expansion inhibits bacterial aggregation and host cell binding: 1) Variation in the number of 28-amino acid repeats may influence the conformation of Cha, thus changing the surface accessibility of the Cha binding pocket. 2) Repeat expansion results in the formation of long, flexible Cha fibers on the bacterial cell surface that may have a greater propensity to interact with neighboring Cha trimers at the N-terminus, thereby precluding adherence to other bacteria or host epithelial cells.

In additional studies screening adherent cryptic genospecies isolates for expression of Cha protein, we identified an additional, antigenically-divergent Cha variant that we refer to as Cha2. Amino acid sequence and domain comparison of Cha2 with Cha (now Cha1) revealed that the structurally undefined N-terminal sequences (encompassing the Cha1 adhesive and aggregative domain) are strikingly divergent. Inspite of this, Cha2 mediates efficient adherence to human epithelial cells, similar to Cha1.

Identification of Cha offers insight into the apparent tissue tropism associated with the Haemophilus cryptic genospecies. We speculate that the unique regulation of Cha adhesive activity enhances the adaptive capability of this pathogenic organism in the human host.

Copper is an essential trace element required by all aerobic organisms as a co-factor for enzymes involved in normal growth, development and physiology. Ctr1 proteins are members of a highly conserved family of copper importers responsible for copper uptake across the plasma membrane. Mice lacking Ctr1 die during embryogenesis from widespread developmental defects, demonstrating the need for adequate copper acquisition in the development of metazoan organisms via as yet uncharacterized mechanisms. The early lethality of the Ctr1 knockout mouse has made it difficult to study the functions of copper and Ctr1 proteins in metazoan development and physiology. Drosophila melanogaster, a genetically tractable system expresses three Ctr1 genes, Ctr1A, Ctr1B and Ctr1C, and may help to further understand the roles of copper and Ctr1 proteins in metazoan development and physiology. Described here is the characterization of Drosophila Ctr1A.

Localization studies using an affinity purified anti-Ctr1A peptide antibody show Ctr1A is predominantly expressed at the plasma membrane in whole embryos and in larval tissues. Ctr1A is an essential gene in Drosophila as loss-of-function mutants, generated by imprecise p-element excision arrest at early larval stages of development. Inductively coupled plasma mass spectroscopy (ICP-MS) demonstrated that whole body copper levels are reduced in Ctr1A mutants and consequently, a number of copper-dependent enzyme deficiencies were detected by in vitro enzyme and cell biological assays. Ctr1A maternal and zygotic mutants have a more severe developmental phenotype and also showed reductions in heart rate, which could be partially rescued by dietary copper supplementation. Heart-specific Ctr1A knockdown flies were subsequently examined for heart rate defects using optical coherence tomography (OCT) and while they did have reduced heart rate measurements, heart contractility was compromised. While investigating tissue-specific requirements for Ctr1A in the development of Drosophila, a genetic interaction between Ctr1A and Ras was observed. Genetic experiments in Drosophila and cell culture experiments in both Drosophila and mammalian cell lines demonstrate a conserved role for Ctr1 proteins and copper as positive modulators of Ras/MAPK pathway signaling. Immunoblot analysis shows that signal transduction is intact until the point at which MEK1/2 phosphorylates ERK1/2. MEK2 protein levels are reduced in copper deficient cells, while MEK1 is able to bind copper-chelated beads, suggesting that these two proteins may be copper-binding proteins. In summary, this work demonstrates that Ctr1A is an essential gene in Drosophila and through characterization studies of Ctr1A, has uncovered conserved roles for Ctr1 proteins and copper in physiological processes and in an important signaling pathway that controls a number of fundamental biological processes.

Leukotrienes are arachidonic acid derivatives long known for their inflammatory properties and their involvement with a number of human diseases, most notably asthma. Recently, leukotriene-based inflammation has also been implicated in atherosclerosis: ALOX5AP and LTA4H, two genes in the leukotriene biosynthesis pathway, have been associated with various cardiovascular disease (CVD) phenotypes. To assess the role of the leukotriene pathway in CVD pathogenesis, we performed genetic association studies of ALOX5AP and LTA4H in a non-familial data set of early onset coronary artery disease. Our results support a modest role for the leukotriene pathway in atherosclerosis pathogenesis, reveal important genomic interactions within the pathway, and suggest the importance of using pathway-based modeling for evaluating the genomics of atherosclerosis susceptibility. Motivated by this need, we investigated the statistical properties of a class of matrix-based statistics to assess epistasis. We simulated multiple two-variant disease models with haplotypes to gain an understanding of pathway interactions in terms of correlation patterns. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequlibrium (LD) patterns with other haplotype markers. The simulated models can be summarized into three categories: 1. No epistasis in the presence of marginal effects and LD; 2. Epistasis in the presence of LD and no marginal effects; and 3. Epistasis in the presence marginal effects and LD. We then assessed previously introduced single-gene methods that compare whole matrices of Single Nucleotide Polymorphism (SNP) LD between two samples. These methods include comparing two sets of principal components, a sum-of-squared-differences comparing pairwise LD, and a contrast test that controls for background LD. We also considered a partial least-square (PLS) approach for modeling gene-gene interactions. Our results indicate that these measures can be used to assess epistasis as well as marginal effects under certain disease models. Understanding and quantifying whole-gene variation and association to disease using multiple SNPs remains a difficult task. Providing a single statistical measure per gene will facilitate combining multiple types of genomic data at a gene-level and will serve as an alternative approach to assess epistasis in genome-wide association studies. The matrix-based measures can also be used in pathway ascertainment tools that require scores on a gene-level.

Telomeres are the ends of chromosomes which are composed of repetitive DNA sequence and telomere associated proteins. In C. elegans, the protein F39H2.5 was found to associate with the telomere, regulating both telomere length and genomic integrity. F39H2.5 is a member of the β-CASP family of proteins that are known to possess nuclease activity on DNA substrates. I thus sought to address whether any of the human β-CASP family proteins associated with telomeres. Here I show that hSnm1B localized to the telomere indirectly, via interaction with the double-stranded telomere binding protein TRF2. The terminal 37 amino acids of hSnm1B are necessary and sufficient for binding TRF2, and moreover that binding to TRF2 stabilized hSnm1B protein by preventing ubiquitination. In the absence of exogenous TRF2 this domain acted as a degron, promoting protein instability. I thus termed the domain the Protection And INstability (PAIN) domain. I hypothesize that TRF2 binding ensures that hSnm1B will only accumulate at telomeres by preventing the degradation of hSnm1B. However, hSnm1B stability appears to be further regulated, as telomere specific DNA damage stabilized hSnm1B independent of the PAIN domain. Thus, it appears that the telomere associated protein, hSnm1B, is regulated by protein stability in a manner that is both dependent and independent of the PAIN domain.

Plants and their pathogens coevolve, with pathogen infection and host resistance acting in evolutionary antagonism of each other. Plant-pathogen coevolution has been shown to effect genetic divergence between populations and species, resulting in localized or specialized interactions between hosts and pathogens. Because most of the studies to date investigating plant-pathogen coevolution have been carried out in managed systems and have focused on pairwise interactions, we know little about three aspects of plant pathosystems in natural settings: 1) the role in nature of the gene-for-gene paradigm for genetic determination of resistance; 2) the relationship of host community diversity and structure, and host-pathogen interaction structure, to the coevolutionary process; and 3) the factors which underlie and drive local adaptation and specialization of interactions.

This dissertation constitutes the results of research in which I have begun addressing these aspects in a natural plant-fungus pathosystem comprising three Ipomoea host species and a single rust pathogen, Coleosporium ipomoeae. I have expanded previous characterization of the genetics of plant resistance in one constituent host species in the system by genetic crosses to characterize the basis of resistance in two additional species, finding support for the expectation that the gene-for-gene paradigm of interaction is important in natural systems. I conducted a cross-inoculation experiment designed to assess host and pathogen variation in infectivity and resistance, to investigate patterns of community interaction structure, and the role that antagonistic coevolution may play in structuring the communities which compose pathosystems. In these experiments I found that the coevolutionary interaction in this system leads to genetic divergence and the substantial amount of host and pathogen variation I discovered, but that it tends to preserve one pattern of community interaction structure across communities. I expanded my cross-inoculation experimental design to facilitate analysis of quantitative aspects of pathogenesis by measuring the intensity of infections, to test existing hypotheses concerning local adaptation and specialization in pathosystems. In this analysis I found strong host local adaptation and pathogen local maladaptation for the qualitative interaction trait of infectivity, and I found weak host local maladaptation and pathogen local adaptation for the quantitative interaction trait of aggressiveness. I also found host specialization among pathogens, and specialized resistance among hosts, to be common in this system. In light of these results, I hypothesize that the geographic scale of host-pathogen coevolution in this system is that of the local community, and that differences between host species result in persistent but incomplete host specialization in pathogen races.

The cell cycle is a series of ordered events that culminates in a single cell dividing into two daughter cells. These events must be properly coordinated to ensure the faithful passage of genetic material. How cell cycle events are carried out accurately remains a fundamental question in cell biology. In this dissertation, I investigate mechanisms orchestrating cell-cycle events in the yeast, Saccharomyces cerevisiae.

Cyclin dependent kinase (CDK) activity is thought to both form the fundamental cell-cycle oscillator and act as an effector of that oscillator, regulating cell-cycle events. By measuring transcript dynamics over time in cells lacking all CDK activity, I show that transcriptional oscillations are not dependent on CDK activity. This data indicates that CDKs do not form the underlying cell-cycle oscillator. I propose a model in which a transcription factor network rather than CDK activity forms the cell-cycle oscillator. In this model, CDKs are activated by the periodic transcription of cyclin genes and feedback on the network increasing the robustness of network oscillations in addition to regulating cell-cycle events.

I also investigate CDK-dependent and -independent mechanism regulating the duplication of the yeast centrosome, the spindle pole body (SPB). It is critical for the formation of a bipolar spindle in mitosis that the SPB duplicates once and only once per cell cycle. Through a combination of genetic and microscopic techniques I show that three distinct mechanisms regulate SPB duplication, ensuring its restriction to once per cell cycle.

Together, the data presented in this dissertation support a model in which CDKs, periodic transcription, and a TF-network oscillator are all important cell-cycle regulatory mechanisms that collaborate to regulate the intricate collection of events that constitute the cell cycle.

Heterochromatin, or condensed chromatin, is a transcriptionally repressive form of chromatin that occurs in many eukaryotic organisms. At its natural locations, heterochromatin is thought to play important roles in genome organization as well as gene expression. Just as important is the restriction of this repressive form of chromatin to appropriate regions of the genome. In the budding yeast Saccaromyces cerevisiae, domains of condensed, transcriptionally silenced chromatin are found at telomeres and at the silent-mating type cassettes, HML and HMR. At these locations, a complex of Silent Information Regulator (SIR) proteins gets recruited to DNA through discrete silencer elements. Once recruited, the Sir protein complex then spreads along chromosomes in a step-wise manner. This process results in the silencing of gene expression. It is unclear whether silenced chromatin is established in the same manner at different genomic locations. Understanding how silenced chromatin is formed is important for determining how these chromatin structures are regulated.

To better understand how silenced chromatin is established in different genomic contexts, I used chromatin immuoprecipitation to follow the rate of silenced chromatin formation at different locations. The rates of Sir protein assembly were compared at two locations, telomere VI-R and HMR. I discovered that the silencers at these two locations were equally proficient at recruiting Sir proteins. However, the rate of Sir protein assembly onto nucleosomes was far more rapid at HMR than at the telomere VI-R. Furthermore, the rate of Sir protein assembly was more rapid on one side of the HMR-E silencer at HMR than the other. Moreover, insertion of the HMR-E silencer adjacent to the telomere VI-R significantly improved the rate of Sir protein assembly onto nucleosomes. Additionally, observations that the association of Sir protein occurs simultaneously across several kilobases at HMR and that silencing at HMR is insensitive to co-expression of wild-type and catalytically inactive Sir2 proteins suggest that HMR-E enables the assembly of silenced chromatin in a non-linear fashion. These results suggest that HMR-E functions to both recruit Sir proteins and promote their assembly across several kilobases.

In addition to the HMR-E silencer, HMR is also characterized by the presence of a second auxiliary HMR-I silencer and a tRNA gene that functions as a boundary element to restrict the spread of silenced chromatin. I used chromatin immunoprecipitation to determine how each of these regulatory elements contribute to the steady-state levels of Sir protein association with chromatin. Consistent with a role for HMR-E beyond recruitment, I discovered that the HMR-E silencer alone promoted higher levels of Sir proteins on nucleosomes compared to the telomere VI-R. The levels of Sir protein association with HMR were further elevated by the HMR-I silencer, even though this silencer does not recruit Sir proteins on its own and does not contribute to any of the known functions of silenced chromatin at HMR. Additionally, although the tRNA gene did block the spread Sir proteins, I discovered that the capacity for Sir proteins to spread beyond a few kilobases was severely limited even in the absence of the boundary.

The results of this thesis work provide new insights into the mechanisms of silenced chromatin establishment and regulation in budding yeast. I show here that the capacity of Sir proteins to assemble onto nucleosomes is inherently limited. Additionally, silencers vary in their ability to promote this assembly. I conclude that the silencer is a key factor in determining the relative size, efficiency, and location of silenced chromatin domains in the cell.

Evolution can be studied at many levels, from phenotypic to molecular, and from a variety of disciplines. An integrative approach can help provide a more complete understanding of the complexities of evolutionary change. This dissertation examines the ecology, genetics, and molecular mechanisms of the evolution of floral anthocyanin pigmentation in four species of Mimulus native to central Chile. Anthocyanins, which create red and purple colors in many plants, are a valuable model for studying evolutionary processes. They are ecologically important and highly variable both within and between species, and the underlying biosynthetic pathway is well characterized. The focus of this dissertation is dramatic diversification in anthocyanin coloration, in four taxa that are closely related to the genomic model system M. guttatus. I posed three primary questions: (1) Is floral diversification associated with pollinator divergence? (2) What is the genetic basis of the floral diversification? (3) What is the molecular mechanism of the increased production of anthocyanin pigment? The first question was addressed by evaluating patterns of pollinator visitation in natural populations of all four study taxa. The second question was explored using segregation analysis for a series of inter- and intraspecific crosses. One trait, increased petal anthocyanins in M. cupreus, was further dissected at the molecular level, using candidate gene testing and quantitative gene expression analysis. Pollinator studies showed little effect of flower color on pollinator behavior, implying that pollinator preference probably did not drive pigment evolution in this group. However, segregation analyses revealed that petal anthocyanin pigmentation has evolved three times independently in the study taxa, suggesting an adaptive origin. In addition to pollinator attraction, anthocyanins and their biochemical precursors protect against a variety of environmental stressors, and selection may have acted on these additional functions. Molecular analysis of petal anthocyanins in M. cupreus revealed that this single-locus trait maps to a transcription factor, McAn1, which is differentially expressed in high- versus low-pigmented flowers. Expression of the anthocyanin structural genes is tightly correlated with McAn1 expression. The results suggest that M. cupreus pigmentation evolved by a mutation cis to McAn1 that alters the intensity of anthocyanin biosynthesis.

Sexual reproduction in fungi is governed by a specialized genomic region called the mating-type locus (MAT). The ascomycetes, the largest phylum of fungi, primarily possess a bipolar mating system while the basidiomycetes, the second largest group, are mostly tetrapolar. The human fungal pathogen and basidiomycetous yeast Cryptococcus neoformans has evolved a bipolar mating system that encodes homeodomain (HD) and pheromone/receptor (P/R) genes. The MAT locus of C. neoformans is unusually large, spans greater than 100 kb, and encodes more than 20 genes. To understand how the pathogenic Cryptococcus species complex evolved this unique bipolar mating system, we investigated the evolution of MAT in closely and distantly related species and discovered an extant sexual cycle in Cryptococcus amylolentus.

Phylogenetic analysis using a six-gene multi-locus sequencing (MLS) approach identified the most closely related species to the pathogenic Cryptococcus species complex that are currently known. The two non-pathogenic sibling species, Tsuchiyaea wingfieldii and Cryptococcus amylolentus, and the more distantly related species Filobasidiella depauperata define the Filobasidiella clade. We also resolved the phylogeny of the species located in the sister clade, Kwoniella. A comprehensive tree dendrogram revealed that the 15 Tremellales species examined suggests a common saprobic ancestor. Moreover, the pathogenic Cryptococcus species have a saprobic origin but later emerged as pathogens. We further characterized the mating-type locus for T. wingfieldii and C. amylolentus by cloning and sequencing two unlinked genomic loci encoding the HD and P/R genes. Interestingly, linked and likely divergently transcribed homologs for SXI1 and SXI2 are present in T. wingfieldii and C. amylolentus, while the P/R alleles contain many genes also found in the MAT locus of the pathogenic Cryptococcus species. Also, hypothetical genes present in C. neoformans MAT are also MAT-linked in both species and indicate a possible translocation event between chromosomes 4 and 5 of C. neoformans. Our analysis of MAT in the sibling species indicates that T. wingfieldii is likely tetrapolar, and the C. amylolentus sequence comparison of the dimorphic SXI1 and SXI2 region and the pheromone receptor, STE3, suggests that C. amylolentus is also tetrapolar. The examination of MAT in these sibling species confirms the model for MAT evolution previously proposed in which this structure in C. neoformans and C. gattii evolved from an ancestral tetrapolar mating system. Moreover, the organization of MAT in these sibling species mirrors key aspects of the proposed intermediates in the evolution of MAT in the pathogenic Cryptococcus species, and for sex chromosomes in plants, animals, and alga in general.

We discovered an extant sexual cycle for C. amylolentus, a species previously thought to be asexual. Matings between two strains of opposite mating-types produce dikaryotic hyphae with fused clamp connections and uni- and bi-nucleate basidiospores. Genotyping of basidiospores using markers linked and unlinked to MAT revealed that genetic exchange (recombination) occurs during the sexual cycle of C. amylolentus, and it is likely that either aneuploids are generated during sex or more than one meiosis event occurs within each basidium. This is in contrast to C. neoformans, where only one meiotic event per basidium has been observed. Uniparental mitochondrial inheritance has also been observed in C. amylolentus progeny; similar to the pathogenic Cryptococcus species, mtDNA is inherited from the C. amylolentus MATa parent. Analysis of sex in C. amylolentus has provided insight into the mechanisms that phylogenetically related fungi employ in orchestrating sexual reproduction.

We also extended our analysis to include the distantly related tetrapolar basidiomycete Tremella mesenterica. We completed comparisons of MAT-specific genes between five strains of T. mesenterica and identified the regions that define its mating-type system. The HD locus is limited to the SXI1- and SXI2-like genes while the P/R locus is defined by STE3, STE12, STE20, and the pheromone gene, tremerogen a-13. Interestingly, many of the genes associated with the MAT locus of the pathogenic Cryptococcus species flank the HD and P/R locus and are not incorporated in MAT in T. mesenterica. The MAT region includes transposons and C. neoformans hypothetical genes also present in T. wingfieldii and C. amylolentus. The mating-type system in T. mesenterica reflects an ancestral intermediate in the evolution of the MAT locus in the pathogenic Cryptococcus species. In conclusion, this study provides an in-depth analysis on the structure, function, and evolution of an unusual mating-type locus with broader implications for the transitions in modes of sexual reproduction in fungi that impact gene flow in populations.

Hydrogen peroxide is used by animals and plants to deter the growth of microbial invaders by inflicting DNA lesions, protein oxidation and lipid membrane modifications. Pathogens protect themselves with enzymes and scavenging proteins. This study investigated population genetic, biochemical and genetic aspects of peroxide survival in Saccharomyces cerevisiae to address its importance for yeast biology and fungal pathogenicity.

Population genetic analyses of DNA sequences from five loci from 103 strains encompassing the known ecological spectrum of S. cerevisiae show that it is a recombining species with three divergent subgroups, which are associated with soil, fruit, and vineyards. Clinical isolates cluster with fruit isolates but are significantly more resistant to peroxide. Clinical isolates are genetically diverse, indicating multiple origins of the pathogenic lifestyle and eliminating the possibility that peroxide resistance is due to shared ancestry rather than it's importance for than its importance in colonizing the host.

Biochemical aspects of peroxide survival were studied in a resistant (high-survival) clinical isolate, a sensitive (low-survival) laboratory strain and their hybrid. Catalase activity and expression levels are indistinguishable among strains. Co-culture assays and growth curve records indicate that a secreted factor improves survival of the laboratory strain and that the phenotypic difference is most pronounced during exponential growth, excluding mechanisms of the General Stress Response effective during stationary phase. Semi-quantitative expression profiles of stress response candidate genes do not differ, suggesting a novel resistance mechanism.

To elucidate the genetic basis of peroxide survival, the hybrid was sporulated and 200 F1 segregants phenotyped and genotyped for oxidative stress candidate genes. Peroxide survival is a dominant quantitative trait and not linked to catalase, peroxidase or superoxide dismutase genes. 1,246 backcross segregants were phenotyped and 93 segregants selectively genotyped using microarrays. A 14-gene locus on chromosome XVI displayed marker-trait association. One gene, RDS2, encodes a zinc cluster protein acting as a regulator of drug sensitivity and contains a non-synonymous polymorphism whose exchange between the parental strains results a 15% decrease in survival in the clinical strain.

This work establishes a novel function for RDS2 in oxidative stress response and demonstrates the effect a quantitative trait nucleotide has on a clinically relevant phenotype.

Neurons, with their enormous membrane contents, depend heavily on regulated membrane trafficking processes to maintain their morphology and function. The phosphatidylinositol 3-kinase class III, or PIK3C3, plays a critical role in various membrane trafficking processes including both the endocytic and autophagic pathways. The functions of PIK3C3 in the nervous system in vivo are un-characterized. We reasoned that studying PIK3C3 in neurons would provide us an entry point into understanding the regulations and functions of the neuronal membrane trafficking processes and their roles in neuronal morphogenesis and homeostasis.

We generated a conditional allele of Pik3c3 and first deleted it specifically in the peripheral sensory neurons. Mutant large-diameter myelinated sensory neurons accumulated numerous enlarged vacuoles and ubiquitin-positive aggregates and underwent rapid degeneration. By contrast, Pik3c3-deficient small-diameter unmyelinated neurons accumulated excessive numbers of lysosome-like organelles and degenerated slower than large-diameter neurons. These differential degenerative phenotypes are unlikely caused by a disruption of the autophagy pathway, because inhibiting autophagy alone by conditional deletion of Atg7 results in a completely distinct subcellular phenotypes and very slow degenerations of all sensory neurons. More surprisingly, a noncanonical PIK3C3-independent LC3-positive autophagosome formation pathway was activated in Pik3c3-deficient small-diameter neurons. This work uncovered unexpected differences of the endo-lysosomal systems in different types of neurons and discovered a novel autophagy initiation pathway in vivo in neurons.

To examine the role of PIK3C3 in the central nervous system (CNS), we next deleted Pik3c3 in CNS neural progenitor cells using the Nestin-Cre transgenic line. The resulting conditional knockout mice displayed a severe cortical lamination abnormality caused by defective cortical neuron migration. This finding uncovered a previously under-appreciated role of endocytic trafficking in neural migration, which was further confirmed by electron microscopic analyses of the developing cortex. Moreover, overexpressing the dominant negative forms of Dynamin2 or Rab5, two regulators of endocytosis, caused similar migration defects as Pik3c3-deletion. Mechanistically, Pik3c3-deficient cortical neurons drastically reduced surface Reelin binding sites, and showed significantly decreased levels of Dab1 phosphorylation, despite expressing normal total amount of Reelin receptor ApoER2. This work suggests endocytosis and recycling of Reelin receptors are likely to play an important role in cortical migration regulated by the Reelin signaling pathway.

These studies represent the first in vivo characterization of PIK3C3 functions in mammals, and provide insight into the complexity and functional importance of neuronal endo-lysosomal and autophagic pathways.

Although evolution results from differential reproduction and survival at the level of the individual, most research in evolutionary genetics is concerned with comparisons made at the level of divergent populations or species. This is particularly true in work focused on the evolutionary genetics of natural populations. While this level of inquiry is extremely valuable, in order to develop a complete understanding of the evolutionary process we also need to understand how traits evolve within populations, on the level of differences between individuals, and in the context of natural ecological and environmental variation. A major difficulty confronting such work stems from the difficulty of assessing interindividual phenotypic variation and its sources within natural populations. This level of inquiry is, however, the main focus for many long-term field studies. Here, I take advantage of one such field study, centered on the wild baboon population of the Amboseli basin, Kenya, to investigate the possibilities for integrating functional, population, and evolutionary genetic approaches with behavioral, ecological, and environmental data. First, I describe patterns of hybridization and admixture in the Amboseli population, a potentially important component of population structure. Second, I combine field sampling, laboratory measurements of gene expression, and a computational approach to examine the possibility of using allele-specific gene expression as a tool to study functional regulatory variation in natural populations. Finally, I outline an example of how these and other methods can be used to understand the relationship between genetic variation and naturally occurring infection by a malaria-like parasite, Hepatocystis, also in the Amboseli baboons. The results of this work emphasize that developing genetic approaches for nonmodel genetic systems is becoming increasingly feasible, thus opening the door to pursuing such studies in behavioral and ecological model systems that provide a broader framework for genetic results. Integrating behavioral, ecological, and genetic perspectives will allow us to better appreciate the interplay between these different factors, and thus achieve a better understanding of the raw material upon which selection acts.

In Saccharomyces cerevisiae, proper maintenance of haploid cell identity requires the SIR complex to mediate the silenced chromatin found at the cryptic mating-type loci, HML and HMR. This complex consists of Sir2, a histone deacetylase and the histone binding proteins Sir3 and Sir4. Interestingly, both Sir2 and Sir3 have paralogs from a genome duplication that occurred after the divergence of Saccharomyces and Kluyveromyces species. The histone deacetylase HST1 is the paralog of SIR2 and works with the promoter-specific SUM1 complex to repress sporulation and alpha-specific genes. ORC1 is the paralog of SIR3 and is an essential subunit of the Origin Recognition Complex and also recruits SIR proteins to the HM loci. I have investigated the functions of these proteins in the non-duplicated species Kluyveromyces lactis and compared these functions to those found in S. cerevisiae.

I have shown that SIR2 and HST1 subfunctionalized post-duplication via the duplication, degeneration and complementation mechanism. In S. cerevisiae, Sir2 has retained the ability to function like Hst1 when in an hst1Δ strain. I have also shown, with a chimeric Sir2-Hst1 protein, that there are distinct specificity domains for Sir2 interaction with the SIR complex and Hst1 interaction with the SUM1 complex that have diverged between Sir2 and Hst1. Trans-species complementation assays show that the non-duplicated Sir2 from K. lactis can interact with both SIR and SUM1 complexes in S. cerevisiae.

Further analysis into the non-duplicated experimental system of K. lactis has revealed that deletion of KlSir2 de-represses the HM loci as well as sporulation and cell-type specific genes. A physical interaction between KlSir2 and the histone binding protein KlSir4 is conserved in K. lactis, and both proteins spread across the HML locus and associate with telomeres in a manner similar to S. cerevisiae. KlSir2 also physically interacts with the DNA-binding protein, KlSum1, to repress sporulation and cell-type specific genes in a promoter-specific manner and recruitment of KlSir2 to these loci is dependent on KlSum1. Surprisingly, deletion of KlSUM1 also de-represses HML and HMR, a phenotype not observed in S. cerevisiae. I show by chromatin immunoprecipitation that KlSum1 directly regulates the HM loci by spreading across these regions in a mechanism that is distinct from its role in repressing sporulation-specific genes. This result indicates that KlSum1 is a key regulator of not only meiotic, but also mating-type transcriptional programming.

The SIR3-ORC1 gene pair has previously been used as an example of neofunctionalization based on accelerated rates of evolution. However, my studies of KlOrc1 show it is distributed across HML and associates with Sir2 and Sir4 at telomeres, indicative of it having Sir3-like capabilities to spread across chromatin. This ability of KlOrc1 to spread is distinct from its functions with ORC, and is entirely dependent on its BAH domain. These findings demonstrate that prior to the genome duplication there was a silencing complex that contained both KlSir2 and KlOrc1. In addition to their functions at HML and the telomeres, KlOrc1 associates with replication origins and KlSir2 and KlSum1 work in complex to repress sporulation genes in a promoter-specific manner. The multiple functions of both KlOrc1 and KlSir2 in K. lactis indicate that after duplication, these properties were divided among paralogs and subsequently specialized to perform the functions that have been characterized in S. cerevisiae.

Although much research has focused on cognitive ability and the genetic and environmental factors that might influence it, this aspect of human nature is still far from being well understood. It has been well-established that certain factors such as age and education have significant impacts on performance on most cognitive tests, but the effects of variables such as cognitive pastimes and strategies used during testing have generally not been assessed. Additionally, no genetic variant has yet been unequivocally shown to influence the normal variation in cognitive ability of healthy individuals. Candidate gene studies of cognition have produced conflicting results that have not been replicable, and genome-wide association studies have not found common variants with large influences on this trait.

Here, we have recruited a large cohort of healthy volunteers (n=1,887) and administered a brief cognitive battery utilizing diverse, common, and well-known tests. In addition to providing standard demographic information, the subjects also filled out a questionnaire that was designed to assess novel factors such as whether they had seen the test before, in what cognitive pastimes they participated, and what strategies they had used during testing. Linear regression models were built to assess the effects of these variables on the test scores. I found that the addition of novel covariates to standard ones increased the percent of the variation in test score that was explained for all tests; for some tests, the increase was as high as 70%.

Next, I examined the effects of genetic variants on test scores. I first performed a genome-wide association study using the Illumina HumanHap 550 and 610 chips. These chips are designed to directly genotype or tag the vast majority of the common variants in the genome. Despite having 80% power to detect a common variant explaining at least 3-6% (depending on the test) of the variation in the trait, I did not find any genetic variants that were significantly associated after correction for multiple testing. This is in line with the general findings from GWA studies that single common variants have a limited impact on complex traits.

Because of the recent technological advances in next-generation sequencing and the apparently limited role of very common variants, many human geneticists are making a transition from genome-wide association study to whole-genome and whole-exome sequencing, which allow for the identification of rarer variants. Because these methods are currently costly, it is important to utilize study designs that have the best chance of finding causal variants in a small sample size. One such method is the extreme-trait design, where individuals from one or both ends of a trait distribution are sequenced and variants that are enriched in the group(s) are identified. Here, I have sequenced the exomes of 20 young individuals of European ethnicity: 10 that performed at the top of the distribution for the cognitive battery and 10 that performed at the bottom. I identified rare genetic variants that were enriched in one extreme group as compared to the other and performed follow-up genotyping of the best candidate variant that emerged from this analysis. Unfortunately, this variant was not found to be associated in a larger sample of individuals. This pilot study indicates that a larger sample size will be needed to identify variants enriched in cognition extremes.

Finally, I assessed the effect of topiramate, an antiepileptic drug that causes marked side effects in certain cognitive areas in certain individuals, on some of the healthy volunteers (n=158) by giving them a 100 mg dose and then administering the cognitive test two hours later. I compared their scores at this testing session to those at the previous session and calculated the overall level to which they were affected by topiramate. I found that the topiramate blood levels, which were highly dependent on weight and the time from dosing to testing, varied widely between individuals after this acute dose, and that this variation explained 35% of the variability in topiramate response. A genome-wide association study of the remaining variability in topiramate response did not identify a genome-wide significant association.

In sum, I studied the contributions of both environmental and genetic variables to cognitive ability and cognitive response to topiramate. I found that I could identify environmental variables explaining large proportions of the variation in these traits, but that I could not identify genetic variants that influenced the traits. My analysis of genetic variants was for the most part restricted to the very common ones found on genotyping chips, and this and other studies have generally found that single common genetic variants do not have large affects on complex traits. As we move forward into studies that involve the sequencing of whole exomes and genomes, genetic variants with large effects on these complex traits may finally be found.