Analysing recombination in nucleotide sequences

Throughout the living world, genetic recombination and nucleotide substitution are the primary processes that create the genetic variation upon which natural selection acts. Just as analyses of substitution patterns can reveal a great deal about evolution, so too can analyses of recombination. Evidence of genetic recombination within the genomes of apparently asexual species can equate with evidence of cryptic sexuality. In sexually reproducing species, nonrandom patterns of sequence exchange can provide direct evidence of population subdivisions that prevent certain individuals from mating. Although an interesting topic in its own right, an important reason for analysing recombination is to account for its potentially disruptive influences on various phylogenetic-based molecular evolution analyses. Specifically, the evolutionary histories of recombinant sequences cannot be accurately described by standard bifurcating phylogenetic trees. Taking recombination into account can therefore be pivotal to the success of selection, molecular clock and various other analyses that require adequate modelling of shared ancestry and draw increased power from accurately inferred phylogenetic trees. Here, we review various computational approaches to studying recombination and provide guidelines both on how to gain insights into this important evolutionary process and on how it can be properly accounted for during molecular evolution studies.     

Fine-Scale Phylogenetic Discordance across the House Mouse Genome

Population genetic theory predicts discordance in the true phylogeny of different genomic regions when studying recently diverged species. Despite this expectation, genome-wide discordance in young species groups has rarely been statistically quantified. The house mouse subspecies group provides a model system for examining phylogenetic discordance. House mouse subspecies are recently derived, suggesting that even if there has been a simple tree-like population history, gene trees could disagree with the population history due to incomplete lineage sorting. Subspecies of house mice also hybridize in nature, raising the possibility that recent introgression might lead to additional phylogenetic discordance. Single-locus approaches have revealed support for conflicting topologies, resulting in a subspecies tree often summarized as a polytomy. To analyze phylogenetic histories on a genomic scale, we applied a recently developed method, Bayesian concordance analysis, to dense SNP data from three closely related subspecies of house mice: Mus musculus musculus, M. m. castaneus, and M. m. domesticus. We documented substantial variation in phylogenetic history across the genome. Although each of the three possible topologies was strongly supported by a large number of loci, there was statistical evidence for a primary phylogenetic history in which M. m. musculus and M. m. castaneus are sister subspecies. These results underscore the importance of measuring phylogenetic discordance in other recently diverged groups using methods such as Bayesian concordance analysis, which are designed for this purpose.     

Phylogeny of the Avian Family Ciconiidae (Storks) Based on CytochromebSequences and DNA–DNA Hybridization Distances

This study is a phylogenetic analysis of the avian family Ciconiidae, the storks, based on two molecular data sets: 1065 base pairs of sequence from the mitochondrial cytochromebgene and a complete matrix of single-copy nuclear DNA–DNA hybridization distances. Sixteen of the nineteen stork species were included in the cytochromebdata matrix, and fifteen in the DNA–DNA hybridization matrix. Both matrices included outgroups from the families Cathartidae (New World vultures) and Threskiornithidae (ibises, spoonbills). Optimal trees based on the two data sets were congruent in those nodes with strong bootstrap support. In the best-fit tree based on DNA–DNA hybridization distances, nodes defining relationships among very recently diverged species had low bootstrap support, while nodes defining more distant relationships had strong bootstrap support. In the optimal trees based on the sequence data, nodes defining relationships among recently diverged species had strong bootstrap support, while nodes defining basal relationships in the family had weak support and were incongruent among analyses. A combinable-component consensus of the best-fit DNA–DNA hybridization tree and a consensus tree based on different analyses of the cytochromebsequences provide the best estimate of relationships among stork species based on the two data sets.     

Molecular evolution and phylogenetic utility of Wolbachia ftsZ and wsp gene sequences with special reference to the origin of male-killing

