The partition matrix: exploring variable phylogenetic signals along nucleotide sequence alignments

Gene Genealogies Reveal High Nucleotide Diversity and Admixture Haplotypes Within Three Alternaria Species Associated with Tomato and Potato

Early blight (EB) and leaf blight are two destructive diseases of tomato in North Carolina (NC), caused by Alternaria linariae and A. alternata, respectively. During the last decade, EB caused by A. solani has increased in potato-producing areas in Wisconsin (WI). We collected 152 isolates of three Alternaria spp. associated with tomato and potato in NC and WI and used the gene genealogical approach to compare the genetic relationships among them. Two nuclear genes: the glyceraldehyde-3-phosphate dehydrogenase (GPDH), RNA polymerase second largest subunit (RPB2), and the rDNA internal transcribed spacer (ITS) region of these isolates were sequenced. Besides, sequences of the GPDH locus from international isolates described in previous studies were included for comparison purposes. A set of single nucleotide polymorphisms was assembled to identify locus-specific and species-specific haplotypes. Nucleotide diversity varied among gene sequences and species analyzed. For example, the estimates of nucleotide diversity and Watterson's theta were higher in A. alternata than in A. linariae and A. solani. There was little or no polymorphisms in the ITS sequences and thus restricted haplotype placement. The RPB2 sequences were less informative to detect haplotype diversity in A. linariae and A. solani, yet six haplotypes were detected in A. alternata. The GPDH sequences enabled strongly supported phylogenetic inferences with the highest haplotype diversity and belonged to five haplotypes (AaH1 to AaH5), which consisted of only A. alternata from NC. However, 13 haplotypes were identified within and among A. linariae and A. solani sequences. Among them, six (AsAlH1 to AsAlH6) were identical to previously reported haplotypes in global samples and the remaining were new haplotypes. The most divergent haplotypes were AaH1, AsAlH2/AsAlH3, and AsAlH4 and consisted exclusively of A. alternata, A. linariae, and A. solani, respectively. Neutrality tests suggested an excess of mutations and population expansion, and selection may play an important role in nucleotide diversity of Alternaria spp.

Frequent nonallelic gene conversion on the human lineage and its effect on the divergence of gene duplicates

Gene conversion is the copying of a genetic sequence from a "donor" region to an "acceptor." In nonallelic gene conversion (NAGC), the donor and the acceptor are at distinct genetic loci. Despite the role NAGC plays in various genetic diseases and the concerted evolution of gene families, the parameters that govern NAGC are not well characterized. Here, we survey duplicate gene families and identify converted tracts in 46% of them. These conversions reflect a large GC bias of NAGC. We develop a sequence evolution model that leverages substantially more information in duplicate sequences than used by previous methods and use it to estimate the parameters that govern NAGC in humans: a mean converted tract length of 250 bp and a probability of [Formula: see text] per generation for a nucleotide to be converted (an order of magnitude higher than the point mutation rate). Despite this high baseline rate, we show that NAGC slows down as duplicate sequences diverge-until an eventual "escape" of the sequences from its influence. As a result, NAGC has a small average effect on the sequence divergence of duplicates. This work improves our understanding of the NAGC mechanism and the role that it plays in the evolution of gene duplicates.

The Rate and Tract Length of Gene Conversion between Duplicated Genes

Interlocus gene conversion occurs such that a certain length of DNA fragment is non-reciprocally transferred (copied and pasted) between paralogous regions. To understand the rate and tract length of gene conversion, there are two major approaches. One is based on mutation-accumulation experiments, and the other uses natural DNA sequence variation. In this review, we overview the two major approaches and discuss their advantages and disadvantages. In addition, to demonstrate the importance of statistical analysis of empirical and evolutionary data for estimating tract length, we apply a maximum likelihood method to several data sets.

Population genetic analysis of Tomato spotted wilt virus on peanut in North Carolina and Virginia

