Timothy B. Sackton

A Scan for Positively Selected Genes in the Genomes of Humans and Chimpanzees

Since the divergence of humans and chimpanzees about 5 million years ago, these species have undergone a remarkable evolution with drastic divergence in anatomy and cognitive abilities. At the molecular level, despite the small overall magnitude of DNA sequence divergence, we might expect such evolutionary changes to leave a noticeable signature throughout the genome. We here compare 13,731 annotated genes from humans to their chimpanzee orthologs to identify genes that show evidence of positive selection. Many of the genes that present a signature of positive selection tend to be involved in sensory perception or immune defenses. However, the group of genes that show the strongest evidence for positive selection also includes a surprising number of genes involved in tumor suppression and apoptosis, and of genes involved in spermatogenesis. We hypothesize that positive selection in some of these genes may be driven by genomic conflict due to apoptosis during spermatogenesis. Genes with maximal expression in the brain show little or no evidence for positive selection, while genes with maximal expression in the testis tend to be enriched with positively selected genes. Genes on the X chromosome also tend to show an elevated tendency for positive selection. We also present polymorphism data from 20 Caucasian Americans and 19 African Americans for the 50 annotated genes showing the strongest evidence for positive selection. The polymorphism analysis further supports the presence of positive selection in these genes by showing an excess of high-frequency derived nonsynonymous mutations.

Dynamic evolution of the innate immune system in Drosophila

The availability of complete genome sequence from 12 Drosophila species presents the opportunity to examine how natural selection has affected patterns of gene family evolution and sequence divergence among different components of the innate immune system. We have identified orthologs and paralogs of 245 Drosophila melanogaster immune-related genes in these recently sequenced genomes. Genes encoding effector proteins, and to a lesser extent genes encoding recognition proteins, are much more likely to vary in copy number across species than genes encoding signaling proteins. Furthermore, we can trace the apparent recent origination of several evolutionarily novel immune-related genes and gene families. Using codon-based likelihood methods, we show that immune-system genes, and especially those encoding recognition proteins, evolve under positive darwinian selection. Positively selected sites within recognition proteins cluster in domains involved in recognition of microorganisms, suggesting that molecular interactions between hosts and pathogens may drive adaptive evolution in the Drosophila immune system.

Evolution of protein-coding genes in Drosophila

Several contributing factors have been implicated in evolutionary rate heterogeneity among proteins, but their evolutionary mechanisms remain poorly characterized. The recently sequenced 12 Drosophila genomes provide a unique opportunity to shed light on these unresolved issues. Here, we focus on the role of natural selection in shaping evolutionary rates. We use the Drosophila genomic data to distinguish between factors that increase the strength of purifying selection on proteins and factors that affect the amount of positive selection experienced by proteins. We confirm the importance of translational selection in shaping protein evolution in Drosophila and show that factors such as tissue bias in expression, gene essentiality, intron number, and recombination rate also contribute to evolutionary rate variation among proteins.

CYTONUCLEAR COADAPTATION IN DROSOPHILA: DISRUPTION OF CYTOCHROME C OXIDASE ACTIVITY IN BACKCROSS GENOTYPES

Abstract The cytochrome c oxidase enzyme (COX) is comprised of 10 nuclear‐encoded subunits and three mito‐chondrial‐encoded subunits in close physical association in the inner mitochondrial membrane. COX passes electrons from cytochrome c to molecular oxygen and pumps protons into the inner mitochondrial space for ATP production. Selection on nuclear‐mitochondrial interactions within species should lead to coadaptation of the proteins comprising this important enzyme. Under this model, there should be relatively little disruption of COX activity when mitochondrial genomes are crossed among strains within species. A more pronounced disruption of activity is expected when the mitochondrial genome is expressed in the nuclear background of a different species. We test these hypotheses in Drosophila using hybridization and backcrossing among lines of D. simulans and D. mauritiana. Disrupted cytonuclear genotypes were constructed using backcrosses between two lines of D. simulans (siI and si II) that introduced each divergent mitochondrial DNA (mtDNA) into each nuclear background due to maternal inheritance of mtDNA. Similar crosses were used to introduce eachD. simulans mtDNA into the D. mauritiana maI nuclear background. Reconstituted cytonuclear control genotypes were constructed by backcrossing the initial F1 females to males of the maternal genotype. COX enzyme activities were compared among these disrupted and reconstituted backcross genotypes within and between species. The disruption effect on COX activity was restricted to males of interspecific genotypes. These data support the coadaptation hypothesis and are consistent with predictions that the evolution of modifiers of male mitochondrial dysfunction is hindered by the maternal inheritance of mtDNA. New sequence data for nuclear encoded subunits of COX identified amino acids that may play a role in the disruption effect.