A detailed assessment of the evolution and phylogenetic utility of two genes, ftsZ and wsp, was used to investigate the origin of male-killing Wolbachia, previously isolated from the ladybird Adalia bipunctata and the butterfly Acraea encedon. The analysis included almost all available sequences of B-group Wolbachia and two outgroup taxa and showed that (1) the two gene regions differ in phylogenetic utility, (2) sequence variation is here correlated with phylogenetic information content, (3) both genes show significant rate heterogeneity between lineages, (4) increased substitution rates are associated with homoplasy in the data, (5) wsp sequences of some taxa appear to be subject to positive selection, and (6) only a limited number of clades can be inferred with confidence due to either lack of phylogenetic information or the presence of homoplasy. With respect to the evolution of male-killing, the two genes nevertheless seemed to provide unbiased information. However, they consistently produce contradictory results. Current data therefore do not permit clarification of the origin of this behavior. In addition, A. bipunctata was found to be a host to two recently diverged strains of male-killing Wolbachia that showed increased substitution rates for both genes. Moreover, the wsp gene, which codes for an outer membrane protein, was found to be subject to positive selection in these taxa. These findings were postulated to be the product of high selection pressures due to antagonistic host-symbiont interactions in this ladybird species. In conclusion, our study demonstrates that the results of a detailed phylogenetic analysis, including characterization of the limitations of such an approach, can serve as a valuable basis for an understanding of the evolution of Wolbachia bacteria. Moreover, particular features of gene evolution, such as elevated substitution rates or the presence of positive selection, may provide information about the dynamics of Wolbachia-host associations.     

Genetic delineation between and within the widespread coccolithophore morpho‐species Emiliania huxleyi and Gephyrocapsa oceanica (Haptophyta)

Emiliania huxleyi and Gephyrocapsa oceanica are abundant coccolithophore morpho‐species that play key roles in ocean carbon cycling due to their importance as both primary producers and cal‐cifiers. Global change processes such as ocean acidification impact these key calcifying species. The physiology of E. huxleyi, a developing model species, has been widely studied, but its genetic delineation from G. oceanica remains unclear due to a lack of resolution in classical genetic markers. Using nuclear (18S rDNA and 28S rDNA), mitochondrial (cox1, cox2, cox3, rpl16, and dam), and plastidial (16S rDNA, rbcL, tufA, and petA) DNA markers from 99 E. huxleyi and 44 G. oceanica strains, we conducted a multigene/multistrain survey to compare the suitability of different markers for resolving phylogenetic patterns within and between these two morpho‐species. The nuclear genes tested did not provide sufficient resolution to discriminate between the two morpho‐species that diverged only 291Kya. Typical patterns of incomplete lineage sorting were generated in phylogenetic analyses using plastidial genes. In contrast, full morpho‐species delineation was achieved with mitochondrial markers and common intra‐morpho‐species phylogenetic patterns were observed despite differing rates of DNA substitution. Mitochondrial genes are thus promising barcodes for distinguishing these coccolithophore morpho‐species, in particular in the context of environmental monitoring.     

Phylogeny of Campanuloideae (Campanulaceae) with Emphasis on the Utility of Nuclear Pentatricopeptide Repeat (PPR) Genes

Background

The Campanuloideae (Campanulaceae) are a highly diverse clade of angiosperms found mostly in the Northern Hemisphere, with the highest diversity in temperate areas of the Old World. Chloroplast markers have greatly improved our understanding of this clade but many relationships remain unclear primarily due to low levels of molecular evolution and recent and rapid divergence. Furthermore, focusing solely on maternally inherited markers such as those from the chloroplast genome may obscure processes such as hybridization. In this study we explore the phylogenetic utility of two low-copy nuclear loci from the pentatricopeptide repeat gene family (PPR). Rapidly evolving nuclear loci may provide increased phylogenetic resolution in clades containing recently diverged or closely related taxa. We present results based on both chloroplast and low-copy nuclear loci and discuss the utility of such markers to resolve evolutionary relationships and infer hybridization events within the Campanuloideae clade.

Results

