Substantial Sequence Divergence Research Articles

Beyond mtDNA: nuclear gene flow suggests taxonomic oversplitting in the little brown bat (Myotis lucifugus)

In southern Alberta and north-central Montana, there is substantial mtDNA sequence divergence between two groups of the little brown bat, Myotis lucifugus (LeConte, 1831), previously thought to be subspecies ( Myotis lucifugus lucifugus and Myotis lucifugus carissima ) but recently hypothesized to be species. We tested this hypothesis using population genetic techniques. Using nuclear microsatellite markers (10 loci), we found a lack of differentiation between these two groups of bats (ascribed based on mitochondrial hypervariable region II sequence), suggesting interbreeding was sympatric. Our findings add to the recent discovery that M. l. lucifugus haplotypes are found throughout the range previously thought to be only M. l. carissima, suggesting widespread sympatry and extensively mixed gene pools, thus refuting the cryptic species hypothesis. Clinal morphology and individual variation demonstrated the impossibility to differentiate groups based on original subspecies definitions. The lack of geographic and morphological boundaries, in addition to the likelihood that the interbreeding observed in this study is occurring across western North America, suggests that no line can be drawn between these two groups. We thus suggest that the carissima subspecies designation be dropped. This study highlights the importance of investigating nuclear gene flow in widely sympatric animals suspected of being cryptic genetic species, and has important implications for applications of the DNA Barcoding Project.

Modular Organization and Combinatorial Energetics of Proline–Tyrosine Nuclear Localization Signals

Proline–tyrosine nuclear localization signals (PY-NLSs) are recognized and transported into the nucleus by human Karyopherin (Kap) β2/Transportin and yeast Kap104p. Multipartite PY-NLSs are highly diverse in sequence and structure, share a common C-terminal R/H/KX2–5PY motif, and can be subdivided into hydrophobic and basic subclasses based on loose N-terminal sequence motifs. PY-NLS variability is consistent with weak consensus motifs, but such diversity potentially renders comprehensive genome-scale searches intractable. Here, we use yeast Kap104p as a model system to understand the energetic organization of this NLS. First, we show that Kap104p substrates contain PY-NLSs, demonstrating their generality across eukaryotes. Previously reported Kapβ2–NLS structures explain Kap104p specificity for the basic PY-NLS. More importantly, thermodynamic analyses revealed physical properties that govern PY-NLS binding affinity: (1) PY-NLSs contain three energetically significant linear epitopes, (2) each epitope accommodates substantial sequence diversity, within defined limits, (3) the epitopes are energetically quasi-independent, and (4) a given linear epitope can contribute differently to total binding energy in different PY-NLSs, amplifying signal diversity through combinatorial mixing of energetically weak and strong motifs. The modular organization of the PY-NLS coupled with its combinatorial energetics lays a path to decode this diverse and evolvable signal for future comprehensive genome-scale identification of nuclear import substrates.

A Trend of Central Versus Peripheral Structuring in Mitochondrial and Nuclear Gene Sequences of the Japanese Wood Mouse, Apodemus speciosus

A phylogeographic analysis was performed on Japanese endemic wood mice (Apodemus speciosus) using nuclear interphotoreceptor retinol binding protein (IRBP) gene sequences (1,152 bp), together with previously published mitochondrial cytochrome b (cyt b) data. In the IRBP analysis, 40 haplotypes were recovered from 84 individuals by statistical and subcloning methods. Substantial sequence variation was determined from the IRBP data (pi=0.0047), and no significant evidence of recombination was detected. From the phylogenetic analysis, the 40 haplotypes fell into two major groups with geographic associations, irrespective of the karyotype groups (2n=46 and 2n=48), yielding a trend of central (Hokkaido, Honshu, Shikoku, Kyushu, and Sado) and peripheral (Izu, Oki, Tsushima, and Satsunan Is.) groupings. This geographic pattern is similar to that observed in the cyt b data, with a different insular grouping of Sado, Hokkaido, Izu, and Satsunan islands, and also to that of morphological features. In both gene data sets, nested clade analyses revealed allopatric fragmentation in the "peripheral island clades" and range expansion in the "central island clades." A mismatch analysis using cyt b data also suggested expansion of the central islands clade. Thus, the trend of central vs. peripheral structuring may be attributable to past demographic dynamics in the two distinct haplotype clades, such as range expansion of one clade in the central area of the Japanese Islands, leaving the other clade in the periphery.

