Patterns of variance in stage-structured populations: Evolutionary predictions and ecological implications

Variability in population growth rate is thought to have negative consequences for organism fitness. Theory for matrix population models predicts that variance in population growth rate should be the sum of the variance in each matrix entry times the squared sensitivity term for that matrix entry. I analyzed the stage-specific demography of 30 field populations from 17 published studies for pattern between the variance of a demographic term and its contribution to population growth. There were no instances in which a matrix entry both was highly variable and had a large effect on population growth rate; instead, correlations between estimates of temporal variance in a term and contribution to population growth (sensitivity or elasticity) were overwhelmingly negative. In addition, survivorship or growth sensitivities or elasticities always exceeded those of fecundity, implying that the former two terms always contributed more to population growth rate. These results suggest that variable life history stages tend to contribute relatively little to population growth rates because natural selection may alter life histories to minimize stages with both high sensitivity and high variation.

Evolutionary EST analysis identifies rapidly evolving male reproductive proteins in Drosophila

Sequence comparisons of genomes or expressed sequence tags (ESTs) from related organisms provide insight into functional conservation and diversification. We compare the sequences of ESTs from the male accessory gland of Drosophila simulans to their orthologs in its close relative Drosophila melanogaster, and demonstrate rapid divergence of many of these reproductive genes. Nineteen (∼11%) of 176 independent genes identified in the EST screen contain protein-coding regions with an excess of nonsynonymous over synonymous changes, suggesting that their divergence has been accelerated by positive Darwinian selection. Genes that encode putative accessory gland-specific seminal fluid proteins had a significantly elevated level of nonsynonymous substitution relative to nonaccessory gland-specific genes. With the 57 new accessory gland genes reported here, we predict that ∼90% of the male accessory gland genes have been identified. The evolutionary EST approach applied here to identify putative targets of adaptive evolution is readily applicable to other tissues and organisms.

Evolutionary genomics of Staphylococcus aureus: Insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic

An emerging theme in medical microbiology is that extensive variation exists in gene content among strains of many pathogenic bacterial species. However, this topic has not been investigated on a genome scale with strains recovered from patients with well-defined clinical conditions. Staphylococcus aureus is a major human pathogen and also causes economically important infections in cows and sheep. A DNA microarray representing >90% of the S. aureus genome was used to characterize genomic diversity, evolutionary relationships, and virulence gene distribution among 36 strains of divergent clonal lineages, including methicillin-resistant strains and organisms causing toxic shock syndrome. Genetic variation in S. aureus is very extensive, with ≈22% of the genome comprised of dispensable genetic material. Eighteen large regions of difference were identified, and 10 of these regions have genes that encode putative virulence factors or proteins mediating antibiotic resistance. We find that lateral gene transfer has played a fundamental role in the evolution of S. aureus. The mec gene has been horizontally transferred into distinct S. aureus chromosomal backgrounds at least five times, demonstrating that methicillin-resistant strains have evolved multiple independent times, rather than from a single ancestral strain. This finding resolves a long-standing controversy in S. aureus research. The epidemic of toxic shock syndrome that occurred in the 1970s was caused by a change in the host environment, rather than rapid geographic dissemination of a new hypervirulent strain. DNA microarray analysis of large samples of clinically characterized strains provides broad insights into evolution, pathogenesis, and disease emergence.

Evolutionary reconstruction of a hairpin deleted from the genome of an RNA virus.

The intercistronic region between the maturation and coat-protein genes of RNA phage MS2 contains important regulatory and structural information. The sequence participates in two adjacent stem-loop structures, one of which, the coat-initiator hairpin, controls coat-gene translation and is thus under strong selection pressure. We have removed 19 out of the 23 nucleotides constituting the intercistronic region, thereby destroying the capacity of the phage to build the two hairpins. The deletion lowered coat-protein yield more than 1000-fold, and the titer of the infectious clone carrying the deletion dropped 10 orders of magnitude as compared with the wild type. Two types of revertants were recovered. One had, in two steps, recruited 18 new nucleotides that served to rebuild the two hairpins and the lost Shine-Dalgarno sequence. The other type had deleted an additional six nucleotides, which allowed the reconstruction of the Shine-Dalgarno sequence and the initiator hairpin, albeit by sacrificing the remnants of the other stem-loop. The results visualize the immense genetic repertoire created by, what appears as, random RNA recombination. It would seem that in this genetic ensemble every possible new RNA combination is represented.

Two new forms of gonadotropin-releasing hormone in a protochordate and the evolutionary implications.