Genome of the house fly, Musca domestica L., a global vector of diseases with adaptations to a septic environment

BackgroundAdult house flies, Musca domestica L., are mechanical vectors of more than 100 devastating diseases that have severe consequences for human and animal health. House fly larvae play a vital role as decomposers of animal wastes, and thus live in intimate association with many animal pathogens.ResultsWe have sequenced and analyzed the genome of the house fly using DNA from female flies. The sequenced genome is 691 Mb. Compared with Drosophila melanogaster, the genome contains a rich resource of shared and novel protein coding genes, a significantly higher amount of repetitive elements, and substantial increases in copy number and diversity of both the recognition and effector components of the immune system, consistent with life in a pathogen-rich environment. There are 146 P450 genes, plus 11 pseudogenes, in M. domestica, representing a significant increase relative to D. melanogaster and suggesting the presence of enhanced detoxification in house flies. Relative to D. melanogaster, M. domestica has also evolved an expanded repertoire of chemoreceptors and odorant binding proteins, many associated with gustation.ConclusionsThis represents the first genome sequence of an insect that lives in intimate association with abundant animal pathogens. The house fly genome provides a rich resource for enabling work on innovative methods of insect control, for understanding the mechanisms of insecticide resistance, genetic adaptation to high pathogen loads, and for exploring the basic biology of this important pest. The genome of this species will also serve as a close out-group to Drosophila in comparative genomic studies.

Natural Selection Constrains Neutral Diversity across A Wide Range of Species

The neutral theory of molecular evolution predicts that the amount of neutral polymorphisms within a species will increase proportionally with the census population size (Nc). However, this prediction has not been borne out in practice: while the range of Nc spans many orders of magnitude, levels of genetic diversity within species fall in a comparatively narrow range. Although theoretical arguments have invoked the increased efficacy of natural selection in larger populations to explain this discrepancy, few direct empirical tests of this hypothesis have been conducted. In this work, we provide a direct test of this hypothesis using population genomic data from a wide range of taxonomically diverse species. To do this, we relied on the fact that the impact of natural selection on linked neutral diversity depends on the local recombinational environment. In regions of relatively low recombination, selected variants affect more neutral sites through linkage, and the resulting correlation between recombination and polymorphism allows a quantitative assessment of the magnitude of the impact of selection on linked neutral diversity. By comparing whole genome polymorphism data and genetic maps using a coalescent modeling framework, we estimate the degree to which natural selection reduces linked neutral diversity for 40 species of obligately sexual eukaryotes. We then show that the magnitude of the impact of natural selection is positively correlated with Nc, based on body size and species range as proxies for census population size. These results demonstrate that natural selection removes more variation at linked neutral sites in species with large Nc than those with small Nc and provides direct empirical evidence that natural selection constrains levels of neutral genetic diversity across many species. This implies that natural selection may provide an explanation for this longstanding paradox of population genetics.

Genetic Variation in Drosophila melanogaster Resistance to Infection: A Comparison Across Bacteria

Insects use a generalized immune response to combat bacterial infection. We have previously noted that natural populations of D. melanogaster harbor substantial genetic variation for antibacterial immunocompetence and that much of this variation can be mapped to genes that are known to play direct roles in immunity. It was not known, however, whether the phenotypic effects of variation in these genes are general across the range of potentially infectious bacteria. To address this question, we have reinfected the same set of D. melanogaster lines with Serratia marcescens, the bacterium used in the previous study, and with three additional bacteria that were isolated from the hemolymph of wild-caught D. melanogaster. Two of the new bacteria, Enterococcus faecalis and Lactococcus lactis, are gram positive. The third, Providencia burhodogranaria, is gram negative like S. marcescens. Drosophila genotypes vary highly significantly in bacterial load sustained after infection with each of the four bacteria, but mean loads are largely uncorrelated across bacteria. We have tested statistical associations between immunity phenotypes and nucleotide polymorphism in 21 candidate immunity genes. We find that molecular variation in some genes, such as Tehao, contributes to phenotypic variation in the suppression of only a subset of the pathogens. Variation in SR-CII and 18-wheeler, however, has effects that are more general. Although markers in SR-CII and 18-wheeler explain >20% of the phenotypic variation in resistance to L. lactis and E. faecalis, respectively, most of the molecular polymorphisms tested explain <10% of the total variance in bacterial load sustained after infection.