Molecular genetic variation within and among isolates of QPX (Thraustochytridae), a parasite of the hard clam Mercenaria mercenaria

The thraustochytrid known as QPX (Quahog Parasite Unknown) has sporadically caused disease in the hard clam Mercenaria mercenaria along the east coast of North America since the 1960s. We hypothesized that genetically distinct QPX strains might be responsible for outbreaks of QPX disease in different areas and tested this hypothesis by comparing several QPX isolates recovered from the recent outbreak in Raritan Bay, New York with QPX strains isolated from 2 outbreaks in Massachusetts, USA. There was no variation in small subunit rDNA (SSU rDNA), 5.8S rDNA, or 4 mitochondrial gene sequences. In contrast, both of the ribosomal ribonucleic acid (rRNA) operon intergenic spacers, internal transcribed spacers 1 and 2 (ITS1 and ITS2), revealed substantial sequence variation. However, strain-specific sequences were not detected because the ITS sequence variation within QPX isolates was comparable to the variation between isolates. ITS1 sequences recovered from an infected clam by amplification with a QPX ITS2-specific primer were identical to those recovered from the QPX isolates.

Helicobacter pylori genome variability in a framework of familial transmission

BackgroundHelicobacter pylori infection is exceptionally prevalent and is considered to be acquired primarily early in life through person-to-person transmission within the family. H. pylori is a genetically diverse bacterial species, which may facilitate adaptation to new hosts and persistence for decades. The present study aimed to explore the genetic diversity of clonal isolates from a mother and her three children in order to shed light on H. pylori transmission and host adaptation.ResultsTwo different H. pylori strains and strain variants were identified in the family members by PCR-based molecular typing and sequencing of five loci. Genome diversity was further assessed for 15 isolates by comparative microarray hybridizations. The microarray consisted of 1,745 oligonucleotides representing the genes of two previously sequenced H. pylori strains. The microarray analysis detected a limited mean number (± standard error) of divergent genes between clonal isolates from the same and different individuals (1 ± 0.4, 0.1%, and 3 ± 0.3, 0.2%, respectively). There was considerable variability between the two different strains in the family members (147 ± 4, 8%) and for all isolates relative to the two sequenced reference strains (314 ± 16, 18%). The diversity between different strains was associated with gene functional classes related to DNA metabolism and the cell envelope.ConclusionThe present data from clonal H. pylori isolates of family members do not support that transmission and host adaptation are associated with substantial sequence diversity in the bacterial genome. However, important phenotypic modifications may be determined by additional genetic mechanisms, such as phase-variation. Our findings can aid further exploration of H. pylori genetic diversity and adaptation.

Genotyping of Hepatitis C virus (HCV) in infected patients from South India

Hepatitis C virus shows substantial nucleotide sequence diversity distributed throughout the viral genome. In the present study genotyping for Hepatitis C virus (HCV) infected patients was based on RFLP analysis of 5′ UTR and using type specific primers of NS5B regions. It was observed that 60% of the patients (30 patients with chronic hepatitis) were infected with variants of genotype 1 and 40% of the patients (4 chronic hepatitis patients, 12 patients with chronic renal failure and 4 cirrhosis) were infected with variants of type 3 of HCV. None of the cirrhotic patients and patients with chronic renal failure, in the present study, were infected with type1 of HCV. While PCR-RFLP, typing was rapid in conjunction with the primers used for RT-PCR, NS5 typing was helpful in determining the subtype. There was good correlation between the two typing methods and this method can be used as a cost-effective method for studying large number of samples. The study shows that predominant genotypes of HCV in South India include type 1 and 3. Type 3 seems to be transmitted nosocomially as suggested by the results in patients with chronic renal failure, as these patients are exposed to multiple medical interventions.

