Search

To study the difference in transcriptome level of fatty acid metabolism pathway in Bamei pork and the difference of pork quality caused by the difference. In this study, Bamei pigs breeding in Huzhu farm of QingHai province were selected as the test object, compared with Gansu Black pigs. Four indexes of nutmeg acid (DX1), palmitic acid (DX2), stearic acid (DX3) and linoleic acid (DX4) were set. The expression profiles of fat metabolism related genes between the two groups samples were analysed by GCMS metabolomics and transcriptomics, then coexpression network analysis were conducted to obtain phenotypic related genes. The results showed that the metabolic levels of DX3 and DX4 were significantly higher than those of other fatty acids. Among these differences, the ENSSSCG00000024681 (G1) and ENSSSCG00000036883 (G2) genes play important regulatory roles in fatty acid metabolism, and the upregulated expression of their gene obviously affects the level of fatty acid metabolism, thereby affecting the quality and taste of pork. In addition, we found that there was a good correlation between the same lines, and the genetic traits of the hybrid lines of Bamei pig and Black pig are more inclined to Bamei pig. In the independent fatty acid metabolism, "Mg 2+ "and flavin adenine dinucleotide are more active, which plays an important role in energy utilization. Therefore, we can be inferred that the metabolism of stearic acid and linoleic acid are important fatty acids for pork quality. It also further confirms that the research method of combined omics is of great significance for the study of species traits and gene functions.

A method is described for the detection of Cu(II). It is based on the use of a room-temperature phosphorescent probe consisting of alginate-capped and manganese(II)-doped ZnS quantum dots. The carboxy groups at the surface of the probe strongly coordinate Cu(II) to form a complex. As a result, the 4 T 1 - 6 A 1 transition of the Mn(II) ions in the probe is quenched, and the long decay time (~2.1 ms in the unquenched state) is accordingly reduced. At excitation/emission wavelengths of 316/590 nm and a delay time of 0.1 ms, the probe shows a linear response in the 0.01 to 12 μM Cu(II) concentration range. The detection limit is 6.0 nM and the RSD is 3.2% (for n = 5). Graphical Abstract A two-step procedure is described to synthesize alginate capped manganese doped ZnS QDs. These coordinate with Cu(II) to form an absorbent complex and can be used as a phosphorescent probe for time-resolved detection of Cu(II).

Purpose: To better understand the gene mutational status and heterogeneity between primary and metastatic CRC (mCRC) using a sensitive sequencing method., Methods: The mutational status of EGFR, KRAS, NRAS, PIK3CA, ERBB2, BRAF, KIT, and PDGFRA was analyzed in 65 patients, with 147 samples of primary and paired live or lung metastatic CRC, using next-generation sequencing (NGS), quantitative RT-PCR (qPCR), and Sanger sequencing., Results: Fifteen cases (15/22, 68.2%) of lung mCRC and thirteen cases (13/20, 65%) of liver mCRC harboured the same mutation profiles of KRAS, NRAS, or BRAF in the primary lesions. To all detected genes, 11 cases (11/22, 50%) of lung mCRC and 11 cases (11/20, 55%) of liver mCRC showed different mutational genes in the primary tumours. KRAS and BRAF mutations were more frequent in lung metastatic lesions (p = 0.004 and 0.003, respectively). The gene mutations in KRAS, NRAS, BRAF, and PIK3CA in the lung metastatic sites were more frequent than those in the liver metastatic sites (86.7 vs. 44%, respectively, p = 0.000). Some new mutations were not covered in the qPCR ranges but were detected by NGS., Conclusion: The study demonstrated that the discordance of gene mutational status between paired primary and metastatic tumours is rather high when detected by NGS. Evaluating the mutational status of both the primary and metastatic tumours should be considered in clinical mutation testing.

A method based on ultrasound-assisted emulsification-microextraction (USAEME) was proposed in this contribution for the determination of ethyl carbamate (EC) in alcoholic beverages using gas chromatography coupled to triple quadrupole mass spectrometry. To achieve the determination of EC in alcoholic beverages, the influences on the extraction efficiency of type and volume of extraction solvent, temperature, ionic strength, alcohol content, and extraction time were studied, once the extraction solvent had been selected. The optimized conditions were 200.0 μL of chloroform at 30 °C during 5 min with 15% (m/v) sodium chloride addition. The detection limit, relative standard deviations, linear range, and recoveries under the optimized conditions were 0.03 μg L(-1), 4.2-6.1%, 0.1-50.0 μg L(-1), and 80.5-87.9%, respectively. Moreover, the feasibility of the present method was also validated by real samples. To the best of our knowledge, this is the first time that USAEME has been applied to determine a strongly hydrophilic compound in alcoholic beverages.

Purpose: This study investigated the in-vivo formation process of laser-induced choroidal neovascularization (CNV) in rat using high-resolution spectral-domain optical coherence tomography (SD-OCT), and compared the results to histological methods., Methods: Brown Norway rats (n = 60, 6-8 weeks of age) received 532-nm diode laser photocoagulation. SD-OCT and fluorescein angiography (FA) were performed in vivo 2, 5, 7, 14, and 21 days post-laser application. Haematoxylin and eosin (H&E) staining and immunohistochemistry for CD31, phosphorylated vascular endothelial factor receptor 2 (pVEGFR2) were conducted at each time point to observe the CNV in vitro. Choroidal flatmount preparations were observed using a confocal laser scanning microscope (CLSM) and a scanning electron microscope (SEM)., Results: SD-OCT monitored the longitudinal morphological changes of laser-induced CNV. CNV reached its maximal size on day 7, and began a gradual reduction on day 14. FA revealed similar dynamic changes in leakage. CNV thickness, as assessed by SD-OCT, was consistent with H&E-stained sections at each time point. CLSM and SEM revealed the details of the fibrovascular membrane. CD31 and pVEGFR2 expression supported the results of SD-OCT and histology., Conclusions: SD-OCT was a convenient and reliable tool for the imaging of the CNV formation process and quantification of the lesion size in vivo.

