Comparison of DNA sequences with protein sequences

Integration of Ixodes ricinus genome sequencing with transcriptome and proteome annotation of the naïve midgut

Background

In Europe, Ixodes ricinus ticks are the most important vectors of diseases threatening humans, livestock, wildlife and companion animals. Nevertheless, genomic sequence information is missing and functional annotation of transcripts and proteins is limited. This lack of information is restricting studies of the vector and its interactions with pathogens and hosts. Here we present and integrate the first analysis of the I. ricinus genome with the transcriptome and proteome of the unfed I. ricinus midgut.

Methods

Whole genome sequencing was performed on I. ricinus ticks and the sequences were de novo assembled. In parallel, I. ricinus ticks were dissected and the midgut transcriptome sequenced. Both datasets were integrated by transcript discovery analysis to identify putative genes and genome contigs were screened for homology. An alignment-based and a motif-search-based approach were combined for the annotation of the midgut transcriptome. Additionally, midgut proteins were identified and annotated by mass spectrometry with public databases and the in-house built transcriptome database as references and results were cross-validated.

Results

The de novo assembly of 1 billion DNA sequences to a reference genome of 393 Mb length provides an unprecedented insight into the I. ricinus genome. A homology search revealed sequences in the assembled genome contigs homologous to 89 % of the I. scapularis genome scaffolds indicating coverage of most genome regions. We identified moreover 6,415 putative genes. More than 10,000 transcripts from naïve midgut were annotated with respect of predicted function and/or cellular localization. By combining an alignment-based with a motif-search-based annotation approach, we doubled the number of annotations throughout all functional categories. In addition, 574 gel spots were significantly identified by mass spectrometry (p < 0.05) and 285 distinct proteins expressed in the naïve midgut were annotated functionally and/or for cellular localization. Our systems approach reveals a midgut metabolism of the unfed tick that is prepared to sense and process an anticipated blood meal.

Conclusions

This multiple-omics study vastly extends the publicly available DNA and RNA databases for I. ricinus, paving the way for further in-depth analysis of the most important European disease vector and its interactions with pathogens and hosts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1981-7) contains supplementary material, which is available to authorized users.

Characterization of the complete mitochondrial genome of Conus tribblei Walls, 1977

The genus Conus sensu lato consists of 500-700 species. However, the mitochondrial genomes of only few species have been fully sequenced and reported so far. In this study, the complete mitochondrial genome of Conus tribblei, a member of the poorly known subgenus Splinoconus is sequenced with the mean coverage of 604×. The mitochondrial genome is 15 570 bp long and consists of genes encoding for 13 respiratory chain proteins, 22 tRNA and 2 rRNA. The gene organization is highly conserved among the Conus species. The longest intergenic region between tRNA-Phe and cytochrome c oxidase subunit III (cox3), which in C. tribblei is 169 bp long and contains a 112 bp long segment of inverted repeat, represents the putative control region. The control regions of Conus species exhibited variability in the length and position of the inverted repeats. Therefore, this region may have the potential to be used as a genetic marker for species discrimination.

Structural features of conopeptide genes inferred from partial sequences of the Conus tribblei genome

The evolvability of venom components (in particular, the gene-encoded peptide toxins) in venomous species serves as an adaptive strategy allowing them to target new prey types or respond to changes in the prey field. The structure, organization, and expression of the venom peptide genes may provide insights into the molecular mechanisms that drive the evolution of such genes. Conus is a particularly interesting group given the high chemical diversity of their venom peptides, and the rapid evolution of the conopeptide-encoding genes. Conus genomes, however, are large and characterized by a high proportion of repetitive sequences. As a result, the structure and organization of conopeptide genes have remained poorly known. In this study, a survey of the genome of Conus tribblei was undertaken to address this gap. A partial assembly of C. tribblei genome was generated; the assembly, though consisting of a large number of fragments, accounted for 2160.5 Mb of sequence. A large number of repetitive genomic elements consisting of 642.6 Mb of retrotransposable elements, simple repeats, and novel interspersed repeats were observed. We characterized the structural organization and distribution of conotoxin genes in the genome. A significant number of conopeptide genes (estimated to be between 148 and 193) belonging to different superfamilies with complete or nearly complete exon regions were observed, ~60 % of which were expressed. The unexpressed conopeptide genes represent hidden but significant conotoxin diversity. The conotoxin genes also differed in the frequency and length of the introns. The interruption of exons by long introns in the conopeptide genes and the presence of repeats in the introns may indicate the importance of introns in facilitating recombination, evolution and diversification of conotoxins. These findings advance our understanding of the structural framework that promotes the gene-level molecular evolution of venom peptides.

Genome Sequence of a Novel Iflavirus from mRNA Sequencing of the Pupa of Bombyx mori Inoculated with Cordyceps militaris

We discovered a novel iflavirus from the transcriptome of the Bombyx mori pupa inoculated with the insect-pathogenic fungus Cordyceps militaris. The assembled iflavirus genome has 10,119 nucleotides, with a 3'-polyadenylated tail, and it encodes a polyprotein composed of 3,004 amino acids.

Proteogenomics Dashboard for the Human Proteome Project

dasHPPboard is a novel proteomics-based dashboard that collects and reports the experiments produced by the Spanish Human Proteome Project consortium (SpHPP) and aims to help HPP to map the entire human proteome. We have followed the strategy of analog genomics projects like the Encyclopedia of DNA Elements (ENCODE), which provides a vast amount of data on human cell lines experiments. The dashboard includes results of shotgun and selected reaction monitoring proteomics experiments, post-translational modifications information, as well as proteogenomics studies. We have also processed the transcriptomics data from the ENCODE and Human Body Map (HBM) projects for the identification of specific gene expression patterns in different cell lines and tissues, taking special interest in those genes having little proteomic evidence available (missing proteins). Peptide databases have been built using single nucleotide variants and novel junctions derived from RNA-Seq data that can be used in search engines for sample-specific protein identifications on the same cell lines or tissues. The dasHPPboard has been designed as a tool that can be used to share and visualize a combination of proteomic and transcriptomic data, providing at the same time easy access to resources for proteogenomics analyses. The dasHPPboard can be freely accessed at: http://sphppdashboard.cnb.csic.es.

Comparative analysis of multiple inducible phages from Mannheimia haemolytica

BACKGROUND

Mannheimia haemolytica is a commensal bacterium that resides in the upper respiratory tract of cattle that can play a role in bovine respiratory disease. Prophages are common in the M. haemolytica genome and contribute significantly to host diversity. The objective of this research was to undertake comparative genomic analysis of phages induced from strains of M. haemolytica serotype A1 (535A and 2256A), A2 (587A and 1127A) and A6 (1152A and 3927A).

RESULTS

Overall, four P2-like (535AP1, 587AP1, 1127AP1 and 2256AP1; genomes: 34.9-35.7 kb; G+C content: 41.5-42.1 %; genes: 51-53 coding sequences, CDSs), four λ-like (535AP2, 587AP2, 1152AP2 and 3927AP1; genomes: 48.6-52.1 kb; 41.1-41.4 % mol G+C; genes: 77-83 CDSs and 2 tRNAs) and one Mu-like (3927AP2; genome: 33.8 kb; 43.1 % mol G+C; encoding 50 CDSs) phages were identified. All P2-like phages are collinear with the temperate phage φMhaA1-PHL101 with 535AP1, 2256AP1 and 1152AP1 being most closely related, followed by 587AP1 and 1127AP1. Lambdoid phages are not collinear with any other known λ-type phages, with 587AP2 being distinct from 535AP2, 3927AP1 and 1152AP2. All λ-like phages contain genes encoding a toxin-antitoxin (TA) system and cell-associated haemolysin XhlA. The Mu-like phage induced from 3927A is closely related to the phage remnant φMhaMu2 from M. haemolytica PHL21, with similar Mu-like phages existing in the genomes of M. haemolytica 535A and 587A.

