20: NASOPHARYNX MICROBIOME COMPOSITION VARIES OVER TIME IN PEDIATRIC ASTHMA

Purpose of Study

The application of next-generation sequencing (NGS) technology has shown that microbial communities in the respiratory airways (i.e., the microbiome) play a significant role in the onset, development and severity of asthma. However, little is known about their temporal dynamics (i.e., microbial succession), which poses a significant obstacle to identifying pulmotypes of disease and assessing inter-patient variation. Here, we couple NGS and 16S rRNA data to characterize the nasopharynx microbiome of children with asthma and determine its stability over time.

Methods Used

We collected nasal washes from 40 children with asthma enrolled in the AsthMaP-2 Project from two consecutive visits, six months apart. Total DNA was extracted and sequenced for the 16S-V4 rRNA gene region (∼250 bp) using the MySeq Illumina platform. Reads were analyzed in Mothur using the SILVAv119 reference database. Alpha diversity metrics and phylogenetic and count-base distance community indexes of beta diversity were compared across samples and time points. PCoA and NJ clustering analysis were used to assess community relatedness. Differences in alpha diversity and OTU abundance between sample pairs across time points were also compared.

Summary of Results

A mean of 27,479 clean 16S sequences corresponding to an average of 173 OTUs were sequenced and detected per sample, respectively. Representatives of Moraxella, Corynebacterium, Prevotella, Staphylococcus, Alloiococcus, Streptococcus, Peptoniphilus, Fusobacterium, and Haemophilus accounted for 36 to 99% of the reads across samples. These genera have been previously found in the nasopharynx of asthmatic and healthy children. A total of 61 OTUs from these genera were present in at least 50% of the samples (i.e., the nasal core microbiome). Significant differences in core microbiome composition were detected between sample pairs, but no directional trend (increase or decrease) was observed across sample pairs. Samples were randomly ordinated and did not cluster together.

Conclusions

Our analysis of nasal microbiomes in 40 asthmatic children revealed significant differences in composition within individuals over six months. Future cross-sectional microbiome studies need to be aware of short span temporal dynamics in nasal microbiota.

Conservation assessment of southern South American freshwater ecoregions on the basis of the distribution and genetic diversity of crabs from the genus Aegla

We assessed the conservation priority of 18 freshwater ecoregions in southern South America on the basis of Aegla (genus of freshwater crabs) genetic diversity and distribution. Geographical distributions for 66 Aegla species were taken from the literature and plotted against ecoregions and main river basins of southern South America. Species richness and number of threatened and endemic species were calculated for each area. To assess taxonomic and phylogenetic diversity, we generated a molecular phylogeny based on DNA sequences for one nuclear (28S) and 4 mitochondrial (12S, 16S, COI, and COII) genes. All species richness and phylogenetic methods agreed, to a large extent, in their rankings of the importance of conservation areas, as indicated by the Spearman's rank correlation coefficient (p < 0.01); nonetheless, some of the lowest correlations were observed between taxonomic and phylogenetic diversity indices. The 5 ecoregions of the Laguna dos Patos Basin (Eastern Brazil), Central Chile, South Brazilian Coast, Chilean Lakes, and Subtropical Potamic Axis (northern Argentina and southern Uruguay and Paraguay) had the highest biodiversity scores. Conservation of these regions will preserve the largest number of species and the greatest amount of genetic diversity within the South American freshwater Aegla fauna. Biodiversity across rivers and within areas was heterogeneously distributed in the ecoregions of Upper Paraná, Ribeira do Iguape, Upper Uruguay, and South Brazilian Coast (i.e., one river showed significantly more biodiversity than any other river from the same ecoregion), but homogeneously distributed in the other ecoregions. Hence, conservation plans in the former regions will potentially require less effort than plans in the latter regions.

A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints

We have developed a modified BOOTSCAN algorithm that may be used to screen nucleotide sequence alignments for evidence of recombination without prior identification of nonrecombinant reference sequences. The algorithm is fast and includes a Bonferroni corrected statistical test of recombination to circumvent the multiple testing problems encountered when using the BOOTSCAN method to explore alignments for evidence of recombination. Using both simulated and real datasets we demonstrate that the modified algorithm is more powerful than other phylogenetic recombination detection methods and performs almost as well as one of the best substitution distribution recombination detection methods.

