It has been theorized that most proteins are under selection pressure to be soluble in crowded cellular spaces. To maintain solubility a proteins’ aggregation propensity should be inversely proportional to their maximum likely concentration. This theory was examined by comparing the proteome of the model organism S. cerevisiae, which has previously undergone a Whole Genome Duplication (WGD) event to the proteome of the closely related yeast K. waltii, which has not undergone WGD. This comparison revealed the following: 1) Predicted aggregation propensities are higher in S. cerevisiae than K. waltii. 2) Aggregation propensity does not predict which genes reverted to a single copy after WGD. 3) In genes which were retained as duplicates in S. cerevisiae after WGD, aggregation propensities rose from the inferred common ancestral protein. 4) Genes retained as duplicates showed less of an increase relative to their homologues in K. waltii than genes which were not retained as duplicates. 5) The relationship between the log predicted aggregation propensity and log mRNA expression level or log protein abundance was not linear as previously predicted. These results suggest that while there is broad selection pressure for reduced aggregation pressure for genes which have been duplicated, the precise relationship between aggregation propensity and gene expression is more complicated than previously predicted. These results also allow speculation that the whole genome duplication in S.cerevisiae may have been made possible by a general relaxation of aggregation-related selection pressure. / text

Inorganic phosphate is an essential nutrient required by all organisms for optimal growth. During phosphate starvation, Saccharomyces cerevisiae induces a set of genes responsible for the regulation of inorganic phosphate acquisition. The phosphate-responsive signaling (PHO) pathway controls this response, with the CDK-cyclin complex Pho85-Pho80 playing a prominent role. Here we report the X-ray structure of the Pho85-Pho80 complex, identifying the unique structural features that distinguish it from other cell cycle associated CDK-cyclin complexes. The structure reveals a specific salt bridge between a Pho85 arginine and a Pho80 aspartate that maintains a Pho80 loop confirmation important for substrate recognition and makes phosphorylation of the Pho85 activation loop dispensible. We show that a cluster of residues distal to the kinase active site are involved in a high affinity interaction between the Pho80 cyclin and the transcription factor substrate (Pho4). The structure also reveals a separate high affinity binding site for the CDK inhibitor (Pho81). The fission yeast, Schizosaccharomyces pombe, regulates expression of the secreted acid phosphatase \((pho1^+)\) via a non-orthologous PHO pathway. The genes induced by phosphate limitation and the molecular mechanism by which the genetically identified positive \((pho7^+)\) and negative \((csk1^+)\) regulators function are not known. Here we use a combination of molecular biology, expression microarrays, chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq), and global transcriptome sequencing (RNA-Seq) to characterize the role of \(pho7^+\) and \(csk1^+\) in the PHO response. We show that there is a fast and slow response to phosphate starvation, each with defined regulatory roles. We use ChIP-Seq to identify members of the Pho7 regulon and characterize Pho7 binding dynamics in response to phosphate-limitation and Csk1 activity. We identify a conserved PHO response for the PHO5 \((pho1^+)\), PHO84 \((spbc8e4.01c^+)\), and GIT1 \((spbc1271.09^+)\) orthologs. We show that activation of \(pho1^+\) requires Pho7 binding to a UAS in the \(pho1^+\) promoter and that a URS is necessary for Csk1 repression. We find that Pho7-dependent activation is not limited to phosphate-starvation, as additional environmental stress response pathways require \(pho7^+\) for maximal induction. Using RNA-Seq we show that Pho7 is also involved in regulating non-coding transcription and is a bi-functional transcription factor.

No description available.