CONCLUSIONS

This is among the first reports of both λ- and Mu-type phages being induced from M. haemolytica. Compared to phages induced from commensal strains of M. haemolytica serotype A2, those induced from the more virulent A1 and A6 serotypes are more closely related. Moreover, when P2-, λ- and Mu-like phages co-existed in the M. haemolytica genome, only P2- and λ-like phages were detected upon induction, suggesting that Mu-type phages may be more resistant to induction. Toxin-antitoxin gene cassettes in λ-like phages may contribute to their genomic persistence or the establishment of persister subpopulations of M. haemolytica. Further work is required to determine if the cell-associated haemolysin XhlA encoded by λ-like phages contributes to the pathogenicity and ecological fitness of M. haemolytica.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel

BACKGROUND

Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels.

RESULTS

We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle.

CONCLUSIONS

The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Transcriptomic Profiling and H3K27me3 Distribution Reveal Both Demethylase-Dependent and Independent Regulation of Developmental Gene Transcription in Cell Differentiation

The removal of histone H3 trimethylation at lysine residue 27 (H3K27me3) plays a critical role in the transcriptional initiation of developmental genes. The H3K27me3-specific KDM6 demethylases JMJD3 and UTX are responsible for the transcriptional initiation of various developmental genes, but some genes are expressed in a KDM6 demethylase-independent manner. To address the role of H3K27me3 in the retinoic acid (RA)-induced differentiation of the human carcinoma NCCIT cell line, we inhibited JMJD3 and UTX using the H3K27me3 demethylase inhibitor GSK-J4. The commitment of JMJD3/UTX-inhibited cells to a specific fate was delayed, and transcriptome profiling also revealed the differential expression of genes related to cell fate specification in demethylase-inactivated cells; the expression levels of RA metabolism and HOX family genes significantly decreased. We observed a weak correlation between H3K27me3 enrichment and transcriptional repression in the control and JMJD/UTX-inhibited cells, except for a few sets of developmental genes that are indispensable for cell fate specification. Taken together, these results provide the H3K27me3 landscape of a differentiating cell line and suggest that both demethylase-dependent and demethylase-independent transcriptional regulation play a role in early differentiation and developmental gene expression activated by H3K27me3 demethylation.

C-It-Loci: a knowledge database for tissue-enriched loci

MOTIVATION

Increasing evidences suggest that most of the genome is transcribed into RNAs, but many of them are not translated into proteins. All those RNAs that do not become proteins are called 'non-coding RNAs (ncRNAs)', which outnumbers protein-coding genes. Interestingly, these ncRNAs are shown to be more tissue specifically expressed than protein-coding genes. Given that tissue-specific expressions of transcripts suggest their importance in the expressed tissue, researchers are conducting biological experiments to elucidate the function of such ncRNAs. Owing greatly to the advancement of next-generation techniques, especially RNA-seq, the amount of high-throughput data are increasing rapidly. However, due to the complexity of the data as well as its high volume, it is not easy to re-analyze such data to extract tissue-specific expressions of ncRNAs from published datasets.

RESULTS

Here, we introduce a new knowledge database called 'C-It-Loci', which allows a user to screen for tissue-specific transcripts across three organisms: human, mouse and zebrafish. C-It-Loci is intuitive and easy to use to identify not only protein-coding genes but also ncRNAs from various tissues. C-It-Loci defines homology through sequence and positional conservation to allow for the extraction of species-conserved loci. C-It-Loci can be used as a starting point for further biological experiments.

AVAILABILITY AND IMPLEMENTATION

C-It-Loci is freely available online without registration at http://c-it-loci.uni-frankfurt.de.

CONTACT

uchida@med.uni-frankfurt.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNAseq Workflow in a Diagnostic Context and an Example of a User Friendly Implementation

Over recent years next generation sequencing (NGS) technologies evolved from costly tools used by very few, to a much more accessible and economically viable technology. Through this recently gained popularity, its use-cases expanded from research environments into clinical settings. But the technical know-how and infrastructure required to analyze the data remain an obstacle for a wider adoption of this technology, especially in smaller laboratories. We present GensearchNGS, a commercial DNAseq software suite distributed by Phenosystems SA. The focus of GensearchNGS is the optimal usage of already existing infrastructure, while keeping its use simple. This is achieved through the integration of existing tools in a comprehensive software environment, as well as custom algorithms developed with the restrictions of limited infrastructures in mind. This includes the possibility to connect multiple computers to speed up computing intensive parts of the analysis such as sequence alignments. We present a typical DNAseq workflow for NGS data analysis and the approach GensearchNGS takes to implement it. The presented workflow goes from raw data quality control to the final variant report. This includes features such as gene panels and the integration of online databases, like Ensembl for annotations or Cafe Variome for variant sharing.

Complete Genome Sequence of Escherichia coli O145:NM Bacteriophage vB_EcoM_AYO145A, a New Member of O1-Like Phages

Previously, bacteriophage vB_EcoM_AYO145A, which lyses Shiga toxin-producing Escherichia coli O145:NM, was classified as an O1-like virus of the Myoviridae family. Here, we report the complete genome sequence of this phage and a comparative genomic analysis with other known O1-like phages.

Comparison of the Venom Peptides and Their Expression in Closely Related Conus Species: Insights into Adaptive Post-speciation Evolution of Conus Exogenomes

Genes that encode products with exogenous targets, which comprise an organism's "exogenome," typically exhibit high rates of evolution. The genes encoding the venom peptides (conotoxins or conopeptides) in Conus sensu lato exemplify this class of genes. Their rapid diversification has been established and is believed to be linked to the high speciation rate in this genus. However, the molecular mechanisms that underlie venom peptide diversification and ultimately emergence of new species remain poorly understood. In this study, the sequences and expression levels of conotoxins from several specimens of two closely related worm-hunting species, Conus tribblei and Conus lenavati, were compared through transcriptome analysis. Majority of the identified putative conopeptides were novel, and their diversity, even in each specimen, was remarkably high suggesting a wide range of prey targets for these species. Comparison of the interspecific expression patterns of conopeptides at the superfamily level resulted in the discovery of both conserved as well as species-specific expression patterns, indicating divergence in the regulatory network affecting conotoxin gene expression. Comparison of the transcriptomes of the individual snails revealed that each specimen produces a distinct set of highly expressed conopeptides, reflecting the capability of individual snails to fine-tune the composition of their venoms. These observations reflect the role of sequence divergence and divergence in the control of expression for specific conopeptides in the evolution of the exogenome and hence venom composition in Conus.

Epithelial Cadherin Determines Resistance to Infectious Pancreatic Necrosis Virus in Atlantic Salmon

