核酸起始序列、终止序列和插入序列的统计分析

对八百个核酸序列进行了统计分析,讨论了起始序列、终止序列、插入序列、间隙序列以及序列的氨基酸编码区的各项信息指标,进行了相互比较,并讨论了这些指标和分子进化的统计相关性.     

Molecular and morphometric evidence for separate species of Uncinaria (Nematoda: Ancylostomatidae) in California sea lions and northern fur seals: hypothesis testing supplants verification

California sea lions (Zalophus californianus) and northern fur seals (Callorhinus ursinus) are each believed to host distinct hookworm species (Uncinaria spp.). However, a recent morphometric analysis suggested that a single species parasitizes multiple pinniped hosts, and that the observed differences are host-induced. To explore the systematics of these hookworms and test these competing hypotheses, we obtained nucleotide sequences of nuclear ribosomal DNA (D2/D3 28S, D18/D19 28S, and internal transcribed spacer [ITS] regions) from 20 individual hookworms parasitizing California sea lion and northern fur seal pups where their breeding grounds are sympatric. Five individuals from an allopatric population of California sea lions were also sampled for ITS-1 and D18/D19 28S sequences. The 28S D2/D3 sequences showed no diagnostic differences among hookworms sampled from individual sea lions and fur seals, whereas the 28S D18/D19 sequences had one derived (apomorphic) character demarcating hookworms from northern fur seals. ITS sequences were variable for 7 characters, with 4 derived (apomorphic) states in ITS-1 demarcating hookworms from California sea lions. Multivariate analysis of morphometric data also revealed significant differences between nematodes representing these 2 host-associated lineages. These results indicate that these hookworms represent 2 species that are not distributed indiscriminately between these host species, but instead exhibit host fidelity, evolving independently with each respective host species. This evolutionary approach to analyzing sequence data for species delimitation is contrasted with similarity-based methods that have been applied to numerous diagnostic studies of nematode parasites.     

Diversity of staphylocoagulase and identification of novel variants of staphylocoagulase gene in Staphylococcus aureus

Staphylocoagulase (SC) is a major phenotypic determinant of Staphylococcus aureus. Serotype of SC (coagulase type) is used as an epidemiological marker and 10 types (I-X) have been discriminated so far. To clarify genetic diversity of SC within a single and among different serotype(s), we determined approximately 1500 bp-nucleotide sequences of SC gene encoding D1, D2, and central regions (N-terminal half and central regions of SC; SC(NC)) for a total of 33 S. aureus strains comprising two to three strains from individual coagulase types (I-VIII, X) and 10 strains which were not determined as previously known SC serotypes (ND-strains). Amino acid sequence identities of SC(NC) among strains with a single coagulase type of II, III, IV, V, VI and X were extremely high (more than 99%), whereas lower identity (56-87%) was observed among different types. In contrast, within a single coagulase type of I, VII, or VIII, sequence divergence was found (lowest identity; 82%). SC(NC) sequences from the ND-strains were discriminated into two genetic groups with an identity of 71% to each other (tentatively assigned to genotypes [XI] and [XII]), and exhibited less than 86% sequence identities to those of most known coagulase types. All the types [XI] and [XII] strains were methicillin susceptible and belonged to different sequence types from those of coagulase types I-X strains reported so far by multilocus sequence typing. These findings indicated genetic heterogeneity of SC in coagulase types I, VII, and VIII strains, and the presence of two novel SC genotypes related to antigenicity of SC serotypes.     

Insulin Binding to Human Astrocytoma Cells and Its Effect on Uridine Incorporation into Nucleic Acid

Binding of [125I]monoiodoinsulin to human astrocytoma cells (U-373 MG) was time dependent, reaching equilibrium after 1 h at 22 degrees C with equilibrium binding corresponding to 2.2 fmol/mg protein: this represents approximately 2,000 occupied binding sites per cell. The t1/2 of 125I-insulin dissociation at 22 degrees C was 10 min; the dissociation rate constant of 1.1 X 10(-2) s-1 was unaffected by a high concentration of unlabeled insulin (16.7 microM). Porcine insulin competed for specific 125I-insulin binding in a dose-dependent manner and Scatchard analysis suggested multiple affinity binding sites (higher affinity Ka = 4.4 X 10(8) M-1 and lower affinity Ka = 7.4 X 10(6) M-1). Glucagon and somatostatin did not compete for specific insulin binding. Incubation of cells with insulin (0.5 microM) for 2 h at 37 degrees C increased [2-14C]uridine incorporation into nucleic acid by 62 +/- 2% (n = 3) above basal. Cyclic AMP, in the absence of insulin, also stimulated nucleoside incorporation into nucleic acid [65 +/- 1% (n = 3)] above basal. Preincubation with cyclic AMP followed by insulin had an additive effect on nucleoside incorporation [160 +/- 4% (n = 3) above basal]. Dipyridamole (50 microM), a nucleoside transport inhibitor, blocked both basal and stimulated uridine incorporation. These studies confirm that human astrocytoma cells possess specific insulin receptors with a demonstrable effect of ligand binding on uridine incorporation into nucleic acid.     