Phylogeography and speciation of colour morphs in the colonial ascidian Pseudodistoma crucigaster

Variation in pigmentation is common in marine invertebrates, although few studies have shown the existence of genetic differentiation of chromatic forms in these organisms. We studied the genetic structure of a colonial ascidian with populations of different colour morphs in the northwestern Mediterranean. A fragment of the c oxidase subunit 1 (COI) mitochondrial gene was sequenced in seven populations of Pseudodistoma crucigaster belonging to three different colour morphs (orange, yellow and grey). Maximum likelihood analyses showed two well-supported clades separating the orange morph from the yellow-grey morphotypes. Genetic divergence between these clades was 2.12%, and gamma(ST) values between populations of the two clades were high (average 0.936), pointing to genetic isolation. Nested clade and coalescence analyses suggest that a past fragmentation event may explain the phylogeographical origin of these two clades. Non-neutral mtDNA evolution is observed in our data when comparing the two clades, showing a significant excess of nonsynonymous polymorphism within the yellow-grey morphotype using the McDonald-Kreitman test, which is interpreted as further support of reproductive isolation. We conclude that the two clades might represent separate species. We compare the population genetic differentiation found with that estimated for other colonial and solitary ascidian species, and relate it to larval dispersal capabilities and other life-history traits.

Evaluation of methods for detecting recombination from DNA sequences: computer simulations

Recombination is a key evolutionary process that shapes the architecture of genomes and the genetic structure of populations. Although many statistical methods are available for the detection of recombination from DNA sequences, their absolute and relative performance is still unknown. Here we evaluated the performance of 14 different recombination detection algorithms. We used the coalescent with recombination to simulate DNA sequences with different levels of recombination, genetic diversity, and rate variation among sites. Recombination detection methods were applied to these data sets, and whether they detected or not recombination was recorded. Different recombination methods showed distinct performance depending on the amount of recombination, genetic diversity, and rate variation among sites. The model of nucleotide substitution under which the data were generated did not seem to have a significant effect. Most methods increase power with more sequence divergence. In general, recombination detection methods seem to capture the presence of recombination, but they are not very powerful. Methods that use substitution patterns or incompatibility among sites were more powerful than methods based on phylogenetic incongruence. Most methods do not seem to infer more false positives than expected by chance. Especially depending on the amount of diversity in the data, different methods could be used to attain maximum power while minimizing false positives. Results shown here will provide some guidance in the selection of the most appropriate method/s for the analysis of the particular data at hand.

Selecting the best-fit model of nucleotide substitution

Despite the relevant role of models of nucleotide substitution in phylogenetics, choosing among different models remains a problem. Several statistical methods for selecting the model that best fits the data at hand have been proposed, but their absolute and relative performance has not yet been characterized. In this study, we compare under various conditions the performance of different hierarchical and dynamic likelihood ratio tests, and of Akaike and Bayesian information methods, for selecting best-fit models of nucleotide substitution. We specifically examine the role of the topology used to estimate the likelihood of the different models and the importance of the order in which hypotheses are tested. We do this by simulating DNA sequences under a known model of nucleotide substitution and recording how often this true model is recovered by the different methods. Our results suggest that model selection is reasonably accurate and indicate that some likelihood ratio test methods perform overall better than the Akaike or Bayesian information criteria. The tree used to estimate the likelihood scores does not influence model selection unless it is a randomly chosen tree. The order in which hypotheses are tested, and the complexity of the initial model in the sequence of tests, influence model selection in some cases. Model fitting in phylogenetics has been suggested for many years, yet many authors still arbitrarily choose their models, often using the default models implemented in standard computer programs for phylogenetic estimation. We show here that a best-fit model can be readily identified. Consequently, given the relevance of models, model fitting should be routine in any phylogenetic analysis that uses models of evolution.

Selecting models of nucleotide substitution: an application to human immunodeficiency virus 1 (HIV-1)