Infectious pancreatic necrosis virus (IPNV) is the cause of one of the most prevalent diseases in farmed Atlantic salmon (Salmo salar). A quantitative trait locus (QTL) has been found to be responsible for most of the genetic variation in resistance to the virus. Here we describe how a linkage disequilibrium-based test for deducing the QTL allele was developed, and how it was used to produce IPN-resistant salmon, leading to a 75% decrease in the number of IPN outbreaks in the salmon farming industry. Furthermore, we describe how whole-genome sequencing of individuals with deduced QTL genotypes was used to map the QTL down to a region containing an epithelial cadherin (cdh1) gene. In a coimmunoprecipitation assay, the Cdh1 protein was found to bind to IPNV virions, strongly indicating that the protein is part of the machinery used by the virus for internalization. Immunofluorescence revealed that the virus colocalizes with IPNV in the endosomes of homozygous susceptible individuals but not in the endosomes of homozygous resistant individuals. A putative causal single nucleotide polymorphism was found within the full-length cdh1 gene, in phase with the QTL in all observed haplotypes except one; the absence of a single, all-explaining DNA polymorphism indicates that an additional causative polymorphism may contribute to the observed QTL genotype patterns. Cdh1 has earlier been shown to be necessary for the internalization of certain bacteria and fungi, but this is the first time the protein is implicated in internalization of a virus.

The complete mitochondrial genome of Chaenogobius gulosus (Gobiidae, Perciformes) from the South Sea, Korea

The morphology-based taxonomy of Chaenogobius is so complex that some species have been transferred to other genera. To obtain basic information on the phylogenetic relationships among gobies, the complete mitochondrial DNA sequence of C. gulosus was determined using next-generation sequencing. The complete mitochondrial genome is 16,477 bp in length and consists of 22 tRNAs, two ribosomal RNAs, 13 protein-coding genes, and two non-coding regions (the light strand replication origin and the control region). The gene order and locations are similar to those in most other vertebrates. The nucleotide composition of the mitogenome contains 29.1% A, 29.6% T, 15.9% G, and 25.4% C, with a slight AT bias (58.7%). The genus most related to Chaenogobius was Luciogobius, as revealed by phylogenetic relationships derived using a maximum likelihood tree. This mitogenome will help to elucidate the morphological systematic complexity and phylogenetic structure of Chaenogobius and related species.

Characterization of the kidney transcriptome of the long-haired mouse Abrothrix hirta (Rodentia, Sigmodontinae) and comparison with that of the olive mouse A. olivacea

To understand how small mammals cope with the challenge of water homeostasis is a matter of interest for ecologists and evolutionary biologists. Here we take a step towards the understanding of the transcriptomic functional response of kidney using as a model the long-haired mouse (Abrothrix hirta) a species that distributes across Patagonian steppes and Austral temperate rainforests in Argentina and Chile. Specifically, we i) characterize the renal transcriptome of A. hirta, and ii) compare it with that-already available-of the co-generic and co-distributed A. olivacea. Renal mRNA transcripts from 16 specimens of A. hirta from natural populations were analyzed. Over 500 million Illumina paired-end reads were assembled de novo under two approaches, an individual assembly for each specimen, and a single in-silico normalized joint assembly including all reads from all specimens. The total number of annotated genes was similar with both strategies: an average of 14,956 in individual assemblies and 14,410 in the joint assembly. Overall, 15,463 distinct genes express in the kidney of A. hirta. Transcriptomes of A. hirta and A. olivacea were similar in terms of gene abundance and composition: 95.6% of the genes of A. hirta were also found in A. olivacea making their functional profiles also similar. However, differences in the transcriptome of these two species were observed in the set of highly expressed genes, in terms of private genes for each species and the functional profiles of highly expressed genes. As part of the novel transcriptome characterization, we provide distinct gene lists with their functional annotation that would constitute the basis for further research on these or any other species of the subfamily Sigmodontinae, which includes about 400 living species distributed from Tierra del Fuego to southern United States.

Aquatic adaptation and the evolution of smell and taste in whales

INTRODUCTION

While olfaction is one of the most important senses in most terrestrial mammals, it is absent in modern toothed whales (Odontoceti, Cetacea). Furthermore, behavioral evidence suggests that gustation is very limited. In contrast, their aquatic sistergroup, baleen whales (Mysticeti) retain small but functional olfactory organs, and nothing is known about their gustation. It is difficult to investigate mysticete chemosensory abilities because experiments in a controlled setting are impossible.

RESULTS

Here, we use the functional regionalization of the olfactory bulb (OB) to identify the loss of specific olfactory functions in mysticetes. We provide the whole-genome sequence of a mysticete and show that mysticetes lack the dorsal domain of the OB, an area known to induce innate avoidance behavior against odors of predators and spoiled foods. Genomic and fossil data suggest that mysticetes lost the dorsal domain of the OB before the Odontoceti-Mysticeti split. Furthermore, we found that all modern cetaceans are revealed to have lost the functional taste receptors.

CONCLUSION

These results strongly indicate that profound changes in the chemosensory capabilities had occurred in the cetacean lineage during the period when ancestral whales migrated from land to water.

Genomic and phenotypic insights into the ecology of Arthrobacter from Antarctic soils

BACKGROUND

Members of the bacterial genus Arthrobacter are both readily cultured and commonly identified in Antarctic soil communities. Currently, relatively little is known about the physiological traits that allow these bacteria to survive in the harsh Antarctic soil environment. The aim of this study is to investigate if Antarctic strains of Arthrobacter owe their resilience to substantial genomic changes compared to Arthrobacter spp. isolated from temperate soil environments.

RESULTS

Quantitative PCR-based analysis revealed that up to 4% of the soil bacterial communities were comprised of Arthrobacter spp. at four locations in the Ross Sea Region. Genome analysis of the seven Antarctic Arthrobacter isolates revealed several features that are commonly observed in psychrophilic/psychrotolerant bacteria. These include genes primarily associated with sigma factors, signal transduction pathways, the carotenoid biosynthesis pathway and genes induced by cold-shock, oxidative and osmotic stresses. However, these genes were also identified in genomes of seven temperate Arthrobacter spp., suggesting that these mechanisms are beneficial for growth and survival in a range of soil environments. Phenotypic characterisation revealed that Antarctic Arthrobacter isolates demonstrate significantly lower metabolic versatility and a narrower salinity tolerance range compared to temperate Arthrobacter species. Comparative analyses also revealed fewer protein-coding sequences and a significant decrease in genes associated with transcription and carbohydrate transport and metabolism in four of the seven Antarctic Arthrobacter isolates. Notwithstanding genome incompleteness, these differences together with the decreased metabolic versatility are indicative of genome content scaling.

CONCLUSIONS

The genomes of the seven Antarctic Arthrobacter isolates contained several features that may be beneficial for growth and survival in the Antarctic soil environment, although these features were not unique to the Antarctic isolates. These genome sequences allow further investigations into the expression of physiological traits that enable survival under extreme conditions and, more importantly, into the ability of these bacteria to respond to future perturbations including climate change and human impacts.

High conopeptide diversity in Conus tribblei revealed through analysis of venom duct transcriptome using two high-throughput sequencing platforms

