Probing Evolutionary Repeatability: Neutral and Double Changes and the Predictability of Evolutionary Adaptation

The question of how organisms adapt is among the most fundamental in evolutionary biology. Two recent studies investigated the evolution of Escherichia coli in response to challenge with the antibiotic cefotaxime. Studying five mutations in the beta-lactamase gene that together confer significant antibiotic resistance, the authors showed a complex fitness landscape that greatly constrained the identity and order of intermediates leading from the initial wildtype genotype to the final resistant genotype. Out of 18 billion possible orders of single mutations leading from non-resistant to fully-resistant form, they found that only 27 (1.5x10(-7)%) pathways were characterized by consistently increasing resistance, thus only a tiny fraction of possible paths are accessible by positive selection. I further explore these data in several ways.Allowing neutral changes (those that do not affect resistance) increases the number of accessible pathways considerably, from 27 to 629. Allowing multiple simultaneous mutations also greatly increases the number of accessible pathways. Allowing a single case of double mutation to occur along a pathway increases the number of pathways from 27 to 259, and allowing arbitrarily many pairs of simultaneous changes increases the number of possible pathways by more than 100 fold, to 4800. I introduce the metric 'repeatability,' the probability that two random trials will proceed via the exact same pathway. In general, I find that while the total number of accessible pathways is dramatically affected by allowing neutral or double mutations, the overall evolutionary repeatability is generally much less affected.These results probe the conceivable pathways available to evolution. Even when many of the assumptions of the analysis of Weinreich et al. (2006) are relaxed, I find that evolution to more highly cefotaxime resistant beta-lactamase proteins is still highly repeatable

Transcriptomic analysis of diplomonad parasites reveals a trans-spliced intron in a helicase gene in Giardia

ABSTRACT Background. The mechanisms by which DNA sequences are expressed is the central preoccupation of molecular genetics. Recently, ourselves and others reported that in the diplomonad protist Giardia lamblia, the coding regions of several mRNAs are produced by ligation of independent RNA species expressed from distinct genomic loci. Such trans-splicing of introns was found to affect nearly as many genes in this organism as does classical cis-splicing of introns. These findings raised questions about the incidence of intron trans-splicing both across the G. lamblia transcriptome and across diplomonad diversity in general, however a dearth of transcriptomic data at the time prohibited systematic study of these questions. Methods. I leverage newly available transcriptomic data from G. lamblia and the related diplomonad Spironucleus salmonicida to search for trans-spliced introns. My computational pipeline recovers all four previously reported trans-spliced introns in G. lamblia, suggesting good sensitivity. Results. Scrutiny of thousands of potential cases revealed only a single additional trans-spliced intron in G. lamblia, in the p68 helicase gene, and no cases in S. salmonicida. The p68 intron differs from the previously reported trans-spliced introns in its high degree of streamlining: the core features of G. lamblia trans-spliced introns are closely packed together, revealing striking economy in the implementation of a seemingly inherently uneconomical molecular mechanism. Discussion. These results serve to circumscribe the role of trans-splicing in diplomonads both in terms of the number of genes effected and taxonomically. Future work should focus on the molecular mechanisms, evolutionary origins and phenotypic implications of this intriguing phenomenon

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions

Intron length distributions and gene prediction

Accurate gene prediction in eukaryotes is a difficult and subtle problem. Here we point out a useful feature of expected distributions of spliceosomal intron lengths. Since introns are removed from transcripts prior to translation, intron lengths are not expected to respect coding frame, thus the number of genomic introns that are a multiple of three bases (‘3n introns’) should be similar to the number that are a multiple of three plus one bases (or plus two bases). Skewed predicted intron length distributions thus suggest systematic errors in intron prediction. For instance, a genome-wide excess of 3n introns suggests that many internal exonic sequences have been incorrectly called introns, whereas a deficit of 3n introns suggests that many 3n introns that lack stop codons have been mistaken for exonic sequence. A survey of genomic annotations for 29 diverse eukaryotic species showed that skew in intron length distributions is a common problem. We discuss several examples of skews in genome-wide intron length distributions that indicate systematic problems with gene prediction. We suggest that evaluation of length distributions of predicted introns is a fast and simple method for detecting a variety of possible systematic biases in gene prediction or even problems with genome assemblies, and discuss ways in which these insights could be incorporated into genome annotation protocols

Functional and evolutionary analysis of alternatively spliced genes is consistent with an early eukaryotic origin of alternative splicing

<p>Abstract</p> <p>Background</p> <p>Alternative splicing has been reported in various eukaryotic groups including plants, apicomplexans, diatoms, amoebae, animals and fungi. However, whether widespread alternative splicing has evolved independently in the different eukaryotic groups or was inherited from their last common ancestor, and may therefore predate multicellularity, is still unknown. To better understand the origin and evolution of alternative splicing and its usage in diverse organisms, we studied alternative splicing in 12 eukaryotic species, comparing rates of alternative splicing across genes of different functional classes, cellular locations, intron/exon structures and evolutionary origins.</p> <p>Results</p> <p>For each species, we find that genes from most functional categories are alternatively spliced. Ancient genes (shared between animals, fungi and plants) show high levels of alternative splicing. Genes with products expressed in the nucleus or plasma membrane are generally more alternatively spliced while those expressed in extracellular location show less alternative splicing. We find a clear correspondence between incidence of alternative splicing and intron number per gene both within and between genomes. In general, we find several similarities in patterns of alternative splicing across these diverse eukaryotes.</p> <p>Conclusion</p> <p>Along with previous studies indicating intron-rich genes with weak intron boundary consensus and complex spliceosomes in ancestral organisms, our results suggest that at least a simple form of alternative splicing may already have been present in the unicellular ancestor of plants, fungi and animals. A role for alternative splicing in the evolution of multicellularity then would largely have arisen by co-opting the preexisting process.</p