We report the imaging findings of a mature thyroid teratoma in a 5-year-old girl. Nuclear imaging showed a decrease in (99)Tcm uptake in the right lobe of the thyroid gland. CT scan showed a slightly lobulated soft-tissue mass without calcification, fat or cystic components. Histological analysis showed that the tumor was composed of mature neural tissue, cartilaginous, and epithelial elements. This case study provides new insights into the CT appearance of mature thyroid teratomas.

The enzymes catalyzing the first two reactions in the sulfate activation pathway, ATP-sulfurylase (S) and APS-kinase (K), are fused as 'KS' in animals but are fused as 'SK' in select bacteria and fungi. We have discovered a novel triple fusion protein of K, S, and pyrophosphatase (P) in several protozoan genomes within the Stramenopile lineage. These triple domain fusion proteins led us to hypothesize that pyrophosphatase enzymes and sulfate activation enzymes physically interact to impact the thermodynamics of the sulfate activation pathway. In support of this hypothesis, we demonstrate through biochemical assays that separately encoded KS and P proteins physically interact and that KS/P complexes activate more sulfate than KS alone. We also conclude on the basis of phylogenetic analyses that all known KS fusions originate from a single fusion event early in the eukaryotic lineage. Strikingly, our analyses support the same conclusion for all known SK fusions. These observations indicate that the patchwork of fused and nonfused S and K genes observed in modern-day eukaryotes and prokaryotes are the result of the two ancestral fusion genes evolving by an assortment of gene fissions, duplications, deletions, and horizontal transfers in different lineages. Our integrative use of genomics, phylogenetics, and biochemistry to characterize pyrophosphatase as a new member of the sulfate activation pathway should be effective at identifying new protein members and connections in other molecular pathways.

Although the evolutionary significance of gene duplication has long been appreciated, it remains unclear what factors determine gene duplicability. In this study we investigated whether metabolism is an important determinant of gene duplicability because cellular metabolism is crucial for the survival and reproduction of an organism. Using genomic data and metabolic pathway data from the yeast (Saccharomyces cerevisiae) and Escherichia coli, we found that metabolic proteins indeed tend to have higher gene duplicability than nonmetabolic proteins. Moreover, a detailed analysis of metabolic pathways in these two organisms revealed that genes in the central metabolic pathways and the catabolic pathways have, on average, higher gene duplicability than do other genes and that most genes in anabolic pathways are single-copy genes.

We estimated the dates of the monocot-dicot split and the origin of core eudicots using a large chloroplast (cp) genomic dataset. Sixty-one protein-coding genes common to the 12 completely sequenced cp genomes of land plants were concatenated and analyzed. Three reliable split events were used as calibration points and for cross references. Both the method based on the assumption of a constant rate and the Li-Tanimura unequal-rate method were used to estimate divergence times. The phylogenetic analyses indicated that nonsynonymous substitution rates of cp genomes are unequal among tracheophyte lineages. For this reason, the constant-rate method gave overestimates of the monocot-dicot divergence and the age of core eudicots, especially when fast-evolving monocots were included in the analysis. In contrast, the Li-Tanimura method gave estimates consistent with the known evolutionary sequence of seed plant lineages and with known fossil records. Combining estimates calibrated by two known fossil nodes and the Li-Tanimura method, we propose that monocots branched off from dicots 140-150 Myr ago (late Jurassic-early Cretaceous), at least 50 Myr younger than previous estimates based on the molecular clock hypothesis, and that the core eudicots diverged 100-115 Myr ago (Albian-Aptian of the Cretaceous). These estimates indicate that both the monocot-dicot divergence and the core eudicot's age are older than their respective fossil records.

The hominoid primates (apes and humans) exhibit remarkable diversity in their social and sexual behavioral systems. This is reflected in many ways in their anatomy and physiology. For example, the testes and seminal vesicles are relatively large in species with high sperm competition like the chimpanzee and small in species with low or no sperm competition like the gorilla. Additionally, the chimpanzee is the only hominoid primate known to produce a firm copulatory plug, which presumably functions in sperm competition by blocking insemination of subsequent males. Here we examine the molecular evolution of the semenogelin genes (SEMG1 and SEMG2), which code for the predominant structural proteins in human semen. High molecular weight complexes of these proteins are responsible for the viscous gelatinous consistency of human semen; their rodent homologs are responsible for the formation of a firm copulatory plug. Chimpanzees have an expanded SEMG1 gene caused by duplications of tandem repeats, each encoding 60 amino acids, resulting in a protein nearly twice as long as that of humans. In contrast, at both SEMG1 and SEMG2 we observed several gorilla haplotypes that contain at least one premature stop codon. We suggest that these structural changes in the semenogelin proteins that have arisen since the human-chimpanzee-gorilla split may be responsible for the physiological differences between these species ejaculated semen that correlate with their sociosexual behavior.