The blind use of models of nucleotide substitution in evolutionary analyses is a common practice in the viral community. Typically, a simple model of evolution like the Kimura two-parameter model is used for estimating genetic distances and phylogenies, either because other authors have used it or because it is the default in various phylogenetic packages. Using two statistical approaches to model fitting, hierarchical likelihood ratio tests and the Akaike information criterion, we show that different viral data sets are better explained by different models of evolution. We demonstrate our results with the analysis of HIV-1 sequences from a hierarchy of samples; sequences within individuals, individuals within subtypes, and subtypes within groups. We also examine results for three different gene regions: gag, pol, and env. The Kimura two-parameter model was not selected as the best-fit model for any of these data sets, despite its widespread use in phylogenetic analyses of HIV-1 sequences. Furthermore, the model complexity increased with increasing sequence divergence. Finally, the molecular-clock hypothesis was rejected in most of the data sets analyzed, throwing into question clock-based estimates of divergence times for HIV-1. The importance of models in evolutionary analyses and their repercussions on the derived conclusions are discussed.

Human immunodeficiency virus type 1 quasi species that rebound after discontinuation of highly active antiretroviral therapy are similar to the viral quasi species present before initiation of therapy

In an effort to identify the sources of the viruses that emerge after discontinuation of therapy, analyses of human immunodeficiency virus (HIV) quasi species were done for 3 patients with sustained levels of HIV RNA of <50 copies/mL for 1-3 years. The sequences found in the rebounding plasma virus were closely related to those of the actively replicating form of viruses present before the initiation of combination therapy. All quasi species found in the rebounding plasma virus were also present in proviral DNA, cell-associated RNA in peripheral blood mononuclear cells (PBMC), and virion RNA derived from PBMC coculture during periods when plasma HIV RNA levels were <50 copies/mL. These findings suggest that the rapid resurgence of plasma viremia observed after discontinuation of therapy and the viruses cocultured from PBMC are derived from a relatively stable pool of the replicating form of virus rather than from activation of a previously latent pool.

Limited heterogeneity of HIV type 1 in infected mothers correlates with lack of vertical transmission

The human immunodeficiency virus type 1 (HIV-1) envelope V3 region sequences of peripheral blood mononuclear cell DNA were analyzed from three nontransmitting mothers (infected mothers who failed to transmit HIV-1 to their infants in the absence of antiretroviral therapy), including one mother with two deliveries, and compared with the sequences of seven previously analyzed transmitting mothers. The coding potential of the envelope open reading frame, including several patient-specific amino acid motifs and previously described molecular features across the V3 region, were highly conserved. There was a low degree of heterogeneity within the sequences of each nontransmitting mother compared with the sequences of transmitting mothers. In addition, the estimates of genetic diversity of nontransmitting mother sequences were significantly lower compared with transmitting mother sequences. Phylogenetic analysis showed that the sequences of each nontransmitting mother formed distinct clusters that were well discriminated from each other and the sequences of seven transmitting mothers. In conclusion, a low degree of HIV-1 genetic heterogeneity in these infected mothers correlates with lack of vertical transmission; this finding may be useful in developing strategies for further prevention of maternal-fetal transmission.

The monophyletic origin of freshwater crayfish estimated from nuclear and mitochondrial DNA sequences

Despite their widespread use as model organisms, the phylogenetic status of the around 520 species of freshwater crayfish is still in doubt. One hypothesis suggests two distinct origins of freshwater crayfish as indicated by their geographical distribution, with two centres of origin near the two present centres of diversity; one in south-eastern United States and the other in Victoria, Australia. An alternative theory proposes a single (monophyletic) origin of freshwater crayfish. Here we use over 3000 nucleotides from three different gene regions in estimating phylogenetic relationships among freshwater crayfish and related Crustacea. We show clear evidence for monophyly of freshwater crayfish and for the sister-group relationship between crayfish and clawed lobsters. Monophyly of the superfamilies Astacoidea and Parastacoidea is also supported. However, the monophyly of the family Cambaridae is questioned with the genus Cambaroides being associated with the Astacidae.

HTLV type I/II in British Columbia Amerindians: a seroprevalence study and sequence characterization of an HTLV type IIa isolate