Three-dimensional cluster analysis identifies interfaces and functional residue clusters in proteins

Three-dimensional cluster analysis offers a method for the prediction of functional residue clusters in proteins. This method requires a representative structure and a multiple sequence alignment as input data. Individual residues are represented in terms of regional alignments that reflect both their structural environment and their evolutionary variation, as defined by the alignment of homologous sequences. From the overall (global) and the residue-specific (regional) alignments, we calculate the global and regional similarity matrices, containing scores for all pairwise sequence comparisons in the respective alignments. Comparing the matrices yields two scores for each residue. The regional conservation score (C(R)(x)) defines the conservation of each residue x and its neighbors in 3D space relative to the protein as a whole. The similarity deviation score (S(x)) detects residue clusters with sequence similarities that deviate from the similarities suggested by the full-length sequences. We evaluated 3D cluster analysis on a set of 35 families of proteins with available cocrystal structures, showing small ligand interfaces, nucleic acid interfaces and two types of protein-protein interfaces (transient and stable). We present two examples in detail: fructose-1,6-bisphosphate aldolase and the mitogen-activated protein kinase ERK2. We found that the regional conservation score (C(R)(x)) identifies functional residue clusters better than a scoring scheme that does not take 3D information into account. C(R)(x) is particularly useful for the prediction of poorly conserved, transient protein-protein interfaces. Many of the proteins studied contained residue clusters with elevated similarity deviation scores. These residue clusters correlate with specificity-conferring regions: 3D cluster analysis therefore represents an easily applied method for the prediction of functionally relevant spatial clusters of residues in proteins.     

Titration of integrated simian virus 40 DNA sequences, using highly radioactive, single-stranded DNA probes.

Nick-translated simian virus 40 (SV40) [32P]DNA fragments (greater than 2 X 10(8) cpm/micrograms) were resolved into early- and late-strand nucleic acid sequences by hybridization with asymmetric SV40 complementary RNA. Both single-stranded DNA fractions contained less than 0.5% self-complementary sequences; both included [32P]-DNA sequences that derived from all regions of the SV40 genome. In contrast to asymmetric SV40 complementary RNA, both single-stranded [32P]DNAs annealed to viral [3H]DNA at a rate characteristic of SV40 DNA reassociation. Kinetics of reassociation between the single-stranded [32P]DNAs indicated that the two fractions contain greater than 90% of the total nucleotide sequences comprising the SV40 genome. These preparations were used as hybridization probes to detect small amounts of viral DNA integrated into the chromosomes of Chinese hamster cells transformed by SV40. Under the conditions used for hybridization titrations in solution (i.e., 10- to 50-fold excess of radioactive probe), as little as 1 pg of integrated SV40 DNA sequence was assayed quantitatively. Among the transformed cells analyzed, three clones contained approximately one viral genome equivalent of SV40 DNA per diploid cell DNA complement; three other clones contained between 1.2 and 1.6 viral genome equivalents of SV40 DNA; and one clone contained somewhat more than two viral genome equivalents of SV40 DNA. Preliminary restriction endonuclease maps of the integrated SV40 DNAs indicated that four clones contained viral DNA sequences located at a single, clone-specific chromosomal site. In three clones, the SV40 DNA sequences were located at two distinct chromosomal sites.     

New statistical approach to discriminate between protein coding and non-coding regions in DNA sequences and its evaluation

We propose a new approach to study protein coding and non-coding regions in DNA sequences, by making use of two complementary statistical methods. The principal component analysis (PCA) is a graphical method to represent DNA sequences which are characterized by some quantitative parameters: it is a help to the intuition. The discriminating analysis (DA) is a quantitative method which permits to classify the DNA sequences. It leads to an evaluation of the first method and to a decision. The value of this approach has been confirmed since we also have found some results which had been described recently in the literature. Furthermore, this general methodology has permitted us to show the existence of parameters which identify the nucleic acid sequence functional domains, without having to make use of the properties of the genetic code.     

