PILER-CR: Fast and accurate identification of CRISPR repeats

BACKGROUND: Sequencing of prokaryotic genomes has recently revealed the presence of CRISPR elements: short, highly conserved repeats separated by unique sequences of similar length. The distinctive sequence signature of CRISPR repeats can be found using general-purpose repeat- or pattern-finding software tools. However, the output of such tools is not always ideal for studying these repeats, and significant effort is sometimes needed to build additional tools and perform manual analysis of the output. RESULTS: We present PILER-CR, a program specifically designed for the identification and analysis of CRISPR repeats. The program executes rapidly, completing a 5 Mb genome in around 5 seconds on a current desktop computer. We validate the algorithm by manual curation and by comparison with published surveys of these repeats, finding that PILER-CR has both high sensitivity and high specificity. We also present a catalogue of putative CRISPR repeats identified in a comprehensive analysis of 346 prokaryotic genomes. CONCLUSION: PILER-CR is a useful tool for rapid identification and classification of CRISPR repeats. The software is donated to the public domain. Source code and a Linux binary are freely available at

Comparing orbiter and rover image-based mapping of an ancient sedimentary environment, Aeolis Palus, Gale crater, Mars

This study provides the first systematic comparison of orbital facies maps with detailed ground-based geology observations from the Mars Science Laboratory (MSL) Curiosity rover to examine the validity of geologic interpretations derived from orbital image data. Orbital facies maps were constructed for the Darwin, Cooperstown, and Kimberley waypoints visited by the Curiosity rover using High Resolution Imaging Science Experiment (HiRISE) images. These maps, which represent the most detailed orbital analysis of these areas to date, were compared with rover image-based geologic maps and stratigraphic columns derived from Curiosity's Mast Camera (Mastcam) and Mars Hand Lens Imager (MAHLI). Results show that bedrock outcrops can generally be distinguished from unconsolidated surficial deposits in high-resolution orbital images and that orbital facies mapping can be used to recognize geologic contacts between well-exposed bedrock units. However, process-based interpretations derived from orbital image mapping are difficult to infer without known regional context or observable paleogeomorphic indicators, and layer-cake models of stratigraphy derived from orbital maps oversimplify depositional relationships as revealed from a rover perspective. This study also shows that fine-scale orbital image-based mapping of current and future Mars landing sites is essential for optimizing the efficiency and science return of rover surface operations

The genome and transcriptome of Trichormus sp NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP

Accurate reconstruction of insertion-deletion histories by statistical phylogenetics

The Multiple Sequence Alignment (MSA) is a computational abstraction that represents a partial summary either of indel history, or of structural similarity. Taking the former view (indel history), it is possible to use formal automata theory to generalize the phylogenetic likelihood framework for finite substitution models (Dayhoff's probability matrices and Felsenstein's pruning algorithm) to arbitrary-length sequences. In this paper, we report results of a simulation-based benchmark of several methods for reconstruction of indel history. The methods tested include a relatively new algorithm for statistical marginalization of MSAs that sums over a stochastically-sampled ensemble of the most probable evolutionary histories. For mammalian evolutionary parameters on several different trees, the single most likely history sampled by our algorithm appears less biased than histories reconstructed by other MSA methods. The algorithm can also be used for alignment-free inference, where the MSA is explicitly summed out of the analysis. As an illustration of our method, we discuss reconstruction of the evolutionary histories of human protein-coding genes.Comment: 28 pages, 15 figures. arXiv admin note: text overlap with arXiv:1103.434

Phylogenetic Codivergence Supports Coevolution of Mimetic Heliconius Butterflies

The unpalatable and warning-patterned butterflies _Heliconius erato_ and _Heliconius melpomene_ provide the best studied example of mutualistic Müllerian mimicry, thought – but rarely demonstrated – to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and this was initially hailed as the most striking known case of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of _H. erato_ and _H. melpomene_, supporting repeated codivergence of mimetic populations. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of _H. erato_ and _H. melpomene_ occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. This evidence supports a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change

Oil Palm Research in Context: Identifying the Need for Biodiversity Assessment

