Pattern Of Nucleotide Substitution Research Articles

MHC (major histocompatibility complex)-DRB genes and polymorphisms in common marmoset.

A New World monkey, the common marmoset (Callithrix jacchus), will be used as a preclinical animal model to study the feasibility of cell and gene therapy targeting immunological and hematological disorders. For elucidating the immunogenetic background of common marmoset to further studies, in the present study, polymorphisms of MHC-DRB genes in this species were examined. Twenty-one Caja-DRB exon 2 alleles, including seven new ones, were detected by means of subcloning and the polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) methods followed by nucleotide sequencing. Based on the alignment of these allele sequences, we designed two pairs of specific primers and established a PCR-SSCP method for DNA-based histocompatibility typing of the common marmoset. According to the family segregation data and phylogenetic analyses, we presumed that Caja-DRB alleles could be classified into five different loci. Southern blotting analysis also supported the existence of multiple DRB loci. The patterns of nucleotide substitutions suggests that positive selection operates in the antigen-recognition sites of Caja-DRB genes.

Molecular Evolution of the Histone 3 Multigene Family in the Drosophila melanogaster Species Subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3′ region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions.

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-<T and G<->C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.

Comparative Genomics of Mitochondrial DNA in Members of the Drosophila melanogaster Subgroup

In this study, a comparative genomics approach is employed to investigate the forces that shape evolutionary change in the mitochondrial DNA (mtDNA) of members of the Drosophila melanogaster subgroup. This approach facilitates differentiation of the patterns of variation resulting from processes acting at a higher level from those acting on a single gene. The mitochondrial genomes of three isofemale lines of D. simulans (siI, -II, and -III), two of D. melanogaster (Oregon R and a line from Zimbabwe), and D. mauritiana (maI and -II), and one of D. sechellia were sequenced and compared with that derived from D. yakuba. Data presented here indicate that at least three broad mechanisms shape the evolutionary dynamics of mtDNA in these taxa. The first set of mechanisms is intrinsic to the molecule. Dominant processes may be interpreted as selection for an increased rate of replication of the mtDNA molecule, biases in DNA repair, and differences in the pattern of nucleotide substitution among strands. In the genes encoded on the major strand (62% of the coding DNA) changes to or from C predominate, whereas on the minor changes to or from G predominate. The second set of mechanisms affects distinct lineages. There are evolutionary rate differences among lineages, possibly owing to population demographic changes or changes in mutational biases. This is supported by the heterogeneity found in synonymous, nonsynonymous, and silent substitutions. The third set of mechanisms differentially affects distinct genes. A maximum-likelihood sliding-window analysis detected four disjunct regions that have a significantly different nucleotide substitution process from that derived from the complete sequence. These data show the potential for comparative genomics to tease apart subtle forces that shape the evolution of DNA.

A MOLECULAR PHYLOGENY OF NORTH AMERICAN PLEUROCERIDAE (GASTROPODA: CERITHIOIDEA) BASED ON MITOCHONDRIAL 16S rDNA SEQUENCES

The Pleuroceridae Fischer, 1885, is one of three freshwater gastropod families currently recognized in the superfamily Cerithioidea Ferussac, 1819 (Mollusca: Caenogastropoda Cox, 1960). Despite considerable literature justifying various proposed generic names of North American pleurocerids, no study has been conducted examining phylogenetic relationships of the recognized genera. In an effort to expand our understanding of evolutionary relationships of North American pleurocerid genrea, we examined a large portion of the mitochondrial 16S rRNA gene among 32 extant North American taxa. Multiple sequence alignment of the amplified region for our taxa resulted in a matrix consisting of 900 nucleotides including insertions and deletions. Based on analysis of nucleotide substitution patterns, we employed two approaches in our phylogenetic analysis: (1) all substitutions received equal weighting and (2) transversions were weighted 2X and 4X transitions to compensate for transition saturation among distantly related taxa. The molecular phylogeny based on the mitochondrial 16S rDNA sequences supports the monophyly of Pleurocera Rafinesque, 1819, Elimia H.&A. Adams, 1854, and Juga H.&A. Adams, 1854, but depicts the genera Lithasia Haldeman, 1840, and Leptoxis Rafinesque, 1819, as polyphyletic. The genus Pleurocera is sister to Elimia, which in turn is sister to a paraphyletic assemblage including representatives of Leptoxis, Lithasia, and the monotypic genus Io Lea, 1831. Juga, a genus restricted to west of the North American continental divide is the basal-most clade and is sister to all the aforementioned genera found east of the continental divide.

