Evolution and historical biogeography of the Eastern Atlantic skates

In my thesis, I tested the hypothesis that the diversification of the Eastern Atlantic skate faunas arose through vicariance rather than dispersal, using combined approach of molecular phylogeny reconstruction and zoogeography (namely historical biogeography). This analyses have been carried out independently on four Rajidae genera belonging to two different tribes: Rajini (Raja and Dipturus) and Amblyrajini (Rajella and Leucoraja). These taxa were selected because they displayed high species diversity and richness of endemic species in the Eastern Atlantic and Mediterranean. The verification of this hypothesis was carried out by reconstructing the best phylogenetic relationships among four genera and 26 species (including several endemism) based on mtDNA and nuDNA gene variation and several statistical approaches. Divergence times of taxa have been estimated based on molecular clock and fossil calibration to explain evolutionary patterns in the context of geological framework. Main issues are (i) the evidence that Eastern Atlantic skate evolution and displacement of species diversity occurred from pulsed geographical speciation (i.e. repeated series of parallel and independent speciation events) started in the Late Eocene-Early Miocene and they have occurred prevalently during Miocene; (ii) such relatively ancient origin of diversification has been allowed the sympatric displacement and evolution of several congeneric taxa likely because they have accumulated huge differences in the genomic and physiological/behavioural phenotypic traits; (iii) recently diverged sister species and taxa showed allopatric or parapatric evolution by the presence of oceanographic or hydrogeographical barriers which likely prevent large mixing between parapatric sister species.

Molekularbiologische Charakterisierung der beiden Untereinheiten des Mega-Hämocyanins der Schnecke Melanoides tuberculata

Bei dem 2010 von unserer Arbeitsgruppe entdeckten Mega-Hämocyanin handelt es sich um einen stark abgewandelten Typ des respiratorischen Proteins Hämocyanin, bestehend aus zwei flankierenden regulären Dekameren und einem zentralen Mega-Dekamer. Diese sind aus zwei immunologisch verschiedenen Untereinheiten mit ~400 bzw. ~550 kDa aufgebaut, die in unserer Arbeitsgruppe bereits proteinbiochemisch charakterisiert wurden. Im Zuge dieser Untersuchungen konnte zudem eine 3D-Rekonstruktion des Oligomers (13,5 MDa) mit einer Auflösung von 13Å erstellt werden. Das Ziel der vorliegenden Arbeit war die Aufklärung der Primärstruktur beider Polypeptide bei der Schnecke Melanoides tuberculata (MtH). Es gelang, die cDNAs der beiden Untereinheiten vollständig zu sequenzieren. Die zu typischen Dekameren assemblierende MtH400-Untereinheit umfasst 3445 Aminosäuren und besitzt eine theoretische Molekularmasse von 390 kDa. Nach dem Signalpeptid von 23 Aminosäuren Länge folgen die für Gastropoden-Hämocyanine typischen funktionellen Einheiten FU-a bis FU-h. Insgesamt verfügt die MtH400-Untereinheit über sechs potentielle N-Glykosylierungsstellen. Die MtH550-Untereinheit, welche mit 10 Kopien das Mega-Dekamer bildet, umfasst 4999 Aminosäuren und besitzt eine theoretische Molekularmasse von 567 kDa. Damit handelt es sich bei dieser Untereinheit um die zweitgrößte jemals bei einem Protein detektierte Polypeptidkette. Die MtH550-Untereinheit besteht aus einem Signalpeptid von 20 Aminosäuren Länge und den typischen Wand-FUs (FU-a bis FU-f). Daran anschließend folgen sechs weitere Varianten der FU-f (FU-f1 bis FU-f6). Die MtH550-Untereinheit verfügt über insgesamt zwölf potentielle N-Glykosylierungsstellen. Anhand der ermittelten Primärstrukturdaten wird klar, dass der auffällig vergrößerte Kragenbereich des Mega-Dekamers aus je 10 Kopien der FU-f1 bis FU-f6 besteht. Die ermittelten Sequenzdaten der beiden MtH-Untereinheiten weisen im Vergleich zu anderen Hämocyanin Sequenzen einige sehr charakteristische Indels sowie unübliche N-Glykosylierungsstellen auf. Es war zudem möglich, anhand einer molekularen Uhr den Entstehungszeitpunkt des Mega-Hämocyanins zu datieren (145 ± 35 MYA). Sowohl die Topologie als auch die berechneten Trennungszeitpunkte des an allen Verzweigungen gut unterstützten Stammbaums stimmen mit den bisher publizierten und auf Hämocyanindaten basierenden molekularen Uhren überein.

Evolution von Miniaturisierung in arktisch-alpinen Lebensräumen in Petasites Mill., Endocellion Turcz. ex Herder, Homogyne Cass. und Tussilago L. (Asteraceae) sowie Soldanella L. (Primulaceae)