Exam Setting and Moderation in a Transnational Engineering Programme

The practice of writing effective exam questions has attracted plenty of interest in the literature. However, the process of exam setting and moderation has received little attention. Thus, the aim of this manuscript is to describe our approach in exam setting and moderation for a transnational education (TNE) programme between the University of Glasgow and the University of Electronic Science and Technology of China (UESTC). Each year, Glasgow University staff members prepare main and resit exam papers for 32 courses in the electronic and electrical engineering subject group. To ensure a smooth workflow and a relatively quick turnover of high-quality exam papers each year, we describe the best practices in exam setting, as well as our approach in exam moderation using a new moderation form and a new system for uploading exam papers online. This new system was designed to facilitate a solid dialogue between exam setters and examiners, so that all summative assessments are fair, consistent and reliable. These best practices can be used by engineering faculty to make informed decisions regarding the process of exam setting and moderation in a transnational engineering programme

Arsenic Speciation in Arsenic-Rich Brazilian Soils from Gold Mining Sites under Anaerobic Incubation

Arsenic enters the environment naturally though rock weathering but its concentrations can be dramatically increased by human activity such as mining. ISTC’s senior chemist, John Scott, and collaborators from the Federal University of Viçosa (Brazil), Illinois State Geological Survey, and the University of Illinois at Urbana-Champaign teamed up to find the arsenic release pathways in arsenic rich soils at gold mining sites in Brazil. Full research results available in de Mello, J.W.V., Talbott, J.L., Scott, J. et al. (2007). "Arsenic speciation in arsenic-rich Brazilian soils from gold mining sites under anaerobic incubation." Environmental Science and Pollution Research - International 14(6), 388-396. https://doi.org/10.1065/espr2006.08.330Ope

SPR-measured dissociation kinetics of PROTAC ternary complexes influence target degradation rate

Bifunctional degrader molecules, known as proteolysis-targeting chimeras (PROTACs), function by recruiting a target to an E3 ligase, forming a target/PROTAC/ligase ternary complex. Despite the importance of this key intermediate species, no detailed validation of a method to directly determine binding parameters for ternary complex kinetics has been reported, and it remains to be addressed whether tuning the kinetics of PROTAC ternary complexes may be an effective strategy to improve the efficiency of targeted protein degradation. Here, we develop an SPR-based assay to quantify the stability of PROTAC-induced ternary complexes by measuring for the first time the kinetics of their formation and dissociation <i>in vitro</i> using purified proteins. We benchmark our assay using four PROTACs that target the bromodomains (BDs) of bromodomain and extraterminal domain proteins Brd2, Brd3, and Brd4 to the von Hippel–Lindau E3 ligase (VHL). We reveal marked differences in ternary complex off-rates for different PROTACs that exhibit either positive or negative cooperativity for ternary complex formation relative to binary binding. The positively cooperative degrader MZ1 forms comparatively stable and long-lived ternary complexes with either Brd4^BD2 or Brd2^BD2 and VHL. Equivalent complexes with Brd3^BD2 are destabilized due to a single amino acid difference (Glu/Gly swap) present in the bromodomain. We observe that this difference in ternary complex dissociative half-life correlates to a greater initial rate of intracellular degradation of Brd2 and Brd4 relative to Brd3. These findings establish a novel assay to measure the kinetics of PROTAC ternary complexes and elucidate the important kinetic parameters that drive effective target degradation

Complex selection on 5' splice sites in intron-rich organisms

In contrast to the typically streamlined genomes of prokaryotes, many eukaryotic genomes are riddled with long intergenic regions, spliceosomal introns, and repetitive elements. What explains the persistence of these and other seemingly suboptimal structures? There are three general hypotheses: (1) the structures in question are not actually suboptimal but optimal, being favored by selection, for unknown reasons; (2) the structures are not suboptimal, but of (essentially) equal fitness to 'optimal' ones; or (3) the structures are truly suboptimal, but selection is too weak to systematically eliminate them. The 5' splice sites of introns offer a rare opportunity to directly test these hypotheses. Intron-poor species show a clear consensus splice site; most introns begin with the same six nucleotide sequence (typically GTAAGT or GTATGT), indicating efficient selection for this consensus sequence. In contrast, intron-rich species have much less pronounced boundary consensus sequences, and only small minorities of introns in intron-rich species share the same boundary sequence. We studied rates of evolutionary change of 5' splice sites in three groups of closely related intron-rich species--three primates, five Drosophila species, and four Cryptococcus fungi. Surprisingly, the results indicate that changes from consensus-to-variant nucleotides are generally disfavored by selection, but that changes from variant to consensus are neither favored nor disfavored. This evolutionary pattern is consistent with selective differences across introns, for instance, due to compensatory changes at other sites within the gene, which compensate for the otherwise suboptimal consensus-to-variant changes in splice boundaries