The venom of each species of Conus contains different kinds of pharmacologically active peptides which are mostly unique to that species. Collectively, the ~500-700 species of Conus produce a large number of these peptides, perhaps exceeding 140,000 different types in total. To date, however, only a small fraction of this diversity has been characterized via transcriptome sequencing. In addition, the sampling of this chemical diversity has not been uniform across the different lineages in the genus. In this study, we used high-throughput transcriptome sequencing approach to further investigate the diversity of Conus venom peptides. We chose a species, Conus tribblei, as a representative of a poorly studied clade of Conus. Using the Roche 454 and Illumina platforms, we discovered 136 unique and novel putative conopeptides belonging to 30 known gene superfamilies and 6 new conopeptide groups, the greatest diversity so far observed from a transcriptome. Most of the identified peptides exhibited divergence from the known conopeptides, and some contained cysteine frameworks observed for the first time in cone snails. In addition, several enzymes involved in posttranslational modification of conopeptides and also some proteins involved in efficient delivery of the conopeptides to prey were identified as well. Interestingly, a number of conopeptides highly similar to the conopeptides identified in a phylogenetically distant species, the generalist feeder Conus californicus, were observed. The high diversity of conopeptides and the presence of conopeptides similar to those in C. californicus suggest that C. tribblei may have a broad range of prey preferences.

Methanobacterium enables high rate electricity-driven autotrophic sulfate reduction

The autotrophic reduction of sulfate can be sustained with a cathode as the only electron donor in bioelectrochemical systems (BES).

Next-Generation Sequencing: A Review of Technologies and Tools for Wound Microbiome Research

Significance: The colonization of wounds by specific microbes or communities of microbes may delay healing and/or lead to infection-related complication. Studies of wound-associated microbial communities (microbiomes) to date have primarily relied upon culture-based methods, which are known to have extreme biases and are not reliable for the characterization of microbiomes. Biofilms are very resistant to culture and are therefore especially difficult to study with techniques that remain standard in clinical settings. Recent Advances: Culture-independent approaches employing next-generation DNA sequencing have provided researchers and clinicians a window into wound-associated microbiomes that could not be achieved before and has begun to transform our view of wound-associated biodiversity. Within the past decade, many platforms have arisen for performing this type of sequencing, with various types of applications for microbiome research being possible on each. Critical Issues: Wound care incorporating knowledge of microbiomes gained from next-generation sequencing could guide clinical management and treatments. The purpose of this review is to outline the current platforms, their applications, and the steps necessary to undertake microbiome studies using next-generation sequencing. Future Directions: As DNA sequencing technology progresses, platforms will continue to produce longer reads and more reads per run at lower costs. A major future challenge is to implement these technologies in clinical settings for more precise and rapid identification of wound bioburden.

High opsin diversity in a non-visual infaunal brittle star

BACKGROUND

In metazoans, opsins are photosensitive proteins involved in both vision and non-visual photoreception. Echinoderms have no well-defined eyes but several opsin genes were found in the purple sea urchin (Strongylocentrotus purpuratus) genome. Molecular data are lacking for other echinoderm classes although many species are known to be light sensitive.

RESULTS

In this study focused on the European brittle star Amphiura filiformis, we first highlighted a blue-green light sensitivity using a behavioural approach. We then identified 13 new putative opsin genes against eight bona fide opsin genes in the genome of S. purpuratus. Six opsins were included in the rhabdomeric opsin group (r-opsins). In addition, one putative ciliary opsin (c-opsin), showing high similarity with the c-opsin of S. purpuratus (Sp-opsin 1), one Go opsin similar to Sp-opsins 3.1 and 3.2, two basal-branch opsins similar to Sp-opsins 2 and 5, and two neuropsins similar to Sp-opsin 8, were identified. Finally, two sequences from one putative RGR opsin similar to Sp-opsin 7 were also detected. Adult arm transcriptome analysis pinpointed opsin mRNAs corresponding to one r-opsin, one neuropsin and the homologue of Sp-opsin 2. Opsin phylogeny was determined by maximum likelihood and Bayesian analyses. Using antibodies designed against c- and r-opsins from S. purpuratus, we detected putative photoreceptor cells mainly in spines and tube feet of A. filiformis, respectively. The r-opsin expression pattern is similar to the one reported in S. purpuratus with cells labelled at the tip and at the base of the tube feet. In addition, r-opsin positive cells were also identified in the radial nerve of the arm. C-opsins positive cells, expressed in pedicellariae, spines, tube feet and epidermis in S. purpuratus were observed at the level of the spine stroma in the brittle star.

CONCLUSION

Light perception in A. filiformis seems to be mediated by opsins (c- and r-) in, at least, spines, tube feet and in the radial nerve cord. Other non-visual opsin types could participate to the light perception process indicating a complex expression pattern of opsins in this infaunal brittle star.

Genome sequences of parvovirus b19 reference strains

We report here the sequences of two reference strains of parvovirus B19 (B19V) used for quantitation of B19V DNA. One reference strain has been established by the World Health Organization (WHO) and the other by the European Pharmacopeia (Ph. Eur.) and belong to B19V genotype 1a1 and 1a2, respectively.

Genome-based comparative analyses of Antarctic and temperate species of Paenibacillus

Antarctic soils represent a unique environment characterised by extremes of temperature, salinity, elevated UV radiation, low nutrient and low water content. Despite the harshness of this environment, members of 15 bacterial phyla have been identified in soils of the Ross Sea Region (RSR). However, the survival mechanisms and ecological roles of these phyla are largely unknown. The aim of this study was to investigate whether strains of Paenibacillus darwinianus owe their resilience to substantial genomic changes. For this, genome-based comparative analyses were performed on three P. darwinianus strains, isolated from gamma-irradiated RSR soils, together with nine temperate, soil-dwelling Paenibacillus spp. The genome of each strain was sequenced to over 1,000-fold coverage, then assembled into contigs totalling approximately 3 Mbp per genome. Based on the occurrence of essential, single-copy genes, genome completeness was estimated at approximately 88%. Genome analysis revealed between 3,043-3,091 protein-coding sequences (CDSs), primarily associated with two-component systems, sigma factors, transporters, sporulation and genes induced by cold-shock, oxidative and osmotic stresses. These comparative analyses provide an insight into the metabolic potential of P. darwinianus, revealing potential adaptive mechanisms for survival in Antarctic soils. However, a large proportion of these mechanisms were also identified in temperate Paenibacillus spp., suggesting that these mechanisms are beneficial for growth and survival in a range of soil environments. These analyses have also revealed that the P. darwinianus genomes contain significantly fewer CDSs and have a lower paralogous content. Notwithstanding the incompleteness of the assemblies, the large differences in genome sizes, determined by the number of genes in paralogous clusters and the CDS content, are indicative of genome content scaling. Finally, these sequences are a resource for further investigations into the expression of physiological attributes that enable survival under extreme conditions and selection processes that affect prokaryotic genome evolution.

Detection of variants of the pRAS3, pAB5S9, and pSN254 plasmids in Aeromonas salmonicida subsp. salmonicida: multidrug resistance, interspecies exchanges, and plasmid reshaping