Tandem gene arrays in Trypanosoma brucei: Comparative phylogenomic analysis of duplicate sequence variation

BackgroundThe genome sequence of the protistan parasite Trypanosoma brucei contains many tandem gene arrays. Gene duplicates are created through tandem duplication and are expressed through polycistronic transcription, suggesting that the primary purpose of long, tandem arrays is to increase gene dosage in an environment where individual gene promoters are absent. This report presents the first account of the tandem gene arrays in the T. brucei genome, employing several related genome sequences to establish how variation is created and removed.ResultsA systematic survey of tandem gene arrays showed that substantial sequence variation existed across the genome; variation from different regions of an array often produced inconsistent phylogenetic affinities. Phylogenetic relationships of gene duplicates were consistent with concerted evolution being a widespread homogenising force. However, tandem duplicates were not usually identical; therefore, any homogenising effect was coincident with divergence among duplicates. Allelic gene conversion was detected using various criteria and was apparently able to both remove and introduce sequence variation. Tandem arrays containing structural heterogeneity demonstrated how sequence homogenisation and differentiation can occur within a single locus.ConclusionThe use of multiple genome sequences in a comparative analysis of tandem gene arrays identified substantial sequence variation among gene duplicates. The distribution of sequence variation is determined by a dynamic balance of conservative and innovative evolutionary forces. Gene trees from various species showed that intraspecific duplicates evolve in concert, perhaps through frequent gene conversion, although this does not prevent sequence divergence, especially where structural heterogeneity physically separates a duplicate from its neighbours. In describing dynamics of sequence variation that have consequences beyond gene dosage, this survey provides a basis for uncovering the hidden functionality within tandem gene arrays in trypanosomatids.

Molecular species identification boosts bat diversity

The lack of obvious morphological differences between species impedes the identification of species in many groups of organisms. Meanwhile, DNA-based approaches are increasingly used to survey biological diversity. In this study we show that sequencing the mitochondrial protein-coding gene NADH dehydrogenase, subunit 1 (nd1) from 534 bats of the Western Palaearctic region corroborates the promise of DNA barcodes in two major respects. First, species described with classical taxonomic tools can be genetically identified with only a few exceptions. Second, substantial sequence divergence suggests an unexpected high number of undiscovered species.

Characterization of a Highly Conserved Island in the Otherwise Divergent Bordetella holmesii and Bordetella pertussis Genomes

The recently discovered pathogen Bordetella holmesii has been isolated from the airways and blood of diseased humans. Genetic events contributing to the emergence of B. holmesii are not understood, and its phylogenetic position among the bordetellae remains unclear. To address these questions, B. holmesii strains were analyzed by comparative genomic hybridization (CGH) to a Bordetella pertussis microarray and by multilocus sequence typing. Both methods indicated substantial sequence divergence between B. pertussis and B. holmesii. However, CGH identified a putative pathogenicity island of 66 kb that is highly conserved between these species and contains several IS481 elements that may have been laterally transferred from B. pertussis to B. holmesii. This island contains, among other genes, a functional, iron-regulated locus encoding the biosynthesis, export, and uptake of the siderophore alcaligin. The acquisition of this genomic island by B. holmesii may have significantly contributed to its emergence as a human pathogen. Horizontal gene transfer between B. pertussis and B. holmesii may also explain the unusually high sequence identity of their 16S rRNA genes.

Spectroscopic Analysis of Benzylidene Anabaseine Complexes with Acetylcholine Binding Proteins as Models for Ligand−Nicotinic Receptor Interactions