It has been established that the human T cell lymphotropic viruses type I and II (HTLV-I and HTLV-II) are both present in some indigenous peoples of the Americas. While HTLV-I has been identified in coastal British Columbia Indians (BCIs), HTLV-II has not been previously reported in the BCIs or other Canadian Amerindians. The prevalence of HTLV-I and HTLV-II in these populations has not been extensively studied. In this article, we examine a group of BCIs from Vancouver Island who belong to the Nuu-Chah-Nulth and are known to have an increased incidence of rheumatic disease. In 494 serum samples from this tribe, the levels of prevalence of HTLV-I and HTLV-II were 2.8 and 1.6%, respectively. No association could be made between arthropathy and HTLV-I infection. In addition, we characterized an HTLV-II isolate of a BCI from the coastal mainland of British Columbia and with a history of intravenous drug abuse. This case represents the first molecular characterization of a Canadian Amerindian HTLV-II isolate: a subtype IIa virus with phylogenetic affinity for intravenous drug user isolates and containing an extended form of the Tax protein. These results are consistent either with this strain having been sampled from a polymorphic ancestral pool of HTLV-II that gave rise to the current epidemic spread of this virus by intravenous drug use and sexual transmission, or with its being "back-transmitted" into the BC Amerindian population in association with intravenous drug use.

PCR detection of host and HIV-1 sequences from archival brain tissue

Mutations in CCR5 and CCR2b have been recently shown to affect disease progression towards AIDS. A role for these host genotypes in AIDS dementia complex (ADC) has also been postulated but remains unclear. Additionally, brain-derived envelope sequences from HIV-1 have been associated with ADC but their specific contribution to pathogenesis remains uncertain. This study demonstrates the successful use of PCR techniques to isolate host CCR5 and CCR2b, and HIV-1 V3 sequences from paraffin embedded tissues from patients with and without ADC. PCR amplification from archival tissue offers a novel approach for studying the interactions between potential neuroprotective elements in the host and virulence determinants in HIV that may contribute to differences in susceptibility to ADC.

Population genetics of the porB gene of Neisseria gonorrhoeae: different dynamics in different homology groups

The porB locus codes for the major outer membrane protein of Neisseria gonorrhoeae. Alleles of this locus have been assigned to two homology groups based on close sequence and immunological relationships and are designated as either PIA or PIB. Several population parameters were estimated and compared among these two groups using a data set of 22 PIA sequences and 91 PIB sequences obtained from diverse geographic localities and from time periods spanning approximately 50 years. Recombination appears to be extensive in the porB gene. While the recombination rates are similar for the PIA and PIB sequences, the relative contribution of recombination to genetic diversity is higher for the PIA sequences. Alleles belonging to the PIB group show greater genetic diversity than do those in the PIA group. Although phylogenetic analysis did not reveal temporal or geographic clustering of sequences, estimates of gene flow and the fixation index suggested that PIB sequences exhibit population substructure based on geographic locality. Selection acts in these homology groups in a different way. While positive Darwinian selection is the dominant force driving the evolution of the PIA sequences, purifying selection operates also on the PIB sequences. These differences may be attributable to the greater propensity of PIA strains, as compared with PIB strains, to cause disseminated gonococcal infection, which would expose the former to intense selection pressure from the host immune system. The molecular evolution of Neisseria gonorrhoeae seems to be driven by the simultaneous action of selection and recombination, but under different rates and selection pressures for the PIA and PIB homology groups.

Intragenomic variation within ITS1 and ITS2 of freshwater crayfishes (Decapoda: Cambaridae): implications for phylogenetic and microsatellite studies

Intragenomic variation in ITS1 and ITS2 is known to exit but is widely ignored in phylogenetic studies using these gene regions. The amount of variation in seven crayfish species, including three populations of Orconectes luteus and two of Procambarus clarkii, was assessed by sequencing 3, 5, or 10 clones from the same individuals, for a total of 77 sequences. The ITS1 and ITS2 sequences reported here are some of the longest known, with aligned lengths of 760 and 1,300 bp, respectively. They contain multiple microsatellite insertions, all of which show considerable intragenomic variation in the number of repeat elements. This variation is enough to obscure phylogenetic relationships at the population level, although relationships between species can be estimated. Given the hybridization techniques used to locate microsatellites, multiple-copy regions like ITS1 and ITS2 will be preferentially found if they contain microsatellites, and in these cases the microsatellites will not behave as typical Mendelian markers and could give spurious results.

Mitochondrial DNA phylogeography and population history of the grey wolf canis lupus