The inclusion of low-copy nuclear genes into the analyses provides increased phylogenetic resolution in two species-rich clades containing recently diverged taxa. We also obtain support for the placement of two early diverging lineages (Jasione and Musschia-Gadellia clades) that have previously been unresolved. Furthermore, phylogenetic analyses of PPR loci revealed potential hybridization events for a number of taxa (e.g., Campanula pelviformis and Legousia species). These loci offer greater overall topological support than obtained with plastid DNA alone.

Conclusion

This study represents the first inclusion of low-copy nuclear genes for phylogenetic reconstruction in Campanuloideae. The two PPR loci were easy to sequence, required no cloning, and the sequence alignments were straightforward across the entire Campanuloideae clade. Although potentially complicated by incomplete lineage sorting, these markers proved useful for understanding the processes of reticulate evolution and resolving relationships at a wide range of phylogenetic levels. Our results stress the importance of including multiple, independent loci in phylogenetic analyses.     

Genetic,ecological, and behavioural divergence between two sibling snapping shrimp species (Crustacea: Decapoda: Alpheus)

Examination of genetic and ecological relationships within sibling species complexes can provide insights into species diversity and speciation processes. Alpheus angulatus and A. armillatus, two snapping shrimp species with overlapping ranges in the north-western Atlantic, are similar in morphology, exploit similar ecological niches and appear to represent recently diverged sibling species. We examined phylogenetic and ecological relationships between these two species with: (i) sequence data from two mitochondrial genes (16S rRNA and COI); (ii) data on potential differences in microhabitat distribution for A. armillatus and A. angulatus; and (iii) data from laboratory experiments on the level of reproductive isolation between the two species. DNA sequence data suggest A. armillatus and A. angulatus are sister species that diverged subsequent to the close of the Isthmus of Panama, and that haplotype diversity is lower in A. armillatus than in A. angulatus. Both species are distantly related to A. heterochaelis and A. estuariensis, two species with which A. angulatus shares some similarities in coloration. Ecological data on the distribution of A. angulatus and A. armillatus from two locations revealed differences in distribution of the two species between habitat patches, with each patch dominated by one or the other species. However, there was no apparent difference in distribution of the two species within habitat patches with respect to microhabitat location. Ecological data also revealed that heterospecific individuals often occur in close proximity (i.e. within metres or centimetres) where sympatric. Behavioural data indicated that these species are reproductively isolated, which is consistent with speciation in transient allopatry followed by post-divergence secondary contact. Our data further resolve taxonomic confusion between the sibling species, A. armillatus and A. angulatus, and suggest that sympatry in areas of range overlap and exploitation of similar ecological niches by these two recently diverged species have selected for high levels of behavioural incompatibility.     

Nested gene fusions as markers of phylogenetic branchpoints in prokaryotes

Phylogenetic trees for prokaryotic microorganisms are being assembled at a rapid pace, primarily through sequence comparisons of ribosomal RNA genes. For lineages that diverged from the ancestral stem at nearly the same time, the order of branching may be uncertain. The problem applies both to minor branches that separated very recently and to major branches that diverged long ago. Bifunctional proteins produced by gene fusion provide the clarity of a plus-or-minus character state, and analysis of the distribution of genefusion patterns can reveal the order of phylogenetic branching.     

Consequences of stop codon reassignment on protein evolution in ciliates with alternative genetic codes

Tetrahymena thermophila and Paramecium tetraurelia are ciliates that reassign TAA and TAG from stop codons to glutamine codons. Because of the lack of full genome sequences, few studies have concentrated on analyzing the effects of codon reassignment in protein evolution. We used the recently sequenced genome of these species to analyze the patterns of amino acid substitution in ciliates that reassign the code. We show that, as expected, the codon reassignment has a large impact on amino acid substitutions in closely related proteins; however, contrary to expectations, these effects also hold for very diverged proteins. Previous studies have used amino acid substitution data to calculate the minimization of the genetic code; our results show that because of the lasting influence of the code in the patterns of substitution, such studies are tautological. These different substitution patterns might affect alignment of ciliate proteins, as alignment programs use scoring matrices based on substitution patterns of organisms that use the standard code. We also show that glutamine is used more frequently in ciliates than in other species, as often as expected based on the presence of the 2 new reassigned codons, indicating that the frequencies of amino acids in proteomes is mostly determined by neutral processes based on their number of codons.     