The neuropeptide gonadotropin-releasing hormone (GnRH) is the major regulator of reproduction in vertebrates. Our goal was to determine whether GnRH could be isolated and identified by primary structure in a protochordate and to examine its location by immunocytochemistry. The primary structure of two novel decapeptides from the tunicate Chelyosoma productum (class Ascidiacea) was determined. Both show significant identity with vertebrate GnRH. Tunicate GnRH-I (pGlu-His-Trp-Ser-Asp-Tyr-Phe-Lys-Pro-Gly-NH2) has 60% of its residues conserved, compared with mammalian GnRH, whereas tunicate GnRH-II (pGlu-His-Trp-Ser-Leu-Cys-His-Ala-Pro-Gly-NH2) is unusual in that it was isolated as a disulfide-linked dimer. Numerous immunoreactive GnRH neurons lie within blood sinuses close to the gonoducts and gonads in both juveniles and adults, implying that the neuropeptide is released into the bloodstream. It is suggested that in ancestral chordates, before the evolution of the pituitary, the hormone was released into the bloodstream and acted directly on the gonads.

DNA sequence insertion and evolutionary variation in gene regulation.

Current evidence on the long-term evolutionary effect of insertion of sequence elements into gene regions is reviewed, restricted to cases where a sequence derived from a past insertion participates in the regulation of expression of a useful gene. Ten such examples in eukaryotes demonstrate that segments of repetitive DNA or mobile elements have been inserted in the past in gene regions, have been preserved, sometimes modified by selection, and now affect control of transcription of the adjacent gene. Included are only examples in which transcription control was modified by the insert. Several cases in which merely transcription initiation occurred in the insert were set aside. Two of the examples involved the long terminal repeats of mammalian endogenous retroviruses. Another two examples were control of transcription by repeated sequence inserts in sea urchin genomes. There are now six published examples in which Alu sequences were inserted long ago into human gene regions, were modified, and now are central in control/enhancement of transcription. The number of published examples of Alu sequences affecting gene control has grown threefold in the last year and is likely to continue growing. Taken together, all of these examples show that the insertion of sequence elements in the genome has been a significant source of regulatory variation in evolution.

A common evolutionary origin for mitochondria and hydrogenosomes.

Trichomonads are among the earliest eukaryotes to diverge from the main line of eukaryotic descent. Keeping with their ancient nature, these facultative anaerobic protists lack two "hallmark" organelles found in most eukaryotes: mitochondria and peroxisomes. Trichomonads do, however, contain an unusual organelle involved in carbohydrate metabolism called the hydrogenosome. Like mitochondria, hydrogenosomes are double-membrane bounded organelles that produce ATP using pyruvate as the primary substrate. Hydrogenosomes are, however, markedly different from mitochondria as they lack DNA, cytochromes and the citric acid cycle. Instead, they contain enzymes typically found in anaerobic bacteria and are capable of producing molecular hydrogen. We show here that hydrogenosomes contain heat shock proteins, Hsp70, Hsp60, and Hsp10, with signature sequences that are conserved only in mitochondrial and alpha-Gram-negative purple bacterial Hsps. Biochemical analysis of hydrogenosomal Hsp60 shows that the mature protein isolated from the organelle lacks a short, N-terminal sequence, similar to that observed for most nuclear-encoded mitochondrial matrix proteins. Moreover, phylogenetic analyses of hydrogenosomal Hsp70, Hsp60, and Hsp10 show that these proteins branch within a monophyletic group composed exclusively of mitochondrial homologues. These data establish that mitochondria and hydrogenosomes have a common eubacterial ancestor and imply that the earliest-branching eukaryotes contained the endosymbiont that gave rise to mitochondria in higher eukaryotes.

Evolutionary conservation of sequence elements controlling cytoplasmic polyadenylylation.

Cytoplasmic polyadenylylation is an evolutionarily conserved mechanism involved in the translational activation of a set of maternal messenger RNAs (mRNAs) during early development. In this report, we show by interspecies injections that Xenopus and mouse use the same regulatory sequences to control cytoplasmic poly(A) addition during meiotic maturation. Similarly, Xenopus and Drosophila embryos exploit functionally conserved signals to regulate polyadenylylation during early post-fertilization development. These experiments demonstrate that the sequence elements that govern cytoplasmic polyadenylylation, and hence one form of translational activation, function across species. We infer that the requisite regulatory sequence elements, and likely the trans-acting components with which they interact, have been conserved since the divergence of vertebrates and arthropods.

Evolutionary changes in the expression pattern of a developmentally essential gene in three Drosophila species.