The ubiquitous water-borne Gram-negative bacterium Aeromonas salmonicida subsp. salmonicida is the causative agent of furunculosis, a worldwide disease in fish farms. Plasmids carrying antibiotic resistance genes have already been described for this bacterium. The aim of the present study was to identify and characterize additional multidrug resistance plasmids in A. salmonicida subsp. salmonicida. We sequenced the plasmids present in two multiple antibiotic-resistant isolates using high-throughput technologies. We also investigated 19 other isolates with various multidrug resistance profiles by genotyping PCR and assessed their resistance to tetracycline. We identified variants of the pAB5S9 and pSN254 plasmids that carry several antibiotic resistance genes and that have been previously reported in bacteria other than A. salmonicida subsp. salmonicida, which suggests a high level of interspecies exchange. Genotyping analyses and the antibiotic resistance profiles of the 19 other isolates support the idea that multiple versions of pAB5S9 and pSN254 exist in A. salmonicida subsp. salmonicida. We also identified variants of the pRAS3 plasmid. The present study revealed that A. salmonicida subsp. salmonicida harbors a wide variety of plasmids, which suggests that this ubiquitous bacterium may contribute to the spread of antibiotic resistance genes in the environment.

Complementary symbiont contributions to plant decomposition in a fungus-farming termite

Termites normally rely on gut symbionts to decompose organic matter but the Macrotermitinae domesticated Termitomyces fungi to produce their own food. This transition was accompanied by a shift in the composition of the gut microbiota, but the complementary roles of these bacteria in the symbiosis have remained enigmatic. We obtained high-quality annotated draft genomes of the termite Macrotermes natalensis, its Termitomyces symbiont, and gut metagenomes from workers, soldiers, and a queen. We show that members from 111 of the 128 known glycoside hydrolase families are represented in the symbiosis, that Termitomyces has the genomic capacity to handle complex carbohydrates, and that worker gut microbes primarily contribute enzymes for final digestion of oligosaccharides. This apparent division of labor is consistent with the Macrotermes gut microbes being most important during the second passage of comb material through the termite gut, after a first gut passage where the crude plant substrate is inoculated with Termitomyces asexual spores so that initial fungal growth and polysaccharide decomposition can proceed with high efficiency. Complex conversion of biomass in termite mounds thus appears to be mainly accomplished by complementary cooperation between a domesticated fungal monoculture and a specialized bacterial community. In sharp contrast, the gut microbiota of the queen had highly reduced plant decomposition potential, suggesting that mature reproductives digest fungal material provided by workers rather than plant substrate.

Mammalian keratin associated proteins (KRTAPs) subgenomes: disentangling hair diversity and adaptation to terrestrial and aquatic environments

Background

Adaptation of mammals to terrestrial life was facilitated by the unique vertebrate trait of body hair, which occurs in a range of morphological patterns. Keratin associated proteins (KRTAPs), the major structural hair shaft proteins, are largely responsible for hair variation.

Results

We exhaustively characterized the KRTAP gene family in 22 mammalian genomes, confirming the existence of 30 KRTAP subfamilies evolving at different rates with varying degrees of diversification and homogenization. Within the two major classes of KRTAPs, the high cysteine (HS) subfamily experienced strong concerted evolution, high rates of gene conversion/recombination and high GC content. In contrast, high glycine-tyrosine (HGT) KRTAPs showed evidence of positive selection and low rates of gene conversion/recombination. Species with more hair and of higher complexity tended to have more KRATP genes (gene expansion). The sloth, with long and coarse hair, had the most KRTAP genes (175 with 141 being intact). By contrast, the “hairless” dolphin had 35 KRTAPs and the highest pseudogenization rate (74% relative to the 19% mammalian average). Unique hair-related phenotypes, such as scales (armadillo) and spines (hedgehog), were correlated with changes in KRTAPs. Gene expression variation probably also influences hair diversification patterns, for example human have an identical KRTAP repertoire as apes, but much less hair.

Conclusions

We hypothesize that differences in KRTAP gene repertoire and gene expression, together with distinct rates of gene conversion/recombination, pseudogenization and positive selection, are likely responsible for micro and macro-phenotypic hair diversification among mammals in response to adaptations to ecological pressures.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-779) contains supplementary material, which is available to authorized users.

Frameshift alignment: statistics and post-genomic applications

MOTIVATION

The alignment of DNA sequences to proteins, allowing for frameshifts, is a classic method in sequence analysis. It can help identify pseudogenes (which accumulate mutations), analyze raw DNA and RNA sequence data (which may have frameshift sequencing errors), investigate ribosomal frameshifts, etc. Often, however, only ad hoc approximations or simulations are available to provide the statistical significance of a frameshift alignment score.

RESULTS

We describe a method to estimate statistical significance of frameshift alignments, similar to classic BLAST statistics. (BLAST presently does not permit its alignments to include frameshifts.) We also illustrate the continuing usefulness of frameshift alignment with two 'post-genomic' applications: (i) when finding pseudogenes within the human genome, frameshift alignments show that most anciently conserved non-coding human elements are recent pseudogenes with conserved ancestral genes; and (ii) when analyzing metagenomic DNA reads from polluted soil, frameshift alignments show that most alignable metagenomic reads contain frameshifts, suggesting that metagenomic analysis needs to use frameshift alignment to derive accurate results.

MIRPIPE: quantification of microRNAs in niche model organisms

MicroRNAs (miRNAs) represent an important class of small non-coding RNAs regulating gene expression in eukaryotes. Present algorithms typically rely on genomic data to identify miRNAs and require extensive installation procedures. Niche model organisms lacking genomic sequences cannot be analyzed by such tools. Here we introduce the MIRPIPE application enabling rapid and simple browser-based miRNA homology detection and quantification. MIRPIPE features automatic trimming of raw RNA-Seq reads originating from various sequencing instruments, processing of isomiRs and quantification of detected miRNAs versus public- or user-uploaded reference databases.

AVAILABILITY AND IMPLEMENTATION

The Web service is freely available at http://bioinformatics.mpi-bn.mpg.de. MIRPIPE was implemented in Perl and integrated into Galaxy. An offline version for local execution is also available from our Web site.

Draft Genome Sequence of Streptomyces iranensis

Streptomyces iranensis HM 35 has been shown to exhibit 72.7% DNA-DNA similarity to the important drug rapamycin (sirolimus)-producing Streptomyces rapamycinicus NRRL5491. Here, we report the genome sequence of HM 35, which represents a partially overlapping repertoire of secondary metabolite gene clusters with S. rapamycinicus, including the gene cluster for rapamycin biosynthesis.

Cystatin C: a new renal marker and its potential use in small animal medicine

The occurrence of chronic kidney disease is underestimated in both human and veterinary medicine. Glomerular filtration rate (GFR) is considered the gold standard for evaluating kidney function. However, GFR assessment is time-consuming and labor-intensive and therefore not routinely used in practice. The commonly used indirect GFR markers, serum creatinine (sCr) and urea, are not sufficiently sensitive or specific to detect early renal dysfunction. Serum cystatin C (sCysC), a proteinase inhibitor, has most of the properties required for an endogenous GFR marker. In human medicine, numerous studies have evaluated its potential use as a GFR marker in several populations. In veterinary medicine, this marker is gaining interest. The measurement is easy, which makes it an interesting parameter for clinical use. This review summarizes current knowledge about cystatin C (CysC) in humans, dogs, and cats, including its history, assays, relationship with GFR, and biological and clinical variations in both human and veterinary medicine.

Genome editing assessment using CRISPR Genome Analyzer (CRISPR-GA)

SUMMARY

