Codon Usage Methods for Horizontal Gene Transfer Detection Generate an Abundance of False Positive and False Negative Results

Bacteria acquire new DNA in a process known as horizontal gene transfer (HGT). To investigate the evolutionary impact of this transfer of DNA, various methods have been developed to detect past HGT events. For example, codon usage-based methods detect the presence of transferred genes by identifying atypical patterns of codon usage. However, some inherited genes exhibit atypical codon usage and some transferred genes have codon usage patterns similar to those of the inherited genes. In this study, we used a comparative phylogenetic approach with Methylobacterium and Caulobacter species to demonstrate that even well-designed codon usage methods fail to detect many HGT events and generate a high rate of false positives (60–75 %) and false negatives (23–61 %). Therefore, we recommend caution when employing codon usage methods to identify transferred genes and suggest that the rapidly increasing availability of bacterial genome sequences makes the phylogenetic approach the method of choice.     

A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5′ end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of ∼250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

Positive Darwinian Selection Promotes Heterogeneity Among Members of the Antifreeze Protein Multigene Family

A variety of organisms have independently evolved proteins exhibiting antifreeze activity that allows survival at subfreezing temperatures. The antifreeze proteins (AFPs) bind ice nuclei and depress the freezing point by a noncolligative absorption–inhibition mechanism. Many organisms have a heterogeneous suite of AFPs with variation in primary sequence between paralogous loci. Here, we demonstrate that the diversification of the AFP paralogues is promoted by positive Darwinian selection in two independently evolved AFPs from fish and beetle. First, we demonstrate an elevated rate of nonsynonymous substitutions compared to synonymous substitutions in the mature protein coding region. Second, we perform phylogeny-based tests of selection to demonstrate a subset of codons is subjected to positive selection. When mapped onto the three-dimensional structure of the fish antifreeze type III antifreeze structure, these codons correspond to amino acid positions that surround but do not interrupt the putative ice-binding surface. The selective agent may be related to efficient binding to diverse ice surfaces or some other aspect of AFP function. Received: 27 February 2001 / Accepted: 12 September 2001     

Assessment of Anti-recombination and Double-strand Break-induced Gene Conversion in Human Cells by a Chromosomal Reporter

Gene conversion is one of the frequent end results of homologous recombination, and it often underlies the inactivation of tumor suppressor genes in cancer cells. Here, we have developed an integrated assay system that allows simultaneous examination of double-strand break (DSB)-induced gene conversion events at the site of a DSB (proximal region) and at a surrounding region ～1 kb away from the break (distal region). Utilizing this assay system, we find that gene conversion events at the proximal and distal regions are relatively independent of one another. The results also indicate that synthesis-dependent strand annealing (SDSA) plays a major role in DSB-induced gene conversion. In addition, our current study has demonstrated that hMLH1 plays an essential role in anti-recombination and gene conversion. Specifically, the anti-recombination activity of hMLH1 is partially dependent on its interaction with hMRE11. Our data suggests that the role of hMLH1 and hMRE11 in the process of gene conversion is complex, and these proteins play different roles in DSB-induced proximal and distal gene conversions. In particular, the involvement of hMLH1 and hMRE11 in the distal gene conversion requires both hMSH2 and heteroduplex formation.     

Positive Selection on HIV Accessory Proteins and the Analysis of Molecular Adaptation After Interspecies Transmission

Studies examining positive selection on accessory proteins of HIV are rare, although these proteins play an important role in pathogenesis in vivo. Moreover, despite the biological relevance of analyses of molecular adaptation after viral transmission between species, the issue is still poorly studied. Here we present evidence that accessory proteins are subjected to positive selective forces exclusively in HIV. This scenario suggests that accessory protein genes are under adaptive evolution in HIV clades, while in SIVcpz such a phenomenon could not be detected. As a result, we show that comparative studies are critical to carry out functional investigation of positively selected protein sites, as they might help to achieve a better comprehension of the biology of HIV pathogenesis.     

酵母PH081在酸性磷酸酯酶基因表达调控中的作用

利用酵母ＰＨＯ８１与大肠杆菌产β－半乳糖苷酶的融合基因,系统地研究了ＰＨＯ８１基因在酵母酸性磷酸醋酶的不同调控因子的单缺失株和双缺失株中的表达规律,它与酸性磷酸酯酶基因ＰＨＯ５和ＰＨＯ１１的控制机制相似。构建了ＡＤＨ１启动子控制下ＰＨＯ８１表达质粒。在ＰＨＯ８１在细胞内组成型表达时,酸性磷酸酯酶基因的表达仍受无机磷的控制,揭示ＰＨＯ８１蛋白可能在不同浓度无机磷条件发生变构。ＰＨＯ８１在酸性磷酸酯酶基因表达系统中是一个中介因子。在此基础提出了一个酸性磷酸酯酶基因调控的分子模型。     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Molecular Mechanisms for the Variation of Mitochondrial Gene Content and Gene Arrangement Among Chigger Mites of the Genus Leptotrombidium (Acari: Acariformes)