Conformational transitions in polyelectrolyte adsorption

The mole fraction X of nucleic acid chromophores, adsorbed on cationic and anionic exchange chromatography paper disks, was calculated [O.L. H?rer and A. Grigorescu, Rev. Roum. Med. Virol., 43 (1992) 33] in terms of the standard Debye screening length D (in the 0.01 to 1.2 M ionic strength range), using data from MLSEA (multiple light scattering enhanced absorption) measurements. In contrast to the cases of mononucleotides or double-stranded and relatively rigid DNA structures, highly polymerized and single-stranded RNA shows sigmoidal-shaped X curves, suggesting that polymer adsorption may be regarded as a conformational transition from a three-dimensional statistical coil to a two-dimensional one.     

Host-specific modulation of the selective constraints driving human immunodeficiency virus type 1 env gene evolution

To address the evolution of human immunodeficiency virus type 1 (HIV-1) within a single host, we analyzed the HIV-1 C2-V5 env regions of both cell-free genomic-RNA- and proviral-DNA-derived clones. Sequential samples were collected over a period of 3 years from six untreated subjects (three typical progressors [TPs] and three slow progressors [SPs], all with a comparable length of infection except one. The evolutionary analysis of the C2-V5 env sequences performed on 506 molecular clones (253 RNA- and 253 DNA-derived sequences) highlighted a series of differences between TPs and SPs. In particular, (i) clonal sequences from SPs (DNA and RNA) showed lower nucleotide similarity than those from TPs (P = 0. 0001), (ii) DNA clones from SPs showed higher intra- and intersample nucleotide divergence than those from TPs (P < 0.05), (iii) higher host-selective pressure was generally detectable in SPs (DNA and RNA sequences), and (iv) the increase in the genetic distance of DNA and RNA sequences over time was paralleled by an increase in both synonymous (Ks) and nonsynonymous (Ka) substitutions in TPs but only in nonsynonymous substitutions in SPs. Several individual peculiarities of the HIV-1 evolutionary dynamics emerged when the V3, V4, and V5 env regions of both TPs and SPs were evaluated separately. These peculiarities, probably reflecting host-specific features of selective constraints and their continuous modulation, are documented by the dynamics of Ka/Ks ratios of hypervariable env domains.     

Metabolic fates of 2′-deoxy-5-[18F]fluorouridine in tumor-bearing mice and human plasma

The metabolic fate of 2′-deoxy-5-[¹⁸F]fluorouridine ([¹⁸F]FdUrd), a useful positron emission tomography (PET) tracer of nucleic acid metabolism in tumors, was investigated in mice and humans. A rapid increase in labeled catabolites was found in mouse and human plasma. In mouse FM3A mammary carcinoma, the corresponding catabolites were also detected in addition to metabolites which were activated by the nucleic acid metabolism. From a biodistribution study of β-[³H]alanine, α-[¹⁸F]fluoro-β-alanine, a major catabolite, was assumed to be taken up twice as much by tumor than by the brain. Nucleic acid metabolism in brain tumors by [¹⁸F]FdUrd-PET may be assessed using normal brain regions as a reference.     

Evolutionary breakpoints are co-localized with fragile sites and intrachromosomal telomeric sequences in primates

The concentration of evolutionary breakpoints in primate karyotypes in some particular regions or chromosome bands suggests that these chromosome regions are more prone to breakage. This is the first extensive comparative study which investigates a possible relationship of two genetic markers (intrachromosomal telomeric sequences [TTAGGG]n, [ITSs] and fragile sites [FSs]), which are implicated in the evolutionary process as well as in chromosome rearrangements. For this purpose, we have analyzed: (a) the cytogenetic expression of aphidicolin-induced FSs in Cebus apella and Cebus nigrivittatus (F. Cebidae, Platyrrhini) and Mandrillus sphinx (F. Cercopithecidae, Catarrhini), and (b) the intrachromosomal position of telomeric-like sequences by FISH with a synthetic (TTAGGG)n probe in C. apella chromosomes. The multinomial FSM statistical model allowed us to determinate 53 FSs in C. apella, 16 FSs in C. nigrivittatus and 50 FSs in M. sphinx. As expected, all telomeres hybridized with the probe, and 55 intrachromosomal loci were also detected in the Cebus apella karyotype. The chi(2) test indicates that the coincidence of the location of Cebus and Mandrillus FSs with the location of human FSs is significant (P < 0.005). Based on a comparative cytogenetic study among different primate species we have identified (or described) the chromosome bands in the karyotypes of Papionini and Cebus species implicated in evolutionary reorganizations. More than 80% of these evolutionary breakpoints are located in chromosome bands that express FSs and/or contain ITSs.     

羊轮状病毒NT株VP1基因的测序和分子进化分析

摘要：【目的】为了研究羊轮状病毒NT株VP1基因的遗传进化规律,【方法】根据GenBank中相关VP1基因的保守序列,设计合成引物,扩增NT株VP1基因并进行克隆测序和序列分析。【结果】 氨基酸序列比较表明NT株与其他毒株VP1基因的相似性为77.3%～98.4%,且氨基酸突变多发生在VP1蛋白的非功能区。VP1蛋白进化树表明NT株与牛轮状病毒处于同一进化分支,有较近的亲缘关系。结合26株具有代表性的轮状病毒,计算毒株间VP1基因的核苷酸和氨基酸进化距离,并对核苷酸的同义突变率（dS）和非同义突变率（dN）进行研究,发现dN/dS的比值小于1,说明同义替代是VP1基因在进化过程中的主要变异。【结论】本文首次对羊轮状病毒NT株进行了VP1基因的测序,并对VP1基因的进化距离和进化规律进行深入探讨。     

Analysis of bacteriophage phi X174 gene A protein-mediated termination and reinitiation of phi X DNA synthesis. II. Structural characterization of the covalent phi X A protein-DNA complex

In the preceeding paper (Brown, D. R., Roth, M. J., Reinberg, D., and Hurwitz, J. (1984) J. Biol. Chem. 259, 10545-10555), it was shown that following bacteriophage phi X174 (phi X) DNA synthesis in vitro using purified proteins, the phi X A protein could be detected covalently linked to nascent 32P-labeled DNA. This phi X A protein-[32P]DNA complex was the product of the reinitiation reaction. The phi X A protein-[32P]DNA complex could be trapped as a protein-32P-oligonucleotide complex by the inclusion of ddGTP in reaction mixtures. In this report, the structure of the phi X A protein-32P-oligonucleotide complex has been analyzed. The DNA sequence of the oligonucleotide bound to the phi X A protein has been determined and shown to be homologous to the phi X (+) strand sequence immediately adjacent (3') to the replication origin. The phi X A protein was directly linked to the 5' position of a dAMP residue of the oligonucleotide; this residue corresponded to position 4306 of the phi X DNA sequence. The phi X A protein-32P-oligonucleotide complex was exhaustively digested with either trypsin or proteinase K and the 32P-labeled proteolytic fragments were analyzed. Each protease yielded two different 32P-labeled peptides in approximately equimolar ratios. The two 32P-labeled peptides formed after digestion with trypsin (designated T1 and T2) and with proteinase K (designated PK1 and PK2) were isolated and characterized. Digestion of peptide T1 with proteinase K yielded a product which co-migrated with peptide PK2. In contrast, peptide T2 was unaffected by digestion with proteinase K. These results suggest that the phi X A protein contains two active sites that are each capable of binding covalently to DNA. The peptide-mononucleotide complexes T1-[32P]pdA and T2-[32P]pdA were isolated and subjected to acid hydrolysis in 6.0 N HCl. In each case, the major 32P-labeled products were identified as [32P] phosphotyrosine and [32P]Pi. This indicates that each active site of the phi X A protein participates in a phosphodiester linkage between a tyrosyl moiety of the protein and the 5' position of dAMP.     

Partitioning of rearranged Ig genes by mutation analysis demonstrates D-D fusion and V gene replacement in the expressed human repertoire

The accurate partitioning of Ig H chain V(H)DJ(H) junctions and L chain V(L)J(L) junctions is problematic. We have developed a statistical approach for the partitioning of such sequences, by analyzing the distribution of point mutations between a determined V gene segment and putative Ig regions. The establishment of objective criteria for the partitioning of sequences between V(H), D, and J(H) gene segments has allowed us to more carefully analyze intervening putative nontemplated (N) nucleotides. An analysis of 225 IgM H chain sequences, with five or fewer V mutations, led to the alignment of 199 sequences. Only 5.0% of sequences lacked N nucleotides at the V(H)D junction (N1), and 10.6% at the DJ(H) junction (N2). Long N regions (>9 nt) were seen in 20.6% of N1 regions and 17.1% of N2 regions. Using a statistical analysis based upon known features of N addition, and mutation analysis, two of these N regions aligned with D gene segments, and a third aligned with an inverted D gene segment. Nine additional sequences included possible alignments with a second D segment. Four of the remaining 40 long N1 regions included 5' sequences having six or more matches to V gene end motifs, which may be the result of V gene replacement. Such sequences were not seen in long N2 regions. The long N regions frequently seen in the expressed repertoire of human Ig gene rearrangements can therefore only partly be explained by V gene replacement and D-D fusion.     

2009−2015年中国甲型H1N1流感病毒血凝素基因进化分析

【目的】为了解中国地区2009?2015年甲型H1N1流感病毒流行态势,分析血凝素(Hemagglutinin,HA)基因的变异情况及其遗传进化特征。【方法】汇集国家流感中心2009?2015年流感周报的流感流行数据,分析甲型H1N1流感的流行病学特征;从全球共享禽流感数据倡议组织数据库及美国国家生物技术中心数据库下载甲型H1N1流感病毒HA基因序列,采用生物学软件进行系统进化和遗传特性的分析。【结果】2009?2015年全国共发生4次甲型H1N1流感的流行高峰。2009?2015年毒株与参考毒株A/California/07/2009(H1N1)的HA基因同源性逐年降低。遗传进化分析显示同一年份的毒株在系统进化树上基本呈现集中分布,2011年的毒株独立形成2个分支。分子特征表现为HA基因的4个抗原决定簇氨基酸位点均有变异,其中Ca区的203位、Sa区的163位和Sb区的185位氨基酸位点逐渐替换为新的氨基酸。除2010年与2012年,其他年份的毒株通过不同模型均得到正向压力选择HA氨基酸位点240。【结论】甲型H1N1流感在中国地区成为主要流行的亚型之一。HA基因与其编码的氨基酸逐年变异,未来进一步的流感监测能力还需加强。     

核酸序列重复片段的统计分析

对核酸序列编码区的重复片段进行了统计分析。发现不同物种的序列具有近似相等的约化重复长度和次数,研究了重复区内碱基分布的特点,并讨论了此结果可能的生物学涵义。     

HPLC analysis of benzo[a]pyrene-albumin adducts in benzo[a]pyrene exposed rats. Detection of cis-tetrols arising from hydrolysis of adducts of anti- and syn-BPDE-III with proteins

Quantitation of protein-benzo[a]pyrene adducts represent a more sensitive analysis method than quantitation of benzo[a]pyrene-DNA adducts. By accurate analysis of benzo[a]pyrene-protein adducts several different molecular adduct forms can be studied. Male Wistar rats were injected i.p. with benzo[a]pyrene, and serum albumin was isolated and subjected to acid hydrolysis at 90 degrees C for 3 h. The hydrolysate was analyzed by HPLC with fluorescence detection. The HPLC profiles obtained after albumin hydrolysis from benzo[a]pyrene exposed animals were compared to similar HPLC profiles from in vitro adducted bovine serum albumin (BSA) and direct hydrolysis of both r-10,t-9-dihydrodiol-c-7,8-oxy-7,8,9,10-tetrahydrobenzo[a]pyrene (syn-BPDE-III) and r-10,t-9-t-dihydrodiol-t-7,8-oxy-7,8,9,10-tetrahydrobenzo[a]pyrene (anti-BPDE-III). After acid hydrolysis of albumin from benzo[a]pyrene exposed rats, 6 fluorescent peaks were separated. Four of the peaks were isomers of benzo[a]pyrene-tetrahydrotetrols, (+/-)-benzo[a]pyrene-r-7,t-8,9,10-tetrahydrotetrol, (+/-)-benzo[a]pyrene-r-7,t-8,9,c-10-tetrahydrotetrol, (+/-)-benzo[a]pyrene-r-7,t-8,c-9,t-10-tetrahydrotetrol and (+/-)-benzo[a]pyrene-r-7,t-8,c-9,10-tetrahydrotetrol. In addition we found two fluorescent peaks, named X1 and X2 with retention times similar to the benzo[a]pyrene-tetrols. The unknown fluorescent peaks reacted similar to the four known tetrols in both dose response experiments and time course experiments. Fluorescent material with retention times equal to X1 and X2 were found after acid hydrolysis of syn-BPDE-III and anti-BPDE-III in acid and in hydrolysates from BSA treated in vitro with syn-BPDE-III and anti-BPDE-III. The ratio X1/X2 was relatively constant indicating epimerization equilibrium between these to species. Synchronous fluorescence analysis of fractions containing X1 or X2 from both in vivo and in vitro experiments showed fluorescence spectra characteristic of benzo[a]pyrene tetrols using a wavelength difference of 34 nm.