Quartet mapping and the extent of lateral transfer in bacterial genomes

Several recent analyses have used quartet-based methods to assess the congruence among phylogenies derived for large sets of genes from prokaryotic genomes. The principal conclusion from these studies is that lateral gene transfer (LGT) has blurred prokaryotic phylogenies to such a degree that the darwinian scheme of treelike evolution might be abandoned in favor of a net or web. Here, we focus on one of these methods, quartet mapping, and show that its application can lead to overestimation of the extent of inferred LGT in prokaryotes, particularly when applied to distantly related taxa.     

Relationship between digital information and thermodynamic stability in bacterial genomes

Ever since the introduction of the Watson-Crick model, numerous efforts have been made to fully characterize the digital information content of the DNA. However, it became increasingly evident that variations of DNA configuration also provide an “analog” type of information related to the physicochemical properties of the DNA, such as thermodynamic stability and supercoiling. Hence, the parallel investigation of the digital information contained in the base sequence with associated analog parameters is very important for understanding the coding capacity of the DNA. In this paper, we represented analog information by its thermodynamic stability and compare it with digital information using Shannon and Gibbs entropy measures on the complete genome sequences of several bacteria, including Escherichia coli (E. coli), Bacillus subtilis (B. subtilis), Streptomyces coelicolor (S. coelicolor), and Salmonella typhimurium (S. typhimurium). Furthermore, the link to the broader classes of functional gene groups (anabolic and catabolic) is examined. Obtained results demonstrate the couplings between thermodynamic stability and digital sequence organization in the bacterial genomes. In addition, our data suggest a determinative role of the genome-wide distribution of DNA thermodynamic stability in the spatial organization of functional gene groups.     

Analysis of septins across kingdoms reveals orthology and new motifs

Background  

Septins are cytoskeletal GTPase proteins first discovered in the fungus Saccharomyces cerevisiae where they organize the septum and link nuclear division with cell division. More recently septins have been found in animals where they are important in processes ranging from actin and microtubule organization to embryonic patterning and where defects in septins have been implicated in human disease. Previous studies suggested that many animal septins fell into independent evolutionary groups, confounding cross-kingdom comparison.     

Disease association mapping in Drosophila can be replicated in the wild

Association and linkage mapping have become important tools in understanding the genetics of complex traits, including diseases in humans. As the success of association mapping is reduced by small effect sizes and limited power, linkage studies in laboratory-based model systems are still heavily used. But whether the results of these studies can be replicated in natural populations has been questioned. Here, we show that a polymorphism in the gene ref(2)P, which had previously been linked to sigma virus resistance in Drosophila melanogaster under laboratory conditions, also provides resistance against the virus in female flies in a wild population in the field. This genetic association is thus upheld in spite of a known genotype-by-genotype interaction and environmental variation.     

Stochastic problems in information transfer across the plasma membrane

The surfaces of many cells are viscous fluids; consequently, most membrane proteins are able to diffuse laterally, in a more or less random fashion, with diffusion coefficients typically of order 10⁻¹⁰ cm²/sec. If a molecule (ligand) in solution outside the cell and a protein molecule on the surface (receptor) each have two or more sites at which they can interact with one another, large, branched receptor-ligand networks can form on the cell surface by virtue of the chemical interactions that surface fluidity permits. Evidence from a variety of systems indicates that such receptor clustering plays a role in the sequence of events leading to cellular activity. This paper describes a number of mathematical problems that arise in the analysis of experiments in which clustering occurs. I begin by reviewing methods for finding the time evolution of the cluster size distribution function in terms of reaction rate constants. The methods solve an essentially infinite system of coupled nonlinear differential equations. Next, the rate constants are analyzed, the Brownian motion problems that arise in attempting to understand ligand recognition are described and relevant experimental systems are discussed. Finally the notion of ligand as a signal amplifier is introduced—an idea that emerges naturally from the requirement that receptors be clustered for a finite amount of time before a signal can be transmitted.     

Order and disorder in bacterial genomes

The availability of sequenced bacterial genomes allows a deeper understanding of their organizational features that are related with fundamental cellular processes such as coordinated gene expression, chromosome replication and cell division. Nevertheless, recent genome comparisons and experimental work highlighted the fluidity of bacterial chromosomes, including genome rearrangements that imperil the selective features of chromosome order. As a result, the clash between elements generating rearrangements and chromosome organization is a classic case of evolutionary conflict.     

Identification and functional characterization of gene components of Type VI Secretion system in bacterial genomes

