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.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Classification of urban multi-angular image sequences by aligning their manifolds

When dealing with multi-angular image sequences, problems of reflectance changes due either to illumination and acquisition geometry, or to interactions with the atmosphere, naturally arise. These phenomena interplay with the scene and lead to a modification of the measured radiance: for example, according to the angle of acquisition, tall objects may be seen from top or from the side and different light scatterings may affect the surfaces. This results in shifts in the acquired radiance, that make the problem of multi-angular classification harder and might lead to catastrophic results, since surfaces with the same reflectance return significantly different signals. In this paper, rather than performing atmospheric or bi-directional reflection distribution function (BRDF) correction, a non-linear manifold learning approach is used to align data structures. This method maximizes the similarity between the different acquisitions by deforming their manifold, thus enhancing the transferability of classification models among the images of the sequence.

Hepatic hemangiomas: Factors associated with T2 shine-through effect on diffusion-weighted MR sequences.

PURPOSE: To determine the frequency and factors associated with the presence of T2 shine-through effect in hepatic hemangiomas on diffusion-weighted (DW) magnetic resonance (MR) sequences. MATERIALS AND METHODS: This retrospective study was approved by institutional review board with waiver of informed consent. One hundred forty-nine consecutive patients with 388 hepatic hemangiomas who underwent a liver MR between January 2010 and November 2011 were included. MR analysis evaluated the lesion characteristics (signal intensities and enhancement patterns (classical, rapidly filling, delayed filling)), the presence of T2 shine-through effect on DW sequences (b values of 0, 150, and 600s/mm(2)), and apparent diffusion coefficient (ADC) values. Multivariate analysis was performed to study the factors associated with the T2 shine-through effect. RESULTS: T2 shine-through effect was observed in 204/388 (52.6%) of hepatic hemangiomas and in 100 (67.1%) patients. Mean ADC value of hemangiomas with T2 shine-through effect was significantly lower than hemangiomas without (2.0±0.48 vs 2.38±0.45, P<.0001). On multivariate analysis, high signal intensity on fat-suppressed T2-weighted fast spin-echo images, hemangiomas with classical or delayed enhancement, and the ADC of the liver were the only significant factors associated with T2 shine-through effect. CONCLUSION: T2 shine-through effect is commonly observed in hepatic hemangiomas and is related to hemangiomas characteristics. Radiologists should be aware of this phenomenon which could lead to misdiagnosis. Its presence should not question the diagnosis of hemangiomas when typical MR findings are found.

Molecular phylogenetics of soricid shrews (Mammalia) based on mitochondrial cytochrome b gene sequences: with special reference to the Soricinae

Molecular phylogeny of soricid shrews (Soricidae, Eulipotyphla, Mammalia) based on 1140 bp mitochondrial cytochrome b gene (cytb) sequences was inferred by the maximum likelihood (ML) method. All 13 genera of extant Soricinae and two genera of Crocidurinae were included in the analyses. Anourosorex was phylogenetically distant from the main groupings within Soricinae and Crocidurinae in the ML tree. Thus, it could not be determined to which subfamily Anourosorex should be assigned: Soricinae, Crocidurinae or a new subfamily. Soricinae (excluding Anourosorex) should be divided into four tribes: Neomyini, Notiosoricini, Soricini and Blarinini. However, monophyly of Blarinini was not robust in the present data set. Also, branching orders among tribes of Soricinae and those among genera of Neomyini could not be determined because of insufficient phylogenetic information of the cytb sequences. For water shrews of Neomyini (Chimarrogale, Nectogale and Neomys), monophyly of Neomys and the Chimarrogale-Nectogale group could not be verified, which implies the possibility of multiple origins for the semi-aquatic mode of living among taxa within Neomyini. Episoriculus may contain several separate genera. Blarinella was included in Blarinini not Soricini, based on the cytb sequences, but the confidence level was rather low; hence more phylogenetic information is needed to determine its phylogenetic position. Furthermore, some specific problems of taxonomy of soricid shrews were clarified, for example phylogeny of local populations of Notiosorex crawfordi, Chimarrogale himalayica and Crocidura attenuata.

Optimization Of Promoter Sequences In Order To Mimic Physiological Pde6b Expression

Purpose: Many retinal degenerations result from defective retina-specific gene expressions. Thus, it is important to understand how the expression of a photoreceptor-specific gene is regulated in vivo in order to achieve successful gene therapy. The present study aims to design an AAV2/8 vector that can regulate the transcript level in a physiological manner to replace missing PDE6b in Rd1 and Rd10 mice. In previous studies (Ogieta, et al., 2000), the short 5' flanking sequence of the human PDE6b gene (350 bp) was shown to be photoreceptor-specific in transgenic mice. However, the efficiency and specificity of the 5' flanking region of the human PDE6b was not investigated in the context of gene therapy during retinal degeneration. In this study, two different sequences of the 5' flanking region of the human PDE6b gene were studied as promoter elements and their expression will be tested in wild type and diseased retinas (Rd 10 mice).Methods: Two 5' flanking fragments of the human PDE6b gene: (-93 to +53 (150 bp) and -297 to +53 (350 bp)) were cloned in different plasmids in order to check their expression in vitro and in vivo by constructing an AAV2/8 vector. These elements drove the activity of either luciferase (pGL3 plasmids) or EGFP. jetPEI transfection in Y 79 cells was used to evaluate gene expression through luciferase activity. Constructs encoding EGFP under the control of the two promoters were performed in AAV2.1-93 (or 297)-EGFP plasmids to produce AAV2/8 vectors.Results: When pGL3-93 (150 bp) or pGL3-297 (350 bp) were transfected in the Y-79 cells, the smaller fragment (150 bp) showed higher gene expression compared to the 350 bp element and to the SV40 control, as previously reported. The 350 bp drove similar levels of expression when compared to the SV40 promoter. In view of these results, the fragments (150 bp or 350 bp) were integrated into the AAV2.1-EGFP plasmid to produce AAV2/8 vector, and we are currently evaluating the efficiency and specificity of the produced constructs in vivo in normal and diseased retinas.Conclusions: Comparisons of these vectors with vectors bearing ubiquitous promoters should reveal which construct is the most suitable to drive efficient and specific gene expression in diseased retinas in order to restore a normal function on the long term.