The gene content of a mitochondrial (mt) genome, i.e., 37 genes and a large noncoding region (LNR), is usually conserved in Metazoa. The arrangement of these genes and the LNR is generally conserved at low taxonomic levels but varies substantially at high levels. We report here a variation in mt gene content and gene arrangement among chigger mites of the genus Leptotrombidium. We found previously that the mt genome of Leptotrombidium pallidum has an extra gene for large-subunit rRNA (rrnL), a pseudo-gene for small-subunit rRNA (PrrnS), and three extra LNRs, additional to the 37 genes and an LNR typical of Metazoa. Further, the arrangement of mt genes of L. pallidum differs drastically from that of the hypothetical ancestor of the arthropods. To find to what extent the novel gene content and gene arrangement occurred in Leptotrombidium, we sequenced the entire or partial mt genomes of three other species, L. akamushi, L. deliense, and L. fletcheri. These three species share the arrangement of all genes with L. pallidum, except trnQ (for tRNA-glutamine). Unlike L. pallidum, however, these three species do not have extra rrnL or PrrnS and have only one extra LNR. By comparison between Leptotrombidium species and the ancestor of the arthropods, we propose that (1) the type of mt genome present in L. pallidum evolved from the type present in the other three Leptotrombidium species, and (2) three molecular mechanisms were involved in the evolution of mt gene content and gene arrangement in Leptotrombidium species. [Reviewing Editor: Dr. Martin Kreitman]     

Selection Conflicts,Gene Expression,and Codon Usage Trends in Yeast

Synonymous codon usage in yeast appears to be influenced by natural selection on gene expression, as well as regional variation in compositional bias. Because of the large number of potential targets of selection (i.e., most of the codons in the genome) and presumed small selection coefficients, codon usage is an excellent model for studying factors that limit the effectiveness of selection. We use factor analysis to identify major trends in codon usage for 5836 genes in Saccharomyces cerevisiae. The primary factor is strongly correlated with gene expression, consistent with the model that a subset of codons allows for more efficient translation. The secondary factor is very strongly correlated with third codon position GC content and probably reflects regional variation in compositional bias. We find that preferred codon usage decreases in the face of three potential limitations on the effectiveness of selection: reduced recombination rate, increased gene length, and reduced intergenic spacing. All three patterns are consistent with the Hill–Robertson effect (reduced effectiveness of selection among linked targets). A reduction in gene expression in closely spaced genes may also reflect selection conflicts due to antagonistic pleiotropy.     

Widespread Gene Conversion of Alpha-2-Fucosyltransferase Genes in Mammals

The alpha-2-fucosyltransferases (α2FTs) are enzymes involved in the biosynthesis of α2fucosylated glycan structures. In mammalian genomes, there are three α2FT genes located in tandem—FUT1, FUT2, and Sec1—each contained within a single exon. It has been suggested that these genes originated from two successive duplications, with FUT1 being generated first and FUT2 and Sec1 second. Despite gene conversion being considered the main mechanism of concerted evolution in gene families, previous studies of primates α2FTs failed to detect it, although the occurrence of gene conversion between FUT2 and Sec1 was recently reported in a human allele. The primary aim of our work was to initiate a broader study on the molecular evolution of mammalian α2FTs. Sequence comparison leads us to confirm that the three genes appeared by two rounds of duplication. In addition, we were able to detect multiple gene-conversion events at the base of primates and within several nonprimate species involving FUT2 and Sec1. Gene conversion involving FUT1 and either FUT2 or Sec1 was also detected in rabbit. The extent of gene conversion between the α2FTs genes appears to be species-specific, possibly related to functional differentiation of these genes. With the exception of rabbits, gene conversion was not observed in the region coding the C-terminal part of the catalytic domain. In this region, the number of amino acids that are identical between FUT1 and FUT2, but different in Sec1, is higher than in other parts of the protein. The biologic meaning of this observation may be related to functional constraints. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Roles of Gene Duplication,Gene Conversion and Positive Selection in Rodent Esp and Mup Pheromone Gene Families with Comparison to the Abp Family