Integrating phylogenetic and population genetic analyses of multiple loci to test species divergence hypotheses in Passerina buntings

Phylogenetic and population genetic analyses of DNA sequence data from 10 nuclear loci were used to test species divergence hypotheses within Passerina buntings, with special focus on a strongly supported, but controversial, sister relationship between Passerina amoena and P. caerulea inferred from a previous mitochondrial study. Here, a maximum-likelihood analysis of a concatenated 10-locus data set, as well as minimize-deep-coalescences and maximum-likelihood analyses of the locus-specific gene trees, recovered the traditional sister relationship between P. amoena and P. cyanea. In addition, a more recent divergence time estimate between P. amoena and P. cyanea than between P. amoena and P. caerulea provided evidence for the traditional sister relationship. These results provide a compelling example of how lineage sorting stochasticity can lead to incongruence between gene trees and species trees, and illustrate how phylogenetic and population genetic analyses can be integrated to investigate evolutionary relationships between recently diverged taxa.     

GENEALOGICAL DISCORDANCE AND PATTERNS OF INTROGRESSION AND SELECTION ACROSS A CRICKET HYBRID ZONE

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus , offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.     

Sliding MinPD: building evolutionary networks of serial samples via an automated recombination detection approach

MOTIVATION: Traditional phylogenetic methods assume tree-like evolutionary models and are likely to perform poorly when provided with sequence data from fast-evolving, recombining viruses. Furthermore, these methods assume that all the sequence data are from contemporaneous taxa, which is not valid for serially-sampled data. A more general approach is proposed here, referred to as the Sliding MinPD method, that reconstructs evolutionary networks for serially-sampled sequences in the presence of recombination. RESULTS: Sliding MinPD combines distance-based phylogenetic methods with automated recombination detection based on the best-known sliding window approaches to reconstruct serial evolutionary networks. Its performance was evaluated through comprehensive simulation studies and was also applied to a set of serially-sampled HIV sequences from a single patient. The resulting network organizations reveal unique patterns of viral evolution and may help explain the emergence of disease-associated mutants and drug-resistant strains with implications for patient prognosis and treatment strategies.     

Phylogenetic, morphological, and chemotaxonomic incongruence in the North American endemic genus Echinacea

The study of recently formed species is important because it can help us to better understand organismal divergence and the speciation process. However, these species often present difficult challenges in the field of molecular phylogenetics because the processes that drive molecular divergence can lag behind phenotypic divergence. In the current study we show that species of the recently diverged North American endemic genus of purple coneflower, Echinacea, have low levels of molecular divergence. Data from three nuclear loci and two plastid loci provide neither resolved topologies nor congruent hypotheses about species-level relationships. This lack of phylogenetic resolution is likely due to the combined effects of incomplete lineage sorting, hybridization, and backcrossing following secondary contact. The poor resolution provided by molecular markers contrasts previous studies that found well-resolved and taxonomically supported relationships from metabolic and morphological data. These results suggest that phenotypic canalization, resulting in identifiable morphological species, has occurred rapidly within Echinacea. Conversely, molecular signals have been distorted by gene flow and incomplete lineage sorting. Here we explore the impact of natural history on the genetic organization and phylogenetic relationships of Echinacea.     

Phylogenetic networks from multi-labelled trees

It is now quite well accepted that the evolutionary past of certain species is better represented by phylogenetic networks as opposed to trees. For example, polyploids are typically thought to have resulted through hybridization and duplication, processes that are probably not best represented as bifurcating speciation events. Based on the knowledge of a multi-labelled tree relating collection of polyploids, we present a canonical construction of a phylogenetic network that exhibits the tree. In addition, we prove that the resulting network is in some well-defined sense a minimal network having this property.     