The discovery of the acetylcholine binding proteins (AChBPs) has provided critical soluble surrogates for examining structure and ligand interactions with nicotinic receptors and related pentameric ligand-gated ion channels. The multiple marine and freshwater sources of AChBP constitute a protein family with substantial sequence divergence and selectivity in ligand recognition for analyzing structure-activity relationships. The purification of AChBP in substantial quantities in the absence of a detergent enables one to conduct spectroscopic studies of the ligand-AChBP complexes. To this end, we have examined the interaction of a congeneric series of benzylidene-ring substituted anabaseines with AChBPs from Lymnaea, Aplysia, and Bulinus species and correlated their binding energetics with spectroscopic changes associated with ligand binding. The anabaseines display agonist activity on the alpha7 nicotinic receptor, a homomeric receptor with sequences similar to those of the AChBPs. Substituted anabaseines show absorbance and fluorescence properties sensitive to the protonation state, relative permittivity (dielectric constant), and the polarizability of the surrounding solvent or the proximal residues in the binding site. Absorbance difference spectra reveal that a single protonation state of the ligand binds to AChBP and that the bound ligand experiences a solvent environment with a high degree of polarizability. Changes in the fluorescence quantum yield of the bound ligand reflect the rigidification of the ring system of the bound ligand. Hence, the spectral properties of the bound ligand allow a description of the electronic character of the bound state of the ligand within its aromatic binding pocket and provide information complementary to that of crystal structures in defining the determinants of interaction.

Occurrence, hosts, morphology, and molecular characterisation of Pasteuria bacteria parasitic in nematodes of the family Plectidae

Parasitic bacteria of the genus Pasteuria are reported for three Anaplectus and four identified and several unidentified Plectus species found in eight countries in various habitats. The pasteurias from plectids agree in essential morphological characters of sporangia and endospores as well as in developmental cycle with those of the Pasteuria species and strains described from tylenchid nematodes, but appear to be mainly distinguished from these by absence of a distinct perisporium in the spores and the endospores obviously not being cup- or saucer-shaped. The wide range of measurements and morphological peculiarities of sporangia and endospores suggest that probably several Pasteuria species have to be distinguished as parasites in Plectidae. From an infected juvenile of an unidentified plectid species the 16S rRNA gene sequence of Pasteuria sp. was obtained. Substantial sequence divergence from described Pasteuria species and its phylogenetic position on molecular trees indicate that this Pasteuria sp. could be considered as a new species. Preliminary results of the analysis of DNA phylogeny of Pasteuria spp. and their nematode hosts provide evidence for incongruence of their phylogenetic history and of host switching events during evolution of the bacterial parasites.

A Shallow BTEX and MTBE Contaminated Aquifer Supports a Diverse Microbial Community

Microbial communities in subsurface environments are poorly characterized and the impacts of anthropogenic contamination on their structure and function have not been adequately addressed. The release of contaminant(s) to a previously unexposed environment is often hypothesized to decrease the diversity of the affected community. We characterized the structure of microbial communities along a gradient of benzene, toluene, ethylbenzene, and xylene (BTEX) and methyl-tert-butyl-ether (MTBE) contamination, resulting from a petroleum spill, within a shallow sandy aquifer at Vandenberg Air Force Base (VAFB) in Lompoc, CA. Differences in microbial community composition along the contaminant plume were assessed via a combinatorial approach utilizing denaturing gradient gel electrophoresis (DGGE), cloning and sequencing, intergenic transcribed spacer analysis (ITS), and comparative phylogenetic analysis of partial 16S rDNA sequences. Substantial bacterial sequence diversity, similar levels of species richness, and similar phylo-groups (including the Cytophaga-Flavobacterium-Bacteroidetes group and numerous members of the alpha-, beta-, gamma-, delta-, and epsilon-groups of the proteobacteria) were observed in both uncontaminated and contaminated regions of the aquifer. High-resolution measures (ITS fingerprinting and phylogenetic inference) readily separated communities impacted by the original petroleum spill (in source zone) from those in other parts of the aquifer and indicated that communities exposed to MTBE only were similar to communities in uncontaminated regions. Collectively, these data suggest that petroleum contamination alters microbial community structure at the species and subspecies level. Further study is required to determine whether these changes have an impact on the functioning of this subsurface ecosystem.