Evolution of nucleotide substitutions and gene regulation in the amylase multigenes in Drosophila kikkawai and its sibling species.

In order to determine evolutionary changes in gene regulation and the nucleotide substitution pattern in a multigene family, the amylase multigenes were characterized in Drosophila kikkawai and its sibling species. The nucleotide substitution pattern was investigated. Drosophila kikkawai has four amylase genes. The Amy1 and Amy2 genes are a head-to-head duplication in the middle of the B arm of the second chromosome, while the Amy3 and Amy4 genes are a tail-to-tail duplication near the centromere of the same chromosome. In the sibling species of D. kikkawai (Drosophila bocki, Drosophila leontia, and Drosophila lini), sequencing of the Amy1, Amy2, Amy3, and Amy4 genes revealed that the Amy1 and Amy2 gene group diverged from Amy3 and Amy4 after duplication. In the Amy1 and Amy2 genes, the divergent evolution occurred in the flanking regions; in contrast, the coding regions have evolved in concerted fashion. The electrophoretic pattern of AMY isozymes was also examined. In D. kikkawai and its siblings, two or three electrophoretically different isozymes are encoded by the Amy1 and Amy2 genes (S isozyme) and by the Amy3 and Amy4 genes (F (M) isozymes). The S and F (M) isozymes show different patterns of band intensity when larvae and flies were fed in different media. Amy1 and Amy2, which encode the S isozyme, are more strikingly regulated than Amy3 and Amy4, which encode the F (M) isozyme. The GC content and codon usage bias were higher for the Amy1 and Amy2 genes than for the Amy3 and Amy4 genes. Although the ratio of synonymous and replacement substitutions within the Amy1 and Amy2 gene group was not significantly different from that within the Amy3 and Amy4 gene group, the synonymous substitution rate in the lineage of Amy1 and Amy2 was lower than that of Amy3 and Amy4. In conclusion, after the first duplication but before speciation of four species, the synonymous substitution rate between the two lineages and the electrophoretic pattern of the isozymes encoded by them changed, although we do not know whether there was any evolutionary relationship between the two.

Adaptive diversification within a large family of recently duplicated, placentally expressed genes.

The pregnancy-associated glycoproteins (PAG) are putative peptide-binding proteins and products of a large family of genes whose expression is localized to the placental surface epithelium of artiodactyl species. We have tested the hypothesis that natural selection has favored diversification of these genes by examining patterns of nucleotide substitution in a sample of 28 closely related bovine, caprine, and ovine family members that are expressed only in trophoblast binucleate cells. Three observations were made. First, in codons encoding highly variable domains of the proteins, there was a greater accumulation of both synonymous and nonsynonymous mutations than in the more conserved regions of the genes. Second, in the variable regions, the mean number of nonsynonymous nucleotide substitutions per site was significantly greater than the mean number of synonymous substitutions per site. Third, nonsynonymous changes affecting amino acid charge occurred more frequently than expected under random substitution. This unusual pattern of nucleotide substitution implies that natural selection has acted to diversify these PAG molecules at the amino acid level, which in turn suggests that these molecules have undergone functional diversification. We estimate that the binucleate cell-expressed PAG originated 52 +/- 6 million years ago, soon after the divergence of the ruminant lineage. Thus, rapid functional diversification of PAG expressed in trophoblast binucleate cells seems to have been associated with the origin of this unique placental adaptation.

Adaptive diversification within a large family of recently duplicated, placentally expressed genes