Beyond barcodes: complex DNA taxonomy of a South Pacific Island radiation

DNA barcodes can provide rapid species identification and aid species inventories in taxonomically unstudied groups. However, the approach may fail in recently diverged groups with complex gene histories, such as those typically found on oceanic islands. We produced a DNA-based inventory of taxonomically little known diving beetles (genus Copelatus) in the Fiji archipelago, where they are a dominant component of the aquatic invertebrate fauna. Sampling from 25 localities on five islands and analysis of sequences from one nuclear (328bp histone 3) and three mitochondrial (492bp rrnL, 786bp cox1, 333bp cob) gene regions revealed high haplotype diversity, mainly originated since the Pleistocene, and subdivided into three major phylogenetic lineages and 22 statistical parsimony networks. A traditional taxonomic study delineated 25 morphologically defined species that were largely incongruent with the DNA-based groups. Haplotype diversity and their spatial arrangement demonstrated a continuum of relatedness in Fijian Copelatus, with evidence for introgression at various hierarchical levels. The study illustrates the difficulties for formal classification in evolutionarily complex lineages, and the potentially misleading conclusions obtained from either DNA barcodes or morphological traits alone. However, the sequence profile of Fijian Copelatus provides an evolutionary framework for the group and a DNA-based reference system for the integration of ecological and other biodiversity data, independent of the Linnaean naming system.     

Elevated substitution rates estimated from ancient DNA sequences

Ancient DNA sequences are able to offer valuable insights into molecular evolutionary processes, which are not directly accessible via modern DNA. They are particularly suitable for the estimation of substitution rates because their ages provide calibrating information in phylogenetic analyses, circumventing the difficult task of choosing independent calibration points. The substitution rates obtained from such datasets have typically been high, falling between the rates estimated from pedigrees and species phylogenies. Many of these estimates have been made using a Bayesian phylogenetic method that explicitly accommodates heterochronous data. Stimulated by recent criticism of this method, we present a comprehensive simulation study that validates its performance. For datasets of moderate size, it produces accurate estimates of rates, while appearing robust to assumptions about demographic history. We then analyse a large collection of 749 ancient and 727 modern DNA sequences from 19 species of animals, plants and bacteria. Our new estimates confirm that the substitution rates estimated from ancient DNA sequences are elevated above long-term phylogenetic levels.     

Vestige: Maximum likelihood phylogenetic footprinting

Background  

Phylogenetic footprinting is the identification of functional regions of DNA by their evolutionary conservation. This is achieved by comparing orthologous regions from multiple species and identifying the DNA regions that have diverged less than neutral DNA. Vestige is a phylogenetic footprinting package built on the PyEvolve toolkit that uses probabilistic molecular evolutionary modelling to represent aspects of sequence evolution, including the conventional divergence measure employed by other footprinting approaches. In addition to measuring the divergence, Vestige allows the expansion of the definition of a phylogenetic footprint to include variation in the distribution of any molecular evolutionary processes. This is achieved by displaying the distribution of model parameters that represent partitions of molecular evolutionary substitutions. Examination of the spatial incidence of these effects across regions of the genome can identify DNA segments that differ in the nature of the evolutionary process.     

Evidence for widespread reticulate evolution within human duplicons