A new secretion system, called the Type VI Secretion system (T6SS), was recently reported in Vibrio cholerae, Pseudomonas aeruginosa and Burkholderia mallei. A total of 18 genes have been identified to be belonging to this secretion system in V. cholerae. Here we attempt to identify presence of T6SS in other bacterial genomes. This includes identification of orthologous sequences, conserved motifs, domains, families, 3D folds, genomic islands containing T6SS components, phylogenetic profiles and protein-protein association of these components. Our analysis indicates presence of T6SS in 42 bacteria and its absence in most of their non-pathogenic species, suggesting the role of T6SS in imparting pathogenicity to an organism. Analysis of genomic regions containing T6SS components, phylogenetic profiles and protein-protein association of T6SS components indicate few additional genes which could be involved in this secretion system. Based on our studies, functional annotations were assigned to most of the components. Except one of the genes, we could group all the other genes of T6SS into those belonging to the puncturing device, and those located in the outer membrane, transmembrane and inner membrane. Based on our analysis, we have proposed a model of T6SS and have compared the same with the other bacterial secretion systems.     

Denitrifying genes in bacterial and Archaeal genomes

Denitrification, the reduction of nitrate or nitrite to nitrous oxide or dinitrogen, is the major mechanism by which fixed nitrogen returns to the atmosphere from soil and water. Although the denitrifying ability has been found in microorganisms belonging to numerous groups of bacteria and Archaea, the genes encoding the denitrifying reductases have been studied in only few species. Recent investigations have led to the identification of new classes of denitrifying reductases, indicating a more complex genetic basis of this process than previously recognized. The increasing number of genome sequencing projects has opened a new way to study the genetics of the denitrifying process in bacteria and Archaea. In this review, we summarized the current knowledge on denitrifying genes and compared their genetic organizations by using new sequences resulting from the analysis of finished and unfinished microbial genomes with a special attention paid to the clustering of genes encoding different classes of reductases. In addition, some evolutionary relationships between the structural genes are presented.     

基因组信息的整合与植物功能基因组学研究的策略

近年来,随着许多植物基因组测序和可利用序列的增加,相继建立了一些基于靶基因诱变的“反向”遗传学研究策略,如T—DNA诱变、基因敲除、基因沉默和超表达分析等。同时,DNA微阵列和基因芯片技术的发展使得快速、定量检测植物发育不同时期和不同组织器官的基因转录时空变化成为现实。作图技术的改进和来自不同物种基因组信息的整合也正在加速图谱克隆程序的简化和发展。因此,随着生物基因组测序工作日益增多,整合不同类群植物基因组的信息和资源,在植物功能基因组学研究中的重要性日趋显著。     

Structural and functional comparative mapping between the Brassica A genomes in allotetraploid Brassica napus and diploid Brassica rapa

Brassica napus (AACC genome) is an important oilseed crop that was formed by the fusion of the diploids B. rapa (AA) and B. oleracea (CC). The complete genomic sequence of the Brassica A genome will be available soon from the B. rapa genome sequencing project, but it is not clear how informative the A genome sequence in B. rapa (A(r)) will be for predicting the structure and function of the A subgenome in the allotetraploid Brassica species B. napus (A(n)). In this paper, we report the results of structural and functional comparative mapping between the A subgenomes of B. napus and B. rapa based on genetic maps that were anchored with bacterial artificial chromosomes (BACs)-sequence of B. rapa. We identified segmental conservation that represented by syntenic blocks in over one third of the A genome; meanwhile, comparative mapping of quantitative trait loci for seed quality traits identified a dozen homologous regions with conserved function in the A genome of the two species. However, several genomic rearrangement events, such as inversions, intra- and inter-chromosomal translocations, were also observed, covering totally at least 5% of the A genome, between allotetraploid B. napus and diploid B. rapa. Based on these results, the A genomes of B. rapa and B. napus are mostly functionally conserved, but caution will be necessary in applying the full sequence data from B. rapa to the B. napus as a result of genomic rearrangements in the A genome between the two species.     

Paralogy and orthology of tyrosine kinases that can extend the life span of Caenorhabditis elegans

Modification of any one of three transmembrane protein tyrosine kinase (PTK) genes, old-1, old-2 (formerly tkr-1 and tkr-2, respectively), and daf-2 can extend the mean and maximum life span of the nematode Caenorhabditis elegans. To identify paralogs and orthologs, we delineated relationships between these three PTKs and all known transmembrane PTKs and all known mammalian nontransmembrane PTKs using molecular phylogenetics. The tree includes a number of invertebrate receptor PTKs and a novel mammalian receptor PTK (inferred from the expressed-sequence tag database) that have not previously been analyzed. old-1 and old-2 were found to be members of a surprisingly large C. elegans PTK family having 16 members. Interestingly, only four members of this transmembrane family appeared to have receptor domains (immunoglobulin-like in each case). The C-terminal domain of this family was found to have a unique sequence motif that could be important for downstream signaling. Among mammalian PTKs, the old-1/old-2 family appeared to be most closely related to the Pdgfr, Fgfr, Ret, and Tie/Tek families. However, these families appeared to have split too early from the old-1/old-2 family to be orthologs, suggesting that a mammalian ortholog could yet be discovered. An extensive search of the expressed-sequence tag database suggested no additional candidate orthologs. In contrast to old-1 and old-2, daf-2 had no C. elegans paralogs. Although daf-2 was most closely related to the mammalian insulin receptor family, a hydra insulin receptor-like sequence suggested that daf-2 might not be an ortholog of the insulin receptor family. Among PTKs, the old-1/old-2 family and daf-2 were not particularly closely related, raising the possibility that other PTK families might extend life span. On a more general note, our survey of the expressed-sequence tag database suggested that few, if any, additional mammalian PTK families are likely to be discovered. The one novel family that was discovered could represent a novel oncogene family, given the prevalence of oncogenes among PTKs. Finally, the PTK tree was consistent with nematodes and fruit flies being as divergent as nematodes and mammals, suggesting that life extension mechanisms shared by nematodes and fruit flies would be reasonable candidates for extending mammalian life spans.     