Exploring the genetic diversity and evolutionary history of plant viruses is critical to understanding their ecology and epidemiology. In this study, maximum-likelihood and population genetics-based methods were used to investigate the population structure, genetic diversity, and sources of genetic variation in field isolates of Tomato spotted wilt virus (TSWV) from peanut in North Carolina and Virginia. Selected regions of the nucleocapsid, movement, and RNA-dependent RNA polymerase genes were amplified and sequenced to identify haplotypes and infer genetic relationships between isolates of TSWV with heuristic methods. The haplotype structure of each locus consisted of 1 or 2 predominant haplotypes and >100 haplotypes represented by a single isolate. No specific haplotypes were associated with geographic area, peanut cultivar, or year of isolation. The population was panmictic at the regional level and high levels of genetic diversity were observed among isolates. There was evidence for positive selection on single amino acids in each gene on a background of predominant purifying selection acting upon each locus. The results of compatibility analyses and the persistence of specific gene sequences in isolates collected over three field seasons suggest that recombination was occurring in the population. Estimates of the population mutation rate suggest that mutation has had a significant effect on the shaping of this population and, together with purifying selection, these forces have been the predominant evolutionary forces influencing the TSWV population in peanut in North Carolina and Virginia.

A pattern analysis of gene conversion literature

Gene conversion is an important biological process that involves the transfer of genetic (sequence) information from one gene to another. This can have a variety of effects on an organism, both short-term and long-term and both positive and detrimental. In an effort to better understand this process, we searched through over 3,000 abstracts that contain research on gene conversions, tagging the important data and performing an analysis on what we extract. Through this we established trends that give a better insight into gene conversion research and genetic research in general. Our results show the importance of the process and the importance of continuing gene conversion research.

The power of the methods for detecting interlocus gene conversion

Interlocus gene conversion can homogenize DNA sequences of duplicated regions with high homology. Such nonvertical events sometimes cause a misleading evolutionary interpretation of data when the effect of gene conversion is ignored. To avoid this problem, it is crucial to test the data for the presence of gene conversion. Here, we performed extensive simulations to compare four major methods to detect gene conversion. One might expect that the power increases with increase of the gene conversion rate. However, we found this is true for only two methods. For the other two, limited power is expected when gene conversion is too frequent. We suggest using multiple methods to minimize the chance of missing the footprint of gene conversion.

Lack of global population genetic differentiation in the arbuscular mycorrhizal fungus Glomus mosseae suggests a recent range expansion which may have coincided with the spread of agriculture

The arbuscular mycorrhizal fungus Glomus mosseae is commonly found in agricultural fields. The cosmopolitan species is found in Africa, Europe, America, Asia and Australia. Three hypotheses may explain this worldwide distribution: First, speciation occurred before the continents separated 120 Ma; second, the distribution is a result of human-mediated dispersal related to agriculture and finally, the morphologically defined species may encompass several local endemic species. To test these hypotheses, three genes were sequenced from 82 isolates of G. mosseae originating from six continents and the resulting sequences analysed for geographical subdivision and estimation of migration between continents. Coalescent analyses estimated divergence and age of mutations. Bayesian coalescent modelling was used to reveal important past population changes in the global population. The sequence data showed no geographical structure, with identical genotypes found on different continents. Coalescence analyses indicated a recent diversification in the species, and the data could be explained by a recent population expansion in G. mosseae. The results of this study suggest that speciation and the range expansion happened much later than continental spread and that human activity may have had a major impact on the dispersal and the population structure of the fungus.

Molecular diversity and evolutionary processes of Alternaria solani in Brazil inferred using genealogical and coalescent approaches

Alternaria spp. form a heterogeneous group of saprophytic and plant-pathogenic fungi widespread in temperate and tropical regions. However, the relationship between evolutionary processes and genetic diversity with epidemics is unknown for several plant-pathogenic Alternaria spp. The interaction of Alternaria solani populations with potato and tomato plants is an interesting case study for addressing questions related to molecular evolution of an asexual fungus. Gene genealogies based on the coalescent process were used to infer evolutionary processes that shape the A. solani population. Sequences of the rDNA internal transcribed spacer (ITS) region and the genes which encode the allergenic protein alt a 1 (Alt a 1) and glyceraldehyde-3-phosphate dehydrogenase (Gpd) were used to estimate haplotype and nucleotide diversity as well as for the coalescent analyses. The highest number of parsimony informative sites (n = 14), nucleotide diversity (0.007), and the average number of nucleotide differences (3.20) were obtained for Alt a 1. Although the highest number of haplotypes (n = 7) was generated for ITS, haplotype diversity was the lowest (0.148) for this region. Recombination was not detected. Subdivision was inferred from populations associated with hosts but there was no evidence of geographic subdivision, and gene flow is occurring among subpopulations. In the analysis of the Alt a 1, balancing selection and population expansion or purifying selection could have occurred in A. solani subpopulations associated with potato and tomato plants, respectively. There is strong evidence that the subpopulation of A. solani that causes early blight in potato is genetically distinct from the subpopulation that causes early blight in tomato. The population of A. solani is clonal, and gene flow and mutation are the main evolutionary processes shaping its genetic structure.