The yeast Saccharomyces cerevisiae genome is descendant from a whole genome duplication event approximately 150 million years ago. Following this duplication many genes were lost however, a certain class of genes, termed ohnologs, persist in duplicate. In this thesis we investigate network centrality as it relates to ohnolog re- tention with the goal of determining why only certain genes were retained. With this in mind, we compare physical and genetic interaction networks and genetic and pro- tein sequence data in order to reveal how network characteristics and post-duplication retention are related. We show that there are two subclasses of ohnologs, those that interact with their duplication sister and those that do not and that these two classes have distinct characteristics that provide insight into the evolutionary mechanisms that affected their retention following whole genome duplication. Namely, a very low ratio of non-synonymous mutations per non-synonymous site for ohnologs that retain an interaction with their duplicate. The opposite observation is seen for ohnologs that have lost their interaction with their duplicate. We interpret this in the fol- lowing way: ohnologs that have retained their interaction with their duplicate are functionally constrained to buffer for the other ohnolog. For this reason they are retained; ohnologs that have lost their interaction with their duplicate are retained because they are functionally divergent to the point of being individually essential. Additionally we investigate small scale duplications and show that, generally, the mechanism of duplication (smale scale or whole genomes) does not affect the distri- bution of network characteristics. Nor do these network characteristics correlate to the selective pressure observed by retained paralogous genes, including both ohnologs and small scale duplicates. In contrast, we show that the network characteristics of individual genes, particularly the magnitude of their physical and genetic network centralities, do influence their retention following whole genome duplication. / Graduate / mjrimrie@gmail.com

Paralog reduction, the loss of duplicate genes after whole genome duplication (WGD) is a pervasive process. Whether this loss proceeds gene by gene or through deletion of multi-gene DNA segments is controversial, as is the question of fractionation bias, namely whether one homeologous chromosome is more vulnerable to gene deletion than the other. As a null hypothesis, we first assume deletion events, on one homeolog only, excise a geometrically distributed number of genes with unknown mean mu, and these events combine to produce deleted runs of length l, distributed approximately as a negative binomial with unknown parameter r; itself a random variable with distribution pi(.). A biologically more realistic model requires deletion events on both homeologs distributed as a truncated geometric. We simulate the distribution of run lengths l in both models, as well as the underlying pi(r), as a function of mu, and show how sampling l allows us to estimate mu. We apply this to data on a total of 15 genomes descended from 6 distinct WGD events and show how to correct the bias towards shorter runs caused by genome rearrangements. Because of the difficulty in deriving pi(.) analytically, we develop a deterministic recurrence to calculate each pi(r) as a function of mu and the proportion of unreduced paralog pairs. This is based on a computing formula containing nested sums. The parameter mu can be estimated based on run lengths of single-copy regions. We then reduce the computing formulae, at least in the one-sided case, to closed form. This virtually eliminates computing time due to highly nested summations. We formulate a continuous version of the fractionation process, deleting line segments of exponentially distributed lengths in analogy to geometric distributed numbers of genes. We derive nested integrals and discover that the number of previously deleted regions to be skipped by a new deletion event is exactly geometrically distributed. We undertook a large simulation experiment to show how to discriminate between the gene-by-gene duplicate deletion model and the deletion of a geometrically distributed number of genes. This revealed the importance of the effects of genome size N, the mean of the geometric distribution, the progress towards completion of the fractionation process, and whether the data are based on runs of deleted genes or undeleted genes.

The study of aflatoxin in Aspergillus spp. has garnered the attention of many researchers due to aflatoxin's carcinogenic properties and frequency as a food and feed contaminant. Significant progress has been made by utilizing the model organism Aspergillus nidulans to characterize the regulation of sterigmatocystin (ST), the penultimate precursor of aflatoxin. A previous forward genetic screen identified 23 A. nidulans mutants involved in regulating ST production. Six mutants were characterized from this screen using classical mapping (five mutations in mcsA) and complementation with a cosmid library (one mutation in laeA). The remaining mutants were backcrossed and sequenced using Illumina and Ion Torrent sequencing platforms. All but one mutant contained one or more sequence variants in predicted open reading frames. Deletion of these genes resulted in identification of mutant alleles responsible for the loss of ST production in 12 of the 17 remaining mutants. Eight of these mutations were in genes already known to affect ST synthesis (laeA, mcsA, fluG, and stcA), while the remaining four mutations (in laeB, sntB, and hamI) were in previously uncharacterized genes not known to be involved in ST production. Deletion of laeB, sntB, and hamI in A. flavus results in loss of aflatoxin production, confirming that these regulators are conserved in the aflatoxigenic aspergilli. This report highlights the multifaceted regulatory mechanisms governing secondary metabolism in Aspergillus. Additionally, these data contribute to the increasing number of studies showing that forward genetic screens of fungi coupled with whole-genome resequencing is a robust and cost-effective technique. IMPORTANCE In a postgenomic world, reverse genetic approaches have displaced their forward genetic counterparts. The techniques used in forward genetics to identify loci of interest were typically very cumbersome and time-consuming, relying on Mendelian traits in model organisms. The current work was pursued not only to identify alleles involved in regulation of secondary metabolism but also to demonstrate a return to forward genetics to track phenotypes and to discover genetic pathways that could not be predicted through a reverse genetics approach. While identification of mutant alleles from whole-genome sequencing has been done before, here we illustrate the possibility of coupling this strategy with a genetic screen to identify multiple alleles of interest. Sequencing of classically derived mutants revealed several uncharacterized genes, which represent novel pathways to regulate and control the biosynthesis of sterigmatocystin and of aflatoxin, a societally and medically important mycotoxin.