Stratigraphy of the Cenomanian-Turonian Oceanic Anoxic Event OAE2 in shallow shelf sequences of NE Egypt

Two shallow water late Cenomanian to early Turonian sequences of NE Egypt have been investigated to evaluate the response to OAE2. Age control based on calcareous nannoplankton, planktic foraminifera and ammonite biostratigraphies integrated with delta(13)C stratigraphy is relatively good despite low diversity and sporadic occurrences. Planktic and benthic foraminiferal faunas are characterized by dysoxic, brackish and mesotrophic conditions, as indicated by low species diversity, low oxygen and low salinity tolerant planktic and benthic species, along with oyster-rich limestone layers. In these subtidal to inner neritic environments the OAE2 delta(13)C excursion appears comparable and coeval to that of open marine environments. However, in contrast to open marine environments where anoxic conditions begin after the first delta(13)C peak and end at or near the Cenomanian-Turonian boundary, in shallow coastal environments anoxic conditions do not appear until the early Turonian. This delay in anoxia appears to be related to the sea-level transgression that reached its maximum in the early Turonian, as observed in shallow water sections from Egypt to Morocco. (C) 2011 Elsevier Ltd. All rights reserved.

trome, trEST and trGEN: databases of predicted protein sequences.

We previously introduced two new protein databases (trEST and trGEN) of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Here, we present the updates made on these two databases plus a new database (trome), which uses alignments of EST data to HTG or full genomes to generate virtual transcripts and coding sequences. This new database is of higher quality and since it contains the information in a much denser format it is of much smaller size. These new databases are in a Swiss-Prot-like format and are updated on a weekly basis (trEST and trGEN) or every 3 months (trome). They can be downloaded by anonymous ftp from ftp://ftp.isrec.isb-sib.ch/pub/databases.

Examination of Heterogeneous Crossing Sequences Between Toner and Rollerball Pen Strokes by Digital Microscopy and 3-D Laser Profilometry

The determination of line crossing sequences between rollerball pens and laser printers presents difficulties that may not be overcome using traditional techniques. This research aimed to study the potential of digital microscopy and 3-D laser profilometry to determine line crossing sequences between a toner and an aqueous ink line. Different paper types, rollerball pens, and writing pressure were tested. Correct opinions of the sequence were given for all case scenarios, using both techniques. When the toner was printed before the ink, a light reflection was observed in all crossing specimens, while this was never observed in the other sequence types. The 3-D laser profilometry, more time-consuming, presented the main advantage of providing quantitative results. The findings confirm the potential of the 3-D laser profilometry and demonstrate the efficiency of digital microscopy as a new technique for determining the sequence of line crossings involving rollerball pen ink and toner. With the mass marketing of laser printers and the popularity of rollerball pens, the determination of line crossing sequences between such instruments is encountered by forensic document examiners. This type of crossing presents difficulties with optical microscopic line crossing techniques involving ballpoint pens or gel pens and toner (1-4). Indeed, the rollerball's aqueous ink penetrates through the toner and is absorbed by the fibers of the paper, leaving the examiner with the impression that the toner is above the ink even when it is not (5). Novotny and Westwood (3) investigated the possibility of determining aqueous ink and toner crossing sequences by microscopic observation of the intersection before and after toner removal. A major disadvantage of their study resides in destruction of the sample by scraping off the toner line to see what was underneath. The aim of this research was to investigate the ways to overcome these difficulties through digital microscopy and three-dimensional (3-D) laser profilometry. The former was used as a technique for the determination of sequences between gel pen and toner printing strokes, but provided less conclusive results than that of an optical stereomicroscope (4). 3-D laser profilometry, which allows one to observe and measure the topography of a surface, has been the subject of a number of recent studies in this area. Berx and De Kinder (6) and Schirripa Spagnolo (7,8) have tested the application of laser profilometry to determine the sequence of intersections of several lines. The results obtained in these studies overcome disadvantages of other methods applied in this area, such as scanning electron microscope or the atomic force microscope. The main advantages of 3-D laser profilometry include the ease of implementation of the technique and its nondestructive nature, which does not require sample preparation (8-10). Moreover, the technique is reproducible and presents a high degree of freedom in the vertical axes (up to 1000 μm). However, when the paper surface presents a given roughness, if the pen impressions alter the paper with a depth similar to the roughness of medium, the results are not always conclusive (8). It becomes difficult in this case to distinguish which characteristics can be imputed to the pen impressions or the quality of the paper surface. This important limitation is assessed by testing different types of paper of variable quality (of different grammage and finishing) and the writing pressure. The authors will therefore assess the limits of 3-D laser profilometry technique and determine whether the method can overcome such constraints. Second, the authors will investigate the use of digital microscopy because it presents a number of advantages: it is efficient, user-friendly, and provides an objective evaluation and interpretation.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.