Clustered regularly interspaced short palindromic repeats (CRISPR)-based technologies have revolutionized human genome engineering and opened countless possibilities to basic science, synthetic biology and gene therapy. Albeit the enormous potential of these tools, their performance is far from perfect. It is essential to perform a posterior careful analysis of the gene editing experiment. However, there are no computational tools for genome editing assessment yet, and current experimental tools lack sensitivity and flexibility. We present a platform to assess the quality of a genome editing experiment only with three mouse clicks. The method evaluates next-generation data to quantify and characterize insertions, deletions and homologous recombination. CRISPR Genome Analyzer provides a report for the locus selected, which includes a quantification of the edited site and the analysis of the different alterations detected. The platform maps the reads, estimates and locates insertions and deletions, computes the allele replacement efficiency and provides a report integrating all the information.

AVAILABILITY AND IMPLEMENTATION

CRISPR-GA Web is available at http://crispr-ga.net. Documentation on CRISPR-GA instructions can be found at http://crispr-ga.net/documentation.html

CONTACT

mguell@genetics.med.harvard.edu.

RNA-Seq transcriptome profiling identifies CRISPLD2 as a glucocorticoid responsive gene that modulates cytokine function in airway smooth muscle cells

Asthma is a chronic inflammatory respiratory disease that affects over 300 million people worldwide. Glucocorticoids are a mainstay therapy for asthma because they exert anti-inflammatory effects in multiple lung tissues, including the airway smooth muscle (ASM). However, the mechanism by which glucocorticoids suppress inflammation in ASM remains poorly understood. Using RNA-Seq, a high-throughput sequencing method, we characterized transcriptomic changes in four primary human ASM cell lines that were treated with dexamethasone--a potent synthetic glucocorticoid (1 µM for 18 hours). Based on a Benjamini-Hochberg corrected p-value <0.05, we identified 316 differentially expressed genes, including both well known (DUSP1, KLF15, PER1, TSC22D3) and less investigated (C7, CCDC69, CRISPLD2) glucocorticoid-responsive genes. CRISPLD2, which encodes a secreted protein previously implicated in lung development and endotoxin regulation, was found to have SNPs that were moderately associated with inhaled corticosteroid resistance and bronchodilator response among asthma patients in two previously conducted genome-wide association studies. Quantitative RT-PCR and Western blotting showed that dexamethasone treatment significantly increased CRISPLD2 mRNA and protein expression in ASM cells. CRISPLD2 expression was also induced by the inflammatory cytokine IL1β, and small interfering RNA-mediated knockdown of CRISPLD2 further increased IL1β-induced expression of IL6 and IL8. Our findings offer a comprehensive view of the effect of a glucocorticoid on the ASM transcriptome and identify CRISPLD2 as an asthma pharmacogenetics candidate gene that regulates anti-inflammatory effects of glucocorticoids in the ASM.

Myocyte-mediated arginase expression controls hyperargininemia but not hyperammonemia in arginase-deficient mice

Human arginase deficiency is characterized by hyperargininemia and infrequent episodes of hyperammonemia that cause neurological impairment and growth retardation. We previously developed a neonatal mouse adeno-associated viral vector (AAV) rh10-mediated therapeutic approach with arginase expressed by a chicken β-actin promoter that controlled plasma ammonia and arginine, but hepatic arginase declined rapidly. This study tested a codon-optimized arginase cDNA and compared the chicken β-actin promoter to liver- and muscle-specific promoters. ARG1(-/-) mice treated with AAVrh10 carrying the liver-specific promoter also exhibited long-term survival and declining hepatic arginase accompanied by the loss of AAV episomes during subsequent liver growth. Although arginase expression in striated muscle was not expected to counteract hyperammonemia, due to muscle's lack of other urea cycle enzymes, we hypothesized that the postmitotic phenotype in muscle would allow vector genomes to persist, and hence contribute to decreased plasma arginine. As anticipated, ARG1(-/-) neonatal mice treated with AAVrh10 carrying a modified creatine kinase-based muscle-specific promoter did not survive longer than controls; however, their plasma arginine levels remained normal when animals were hyperammonemic. These data imply that plasma arginine can be controlled in arginase deficiency by muscle-specific expression, thus suggesting an alternative approach to utilizing the liver for treating hyperargininemia.

Analytical validation of a human particle-enhanced nephelometric assay for cystatin C measurement in feline serum and urine

BACKGROUND

In people and dogs, Cystatin C (CysC), a renal glomerular and tubular marker, seems superior to serum creatinine to estimate the glomerular filtration rate (GFR). A particle-enhanced nephelometric immunoassay is available to measure human CysC, but there are no reports in cats.

OBJECTIVE

The goal of this study was the validation of the human CysC nephelometric assay with feline serum and urine, and to perform a pilot study comparing serum and urine CysC between healthy cats and cats with chronic kidney disease (CKD).

METHODS

Western blot analysis was used to assess cross-reactivity between the polyclonal rabbit anti-human CysC antibody and feline CysC. Imprecision and linearity were determined for feline serum and urine CysC. Serum and urine CysC were measured in 10 healthy and 10 CKD cats.

RESULTS

Cross-reactivity between the polyclonal rabbit anti-human CysC antibody and feline CysC was demonstrated. Intra- and inter-assay coefficients of variation in feline serum and urine were 1.3% and 0.4%, and 12.5%, and 4.1%, respectively. Cats with CKD had a significantly higher serum CysC concentration (1.24 [0.63-2.99] vs 0.79 [0.43-1.05] mg/L; P = .02) and urine CysC/urinary Creatinine (uCr) ratio (565.6 [0-1311] vs < 0.049/uCr mg/mol; P = .005) compared with healthy cats.

CONCLUSIONS

The human nephelometric assay showed satisfactory validation results for feline CysC. Cats with CKD had a significantly higher sCysC concentration and uCysC/uCr ratio compared with healthy cats. Additional studies are necessary to evaluate CysC as an early marker of renal damage in cats.

Cloning and transcriptional expression of mouse mannosyltransferase IV/V cDNA, which is involved in the synthesis of lipid-linked oligosaccharides

We cloned the mouse mannosyltransferase IV/V gene (mALG11) from FM3A cells by a bioinformatic approach. The ORF contained 1476 bp encoding 492 amino acids. The cloned mALG11 complemented the growth defect of the Saccharomyces cerevisiae ALG11Δ mutant. In addition, we detected a variant cDNA by alternate splicing that had an additional four-nucleotide ATGC insertion at base 276 of the ORF. Consequently the variant cDNA encoded a truncated protein with 92 amino acids, lacking the glycosyltransferase group-1 domain. The variant cDNA occurs in many mouse strains according to EST database searches. Moreover, we detected it in FM3A cDNA, but we did not detect any such variants in the human EST database or in HeLa cDNA, although human ALG11 (hALG11) genomic DNA has the same sequence around the intron-exon boundaries as those of mALG11 genomic DNA. Hence, we concluded that there is different transcriptional control mechanism between mALG11 and hALG11.

DNA methylation and evolution of duplicate genes

