Jon P. Anderson

Likelihood-Based Tests of Topologies in Phylogenetics

Likelihood-based statistical tests of competing evolutionary hypotheses (tree topologies) have been available for approximately a decade. By far the most commonly used is the Kishino-Hasegawa test. However, the assumptions that have to be made to ensure the validity of the Kishino-Hasegawa test place important restrictions on its applicability. In particular, it is only valid when the topologies being compared are specified a priori. Unfortunately, this means that the Kishino-Hasegawa test may be severely biased in many cases in which it is now commonly used: for example, in any case in which one of the competing topologies has been selected for testing because it is the maximum likelihood topology for the data set at hand. We review the theory of the Kishino-Hasegawa test and contend that for the majority of popular applications this test should not be used. Previously published results from invalid applications of the Kishino-Hasegawa test should be treated extremely cautiously, and future applications should use appropriate alternative tests instead. We review such alternative tests, both nonparametric and parametric, and give two examples which illustrate the importance of our contentions.

Multiple mutations and cancer

Most human tumors are highly heterogenous. We have hypothesized that this heterogeneity results from a mutator phenotype. Our premise is that normal mutation rates are insufficient to account for the multiple mutations found in human cancers, and, instead, that cancers must exhibit a mutator phenotype early during their evolution. Here, we examine the current status and implications of the mutator phenotype hypothesis for the prognosis, treatment, and prevention of human cancers.

Testing the Hypothesis of a Recombinant Origin of Human Immunodeficiency Virus Type 1 Subtype E

ABSTRACT The human immunodeficiency virus type 1 (HIV-1) epidemic in Southeast Asia has been largely due to the emergence of clade E (HIV-1E). It has been suggested that HIV-1E is derived from a recombinant lineage of subtype A (HIV-1A) and subtype E, with multiple breakpoints along the E genome. We obtained complete genome sequences of clade E viruses from Thailand (93TH057 and 93TH065) and from the Central African Republic (90CF11697 and 90CF4071), increasing the total number of HIV-1E complete genome sequences available to seven. Phylogenetic analysis of complete genomes showed that subtypes A and E are themselves monophyletic, although together they also form a larger monophyletic group. The apparent phylogenetic incongruence at different regions of the genome that was previously taken as evidence of recombination is shown to be not statistically significant. Furthermore, simulations indicate that bootscanning and pairwise distance results, previously used as evidence for recombination, can be misleading, particularly when there are differences in substitution or evolutionary rates across the genomes of different subtypes. Taken jointly, our analyses suggest that there is inadequate support for the hypothesis that subtype E variants are derived from a recombinant lineage. In contrast, many other HIV strains claimed to have a recombinant origin, including viruses for which only a single parental strain was employed for analysis, do indeed satisfy the statistical criteria we propose. Thus, while intersubtype recombinant HIV strains are indeed circulating, the criteria for assigning a recombinant origin to viral structures should include statistical testing of alternative hypotheses to avoid inappropriate assignments that would obscure the true evolutionary properties of these viruses.

Incorporation of reporter-labeled nucleotides by DNA polymerases

The incorporation of fluorescently labeled nucleotides into DNA by DNA polymerases has been used extensively for tagging genes and for labeling DNA. However, we lack studies comparing polymerase efficiencies for incorporating different fluorescently labeled nucleotides. We analyzed the incorporation of fluorescent deoxynucleoside triphosphates by 10 different DNA polymerases, representing a cross-section of DNA polymerases from families A, B, and reverse transcriptase. The substitution of one or more different reporter-labeled nucleotides for the cognate nucleotides was initially investigated by using an in vitro polymerase extension filter-binding assay with natural DNA as a template. Further analysis on longer DNA fragments containing one or more nucleotide analogs was performed using a newly developed extension cut assay. The results indicate that incorporation of fluorescent nucleotides is dependent on the DNA polymerase, fluorophore, linker between the nucleotide and the fluorophore, and position for attachment of the linker and the cognate nucleotide. Of the polymerases tested, Taq and Vent exo DNA polymerases were most efficient at incorporating a variety of fluorescently labeled nucleotides. This study suggests that it should be feasible to copy DNA with reactions mixtures that contain all four fluorescently labeled nucleotides allowing for high-density labeling of DNA.

FVB/N (H2(q)) mouse is resistant to arthritis induction and exhibits a genomic deletion of T-cell receptor V beta gene segments.

Abstract Animal models of autoimmune diseases have been instrumental in advancing our understanding of autoimmunity in humans. Collagen-induced arthritis (CIA) in mice is an autoimmune disease model of rheumatoid arthritis. Susceptibility to CIA in mice is linked to genes of the major histocompatibility complex (MHC). CD4+ T cells that express the T-cell receptor (TCR) Tcra-V11.1 and/or Tcrb-V8.2 play a key role in the pathogenesis of arthritis in the DBA/1 mouse (H2q). We identified an inbred mouse strain, FVB/NJ (H2q), that is resistant to arthritis induction and exhibits a genomic deletion of certain Tcrb-V gene segments. We report a novel polymerase chain reaction-based method for the rapid identification of new mouse strains that exhibit germline Tcrb-V gene deletions. We mapped for the first time both the 5′ and 3′ breakpoints of the Tcrb-V deletion in the FVB/NJ, SWR, SJL, C57L, and C57BR strains to within 1.1 kilobases. Since there is an association between a particular Tcra-V allele (Tcra-V11.1d) and arthritis susceptibility in H2q mouse strains, we examined the allelic polymorphisms of the Tcra-V11 gene subfamily members between the arthritis-susceptible DBA/1 mouse and the arthritis-resistant FVB/NJ mouse strain. The amino acid sequences of the Tcra-V11.1 alleles differ at two positions (codons 18 and 68). Therefore, the resistance of FVB/NJ mouse to arthritis induction may be due in part to Tcra-V11.1 coding sequence polymorphism and Tcrb-V8.2 gene segment deletion, as we have recently demonstrated in the case of SWR mouse strain.

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Abstract. Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences.

T-cell receptor vbeta deletion and valpha polymorphism are responsible for the resistance of SWR mouse to arthritis induction.

Abstract Collagen type II-induced arthritis (CIA) develops in susceptible mouse strains after intradermal injections of type II collagen (CII) in complete Freunds adjuvant (CFA). Susceptibility to CIA in mice is linked to genes of the major histocompatibility complex (MHC). Although the SWR mouse has a susceptible MHC haplotype (H2q), it is resistant to CIA. SWR exhibits at least two known immunological defects: (1) it contains a germline deletion of about 50% of T-cell receptor (TCR) Vβ-chain gene segments, and (2) SWR is deficient in complement component C5. It has been shown that T cells that express TCRVα11.1 and TCRVβ8.2 play a substantial role in the pathogenesis of arthritis in the DBA/1 mouse (H2q). We generated SWR transgenic (tg) mice to determine whether the expression of pathogenic Vα11.1 and/or Vβ8.2 transgenes would confer arthritis susceptibility. Arthritis was induced in the SWR TCRαβ tg mice, but not in SWR TCRβ tg mice. To address the role of Vα11.1 in arthritis susceptibility, we examined the allelic polymorphisms of the Tcra-V11-gene subfamily members between the arthritis susceptible DBA/1 mouse and the arthritis-resistant SWR mouse strain. The amino acid sequences of the Vα11.1 alleles differ at two positions (codons 18 and 68). Accordingly, these two amino acid changes are sufficient to allow the production of pathogenic T cells in SWR mice. This is the first demonstration of the association of a particular Tcra-V allele and arthritis susceptibility in mice.