Evidence for effect of random genetic drift on G+C content after lateral transfer of fucose pathway genes to Escherichia coli K-12

The cps cluster of Escherichia coli K-12 comprises genes involved in synthesis of capsular polysaccharide colanic acid. Part of the E. coli K-12 cps region has been cloned and sequenced and compared to its Salmonella enterica LT2 counterpart. The cps genes from the two organisms are homologous; in the case of the LT2 genes, with G+C content of 0.61 and codons characteristic of high G+C species, it seems clear that they have been acquired relatively recently by lateral transfer from a high G+C species. The K-12 form of these cps genes is closely related to those of LT2 so must derive from the same high G+C species, but it appears to have transferred much earlier such that random genetic drift has brought P3 (the corrected G+C content of codon base 3) down from 0.77 to 0.64, more than halfway to the E. coli average of 0.57. We estimate, using an equation developed by Sueoka, that the lateral transfer to E. coli took place approximately 45 million years ago. This is the first report we are aware of demonstrating the expected adjustment of P3 after lateral transfer between species with different G+C content DNA.      

Molecular evolution of centromere-associated nucleotide sequences in two species of canids

The major centromeric satellite nt sequences present in the domestic dog (Canis familiaris) and in the grey fox (Urocyon cineroargenteus) have been examined. The dog satellite monomer is 737 bp long and contains 51% G + C; the grey fox satellite monomer is 880 b long and contains 54% G + C. The two satellites share three regions of 78, 92 and 314 bp with 70-80% sequence similarity. Sequence data from 16 monomers of dog satellite and 19 monomers of grey fox satellite demonstrate that the substitution spectra are different in the two canid species. For example, substitutions involving G or C residues are much more common in the grey fox satellite than in the domestic dog satellite despite their similar G + C contents.     

The mitochondrial genome of wild-type yeast cells. IV. Genes and spacers

The organization of the mitochondrial genome of wild-type Saccharomyces cerevisiae cells has been investigated further, by degrading mitochondrial DNA with micrococcal nuclease. Under the conditions used, this enzyme very strongly degrades the A + T-rich stretches (spacers) whereas it only inflicts a limited number of breaks into the G + C-rich stretches (genes). The macromolecular fragments derived from the “genes” have been separated from the oligonucleotides originating from the “spacers” by gel filtration, and both sorts of products have been investigated. It has been shown (a) that the spacers are very homogeneous in base composition and have a G + C content lower than 5% (mitochondrial DNA has a G + C content of 18%); (b) that the genes are very heterogeneous in base composition, the G + C content ranging from about 25% to 50%, when the average size of the fragments is 1·2 × 10⁵; smaller fragments, molecular weight 4 × 10⁴, having a G + C level as high as 65%, have been isolated in a yield of 10%; the average G + C content of genes is about 32%; (c) that genes and spacers are present in about equal amounts in the mitochondrial genome and that they have comparable average sizes.     

Isochore structures in the mouse genome

The distribution of the G+C content in the mouse genome has been studied using a windowless technique. We have found that: (i). Abrupt variations of the G+C content from a GC-rich region to a GC-poor region, and vice versa, occur frequently at some sites along the sequence of the mouse genome. (ii). Long domains with relatively homogeneous G+C content (isochores) exist, which usually have sharp boundaries. Consequently, 28 isochores longer than 1 Mb have been identified in the mouse genome. A homogeneity index was used to quantify the variations of the G+C content within isochores. The precise boundaries, sizes, and G+C contents of these isochores have been determined. The windowless technique for the G+C content computation was also used to analyze the DNA sequence containing the mouse MHC region, which has a GC-poor isochore. This isochore is located at the central part of the sequence with boundaries at 468459 and 812716 bp, where the sequence is extended from the centromeric end to the telomeric end. In addition, the analysis of a segment of the rat genome shows that the rat genome also has clear isochore structures.     

Association of interleukin 23 receptor gene polymorphisms (rs10489629, rs7517847) with rheumatoid arthritis in European population: a meta-analysis