Approximately 5% of the human genome consists of segmental duplications that can cause genomic mutations and may play a role in gene innovation. Reticulate evolutionary processes, such as unequal crossing-over and gene conversion, are known to occur within specific duplicon families, but the broader contribution of these processes to the evolution of human duplications remains poorly characterized. Here, we use phylogenetic profiling to analyze multiple alignments of 24 human duplicon families that span >8 Mb of DNA. Our results indicate that none of them are evolving independently, with all alignments showing sharp discontinuities in phylogenetic signal consistent with reticulation. To analyze these results in more detail, we have developed a quartet method that estimates the relative contribution of nucleotide substitution and reticulate processes to sequence evolution. Our data indicate that most of the duplications show a highly significant excess of sites consistent with reticulate evolution, compared with the number expected by nucleotide substitution alone, with 15 of 30 alignments showing a >20-fold excess over that expected. Using permutation tests, we also show that at least 5% of the total sequence shares 100% sequence identity because of reticulation, a figure that includes 74 independent tracts of perfect identity >2 kb in length. Furthermore, analysis of a subset of alignments indicates that the density of reticulation events is as high as 1 every 4 kb. These results indicate that phylogenetic relationships within recently duplicated human DNA can be rapidly disrupted by reticulate evolution. This finding has important implications for efforts to finish the human genome sequence, complicates comparative sequence analysis of duplicon families, and could profoundly influence the tempo of gene-family evolution.     

Substitution bias, rapid saturation, and the use of mtDNA for nematode systematics

Only relatively recently have researchers turned to molecular methods for nematode phylogeny reconstruction. Thus, we lack the extensive literature on evolutionary patterns and phylogenetic usefulness of different DNA regions for nematodes that exists for other taxa. Here, we examine the usefulness of mtDNA for nematode phylogeny reconstruction and provide data that can be used for a priori character weighting or for parameter specification in models of sequence evolution. We estimated the substitution pattern for the mitochondrial ND4 gene from intraspecific comparisons in four species of parasitic nematodes from the family Trichostrongylidae (38-50 sequences per species). The resulting pattern suggests a strong mutational bias toward A and T, and a lower transition/transversion ratio than is typically observed in other taxa. We also present information on the relative rates of substitution at first, second, and third codon positions and on relative rates of saturation of different types of substitutions in comparisons ranging from intraspecific to interordinal. Silent sites saturate extremely quickly, presumably owing to the substitution bias and, perhaps, to an accelerated mutation rate. Results emphasize the importance of using only the most closely related sequences in order to infer patterns of substitution accurately for nematodes or for other taxa having strongly composition-biased DNA. ND4 also shows high amino acid polymorphism at both the intra- and interspecific levels, and in higher level comparisons, there is evidence of saturation at variable amino acid sites. In general, we recommend using mtDNA coding genes only for phylogenetics of relatively closely related nematode species and, even then, using only nonsynonymous substitutions and the more conserved mitochondrial genes (e.g., cytochrome oxidases). On the other hand, the high substitution rate in genes such as ND4 should make them excellent for population genetics studies, identifying cryptic species, and resolving relationships among closely related congeners when other markers show insufficient variation.      

Utility of nuclear allele networks for the analysis of closely related species in the genus Carabus, subgenus Ohomopterus

Nuclear DNA sequence data for diploid organisms are potentially a rich source of phylogenetic information for disentangling the evolutionary relationships of closely related organisms, but present special phylogenetic problems owing to difficulties arising from heterozygosity and recombination. We analyzed allelic relationships for two nuclear gene regions (phosphoenolpyruvate carboxykinase and elongation factor-1a), along with a mitochondrial gene region (NADH dehydrogenase subunit 5), for an assemblage of closely related species of carabid beetles (Carabus subgenus Ohomopterus). We used a network approach to examine whether the nuclear gene sequences provide substantial phylogenetic information on species relationships and evolutionary history. The mitochondrial gene genealogy strongly contradicted the morphological species boundary as a result of introgression of heterospecific mitochondria. Two nuclear gene regions showed high allelic diversity within species, and this diversity was partially attributable to recombination between various alleles and high variability in the intron region. Shared nuclear alleles among species were rare and were considered to represent shared ancestral polymorphism. Despite the presence of recombination, nuclear allelic networks recovered species monophyly more often and presented genetic differentiation patterns (low to high) among species more clearly. Overall, nuclear gene networks provide clear evidence for separate biological species and information on the phylogenetic relationships among closely related carabid beetles.