Three proteinaceous pheromone families, the androgen-binding proteins (ABPs), the exocrine-gland secreting peptides (ESPs) and the major urinary proteins (MUPs) are encoded by large gene families in the genomes of Mus musculus and Rattus norvegicus. We studied the evolutionary histories of the Mup and Esp genes and compared them with what is known about the Abp genes. Apparently gene conversion has played little if any role in the expansion of the mouse Class A and Class B Mup genes and pseudogenes, and the rat Mups. By contrast, we found evidence of extensive gene conversion in many Esp genes although not in all of them. Our studies of selection identified at least two amino acid sites in β-sheets as having evolved under positive selection in the mouse Class A and Class B MUPs and in rat MUPs. We show that selection may have acted on the ESPs by determining K_a/K_s for Exon 3 sequences with and without the converted sequence segment. While it appears that purifying selection acted on the ESP signal peptides, the secreted portions of the ESPs probably have undergone much more rapid evolution. When the inner gene converted fragment sequences were removed, eleven Esp paralogs were present in two or more pairs with K_a/K_s >1.0 and thus we propose that positive selection is detectable by this means in at least some mouse Esp paralogs. We compare and contrast the evolutionary histories of all three mouse pheromone gene families in light of their proposed functions in mouse communication.     

Gene Conversion and Positive Selection Driving the Evolution of the <Emphasis Type="Italic">Caenorhabditis</Emphasis> ssp. ZIM/HIM-8 Protein Family

In C. elegans, four C2H2 zinc-finger proteins (ZIM-1, ZIM-2, ZIM-3, and HIM-8), which are arranged in tandem, mediate chromosome-specific pairing and synapsis during meiosis. The zim/him-8 genes from three Caenorhabditis species were contrasted in an effort to investigate the mechanisms driving their evolution. Here it is shown that the preservation of higher degree of sequence similarity in the N-terminal portion, particularly in several regions within the second exon between paralogous zim genes (especially between zim-1 and zim-3), is due to independent interparalogue gene conversions. However, the evolutionary force is not uniformly strong across species. The present data reveal that more frequent gene conversion events have occurred in C. elegans, whereas only gene conversions between zim-1 and zim-3 are detected in C. remanei. Although gene conversions are predicted to be present among zim-1, zim-2, and zim-3 in C. briggsae, the conversion tracts between zim-1/zim-2 and zim-2/zim-3 are very short. Moreover, positive selection analysis was performed on the basis of the significantly discordant phylogenies reconstructed using the N- and C-terminal sequences, respectively. Several codon sites located in the regions that are supposed not to have experienced gene conversions are predicted to be under the influence of positive selection. In comparison, stronger positive selection has acted on the C-terminal region relative to the N-terminal region. Thus, the zim/him-8 genes that evolve concertedly have also been shown to undergo adaptive diversifying selection. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Shedding Genomic Ballast: Extensive Parallel Loss of Ancestral Gene Families in Animals

Loss of ancestral gene families has played an important role in genomic specialization in animals. An examination of the pattern of gene family loss in completely sequenced animal genomes revealed that the same gene families have been lost independently in different lineages to a far greater extent than expected if gene loss occurred at random. This result implies that certain ancestral gene families—and thus the biological functions they encode—have been more expendable than others over the radiation of the animal phyla.Reviewing Editor: Dr. Manyuan Long     

Positive Selection Results in Frequent Reversible Amino Acid Replacements in the G Protein Gene of Human Respiratory Syncytial Virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a “flip-flop” phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.     

Evolution of the bovine lysozyme gene family: Changes in gene expression and reversion of function

Recruitment of lysozyme to a digestive function in ruminant artiodactyls is associated with amplification of the gene. At least four of the approximately ten genes are expressed in the stomach, and several are expressed in nonstomach tissues. Characterization of additional lysozymelike sequences in the bovine genome has identified most, if not all, of the members of this gene family. There are at least six stomachlike lysozyme genes, two of which are pseudogenes. The stomach lysozyme pseudogenes show a pattern of concerted evolution similar to that of the functional stomach genes. At least four nonstomach lysozyme genes exist. The nonstomach lysozyme genes are not monophyletic. A gene encoding a tracheal lysozyme was isolated, and the stomach lysozyme of advanced ruminants was found to be more closely related to the tracheal lysozyme than to the stomach lysozyme of the camel or other nonstomach lysozyme genes of ruminants. The tracheal lysozyme shares with stomach lysozymes of advanced ruminants the deletion of amino acid 103, and several other adaptive sequence characteristics of stomach lysozymes. I suggest here that tracheal lysozyme has reverted from a functional stomach lysozyme. Tracheal lysozyme then represents a second instance of a change in lysozyme gene expression and function within ruminants. Correspondence to: D.M. Irwin     