It has been hypothesized that the length of an exon tends to increase with the GC content because stop codons are AT-rich and should occur less frequently in GC-rich exons. This prediction assumes that mutation pressure plays a significant role in the occurrence and distribution of stop codons. However, the prediction is applicable not to all exons, but only to the last coding exon of a gene and to single-exon CDS sequences. We classified exons in multiexon genes in eight eukaryotic species into three groups-the first exon, the internal, and the last exon-and computed the Spearman correlation between the exon length and the percentage GC (%GC) for each of the three groups. In only five of the species studied is the correlation for the last coding exon greater than that for the first or internal exons. For the single-exon CDS sequences, the correlation between CDS length and %GC is mostly negative. Thus, eukaryotic genomes do not support the predicted relationship between exon length and %GC. In prokaryotic genomes, CDS length and %GC are positively correlated in each of the 68 completely sequenced prokaryotic genomes in GenBank with genomic GC contents varying from 25 to 68%, except for the wall-less Mycoplasma genitalium and the syphilis pathogen Treponema pallidum. Moreover, the average CDS length and the genomic GC content are also positively correlated. After correcting for genome size, the partial correlation between the average CDS length and the genomic GC content is 0.3217 ( p < 0.025).

In this paper we present a new method for detecting block duplications in a genome. It is more stringent than previous ones in that it requires a more rigorous definition of paralogous genes and that it requires the paralogous proteins on the two blocks to be contiguous. In addition, it provides three criterion choices: (1) the same composition (i.e., having the same paralogues in the two windows), (2) the same composition and gene order, and (3) the same composition, gene order, and gene orientation. The method is completely automated, requiring no visual inspection as in previous methods. We applied it to analyze the complete genomes of S. cerevisiae and C. elegans. In yeast we detected fewer duplicated blocks than previously reported. In C. elegans, however, we detected more block duplications than previously reported, indicating that although our method has a more stringent definition of block duplication than previous ones, it may be more sensitive in detection because it considers every possible window rather than only fixed nonoverlapping windows. Our results show that block duplication is a common phenomenon in both organisms. The patterns of block duplication in the two species are, however, markedly different. The yeast shows much more extensive block duplication than the nematode, with some chromosomes having more than 40% of the duplications derived from block duplications. Moreover, in the yeast the majority of block duplications occurred between chromosomes, while in the nematode most block duplications occurred within chromosomes.

It has been known that in noncoding regions of the chloroplast genome, the pattern of nucleotide substitution is influenced by the two nucleotides flanking the substitution site. In a GC-rich environment, a bias toward transition was observed, whereas in an AT-rich environment, a bias toward transversion was observed. In this study, the influence of the two adjacent neighbors on the substitution pattern was observed in the first intron of the mitochondrial nad4 gene, although the AT content of this intron is only 48%. The proportion of transversions increases from 0.32 to 0.75 as the A + T content (number of A's + T's) of the two nearest neighbors increases from 0 to 2. This trend was also observed in another mitochondrial group I intron with an AT content of 64%. In addition, a similar, though weaker, effect was observed in vertebrate pseudogenes. So this effect is present in all three types of genomes. Furthermore, in contrast to the situation in the noncoding regions of chloroplast DNA, where most nucleotide substitutions occurred in the categories with an A + T content of either 1 or 2, nucleotide substitutions in the mitochondrial first nad4 intron occurred more evenly in three categories of different A + T contents. This might be due largely to the difference in the AT content (0.48 vs. 0.72) between the mitochondrial first nad4 intron and the chloroplast DNA regions studied.

Most Neotropical primate species possess a polymorphic X-linked and a monomorphic autosomal color vision gene. Consequently, populations are composed of both dichromatics and trichromatics. Most theories on the maintenance of this genetic system revolve around possible advantages for foraging ecology. To examine the issue from a different angle, we compared the numbers and relative frequencies of alleles at the X-linked locus among three species of Saimiri representing a wide range of geographical and behavioral variation in the genus. Exons 3, 4, and 5 of the X-linked opsin gene were sequenced for a large number of X chromosomes for all three species. Several synonymous mutations were detected in exons 4 and 5 for the originally reported alleles but only a single nonsynonymous change was detected. Two alleles were found that appeared to be the result of recombination events. The low occurrence of recombinant alleles and absence of mutations in the amino acids critical for spectral tuning indicates that stabilizing selection acts to maintain the combinations of critical sites specific to each allele. Allele frequencies were approximately the same for all Saimiri species, with a slight but significant difference between S. boliviensis and S. oerstedii. No apparent correlation exists between allele frequencies and behavioral or biogeographical differences between species, casting doubt on the speculation that the spectral sensitivities of the alleles have been maintained because they are specifically well-tuned to Saimiri visual ecology. Rather, the spectral tuning peaks might have been maintained because they are as widely spaced as possible within the limited range of middlewave to longwave spectra useful to all primates. This arrangement creates a balance between maximizing the distance between spectral tuning peaks (allowing the color opponency of the visual system to distinguish between peaks) and maximizing the number of alleles within a limited range (yielding the greatest possible frequency of heterozygotes).