Paralog reduction, the loss of duplicate genes after whole genome duplication (WGD) is a pervasive process. Whether this loss proceeds gene by gene or through deletion of multi-gene DNA segments is controversial, as is the question of fractionation bias, namely whether one homeologous chromosome is more vulnerable to gene deletion than the other. As a null hypothesis, we first assume deletion events, on one homeolog only, excise a geometrically distributed number of genes with unknown mean mu, and these events combine to produce deleted runs of length l, distributed approximately as a negative binomial with unknown parameter r; itself a random variable with distribution pi(.). A biologically more realistic model requires deletion events on both homeologs distributed as a truncated geometric. We simulate the distribution of run lengths l in both models, as well as the underlying pi(r), as a function of mu, and show how sampling l allows us to estimate mu. We apply this to data on a total of 15 genomes descended from 6 distinct WGD events and show how to correct the bias towards shorter runs caused by genome rearrangements. Because of the difficulty in deriving pi(.) analytically, we develop a deterministic recurrence to calculate each pi(r) as a function of mu and the proportion of unreduced paralog pairs. This is based on a computing formula containing nested sums. The parameter mu can be estimated based on run lengths of single-copy regions. We then reduce the computing formulae, at least in the one-sided case, to closed form. This virtually eliminates computing time due to highly nested summations. We formulate a continuous version of the fractionation process, deleting line segments of exponentially distributed lengths in analogy to geometric distributed numbers of genes. We derive nested integrals and discover that the number of previously deleted regions to be skipped by a new deletion event is exactly geometrically distributed. We undertook a large simulation experiment to show how to discriminate between the gene-by-gene duplicate deletion model and the deletion of a geometrically distributed number of genes. This revealed the importance of the effects of genome size N, the mean of the geometric distribution, the progress towards completion of the fractionation process, and whether the data are based on runs of deleted genes or undeleted genes.

Naturally competent bacterial species, which self-induce the recombination mechanism of transformation, are wide spread across the bacterial tree-of-life. However, it remains unclear why competence has evolved in these bacteria. Although it is likely that exact explanations will be different for each species, a common selective factor cannot be excluded. Currently, three dominant hypotheses, which focus on the transformation function, try to explain the benefits of competence. Firstly, competence is thought to increase the rate of adaptation by combining beneficial alleles in single genotypes. Secondly, competence can repair DNA-damage by replacing the damaged DNA fragments with undamaged ones. Thirdly, the DNA uptake during competence is used to recycle environmental DNA fragments for nutrients. One of the naturally competent species is the Gram-positive Streptococcus pneumoniae, which is an opportunistic pathogen generally inhabiting the naso-pharyngeal area of young children. Competence in S. pneumoniae is regulated via density dependent extracellular signaling peptide. Here I use a combination of experiments designed around knockout mutants of the signaling mechanism and next-generation sequencing methods to test the first two hypotheses in S. pneumoniae. First, I extend on the DNA-for-repair hypothesis by showing that competent populations of S. pneumoniae are better protected not only against a DNA-damaging agent, but also against protein synthesis inhibitors. However, the mechanisms underlying this protection differ between types of stress. DNA-damage requires the full process of transformation, while protection against protein synthesis inhibitors only requires the activation of the competent cell state. This shows that benefits of competence cannot be totally explained by the benefits of transformation. Second, I use a long-term evolution experiment, where competent and non-competent strains are kept in the presence and absence of periodic stress, to determine the importance of competence for the generation of genetic variation. I find that competence does not increase the rate of adaptation in S. pneumoniae. The fitness of evolved competent populations was significantly lower than those of non-competent populations evolved over the same period of time. However, the intrinsic costs of competence are mitigated by the addition of short periods of stress exposure. These results confirm the prediction of the fitness associated recombination (FAR) hypothesis that competence is favoured in low-fitness situations. Thirdly, whole genome re-sequencing of the evolved populations allowed me to explore genomic evolution next to fitness changes. The genomic data revealed that competence reduces the mutational load of deleterious mutations rather than generating combinations of beneficial alleles. In addition I show several case of parallel genomic evolution within each treatment and across treatments. This shows that parallel evolution is not restricted by genotypic background (competence) or environment (periodic stress). Finally, these results show that competence has evolved in populations of S. pneumoniae as a mechanism to deal with various forms of stress.