Parallel adaptive radiations in arid and temperate Australia: molecular phylogeography and systematics of the Egernia whitii (Lacertilia: Scincidae) species group

It has been an enduring belief that increasing aridity combined with Pliocene-Pleistocene glacial cycles resulted in the formation of distinct arid zone and temperate zone faunas within Australia. We assembled a molecular phylogeny for the Egernia whitii species group, an endemic group of skinks that comprises representatives from arid and temperate Australia, in order to test several biogeographical hypotheses regarding the origin of the Australian arid zone fauna. Sequence data were obtained from ten of the 11 species within the species group, plus three other Egernia species and an outgroup (Eulamprus heatwolei). We targeted portions of the ND4 (696 base pairs) and 16S rRNA (500 bp) mitochondrial genes and the b-Fibrinogen 7th Intron nuclear gene (648 bp). The edited alignment comprised 1844 characters, of which 551 (30%) were variable and 382 (69%) were parsimony informative. We analysed the data using maximum likelihood and Bayesian techniques and produced a single optimal tree. Our phylogeny strongly supports two major clades within the species group, corresponding to temperate-adapted rock-dwelling species and arid-adapted obligate burrowing species. However, the phylogenetic affinities of E. pulchra were not resolved. Our topology indicates that the New South Wales population of E. margaretae is actually E. whitii and reveals that E. margaretae margaretae and E. m. personata are distinct species. There also appears to be a major phylogeographical break within E. whitii occurring in eastern Victoria. Although our data supported several previously proposed phylogenetic relationships, Shimodaira‐Hasegawa tests soundly rejected several suggested affinities between certain species. The arid zone members of the E. whitii species group had been suggested to have originated as a result of multiple periods of colonization during the Pleistocene glaciation cycles. However, our genetic data suggest a single origin (presumably from a semiarid E. multiscutata-like ancestor) for the arid zone members of the group prior to the Plio-Pleistocene, probably during the late Miocene to early Pliocene. Our topology displays substantial sequence divergence between species with short internodes and long terminal branches, indicating rapid adaptive radiations within the arid and temperate zones. The presence of temperate-adapted species within more mesic refugia of the arid zone suggests that the necessary adaptations to aridity for colonizing the dry interior of the continent have not evolved since the initial period of adaptive radiation, despite the long evolutionary history of the species group. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 83, 157‐173. ADDITIONAL KEYWORDS: 16S rRNA ‐ b-fibrinogen ‐ biogeography ‐ mtDNA ‐ ND4 ‐ nuclear gene ‐ phylogeny ‐ skink ‐ taxonomy.

Genetic variation at the mitochondrial large-subunit rRNA locus of Pneumocystis isolates from simian immunodeficiency virus-infected rhesus macaques.

The nucleotide sequences of a segment of the Pneumocystis mitochondrial large-subunit (mt LSU) rRNA gene from rhesus macaques coinfected with simian immunodeficiency virus (SIV) and Pneumocystis carinii were examined. Of 12 isolates examined, 3 were found to be identical and the others showed substantial sequence variation, with up to 13% divergence among variants. We identified two general sequence types that differed at several sites, including a conserved 26-nucleotide insertion. Four monkeys had evidence of two Pneumocystis variants present simultaneously, indicative of a mixed infection. There was a high degree of variance between the rhesus macaque-derived Pneumocystis mt LSU rRNA gene sequence and the cognate sequences in Pneumocystis organisms derived from other hosts. Analysis of the mt LSU rRNA genes of Pneumocystis organisms derived from rhesus macaques and several other mammalian hosts supports the observation that rhesus macaque-derived Pneumocystis is most closely related to human-derived PNEUMOCYSTIS: In addition, the data identify the mt LSU rRNA gene as an informative locus for transmission and epidemiological studies of the SIV-rhesus macaque model of Pneumocystis infection.

