Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Phylogenetic structures of the Holarctic Sorex araneus group and its relationships with S-samniticus, as inferred from mtDNA sequences

The shrews of the Sorex araneus group, characterized by the sexual chromosome complex XY1, Y2 have been intensively studied by morphological, karyotypical, and biochemical analyses. Nevertheless, the phylogenetic relationships among the species belonging to the araneus complex are still under debate, as different approaches gave often contradictory results. In this paper, partial nucleotide sequences of the mitochondrial DNA cytochrome b gene (1011 bp) were determined for 6 species of the araneus group from Eurasia and North America. We also included in the data set the sequences of Sorex samniticus, whose relationships with the araneus group remain controversial. Three other species representing two major karyological groups were also examined. Both parsimony and distance trees strongly support the monophyly of the araneus group. Sorex sumniticus is significantly more closely related to the araneus complex than to the other species included in the analysis. Based on the branching pattern within the araneus group, an attempt has been made to reconstruct the colonization history of the Holarctic region.

MyHits: improvements to an interactive resource for analyzing protein sequences.

The MyHits web site (http://myhits.isb-sib.ch) is an integrated service dedicated to the analysis of protein sequences. Since its first description in 2004, both the user interface and the back end of the server were improved. A number of tools (e.g. MAFFT, Jacop, Dotlet, Jalview, ESTScan) were added or updated to improve the usability of the service. The MySQL schema and its associated API were revamped and the database engine (HitKeeper) was separated from the web interface. This paper summarizes the current status of the server, with an emphasis on the new services.

Reconstruction of metal pollution and recent sedimentation processes in Havana Bay (Cuba): A tool for coastal ecosystem management.

Since 1998 the highly polluted Havana Bay ecosystem has been the subject of a mitigation program. In order to determine whether pollution-reduction strategies were effective, we have evaluated the historical trends of pollution recorded in sediments of the Bay. A sediment core was dated radiometrically using natural and artificial fallout radionuclides. An irregularity in the (210)Pb record was caused by an episode of accelerated sedimentation. This episode was dated to occur in 1982, a year coincident with the heaviest rains reported in Havana over the XX century. Peaks of mass accumulation rates (MAR) were associated with hurricanes and intensive rains. In the past 60 years, these maxima are related to strong El Niño periods, which are known to increase rainfall in the north Caribbean region. We observed a steady increase of pollution (mainly Pb, Zn, Sn, and Hg) since the beginning of the century to the mid 90s, with enrichment factors as high as 6. MAR and pollution decreased rapidly after the mid 90s, although some trace metal levels remain high. This reduction was due to the integrated coastal zone management program introduced in the late 90s, which dismissed catchment erosion and pollution.

Comparison of human chromosome 21 conserved nongenic sequences (CNGs) with the mouse and dog genomes shows that their selective constraint is independent of their genic environment.

The analysis of conservation between the human and mouse genomes resulted in the identification of a large number of conserved nongenic sequences (CNGs). The functional significance of this nongenic conservation remains unknown, however. The availability of the sequence of a third mammalian genome, the dog, allows for a large-scale analysis of evolutionary attributes of CNGs in mammals. We have aligned 1638 previously identified CNGs and 976 conserved exons (CODs) from human chromosome 21 (Hsa21) with their orthologous sequences in mouse and dog. Attributes of selective constraint, such as sequence conservation, clustering, and direction of substitutions were compared between CNGs and CODs, showing a clear distinction between the two classes. We subsequently performed a chromosome-wide analysis of CNGs by correlating selective constraint metrics with their position on the chromosome and relative to their distance from genes. We found that CNGs appear to be randomly arranged in intergenic regions, with no bias to be closer or farther from genes. Moreover, conservation and clustering of substitutions of CNGs appear to be completely independent of their distance from genes. These results suggest that the majority of CNGs are not typical of previously described regulatory elements in terms of their location. We propose models for a global role of CNGs in genome function and regulation, through long-distance cis or trans chromosomal interactions.

Complete Nucleotide Sequence of Mouse Mammary Tumor Virus from JYG Chinese Wild Mice: Absence of Bacterial Insertion Sequences in the Cloned Viral gag Gene.

Mammary tumors of a newly isolated strain of Chinese wild mouse (JYG mouse) harbor exogenous mouse mammary tumor virus (MMTV). The complete nucleotide sequence of exogenous JYG-MMTV was determined on the proviral 5' long terminal repeat (LTR)(partial)-gag-pol-env-3' LTR (partial) fragment cloned into a plasmid vector and the cDNA sequence from JYG-MMTV producing cells. Similarly to the other MMTV species the LTR of JYG-MMTV contains an open reading frame (ORF). The amino acid sequence of the JYG-MMTV ORF resembles that of SW-MMTV (92% identity) and endogenous Mtv-7 (93% identity) especially at the C-terminal region. Thus, a functional similarity in T-cell receptor V beta recognition as a superantigen is implicated among these MMTV species. Analysis of the viral gag nucleotide sequence revealed that this gene is not disrupted by the bacterial insertion sequence IS1 or IS2, which have been reported to be present in the majority of the plasmids containing the gag region. Comparison of amino acid sequences of JYG-MMTV with those of BR6-MMTV showed that over 96% of the amino acids of gag, pol, protease and env products are identical. These results suggest the intact nature of the nucleotide sequence of the near full-length MMTV genome cloned in the plasmid.

The complete amino acid sequence for brain beta spectrin (beta fodrin): relationship to globin sequences.

The amino acid sequence of mouse brain beta spectrin (beta fodrin), deduced from the nucleotide sequence of complementary DNA clones, reveals that this non-erythroid beta spectrin comprises 2363 residues, with a molecular weight of 274,449 Da. Brain beta spectrin contains three structural domains and we suggest the position of several functional domains including f-actin, synapsin I, ankyrin and spectrin self association sites. Analysis of deduced amino acid sequences indicated striking homology and similar structural characteristics of brain beta spectrin repeats beta 11 and beta 12 to globins. In vitro analysis has demonstrated that heme is capable of specific attachment to brain spectrin, suggesting possible new functions in electron transfer, oxygen binding, nitric oxide binding or heme scavenging.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Localized in vivo hyperpolarization transfer sequences.

In vivo localized and fully adiabatic homonuclear and heteronuclear polarization transfer experiments were designed and performed in the rat brain at 9.4 T after infusion of hyperpolarized sodium [1,2-(13)C(2)] and sodium [1-(13)C] acetate. The method presented herein leads to highly enhanced in vivo detection of short-T(1) (13)C as well as attached protons. This indirect detection scheme allows for probing additional molecular sites in hyperpolarized substrates and their metabolites and can thus lead to improved spectral resolution such as in the case of (13)C-acetate metabolism.

Temporal and spatial origin of Gesneriaceae in the New World inferred from plastid DNA sequences

Gesneriaceae are represented in the New World (NW) by a major clade (c. 1000 species) currently recognized as subfamily Gesnerioideae. Radiation of this group occurred in all biomes of tropical America and was accompanied by extensive phenotypic and ecological diversification. Here we performed phylogenetic analyses using DNA sequences from three plastid loci to reconstruct the evolutionary history of Gesnerioideae and to investigate its relationship with other lineages of Gesneriaceae and Lamiales. Our molecular data confirm the inclusion of the South Pacific Coronanthereae and the Old World (OW) monotypic genus Titanotrichum in Gesnerioideae and the sister-group relationship of this subfamily to the rest of the OW Gesneriaceae. Calceolariaceae and the NW genera Peltanthera and Sanango appeared successively sister to Gesneriaceae, whereas Cubitanthus, which has been previously assigned to Gesneriaceae, is shown to be related to Linderniaceae. Based on molecular dating and biogeographical reconstruction analyses, we suggest that ancestors of Gesneriaceae originated in South America during the Late Cretaceous. Distribution of Gesneriaceae in the Palaeotropics and Australasia was inferred as resulting from two independent long-distance dispersals during the Eocene and Oligocene, respectively. In a short time span starting at 34 Mya, ancestors of Gesnerioideae colonized several Neotropical regions including the tropical Andes, Brazilian Atlantic forest, cerrado, Central America and the West Indies. Subsequent diversification within these areas occurred largely in situ and was particularly extensive in the mountainous systems of the Andes, Central America and the Brazilian Atlantic forest. Only two radiations account for 90% of the diversity of Gesneriaceae in the Brazilian Atlantic forest, whereas half of the species richness in the northern Andes and Central America originated during the last 10 Myr from a single radiation.

Scattering of repetitive DNA sequences in the albumin and vitellogenin gene loci of Xenopus laevis.

We have analyzed middle repetitive DNA in the albumin and vitellogenin gene families of Xenopus laevis. Mapping specific repetitive DNA sequences derived from introns of the A1 vitellogenin gene reveals that these sequences are scattered within and around the four vitellogenin genes (A1, A2, B1 and B2) and the two albumin genes (74 kd and 68 kd). Three repetitive DNA elements present in the A1 vitellogenin transcriptional unit are also located in introns of the 74 kd albumin gene. This apparently random distribution of middle repetitive DNA in the two gene families suggests that the analyzed sequences are not involved in gene regulation, but rather that they might represent unstable genetic elements. This hypothesis is further supported by the finding that size polymorphism in the A1 vitellogenin gene and in the 74 kd albumin gene is correlated with the presence or absence of repetitive DNA.