The pregnancy-associated glycoproteins (PAG) are putative peptide-binding proteins and products of a large family of genes whose expression is localized to the placental surface epithelium of artiodactyl species. We have tested the hypothesis that natural selection has favored diversification of these genes by examining patterns of nucleotide substitution in a sample of 28 closely related bovine, caprine, and ovine family members that are expressed only in trophoblast binucleate cells. Three observations were made. First, in codons encoding highly variable domains of the proteins, there was a greater accumulation of both synonymous and nonsynonymous mutations than in the more conserved regions of the genes. Second, in the variable regions, the mean number of nonsynonymous nucleotide substitutions per site was significantly greater than the mean number of synonymous substitutions per site. Third, nonsynonymous changes affecting amino acid charge occurred more frequently than expected under random substitution. This unusual pattern of nucleotide substitution implies that natural selection has acted to diversify these PAG molecules at the amino acid level, which in turn suggests that these molecules have undergone functional diversification. We estimate that the binucleate cell-expressed PAG originated 52 +/- 6 million years ago, soon after the divergence of the ruminant lineage. Thus, rapid functional diversification of PAG expressed in trophoblast binucleate cells seems to have been associated with the origin of this unique placental adaptation.

Examining rates and patterns of nucleotide substitution in plants.

Driven by rapid improvements in affordable computing power and by the even faster accumulation of genomic data, the statistical analysis of molecular sequence data has become an active area of interdisciplinary research. Maximum likelihood methods have become mainstream because of their desirable properties and, more importantly, their potential for providing statistically sound solutions in complex data analysis settings. In this chapter, a review of recent literature focusing on rates and patterns of nucleotide substitution rates in the nuclear, chloroplast, and mitochondrial genomes of plants demonstrates the power and flexibility of these new methods. The emerging picture of the nucleotide substitution process in plants is a complex one. Evolutionary rates are seen to be quite variable, both among genes and among plant lineages. However, there are hints, particularly in the chloroplast, that individual factors can have important effects on many genes simultaneously.

Evolutionary diversification of the mammalian defensins.

Defensins are cysteine-rich cationic peptides that function in antimicrobial defense in both invertebrates and vertebrates. Three main groups of animal defensins are known: insect defensins; mammalian alpha-defensins and vertebrate beta-defensins. It has been difficult to determine whether these molecules are homologous or have independently evolved similar features, but overall the evidence favors a distant relationship. The best evidence of this relationship is structural, particularly from their overall three-dimensional structure and from the spacing of half-cystine residues involved in intra-chain disulfide bonds. Some evidence favors a closer relationship between vertebrate beta-defensins and insect defensins than between the two groups of vertebrate defensins. Examination of nucleotide substitutions between recently duplicated mammalian defensin genes shows that the rate of nonsynonymous (amino-acid-altering) substitution exceeds that of synonymous substitution in the region of the gene encoding the mature defensin. This highly unusual pattern of nucleotide substitution is evidence that natural selection has acted to diversify defensins at the amino acid level. The resulting rapid evolution explains why it is difficult to reconstruct the evolutionary history of these molecules.

Pattern of nucleotide substitution and rate heterogeneity in the hypervariable regions I and II of human mtDNA.

This study provides a comprehensive survey of the complex pattern of nucleotide substitution in the control region of human mtDNA, which is of central importance to the studies of human evolution. A total of 1229 different hypervariable region I (HVRI) and 385 different hypervariable region II (HVRII) sequences were analyzed using a complex substitution model. Moreover, we suggest a new method to assign relative rates to each site in the sequence. Estimates are based on maximum-likelihood methods applied to randomly selected subsets of sequences. Our results indicate that the rate of substitution in HVRI is approximately twice as high as in HVRII and that this difference is mainly due to a higher frequency of pyrimidine transitions in HVRI. However, rate heterogeneity is more pronounced in HVRII.

The complete sequence of the mitochondrial genome of Daphnia pulex (Cladocera: Crustacea).