Loss of DNA recombinational repair enzymes in the initial stages of genome degeneration.

Many obligate intracellular pathogens and symbionts undergo genome degeneration during long-term association with eukaryotic hosts; however, very little is known about genome changes that occur in the initial stages of such intracellular associations. By focusing on a clade of bacteria that have recently established symbiotic associations with insect hosts, we have identified events that may contribute to the reduction and degeneration of symbiont genomes. Unlike virtually all other bacteria, the obligate symbionts of maize and rice weevils each display substantial sequence divergence between multiple copies of their rDNA genes, resulting from a reduction in the efficacy of recombinational gene conversion, coincident with the inactivation of the recombinational repair gene recF in the common ancestor of both symbionts. The maize weevil endosymbiont also lacks a functional recA, resulting in further reduction in the efficacy of gene conversion between paralogous rDNAs and in a novel IS-mediated deletion in a 23S rDNA gene. Similar events may be pervasive during the evolution of symbiosis because symbiont genomes typically lack recombinational repair genes and have reduced numbers of ribosomal operons.

Intragenomic and interspecific 5S rDNA sequence variation in five Asian pines

Patterns of intragenomic and interspecific variation of 5S rDNA in Pinus (Pinaceae) were studied by cloning and sequencing multiple 5S rDNA repeats from individual trees. Five pines, from both subgenera, Pinus and Strobus, were selected. The 5S rDNA repeat in pines has a conserved 120-base pair (bp) transcribed region and an intergenic spacer region of variable length (382-608 bp). The evolutionary rate in the spacer region is three- to sevenfold higher than in the genic region. We found substantial sequence divergence between the two subgenera. Intragenomic sequence heterogeneity was high for all species, and more than 86% of the clones within each individual were unique. The 5S gene tree revealed that different 5S repeats within individuals are polyphyletic, indicating that their ancestral divergence preceded the speciation events. The degrees of interspecific and intragenomic divergence among diploxylon pines are similar. The observed sequence patterns suggest that concerted evolution has been acting after the diversification of the two subgenera but very weak after the speciation of the four diploxylon pines. Sequence patterns in P. densata are consistent with hybrid origin. It had higher intragenomic diversity and maintained polymorphic copies of the parental types in addition to new and recombinant types unique to the hybrid.

Inhibition of hepatitis delta virus RNA editing by short inhibitory RNA-mediated knockdown of ADAR1 but not ADAR2 expression.

Hepatitis delta virus (HDV) requires host RNA editing at the viral RNA amber/W site. Of the two host genes responsible for RNA editing via deamination of adenosines in double-stranded RNAs, short inhibitory RNA-mediated knockdown of host ADAR1 expression but not that of ADAR2 led to decreased HDV amber/W editing and virus production. Despite substantial sequence and structural variation among the amber/W sites of the three HDV genotypes, ADAR1a was primarily responsible for editing all three. We conclude that ADAR1 is primarily responsible for editing HDV RNA at the amber/W site during HDV infection.

Patterns of sequence divergence in Daphniid hemoglobin genes.

In common with other multigene families, sequence diversity in the hemoglobin genes of cladoceran crustaceans has been heavily impacted by gene conversion events. Because of their structural complexity (six exons, five introns), these genes provide a good opportunity to study the influence of intron length and position on the conversion process. This study surveys the patterns of divergence in variants of one hemoglobin gene (H1) from two closely related species of Daphnia using a PCR-based approach. Although its effects were most pronounced at their 5' ends, intron and exon regions of these genes showed similar exposure to gene conversion, excepting intron 2. This intron, which was the only one with a marked length difference among variants, showed substantial sequence divergence, suggesting that gene conversion was disrupted. These results, together with those on hemoglobin gene families in other organisms, indicate that sequence tracts showing gene conversion are often distributed in a mosaic fashion. The reactivation of gene conversion downstream of a block protected from its effects suggests that there are multiple initiation points, and the distribution of conversion tracts suggests that exon/intron splice sites are important in this regard.