The interleukin 23 receptor (IL-23R) polymorphisms have been already discussed in rheumatoid arthritis (RA) repeatedly, but the results are conflict. The purpose of this meta-analysis was to assess whether IL-23R gene polymorphisms are associated with RA. We retrieved the available data from Pubmed, Medline, CNKI and CBM. Our study evaluated the effects of two polymorphisms (rs10489629, rs7517847) in European population. Pooling all the subjects, we found significant associations between the two polymorphisms and RA. For rs10489629, the pooled ORs (95 % CI) of C versus T, C/C+C/T versus T/T and C/C versus C/T+T/T were 1.092 (1.038-1.149), 1.146 (1.059-1.240) and 1.099 (1.008-1.199), respectively. For rs7517847, the combined ORs (95 % CI) of G versus T, G/G+G/T versus T/T and G/G versus G/T+T/T were 1.121 (1.063-1.183), 1.184 (1.092-1.283) and 1.133 (1.030-1.246), respectively. In conclusion, this meta-analysis demonstrates that the polymorphisms rs10489629 and rs7517847 of the IL-23R gene may be considered as risk factors for developing RA in European population.     

A novel method to calculate the G+C content of genomic DNA sequences.

The base composition of a DNA fragment or genome is usually measured by the proportion of A+T or G+C in the sequence. The G+C content along genomic sequences is usually calculated using an overlapping or non-overlapping sliding window method. The result and accuracy of such an approach depends on the size of the window and the moving distance adopted. In this paper, a novel windowless technique to calculate the G+C content of genomic sequences is proposed. By this method, the G+C content can be calculated at different "resolution". In an extreme case, the G+C content may be computed at a specific point, rather than in a window of finite size. This is particularly useful to analyze the fine variation of base composition along genomic sequences. As the first example, the variation of G+C content along each of 16 yeast chromosomes is analyzed. The G+C-rich regions with length larger than 5 kb sequences are detected and listed in details. It is found that each chromosome consists of several G+C-rich and G+C-poor regions alternatively, i.e., a mosaic structure. Another example is to analyze the G+C content for each of the two chromosomes of the Vibrio cholerae genome. Based on the variations of the G+C content in each chromosome, it is shown that some fragments in the Vibrio cholerae genome may have been transferred from other species. Especially, the position and size of the large integron island on the smaller chromosome was precisely predicted. This method would be a useful tool for analyzing genomic sequences.     

Determination of guanine-plus-cytosine content of bacterial DNA by dual-laser flow cytometry

A dual-laser flow cytometer was used to analyse different species of bacteria for the molar percentage of guanine-plus-cytosine (% G + C) without the need for DNA extraction or purification. Ethanol-fixed bacterial cells were stained with a combination of DNA-specific fluorochromes, Hoechst 33258 and chromomycin A3, which bind to AT- and GC-rich regions of DNA, respectively. A linear relationship (r = 0.99) was demonstrated between the log of the ratio of chromomycin A3 to Hoechst 33258 fluorescence and the log of the % G + C as determined by thermal denaturation (Tm) or buoyant density centrifugation (Bd) methods. Linearity was maintained for all bacterial species tested over the range of 28-67% G + C. A standard curve was constructed using five strains whose % G + C had been determined by other methods. From the equation describing this line, the % G + C values of nine other strains with known DNA base composition, together with the five strains used to construct the curve, were calculated using the chromomycin A3 to Hoechst 33258 ratio and were in agreement with values obtained by Tm, Bd or HPLC. The reproducibility of flow cytometric analysis (mean error 0.7% G + C) compared well with the reproducibility of other methods. Mixtures containing two species were also analysed. Two cell populations could be discerned in mixtures containing two species which differed in base composition by as little as 4% G + C.(ABSTRACT TRUNCATED AT 250 WORDS)     

An isochore map of the human genome based on the Z curve method