Hintergrund: Miniaturisierung ist ein häufig beobachtetes Phänomen bei Pflanzen in arktisch-alpinen Lebensräumen und wird als Anpassung an niedrige Jahresmitteltemperaturen und eine kurze Vegetationsperiode interpretiert. Ziele: In der vorliegenden Arbeit wird im Petasites-Clade (Petasites Mill., Endocellion Turcz. ex Herder, Homogyne Cass., Tussilago L.; Asteraceae) und in Soldanella (Primulaceae) die Evolution der Miniaturisierung arktisch-alpiner Arten untersucht. Zudem wird innerhalb von Homogyne untersucht, ob unterschiedliche edaphische Präferenz von H. alpina (variabel) und H. discolor (kalkliebend) genetisch fixiert ist. rnMethoden: Molekulare Phylogenien des Petasites-Clades und von Soldanella wurden mit nukleären und plastidären Markern erstellt, und mit den in den Alpen vorkommenden Soldanella-Arten wurde zudem eine Fingerprint-Studie (AFLPs) gemacht. Zur Datierung der Diversifizierungsereignisse im Petasites-Clade diente eine molekulare Uhr, und die Evolution von Miniaturisierung wurde rekonstruiert. Mit H. alpina und H. discolor wurde ein vergleichendes Kulturexperiment durchgeführt.rnErgebnisse: Miniaturisierung entstand mehrere Male unabhängig voneinander in den arktisch-alpinen Vertretern des Petasites-Clade, aber nicht alle arktisch-alpinen Arten sind klein. Das Alter der arktisch-alpinen Arten deutet darauf hin, dass diese Taxa ihren Ursprung in der arkto-tertiären Flora haben. In Soldanella sind reduzierte Blütenmorphologie sowie Kleinwüchsigkeit der beiden alpinen Arten zweimal parallel entstanden. Homogyne alpina und H. discolor zeigen keine edaphischen Unterschiede hinsichtlich des Keimverhaltens, aber in Kultur zeigt sich, dass die Präferenz von H. discolor für Kalk wahrscheinlich genetisch fixiert ist.rnSchlussfolgerungen: Miniaturisierung von Pflanzen in größerer Höhe und höherer geographischer Breite kann in der Regel beobachtet werden. Allerdings kann die Evolution arktisch-alpiner Arten auch durch Faktoren wie Nährstoffverfügbarkeit, Konkurrenz und Störung beeinflusst werden, die dem Effekt der Temperatur entgegenwirken, so dass nicht alle Pflanzen in arktisch-alpinen Habitaten klein sind. Blütenmorphologische Reduktion in Soldanella kann als Anpassung an einen höheren Grad an Selbstbestäubung interpretiert werden, um eine geringere Bestäuberaktivität im alpinen Lebensraum zu kompensieren.

Rhythms of anorexigenic neuropeptides in senegalese sole: Effects of light, feeding regimes and developmental stage.

Increasing knowledge on the endocrine mechanisms that regulate feeding and growth in cultured fish can contribute to make improvement in fish holding conditions and feeding strategies, supporting the development of new techniques that could ameliorate feeding, food conversion efficiency and growth in aquaculture practice. The main objective of this study was to investigate how daily mRNA expression of three specific anorexigenic hormones, i.e. the corticotropin-releasing hormone (CRH) and the paralogues α- and β- proopiomelanocortin (POMC), is modulated by different photoperiods, light spectra and feeding regimes, in both adult and larvae of Solea senegalensis. In addition, as Senegalese sole exhibits a shift from diurnal to nocturnal in locomotor activity and feeding habits during metamorphic process, we tried to elucidate if this shift is accompanied by relevant daily variations in the expression of these anorexigenic hormones before, during and after the completion of metamorphosis. In order to reach this main objective, three main experiments were developed. In a first experiment, adults were reared under LD (12 h light: 12h dark) cycle and fed at mid-light (ML), mid-dark (MD) and at random (RND). In a second experiment, adult specimens were reared in constant darkness (DD) and fed at subjective mid-light (sML) or at RND. Larvae of Senegalese sole were reared under LD cycle with white, blue or red light for 40 days. Our results show an independence of crh mRNA expression from the feeding time and suggest an endogenous control of crh expression in sole. Both pomc paralogues showed significant daily rhythms under LD conditions. The rhythms were maintained or were even more robust under DD conditions for pomc_a, but were completely abolished for pomc_b. Our results indicate an endogenous control of pomc_a expression by the molecular clock in telencephalon and diencephalon, but not in the pituitary gland. Our findings confirm for the first time the significant influence that ambient lighting has on larval growth and development in Senegalese sole, revealing an important effect of light spectra upon functional elements of this species. Our results also emphasize the importance of maintaining cycling light-dark conditions of the adequate wavelengths in aquaculture practices during early development of sole.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Macroevolutionary patterns in the diversification of parrots: effects of climate change, geological events and key innovations

