Multilocus analysis of extracellular putative virulence proteins made by group A Streptococcus: Population genetics, human serologic response, and gene transcription

Species of pathogenic microbes are composed of an array of evolutionarily distinct chromosomal genotypes characterized by diversity in gene content and sequence (allelic variation). The occurrence of substantial genetic diversity has hindered progress in developing a comprehensive understanding of the molecular basis of virulence and new therapeutics such as vaccines. To provide new information that bears on these issues, 11 genes encoding extracellular proteins in the human bacterial pathogen group A Streptococcus identified by analysis of four genomes were studied. Eight of the 11 genes encode proteins with a LPXTG(L) motif that covalently links Gram-positive virulence factors to the bacterial cell surface. Sequence analysis of the 11 genes in 37 geographically and phylogenetically diverse group A Streptococcus strains cultured from patients with different infection types found that recent horizontal gene transfer has contributed substantially to chromosomal diversity. Regions of the inferred proteins likely to interact with the host were identified by molecular population genetic analysis, and Western immunoblot analysis with sera from infected patients confirmed that they were antigenic. Real-time reverse transcriptase–PCR (TaqMan) assays found that transcription of six of the 11 genes was substantially up-regulated in the stationary phase. In addition, transcription of many genes was influenced by the covR and mga trans-acting gene regulatory loci. Multilocus investigation of putative virulence genes by the integrated approach described herein provides an important strategy to aid microbial pathogenesis research and rapidly identify new targets for therapeutics research.

Population genetics of Ice Age brown bears

The Pleistocene was a dynamic period for Holarctic mammal species, complicated by episodes of glaciation, local extinctions, and intercontinental migration. The genetic consequences of these events are difficult to resolve from the study of present-day populations. To provide a direct view of population genetics in the late Pleistocene, we measured mitochondrial DNA sequence variation in seven permafrost-preserved brown bear (Ursus arctos) specimens, dated from 14,000 to 42,000 years ago. Approximately 36,000 years ago, the Beringian brown bear population had a higher genetic diversity than any extant North American population, but by 15,000 years ago genetic diversity appears similar to the modern day. The older, genetically diverse, Beringian population contained sequences from three clades now restricted to local regions within North America, indicating that current phylogeographic patterns may provide misleading data for evolutionary studies and conservation management. The late Pleistocene phylogeographic data also indicate possible colonization routes to areas south of the Cordilleran ice sheet.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.

Molecular markers reveal cryptic sex in the human pathogen Coccidioides immitis.

Coccidioides immitis, cause of a recent epidemic of "Valley fever" in California, is typical of many eukaryotic microbes in that mating and meiosis have yet to be reported, but it is not clear whether sex is truly absent or just cryptic. To find out, we have undertaken a population genetic study using PCR amplification, screening for single-strand conformation polymorphisms, and direct DNA sequencing to find molecular markers with nucleotide-level resolution. Both population genetic and phylogenetic analyses indicate that C. immitis is almost completely recombining. To our knowledge, this study is the first to find molecular evidence for recombination in a fungus for which no sexual stage has yet been described. These results motivate a directed search for mating and meiosis and illustrate the utility of single-strand conformation polymorphism and sequencing with arbitrary primer pairs in molecular population genetics.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

ALFRED: an allele frequency database for diverse populations and DNA polymorphisms—an update

ALFRED (the ALelle FREquency Database) is designed to store and disseminate frequencies of alleles at human polymorphic sites for multiple populations, primarily for the population genetics and molecular anthropology communities. Currently ALFRED has information on over 180 polymorphic sites for more than 70 populations. Since our initial release of the database we have focussed on increasing the quantity and quality of data, making reciprocal links between ALFRED and other related databases, and providing useful tools to make the data more comprehensible to the end user. ALFRED is accessible from the Kidd Lab home page (http://info.med.yale.edu/genetics/kkidd/) or from ALFRED directly (http://alfred.med.yale.edu/alfred/index.asp).

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd′ of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd′ suggests that the inversion may be the result of the transposable element’s activity. Characteristics of Odysseus’ terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

Concordance of gene genealogies reveals reproductive isolation in the pathogenic fungus Coccidioides immitis

Simple cladogenetic theory suggests that gene genealogies can be used to detect mixis in a population and delineate reproductively isolated groups within sexual taxa. We have taken this approach in a study of Coccidioides immitis, an ascomycete fungus responsible for a recent epidemic of coccidioidomycosis (Valley fever) in California. To test whether this fungus represents a single sexual species throughout its entire geographic range, we have compared genealogies from fragments of five nuclear genes. The five genealogies show multiple incompatibilities indicative of sex, but also share a branch that partitions the isolates into two reproductively isolated taxa, one centered in California and the other outside California. We conclude that coccidioidomycosis can be caused by two distinct noninterbreeding taxa. This result should aid the future study of the disease and illustrates the utility of the genealogical approach in population genetics.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

Social evolution in a new environment: the case of introduced fire ants.

The inadvertent introduction of the fire ant Solenopsis invicta to the United States from South America provides the opportunity to study recent social evolution by comparing social organization in native and introduced populations. We report that several important elements of social organization in multiple-queen nests differ consistently and dramatically between ants in Argentina and the United States. Colonies in Argentina contain relatively few queens and they are close relatives, whereas colonies in the United States contain high numbers of unrelated queens. A corollary of these differences is that workers in the native populations are significantly related to the new queens that they rear in contrast to the zero relatedness between workers and new queens in the introduced populations. The observed differences in queen number and relatedness signal a shift in the breeding biology of the introduced ants that is predicted on the basis of the high population densities in the new range. An additional difference in social organization that we observed, greater proportions of permanently unmated queens in introduced than in native populations, is predicted from the loss of alleles at the sex-determining locus and consequent skewing of operational sex ratios in the colonizing ants. Thus, significant recent social evolution in fire ants is consistent with theoretical expectations based on the altered ecology and population genetics of the introduced populations.

Inferring identify from DNA profile evidence.

The controversy over the interpretation of DNA profile evidence in forensic identification can be attributed in part to confusion over the mode(s) of statistical inference appropriate to this setting. Although there has been substantial discussion in the literature of, for example, the role of population genetics issues, few authors have made explicit the inferential framework which underpins their arguments. This lack of clarity has led both to unnecessary debates over ill-posed or inappropriate questions and to the neglect of some issues which can have important consequences. We argue that the mode of statistical inference which seems to underlie the arguments of some authors, based on a hypothesis testing framework, is not appropriate for forensic identification. We propose instead a logically coherent framework in which, for example, the roles both of the population genetics issues and of the nonscientific evidence in a case are incorporated. Our analysis highlights several widely held misconceptions in the DNA profiling debate. For example, the profile frequency is not directly relevant to forensic inference. Further, very small match probabilities may in some settings be consistent with acquittal. Although DNA evidence is typically very strong, our analysis of the coherent approach highlights situations which can arise in practice where alternative methods for assessing DNA evidence may be misleading.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.