We obtained 16 nucleotide sequences ( approximately 1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16-0.23 x 10(-9) substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5-0.7 x 10(-9) substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170-235 million years (Myr) for the monocot-dicot split, 112-156 Myr for the Brassicaceae-Lettuce split, 14.5-20.4 Myr for the Brassica-Arabidopsis split, and 14.5-20.4 Myr for the Arabidopsis-Arabideae split.

We studied 10 protein-coding mitochondrial genes from 19 mammalian species to evaluate the effects of 10 amino acid properties on the evolution of the genetic code, the amino acid composition of proteins, and the pattern of nonsynonymous substitutions. The 10 amino acid properties studied are the chemical composition of the side chain, two polarity measures, hydropathy, isoelectric point, volume, aromaticity, aliphaticity, hydrogenation, and hydroxythiolation. The genetic code appears to have evolved toward minimizing polarity and hydropathy but not the other seven properties. This can be explained by our finding that the presumably primitive amino acids differed much only in polarity and hydropathy, but little in the other properties. Only the chemical composition (C) and isoelectric point (IE) appear to have affected the amino acid composition of the proteins studied, that is, these proteins tend to have more amino acids with typical C and IE values, so that nonsynonymous mutations tend to result in small differences in C and IE. All properties, except for hydroxythiolation, affect the rate of nonsynonymous substitution, with the observed amino acid changes having only small differences in these properties, relative to the spectrum of all possible nonsynonymous mutations.

Although most New World monkeys have only one X-linked photopigment locus, many species have three polymorphic alleles at the locus. The three alleles in the squirrel monkey and capuchin have spectral peaks near 562, 550, and 535 nm, respectively, and the three alleles in the marmoset and tamarin have spectral peaks near 562, 556, and 543 nm, respectively. To determine the amino acids responsible for the spectral sensitivity differences among these pigment variants, we sequenced all exons of the three alleles in each of these four species. From the deduced amino acid sequences and the spectral peak information and from previous studies of the spectral tuning of X-linked pigments in humans and New World monkeys, we estimated that the Ala --> Ser, Ile --> Phe, Gly --> Ser, Phe --> Tyr, and Ala --> Tyr substitutions at residue positions 180, 229, 233, 277, and 285, respectively, cause spectral shifts of about 5, -2, -1, 8, and 15 nm. On the other hand, the substitutions His --> Tyr, Met --> Val or Leu, and Ala --> Tyr at positions 116, 275, and 276, respectively, have no discernible spectral tuning effect, though residues 275 and 276 are inside the transmembrane domains. Many substitutions between Val and Ile or between Val and Ala have occurred in the transmembrane domains among the New World monkey pigment variants but apparently have no effect on spectral tuning. Our study suggests that, in addition to amino acid changes involving a hydroxyl group, large changes in residue size can also cause a spectral shift in a visual pigment.

Bush babies have had a long history of nocturnal life and it would be interesting to know whether their color vision genes have become degenerate. Therefore, we used PCR techniques to sequence the X-linked pigment gene of two of these nocturnal prosimians: Galago senegalensis and Otolemur garnettii. Southern hybridization of genomic DNA of G. senegalensis showed a single X-linked pigment gene. Interestingly, the deduced pigment sequences of the two bush babies are identical. By comparing the X-linked pigments of bush baby, human, squirrel monkey, and marmoset, 38 variable positions were identified. At those positions that may cause a spectral shift, the bush baby pigment has identical or biochemically similar residues to those of the marmoset cone pigment with a spectral peak of 543 nm. This result is consistent with the estimate of 544-545 nm for the spectral peak of the X-linked pigment of Otolemur crassicaudatus, which is closely related to Otolemur garnettii. The neighbor-joining tree of mammalian X-linked pigments showed a significantly shorter branch in the bush baby lineage than in other primate lineages. A relative rate test showed that the nonsynonymous substitution rate of the bush baby X-linked pigment gene is about three times slower than that of the human red pigment gene, though the synonymous substitution rates of the two genes are similar. The slower nonsynonymous rate in the bush baby lineage suggests that the bush baby X-linked pigment gene is under functional constraints, in spite of its nocturnal life. Two radical changes at positions in the intradiskal surface next to the sixth transmembrane domain were observed in the X-linked cone pigment of bush babies but not in other primates. They are changes from Ala to Ser and from Asn to His, which are similar in function to the corresponding residues in rhodopsins. These two changes may be of importance for dim light sensitivity, which is consistent with our proposal that the evolution of the bush baby X-linked pigment gene is under selective pressure. In addition, the 2.5% divergence in introns 2 and 5 of the X-linked pigment gene between the two bush babies supports their classification into two separate genera.

The sequences of the entire blue opsin gene in the squirrel monkey (Saimiri boliviensis) and the five introns of the human blue opsin gene were obtained. Intron 3 of these genes contains an Alu sequence and intron 4 contains a partial mer13 sequence. A comparison of the squirrel monkey opsin sequence with published mammalian opsin sequences shows that features believed to be functionally critical are all conserved. However, the blue opsin has evolved twice as fast as rhodopsin and is only as conservative as the beta globin, which has evolved at the average rate of mammalian proteins. Interestingly, the interhelical loops are, on average, actually more conservative than the transmembrane alpha helical regions. The introns of the blue opsin gene have evolved at the average rate of introns in primate genes.

Substitution rates were estimated for the coding and noncoding regions of the hepatitis delta virus (HDV). The estimated rates of synonymous substitution in HDV were lower than the rates of substitution at non-synonymous sites and in the noncoding region. HDV has lower synonymous substitution rates than the hepatitis C virus, though both are RNA viruses. The relatively low rate of synonymous substitution in HDV may be due to a strong preference of G and C nucleotides at third codon positions. Variation in substitution rate among HDV lineages may be correlated with the clinical development of the HDV-induced hepatitis. The phylogenetic tree inferred for 24 HDV strains reveals similarities between lineages isolated from the same geographic region.

To determine the evolutionary positions of the conifer genera Amentotaxus, Phyllocladus, and Nageia, we obtained 18S rRNA sequences from 11 new taxa representing the major living orders and families of gymnosperms. With the published Chlamydomonas as an outgroup, phylogenetic analyses of our new data and available sequences indicate that (1) the Gnetales form a monophyletic group, which is an outgroup to the conifers, (2) the conifers are monophyletic, (3) Taxaceae, Cephalotaxaceae, Cupressaceae, and Taxodiaceae form a monophyletic group, (4) Amentotaxus is closer to Torreya than to Cephalotaxus, suggesting that Amentotaxus is better to be classified as a member of Taxaceae, (5) Phyllocladus, Dacrycarpus, Podocarpus, and Nageia form a monophyletic group, and (6) Pinaceae is an out-group to the other families of conifers. Our finding that Phyllocladus is a sister group of the Podocarpaceae disagrees with the suggestion that the phylloclade of the genus is an ancient structure and that the genus is a terminal taxon within the Podocarpaceae. The genus Nageia is more closely related to Podocarpus than to Dacrycarpus and was derived from within the Podocarpaceae. In conclusion, our data indicate that in conifers, the uniovulate cone occurred independently in Taxacaeae and Cephalotaxaceae, and in Podocarpaceae after the three families separated from Pinaceae, and support the hypothesis that the uniovulate cone is derived from reduction of a multiovulate cone.

The size distributions of deletions, insertions, and indels (i.e., insertions or deletions) were studied, using 78 human processed pseudogenes and other published data sets. The following results were obtained: (1) Deletions occur more frequently than do insertions in sequence evolution; none of the pseudogenes studied shows significantly more insertions than deletions. (2) Empirically, the size distributions of deletions, insertions, and indels can be described well by a power law, i.e., fk = Ck-b, where fk is the frequency of deletion, insertion, or indel with gap length k, b is the power parameter, and C is the normalization factor. (3) The estimates of b for deletions and insertions from the same data set are approximately equal to each other, indicating that the size distributions for deletions and insertions are approximately identical. (4) The variation in the estimates of b among various data sets is small, indicating that the effect of local structure exists but only plays a secondary role in the size distribution of deletions and insertions. (5) The linear gap penalty, which is most commonly used in sequence alignment, is not supported by our analysis; rather, the power law for the size distribution of indels suggests that an appropriate gap penalty is wk = a + b ln k, where a is the gap creation cost and blnk is the gap extension cost. (6) The higher frequency of deletion over insertion suggests that the gap creation cost of insertion (ai) should be larger than that of deletion (ad); that is, ai - ad = ln R, where R is the frequency ratio of deletions to insertions.

Using sequence data from the last introns of ZFX and ZFY genes, we previously estimated the male-to-female ratio (alpha) of mutation rate to be close to 6 in higher primates and 1.8 in rodents. As the mutation rate may vary among different regions of the mammalian genome, it is interesting to see whether sequence data from other regions will give similar estimates. In this study, we have determined the partial genomic sequences of the ubiquitin-activating enzyme E1 genes (Ube 1x and Ube 1y for the X-linked and Y-linked homologues, respectively) of mice and rats and two mouse Ube 1y pseudogenes. From the intron sequences of the Ube 1 genes, we calculated the divergence of the Y-linked genes (Y = 0.161) and that of the X-linked genes (X = 0.107) between mouse and rat, and found the Y/X ratio to be 1.50. This ratio led to an estimate of alpha = 2.0 with a 95% confidence interval of (1.0, 3.9). Similar estimates of alpha were obtained if mouse Ube 1y pseudogenes were used instead of the mouse Ube 1y functional gene. These estimates are consistent with our previous estimate for rodents and suggest that the sex ratio of mutation rate in rodents is approximately only one-third of that in higher primates. Our estimate of the divergence time between Ube 1x and Ube 1y supports the view that the two genes separated before the eutherian radiation.

We have sequenced the entire exon (approximately 1.180 bp) encoding the zinc finger domain of the X-linked and Y-linked zinc finger genes (ZFX and ZFY, respectively) in the orangutan, the baboon, the squirrel monkey, and the rat; a total of 9,442 bp were sequenced. The ratio of the rates of synonymous substitution in the ZFY and ZFX genes is estimated to be 2.1 in primates. This is close to the ratio of 2.3 estimated from primate ZFY and ZFX intron sequences and supports the view that the male-to-female ratio of mutation rate in humans in considerably higher than 1 but not extremely large. The ratio of synonymous substitution rates in ZFY and ZFX is estimated to be 1.3 in the rat lineage but 4.2 in the mouse lineage. The former is close to the estimate (1.4) from introns. The much higher ratio in the mouse lineage (not statistically significant) might have arisen from relaxation of selective constraints. The synonymous divergence between mouse and rat ZFX is considerably lower than that between mouse and rat autosomal genes, agreeing with previous observations and providing some evidence for stronger selective constraints on synonymous changes in X-linked genes than in autosomal genes. At the protein level ZFX has been highly conserved in all placental mammals studied while ZFY has been well conserved in primates and foxes but has evolved rapidly in mice and rats, possibly due to relaxation of functional constraints as a result of the development of X-inactivation of ZFX in rodents. The long persistence of the ZFY-ZFX gene pair in mammals provides some insight into the process of degeneration of Y-linked genes.

Based on the biochemical kinetics of DNA replication and mutagenesis, including misincorporation and correction, a model has been developed for studying the relationships among the mutation rate (u), the G+C content of the sequence (f), and the G+C proportion in the nucleotide precursor pool (N). Also a measure for the next-nucleotide effect, called the maximum capacity of the next-nucleotide effect (MC), has been proposed. Under the normal physiological conditions of mammalian germ cells, our results indicate: (1) the equilibrium G+C content in a sequence is approximately equal to the G+C proportion in the nucleotide precursor pool, i.e., f approximately N, which is independent of the next-nucleotide effect; (2) an inverted-V-shaped distribution of mutation rates with respect to G+C contents is predicted, when the next-nucleotide effect is week, i.e., MC approximately 1; (3) the distribution becomes flatter (i.e., inverted-U-shaped) as MC increases, but the peak at 50% GC is still observed when MC < 2; and (4) the peak disappears when MC > 2.8, that is, when the next-nucleotide effect becomes strong. Our results suggest that changes in the relative concentrations of nucleotide precursors can cause variations among genes both in mutation rate and in G+C content and that compositional isochores (DNA segments with a homogeneous G+C content) can arise in a genome due to differences in replication times of DNA segments.

The evolutionary position of the yew family, Taxaceae, has been very controversial. Some plant taxonomists strongly advocate excluding Taxaceae from the conifer order and raising its taxonomic status to a new order or even class because of its absence of seed cones, contrary to the case in the majority of conifers. However, other authors believe that the Taxaceae are not fundamentally different from the rest of the conifers except in that they possess the most reduced solitary-ovule cones. To resolve the controversy, we have sequenced the 18S rRNA genes from representative gymnosperms: Taxus mairei (Taxaceae), Podocarpus nakaii (Podocarpaceae), Pinus luchuensis (Pinaceae), and Ginkgo biloba (Ginkgoales). Our phylogenetic analysis of the new sequence data with the published 18S rRNA sequence of Zamia pumila (a.cycad) as an outgroup strongly indicates that Taxus, Pinus, and Podocarpus form a monophyletic group with the exclusion of Ginkgo and that Taxus is more closely related to Pinus than to Podocarpus. Therefore, Taxaceae should be classified as a family of Coniferales. Our finding that Taxaceae, Pinaceae, and Podocarpaceae form a clade contradicts both the view that the uniovulate seed of Taxaceae is a primitive character and the view that the Taxaceae are descendants of the Podocarpaceae. Rather, the uniovulate seed of Taxaceae and that of some species of Podocarpus appear to have different origins, probably all reduced from multiovulate cones.

It is commonly believed that the rate of mutation is much higher in males than in females because the number of germ-cell divisions per generation is much larger in males than in females. However, the precise magnitude of the male-to-female mutation rate ratio (alpha m) remains unknown. Recently there have been efforts to estimate alpha m by using DNA sequence data from different species. We have studied the potential problems in such an approach. We found that the rate of synonymous substitution varies about fivefold among X-linked genes, as large as the variation among autosomal genes. This large variation makes the assumption of selective neutrality of synonymous changes dubious, so one should be cautious in using the synonymous rates in X-linked and autosomal genes to estimate alpha m. A similar difficulty was also observed in using nonhomologous intron sequences to estimate alpha m. Contrary to the expectation that X-linked sequences should evolve more slowly than autosomal sequences, the Alu repeat in the last intron of the X-linked zinc finger gene has evolved faster than the four autosomal Alu repeats used in this study. It appears that the best way to estimate alpha m is to use homologous sequences. However, such sequences may be involved in gene conversion events. In fact, we found evidence that the Y-linked and X-linked zinc finger genes have been involved in multiple conversion events during primate evolution. Thus, the possibility of gene conversion should be considered when using homologous sequences to estimate alpha m.

We have sequenced the mitochondrial cytochrome b gene from the guinea pig, the African porcupine, and a South American opossum. A phylogenetic analysis, which includes 22 eutherian and four other vertebrate cytochrome b sequences, indicates that the guinea pig and the porcupine constitute a natural clade (Hystricomorpha) that is not a sister group to the clade of mice and rats (Myomorpha). Therefore, the hypothesis that the Rodentia is paraphyletic receives additional support. The artiodactyls, the perissodactyls, and the cetaceans form a group that is separated from the primates and the rodents. The 26 sequences are used to study the structure/function relationships in cytochrome b, whose function is electron transport. Most of the amino acid residues involved in the two reaction centers are well conserved in evolution. The four histidines that are believed to ligate the two hemes are invariant among the 26 sequences, but their nearby residues are not well conserved in evolution. The eight transmembrane domains represent some of the most divergent regions in the cytochrome b sequence. The rate of nonsynonymous substitution is considerably faster in the human and elephant lineages than in other eutherian lineages; the faster rate might be due to coevolution between cytochrome b and cytochrome c.

Extensive DNA sequence data are used to compare the rates of nucleotide substitution in the mouse, rat, and hamster lineages. A relative rate test using hamster sequences as references shows that the rates of synonymous and nonsynonymous substitution in the mouse and rat lineages are nearly equal and a test using human sequences as references shows that the rates in the mouse, rat, and hamster lineages are also nearly equal. Under the assumptions that the guinea pig lineage and the myomorph (mouse, rat, and hamster) lineage diverged 70-100 million years (Myr) ago and that the rate of nucleotide substitution has been constant in all these lineages since their divergence, the date of the mouse-rat split is estimated to be between 20 and 29 Myr ago, which is considerably older than the date (approximately 12 Myr) suggested by available rodent fossils and considerably younger than the date (approximately 35 Myr) suggested by Wilson and colleagues. The murid-hamster split is estimated to be 1.6 times older than the mouse-rat split.

The statistical properties of sample estimation and bootstrap estimation of phylogenetic variability from a sample of nucleotide sequences were studied by considering model trees of three taxa with an outgroup. The cases of constant and varying rates of nucleotide substitution were compared. From sequences obtained by simulation, phylogenetic trees were constructed by using the maximum parsimony (MP) and neighbor-joining (NJ) methods. The effectiveness and consistency of the MP method were studied in terms of proportions of informative sites. The results of simulation showed that bootstrap estimation of the confidence level for an inferred phylogeny can be used even under unequal rates of evolution if the rate differences are not large so that the MP method is not misleading. The condition under which the MP method becomes misleading (inconsistent) is more stringent for slowly evolving sequences than for rapidly evolving ones, and it also depends on the length of the internal branch. If the rate differences are large so that the MP method becomes consistently misleading, then bootstrap estimation will reinforce an erroneous conclusion on topology. Similar conclusions apply to the NJ method with uncorrected distances. The NJ method with corrected distances performs poorly when the sequence length is short but can avoid the inconsistency problem if the sequence length is long and if the distances can be estimated accurately.

Apolipoprotein (apo) E, a major protein component of plasma lipoproteins, is a physiological ligand for the low density lipoprotein (LDL) receptor as well as for a specific apoE receptor; it is therefore an important modulator of lipoprotein metabolism. In this study we cloned and sequenced bovine apoE complementary DNA. Comparison of nucleotide substitution rates shows that apoE is less conservative than apoA-I and evolves about 30% faster than an average mammalian protein. Although apoE is not a conservative protein, several regions have been well conserved among all eight mammalian sequences now available. These include a 33-amino-acid block immediately upsteam from the third intron/exon junction and the LDL receptor binding region. We have also compared published apoC-I and apoC-II sequences. Both proteins are less conservative than apoE. In particular, apoC-I shows no well-conserved region except for a small region in the common 33-amino-acid block, suggesting that the function of apoC-I does not have stringent structural requirements. On the other hand, in apoC-II the region encoded by exon 4, which consists of the last 29 amino acids of the polypeptide, has been rather well conserved, probably because this region is important for the activation of lipoprotein lipase and chylomicron and very low density lipoprotein metabolism.

The population dynamics of nearly neutral mutations are studied using a single-site and a multisite model. In the latter model, the nucleotides in a sequence are completely linked and the selection schemes employed are additive, multiplicative, and additive with a threshold. Although the third selection scheme is very different from the first two, the three schemes produce identical results for a wide range of parameter values. Thus the present study provides a general theory for the population dynamics of nearly neutral mutations because the results can also be used to draw inferences about other selection schemes such as stabilizing selection and synergistic selection. It is shown that the number of slightly deleterious mutations accumulated in a sequence can be considerably larger under the multisite model than under the single-site model, particularly if the sequence is long or if the mutation rate per site is high. The results show that even a very slight selective difference between synonymous codons can produce a strong bias in codon usage. Three alternative explanations for the strong bias in codon usage in bacterial and yeast genes are considered. The implications of the present results for molecular evolution are discussed.

A simulation study has been conducted to check the accuracy of Nei and Li's (1979) formulas for the mean and variance of the proportion (S) of identical restriction sites between two DNA sequences and for estimating the mean and variance of the number (delta) of base substitutions per nucleotide site between two DNA sequences. The results show that these formulas are quite accurate as long as the probability of S becoming zero is negligibly. In addition to the simulation, approximate formulas have also been obtained for the probability for S to become zero at time t and for the contribution to S due to parallel mutation.

Under the assumption of unequal rates of nucleotide substitution among the three positions of codons, mathematical formulas are derived for the probability that a restriction site observed at time 0 in a DNA sequence will be present at time t, the probability that a new restriction site will emerge at a particular place in two descendant sequences, and the proportion of identical restriction sites between two such sequences. All three quantities are shown to be larger for the case of unequal rates than for the case of equal rates. As a consequence, estimates of nucleotide divergence (2 lambda t) between two sequences based on restriction-site data tend to be lower than the actual values if the assumption of equal rates is made but actually does not hold. However, the degree of underestimation is slight if 2 lambda t is 0.10 or smaller. The underestimation can be serious when 2 lambda t becomes 0.20 or larger, particularly if the substitution rates at the three codon positions are very different or if transitional substitution occurs more frequently than transversional substitution. The underestimation is more serious for four-base enzymes than for six-base enzymes.

The relative rates of point nucleotide substitution and accumulation of gap events (deletions and insertions) were calculated for 22 human and 30 rodent processed pseudogenes. Deletion events not only outnumbered insertions (the ratio being 7:1 and 3:1 for human and rodent pseudogenes, respectively), but also the total length of deletions was greater than that of insertions. Compared with their functional homologs, human processed pseudogenes were found to be shorter by about 1.2%, and rodent pseudogenes by about 2.3%. DNA loss from processed pseudogenes through deletion is estimated to be at least seven times faster in rodents than in humans. In comparison with the rate of point substitutions, the abridgment of pseudogenes during evolutionary times is a slow process that probably does not retard the rate of growth of the genome due to the proliferation of processed pseudogenes.

Observed patterns of synonymous codon usage are explained in terms of the joint effects of mutation, selection, and random drift. Examination of the codon usage in 165 Escherichia coli genes reveals a consistent trend of increasing bias with increasing gene expression level. Selection on codon usage appears to be unidirectional, so that the pattern seen in lowly expressed genes is best explained in terms of an absence of strong selection. A measure of directional synonymous-codon usage bias, the Codon Adaptation Index, has been developed. In enterobacteria, rates of synonymous substitution are seen to vary greatly among genes, and genes with a high codon bias evolve more slowly. A theoretical study shows that the patterns of extreme codon bias observed for some E. coli (and yeast) genes can be generated by rather small selective differences. The relative plausibilities of various theoretical models for explaining nonrandom codon usage are discussed.

A statistical analysis of extensive DNA sequence data from primates, rodents, and artiodactyls clearly indicates that no global molecular clock exists in mammals. Rates of nucleotide substitution in rodents are estimated to be four to eight times higher than those in higher primates and two to four times higher than those in artiodactyls. There is strong evidence for lower substitution rates in apes and humans than in monkeys, supporting the hominoid slowdown hypothesis. There is also evidence for lower rates in humans than in apes, suggesting a further rate slowdown in the human lineage after the separation of humans from apes. By contrast, substitution rates are nearly equal in mouse and rat. These results suggest that differences in generation time or, more precisely, in the number of germline DNA replications per year are the primary cause of rate differences in mammals. Further, these differences are more in line with the neutral mutation hypothesis than if the rates are the same for short- and long-living mammals.

A statistical analysis of the nucleotide sequences of cytochrome c genes from four species of animals and two of yeast and of cytochrome c pseudogenes from rat, mouse, and human was conducted. It was estimated that animals and yeast diverged 1.2 billion years ago, that the two duplicated genes DC3 and DC4 in Drosophila diverged 520 million years ago, and that the two duplicated genes Iso-1 and Iso-2 in the yeast Saccharomyces cerevisiae diverged 200 million years ago. DC3 is expressed at a low level and has evolved 3 times faster than DC4. This observation supports the neutralist view that relaxation of functional constraints is a more likely cause of accelerated evolution following gene duplication than is advantageous mutation. All the rodent pseudogenes examined appear to be processed pseudogenes derived directly from the functional genes, and most of them apparently arose after the mouse-rat split. No event of gene conversion could be detected between any pair of the rodent pseudogenes. Our analysis suggests that the human cytochrome c gene has evolved at a rate comparable to the average rate for pseudogenes, whereas some human cytochrome c pseudogenes have evolved at an exceptionally low rate.

The effect of changes in population size on the correlation between mutation rate and heterozygosity was studied by using two models: sudden change in population size and gradual change. It was shown that the results for the two models are close to each other, unless the rate of change for the latter is exceedingly slow. Thus, in many cases, the former model, which is much simpler than the latter, can be used to treat the present problem. Numerical computations showed that the correlation in a population that is expanding or has expanded in the recent past is stronger while the correlation in a population that is decreasing or has experienced a population reduction or bottleneck in the recent past is weaker than that for an equilibrium population with the same mean heterozygosity. However, regardless of whether the population is at equilibrium or not, the proportion of variation in heterozygosity that is attributable to variation in molecular weight over loci is rather small if the mean heterozygosity of the population is low, say of the order 0.05 or smaller.

Analysis of the rate of nucleotide substitution at silent sites in Drosophila genes reveals three main points. First, the silent rate varies (by a factor of two) among nuclear genes; it is inversely related to the degree of codon usage bias, and so selection among synonymous codons appears to constrain the rate of silent substitution in some genes. Second, mitochondrial genes may have evolved only as fast as nuclear genes with weak codon usage bias (and two times faster than nuclear genes with high codon usage bias); this is quite different from the situation in mammals where mitochondrial genes evolve approximately 5-10 times faster than nuclear genes. Third, the absolute rate of substitution at silent sites in nuclear genes in Drosophila is about three times higher than the average silent rate in mammals.

We have obtained a revised estimate of the pattern of point mutation by considering more pseudogene sequences. Compared with our previous estimate, it agrees better with expectations based on the double-strand structure of DNA. The revised pattern, like the previous one, indicates that mutation occurs nonrandomly among the four nucleotides. In particular, the proportion of transitional mutations (59%) is almost twice as high as the value (33%) expected under random mutation. The same high proportion of transitions is observed in synonymous substitutions in genes. The proportion of transitional changes observed among electrophoretic variants of human hemoglobin is about the same as that predicted by the revised pattern of mutation. We also show that nonrandom mutation increases, by about 15%, the proportion of synonymous mutations due to single-nucleotide changes in the codon table, and increases, from 10% to 50%, the rate of synonymous mutation in the seven genes studied. However, nonrandom mutation reduces (by about 10%) the proportion of polar changes among nonsynonymous mutations in a gene. As far as single-nucleotide changes (in the codon table) are concerned, nonrandom mutation only slightly favors relatively conservative amino acid interchanges, and has virtually no effect on the proportions of radical changes and nonsense mutations.

In at least two instances involving serine proteinase inhibitors it has been shown that functionally important sites evolve faster and exhibit more interspecific variability than functionally neutral sites. Because these phenomena are difficult to reconcile with the neutral theory of molecular evolution, it has been suggested that the accelerated rate of amino acid substitution at the reactive sites is brought about by positive Darwinian selection. We show that differences in the amino acid composition in the different regions of proteinase inhibitors can account for the differences in the rates of amino acid substitution. By using an index of protein mutability [D. Graur (1985) J Mol Evol 22:53-62], we show that the amino acid composition of the reactive center in the ovomucoids and Spi-2 gene products is such that, regardless of function, they are expected to evolve more rapidly than any other polypeptide for which the rate of substitution is known. In addition, the reactive region in the Spi-2 proteins is shown to be free of compositional constraint. Positive Darwinian selection need not be invoked at the present time in these cases.

Ubiquitin is remarkable for its ubiquitous distribution and its extreme protein sequence conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences of several ubiquitin repeats from each of humans, chicken, Xenopus, Drosophila, barley, and yeast have recently been determined. By analysis of these data we show that ubiquitin is evolving more slowly than any other known protein, and that this (together with its gene organization) contributes to an ideal situation for the occurrence of concerted evolution of tandem repeats. By contrast, there is little evidence of between-cluster concerted evolution. We deduce that in ubiquitin genes, concerted evolution involves both unequal crossover and gene conversion, and that the average time since two repeated units within the polyubiquitin locus most recently shared a common ancestor is approximately 38 million years (Myr) in mammals, but perhaps only 11 Myr in Drosophila. The extreme conservatism of ubiquitin evolution also allows the inference that certain synonymous serine codons differing at the first two positions were probably mutated at single steps.