The grey wolf (Canis lupus) and coyote (C. latrans) are highly mobile carnivores that disperse over great distances in search of territories and mates. Previous genetic studies have shown little geographical structure in either species. However, population genetic structure is also influenced by past isolation events and population fluctuations during glacial periods. In this study, control region sequence data from a worldwide sample of grey wolves and a more limited sample of coyotes were analysed. The results suggest that fluctuating population sizes during the late Pleistocene have left a genetic signature on levels of variation in both species. Genealogical measures of nucleotide diversity suggest that historical population sizes were much larger in both species and grey wolves were more numerous than coyotes. Currently, about 300 000 wolves and 7 million coyotes exist. In grey wolves, genetic diversity is greater than that predicted from census population size, reflecting recent historical population declines. By contrast, nucleotide diversity in coyotes is smaller than that predicted by census population size, reflecting a recent population expansion following the extirpation of wolves from much of North America. Both species show little partitioning of haplotypes on continental or regional scales. However, a statistical parsimony analysis indicates local genetic structure that suggests recent restricted gene flow.

Different models, different trees: the geographic origin of PTLV-I

Using nucleotide sequences from three genomic regions of the human and simian T-cell lymphotropic virus type I (HTLV-I/STLV-I)-consisting of 69 sequences from a 140-bp segment of the pol region, 98 sequences from a 503-bp segment of the LTR, and 154 sequences from a 386-bp segment of the env region-we tested two hypotheses concerning the geographic origin and evolution of STLV-I and HTLV-I. First, we tested the assumption of equal rates of evolution along STLV-I and HTLV-I lineages using a likelihood ratio test to ascertain whether current levels of genomic diversity can be used to determine ancestry. We demonstrated that unequal rates of evolution along HTLV-I and STLV-I lineages have occurred throughout evolutionary time, thus calling into question the use of pairwise distances to assign ancestry. Second, we constructed phylogenetic trees using multiple phylogenetic techniques to test for the geographic origin of STLV-I and HTLV-I. Using the principle of likelihood, we chose a statistically justified model of evolution for each data set. We demonstrated the utility of the likelihood ratio test to determine which model of evolution should be chosen for phylogenetic analyses, revealing that using different models of evolution produces conflicting results, and neither the hypothesis of an African origin nor the hypothesis of an Asian origin can be rejected statistically. Our best estimates of phylogenetic relationships, however, support an African origin of PTLV for each gene region.

Independent evolution of HIV type 1 in different brain regions

HIV-1-associated brain pathology exhibits regional variability and we therefore studied the genetic differences in the V1-V5 domains of the HIV env gene in up to four regions of brain (frontal lobe, basal ganglia, medial temporal lobe, and nonmedial temporal lobe) from three patients. We found that in each separate brain region HIV-1 forms different quasispecies and that there is little gene flow among these regions. In further support of brain region-specific evolution of HIV-1, we analyzed amino acid signatures in these clones. In addition to known amino acid signatures associated with macrophage tropism and the lack of syncytium formation, we found 15 majority amino acid signature patterns from the V1-V5 env sequences associated with the neuroanatomical regions analyzed from the three individuals. Furthermore, on average, intrabrain genetic distances for the HIV-1 env were estimated to be much smaller than genetic distances between brain regions. Specific strains of HIV-1 may be neurotropic or neuroinvasive (replication preference in brain tissue) and may contribute to pathology, cognitive loss, and neuropsychiatric disease.

Parallel evolution of drug resistance in HIV: failure of nonsynonymous/synonymous substitution rate ratio to detect selection

Parallel or convergent evolution at the molecular level has been difficult to demonstrate especially when rigorous statistical criteria are applied. We present sequence data from the protease gene from eight patients infected with the human immunodeficiency virus (HIV-1). These patients have been on multiple drug therapies for at least 2 years. We present sequence data from two timepoints: time zero--the initiation of drug therapy--and a subsequent timepoint between 59 and 104 weeks after the initiation of drug therapy. In addition to the sequence data, we present viral load data from both initial and final timepoints. Our phylogenetic analyses indicate significant evolution of virus from initial to final time points, even in three of eight patients who show low viral loads. Of the five patients who escaped drug therapy, identical amino acid replacements were seen in all five patients at two different codon positions, an indication of parallel evolution. We also measured genetic diversity for these patients and found no correlation between genetic diversity and viral load. Finally, we calculated the nonsynonymous and synonymous substitution rates and showed that the ratio of nonsynonymous to synonymous substitution compared to the value of one may be a poor indicator of natural selection.