Inferring functional constraints and divergence in protein families using 3D mapping of phylogenetic information

Comparative sequence analysis has been used to study specific questions about the structure and function of proteins for many years. Here we propose a knowledge-based framework in which the maximum likelihood rate of evolution is used to quantify the level of constraint on the identity of a site. We demonstrate that site-rate mapping on 3D structures using datasets of rhodopsin-like G-protein receptors and alpha- and beta-tubulins provides an excellent tool for pinpointing the functional features shared between orthologous and paralogous proteins. In addition, functional divergence within protein families can be inferred by examining the differences in the site rates, the differences in the chemical properties of the side chains or amino acid usage between aligned sites. Two novel analytical methods are introduced to characterize rate- independent functional divergence. These are tested using a dataset of two classes of HMG-CoA reductases for which only one class can perform both the forward and reverse reaction. We show that functionally divergent sites occur in a cluster of sites interacting with the catalytic residues and that this information should facilitate the design of experimental strategies to directly test functional properties of residues.     

Ongoing evolution of strand composition in bacterial genomes

We tried to identify the substitutions involved in the establishment of replication strand bias, which has been recognized as an important evolutionary factor in the evolution of bacterial genomes. First, we analyzed the composition asymmetry of 28 complete bacterial genomes and used it to test the possibility that asymmetric deamination of cytosine might be at the origin of the bias. The model showed significant correlation to the data but left unexplained a significant portion of the variance and indicated a systematic underestimation of GC skews in comparison with TA skews. Second, we analyzed the substitutions acting on the genes from five fully sequenced Chlamydia genomes that had not suffered strand switch since speciation. This analysis showed that substitutions were not at equilibrium in Chlamydia trachomatis or in C. muridarum and that strand bias is still an on-going process in these genes. Third, we identified substitutions involved in the adaptation of genes that had switched strands after speciation. These genes adapted quickly to the skewed composition of the new strand, mostly due to C-->T, A-->G, and C-->G asymmetric substitutions. This observation was reinforced by the analysis of genes that switched strands after divergence between Bacillus subtilis and B. halodurans. Finally, we propose a more extended model based on the analysis of the substitution asymmetries of CHLAMYDIA: This model fits well with the data provided by bacterial genomes presenting strong strand bias.     

Comparative genomics of small RNAs in bacterial genomes

In recent years, various families of small non-coding RNAs (sRNAs) have been discovered by experimental and computational approaches, both in bacterial and eukaryotic genomes. Although most of them await elucidation of their function, it has been reported that some play important roles in gene regulation. Here we carried out comparative genomics analysis of possible sRNAs that are computationally identified in 30 bacterial genomes from gamma- and alpha-proteobacteria and Deinococcus radiodurans. Identified sRNAs are clustered by a complete-linkage clustering method to see conservation among the organisms. On average, sRNAs are found in approximately 30% of intergenic regions of each genome sequence. Of these, 25.7% are conserved among three or more organisms. Approximately 60% of the conserved sRNAs do not locate in orthologous intergenic regions, implying that sRNAs may be shuffled their positions in genomes. The current study implies that sRNAs may be involved in a more extensive range of functions in bacteria.     

Synonymous and nonsynonymous polymorphisms versus divergences in bacterial genomes

Comparison of the ratio of nonsynonymous to synonymous polymorphisms within species with the ratio of nonsynonymous to synonymous substitutions between species has been widely used as a supposed indicator of positive Darwinian selection, with the ratio of these 2 ratios being designated as a neutrality index (NI). Comparison of genome-wide polymorphism within 12 species of bacteria with divergence from an outgroup species showed substantial differences in NI among taxa. A low level of nonsynonymous polymorphism at a locus was the best predictor of NI < 1, rather than a high level of nonsynonymous substitution between species. Moreover, genes with NI < 1 showed a strong tendency toward the occurrence of rare nonsynonymous polymorphisms, as expected under the action of ongoing purifying selection. Thus, our results are more consistent with the hypothesis that a high relative rate of between-species nonsynonymous substitution reflects mainly the action of purifying selection within species to eliminate slightly deleterious mutations rather than positive selection between species. This conclusion is consistent with previous results highlighting an important role of slightly deleterious variants in bacterial evolution and suggests caution in the use of the McDonald-Kreitman test and related statistics as tests of positive selection.