The evolutionary mechanisms underlying duplicate gene maintenance and divergence remain highly debated. Epigenetic modifications, such as DNA methylation, may contribute to duplicate gene evolution by facilitating tissue-specific regulation. However, the role of epigenetic divergence on duplicate gene evolution remains little understood. Here we show, using comprehensive data across 10 diverse human tissues, that DNA methylation plays critical roles in several aspects of duplicate gene evolution. We first demonstrate that duplicate genes are initially heavily methylated, before gradually losing DNA methylation as they age. Within each pair, DNA methylation divergence between duplicate partners increases with evolutionary age. Importantly, tissue-specific DNA methylation of duplicates correlates with tissue-specific expression, implicating DNA methylation as a causative factor for functional divergence of duplicate genes. These patterns are apparent in promoters but not in gene bodies, in accord with the complex relationship between gene-body DNA methylation and transcription. Remarkably, many duplicate gene pairs exhibit consistent division of DNA methylation across multiple, divergent tissues: For the majority (73%) of duplicate gene pairs, one partner is always hypermethylated compared with the other. This is indicative of a common underlying determinant of DNA methylation. The division of DNA methylation is also consistent with their chromatin accessibility profiles. Moreover, at least two sequence motifs known to interact with the Sp1 transcription factor mark promoters of more hypomethylated duplicate partners. These results demonstrate critical roles of DNA methylation, as well as complex interaction between genome and epigenome, on duplicate gene evolution.

A survey of plant and algal genomes and transcriptomes reveals new insights into the evolution and function of the cellulose synthase superfamily

Background

Enzymes of the cellulose synthase (CesA) family and CesA-like (Csl) families are responsible for the synthesis of celluloses and hemicelluloses, and thus are of great interest to bioenergy research. We studied the occurrences and phylogenies of CesA/Csl families in diverse plants and algae by comprehensive data mining of 82 genomes and transcriptomes.

Results

We found that 1) charophytic green algae (CGA) have orthologous genes in CesA, CslC and CslD families; 2) liverwort genes are found in the CesA, CslA, CslC and CslD families; 3) The fern Pteridium aquilinum not only has orthologs in these conserved families but also in the CslB, CslH and CslE families; 4) basal angiosperms, e.g. Aristolochia fimbriata, have orthologs in these families too; 5) gymnosperms have genes forming clusters ancestral to CslB/H and to CslE/J/G respectively; 6) CslG is found in switchgrass and basal angiosperms; 7) CslJ is widely present in dicots and monocots; 8) CesA subfamilies have already diversified in ferns.

Conclusions

We speculate that: (i) ferns and horsetails might both have CslH enzymes, responsible for the synthesis of mixed-linkage glucans and (ii) CslD and similar genes might be responsible for the synthesis of mannans in CGA. Our findings led to a more detailed model of cell wall evolution and suggested that gene loss played an important role in the evolution of Csl families. We also demonstrated the usefulness of transcriptome data in the study of plant cell wall evolution and diversity.

Novel coronavirus and astrovirus in Delaware Bay shorebirds

BACKGROUND

Wild birds are an important but to some extent under-studied reservoir for emerging pathogens. We used unbiased sequencing methods for virus discovery in shorebird samples from the Delaware Bay, USA; an important feeding ground for thousands of migratory birds.

FINDINGS

Analysis of shorebird fecal samples indicated the presence of a novel astrovirus and coronavirus. A sanderling sample yielded sequences with distant homology to avian nephritis virus 1, an astrovirus associated with acute nephritis in poultry. A ruddy turnstone sample yielded sequences with homology to deltacoronaviruses.

CONCLUSIONS

Our findings highlight shorebirds as a virus reservoir and the need to closely monitor wild bird populations for the emergence of novel virus variants.

Using RNA sequencing to characterize female reproductive genes between Z and E Strains of European Corn Borer moth (Ostrinia nubilalis)

Background

Reproductive proteins often evolve rapidly and are thought to be subject to strong sexual selection, and thus may play a key role in reproductive isolation and species divergence. However, our knowledge of reproductive proteins has been largely limited to males and model organisms with sequenced genomes. With advances in sequencing technology, Lepidoptera are emerging models for studies of sexual selection and speciation. By profiling the transcriptomes of the bursa copulatrix and bursal gland from females of two incipient species of moth, we characterize reproductive genes expressed in the primary reproductive tissues of female Lepidoptera and identify candidate genes contributing to a one-way gametic incompatibility between Z and E strains of the European corn borer (Ostrinia nubilalis).

Results

Using RNA sequencing we identified transcripts from ~37,000 and ~36,000 loci that were expressed in the bursa copulatrix or the bursal gland respectively. Of bursa copulatrix genes, 8% were significantly differentially expressed compared to the female thorax, and those that were up-regulated or specific to the bursa copulatrix showed functional biases toward muscle activity and/or organization. In the bursal gland, 9% of genes were differentially expressed compared to the thorax, with many showing reproduction or gamete production functions. Of up-regulated bursal gland genes, 46% contained a transmembrane region and 16% possessed secretion signal peptides. Divergently expressed genes in the bursa copulatrix were exclusively biased toward protease-like functions and 51 proteases or protease inhibitors were divergently expressed overall.

Conclusions

This is the first comprehensive characterization of female reproductive genes in any lepidopteran system. The transcriptome of the bursa copulatrix supports its role as a muscular sac that is the primary site for disruption of the male ejaculate. We find that the bursal gland acts as a reproductive secretory body that might also interact with male ejaculate. In addition, differential expression of proteases between strains supports a potential role for these tissues in contributing to reproductive isolation. Our study provides new insight into how male ejaculate is processed by female Lepidoptera, and paves the way for future work on interactions between post-mating sexual selection and speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-189) contains supplementary material, which is available to authorized users.

BLAST and FASTA similarity searching for multiple sequence alignment

BLAST, FASTA, and other similarity searching programs seek to identify homologous proteins and DNA sequences based on excess sequence similarity. If two sequences share much more similarity than expected by chance, the simplest explanation for the excess similarity is common ancestry-homology. The most effective similarity searches compare protein sequences, rather than DNA sequences, for sequences that encode proteins, and use expectation values, rather than percent identity, to infer homology. The BLAST and FASTA packages of sequence comparison programs provide programs for comparing protein and DNA sequences to protein databases (the most sensitive searches). Protein and translated-DNA comparisons to protein databases routinely allow evolutionary look back times from 1 to 2 billion years; DNA:DNA searches are 5-10-fold less sensitive. BLAST and FASTA can be run on popular web sites, but can also be downloaded and installed on local computers. With local installation, target databases can be customized for the sequence data being characterized. With today's very large protein databases, search sensitivity can also be improved by searching smaller comprehensive databases, for example, a complete protein set from an evolutionarily neighboring model organism. By default, BLAST and FASTA use scoring strategies target for distant evolutionary relationships; for comparisons involving short domains or queries, or searches that seek relatively close homologs (e.g. mouse-human), shallower scoring matrices will be more effective. Both BLAST and FASTA provide very accurate statistical estimates, which can be used to reliably identify protein sequences that diverged more than 2 billion years ago.

N6-methyladenosine-dependent regulation of messenger RNA stability

N(6)-methyladenosine (m(6)A) is the most prevalent internal (non-cap) modification present in the messenger RNA of all higher eukaryotes. Although essential to cell viability and development, the exact role of m(6)A modification remains to be determined. The recent discovery of two m(6)A demethylases in mammalian cells highlighted the importance of m(6)A in basic biological functions and disease. Here we show that m(6)A is selectively recognized by the human YTH domain family 2 (YTHDF2) 'reader' protein to regulate mRNA degradation. We identified over 3,000 cellular RNA targets of YTHDF2, most of which are mRNAs, but which also include non-coding RNAs, with a conserved core motif of G(m(6)A)C. We further establish the role of YTHDF2 in RNA metabolism, showing that binding of YTHDF2 results in the localization of bound mRNA from the translatable pool to mRNA decay sites, such as processing bodies. The carboxy-terminal domain of YTHDF2 selectively binds to m(6)A-containing mRNA, whereas the amino-terminal domain is responsible for the localization of the YTHDF2-mRNA complex to cellular RNA decay sites. Our results indicate that the dynamic m(6)A modification is recognized by selectively binding proteins to affect the translation status and lifetime of mRNA.

