Gene conversion yields novel gene combinations in paralogs of GOT1 in the copepod Tigriopus californicus
© Willett; licensee BioMed Central Ltd. 2013
Received: 26 March 2013
Accepted: 8 July 2013
Published: 12 July 2013
Gene conversion of duplicated genes can slow the divergence of paralogous copies over time but can also result in other interesting evolutionary patterns. Islands of genetic divergence that persist in the face of gene conversion can point to gene regions undergoing selection for new functions. Novel combinations of genetic variation that differ greatly from the original sequence can result from the transfer of genetic variation between paralogous genes by rare gene conversion events. Genetically divergent populations of the copepod Tigriopus californicus provide an excellent model to look at the patterns of divergence among paralogs across multiple independent evolutionary lineages.
In this study the evolution of a set of paralogous genes encoding putative aspartate transaminase proteins (called GOT1 here) are examined in populations of the copepod T. californicus. One pair of duplicated genes, GOT1p1 and GOT1p2, has regions of high divergence between the copies in the face of apparent on-going gene conversion. The GOT1p2 gene also has unique haplotypes in two populations that appear to have resulted from a transfer of genetic variation via inter-paralog gene conversion. A second pair of duplicated genes GOT1Sr and GOT1Sd also shows evidence of gene conversion, but this gene conversion does not appear to have maintained each as a functional copy in all populations.
The patterns of conservation and sequence divergence across this set of paralogous genes among populations of T. californicus suggest that some interesting evolutionary patterns are occurring at these loci. The results for the GOT1p1/GOT1p2 paralogs illustrate how gene conversion can factor in the creation of a mosaic pattern of regions of high divergence and low divergence. When coupled with rare gene conversion events of divergent regions, this pattern can result in the formation of novel proteins differing substantially from either original protein. The evolutionary patterns across these paralogs show how gene conversion can both constrain and facilitate diversification of genetic sequences.
Gene conversion can impact the evolution of duplicated genes in a number of different ways including both impeding sequence divergence between genes and transferring variation between them . Gene conversion is a common mechanism of unidirectional homologous recombination in eukaryotes that results in a cut-and-paste like copying of sequence between similar alleles that are either at the same locus or at another locus in the same genome (reviewed in Chen et al. ). Concerted evolution can result from loci undergoing repeated gene conversion, which causes duplicated genes to evolve in tandem and not diverge from one another over evolutionary time. Not all duplicated genes are subject to gene conversion, in fact, surveys in mammals and fruit flies suggest that only about ten percent of paralogous copies show signs of gene conversion, and only a small fraction of the total sequence length is typically impacted [3, 4].
Duplicated genes that are experiencing concerted evolution typically will go through a series of phases of differential divergence. Rates of gene conversion between sequences go down as the sequences become more dissimilar. For gene duplicates undergoing some level of concerted evolution, divergence between them will not begin to increase markedly until a threshold of sequence divergence is breached (as high as 20 percent ). Models of this process suggest that there will typically be a long period of evolution with only low levels of divergence until a threshold level of divergence is passed at which point the rate of divergence will increase . Selective divergence can counter this homogenization and lead to the establishment and maintenance of regions of higher sequence divergence in the face of gene conversion if differences in specific regions of the gene between the two duplicates are adaptive (e.g. with neofunctionalization ). Teshima and Innan  propose scanning for this specific pattern as a method of identifying such regions undergoing selection. Using this method in a study in yeast, Takuno and Innan  identified two sets of duplicated heat shock proteins that likely fit this model.
In addition to the role outlined above in slowing or countering adaptive divergence between duplicated genes, gene conversion can also play a role in transferring adaptive variation between duplicate genes. Under such a scenario, gene conversion acts to increase the effective population size of the duplicated genes, making selection more efficient. This transfer can spread advantageous variation and remove deleterious mutations [9, 10]. A number of studies have shown that gene conversion between duplicate genes with some degree of initial divergence between them can result in the introduction of high levels of variation at the converted locus [11–18]. For many of these cases, this variation appears to be adaptive with a number of these genes under selection for higher haplotype diversity (e.g. MHC, attacin, and resistance genes in plants).
The copepod Tigriopus californicus has a set of unique features that makes it useful system in which to look at patterns of molecular evolution in duplicated genes. T. californicus exists in a series of extensively genetically divergent populations that have undergone substantial periods of independent evolution from one another. This species occurs in rocky upper intertidal pools along the Pacific coast of North America from Mexico to Alaska. Populations of this species can be highly genetically divergent from one another even over relatively short distances, with mitochondrial DNA (mtDNA) divergences greater than 20 percent between populations [19–21]. Divergence in the nuclear genome is lower but still substantial, likely reflecting a substantially higher rate of mutation for the mtDNA . Even with these higher rates of mtDNA evolution, the levels of divergence among populations suggest that these populations have been evolving fairly independently of one another for long periods of time. Genomic resources are being developed for this species and now include published transcriptomes from a pair of populations, and these resources facilitate the characterization of paralogs .
In this paper the molecular evolution of a set of aspartate transaminase-encoding homologs is examined in populations of T. californicus. A putatively mitochondrially targeted homolog was previously identified from this species  and named after the corresponding allozyme locus (GOT2, the enzyme aspartate transaminase was formerly called glutamate-oxaloacetate transamine; EC 126.96.36.199). Five additional homologs are described in this paper that have originated from a series of gene duplication events in the evolutionary lineage leading to this species. Sequence similarity suggests that these genes are likely to be cytoplasmically targeted GOT1 proteins. Two sets of somewhat more recently duplicated pairs of genes show strong evidence of gene conversion. In this paper the differential impact of gene conversion on the evolution of these two pairs of duplicated genes is examined.
Identification of GOT paralogs
Fixed genetic divergence in coding regions for orthologs and paralogs of GOT1 in populations of T. californicus
Between GOT1p1/GOT1p2 paralogs
SD p1/SD p2
LJS p1/LJS p2
AB p1/AB p2
SCN p1/SCN p2
SD p1/LJS p1
SD p1/AB p1
SD p1/SCN p1
SD p2/LJS p2
SD p2/AB p2
SD p2/SCN p2
Between GOT1Sd/GOT1Sr paralogs
SD Sd/SD Sr
LJS Sd/LJS Sr
AB Sd/AB Sr
SCN Sd/SCN Sr
SD Sd/LJS Sd
SD Sd/AB Sd
SD Sd/SCN Sd
SD Sr/LJS Sr
SD Sr/AB Sr
SD Sr/SCN Sr
SD 6a/LJS 6a
SD 6a/AB 6a
SD 6a/SCN 6a
LJS 6a/AB 6a
LJS 6a/SCN 6a
AB 6a/SCN 6a
In addition to the large amount of amino acid divergence among the more divergent GOT1 paralogs, there are also a number of structural differences at these loci. The GOT1_6a gene and the GOT1Sr genes each have four introns in the same positions in the gene (as assessed by their position in the amino acid alignment). GOT2 also has four introns but only one of these shares a position with those of the GOT1_6a and GOT1Sr genes (the third intron). The size of this third intron varies widely from 152 bp in GOT1_6a to 3894 bp for the SCN population for GOT1Sr (the other three populations each have a 2723 bp for this intron in the GOT1Sr gene). Interestingly, the GOT1p1/GOT1p2 genes lack introns completely. The transcript for each of these genes is between 1257 bp for GOT1Sr and 1532 bp for GOT1_6a, while the coding regions are all close to 1224 bp (with GOT1_6a being 1233 bp). We did not obtain sequence corresponding to the first 618 bp of the coding region for the GOT1Sd gene, but the sequenced portion is consistent with the presence of the final three introns. For the AB and LJS populations the second and third introns respectively have polymorphisms that would alter the predicted splice sequences for the GOT1Sd gene.
Five of these six GOT genes can be found in the published transcriptome dataset derived from the SD and SCN populations of T. californicus with only the GOT1Sd gene missing. These data also give some hints as to the relative expression levels of these genes. Total read numbers per gene are somewhat low overall in this 454 dataset but the highest counts were found for the GOT2 gene with 247 reads and the GOT1p1/GOT1p2 genes with 154 reads summed over both copies. Examination of the proportion of reads from the diagnostic regions of the GOT1p1/GOT1p2 genes suggests that the expression of the GOT1p1 gene is about 6-fold higher than that of the GOT1p2 gene. The GOT1_6a and GOT1Sr genes had fewer than 10 reads each suggesting that they are expressed at a much lower level. Consistent with its absence from the transcriptomes, our lab found no expression of the GOT1Sd gene using qualitative RT-PCR assays in the San Diego (SD) population, but we did find expression of both the GOT1Sr and GOT1p1/GOT1p2 paralogs (Willett CS, unpublished data). Sequences of mRNA obtained from individual copepods from these experiments were identical to haplotypes obtained via direct sequencing from the coding regions.
Divergence in GOT1 paralogs and gene conversion
Polymorphism capture via gene conversion
Levels of polymorphism in GOT paralogs in T. californicus
Indels in coding region
70 bp insert (one haplotype)
(stop codon poly.)
4bp poly., 1bp poly., 4bp fixed
Although there is evidence for inter-paralog gene conversion for the GOT1Sd and GOT1Sr paralogs as well, it does not appear to have been substantial enough to result in both copies retaining their open reading frames in all haplotypes. For the SD, LJS, and particularly SCN populations there are fixed and polymorphic indels in exons in GOT1Sd that should disrupt the reading frame and result in greatly truncated mRNAs (Table 2). In the AB population there appears to be a premature stop codon in the GOT1Sr that is polymorphic in this population. For both the GOT1Sr and GOT1Sd paralogs elevated ka/ks ratios are seen for some comparisons further suggesting reduced functional constraint (Table 1). For the GOT1p2 gene one haplotype in the SD population also had an insertion that would disrupt the reading frame suggesting that non-functional alleles can also be found at this locus. A one bp deletion was found in the coding region for a single haplotype in GOT1_6a in the SD population as well. Only for the GOT1p1 gene copy were no such truncating or frameshift polymorphisms found in any of this set of four populations of T. californicus for these five GOT1 homologs.
I have identified a set of homologous genes from T. californicus that appear to encode aspartate transaminase proteins and these genes display a number of interesting patterns of inter-locus gene conversion. In discussing these results, first, I will discuss the potential deeper level relationships among these duplicates within and between species and then, second, I will look at the interesting patterns of gene conversion in two pairs of more closely related duplicates.
The cytosolic GOT1 proteins have undergone a number of gene duplication events in copepods and in the T. californicus lineage. The GOT1p1/GOT1p2 paralogs cluster phylogenetically with cytosolic GOT1 proteins in other species of arthropods and are their most likely orthologs. The relationships of the other three GOT1 paralogs to other GOT1 proteins are not resolved with the exception of a weakly supported relationship to putative GOT1 paralogs in two other distantly related copepod species (Caligus clemensi and Lepeoophtheirus salmonis). The lack of deeply divergent GOT1 paralogs in other sequenced metazoan genomes suggests that the duplication events producing the GOT1_6a and GOT1Sd/GOT1Sr paralogs may have occurred within copepods and were not the result of an ancient metazoan duplication event. Other examples of older duplicates of aspartate transaminases in animals are restricted to individual clades such as mammals as can be seen in panther gene family trees http://www.pantherdb.org/ for aspartate aminotransferases . If the duplications did occur within copepods, perhaps relatively high levels of amino acid divergence in these paralogs are obscuring their relationship to the other GOT1 proteins. Regardless of the deeper level relationships, it is clear that the duplications that have resulted in the production of the GOT1Sd/GOT1Sr and GOT1p1/GOT1p2 gene pairs occurred more recently than these deeper splits. Most likely these splits occurred in the common ancestor of these four populations of T. californicus given the presence of each copy in each population.
The GOT1p1 paralog is the most conserved of the five paralogs with no evidence for segregating non-functional alleles (Table 2) and it has the highest levels of constraint as measured by ka/ks values (Table 1). The higher expression level of the GOT1p1 copy, coupled with potential matches between predicted amino acid differences and allozyme allele differences among populations together suggest that the GOT1p1 paralog could be the same locus as the GOT1 allozyme used previously to examine genetic variation among T. californicus populations [19, 28, 29] and may be the primary cytosolic aspartate transaminase protein in this species. The GOT1p2 paralog has slightly lower levels of constraint than the GOT1p1 paralog and has one haplotype that contains a frameshift polymorphism in this sample of sequences from the SD population. Of the five paralogs, the GOT1Sd gene is behaving the most like a pseudogene. It does not appear to be expressed at detectable levels and has a series of frameshift substitutions in each of the populations that disrupt the reading frame (with the exception of the AB population).
Turning now to the patterns of gene conversion in the more recently duplicated pairs of paralogs, GOT1Sd/GOT1Sr and GOT1p1/GOT1p2, it is clear that there has been gene conversion in the past within each pair. There is no evidence of gene conversion between the more divergent paralogs, e.g. between GOT1Sr and GOT1_6a. There are numerous likely gene conversion tracks resulting from both inter- and intra-locus events between pairs for both of these sets of paralogs (Figures 1 and 2; Additional file 3: Table S2). For the GOT1Sd/GOT1Sr pair the inter-paralog gene conversion events are largely restricted to the exonic sequences with a large intron becoming largely un-alignable between paralogs. The GOT1Sd gene appears to be evolving as a pseudogene in several populations as discussed above despite evidence for inter-locus gene conversion events with the largely intact GOT1Sr gene. Apparently these gene conversion events are not happening frequently enough to maintain the open reading frame of this GOT1Sd copy in all populations. In contrast to the GOT1Sd/GOT1Sr pair of genes, there are no introns in the coding sequences of the GOT1p1/GOT1p2 paralogs and the regions of elevated divergence between the two paralogs are therefore located within the single exon. Close physical proximity in the genome can facilitate interlocus gene conversion  and in fact, the GOT1p1/GOT1p2 paralogs are tightly linked (and are also located on the same chromosome as GOT2; Willett CS, unpublished data). The allozyme loci GOT1 and GOT2 were previously shown to be linked , lending further credence to the idea that the GOT1p1 and/or GOT1p2 loci might encode the allozyme marker GOT1 that has been previously characterized in this species.
Both pairs of paralogs GOT1Sd/GOT1Sr and GOT1p1/GOT1p2 show islands of genetic divergence amid regions of higher similarity but the evolutionary explanation for this pattern may differ between the two sets of duplicates. For the GOT1Sd/GOT1Sr pair the divergence is restricted to the introns and may be a result of the accumulation of substitutions that can terminate inter-paralog gene conversion in those stretches of the gene. Divergence in sequence similarity that lowers the level of gene conversion could accumulate either via the gradual accumulation of single-base differences or more rapidly by larger changes such as large indels [32, 33]. The GOT1Sd/GOT1Sr paralogs have both very large size differences and low sequence similarity in the intron so that either mode of divergence could have contributed to the absence of gene conversion in these regions. Even small regions of clustered sequence divergence (with multiple substitutions or indels) can dramatically reduce the rate of gene conversion for a region of a gene [34, 35]. The net result of this divergence for the GOT1Sd/GOT1Sr paralogs is that interlocus gene conversion is not likely to occur in this intronic region of the gene and these regions are free to accumulate further differences.
In contrast for the GOT1p1/GOT1p2 paralogs the regions of genetic divergence occur in the exons and there are no fixed indels in these regions that could disrupt interlocus gene conversion. Teshima and Innan  have suggested that such regions of differentiation in the face of on-going gene conversion can be a signal that selection is maintaining divergence in the paralogs (i.e. the paralogs have begun the process of neofunctionalization). Under such a model the width of the divergent region should extend less than the average length of a gene conversion tract from the selected site or sites. A number of duplicated genes show such islands of divergence that are associated with clear functional differences in the resulting proteins (e.g. RH factor and opsin proteins ). Other duplicated genes in yeast and Drosophila show a similar pattern consistent with selection but lack evidence for functional differences [8, 37]. For the GOT1p1/GOT1p2 paralogs, one potential neutral explanation for this pattern could posit that gene conversion initiation is lower in these regions and that these regions have accumulated enough differentiation to begin to suppress gene conversion. An argument against this limited initiation idea is that intralocus gene conversion is common in the region of sequence differentiation between these two paralogs in the first half of the gene. This observation suggests that sequence factors are not completely suppressing the initiation of gene conversion events in the divergent regions of the gene. Other factors that could also suppress interlocus gene conversion such as indel differences are also absent. The loss of fixed divergences between paralogs for one of these islands of genetic divergence in the SD and LJS populations in the second half of the gene (discussed further below) also argues that gene conversion is still possible for these regions. Although these results are suggestive of a selective explanation, further study attempting to identify functional differences between the GOT1p1/GOT1p2 paralogs is needed to confirm or reject this hypothesis.
A region of high polymorphism and lowered divergence between a set of alleles in the GOT1p1/GOT1p2 paralogs in the SD and LJS populations is likely to have been created by inter-paralog gene conversion. The patterns of variation and phylogenetic evidence (Figure 1 and Figure 4) are consistent with one-way transfers of variation from each population’s GOT1p1 locus to the GOT1p2 locus. One-way exchange like this is consistent with other studies where gene conversion shows biased directionality [2, 38]. The net result of this directional gene conversion is to transfer variants from one paralog to the other. In this case this transfer is limited to the second half of the gene resulting in haplotypes that are a chimera of the GOT1p1 and the GOT1p2 paralogs and this transfer also results in an increase in the levels of polymorphism in this region of the gene. The chimeric protein that results is substantially altered from that produced by other GOTp2 alleles, differing by 8 amino acids, while still differing from GOTp1 by 20 amino acid in the first half of the gene.
It is possible that gene conversion events that result in greatly augmented polymorphism in gene duplicates are effectively neutral, but in a number of other cases they appear to be under selection, often occurring in genes undergoing selection for diversification [11–17, 39]. For the GOT1p1 and GOT1p2 genes there is not a clear signal of diversifying selection in comparisons of orthologous copies across populations with Ka/Ks values much lower than one (Table 2). Without any further functional information it is difficult to say whether the gene conversion events that resulted in greatly increased diversity in the GOT1p2 gene in the SD and LJS population are adaptive in nature. Clearly this process has generated a large amount of novel variation at this locus both in DNA and protein sequence.
The two sets of duplicate genes of GOT1 illustrate different patterns of evolution with ongoing gene conversion among duplicated copies. The set of GOT1Sd/GOT1Sr genes appear to be in the process of diverging with gradually decreasing gene conversion given that one copy does not maintain its open reading frame and does not appear to be expressed. The central intron in this gene is already quite divergent. In contrast for the GOT1p1/GOT1p2 pair, gene conversion is maintaining much higher similarity in some regions of the gene but other exonic portions are substantially diverged. The combination of these islands of genetic divergence between paralogs with rare gene conversion events has the ability to construct radically different haplotypes from the combination of variation in both paralogs (as has happened in the SD and LJS GOT1p2 gene). Further work on the function of these two duplicates could help to determine whether there are likely to be adaptive differences between these copies.
Isolation and sequencing of GOT1 homologs
The putative GOT1 homologs were uncovered from T. californicus using an analogous strategy to that used to obtain the GOT2 homolog in this species . Briefly, a cDNA library was screened for putative homologs using a PCR-RACE procedure with primers designed to match conserved regions of GOT proteins from a range of species. Five homologs of GOT1 were eventually identified using this screen after cloning and sequencing the products to separate the more closely related paralogs. Initial work was done for the San Diego population in southern California (SD, 32.7457˚N, 117.2550˚W, San Diego County, CA). Three other sites were used to examine the evolution of these GOT1 paralogs, two more in southern California, La Jolla (LJS, 32.8434˚N, 117.2808˚W, San Diego County, CA), and Abalone Cove (AB, 33.7377˚N, 118.3753˚W, Los Angeles County, CA), and one site in central California, Santa Cruz (SCN, 36.9495˚N, 122.0470˚W, Santa Cruz County, CA). These sites were selected because they capture a number of divergent lineages of T. californicus and have been used extensively in other studies of sequence evolution in this species [22, 24, 26].
Primer sequences and amplification conditions for GOT1 paralogs from T. californicus
Forward Primer (5′ to 3′)
Reverse Primer (5′ to 3′)
1751 SD, 1722 AB, 2855 SCN
All sequences were aligned and edited using Sequencer v4.8 software (Genecodes, Ann Arbor, Michigan). The program DNAsp v.5  was used to perform the polymorphism and divergence analyses for each gene. In addition to calculations of polymorphism and divergence (including analyses over sliding windows), Tajima’s D test  was also implemented. The program GENECONV (version 1.81a http://www.math.wustl.edu/~sawyer/geneconv/) was used to identify regions of the paralogous genes that have sequence patterns consistent with gene conversion . Gene conversion events were identified both within and between paralogs within a single population by setting up the group structure within the file and allowing only gene conversion events within populations. The protein variability server (http://imed.med.ucm.es/PVS/) was used to look at patterns of amino acid conservation across GOT1 proteins of arthropods . Conservation was measured by looking at the diversity of amino acids at each site using the Shannon entropy H value.
Phylogenetic trees were constructed using both parsimony and Bayesian analyses with amino acid sequence data and only with parsimony for DNA sequence data from within Tigriopus. The program PAUP*v4b10-×86  was used for the parsimony reconstructions of relationships among GOT1p1/GOT1p2 haplotypes. Heuristic searches were done with 100 random starting trees using either the first 923 bp of the sequence or the last 282 bp in separate analyses. A similar search approach was used for analyses of the divergent sets of GOT amino acid sequences for parsimony analyses. A variety of search conditions using Bayesian analyses and the program MrBayes v3.1.2  were also performed on these protein alignments but did not provide strong support for unresolved relationships in the parsimony analyses.
Availability of supporting data
Sequences are available in Genbank with the accession numbers [KF135593 to KF135616]. The data sets (sequence alignments) supporting the results of this article are available in the Dryad repository http://dx.doi.org/10.5061/dryad.8r6jp.
E. Burch, A. Craven, H. Kunduru collected much of the data for this paper and helped with the analyses. T. Lima provided helpful comments on the manuscript. This work supported by the National Science Foundation (grant DEB-0821003 and IOS-1155325). Two anonymous reviewers provided useful comments as well on a previous draft.
- Innan H: Population genetic models of duplicated genes. Genetica. 2009, 137 (1): 19-37. 10.1007/s10709-009-9355-1.PubMedView ArticleGoogle Scholar
- Chen JM, Cooper DN, Chuzhanova N, Ferec C, Patrinos GP: Gene conversion: mechanisms, evolution and human disease. Nat Rev Genet. 2007, 8 (10): 762-775. 10.1038/nrg2193.PubMedView ArticleGoogle Scholar
- McGrath CL, Casola C, Hahn MW: Minimal effect of ectopic gene conversion among recent duplicates in four mammalian genomes. Genetics. 2009, 182 (2): 615-622. 10.1534/genetics.109.101428.PubMed CentralPubMedView ArticleGoogle Scholar
- Casola C, Ganote CL, Hahn MW: Nonallelic gene conversion in the genus Drosophila. Genetics. 2010, 185 (1): 95-103. 10.1534/genetics.110.115444.PubMed CentralPubMedView ArticleGoogle Scholar
- Mansai SP, Kado T, Innan H: The rate and tract length of gene conversion between duplicated genes. Genes. 2011, 2 (4): 313-331. 10.3390/genes2020313.PubMed CentralPubMedView ArticleGoogle Scholar
- Teshima KM, Innan H: The effect of gene conversion on the divergence between duplicated genes. Genetics. 2004, 166 (3): 1553-1560. 10.1534/genetics.166.3.1553.PubMed CentralPubMedView ArticleGoogle Scholar
- Teshima KM, Innan H: Neofunctionalization of Duplicated Genes Under the Pressure of Gene Conversion. Genetics. 2008, 178 (3): 1385-1398. 10.1534/genetics.107.082933.PubMed CentralPubMedView ArticleGoogle Scholar
- Takuno S, Innan H: Selection to maintain paralogous amino acid differences under the pressure of gene conversion in the heat-shock protein genes in yeast. Mol Biol Evol. 2009, 26 (12): 2655-2659. 10.1093/molbev/msp211.PubMedView ArticleGoogle Scholar
- Mano S, Innan H: The evolutionary rate of duplicated genes under concerted evolution. Genetics. 2008, 180 (1): 493-505. 10.1534/genetics.108.087676.PubMed CentralPubMedView ArticleGoogle Scholar
- Fawcett J, Innan H: Neutral and non-neutral evolution of duplicated genes with gene conversion. Genes. 2011, 2 (4): 191-209. 10.3390/genes2010191.PubMed CentralPubMedView ArticleGoogle Scholar
- Gorski J, Mach B: Polymorphism of human Ia antigens: gene conversion between two DR β loci results in new HLA-D/DR specificity. Nature. 1986, 322 (6074): 67-70. 10.1038/322067a0.PubMedView ArticleGoogle Scholar
- Winderickx J, Battisti L, Hibiya Y, Motulsky AG, Deeb SS: Haplotype diversity in the human red and green opsin genes: Evidence for frequent sequence exchange in exon 3. Hum Mol Genet. 1993, 2 (9): 1413-1421. 10.1093/hmg/2.9.1413.PubMedView ArticleGoogle Scholar
- Parniske M, Hammond-Kosack KE, Golstein C, Thomas CM, Jones DA, Harrison K, Wulff BBH, Jones JDG: Novel disease resistance specificities result from sequence exchange between tandemly repeasted genes at the Cf-4/9 locus of tomato. Cell. 1997, 91 (6): 821-832. 10.1016/S0092-8674(00)80470-5.PubMedView ArticleGoogle Scholar
- Charlesworth D, Bartolome C, Schierup MH, Mable BK: Haplotype structure of the stigmatic self-incompatibility gene in natural populations of Arabidopsis lyrata. Mol Biol Evol. 2003, 20 (11): 1741-1753. 10.1093/molbev/msg170.PubMedView ArticleGoogle Scholar
- Lazzaro BP, Clark AG: Evidence for recurrent paralogous gene conversion and exceptional allelic divergence in the Attacin genes of Drosophila melanogaster. Genetics. 2001, 159 (2): 659-671.PubMed CentralPubMedGoogle Scholar
- Nielsen KM, Kasper J, Choi M, Bedford T, Kristiansen K, Wirth DF, Volkman SK, Lozovsky ER, Hartl DL: Gene conversion as a source of nucleotide diversity in Plasmodium falciparum. Mol Biol Evol. 2003, 20 (5): 726-734. 10.1093/molbev/msg076.PubMedView ArticleGoogle Scholar
- Storz JF, Baze M, Waite JL, Hoffmann FG, Opazo JC, Hayes JP: Complex signatures of selection and gene conversion in the duplicated globin genes of house mice. Genetics. 2007, 177 (1): 481-500. 10.1534/genetics.107.078550.PubMed CentralPubMedView ArticleGoogle Scholar
- Rane HS, Smith JM, Bergthorsson U, Katju V: Gene conversion and DNA sequence polymorphism in the sex-determination gene fog-2 and its paralog ftr-1 in Caenorhabditis elegans. Mol Biol Evol. 2010, 27 (7): 1561-1569. 10.1093/molbev/msq039.PubMed CentralPubMedView ArticleGoogle Scholar
- Burton RS: Intraspecific phylogeography across the Point Conception biogeographic boundary. Evolution. 1998, 52 (3): 734-745. 10.2307/2411268.View ArticleGoogle Scholar
- Edmands S: Phylogeography of the intertidal copepod Tigriopus californicus reveals substantially reduced population differentiation at northern latitudes. Mol Ecol. 2001, 10 (7): 1743-1750. 10.1046/j.0962-1083.2001.01306.x.PubMedView ArticleGoogle Scholar
- Willett CS, Ladner JT: Investigations of fine-scale phylogeography in Tigriopus californicus reveal historical patterns of population divergence. BMC Evol Biol. 2009, 9 (1): 139-10.1186/1471-2148-9-139.PubMed CentralPubMedView ArticleGoogle Scholar
- Willett CS: Quantifying the elevation of mitochondrial DNA evolutionary substitution rates over nuclear rates in the intertidal copepod Tigriopus californicus. J Mol Evol. 2012, 74 (5–6): 310-318.PubMedView ArticleGoogle Scholar
- Barreto FS, Moy GW, Burton RS: Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus. Mol Ecol. 2011, 20 (3): 560-572. 10.1111/j.1365-294X.2010.04963.x.PubMedView ArticleGoogle Scholar
- Willett CS: Complex deleterious interactions associated with malic enzyme may contribute to reproductive isolation in the copepod Tigriopus californicus. PLoS One. 2011, 6 (6): e21177-10.1371/journal.pone.0021177.PubMed CentralPubMedView ArticleGoogle Scholar
- Sawyer SA: GENECONV: A computer package for the statistical detection of gene conversion. 1999, Distributed by the author, available at http://www.math.wustl.edu/~sawyer/geneconv/Google Scholar
- Willett CS: Evolution of interacting proteins in the mitochondrial electron transport system in a marine copepod. Mol Biol Evol. 2004, 21 (3): 443-453.PubMedView ArticleGoogle Scholar
- Thomas PD, Kejariwal A, Campbell MJ, Mi HY, Diemer K, Guo N, Ladunga I, Ulitsky-Lazareva B, Muruganujan A, Rabkin S, et al: PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 2003, 31 (1): 334-341. 10.1093/nar/gkg115.PubMed CentralPubMedView ArticleGoogle Scholar
- Burton RS, Feldman MW: Population genetics of Tigriopus californicus. II. Differentiation among neighboring populations. Evolution. 1981, 35 (6): 1192-1205. 10.2307/2408132.View ArticleGoogle Scholar
- Burton RS: Genetic evidence for long term persistence of marine invertebrate populations in an ephemeral environment. Evolution. 1997, 51 (3): 993-998. 10.2307/2411174.View ArticleGoogle Scholar
- Lichten M, Haber JE: Position effects in ectopic and allelic mitotic recombination in Saccharomyces cerevisiae. Genetics. 1989, 123 (2): 261-268.PubMed CentralPubMedGoogle Scholar
- Dill MM, Burton RS: Genetics of mitochondrial glutamate-oxaloacetate transaminase (GOT-2) in Tigriopus californicus. Biochem Genet. 1984, 22 (3–4): 339-347.PubMedView ArticleGoogle Scholar
- Walsh JB: Sequence-dependent gene conversion: can duplicated genes diverge fast enough to escape conversion. Genetics. 1987, 117 (3): 543-557.PubMed CentralPubMedGoogle Scholar
- Chen WL, Jinks-Robertson S: The role of the mismatch repair machinery in regulating mitotic and meiotic recombination between diverged sequences in yeast. Genetics. 1999, 151 (4): 1299-1313.PubMed CentralPubMedGoogle Scholar
- Godwin AR, Liskay RM: The effects of insertions on mammalian intrachromosomal recombination. Genetics. 1994, 136 (2): 607-617.PubMed CentralPubMedGoogle Scholar
- Lukacsovich T, Waldman AS: Suppression of intrachromosomal gene conversion in mammalian cells by small degrees of sequence divergence. Genetics. 1999, 151 (4): 1559-1568.PubMed CentralPubMedGoogle Scholar
- Innan H: A two-locus gene conversion model with selection and its application to the human RHCE and RHD genes. P Natl Acad Sci USA. 2003, 100 (15): 8793-8798. 10.1073/pnas.1031592100.View ArticleGoogle Scholar
- Osada N, Innan H: Duplication and gene conversion in the Drosophila melanogaster genome. PLoS Genet. 2008, 4 (12): e1000305-10.1371/journal.pgen.1000305.PubMed CentralPubMedView ArticleGoogle Scholar
- Bosch E, Hurles ME, Navarro A, Jobling MA: Dynamics of a human interparalog gene conversion hotspot. Genome Res. 2004, 14 (5): 835-844. 10.1101/gr.2177404.PubMed CentralPubMedView ArticleGoogle Scholar
- Ohta T: Role of Diversifying Selection and Gene Conversion in Evolution of Major Histocompatibility Complex Loci. P Natl Acad Sci USA. 1991, 88 (15): 6716-6720. 10.1073/pnas.88.15.6716.View ArticleGoogle Scholar
- Hoelzel AR, Green A: Analysis of population-level variation by sequencing PCR-amplified DNA. Practical Approach Series: Molecular Genetic Analysis Of Populations. Edited by: Hoelzel AR. 1992, New York, New York: Oxford University Press, 159-187.Google Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.PubMedView ArticleGoogle Scholar
- Tajima F: Statistical methods for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Garcia-Boronat M, Diez-Rivero CM, Reinherz EL, Reche PA: PVS: a web server for protein sequence variability analysis tuned to facilitate conserved epitope discovery. Nucleic Acids Res. 2008, 36: W35-W41. 10.1093/nar/gkn211.PubMed CentralPubMedView ArticleGoogle Scholar
- Swofford DL: PAUP*. phylogenetic analysis using parsimony (*and other methods). In., 4 edn. 1993, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.