Positive Selection in the N-Terminal Extramembrane Domain of Lung Surfactant Protein C (SP-C) in Marine Mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (ω) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Richard Kliman     

Selection, mutation and sexual reproduction in an infinite haploid population with a genome of finite length

We present a model which describes mutation, selection and sexual reproduction in an infinite haploid population with a finite genome. Each generation is described using an approximation which assures a certain persistent form of the distribution of the number of deleterious elements. The steady state exists and is determined. In addition, we conclude that the introduction of sexual reproduction increases the mean fitness in the equilibrium.     

Molecular Evolution and Positive Selection of the Symbiotic Gene NORK in Medicago truncatula

Understanding the selective constraints of partner specificity in mutually beneficial symbiosis is a significant, yet largely unexplored, prospect of evolutionary biology. These selective constraints can be explored through the study of nucleotide polymorphism at loci controlling specificity. The membrane-anchored receptor NORK (nodulation receptor kinase) of the legume Medicago truncatula controls early steps of root infection by two symbiotic microorganisms: nitrogen-fixing bacteria (rhizobia) and endomycorrhizal fungi (Glomales). We analyzed the diversity of the gene NORK by sequencing 4 kilobases in 28 inbred lines sampled from natural populations. We detected 33 polymorphic sites with only one nonsynonymous change. Analysis based on Tajima’s D and Fay and Wu’s H summary statistics revealed no departure from the neutral model. We analyzed divergence using sequences from the closely related species M. coerulea. The McDonald-Kreitman test indicated a significant excess of nonsynonymous changes contributing to this divergence. Furthermore, maximum-likelihood analysis of a molecular phylogeny of a few legume species indicated that a number of amino acid sites, likely located in the receptor domain of the protein, evolved under the regime of positive selection. Further research should focus on the rate and direction of molecular coevolution between microorganisms’ signaling molecules and legumes’ receptors. [Reviewing Editor: Dr. Deborah Charlesworth] Sequence data were deposited in the GenBank database under accession nos. AY676428 to AY676457 and AJ884582.     

Detection of DNA in Ancient Bones Using Histochemical Methods

We describe histochemical techniques for detecting DNA within the osteocytic lacunae of ancient bones. The bones examined were fragments of femurs from two human individuals found in the Pompeian C. I. Polybius house and fragments of metacarpals from two horses (Equus sp.) found in the Pompeian “Casti Amanti” house. Both buildings were buried by the 79 A. D. Vesuvius eruption. Fragments of femurs from a modern horse, a modern swine and a modern amphibian also were studied as controls. Some bone sections were stained with two different DNA-specific fluorochromes, 4′ -′ 6-diamidino-2-phenylindole (DAPI) and chromomycin A3 (CMA), while others were stained by the Feulgen reaction. All of the techniques gave a positive reaction within the osteocytic lacunae. Histological analysis of the undecalcified, ground and unstained sections agreed well with results of bone sections stained with either the fluorochromes or the Feulgen reaction. Bones showing good histology also were positive by our DNA-specific stain. Histochemical and histological analyses correlated well with the success of DNA extraction and amplification. Using conventional DNA-specific histochemical techniques in conjunction with histological analysis can be useful in the study of DNA extracted from ancient bone remains while reducing both the amount of time and cost.     

Detection of DNA in Ancient Bones Using Histochemical Methods

We describe histochemical techniques for detecting DNA within the osteocytic lacunae of ancient bones. The bones examined were fragments of femurs from two human individuals found in the Pompeian C. I. Polybius house and fragments of metacarpals from two horses (Equus sp.) found in the Pompeian “Casti Amanti” house. Both buildings were buried by the 79 A. D. Vesuvius eruption. Fragments of femurs from a modern horse, a modern swine and a modern amphibian also were studied as controls. Some bone sections were stained with two different DNA-specific fluorochromes, 4' -' 6-diamidino-2-phenylindole (DAPI) and chromomycin A3 (CMA), while others were stained by the Feulgen reaction. All of the techniques gave a positive reaction within the osteocytic lacunae. Histological analysis of the undecalcified, ground and unstained sections agreed well with results of bone sections stained with either the fluorochromes or the Feulgen reaction. Bones showing good histology also were positive by our DNA-specific stain. Histochemical and histological analyses correlated well with the success of DNA extraction and amplification. Using conventional DNA-specific histochemical techniques in conjunction with histological analysis can be useful in the study of DNA extracted from ancient bone remains while reducing both the amount of time and cost.