The distribution of the G+C content in the human genome has been studied by using a windowless technique derived from the Z curve method. The most important findings presented in this paper are twofold. First, abrupt variations of the G+C content along human chromosome sequences are the main variation patterns of G+C content. It is found that at some sites, the G+C content undergoes abrupt changes from a G+C-rich region to a G+C-poor region alternatively and vice versa. Second, it is shown that long domains with relatively homogeneous G+C content along each chromosome do exist. These domains are thought to be isochores, which usually have sharp boundaries. Consequently, 56 isochores longer than 3 Mb have been identified in chromosomes 1-22, X and Y. Boundaries, size and G+C content of each isochore identified are listed in detail. As an example to demonstrate the power of the method, the boundary between the Classes III and II isochores of the MHC sequence has been determined and found to be at 2,477,936, which is in good agreement with the experimental evidence. A homogeneity index is introduced to measure the homogeneity of G+C content in isochores. We emphasize that the homogeneity of G+C content is relative. The isochores in which the G+C content keeps absolutely constant do not exist. Isochore structures appear to be a basic organization of the human genome. Due to the relevance to many important biological functions, the clarification of isochore structures will provide much insight into the understanding of the human genome.     

Base pairing between Escherichia coli RNase P RNA and its substrate.

Base pairing between the substrate and the ribozyme has previously been shown to be essential for catalytic activity of most ribozymes, but not for RNase P RNA. By using compensatory mutations we have demonstrated the importance of Watson-Crick complementarity between two well-conserved residues in Escherichia coli RNase P RNA (M1 RNA), G292 and G293, and two residues in the substrate, +74C and +75C (the first and second C residues in CCA). We suggest that these nucleotides base pair (G292/+75C and G293/+74C) in the ribozyme-substrate complex and as a consequence the amino acid acceptor stem of the precursor is partly unfolded. Thus, a function of M1 RNA is to anchor the substrate through this base pairing, thereby exposing the cleavage site such that cleavage is accomplished at the correct position. Our data also suggest possible base pairing between U294 in M1 RNA and the discriminator base at position +73 of the precursor. Our findings are also discussed in terms of evolution.     

TCRαβ^＋CD^＋4CD^—8胸腺髓质细胞的表型分析

Phenotypic analysis of the medullary-type CD4+CD8- (CD4SP) thymocytes have revealed phenotypic heterogeneity within these cells. The phenotype of mature peripheral T cells is Qa-2+ HSA- CD69-, whereas in the medullary-type CD4SP thymocytes, the expression pattern of many markers were quite different, suggesting that the medullary-type CD4SP thymocytes may undergo phenotypic maturation. According to the results of two-color cytometry, seven discrete phenotypes were defined by the relative expression of Qa-2, HSA, CD69, 3G11 and 6C10: 3G11-6C10+CD69+HSAhi-->3G11+6C10+CD69+ HSAhi-->3G11+6C10-CD69+HSAint-->3G11+6C10- CD69-HSAint Qa-2(-)-->3G11+HSAlo/-Qa-2lo, at the same time, 3G11+6C10-CD69-HSAint Qa-2(-)-->3G11-HSAlo Qa-2(-)-->3G11-HSAlo/- Qa-2hi, the last two Qa-2 positive subsets could exit the thymus and home into periphery.     

Molecular mechanics calculations on a triple stranded DNA involving C+.G-T and T.A+-C mismatched base triples

We have carried out molecular modeling of a triple stranded pyrimidine(Y). purine(R): pyrimidine(Y) (where ':' refers to Watson-Crick and '.' to Hoogsteen bonding) DNA, formed by a homopurine (d-TGAGGAAAGAAGGT) and homo-pyrimidine (d-CTCCTTTCTTCC). Molecular mechanics calculations using NMR constraints have provided a detailed three dimensional structure of the triplex. The entire stretches of purine and the pyrimidine nucleotides have a conformation close to B-DNA. The three strands are held by the canonical C+.G:C and T.A:T hydrogen bonds. The structure also contains two mismatch C+.G-T and T.A+-C base triples which have been characterized for the first time. In the A+-C base-pair of the T.A+-C triple, both hydrogen donors are situated on the purine (A+(1N) and A+(6N)). We observe a unique hydrogen bonding interaction scheme in case of C+.G-T where one acceptor, G(60), is bonded to three donors (C+(3NH), C+(4NH2) and T(3NH)). Though the C+.G-T base triple is less stable than C+.G:C, it is significantly more stable than T.A:T. On the other hand, T.A+-C is as stable as the T.A:T base triad.     