Pore-forming Polypeptides of the Pathogenic Protozoon Naegleria fowleri

The free-living amoeboflagellate and potential human pathogen Naegleria fowleri causes the often fatal disease primary amoebic meningoencephalitis. The molecular repertoire responsible for the cytolytic and tissue-destructive activity of this amoeboid protozoon is largely unknown. We isolated two pore-forming polypeptides from extracts of highly virulent trophozoites of N. fowleri by measuring their membrane-permeabilizing activity. N-terminal sequencing and subsequent molecular cloning yielded the complete primary structures and revealed that the two polypeptides are isoforms. Both polypeptides share similar structural properties with antimicrobial and cytolytic polypeptides of the protozoon Entamoeba histolytica (amoebapores) and of cytotoxic natural killer (NK) and T cells of human (granulysin) and pig (NK-lysin), all characterized by a structure of amphipathic alpha-helices and an invariant framework of cysteine residues involved in disulfide bonds. In contrast to the aforementioned proteins, the Naegleria polypeptides both are processed from large precursor molecules containing additional isoforms of substantial sequence divergence. Moreover, biochemical characterization of the isolated polypeptides in combination with mass determination showed that they are N-glycosylated and variably processed at the C terminus. The biological activity of the purified polypeptides of Naegleria was examined toward human cells and bacteria, and it was found that these factors, named naegleriapores, are active against both types of target cells, which is in good agreement with their proposed biological role as a broad-spectrum effector molecule.

A single ancient origin of brood parasitism in African finches: implications for host-parasite coevolution.

Robust phylogenies for brood-parasitic birds, their hosts, and nearest nesting relatives provide the framework to address historical questions about host-parasite coevolution and the origins of parasitic behavior. We tested phylogenetic hypotheses for the two genera of African brood-parasitic finches, Anomalospiza and Vidua, using mitochondrial DNA sequence data from 43 passeriform species. Our analyses strongly support a sister relationship between Vidua and Anomalospiza, leading to the conclusion that obligate brood parasitism evolved only once in African finches rather than twice, as has been the conventional view. In addition, the parasitic finches (Viduidae) are not recently derived from either weavers (Ploceidae) or grassfinches (Estrildidae), but represent a third distinct lineage. Among these three groups, the parasitic finches and estrildids, which includes the hosts of all 19 Vidua species, are sister taxa in all analyses of our full dataset. Many characters shared by Vidua and estrildids, including elaborate mouth markings in nestlings, unusual begging behavior, and immaculate white eggs, can therefore be attributed to common ancestry rather than convergent evolution. The host-specificity of mouth mimicry in Vidua species, however, is clearly the product of subsequent host-parasite coevolution. The lineage leading to Anomalospiza switched to parasitizing more distantly related Old World warblers (Sylviidae) and subsequently lost these characteristics. Substantial sequence divergence between Vidua and Anomalospiza indicates that the origin of parasitic behavior in this clade is ancient (approximately 20 million years ago), a striking contrast to the recent radiation of extant Vidua. We suggest that the parasitic finch lineage has experienced repeated cycles of host colonization, speciation, and extinction through their long history as brood parasites and that extant Vidua species represent only the latest iterations of this process. This dynamic process may account for a significantly faster rate of DNA sequence evolution in parasitic finches as compared to estrildids and other passerines. Our study reduces by one the tally of avian lineages in which obligate brood parasitism has evolved and suggests an origin of parasitism that involved relatively closely related species likely to accept and provide appropriate care to parasitic young. Given the ancient origin of parasitism in African finches, ancestral estrildids must have been parasitized well before the diversification of extant Vidua, suggesting a long history of coevolution between these lineages preceding more recent interactions between specific hosts and parasites.