Sexy gene conversions: locating gene conversions on the X-chromosome

Gene conversion can have a profound impact on both the short- and long-term evolution of genes and genomes. Here, we examined the gene families that are located on the X-chromosomes of human (Homo sapiens), chimpanzee (Pan troglodytes), mouse (Mus musculus) and rat (Rattus norvegicus) for evidence of gene conversion. We identified seven gene families (WD repeat protein family, Ferritin Heavy Chain family, RAS-related Protein RAB-40 family, Diphosphoinositol polyphosphate phosphohydrolase family, Transcription Elongation Factor A family, LDOC1-related family, Zinc Finger Protein ZIC, and GLI family) that show evidence of gene conversion. Through phylogenetic analyses and synteny evidence, we show that gene conversion has played an important role in the evolution of these gene families and that gene conversion has occurred independently in both primates and rodents. Comparing the results with those of two gene conversion prediction programs (GENECONV and Partimatrix), we found that both GENECONV and Partimatrix have very high false negative rates (i.e. failed to predict gene conversions), which leads to many undetected gene conversions. The combination of phylogenetic analyses with physical synteny evidence exhibits high resolution in the detection of gene conversions.

Evolution of suiform aromatases: ancestral duplication with conservation of tissue-specific expression in the collared peccary (Pecari tayassu)

Aromatase cytochrome P450 (P450arom), the enzyme that catalyzes estrogen synthesis, is required for successful reproduction and is encoded by a single copy gene (CYP19) in most mammals. However, pigs and their distant suiform relatives the peccaries experienced CYP19 duplication. Here, the evolutionary origin of CYP19 duplication, and the evolution of the gene paralogs, was explored further in collared peccaries (Pecari tayassu). Exons IV and V, and the intervening intron, representing duplicated CYP19 genes, were cloned and sequenced from collared peccary, pig, and hippopotamus. Sequence alignment and analysis identified a gene conversion in collared peccary with a breakpoint 102 base pairs (bp) upstream of exon V. Phylogenetic analyses of nucleotide and amino acid sequence upstream of the breakpoint supported a tree in which one peccary sequence was orthologous with the porcine gonadal gene. Cloning and sequencing of tissue transcripts, using reverse-transcriptase polymerase chain reaction techniques (RT-PCR), confirmed that the gonadal ortholog was expressed in collared peccary testis. Orthology of the other genomic sequence with the porcine placental gene was not resolved, but its placenta-specific expression in collared peccary was confirmed by similar transcript analysis. Immunoblot and enzyme activity in collared peccary testes demonstrated much lower levels of P450arom than in pig testis. Collared peccary placental P450arom expression also seemed much lower than pigs. Thus, suiform CYP19 genes arose from an ancestral duplication that has maintained gonad- and placenta-specific expression, but at lower levels in peccaries than pigs, perhaps facilitating the emergence of different reproductive strategies as Suiformes diverged and evolved.

PHYLOGENETIC EVIDENCE FOR MULTIPLE SYMPATRIC ECOLOGICAL DIVERSIFICATION IN A MARINE SNAIL