Dependence of the Raman signature of genomic B-DNA on nucleotide base sequence.

The vibrational spectra of four genomic and two synthetic DNAs, encompassing a wide range in base composition [poly(dA-dT). poly(dA-dT), 0% G + C; Clostridium perfringens DNA, 27% G + C; calf thymus DNA, 42% G + C; Escherichia coli DNA, 50% G + C; Micrococcus luteus DNA, 72% G + C; poly(dG-dC).poly(dG-dC), 100% G + C] (dA: deoxyadenosine; dG: deoxyguanosine; dC: deoxycytidine; dT: thymidine), have been analyzed using Raman difference methods of high sensitivity. The results show that the Raman signature of B DNA depends in detail upon both genomic base composition and sequence. Raman bands assigned to vibrational modes of the deoxyribose-phosphate backbone are among the most sensitive to base sequence, indicating that within the B family of conformations major differences occur in the backbone geometry of AT- and GC-rich domains. Raman bands assigned to in-plane vibrations of the purine and pyrimidine bases-particularly of A and T-exhibit large deviations from the patterns expected for random base distributions, establishing that Raman hypochromic effects in genomic DNA are also highly sequence dependent. The present study provides a basis for future use of Raman spectroscopy to analyze sequence-specific DNA-ligand interactions. The demonstration of sequence dependency in the Raman spectrum of genomic B DNA also implies the capability to distinguish genomic DNAs by means of their characteristic Raman signatures.     

The effects of guanine and cytosine variation on dinucleotide frequency and amino acid composition in the human genome

Summary One hundred twelve human DNA sequences were analyzed with respect to dinucleotide frequency and amino acid composition. The variation in guanine and cytosine (G+C) content revealed: (1) at 2–3 and 3-1 doublet positions CG discrimination is attenuated at high G+C, but TA disfavor is enhanced, and (2) several amino acids are subject to G+C change. These findings have been reported in part for collections of sequences from various species. The present study confirms that in a single organism-the human-the G+C effects do exist. Aspects of the argument that connects G+C with protein thermal stability are also discussed.     

Crystal structure of pullulanase: evidence for parallel binding of oligosaccharides in the active site

The crystal structures of Klebsiella pneumoniae pullulanase and its complex with glucose (G1), maltose (G2), isomaltose (isoG2), maltotriose (G3), or maltotetraose (G4), have been refined at around 1.7-1.9A resolution by using a synchrotron radiation source at SPring-8. The refined models contained 920-1052 amino acid residues, 942-1212 water molecules, four or five calcium ions, and the bound sugar moieties. The enzyme is composed of five domains (N1, N2, N3, A, and C). The N1 domain was clearly visible only in the structure of the complex with G3 or G4. The N1 and N2 domains are characteristic of pullulanase, while the N3, A, and C domains have weak similarity with those of Pseudomonas isoamylase. The N1 domain was found to be a new type of carbohydrate-binding domain with one calcium site (CBM41). One G1 bound at subsite -2, while two G2 bound at -1 approximately -2 and +2 approximately +1, two G3, -1 approximately -3 and +2 approximately 0', and two G4, -1 approximately -4 and +2 approximately -1'. The two bound G3 and G4 molecules in the active cleft are almost parallel and interact with each other. The subsites -1 approximately -4 and +1 approximately +2, including catalytic residues Glu706 and Asp677, are conserved between pullulanase and alpha-amylase, indicating that pullulanase strongly recognizes branched point and branched sugar residues, while subsites 0' and -1', which recognize the non-reducing end of main-chain alpha-1,4 glucan, are specific to pullulanase and isoamylase. The comparison suggested that the conformational difference around the active cleft, together with the domain organization, determines the different substrate specificities between pullulanase and isoamylase.     

Molecular cloning and nucleotide sequence of the gene for the ribosomal protein S11 from the archaebacterium Halobacterium marismortui