Aim Parrots are thought to have originated on Gondwana during the Cretaceous. The initial split within crown group parrots separated the New Zealand taxa from the remaining extant species and was considered to coincide with the separation of New Zealand from Gondwana 82-85 Ma, assuming that the diversification of parrots was mainly shaped by vicariance. However, the distribution patterns of several extant parrot groups cannot be explained without invoking transoceanic dispersal, challenging this assumption. Here, we present a temporal and spatial framework for the diversification of parrots using external avian fossils as calibration points in order to evaluate the relative importance of the influences of past climate change, plate tectonics and ecological opportunity. Location Australasian, African, Indo-Malayan and Neotropical regions. Methods Phylogenetic relationships were investigated using partial sequences of the nuclear genes c-mos, RAG-1 and Zenk of 75 parrot and 21 other avian taxa. Divergence dates and confidence intervals were estimated using a Bayesian relaxed molecular clock approach. Biogeographic patterns were evaluated taking temporal connectivity between areas into account. We tested whether diversification remained constant over time and if some parrot groups were more species-rich than expected given their age. Results Crown group diversification of parrots started only about 58 Ma, in the Palaeogene, significantly later than previously thought. The Australasian lories and possibly also the Neotropical Arini were found to be unexpectedly species-rich. Diversification rates probably increased around the Eocene/Oligocene boundary and in the middle Miocene, during two periods of major global climatic aberrations characterized by global cooling. Main conclusions The diversification of parrots was shaped by climatic and geological events as well as by key innovations. Initial vicariance events caused by continental break-up were followed by transoceanic dispersal and local radiations. Habitat shifts caused by climate change and mountain orogenesis may have acted as a catalyst to the diversification by providing new ecological opportunities and challenges as well as by causing isolation as a result of habitat fragmentation. The lories constitute the only highly nectarivorous parrot clade, and their diet shift, associated with morphological innovation, may have acted as an evolutionary key innovation, allowing them to explore underutilized niches and promoting their diversification.

Age of cichlids: new dates for ancient lake fish radiations

Timing divergence events allow us to infer the conditions under which biodiversity has evolved and gain important insights into the mechanisms driving evolution. Cichlid fishes are a model system for studying speciation and adaptive radiation, yet, we have lacked reliable timescales for their evolution. Phylogenetic reconstructions are consistent with cichlid origins prior to Gondwanan landmass fragmentation 121-165 MYA, considerably earlier than the first known fossil cichlids (Eocene). We examined the timing of cichlid evolution using a relaxed molecular clock calibrated with geological estimates for the ages of 1) Gondwanan fragmentation and 2) cichlid fossils. Timescales of cichlid evolution derived from fossil-dated phylogenies of other bony fishes most closely matched those suggested by Gondwanan breakup calibrations, suggesting the Eocene origins and marine dispersal implied by the cichlid fossil record may be due to its incompleteness. Using Gondwanan calibrations, we found accumulation of genetic diversity within the radiating lineages of the African Lakes Malawi, Victoria and Barombi Mbo, and Palaeolake Makgadikgadi began around or after the time of lake basin formation. These calibrations also suggest Lake Tanganyika was colonized independently by the major radiating cichlid tribes that then began to accumulate genetic diversity thereafter. These results contrast with the widely accepted theory that diversification into major lineages took place within the Tanganyika basin. Together, this evidence suggests that ancient lake habitats have played a key role in generating and maintaining diversity within radiating lineages and also that lakes may have captured preexisting cichlid diversity from multiple sources from which adaptive radiations have evolved.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last similar to 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events. (C) 2009 Elsevier B.V. All rights reserved.

The Accumulation of Reproductive Incompatibilities in African Cichlid Fish

The rate at which different components of reproductive isolation accumulate with divergence time between species has only been studied in a limited, but growing, number of species. We measured premating isolation and hybrid inviability at four different ontogenetic stages from zygotes to adults in interspecific hybrids of 26 pairs of African cichlid species, spanning the entire East African haplochromine radiation. We then used multiple relaxed molecular clock calibrations to translate genetic distances into absolute ages to compare evolutionary rates of different components of reproductive isolation. We find that premating isolation accumulates fast initially but then changes little with increasing genetic distance between species. In contrast, postmating isolation between closely related species is negligible but then accumulates rapidly, resulting in complete hybrid inviability after 4.4/8.5/18.4 million years (my). Thus, the rate at which complete intrinsic incompatibilities arise in this system is orders of magnitude lower than rates of speciation within individual lake radiations. Together these results suggest divergent ecological adaptations may prevent populations from interbreeding and help maintain cichlid species diversity, which may be vulnerable to environmental degradation. By quantifying the capacity to produce viable hybrids between allopatric, distantly related lineages our results also provide an upper divergence time limit for the "hybrid swarm origin" model of adaptive radiation.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.