Daniel L. Hartl

Mitochondrial pseudogenes: evolution's misplaced witnesses

Nuclear copies of mitochondrial DNA (mtDNA) have contaminated PCR-based mitochondrial studies of over 64 different animal species. Since the last review of these nuclear mitochondrial pseudogenes (Numts) in animals, Numts have been found in 53 of the species studied. The recent evidence suggests that Numts are not equally abundant in all species, for example they are more common in plants than in animals, and also more numerous in humans than in Drosophila. Methods for avoiding Numts have now been tested, and several recent studies demonstrate the potential utility of Numt DNA sequences in evolutionary studies. As relics of ancient mtDNA, these pseudogenes can be used to infer ancestral states or root mitochondrial phylogenies. Where they are numerous and selectively unconstrained, Numts are ideal for the study of spontaneous mutation in nuclear genomes.

A genome-wide view of the spectrum of spontaneous mutations in yeast

The mutation process ultimately defines the genetic features of all populations and, hence, has a bearing on a wide range of issues involving evolutionary genetics, inheritance, and genetic disorders, including the predisposition to cancer. Nevertheless, formidable technical barriers have constrained our understanding of the rate at which mutations arise and the molecular spectrum of their effects. Here, we report on the use of complete-genome sequencing in the characterization of spontaneously arising mutations in the yeast Saccharomyces cerevisiae. Our results confirm some findings previously obtained by indirect methods but also yield numerous unexpected findings, in particular a very high rate of point mutation and skewed distribution of base-substitution types in the mitochondrion, a very high rate of segmental duplication and deletion in the nuclear genome, and substantial deviations in the mutational profile among various model organisms.

Missense meanderings in sequence space: a biophysical view of protein evolution

Proteins are finicky molecules; they are barely stable and are prone to aggregate, but they must function in a crowded environment that is full of degradative enzymes bent on their destruction. It is no surprise that many common diseases are due to missense mutations that affect protein stability and aggregation. Here we review the literature on biophysics as it relates to molecular evolution, focusing on how protein stability and aggregation affect organismal fitness. We then advance a biophysical model of protein evolution that helps us to understand phenomena that range from the dynamics of molecular adaptation to the clock-like rate of protein evolution.

Selection for short introns in highly expressed genes.

Transcription is a slow and expensive process: in eukaryotes, approximately 20 nucleotides can be transcribed per second at the expense of at least two ATP molecules per nucleotide. Thus, at least for highly expressed genes, transcription of long introns, which are particularly common in mammals, is costly. Using data on the expression of genes that encode proteins in Caenorhabditis elegans and Homo sapiens, we show that introns in highly expressed genes are substantially shorter than those in genes that are expressed at low levels. This difference is greater in humans, such that introns are, on average, 14 times shorter in highly expressed genes than in genes with low expression, whereas in C. elegans the difference in intron length is only twofold. In contrast, the density of introns in a gene does not strongly depend on the level of gene expression. Thus, natural selection appears to favor short introns in highly expressed genes to minimize the cost of transcription and other molecular processes, such as splicing.

The cost of inbreeding in Arabidopsis

Population geneticists have long sought to estimate the distribution of selection intensities among genes of diverse function across the genome. Only recently have DNA sequencing and analytical techniques converged to make this possible. Important advances have come from comparing genetic variation within species (polymorphism) with fixed differences between species (divergence). These approaches have been used to examine individual genes for evidence of selection. Here we use the fact that the time since species divergence allows combination of data across genes. In a comparison of amino-acid replacements among species of the mustard weed Arabidopsis with those among species of the fruitfly Drosophila, we find evidence for predominantly beneficial gene substitutions in Drosophila but predominantly detrimental substitutions in Arabidopsis. We attribute this difference to the Arabidopsis mating system of partial self-fertilization, which corroborates a prediction of population genetics theory that species with a high frequency of inbreeding are less efficient in eliminating deleterious mutations owing to their reduced effective population size.

A genome-wide map of diversity in Plasmodium falciparum

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (π = 1.16 × 10−3) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.

GeneMerge—post-genomic analysis, data mining, and hypothesis testing

SUMMARY GeneMerge is a web-based and standalone program written in PERL that returns a range of functional and genomic data for a given set of study genes and provides statistical rank scores for over-representation of particular functions or categories in the data set. Functional or categorical data of all kinds can be analyzed with GeneMerge, facilitating regulatory and metabolic pathway analysis, tests of population genetic hypotheses, cross-experiment comparisons, and tests of chromosomal clustering, among others. GeneMerge can perform analyses on a wide variety of genomic data quickly and easily and facilitates both data mining and hypothesis testing. AVAILABILITY GeneMerge is available free of charge for academic use over the web and for download from: http://www.oeb.harvard.edu/hartl/lab/publications/GeneMerge.html.

Selective sweep of a newly evolved sperm-specific gene in Drosophila

The pattern of genetic variation across the genome of Drosophila melanogaster is consistent with the occurrence of frequent ‘selective sweeps’, in which new favourable mutations become incorporated into the species so quickly that linked alleles can ‘hitchhike’ and also become fixed. Because of the hitchhiking of linked genes, it is generally difficult to identify the target of any putative selective sweep. Here, however, we identify a new gene in D. melanogaster that codes for a sperm-specific axonemal dynein subunit. The gene has a new testes-specific promoter derived from a protein-coding region in a gene encoding the cell-adhesion protein annexin X (AnnX), and it contains a new protein-coding exon derived from an intron in a gene encoding a cytoplasmic dynein intermediate chain (Cdic). The new transcription unit, designated Sdic (for sperm-specific dynein intermediate chain), has been duplicated about tenfold in a tandem array. Consistent with the selective sweep of this gene, the level of genetic polymorphism near Sdic is unusually low. The discovery of this gene supports other results that point to the rapid molecular evolution of male reproductive functions.

A single mode of canalization

The evolution of mechanisms underlying the buffering of the phenotype against genetic and environmental influences has received much theoretical and experimental attention, yet many issues remain unresolved. Here, we consider the kinds of biological process that are likely to promote this buffering, or canalization, and the circumstances under which the evolution of these mechanisms will be favored. We conclude that evolution should produce a single mode of canalization that will buffer the phenotype against all kinds of perturbation, and that the major fitness benefit driving the fixation of canalizing alleles derives from a reduction in environmental influences on phenotypic variation.

RADseq underestimates diversity and introduces genealogical biases due to nonrandom haplotype sampling

Reduced representation genome‐sequencing approaches based on restriction digestion are enabling large‐scale marker generation and facilitating genomic studies in a wide range of model and nonmodel systems. However, sampling chromosomes based on restriction digestion may introduce a bias in allele frequency estimation due to polymorphisms in restriction sites. To explore the effects of this nonrandom sampling and its sensitivity to different evolutionary parameters, we developed a coalescent‐simulation framework to mimic the biased recovery of chromosomes in restriction‐based short‐read sequencing experiments (RADseq). We analysed simulated DNA sequence datasets and compared known values from simulations with those that would be estimated using a RADseq approach from the same samples. We compare these ‘true’ and ‘estimated’ values of commonly used summary statistics, π, θw, Tajimas D and FST. We show that loci with missing haplotypes have estimated summary statistic values that can deviate dramatically from true values and are also enriched for particular genealogical histories. These biases are sensitive to nonequilibrium demography, such as bottlenecks and population expansion. In silico digests with 102 completely sequenced Drosophila melanogaster genomes yielded results similar to our findings from coalescent simulations. Though the potential of RADseq for marker discovery and trait mapping in nonmodel systems remains undisputed, our results urge caution when applying this technique to make population genetic inferences.