MODELTEST: testing the model of DNA substitution

SUMMARY

The program MODELTEST uses log likelihood scores to establish the model of DNA evolution that best fits the data.

AVAILABILITY

The MODELTEST package, including the source code and some documentation is available at http://bioag.byu. edu/zoology/crandall_lab/modeltest.html.

Cocaine abuse and HIV-1 infection: epidemiology and neuropathogenesis

The epidemiology of cocaine abuse and potential relationships of cocaine withdrawal to human immunodeficiency virus type 1 (HIV-1)-associated dementia (HAD) are discussed. Neuroendocrinological changes in HIV-1 infection of the central nervous system (CNS) are discussed with the relevant impact of cocaine abuse. HIV-1 load in the brain tissue of infected substance users is described along with possible associations with neuropathology and HAD. Finally, the molecular epidemiology and sequence heterogeneity of HIV-1 and their implications for neuropathogenesis are summarized. The complex context of addressing cocaine abuse in the setting of HIV-1 infection appears more tractable when decomposed into its components.

The molecular evolution of visual pigments of freshwater crayfishes (Decapoda: Cambaridae)

This study examines the diverse maximum wavelength absorption (lambdamax) found in crayfishes (Decapoda: Cambaridae and Parastacidae) and the associated genetic variation in their opsin locus. We measured the wavelength absorption in the photoreceptors of six species that inhabit environments of different light intensities (i.e., burrows, streams, standing waters, and subterranean waters). Our results indicate that there is relatively little variation in lambdamax (522-530 nm) among species from different genera and families. The existing variation did not correlate with the habitat differences of the crayfishes studied. We simultaneously sequenced the rhodopsin gene to identify the amino acid replacements that affect shifts in maximum wavelength absorption. We then related these to changes that correlated with shifts in lambdamax by reconstructing ancestral character states using a maximum-likelihood approach. Using amino acid sequences obtained from five species (all were 301 amino acids in length), we identified a number of candidates for producing shifts of 4 to 8 nm in lambdamax. These amino acid replacements occurred in similar regions to those involved in spectral shifts in vertebrates.

Multiple and ancient origins of the domestic dog

Mitochondrial DNA control region sequences were analyzed from 162 wolves at 27 localities worldwide and from 140 domestic dogs representing 67 breeds. Sequences from both dogs and wolves showed considerable diversity and supported the hypothesis that wolves were the ancestors of dogs. Most dog sequences belonged to a divergent monophyletic clade sharing no sequences with wolves. The sequence divergence within this clade suggested that dogs originated more than 100,000 years before the present. Associations of dog haplotypes with other wolf lineages indicated episodes of admixture between wolves and dogs. Repeated genetic exchange between dog and wolf populations may have been an important source of variation for artificial selection.

Multiple interspecies transmissions of human and simian T-cell leukemia/lymphoma virus type I sequences

Using two sets of nucleotide sequences of the human and simian T-cell leukemia/lymphoma virus type I (HTLV-I/STLV-I), one consisting of 522 bp of the env gene from 70 viral strains and the other a 140-bp segment from the pol gene of 52 viral strains, I estimated cladograms based on a statistical parsimony procedure that was developed specifically to estimate within-species gene trees. An extension of a nesting procedure is offered for sequence data that forms nested clades used in hypothesis testing. The nested clades were used to test three hypotheses relating to transmission of HTLV/STLV sequences: (1) Have cross-species transmissions occurred and, if so, how many? (2) In what direction have they occurred? (3) What are the geographic relationships of these transmission events? The analyses support a range of 11-16 cross-species transmissions throughout the history of these sequences. Additionally, outgroup weights were assigned to haplotypes using arguments from coalescence theory to infer directionality of transmission events. Conclusions on geographic origins of transmission events and particular viral strains are inconclusive due to small samples and inadequate sampling design. Finally, this approach is compared directly to results obtained from a traditional maximum parsimony approach and found to be superior at establishing relationships and identifying instances of transmission.

Intraspecific phylogenetics: support for dental transmission of human immunodeficiency virus