Parallel speciation can occur when traits determining reproductive isolation evolve independently in different populations that experience a similar range of environments. However, a common problem in studies of parallel evolution is to distinguish this hypothesis from an alternative one in which different ecotypes arose only once in allopatry and now share a sympatric scenario with substantial gene flow between them. Here we show that the combination of a phylogenetic approach with life-history data is able to disentangle both hypotheses in the case of the intertidal marine snail Littorina saxatilis on the rocky shores of Galicia in northwestern Spain. In this system, numerous phenotypic and genetic differences have evolved between two sympatric ecotypes spanning a sharp ecological gradient, and as aside effect of the former have produced partial reproductive isolation. A mitochondrial phylogeny of these populations strongly suggests that the two sympatric ecotypes have originated independently several times. Building upon earlier work demonstrating size-based assortative mating as the main contributor to reproductive isolation among ecotypes, our analysis provides strong evidence that divergent selection across a sharp ecological gradient promoted the parallel divergence of body size and shape between two sympatric ecotypes. Thus, divergent selection occurring independently in different populations has produced the marine equivalent of host races, which may represent the first step in speciation.

New Populations of Sclerotinia sclerotiorum from Lettuce in California and Peas and Lentils in Washington

ABSTRACT Four populations of Sclerotinia sclerotiorum in North America were inferred previously, based on analyses of both rapidly evolving markers (DNA fingerprint and mycelial compatiblity), and multilocus DNA sequence spanning the range between fast and slow evolution. Each population was defined as an interbreeding unit of conspecific individuals sharing a common recent ancestor and arising in a unique evolutionary event. The present study applies this standard to extend characterization of S. sclerotiorum populations to the Western United States. Isolates of S. sclerotiorum (N = 294) were determined to represent three genetically differentiated populations: California (CA, lettuce), Washington (WA, pea/lentil), and Ontario (ON, lettuce). CA was the most diverse population yet sampled in North America. Clonality was detected in ON and WA. No DNA fingerprints were common among the populations. The index of association (I(A)), based on fingerprint, was closer to zero (0) for CA than it was for the other populations. High diversity and lack of association of markers in California are consistent either with genetic exchange and recombination, or with large population size and high standing genetic variation. Intra- and interlocus conflict among three DNA sequence loci was consistent with recombination. The coalescent IGS genealogy confirmed subdivision and showed CA to be older than WA or ON. The Nearest Neighbor statistic on combined data confirmed subdivision among all present and previously defined populations. All isolates had both MAT1-1 and MAT1-2, consistent with uniform homothallism.

Reassortment and concerted evolution in banana bunchy top virus genomes

The nanovirus Banana bunchy top virus (BBTV) has six standard components in its genome and occasionally contains components encoding additional Rep (replication initiation protein) genes. Phylogenetic network analysis of coding sequences of DNA 1 and 3 confirmed the two major groups of BBTV, a Pacific and an Asian group, but show evidence of web-like phylogenies for some genes. Phylogenetic analysis of 102 major common regions (CR-Ms) from all six components showed a possible concerted evolution within the Pacific group, which is likely due to recombination in this region. The CR-M of additional Rep genes is close to that of DNA 1 and 2. Comparison of tree topologies constructed with DNA 1 and DNA 3 coding sequences of 14 BBTV isolates showed distinct phylogenetic histories based on Kishino-Hasegawa and Shimodaira-Hasegawa tests. The results of principal component analysis of amino acid and codon usages indicate that DNA 1 and 3 have a codon bias different from that of all other genes of nanoviruses, including all currently known additional Rep genes of BBTV, which suggests a possible ancient genome reassortment event between distinctive nanoviruses.

Codon-usage bias versus gene conversion in the evolution of yeast duplicate genes

Many Saccharomyces cerevisiae duplicate genes that were derived from an ancient whole-genome duplication (WGD) unexpectedly show a small synonymous divergence (K(S)), a higher sequence similarity to each other than to orthologues in Saccharomyces bayanus, or slow evolution compared with the orthologue in Kluyveromyces waltii, a non-WGD species. This decelerated evolution was attributed to gene conversion between duplicates. Using approximately 300 WGD gene pairs in four species and their orthologues in non-WGD species, we show that codon-usage bias and protein-sequence conservation are two important causes for decelerated evolution of duplicate genes, whereas gene conversion is effective only in the presence of strong codon-usage bias or protein-sequence conservation. Furthermore, we find that change in mutation pattern or in tDNA copy number changed codon-usage bias and increased the K(S) distance between K. waltii and S. cerevisiae. Intriguingly, some proteins showed fast evolution before the radiation of WGD species but little or no sequence divergence between orthologues and paralogues thereafter, indicating that functional conservation after the radiation may also be responsible for decelerated evolution in duplicates.