The gene encoding ribosomal protein S11 (Escherichia coli S15 homologue) from Halobacterium marismortui was cloned employing two synthetic oligonucleotide mixtures, 23 and 32 bases in length, as hybridization probes. The nucleotide sequence of the gene and the adjacent 5'- and 3'-flanking regions (1300 base pairs) were then determined by the dideoxy chain termination method. Comparison of the nucleotide sequence of the H. marismortui S11 gene with that of the E. coli S15 gene (rpsO) showed that the 3'-end of the S11 gene can be aligned with the entire E. coli S15 gene, sharing 44% identical nucleotides. It has been found that the S11 gene has a higher G + C content (G + C = 65%) than that of the E. coli S15 gene (G + C = 53%). This increase in G + C content specifically shows up as a preference for G + C in the 3rd position of the codon. Upstream of the S11 gene, an archaebacterial promoter sequence (GGACTTTCA) and a putative ribosomal binding site (GCGGT) have been found, 88 and 15 (or 24) base pairs from the initiation codon of the gene. In addition, an open reading frame could be identified immediately after the stop codon for the S11 gene. Northern blotting analysis using the S11 coding region as probe has shown that the S11 gene is located on a 2.4-kilobase mRNA, suggesting that it is cotranscribed with other downstream gene(s).     

Transforming growth factor-β1 gene +869T/C,but not +915G/C polymorphism is associated with essential hypertension in a Chinese patient cohort

Many studies have suggested that transforming growth factor-β1 (TGF-β1) gene might be involved in the development of hypertension. However, results have been inconsistent. In this study, the authors performed a meta-analysis to investigate the associations of +869T/C and +915G/C polymorphisms in TGF-β1 gene with hypertension risk in Chinese. Published literature from PubMed, EMBASE, CNKI, CBM, and Wanfang Data were searched. Pooled odds ratio (OR) and 95% confidence interval (CI) were calculated using fixed or random-effects model. Nine studies (1,995 cases/1,840 controls) for +869T/C polymorphism and seven studies (1,547 cases/1,577 controls) for +915G/C polymorphism were included in the meta-analysis. The overall result showed that there was a statistically significant association between +869T/C polymorphism and hypertension risk (CC vs. TT: OR = 1.80, 95% CI 1.34–2.44). Similar results were found among two geographic locations and two subgroups with different sample size. However, no significant association was found for +915G/C polymorphism with the risk of hypertension (CC vs. GG: OR = 1.66, 95% CI 0.74–3.74). The meta-analysis indicated the significant association of +869T/C, but not +915G/C polymorphism with hypertension susceptibility. However, given the limited sample size, the associations warrant further investigation.     

Combined segregation-linkage analysis of plasma thrombin activatable fibrinolysis inhibitor (TAFI) antigen levels with TAFI gene polymorphisms

By decreasing plasminogen binding to fibrin surface, the thrombin activatable fibrinolysis inhibitor (TAFI) has been hypothesized to constitute an early marker for atherothrombotic diseases. Previous studies have shown that plasma TAFI levels exhibit a high interindividual variability that is only poorly explained by lifestyle factors. Several polymorphisms of the TAFI gene have been described, and a combination of a C+1542G substitution in the 3' untranslated region and an Ala147Thr amino acid change has been shown to explain 60% of TAFI variability in a sample of unrelated individuals. A segregation-linkage analysis was performed to determine whether these polymorphisms are directly involved in the genetic regulation of TAFI levels, or whether they are only markers in linkage disequilibrium (LD) with unmeasured TAFI-linked quantitative trait loci (QTLs). The sample consisted of 97 healthy nuclear families from the Stanislas Cohort. The C+1542G and Ala147Thr polymorphisms were in complete negative LD, with minor allele frequencies of 0.27 and 0.28, respectively. Results of the segregation-linkage analysis provided evidence of two TAFI-linked QTLs in LD with the two measured polymorphisms, which would explain 78% of the TAFI variance, as compared with 55% explained by the C+1542G and the Ala147Thr polymorphisms combined. The two putative QTLs would have minor allele frequencies of 0.45 and 0.32, respectively. The hypothesis that one of the measured polymorphisms is one of the QTLs was rejected. The putative QTLs also did not seem compatible with the other TAFI gene polymorphisms that we have previously described. More extensive sequencing of the TAFI gene is necessary to identify the functional variants.     