A new method to estimate within-species gene genealogies was used to establish linkages among individuals associated with the Florida dental human immunodeficiency virus transmission case. Phylogenetic relationships were estimated from 103 nucleotide sequences from the V3 region of the env gene representing the Florida dentist, eight of his seropositive patients, and many local controls. The cladogram estimation procedure supports linkages among individuals within the previously described dental clade, whereas local controls and other patients form independent networks or are outliers in the main network, indicating more distant evolutionary relationships. A nested statistical analysis also indicates significant cohesion of the dental clade group.

Opsin phylogeny and evolution: a model for blue shifts in wavelength regulation

The vast diversity in spectral sensitivities in the vision of many organisms is mediated mostly (although not entirely) through variation in the photosensitive visual pigments (opsins) of the eye. Specifically, shifts in absorption maxima of visual pigments are thought to be a result of interactions within the binding pocket of the opsin, between amino acid side chains and the retinal chromophore, However, it has proven difficult to identify specific amino acid residues important in determining wavelength absorption maxima, especially for some of the short wavelength (blue) opsins. In this paper, a comparative phylogenetic approach was applied to opsin protein sequence data to identify residues important in opsin wavelength regulation. In essence, this approach consisted of interpreting evolutionary history as a series of experiments in which natural selection has repeatedly favored amino acid replacements of certain residues to shift the opsin absorption spectra to either shorter or longer wavelengths. Opsin protein sequences were obtained from GenBank, aligned, and used to reconstruct a phylogenetic tree. Amino acid replacements were traced along the branches of this opsin tree, focusing only on residues likely to reside within the chromophore-binding pocket. A number of functionally convergent, nonconservative amino acid replacements in independently evolved opsins with similar shifts in spectral properties were identified. In short, reconstruction of the phylogeny of the opsin molecule allowed us to track amino acid substitutions in specific sites within the opsin and to target those particular substitutions that are repeatedly associated with marked changes in peak absorbance, shifting the spectral sensitivity of the opsin toward shorter or longer wavelengths. Based on these results, we propose a model for blue shifts of opsin absorption spectra. Amino acid replacements of four polar and charged residues near the protonated Schiff base (SBH+) end of the chromophore are proposed to result in blue shifts of the opsin absorption spectra. This model may explain some of the diversity of blue opsins apparent in both vertebrates and invertebrates.

Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction

Empirical data sets of intraspecific restriction site polymorphism in Drosophila have been gathered in order to test hypotheses derived from coalescent theory. Three main ideas are tested: (1) haplotype frequency in the sample contains information on the topological position of a given haplotype in a cladogram, (2) the frequency of a haplotype is related to the number of mutational connections to other haplotypes in the cladogram and (3) geographic location can be used to infer topological positioning of haplotypes in a cladogram. These relationships can then be used to better estimate intraspecific phylogenies in two ways: (1) rooting the phylogeny and (2) resolving ambiguities in a cladogram. This information will allow one to reduce the number of alternative phylogenies and incorporate the uncertainties involved in reconstructing intraspecific phylogenies into subsequent analyses that depend heavily on the topology of the tree.

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation

We previously developed a cladistic approach to identify subsets of haplotypes defined by restriction endonuclease mapping or DNA sequencing that are associated with significant phenotypic deviations. Our approach was limited to segments of DNA in which little recombination occurs. In such cases, a cladogram can be constructed from the restriction site or sequence data that represents the evolutionary steps that interrelate the observed haplotypes. The cladogram is used to define a nested statistical design to identify mutational steps associated with significant phenotypic deviations. The central assumption behind this strategy is that any undetected mutation causing a phenotypic effect is embedded within the same evolutionary history that is represented by the cladogram. The power of this approach depends upon the confidence one has in the particular cladogram used to draw inferences. In this paper, we present a strategy for estimating the set of cladograms that are consistent with a particular sample of either restriction site or nucleotide sequence data and that includes the possibility of recombination. We first evaluate the limits of parsimony in constructing cladograms. Once these limits have been determined, we construct the set of parsimonious and nonparsimonious cladograms that is consistent with these limits. Our estimation procedure also identifies haplotypes that are candidates for being products of recombination. If recombination is extensive, our algorithm subdivides the DNA region into two or more subsections, each having little or no internal recombination. We apply this estimation procedure to three data sets to illustrate varying degrees of cladogram ambiguity and recombination.