The sequence of the mitochondrial DNA (mtDNA) of the branchiopod crustacean Daphnia pulex has been completed. It is 15333bp with an A+T content of 62.3%, and contains the typical complement of 13 protein-coding, 22 transfer RNA (tRNA) and two ribosomal RNA (rRNA) genes. Comparison of this sequence with the sequences of the other eight completely sequenced arthropod mtDNAs showed that gene order and orientation are identical to that of Drosophila but different from Artemia due to the rearrangement of two tRNA genes. Nucleotide composition, codon usage, and amino acid composition are very similar in the crustaceans, but divergent from insects and chelicerates which show a much higher bias towards A+T. However, with few exceptions, the mitochondrial proteins of Daphnia are more similar to those of the dipteran insects (Drosophila and Anopheles) than to those of Artemia, at both the nucleotide and amino acid levels, suggesting that Artemia mtDNA is evolving at an accelerated rate. These results also show that sequence evolution and the evolution of nucleotide composition can be decoupled. Analysis of nucleotide substitution patterns in COII showed that there has been an unbiased acceleration of the overall substitution rate in Artemia. In contrast, the accelerated substitution rate in Apis is due partly to extreme A+T mutation pressure. Secondary structures are proposed for the Daphnia tRNAs and rRNAs. The tRNAs are similar to those of other arthropods but tend to have TPsiC arms that are only 4bp long. The rRNA secondary structures are similar to those proposed for insects except for the absence of a small number of helices in Daphnia. Phylogenetic analysis of second codon positions grouped Daphnia with Artemia, as expected, despite the latter's accelerated divergence rate. In contrast, the unusual pattern of mtDNA divergence in Apis led to a topology in which the holometabolous insects (Anopheles, Drosophila, Apis) appeared to be paraphyletic with respect to the hemimetabolous insect, Locusta, due to the early branching of Apis.

Patterns of nucleotide substitution in Drosophila and mammalian genomes.

To estimate patterns of molecular evolution of unconstrained DNA sequences, we used maximum parsimony to separate phylogenetic trees of a non-long terminal repeat retrotransposable element into either internal branches, representing mainly the constrained evolution of active lineages, or into terminal branches, representing mainly nonfunctional "dead-on-arrival" copies that are unconstrained by selection and evolve as pseudogenes. The pattern of nucleotide substitutions in unconstrained sequences is expected to be congruent with the pattern of point mutation. We examined the retrotransposon Helena in the Drosophila virilis species group (subgenus Drosophila) and the Drosophila melanogaster species subgroup (subgenus Sophophora). The patterns of point mutation are indistinguishable, suggesting considerable stability over evolutionary time (40-60 million years). The relative frequencies of different point mutations are unequal, but the "transition bias" results largely from an approximately 2-fold excess of G.C to A.T substitutions. Spontaneous mutation is biased toward A.T base pairs, with an expected mutational equilibrium of approximately 65% A + T (quite similar to that of long introns). These data also enable the first detailed comparison of patterns of point mutations in Drosophila and mammals. Although the patterns are different, all of the statistical significance comes from a much greater rate of G.C to A.T substitution in mammals, probably because of methylated cytosine "hotspots." When the G.C to A.T substitutions are discounted, the remaining differences are considerably reduced and not statistically significant.

Variation in the pattern of nucleotide substitution across sites.

A model of nucleotide substitution that allows the transition/transversion rate bias to vary across sites was constructed. We examined the fit of this model using likelihood-ratio tests by analyzing 13 protein coding genes and 1 pseudogene. Likelihood-ratio testing indicated that a model that allows variation in the transition/transversion rate bias across sites provided a significant improvement in fit for most protein coding genes but not for the pseudogene. When the analysis was repeated with parameters estimated separately for first, second, and third codon positions, strong heterogeneity was uncovered for the first and second codon positions; the variation in the transition/transversion rate was generally weaker at the third codon position. The transition rate bias and branch lengths are underestimated when variation in the transition/transversion rate was not accommodated, suggesting that it may be important to accommodate variation in the pattern of nucleotide substitution for accurate estimation of evolutionary parameters.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution.

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Patterns of Y and X Chromosome DNA Sequence Divergence During the Felidae Radiation