Unraveling functional and structural interactions between transmembrane domains IV and XI of NhaA Na+/H+ antiporter of Escherichia coli

A functionally important, interface domain between transmembrane segments (TMSs) IV and XI of the NhaA Na+/H+ antiporter of Escherichia coli has been unraveled. Scanning by single Cys replacements identified new mutations (F136C, G125C, and A137C) that cluster in one face of TMS IV and increase dramatically the Km of the antiporter. Whereas G125C, in addition, causes a drastic alkaline shift to the pH dependence of the antiporter, G338C alleviates the pH control of NhaA. Scanning by double Cys replacements (21 pairs of one replacement per TMS) identified genetically eight pairs of residues that showed very strong negative complementation. Cross-linking of the double mutants identified six double mutants (T132C/G338C, D133C/G338C, F136C/S342C, T132C/S342C, A137C/S342C, and A137C/G338C) of which pronounced intramolecular cross-linking defined an interface domain between the two TMSs. Remarkably, cross-linking by a short and rigid reagent (N,N'-o-phenylenedimaleimide) revived the Li+/H+ antiport activity, whereas a shorter reagent (1,2-ethanediyl bismethanethiosulfonate) revived both Na+/H+ and Li+/H+ antiporter activities and even the pH response of the dead mutant T132C/G338C. Hence, cross-linking at this position restores an active conformation of NhaA.     

Identification of Shewanella baltica as the most important H2S-producing species during iced storage of Danish marine fish

Shewanella putrefaciens has been considered the main spoilage bacteria of low-temperature stored marine seafood. However, psychrotropic Shewanella have been reclassified during recent years, and the purpose of the present study was to determine whether any of the new Shewanella species are important in fish spoilage. More than 500 H2S-producing strains were isolated from iced stored marine fish (cod, plaice, and flounder) caught in the Baltic Sea during winter or summer time. All strains were identified as Shewanella species by phenotypic tests. Different Shewanella species were present on newly caught fish. During the warm summer months the mesophilic human pathogenic S. algae dominated the H2S-producing bacterial population. After iced storage, a shift in the Shewanella species was found, and most of the H2S-producing strains were identified as S. baltica. The 16S rRNA gene sequence analysis confirmed the identification of these two major groups. Several isolates could only be identified to the genus Shewanella level and were separated into two subgroups with low (44%) and high (47%) G+C mol%. The low G+C% group was isolated during winter months, whereas the high G+C% group was isolated on fish caught during summer and only during the first few days of iced storage. Phenotypically, these strains were different from the type strains of S. putrefaciens, S. oneidensis, S. colwelliana, and S. affinis, but the high G+C% group clustered close to S. colwelliana by 16S rRNA gene sequence comparison. The low G+C% group may constitute a new species. S. baltica, and the low G+C% group of Shewanella spp. strains grew well in cod juice at 0 degrees C, but three high G+C Shewanella spp. were unable to grow at 0 degrees C. In conclusion, the spoilage reactions of iced Danish marine fish remain unchanged (i.e., trimethylamine-N-oxide reduction and H2S production); however, the main H2S-producing organism was identified as S. baltica.     

新变种上犹杜鹃的ITS分子标记证据分析

对杜鹃花科杜鹃属的新变种上犹杜鹃进行了ITS分子标记分析。结果表明:(1)上犹杜鹃与原变种毛果杜鹃在系统树上没有聚合为一支;(2)上犹杜鹃ITS序列密码子第二位上的G+C含量与原变种毛果杜鹃差异明显,而且编码位上的G+C含量(G+Cc)也与原变种毛果杜鹃有差异;(3)上犹杜鹃与原变种毛果杜鹃的形态特征有明显的区别,并具间断性特征;(4)密码子偏向性指数CBI分析正好反映了上犹杜鹃处于"变种"等级阶段的特点。基于这些证据,认为上犹杜鹃处理为种以下等级"变种"比较合理。