Oil palm cultivation is frequently cited as a major threat to tropical biodiversity as it is centered on some of the world's most biodiverse regions. In this report, Web of Science was used to find papers on oil palm published since 1970, which were assigned to different subject categories to visualize their research focus. Recent years have seen a broadening in the scope of research, with a slight growth in publications on the environment and a dramatic increase in those on biofuel. Despite this, less than 1% of publications are related to biodiversity and species conservation. In the context of global vegetable oil markets, palm oil and soyabean account for over 60% of production but are the subject of less than 10% of research. Much more work must be done to establish the impacts of habitat conversion to oil palm plantation on biodiversity. Results from such studies are crucial for informing conservation strategies and ensuring sustainable management of plantations

Bacteremic complications of intravascular catheter tip colonization with Gram-negative micro-organisms in patients without preceding bacteremia

Although Gram-negative micro-organisms are frequently associated with catheter-related bloodstream infections, the prognostic value and clinical implication of a positive catheter tip culture with Gram-negative micro-organisms without preceding bacteremia remains unclear. We determined the outcomes of patients with intravascular catheters colonized with these micro-organisms, without preceding positive blood cultures, and identified risk factors for the development of subsequent Gram-negative bacteremia. All patients with positive intravascular catheter tip cultures with Gram-negative micro-organisms at the University Medical Center, Utrecht, The Netherlands, between 2005 and 2009, were retrospectively studied. Patients with Gram-negative bacteremia within 48 h before catheter removal were excluded. The main outcome measure was bacteremia with Gram-negative micro-organisms. Other endpoints were length of the hospital stay, in-hospital mortality, secondary complications of Gram-negative bacteremia, and duration of intensive care admission. A total of 280 catheters from 248 patients were colonized with Gram-negative micro-organisms. Sixty-seven cases were excluded because of preceding positive blood cultures, leaving 213 catheter tips from 181 patients for analysis. In 40 (19%) cases, subsequent Gram-negative bacteremia developed. In multivariate analysis, arterial catheters were independently associated with subsequent Gram-negative bacteremia (odds ratio [OR] = 5.00, 95% confidence interval [CI]: 1.20–20.92), as was selective decontamination of the digestive tract (SDD) (OR = 2.47, 95% CI: 1.07–5.69). Gram-negative bacteremia in patients who received SDD was predominantly caused by cefotaxime (part of the SDD)-resistant organisms. Mortality was significantly higher in the group with subsequent Gram-negative bacteremia (35% versus 20%, OR = 2.12, 95% CI: 1.00–4.49). Patients with a catheter tip colonized with Gram-negative micro-organisms had a high chance of subsequent Gram-negative bacteremia from any cause. This may be clinically relevant, as starting antibiotic treatment pre-emptively in high-risk patients with Gram-negative micro-organisms cultured from arterial intravenous catheters may be beneficial

Horizontal DNA transfer mechanisms of bacteria as weapons of intragenomic conflict

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell-cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing "arms race." Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic's effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell-cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated

Increased Genetic Diversity of HIV-1 Circulating in Hong Kong

HIV-1 group M strains are characterized into 9 pure subtypes and 48 circulating recombinant forms (CRFs). Recent studies have identified the presence of new HIV-1 recombinants in Hong Kong and their complexity continues to increase. This study aims to characterize the HIV-1 genetic diversity in Hong Kong. Phylogenetic analyses were performed by using HIV-1 pol sequences including protease and partial reverse transcriptase isolated from 1045 local patients in Hong Kong from 2003 to 2008. For the pol sequences with unassigned genotype, the evidence of recombination was determined by using sliding-window based bootscan plots and their env C2V3 region were also sequenced. Epidemiological background of these patients was further collected. The pol phylogenetic analyses highlighted the extent of HIV-1 genetic diversity in Hong Kong. Subtype B (450/1045; 43.1%) and CRF01_AE (469/1045; 44.9%) variants were clearly predominant. Other genotypes (126/1045; 12.1%) including 3 defined subtypes, 10 CRFs, 1 unassigned subtype and 33 recombinants with 11 different mosaic patterns were observed. Recombinants of subtype B and CRF01_AE were mainly found among local Chinese MSM throughout 2004 to 2008, while the CRF02_AG and subtype G recombinants were circulating among non-Chinese Asian population in Hong Kong through heterosexual transmission starting from 2008. Our study demonstrated the complex recombination of HIV-1 in Hong Kong and the need in developing surveillance system for tracking the distribution of new HIV-1 genetic variants

Virulence related sequences: insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background: Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains. Results: Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection. Conclusion: The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content