Detecting recombination in evolving nucleotide sequences

BACKGROUND

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. These recombination events can be obscured by subsequent residue substitutions, which consequently complicate their detection. While there are many algorithms for the identification of recombination events, little is known about the effects of subsequent substitutions on the accuracy of available recombination-detection approaches.

RESULTS

We assessed the effect of subsequent substitutions on the detection of simulated recombination events within sets of four nucleotide sequences under a homogeneous evolutionary model. The amount of subsequent substitutions per site, prior evolutionary history of the sequences, and reciprocality or non-reciprocality of the recombination event all affected the accuracy of the recombination-detecting programs examined. Bayesian phylogenetic-based approaches showed high accuracy in detecting evidence of recombination event and in identifying recombination breakpoints. These approaches were less sensitive to parameter settings than other methods we tested, making them easier to apply to various data sets in a consistent manner.

CONCLUSION

Post-recombination substitutions tend to diminish the predictive accuracy of recombination-detecting programs. The best method for detecting recombined regions is not necessarily the most accurate in identifying recombination breakpoints. For difficult detection problems involving highly divergent sequences or large data sets, different types of approach can be run in succession to increase efficiency, and can potentially yield better predictive accuracy than any single method used in isolation.

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.

DNA sequence-based subtyping and evolutionary analysis of selected Salmonella enterica serotypes

While serotyping and phage typing have been used widely to characterize Salmonella isolates, sensitive subtyping methods that allow for evolutionary analyses are essential for examining Salmonella transmission, ecology, and evolution. A set of 25 Salmonella enterica isolates, representing five clinically relevant serotypes (serotypes Agona, Heidelberg, Schwarzengrund, Typhimurium, and Typhimurium var. Copenhagen) was initially used to develop a multilocus sequence typing (MLST) scheme for Salmonella targeting seven housekeeping and virulence genes (panB, fimA, aceK, mdh, icdA, manB, and spaN). A total of eight MLST types were found among the 25 isolates sequenced. A good correlation between MLST types and Salmonella serotypes was observed; only one serotype Typhimurium var. Copenhagen isolate displayed an MLST type otherwise typical for serotype Typhimurium isolates. Since manB, fimA, and mdh allowed for the highest subtype discrimination among the initial 25 isolates, we chose these three genes to perform DNA sequencing of an additional 41 Salmonella isolates representing a larger diversity of serotypes. This "three-gene sequence typing scheme" allowed discrimination of 25 sequence types (STs) among a total of 66 isolates; STs correlated well with serotypes and allowed within-serotype differentiation for 9 of the 12 serotypes characterized. Phylogenetic analyses showed that serotypes Kentucky and Newport could each be separated into two distinct, statistically well supported evolutionary lineages. Our results show that a three-gene sequence typing scheme allows for accurate serotype prediction and for limited subtype discrimination among clinically relevant serotypes of Salmonella. Three-gene sequence typing also supports the notion that Salmonella serotypes represent both monophyletic and polyphyletic lineages.

Gene conversion and functional divergence in the beta-globin gene family

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication-degeneration-complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d(N)/ d(S) (omega) rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the beta-globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate beta-globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohta's model and the DDC model by examining selective pressures along lineages in a phylogeny of beta-globin genes in eutherian mammals. We did not find significant evidence for an increase in the omega ratio following major duplication events in this family. However, one exception to this pattern was the duplication of gamma-globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of beta-globin genes evolved under a nonepisodic process of functional divergence.

Widespread Recombination in Published Animal mtDNA Sequences1

Mitochondrial DNA (mtDNA) recombination has been observed in several animal species, but there are doubts as to whether it is common or only occurs under special circumstances. Animal mtDNA sequences retrieved from public databases were unambiguously aligned and rigorously tested for evidence of recombination. At least 30 recombination events were detected among 186 alignments examined. Recombinant sequences were found in invertebrates and vertebrates, including primates. It appears that mtDNA recombination may occur regularly in the animal cell but rarely produces new haplotypes because of homoplasmy. Common animal mtDNA recombination would necessitate a reexamination of phylogenetic and biohistorical inference based on the assumption of clonal mtDNA transmission. Recombination may also have an important role in producing and purging mtDNA mutations and thus in mtDNA-based diseases and senescence.

