horizontal gene transfer ( hgt ) is very common and of great importance in bacterial evolution and is also relatively common in certain eukaryotic lineages . one consequence of hgt is incongruence between phylogenies reconstructed from genes with different histories of transfer , including of course no transfer at all . finding and examining phylogenetically atypical genes has become a standard task in hgt studies . in cases where an entire gene has been transferred , without the complication of recombination / conversion with a native homolog , the donor lineage can in principle gene conversion enters the picture either during the act of gene transfer , when transiently - present foreign dna directly converts ( replaces ) part of a native locusor , after duplicative hgt , via potentially ongoing gene conversion between co - existing native and foreign copies . overall , then , gene conversion can lead to a potentially complex and diverse set of patchwork recombinant sequences , especially if it occurs repeatedly , and differentially , over the course of speciation . each recombinant gene , if analyzed as a whole , might or might not reflect the true evolutionary history of either or both parental sequences . when parental sequences contribute differentially to the number of informative characters in a recombinant sequence , this sequence will tend to resemble the parental sequence that contributes more informative characters ( e.g. , refs . 15 and 16 ) , whereas when parental sequences contribute similar numbers of informative characters , the recombinant will potentially be quite different from both parental sequences , depending of course on the degree of divergence of the two parental sequences from each other . when properly recognized and dealt with , recombination poses few problems for phylogenetic analysis and interpretation . in practice , however , recombination detection is challenging and often subject to failure . first , it is well established that recombination detection programs perform poorly when sequence divergence is low . unfortunately , plant mitochondrial genomes , which collectively constitute a premiere model system for eukaryotic hgt studies , usually have very low rates of nucleotide substitution . second , gene conversion often involves very short tracts of dna , which can make recombination detection very difficult . for instance , ten previously published recombinant regions between plant mitochondrial and chloroplast genes range in length from only 14 to 79 nucleotides . third , existing recombination detection programs are generally designed to identify a single or only a small number of recombination breakpoints . intricate gene conversion during the process of duplicative hgt and differential gene conversion ( dh - dc ) can , however , lead to mosaic gene structures , with multiple foreign regions interspersed with native regions on a fine scale . we recently reported such mosaicism in mitochondrial atp1 and matr genes belonging to different groups of flowering plants ( angiosperms ) . we show that these mosaic genes largely escaped detection by recombination - detection programs and were recognizable only by direct visual inspection of dna sequence alignments . in this report , we explore the effects of chimeric sequences on phylogeny reconstruction by conducting phylogenetic analyses on simulated sequences and on an augmented set of the atp1 sequences analyzed in reference 12 . in a recent study , we reported the presence of three differentially mosaic types of mitochondrial atp1 genes in the angiosperm genus ternstroemia ( pentaphylaceae , ericales ) and concluded that they arose via dh - dc , with the blueberry genus vaccinium ( ericaceae ) the best candidate to be the donor group in the initiating hgt event . 12 ) , each of the three major clades within ternstroemia possesses a differentially mosaic atp1 gene , each with multiple ( four to five ) foreign regions interspersed with native regions . in the one atp1 phylogeny presented in reference 12 , the mosaic ternstroemia genes were all placed within the eriaceae , in a paraphyletic relationship with respect to vaccinium . black vertical lines represent the 38 nucleotide positions inferred12 to have differed between donor and recipient atp1 genes at the time of atp1 transfer from vaccinium to a common ancestor of ternstroemia . lines at the top of the boxes and red shading indicate sites and regions , respectively , of putatively foreign , vaccinium ancestry , while bottom lines and blue shading represent native sites and regions . white lines centered within the boxes represent the only two sites that otherwise differ within the ternstroemia clade . t - ips  refers to the ternstroemia subclade containing t. impressa , t. peduncularis and t. stahlii ) and  t - glj  to the subclade containing t. gymnanthera , t. longipes and t. japonica ( see also fig . ( a ) chronogram showing organismal relationships and divergence times of relevant taxa belonging to the ericales . as described in reference 12 , the chronogram was constructed by the beast program using a eurya reference fossil calibration of 86 myr ago . ( b - f ) maximum likelihood phylogenies of mitochondrial atp1 genes from the taxa shown in ( a ) , with these analyses varying as to which members of ternstroemia ( shown in red ) , whose atp1 gene is differentially mosaic , were included . a total of 1000 bootstrap iterations were performed , with all bootstrap values  50% shown on the trees . phylogenies were rooted using fouquieria , marcgravia and pentamerista as unshown outgroups ( hence the stub branch at the base of each gene tree ) . to better understand how mosaic genes affect phylogeny reconstruction , we sequenced the mitochondrial atp1 gene of chamaedaphne calyculata , a close relative of vaccinium , using the same set of primers and method as in reference 12 , and employed it together with varying subsets of previously sequenced ericales atp1 genes in maximum likelihood phylogenetic analyses . consistent with our recent study , and in contrast to organismal phylogeny ( fig . 2a ) , in an analysis that included all relevant genes , the three types of mosaic atp1 genes in ternstroemia formed a paraphyletic assemblage , with t. fragrans sister to the vaccinium / chamaedaphne clade , the t - glj clade the most distant from the vaccinium / chamaedaphne clade and the t - ips clade in an intermediate position ( fig . remarkably , when only one type of mosaic atp1 gene was included in a given analysis , each of the three types fell in a different phylogenetic position ( fig .  2c  e ) . this shows that the phylogenetic position of a mosaic gene can vary depending on the inclusion of additional , related mosaic genes and emphasizes that this position provides little or no reliable information on the nature of the gene s parental sequences . 2b  f ) involves chamaedaphne , which was weakly placed ( 41% bootstrap support ) within vaccinium when all mosaic genes were included ( fig . 2b ) , but was placed as sister ( with 6892% support ) to a monophyletic vaccinium in the other four gene trees . this result raises the possibility that the inclusion of mosaic sequences in phylogenetic analyses can affect not only the placement of related mosaic sequences , but also the placement of apparently native sequences . we used simulation studies ( conducted using seq - gen ) to further explore the effects of chimeric sequences on phylogeny reconstruction . the following simulation parameters were chosen to be the same as those used in the above analysis of atp1 sequences : ( 1 ) sequence length , 1200 nucleotides ; ( 2 ) substitution model , gtr ; ( 3 ) gamma shape parameter , 0.218 ; ( 4 ) proportion of invariant sites : 0.371 ; ( 5 ) nucleotide frequencies , 0.271 ( a ) , 0.207 ( c ) , 0.261 ( g ) and 0.261 ( t ) ; and ( 6 ) gtr relative rate parameters , a < - > c = 0.818 , a < - > g = 1.938 , a < - > t = 0.244 , c < - > g = 0.884 , c all but the first of these parameters are based on phyml estimates derived from the mitochondrial atp1 alignment shown in figure  s1 , except with the mosaic ternstroemia genes excluded ( as in fig . 2f ) because recombinant genes have been shown to alter the estimation of substitution rate heterogeneity . the remaining simulation parameters were independent of the atp1 data and include the number of sequences , their topology , relative branch lengths and absolute amount of divergence ( fig . 3a ) . for simplicity , only one recombination breakpoint was allowed in each chimeric sequence , with these chimerics constructed to contain varying proportions ( from 50:50 to 10:90 ) of two of the 16  parental  sequences generated by the simulations . the 16 purely - simulated sequences together with the artificially constructed chimeric sequence were used in phylogenetic analyses performed using raxml version 7.0.4 . chimeric sequences were generated by combining the 5 and 3 portions of sequences 1 and 9 , respectively ( both circled in red ) , with different length ratios ( 10:90 , 30:70 or 50:50 ) of the two parental sequences . is based ( for computational ease ) on the concatenated sequences of the first 100 iterations , while the bootstrap support values are from all 1000 trees . two sets of simulation analyses were performed . in the first set ( fig . 3 ) , we varied the parental proportions that comprise the chimeric sequence(s ) and the number of chimeric sequences included in a given analysis . in analyses with a single chimeric sequence , this sequence grouped with 100% bootstrap support with its majority parental sequence when the parental ratio was 10:90 ( i.e. , when the chimeric sequence consisted of 10% of sequence 9 and 90% of sequence 1 ; fig . 3d ) , it went to the base of the tree , in between the two main clades of simulated sequences , and with bootstrap support reduced along the branches leading to both parental sequences . these results thus show that the phylogenetic position of chimeric genes can vary substantially as the proportion of parental sequences that comprise the chimerics varies . the phylogenetic position of the 50:50 chimeric sequence also varied substantially depending on whether it was the only chimeric sequence in the analysis ( fig . 3d ) or whether the 30:70 and/or 10:90 sequences were also included ( fig .  3e  g ) . the 50:50 from the base of the tree to its periphery , together with the 30:70 and parental sequence 1 ( and with the 10:90 when included ) and with strong support ( 96% and 92% , respectively ) . this peripheral attraction was more subdued when the 50:50 was paired with the 10:90 ( fig . 3f ) , presumably because of the greater proportion of sequence length shared by the 30:70 and 50:50 ( 80% ) relative to the 10:90 and 50:50 ( 60% ) . also , there is evidence for a mutual attraction between the 10:90 and 50:50 ( but not between the 30:70 and 50:50 ) , in that the bootstrap support for the 10:90 grouping with parental sequence 1 was reduced from 100% ( fig . 2 ) in showing that the phylogenetic position of a chimeric gene can vary substantially when different additional chimeric genes are sampled . we also found that chimeric genes can perturb phylogenetic analysis by producing branch length bias . first , certain chimeric genes were directly associated with elongated branch lengths . in the single - chimeric analyses , the 10:90 chimera had a notably longer branch length than sequence 1 ( fig . 3c ) , with the ( homoplasious ) substitutions contributed by the 30% of the chimeric sequence originating from sequence 9 presumably responsible for this greatly extended branch length . similar results were obtained in figure  3e  h , where multiple chimeric sequences were included in each analysis . second , the presence of the 50:50 chimeric in the single - chimeric analysis of figure  3d resulted in the loss of a molecular clock among the non - recombinant sequences in this analysis , with the branch leading to the ( 14 ) clade only 60% of the length leading to the ( 58 ) clade and the branch leading to the ( 912 ) clade only 58% of the length leading to the ( 1316 ) clade . there are consistently two additional  albeit much less pronounced  sets of branch length differences in the six trees ( fig . 3b , c and e  h ) in which one or more chimeric sequences are sister to sequence 1 . the second set of simulations showed that chimeric sequences can also confound phylogenetic analysis by altering the placement of non - recombinant ( native ) sequences . in these simulations ( fig . 4 ) , we varied the absolute amount of sequence divergence across the tree , with the 50:50 sequence the only chimeric sequence in each analysis . figure  4a shows the same tree as figure  3d , while figure  4b  d has expanded divergence by a factor of 2  , 5  and 10  , respectively . as the simulated sequences become more divergent , the parental sequences show an increased and pronounced , tendency to be attracted toward the base of the tree by the chimeric gene . also , the branch leading to the chimeric sequence becomes increasingly short , approaching zero in figure  4candd . it is important to note that this effect is not simply the result of there being increasingly more informative characters from figure  4a  d . in analyses with 5  or 10  divergence ( data not shown ) , but only 1/10th the sequence length ( 120 nucleotides rather than the 1200 in the simulations shown in figure  4 ) , the topology and bootstrap values were essentially identical with those in figure  4candd , respectively , with these being substantially different from those in figure  4aandb . the different topologies shown in figure  4 must therefore result from deterministic effects arising from the varying levels of divergence in these trees . these simulations thus show that inclusion of chimeric sequences can distort the branching pattern of nonchimeric , native sequences and that the extent of this distortion varies directly with the absolute level of sequence divergence . ( b  d ) maximum likelihood analysis of the same simulated sequences , but with increasing branch - length scales across the set of sequences : ( b ) 2  the scale in ( a ) ; ( c ) 5  ; and ( d ) 10  . is based ( for computational ease ) on the concatenated sequences of the first 100 iterations , while the bootstrap support values are from all 1,000 trees . in a recent study , we reported the presence of three differentially mosaic types of mitochondrial atp1 genes in the angiosperm genus ternstroemia ( pentaphylaceae , ericales ) and concluded that they arose via dh - dc , with the blueberry genus vaccinium ( ericaceae ) the best candidate to be the donor group in the initiating hgt event . 12 ) , each of the three major clades within ternstroemia possesses a differentially mosaic atp1 gene , each with multiple ( four to five ) foreign regions interspersed with native regions . in the one atp1 phylogeny presented in reference 12 , the mosaic ternstroemia genes were all placed within the eriaceae , in a paraphyletic relationship with respect to vaccinium . black vertical lines represent the 38 nucleotide positions inferred12 to have differed between donor and recipient atp1 genes at the time of atp1 transfer from vaccinium to a common ancestor of ternstroemia . lines at the top of the boxes and red shading indicate sites and regions , respectively , of putatively foreign , vaccinium ancestry , while bottom lines and blue shading represent native sites and regions . white lines centered within the boxes represent the only two sites that otherwise differ within the ternstroemia clade . t - ips  refers to the ternstroemia subclade containing t. impressa , t. peduncularis and t. stahlii ) and  t - glj  to the subclade containing t. gymnanthera , t. longipes and t. japonica ( see also fig . ( a ) chronogram showing organismal relationships and divergence times of relevant taxa belonging to the ericales . as described in reference 12 , the chronogram was constructed by the beast program using a eurya reference fossil calibration of 86 myr ago . ( b - f ) maximum likelihood phylogenies of mitochondrial atp1 genes from the taxa shown in ( a ) , with these analyses varying as to which members of ternstroemia ( shown in red ) , whose atp1 gene is differentially mosaic , were included . a total of 1000 bootstrap iterations were performed , with all bootstrap values  50% shown on the trees . phylogenies were rooted using fouquieria , marcgravia and pentamerista as unshown outgroups ( hence the stub branch at the base of each gene tree ) . to better understand how mosaic genes affect phylogeny reconstruction , we sequenced the mitochondrial atp1 gene of chamaedaphne calyculata , a close relative of vaccinium , using the same set of primers and method as in reference 12 , and employed it together with varying subsets of previously sequenced ericales atp1 genes in maximum likelihood phylogenetic analyses . consistent with our recent study , and in contrast to organismal phylogeny ( fig . 2a ) , in an analysis that included all relevant genes , the three types of mosaic atp1 genes in ternstroemia formed a paraphyletic assemblage , with t. fragrans sister to the vaccinium / chamaedaphne clade , the t - glj clade the most distant from the vaccinium / chamaedaphne clade and the t - ips clade in an intermediate position ( fig . remarkably , when only one type of mosaic atp1 gene was included in a given analysis , each of the three types fell in a different phylogenetic position ( fig .  2c  e ) . this shows that the phylogenetic position of a mosaic gene can vary depending on the inclusion of additional , related mosaic genes and emphasizes that this position provides little or no reliable information on the nature of the gene s parental sequences . 2b  f ) involves chamaedaphne , which was weakly placed ( 41% bootstrap support ) within vaccinium when all mosaic genes were included ( fig . 2b ) , but was placed as sister ( with 6892% support ) to a monophyletic vaccinium in the other four gene trees . this result raises the possibility that the inclusion of mosaic sequences in phylogenetic analyses can affect not only the placement of related mosaic sequences , but also the placement of apparently native sequences . we used simulation studies ( conducted using seq - gen ) to further explore the effects of chimeric sequences on phylogeny reconstruction . the following simulation parameters were chosen to be the same as those used in the above analysis of atp1 sequences : ( 1 ) sequence length , 1200 nucleotides ; ( 2 ) substitution model , gtr ; ( 3 ) gamma shape parameter , 0.218 ; ( 4 ) proportion of invariant sites : 0.371 ; ( 5 ) nucleotide frequencies , 0.271 ( a ) , 0.207 ( c ) , 0.261 ( g ) and 0.261 ( t ) ; and ( 6 ) gtr relative rate parameters , a < - > c = 0.818 , a < - > g = 1.938 , a < - > t = 0.244 , c < - > g = 0.884 , c < - > t = 2.219 and g all but the first of these parameters are based on phyml estimates derived from the mitochondrial atp1 alignment shown in figure  s1 , except with the mosaic ternstroemia genes excluded ( as in fig . 2f ) because recombinant genes have been shown to alter the estimation of substitution rate heterogeneity . the remaining simulation parameters were independent of the atp1 data and include the number of sequences , their topology , relative branch lengths and absolute amount of divergence ( fig . for simplicity , only one recombination breakpoint was allowed in each chimeric sequence , with these chimerics constructed to contain varying proportions ( from 50:50 to 10:90 ) of two of the 16  parental  sequences generated by the simulations . the 16 purely - simulated sequences together with the artificially constructed chimeric sequence were used in phylogenetic analyses performed using raxml version 7.0.4 . chimeric sequences were generated by combining the 5 and 3 portions of sequences 1 and 9 , respectively ( both circled in red ) , with different length ratios ( 10:90 , 30:70 or 50:50 ) of the two parental sequences . is based ( for computational ease ) on the concatenated sequences of the first 100 iterations , while the bootstrap support values are from all 1000 trees . two sets of simulation analyses were performed . in the first set ( fig . 3 ) , we varied the parental proportions that comprise the chimeric sequence(s ) and the number of chimeric sequences included in a given analysis . in analyses with a single chimeric sequence , this sequence grouped with 100% bootstrap support with its majority parental sequence when the parental ratio was 10:90 ( i.e. , when the chimeric sequence consisted of 10% of sequence 9 and 90% of sequence 1 ; fig . 3d ) , it went to the base of the tree , in between the two main clades of simulated sequences , and with bootstrap support reduced along the branches leading to both parental sequences . these results thus show that the phylogenetic position of chimeric genes can vary substantially as the proportion of parental sequences that comprise the chimerics varies . the phylogenetic position of the 50:50 chimeric sequence also varied substantially depending on whether it was the only chimeric sequence in the analysis ( fig . 3d ) or whether the 30:70 and/or 10:90 sequences were also included ( fig . the 50:50 from the base of the tree to its periphery , together with the 30:70 and parental sequence 1 ( and with the 10:90 when included ) and with strong support ( 96% and 92% , respectively ) . this peripheral attraction was more subdued when the 50:50 was paired with the 10:90 ( fig . 3f ) , presumably because of the greater proportion of sequence length shared by the 30:70 and 50:50 ( 80% ) relative to the 10:90 and 50:50 ( 60% ) . also , there is evidence for a mutual attraction between the 10:90 and 50:50 ( but not between the 30:70 and 50:50 ) , in that the bootstrap support for the 10:90 grouping with parental sequence 1 was reduced from 100% ( fig . 2 ) in showing that the phylogenetic position of a chimeric gene can vary substantially when different additional chimeric genes are sampled . we also found that chimeric genes can perturb phylogenetic analysis by producing branch length bias . first , certain chimeric genes were directly associated with elongated branch lengths . in the single - chimeric analyses , the 10:90 chimera had a notably longer branch length than sequence 1 ( fig . 3c ) , with the ( homoplasious ) substitutions contributed by the 30% of the chimeric sequence originating from sequence 9 presumably responsible for this greatly extended branch length . similar results were obtained in figure  3e  h , where multiple chimeric sequences were included in each analysis . second , the presence of the 50:50 chimeric in the single - chimeric analysis of figure  3d resulted in the loss of a molecular clock among the non - recombinant sequences in this analysis , with the branch leading to the ( 14 ) clade only 60% of the length leading to the ( 58 ) clade and the branch leading to the ( 912 ) clade only 58% of the length leading to the ( 1316 ) clade . there are consistently two additional  albeit much less pronounced  sets of branch length differences in the six trees ( fig . 3b , c and e  h ) in which one or more chimeric sequences are sister to sequence 1 . the second set of simulations showed that chimeric sequences can also confound phylogenetic analysis by altering the placement of non - recombinant ( native ) sequences . in these simulations ( fig . 4 ) , we varied the absolute amount of sequence divergence across the tree , with the 50:50 sequence the only chimeric sequence in each analysis . figure  4a shows the same tree as figure  3d , while figure  4b  d has expanded divergence by a factor of 2  , 5  and 10  , respectively . as the simulated sequences become more divergent , the parental sequences show an increased and pronounced , tendency to be attracted toward the base of the tree by the chimeric gene . also , the branch leading to the chimeric sequence becomes increasingly short , approaching zero in figure  4candd . it is important to note that this effect is not simply the result of there being increasingly more informative characters from figure  4a  d . in analyses with 5  or 10  divergence ( data not shown ) , but only 1/10th the sequence length ( 120 nucleotides rather than the 1200 in the simulations shown in figure  4 ) , the topology and bootstrap values were essentially identical with those in figure  4candd , respectively , with these being substantially different from those in figure  4aandb . the different topologies shown in figure  4 must therefore result from deterministic effects arising from the varying levels of divergence in these trees . these simulations thus show that inclusion of chimeric sequences can distort the branching pattern of nonchimeric , native sequences and that the extent of this distortion varies directly with the absolute level of sequence divergence . ( b  d ) maximum likelihood analysis of the same simulated sequences , but with increasing branch - length scales across the set of sequences : ( b ) 2  the scale in ( a ) ; ( c ) 5  ; and ( d ) 10  . is based ( for computational ease ) on the concatenated sequences of the first 100 iterations , while the bootstrap support values are from all 1,000 trees . the analyses reported in this study show that the inclusion in phylogenetic analyses of chimeric sequences arising from hgt and gene conversion can produce a variety of spurious phylogenetic results . these include the misplacement of both chimeric and native sequences , as well as branch length distortions . accordingly , we introduce the term  hgt turbulence  as a general moniker for this category of phylogenetic artifacts . mutually reinforcing evidence for these types of hgt turbulence was apparent in the phylogenetic analyses of both simulated sequences and a naturally occurring set of native and chimeric mitochondrial sequences , with the simulations providing greater it is also important to realize that while these simulations were framed and presented in the context of hgt , their results apply equally well to conversion between paralogs arising from internal gene duplication as to xenologs arising from hgt . hgt turbulence is probably relatively common in bacteria , given the prevalence of hgt and recombination during bacterial genome evolution . for example , such diverse bacteria as neisseria meningitidis , 
 streptococcus pneumoniae , 
 helicobacter pylori and wolbachia have been found to be so recombinogenic that scientists have resorted to using multiple loci ( e.g. , multilocus sequence typing ) as opposed to a single locus to identify clones . surprisingly , however , the phenomenon of hgt turbulence has never been explicitly addressed in the bacterial literature . one reason for this is that these studies have mainly focused on minimizing the effect of recombination and thereby inferring accurate evolutionary relationships , or on quantifying the number of recombinant genes , rather than actually exploring the topological alterations caused by hgt turbulence . also , the precise origins of horizontally transferred genes ( or gene fragments ) in bacterial genomes can be extremely difficult to recover , especially when transfer and/or subsequent recombination have occurred on a fine scale . several other studies have used recombinant sequences in phylogenetic simulations , but most of these have focused on the issue of whether recombinants are detectable and/or how to detect them . perhaps most relevant to our study is the 2002 study by posada and crandall , which also explored the relationship between the location of recombination breakpoints and the phylogenetic placement of recombinant sequences , reaching similar conclusion to ours on this point . however , none of these studies explored the effects of combining multiple , related chimeric sequences in the same analysis , nor did they show that chimeric sequences can , under certain conditions , substantially alter the phylogenetic behavior of native sequences . further simulation studies on the effects of recombination on phylogenetic inference would therefore appear to be called for . recognition of the remarkable frequency and extent of horizontal gene transfer , and its often great evolutionary importance , is arguably the greatest accomplishment of the past 15 years of comparative genomics research . because hgt is so common and important , recognizing and properly dealing with hgt turbulence is likewise important . this is so not only because of the obvious need for obtaining accurate estimates of gene and species phylogeny , but also because otherwise too many cases of chimeric hgt , including complex situations involving dh - dc , will continue to go overlooked .
horizontal gene transfer ( hgt ) often leads to phylogenetic incongruence . when  duplicative hgt 
 introduces a second copy of a pre - existing gene , the two copies may then engage in gene conversion , leading to phylogenetically mosiac genes . 
 when duplicative hgt is followed by differential gene conversion among descendant lineages , as under the dh - dc model , phylogenetic analysis is further complicated . to explore the effects of dh - dc on phylogeny reconstruction , we analyzed two sets of sequences : ( 1 ) an augmented set of plant mitochondrial atp1 sequences for which we recently published evidence of dh - dc ; and ( 2 ) a set of simulated sequences for which we varied the extent of chimerism , the number of chimeric genes and nucleotide substitution rates . 
 we show that the phylogenetic behavior of evolutionarily chimeric genes is highly volatile and depends on both the degree of chimerism and the number of differentially chimeric genes present in the analysis . 
 furthermore , we show that the presence of chimeric genes in gene trees can spuriously affect the phylogenetic position of purely native sequences , especially by attracting these sequences toward basal positions in trees . 
 we propose the term  hgt turbulence  to describe these complex effects of evolutionarily chimeric genes on phylogenetic results .