Synergistic blockade of mitotic exit by two chemical inhibitors of the APC/C

Protein machines are multi-subunit protein complexes that orchestrate highly regulated biochemical tasks. An example is the anaphase-promoting complex/cyclosome (APC/C), a 13-subunit ubiquitin ligase that initiates the metaphase–anaphase transition and mitotic exit by targeting proteins such as securin and cyclin B1 for ubiquitin-dependent destruction by the proteasome. Because blocking mitotic exit is an effective approach for inducing tumour cell death, the APC/C represents a potential novel target for cancer therapy. APC/C activation in mitosis requires binding of Cdc20 (ref. 5), which forms a co-receptor with the APC/C to recognize substrates containing a destruction box (D-box). Here we demonstrate that we can synergistically inhibit APC/C-dependent proteolysis and mitotic exit by simultaneously disrupting two protein–protein interactions within the APC/C–Cdc20–substrate ternary complex. We identify a small molecule, called apcin (APC inhibitor), which binds to Cdc20 and competitively inhibits the ubiquitylation of D-box-containing substrates. Analysis of the crystal structure of the apcin–Cdc20 complex suggests that apcin occupies the D-box-binding pocket on the side face of the WD40-domain. The ability of apcin to block mitotic exit is synergistically amplified by co-addition of tosyl-l-arginine methyl ester, a small molecule that blocks the APC/C–Cdc20 interaction. This work suggests that simultaneous disruption of multiple, weak protein–protein interactions is an effective approach for inactivating a protein machine.

Population Genomic Inferences from Sparse High-Throughput Sequencing of Two Populations of Drosophila melanogaster

Short-read sequencing techniques provide the opportunity to capture genome-wide sequence data in a single experiment. A current challenge is to identify questions that shallow-depth genomic data can address successfully and to develop corresponding analytical methods that are statistically sound. Here, we apply the Roche/454 platform to survey natural variation in strains of Drosophila melanogaster from an African (n = 3) and a North American (n = 6) population. Reads were aligned to the reference D. melanogaster genomic assembly, single nucleotide polymorphisms were identified, and nucleotide variation was quantified genome wide. Simulations and empirical results suggest that nucleotide diversity can be accurately estimated from sparse data with as little as 0.2× coverage per line. The unbiased genomic sampling provided by random short-read sequencing also allows insight into distributions of transposable elements and copy number polymorphisms found within populations and demonstrates that short-read sequencing methods provide an efficient means to quantify variation in genome organization and content. Continued development of methods for statistical inference of shallow-depth genome-wide sequencing data will allow such sparse, partial data sets to become the norm in the emerging field of population genomics.

Genotype and gene expression associations with immune function in Drosophila.

It is now well established that natural populations of Drosophila melanogaster harbor substantial genetic variation associated with physiological measures of immune function. In no case, however, have intermediate measures of immune function, such as transcriptional activity of immune-related genes, been tested as mediators of phenotypic variation in immunity. In this study, we measured bacterial load sustained after infection of D. melanogaster with Serratia marcescens, Providencia rettgeri, Enterococcus faecalis, and Lactococcus lactis in a panel of 94 third-chromosome substitution lines. We also measured transcriptional levels of 329 immune-related genes eight hours after infection with E. faecalis and S. marcescens in lines from the phenotypic tails of the test panel. We genotyped the substitution lines at 137 polymorphic markers distributed across 25 genes in order to test for statistical associations among genotype, bacterial load, and transcriptional dynamics. We find that genetic polymorphisms in the pathogen recognition genes (and particularly in PGRP-LC, GNBP1, and GNBP2) are most significantly associated with variation in bacterial load. We also find that overall transcriptional induction of effector proteins is a significant predictor of bacterial load after infection with E. faecalis, and that a marker upstream of the recognition gene PGRP-SD is statistically associated with variation in both bacterial load and transcriptional induction of effector proteins. These results show that polymorphism in genes near the top of the immune system signaling cascade can have a disproportionate effect on organismal phenotype due to the amplification of minor effects through the cascade.