Recombination and migration of Cryphonectria hypovirus 1 as inferred from gene genealogies and the coalescent

Genealogy-based methods were used to estimate migration of the fungal virus Cryphonectria hypovirus 1 between vegetative compatibility types of the host fungus, Cryphonectria parasitica, as a means of estimating horizontal transmission within two host populations. Vegetative incompatibility is a self/non-self recognition system that inhibits virus transmission under laboratory conditions but its effect on transmission in nature has not been clearly demonstrated. Recombination within and among different loci in the virus genome restricted the genealogical analyses to haplotypes with common mutation and recombinational histories. The existence of recombination necessitated that we also use genealogical approaches that can take advantage of both the mutation and recombinational histories of the sample. Virus migration between populations was significantly restricted. In contrast, estimates of migration between vegetative compatibility types were relatively high within populations despite previous evidence that transmission in the laboratory was restricted. The discordance between laboratory estimates and migration estimates from natural populations highlights the challenges in estimating pathogen transmission rates. Genealogical analyses inferred migration patterns throughout the entire coalescent history of one viral region in natural populations and not just recent patterns of migration or laboratory transmission. This application of genealogical analyses provides markedly stronger inferences on overall transmission rates than laboratory estimates do.

Tracing the Decay of the Historical Signal in Biological Sequence Data

Alignments of nucleotide or amino acid sequences may contain a variety of different signals, one of which is the historical signal that we often try to recover by phylogenetic analysis. Other signals, such as those arising due to compositional heterogeneities, among-lineage and among-site rate heterogeneities, invariant sites, and covariotides, may interfere adversely with the recovery of the historical signal. The effect of the interaction of these signals on phylogenetic inference is not well understood and may, in many cases, even be underappreciated. In this study, we investigate this matter and present results based on Monte Carlo simulations. We explored the success of four phylogenetic methods in recovering the true tree from data that had evolved under conditions where the equilibrium base frequencies and substitution rates were allowed to vary among lineages. Seven scenarios with increasingly complex conditions were investigated. All of the methods tested, with the exception of neighbor-joining using LogDet distances, were sensitive to compositional convergence in nonsister lineages. Maximum parsimony was also susceptible to attraction between long edges. In many cases, however, phylogenetic inference methods can still recover the true tree when misleading signals are present, in some instances even when the historical signal is no longer dominant. These results highlight the growing need for simple methods to detect violation of the phylogenetic assumptions.

Recombination and Migration of Cryphonectria hypovirus 1 as Inferred From Gene Genealogies and the Coalescent

Genealogy-based methods were used to estimate migration of the fungal virus Cryphonectria hypovirus 1 between vegetative compatibility types of the host fungus, Cryphonectria parasitica, as a means of estimating horizontal transmission within two host populations. Vegetative incompatibility is a self/non-self recognition system that inhibits virus transmission under laboratory conditions but its effect on transmission in nature has not been clearly demonstrated. Recombination within and among different loci in the virus genome restricted the genealogical analyses to haplotypes with common mutation and recombinational histories. The existence of recombination necessitated that we also use genealogical approaches that can take advantage of both the mutation and recombinational histories of the sample. Virus migration between populations was significantly restricted. In contrast, estimates of migration between vegetative compatibility types were relatively high within populations despite previous evidence that transmission in the laboratory was restricted. The discordance between laboratory estimates and migration estimates from natural populations highlights the challenges in estimating pathogen transmission rates. Genealogical analyses inferred migration patterns throughout the entire coalescent history of one viral region in natural populations and not just recent patterns of migration or laboratory transmission. This application of genealogical analyses provides markedly stronger inferences on overall transmission rates than laboratory estimates do.

Recombination in evolutionary genomics