Mining of haplotype-based expressed sequence tag single nucleotide polymorphisms in citrus

Background

Single nucleotide polymorphisms (SNPs), the most abundant variations in a genome, have been widely used in various studies. Detection and characterization of citrus haplotype-based expressed sequence tag (EST) SNPs will greatly facilitate further utilization of these gene-based resources.

Results

In this paper, haplotype-based SNPs were mined out of publicly available citrus expressed sequence tags (ESTs) from different citrus cultivars (genotypes) individually and collectively for comparison. There were a total of 567,297 ESTs belonging to 27 cultivars in varying numbers and consequentially yielding different numbers of haplotype-based quality SNPs. Sweet orange (SO) had the most (213,830) ESTs, generating 11,182 quality SNPs in 3,327 out of 4,228 usable contigs. Summed from all the individually mining results, a total of 25,417 quality SNPs were discovered – 15,010 (59.1%) were transitions (AG and CT), 9,114 (35.9%) were transversions (AC, GT, CG, and AT), and 1,293 (5.0%) were insertion/deletions (indels). A vast majority of SNP-containing contigs consisted of only 2 haplotypes, as expected, but the percentages of 2 haplotype contigs varied widely in these citrus cultivars. BLAST of the 25,417 25-mer SNP oligos to the Clementine reference genome scaffolds revealed 2,947 SNPs had “no hits found”, 19,943 had 1 unique hit / alignment, 1,571 had one hit and 2+ alignments per hit, and 956 had 2+ hits and 1+ alignment per hit. Of the total 24,293 scaffold hits, 23,955 (98.6%) were on the main scaffolds 1 to 9, and only 338 were on 87 minor scaffolds. Most alignments had 100% (25/25) or 96% (24/25) nucleotide identities, accounting for 93% of all the alignments. Considering almost all the nucleotide discrepancies in the 24/25 alignments were at the SNP sites, it served well as in silico validation of these SNPs, in addition to and consistent with the rate (81%) validated by sequencing and SNaPshot assay.

Conclusions

High-quality EST-SNPs from different citrus genotypes were detected, and compared to estimate the heterozygosity of each genome. All the SNP oligo sequences were aligned with the Clementine citrus genome to determine their distribution and uniqueness and for in silico validation, in addition to SNaPshot and sequencing validation of selected SNPs.

Genomically recoded organisms expand biological functions

We describe the construction and characterization of a genomically recoded organism (GRO). We replaced all known UAG stop codons in Escherichia coli MG1655 with synonymous UAA codons, which permitted the deletion of release factor 1 and reassignment of UAG translation function. This GRO exhibited improved properties for incorporation of nonstandard amino acids that expand the chemical diversity of proteins in vivo. The GRO also exhibited increased resistance to T7 bacteriophage, demonstrating that new genetic codes could enable increased viral resistance.

Identification of novel vaccine candidates against multidrug-resistant Acinetobacter baumannii

Acinetobacter baumannii is an emerging opportunistic bacterium associated with nosocomial infections in intensive care units. The alarming increase in infections caused by A. baumannii is strongly associated with enhanced resistance to antibiotics, in particular carbapenems. This, together with the lack of a licensed vaccine, has translated into significant economic, logistic and health impacts to health care facilities. In this study, we combined reverse vaccinology and proteomics to identify surface-exposed and secreted antigens from A. baumannii. Using in silico prediction tools and comparative genome analysis in combination with in vitro proteomic approaches, we identified 42 antigens that could be used as potential vaccine targets. Considering the paucity of effective antibiotics available to treat multidrug-resistant A. baumannii infections, these vaccine targets may serve as a framework for the development of a broadly protective multi-component vaccine, an outcome that would have a major impact on the burden of A. baumannii infections in intensive care units across the globe.

The effects of pre-weaning undernutrition on the expression levels of free radical deactivating enzymes in the mouse brain

A mild degree of undernutrition brought about by restricting the amount of food in the diet is known to alter the life span of an animal. It has been hypothesised that this may be related to the effects of undernutrition on an animals anti-oxidant defense system. We have therefore, used real-time PCR (rt-PCR) techniques to determine the levels of mRNA expression for manganese superoxide dismutase (MnSOD), copper/zinc superoxide dismutase (Cu/ZnSOD), glutathione peroxidase 1 (GPx 1) and catalase in the brains of Quackenbush mice undernourished from conception until 21-post-natal days of age. It was found that 21- and 61-day-old undernourished mice had a deficit in the expression of Cu/ZnSOD in both the cerebellum and forebrain regions compared to age-matched controls. The expression of MnSOD was found to be greater in the cerebellum, but not the forebrain region, of 21-day-old undernourished mice. There were no significant differences in the expression of GPx 1 and catalase between control and undernourished or previously undernourished mice. Our results confirm that undernutrition during the early life of a mouse may disrupt some of the enzymes involved in the anti-oxidant defense systems.

Conservation of meningococcal antigens in the genus Neisseria

Neisseria meningitidis, one of the major causes of bacterial meningitis and sepsis, is a member of the genus Neisseria, which includes species that colonize the mucosae of many animals. Three meningococcal proteins, factor H-binding protein (fHbp), neisserial heparin-binding antigen (NHBA), and N. meningitidis adhesin A (NadA), have been described as antigens protective against N. meningitidis of serogroup B, and they have been employed as vaccine components in preclinical and clinical studies. In the vaccine formulation, fHbp and NHBA were fused to the GNA2091 and GNA1030 proteins, respectively, to enhance protein stability and immunogenicity. To determine the possible impact of vaccination on commensal neisseriae, we determined the presence, distribution, and conservation of these antigens in the available genome sequences of the genus Neisseria, finding that fHbp, NHBA, and NadA were conserved only in species colonizing humans, while GNA1030 and GNA2091 were conserved in many human and nonhuman neisseriae. Sequence analysis showed that homologous recombination contributed to shape the evolution and distribution of both NHBA and fHbp, three major variants of which have been defined. fHbp variant 3 was probably the ancestral form of meningococcal fHbp, while fHbp variant 1 from N. cinerea was introduced into N. meningitidis by a recombination event. fHbp variant 2 was the result of a recombination event inserting a stretch of 483 bp from variant 1 into the variant 3 background. These data indicate that a high rate of exchange of genetic material between neisseriae that colonize the human upper respiratory tract exists.

The upper respiratory tract of healthy individuals is a complex ecosystem colonized by many bacterial species. Among these, there are representatives of the genus Neisseria, including Neisseria meningitidis, a major cause of bacterial meningitis and sepsis. Given the close relationship between commensal and pathogenic species, a protein-based vaccine against N. meningitidis has the potential to impact the other commensal species of Neisseria. For this reason, we have studied the distribution and evolutionary history of the antigen components of a recombinant vaccine, 4CMenB, that recently received approval in Europe under the commercial name of Bexsero®. We found that fHbp, NHBA, and NadA can be found in some of the human commensal species and that the evolution of these antigens has been essentially shaped by the high rate of genetic exchange that occurs between strains of neisseriae that cocolonize the same environment.