Mycobacterium tuberculosis (Mtb) remains a global challenge that has been worsened by the emergence of drug resistant strains of Mtb. We used publicly available Next Generation Sequencing (NGS) and drug susceptibility (DST) data to develop “Resistance sniffer”, an online software program for the rapid prediction of lineage and Mtb drug resistance. Based on the distribution of polymorphisms in the genomes of Mtb, we calculated the power of association between the polymorphisms in different clades of Mtb and resistance to 13 anti-TB drugs. Our data suggests that the development of drug resistance in Mtb is a stepwise process that involves the accumulation of polymorphisms in the Mtb genome. We carefully curated the polymorphisms based on their association powers to create a diagnostic key that captures patterns of these polymorphisms that can be used to predict lineage and drug resistance in Mtb. This diagnosis key was incorporated into the Resistance Sniffer tool, an online software program that we developed for the rapid diagnosis of drug resistance in Mtb. The tool was tested using sequence data from the South Africa Medical Research Council (SA-MRC). Our data suggests that the majority of the strains in SA may have been brought by the arrival of European settlers while the more resistant strains may have been introduced in the region by Asian travellers later on. We next sought to determine non-random associations between polymorphic sites in genomes of Mtb. Using the attributable risk (Ra) statistical methods, we distinguished between functional associations and associations that may have been due to genetic drift events for different Mtb clades. We then integrated the (Ra) data with drug susceptibility and annotation data to generate networks in Cytoscape 3.71. These networks were then used to infer evolutionary trajectories that drive the emergence and fixation of the drug resistant phenotype in different clades of Mtb. We demonstrate that strains from the Lineage 1.2 are associated with less complex functional associations compared to the strains from other clades such as the Asian and Euro-American clades. Our data also shows that the predisposition of strains from the Asian clades to develop multi-drug resistance may be attributed to a complex network of functional interactions of mutations in genes that are involved in several aspects of Mtb physiology such as cell wall modelling, lipid metabolism, stress response and DNA repair. / Dissertation (MSc)--University of Pretoria, 2019. / Biochemistry / MSc / Unrestricted

Population dynamics of large ungulates are complex and vary with fluctuations in factors such as predation, resource availability, human disturbance, and weather (Gaillard et al. 1998, Forrester and Wittmer 2013). These regulating factors exhibit similar effects on ungulate populations by changing vital rates such as birthrate, death rate, emigration or immigration (Gaillard et al. 2000). To better understand the mechanisms influencing population change, it is useful to involve tools from multiple disciplines (Krausman et al. 2013). Here we explore population dynamics of mule deer (Odocoileus hemionus) through the lenses of two distinct fields: population ecology and genomics. In the first chapter we examine the influence of maternal effects on offspring fitness. In the second chapter we present a high-quality, chromosome-level reference genome for mule deer. We expect results from each of these studies to provide valuable resources for continued research and conservation of mule deer.