Recombination can be a dominant force in shaping genomes and associated phenotypes. To better understand the impact of recombination on genomic evolution, we need to be able to identify recombination in aligned sequences. We review bioinformatic approaches for detecting recombination and measuring recombination rates. We also examine the impact of recombination on the reconstruction of evolutionary histories and the estimation of population genetic parameters. Finally, we review the role of recombination in the evolutionary history of bacteria, viruses, and human mitochondria. We conclude by highlighting a number of areas for future development of tools to help quantify the role of recombination in genomic evolution.

Genomic analysis of the murine odorant receptor MOR28 cluster: a possible role of gene conversion in maintaining the olfactory map

Genomic analysis was performed for the murine odorant receptor (OR) genes. The MOR28 cluster on chromosome 14 was extensively studied. It contains six OR genes, MOR28, 10, 83, 29A, 29B and 30. The human homolog of this cluster is located on the human chromosome 14, and contains five OR genes, HOR28/10, 83, 29A, 29B and 30. Sequence comparison of these OR gene paralogs and orthologs suggests that the coding homologies are accounted for not only by recent gene duplication, but also by gene conversion among the coding sequences within the cluster. A possible role of gene conversion in the olfactory system is discussed in the context of the olfactory map.

Evaluation of methods for detecting recombination from DNA sequences: empirical data

The performance of 14 different recombination detection methods was evaluated by analyzing several empirical data sets where the presence of recombination has been suggested or where recombination is assumed to be absent. In general, recombination methods seem to be more powerful with increasing levels of divergence, but different methods showed distinct performance. Substitution methods using summary statistics gave more accurate inferences than most phylogenetic methods. However, definitive conclusions about the presence of recombination should not be derived on the basis of a single method. Performance patterns observed from the analysis of real data sets coincided very well with previous computer simulation results. Previous recombination inferences from some of the data sets analyzed here should be reconsidered. In particular, recombination in HIV-1 seems to be much more widespread than previously thought. This finding might have serious implications on vaccine development and on the reliability of previous inferences of HIV-1 evolutionary history and dynamics.

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

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

Comparison of chimpanzee and human leukocyte Ig-like receptor genes reveals framework and rapidly evolving genes

The leukocyte receptor complex (LRC) on human chromosome 19 contains related Ig superfamily killer cell Ig-like receptor (KIR) and leukocyte Ig-like receptor (LIR) genes. Previously, we discovered much difference in the KIR genes between humans and chimpanzees, primate species estimated to have approximately 98.8% genomic sequence similarity. Here, the common chimpanzee LIR genes are identified, characterized, and compared with their human counterparts. From screening a chimpanzee splenocyte cDNA library, clones corresponding to nine different chimpanzee LIRs were isolated and sequenced. Analysis of genomic DNA from 48 unrelated chimpanzees showed 42 to have all nine LIR genes, and six animals to lack just one of the genes. In structural diversity and functional type, the chimpanzee LIRs cover the range of human LIRs. Although both species have the same number of inhibitory LIRs, humans have more activating receptors, a trend also seen for KIRs. Four chimpanzee LIRs are clearly orthologs of human LIRs. Five other chimpanzee LIRs have paralogous relationships with clusters of human LIRs and have undergone much recombination. Like the human genes, chimpanzee LIR genes appear to be organized into two duplicated blocks, each block containing two orthologous genes. This organization provides a conserved framework within which there are clusters of faster evolving genes. Human and chimpanzee KIR genes have an analogous arrangement. Whereas both KIR and LIR genes can exhibit greater interspecies differences than the genome average, within each species the LIR gene family is more conserved than the KIR gene family.

A simulation study of the reliability of recombination detection methods

There exist many methods to detect recombination or mosaic structure in a sample of DNA sequences. But how reliable are they? Four methods were investigated with respect to their power to detect recombination in simulated samples with different amounts of recombination and mutation. In addition, we investigated the impact of the shape of the underlying genealogy on their performances. We found that the methods detected far fewer recombinations than were theoretically possible and that methods based on the principle of incompatibility in general had more power than methods that did not make use of this principle explicitly. This seemed, in particular, to be the case for phylogenies generated under population expansion scenarios which result in long branches at the tips and small deep branches. In addition to the results obtained through simulations, a series of new theoretical results on recombination is presented.

Recombinational and mutational hotspots within the human lipoprotein lipase gene