The 37 species of modern cats have evolved from approximately eight phylogenetic lineages within the past 10 to 15 million years. The Felidae family has been described with multiple measures of morphologic and molecular evolutionary methods that serve as a framework for tracking gene divergence during brief evolutionary periods. In this report, we compare the mode and tempo of evolution of noncoding sequences of a large intron within Zfy (783 bp) and Zfx (854 bp), homologous genes located on the felid Y and X chromosomes, respectively. Zfy sequence variation evolves at about twice the rate of Zfx, and both gene intron sequences track feline hierarchical topologies accurately. As homoplasies are infrequent in patterns of nucleotide substitution, the Y chromosome sequence displays a remarkable degree of phylogenetic consistency among cat species and provides a highly informative glimpse of divergence of sex chromosome sequences in Felidae.

Intralineage variation in the pattern ofrbcL nucleotide substitution

Variation in chloroplastrbcL sequences was studied in representative species of four different lineages: the tribeRubieae (Rubiaceae), and the generaDrosera (Droseraceae),Nothofagus (Nothofagaceae) andIlex (Aquifoliaceae). Each lineage has its particular non-overlapping set ofrbcL polymorphic sites, indicating that common unconstrainedrbcL sites are not shared. Large differences in the rate and pattern of nucleotide substitution are observed among the four lineages. The genusIlex has the lowest rate of substitution, the lowest transition/transversion ratio, the lowest synonymous/replacement ratio and the lowest number of substitutions at the third codon position. An apparent relationship of these measures to the age of the lineages is observed. The A + T content and codon use among the four lineages are very similar and, apparently, cannot account for the observed differences in patterns of nucleotide substitution. However, the A + T content of the two bases immediately flanking the polymorphic sites is higher inIlex than in the other lineages. This could be correlated with the transversion/transition bias observed inIlex. The particularly low synonymous/replacement ratio found inIlex could also be explained by the small population sizes of species in this genus.

Methods for the detection of non-random base substitution in virus genes: models of synonymous nucleotide substitution in picornavirus genes.

A substantial fraction of phylogenetic divergence between closely related RNA virus genes is generally accounted for by synonymous (non-amino acid changing) point mutation. Viral evolution may be a complicated phenomena, governed by many different processes. However in this study we ask whether there are any properties in the patterns of synonymous nucleotide substitutions in three different Picornavirus genes that permit the process of accumulation of synonymous point mutation in these genes to be distinguished from some of the simplest most basic evolutionary models. We conclude that while the observed patterns in the occurrence of synonymous point substitution are consistent with those predicted by a model in which base mutation is equi-probable along a gene, and the probability of synonymous substitution determined only by local codon usage, some patterns in the actual nucleotides exchanged remain to be explained.

Nucleotide Sequence Evolution at the κ-Casein Locus: Evidence for Positive Selection Within the Family Bovidae

Kappa-casein is a mammalian milk protein involved in a number of important physiological processes. In the gut, the ingested protein is split into an insoluble peptide (para kappa-casein) and a soluble hydrophilic glycopeptide (caseinomacropeptide). Caseinomacropeptide is responsible for increased efficiency of digestion, prevention of neonate hypersensitivity to ingested proteins, and inhibition of gastric pathogens. Variation within this peptide has significant effects associated with important traits such as milk production. The nucleotide sequences for regions of kappa-casein exon and intron four were determined for representatives of the artiodactyl family Bovidae. The pattern of nucleotide substitution in kappa-casein sequences for distantly related bovid taxa demonstrates that positive selection has accelerated their divergence at the amino acid sequence level. This selection has differentially influenced the molecular evolution of the two kappa-casein split peptides and is focused within a 34-codon region of caseinomacropeptide.

Bovine IgG repertoire is dominated by a single diversified VH gene family.

Using methods based upon the PCR, we have analyzed sequences that encode the V domain of bovine Ig heavy chain. Pairwise comparison indicates that the heavy chain repertoire in cattle is dominated by a single VH gene family that is homologous to human VHII. The third complementarity-determining region (CDR3) is particularly variable and 60% of the sequences encode a CDR3 in excess of 21 amino acids. Southern blots, single-stranded conformation polymorphism analysis, and sequencing of segments from the germline indicate that the dominant VHII gene family is of small size, comprising as few as 10 unique members, and carries a very low level of diversity in CDRs 1 and 2. Comparison of germline VH sequences with those gathered from adult bovine splenic cDNA demonstrates a pattern of nucleotide substitution that is consistent with diversification through somatic hypermutation.