The hypothesis that morphological evolution may largely result from changes in gene regulation rather than gene structure has been difficult to test. Morphological differences among insects are often apparent in the cuticle structures produced. The dopa decarboxylase (Ddc) and alpha-methyldopa hypersensitive (amd) genes arose from an ancient gene duplication. In Drosophila, they have evolved nonoverlapping functions, including the production of distinct types of cuticle, and for Ddc, the production of the neurotransmitters, dopamine and serotonin. The amd gene is particularly active in the production of specialized flexible cuticles in the developing embryo. We have compared the pattern of amd expression in three Drosophila species. Several regions of expression conserved in all three species but, surprisingly, a unique domain of expression is found in Drosophila simulans that does occur in the closely related (2-5 million years) Drosophila melanogaster or in the more remote (40-50 million years) Drosophila virilis. The "sudden" appearance of a completely new and robust domain of expression provides a glimpse of evolutionary variation resulting from changes in regulation of structural gene expression.

Microsatellite spreading in the human genome: evolutionary mechanisms and structural implications.

Microsatellites are tandem repeat sequences abundant in the genomes of higher eukaryotes and hitherto considered as "junk DNA." Analysis of a human genome representative data base (2.84 Mb) reveals a distinct juxtaposition of A-rich microsatellites and retroposons and suggests their coevolution. The analysis implies that most microsatellites were generated by a 3'-extension of retrotranscripts, similar to mRNA polyadenylylation, and that they serve in turn as "retroposition navigators," directing the retroposons via homology-driven integration into defined sites. Thus, they became instrumental in the preservation and extension of primordial genomic patterns. A role is assigned to these reiterating A-rich loci in the higher-order organization of the chromatin. The disease-associated triplet repeats are mostly found in coding regions and do not show an association with retroposons, constituting a unique set within the family of microsatellite sequences.

Reduction of two functional gamma-globin genes to one: an evolutionary trend in New World monkeys (infraorder Platyrrhini).

Nucleotide sequences were determined for the gamma1- and gamma2-globin loci from representatives of the seven anciently separated clades in the three extant platyrrhine families (Atelidae, Pitheciidae, and Cebidae). These sequences revealed an evolutionary trend in New World monkeys either to inactivate the gamma1 gene or to fuse it with the gamma2 gene, i.e. to have only one functional fetally expressed gamma gene. This trend is clearly evident in six of the seven clades: (i) it occurred in atelids by deletion of most of the gamma1 gene in the basal ancestor of this clade; (ii-iv) in pitheciid titi, saki, and cebid capuchin monkeys by potentially debilitating nucleotide substitutions in the proximal CCAAT box of the gamma1 promoters and (v and vi) in cebid owl and squirrel monkeys by crossovers that fused 5' sequence from gamma1 with 3' sequence from gamma2. In the five clades with gamma1 and gamma2 loci separated by intergenic sequences (the fifth clade being the cebid marmosets), the gamma2 genes retained an unaltered proximal CCAAT motif and their gamma2 promoters accumulated fewer nucleotide substitutions than did the gamma1 promoters. Thus, phylogenetic considerations indicate that the stem platyrrhines, ancestral to all New World monkeys, had gamma2 as the primary fetally expressed gamma gene. A further inference is that when the earlier stem anthropoid gamma gene duplicated, gamma2 (at its greater downstream distance from epsilon) could evade embryonic activation by the locus control region but could be fetally activated once released by regulatory mutations from fetal repressors.

Evolutionary relationships of the coelacanth, lungfishes, and tetrapods based on the 28S ribosomal RNA gene.

The origin of land vertebrates was one of the major transitions in the history of vertebrates. Yet, despite many studies that are based on either morphology or molecules, the phylogenetic relationships among tetrapods and the other two living groups of lobe-finned fishes, the coelacanth and the lungfishes, are still unresolved and debated. Knowledge of the relationships among these lineages, which originated back in the Devonian, has profound implications for the reconstruction of the evolutionary scenario of the conquest of land. We collected the largest molecular data set on this issue so far, about 3,500 base pairs from seven species of the large 28S nuclear ribosomal gene. All phylogenetic analyses (maximum parsimony, neighbor-joining, and maximum likelihood) point toward the hypothesis that lungfishes and coelacanths form a monophyletic group and are equally closely related to land vertebrates. This evolutionary hypothesis complicates the identification of morphological or physiological preadaptations that might have permitted the common ancestor of tetrapods to colonize land. This is because the reconstruction of its ancestral conditions would be hindered by the difficulty to separate uniquely derived characters from shared derived characters in the coelacanth/lungfish and tetrapod lineages. This molecular phylogeny aids in the reconstruction of morphological evolutionary steps by providing a framework; however, only paleontological evidence can determine the sequence of morphological acquisitions that allowed lobe-finned fishes to colonize land.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.