Here an analysis is presented of the roles of recombination and mutation in shaping previously determined haplotype variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase gene (LPL), scored in 71 individuals from three populations: 24 African Americans, 24 Finns, and 23 non-Hispanic whites. Recombination and gene-conversion events inferred from data on 88 haplotypes that were defined by 69 variable sites were tested. The analysis revealed 29 statistically significant recombination events and one gene-conversion event. The recombination events were concentrated in a 1.9-kb region, near the middle of the segment, that contains a microsatellite and a pair of tandem and complementary mononucleotide runs; both the microsatellite and the runs show length variation. An analysis of site variation revealed that 9.6% of the nucleotides at CpG sites were variable, as were 3% of the nucleotides found in mononucleotide runs of >/=5 nucleotides, 3% of the nucleotides found </=3 bp from certain putative polymerase alpha-arrest sites, and 0. 5% of the remaining nucleotides. This nonhomogeneous distribution of variation suggests that multiple mutational hits at certain sites are common, an observation that challenges the fundamental assumption of the infinite-sites-mutation model. The nonrandom patterns of recombination and mutation suggest that randomly chosen single-nucleotide polymorphisms may not be optimal for disequilibrium mapping of this gene. Overall, these results indicate that both recombinational and mutational hotspots have played significant roles in shaping the haplotype variation at the LPL locus.

Possible emergence of new geminiviruses by frequent recombination

Although exchange of genetic information by recombination plays a role in the evolution of viruses, the extent to which it generates diversity is not clear. We analyzed genomes of geminiviruses for recombination using a new statistical procedure developed to detect gene conversions. Geminiviruses (family, Geminiviridae) are a group of plant viruses characterized by a genome of circular single-stranded DNA (approximately 2700 nucleotides in length) encapsidated in twinned quasi-isometric particles. Complete nucleotide sequences of geminiviruses were aligned, and recombination events were detected by searching pairs of viruses for sequences that are significantly more similar than expected based on random distribution of polymorphic sites. The analyses revealed that recombination is very frequent and occurs between species and within and across genera. Tests identified 420 statistically significant recombinant fragments distributed across the genome. The results suggest that recombination is a significant contributor to geminivirus evolution. The high rate of recombination may be contributing to the recent emergence of new geminivirus diseases.

The detection and measurement of recombination from sequence data

There are two types of recombination that we may wish to detect: rare recombinants between members of different populations or species and repeated recombination within a population. Methods appropriate in the former context are inappropriate in the latter because they depend on recognizing the existence of runs of nucleotides with similar ancestry. If recombination is sufficiently frequent, no such runs will be present. Several methods, including the homoplasy test and the incompatibility test, are described that are appropriate for detecting repeated recombination and for measuring its importance, relative to mutation, in causing genetic change. The sensitivity of these tests is investigated by simulating populations with varying frequencies of mutation and recombination and calculating the various statistics on samples.

Haplotype structure and population genetic inferences from nucleotide-sequence variation in human lipoprotein lipase

Allelic variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase gene (LPL) was scored in 71 healthy individuals (142 chromosomes) from three populations: African Americans (24) from Jackson, MS; Finns (24) from North Karelia, Finland; and non-Hispanic Whites (23) from Rochester, MN. The sequences had a total of 88 variable sites, with a nucleotide diversity (site-specific heterozygosity) of .002+/-.001 across this 9.7-kb region. The frequency spectrum of nucleotide variation exhibited a slight excess of heterozygosity, but, in general, the data fit expectations of the infinite-sites model of mutation and genetic drift. Allele-specific PCR helped resolve linkage phases, and a total of 88 distinct haplotypes were identified. For 1,410 (64%) of the 2,211 site pairs, all four possible gametes were present in these haplotypes, reflecting a rich history of past recombination. Despite the strong evidence for recombination, extensive linkage disequilibrium was observed. The number of haplotypes generally is much greater than the number expected under the infinite-sites model, but there was sufficient multisite linkage disequilibrium to reveal two major clades, which appear to be very old. Variation in this region of LPL may depart from the variation expected under a simple, neutral model, owing to complex historical patterns of population founding, drift, selection, and recombination. These data suggest that the design and interpretation of disease-association studies may not be as straightforward as often is assumed.