Tombusvirus p19 Captures RNase III-Cleaved Double-Stranded RNAs Formed by Overlapping Sense and Antisense Transcripts in Escherichia coli.

Antisense transcription is widespread in bacteria. By base pairing with overlapping sense RNAs, antisense RNAs (asRNA) can form double-stranded RNAs (dsRNA), which are cleaved by RNase III, a dsRNA endoribonuclease. The ectopic expression of plant Tombusvirus p19 in Escherichia coli stabilizes ∼21-nucleotide (nt) dsRNA RNase III decay intermediates, which enabled us to characterize otherwise highly unstable asRNA by deep sequencing of p19-captured dsRNA. RNase III-produced small dsRNA were formed at most bacterial genes in the bacterial genome and in a plasmid. We classified the types of asRNA in genomic clusters producing the most abundant p19-captured dsRNA and confirmed RNase III regulation of asRNA and sense RNA decay at three type I toxin-antitoxin loci and at a coding gene, rsd Furthermore, we provide potential evidence for the RNase III-dependent regulation of CspD protein by asRNA. The analysis of p19-captured dsRNA revealed an RNase III sequence preference for AU-rich sequences 3 nucleotides on either side of the cleavage sites and for GC-rich sequences in the 2-nt overhangs. Unexpectedly, GC-rich sequences were enriched in the middle section of p19-captured dsRNA, suggesting some unexpected sequence bias in p19 protein binding. Nonetheless, the ectopic expression of p19 is a sensitive method for identifying antisense transcripts and RNase III cleavage sites in dsRNA formed by overlapping sense and antisense transcripts in bacteria.

A computational framework to explore large-scale biosynthetic diversity.

Genome mining has become a key technology to exploit natural product diversity. Although initially performed on a single-genome basis, the process is now being scaled up to mine entire genera, strain collections and microbiomes. However, no bioinformatic framework is currently available for effectively analyzing datasets of this size and complexity. In the present study, a streamlined computational workflow is provided, consisting of two new software tools: the 'biosynthetic gene similarity clustering and prospecting engine' (BiG-SCAPE), which facilitates fast and interactive sequence similarity network analysis of biosynthetic gene clusters and gene cluster families; and the 'core analysis of syntenic orthologues to prioritize natural product gene clusters' (CORASON), which elucidates phylogenetic relationships within and across these families. BiG-SCAPE is validated by correlating its output to metabolomic data across 363 actinobacterial strains and the discovery potential of CORASON is demonstrated by comprehensively mapping biosynthetic diversity across a range of detoxin/rimosamide-related gene cluster families, culminating in the characterization of seven detoxin analogues.

Total Lipopolysaccharide from the Human Gut Microbiome Silences Toll-Like Receptor Signaling.

Cohabitation of microbial communities with the host enables the formation of a symbiotic relationship that maintains homeostasis in the gut and beyond. One prevailing model suggests that this relationship relies on the capacity of host cells and tissues to remain tolerant to the strong immune stimulation generated by the microbiota such as the activation of Toll-like receptor 4 (TLR4) pathways by lipopolysaccharide (LPS). Indeed, gut microbial LPS is thought to be one of the most potent activators of innate immune signaling and an important mediator of the microbiome's influence on host physiology. In this study, we performed computational and experimental analyses of healthy human fecal samples to examine the TLR4 signaling capacity of the gut microbiota. These analyses revealed that an immunoinhibitory activity of LPS, conserved across the members of the order Bacteroidales and derived from an underacylated structural feature, silences TLR4 signaling for the entire consortium of organisms inhabiting the human gut. Comparative analysis of metagenomic data from the Human Microbiome Project and healthy-donor samples indicates that immune silencing via LPS is a microbe-intrinsic feature in all healthy adults. These findings challenge the current belief that robust TLR4 signaling is a feature of the microbiome and demonstrate that microbiome-derived LPS has the ability to facilitate host tolerance of gut microbes. These findings have broad implications for how we model host-microbe interactions and for our understanding of microbiome-linked disease. IMPORTANCE While the ability for humans to host a complex microbial ecosystem is an essential property of life, the mechanisms allowing for immune tolerance of such a large microbial load are not completely understood and are currently the focus of intense research. This study shows that an important proinflammatory pathway that is commonly triggered by pathogenic bacteria upon interaction with the host is, in fact, actively repressed by the bacteria of the gut microbiome, supporting the idea that beneficial microbes themselves contribute to the immune tolerance in support of homeostasis. These findings are important for two reasons. First, many currently assume that proinflammatory signaling by lipopolysaccharide is a fundamental feature of the gut flora. This assumption influences greatly how host-microbiome interactions are theoretically modeled but also how they are experimentally studied, by using robust TLR signaling conditions to simulate commensals. Second, elucidation of the mechanisms that support host-microbe tolerance is key to the development of therapeutics for both intestinal and systemic inflammatory disorders.

Indoleacrylic Acid Produced by Commensal Peptostreptococcus Species Suppresses Inflammation.

Host factors in the intestine help select for bacteria that promote health. Certain commensals can utilize mucins as an energy source, thus promoting their colonization. However, health conditions such as inflammatory bowel disease (IBD) are associated with a reduced mucus layer, potentially leading to dysbiosis associated with this disease. We characterize the capability of commensal species to cleave and transport mucin-associated monosaccharides and identify several Clostridiales members that utilize intestinal mucins. One such mucin utilizer, Peptostreptococcus russellii, reduces susceptibility to epithelial injury in mice. Several Peptostreptococcus species contain a gene cluster enabling production of the tryptophan metabolite indoleacrylic acid (IA), which promotes intestinal epithelial barrier function and mitigates inflammatory responses. Furthermore, metagenomic analysis of human stool samples reveals that the genetic capability of microbes to utilize mucins and metabolize tryptophan is diminished in IBD patients. Our data suggest that stimulating IA production could promote anti-inflammatory responses and have therapeutic benefits.

mBLAST: Keeping up with the sequencing explosion for (meta)genome analysis.

Recent advances in next-generation sequencing technologies require alignment algorithms and software that can keep pace with the heightened data production. Standard algorithms, especially protein similarity searches, represent significant bottlenecks in analysis pipelines. For metagenomic approaches in particular, it is now often necessary to search hundreds of millions of sequence reads against large databases. Here we describe mBLAST, an accelerated search algorithm for translated and/or protein alignments to large datasets based on the Basic Local Alignment Search Tool (BLAST) and retaining the high sensitivity of BLAST. The mBLAST algorithms achieve substantial speed up over the National Center for Biotechnology Information (NCBI) programs BLASTX, TBLASTX and BLASTP for large datasets, allowing analysis within reasonable timeframes on standard computer architectures. In this article, the impact of mBLAST is demonstrated with sequences originating from the microbiota of healthy humans from the Human Microbiome Project. mBLAST is designed as a plug-in replacement for BLAST for any study that involves short-read sequences and includes high-throughput analysis. The mBLAST software is freely available to academic users at www.multicorewareinc.com.

Metabolic and metagenomic outcomes from early-life pulsed antibiotic treatment.

Mammalian species have co-evolved with intestinal microbial communities that can shape development and adapt to environmental changes, including antibiotic perturbation or nutrient flux. In humans, especially children, microbiota disruption is common, yet the dynamic microbiome recovery from early-life antibiotics is still uncharacterized. Here we use a mouse model mimicking paediatric antibiotic use and find that therapeutic-dose pulsed antibiotic treatment (PAT) with a beta-lactam or macrolide alters both host and microbiota development. Early-life PAT accelerates total mass and bone growth, and causes progressive changes in gut microbiome diversity, population structure and metagenomic content, with microbiome effects dependent on the number of courses and class of antibiotic. Whereas control microbiota rapidly adapts to a change in diet, PAT slows the ecological progression, with delays lasting several months with previous macrolide exposure. This study identifies key markers of disturbance and recovery, which may help provide therapeutic targets for microbiota restoration following antibiotic treatment.

Mining novel effector proteins from the esophageal gland cells of Meloidogyne incognita.

Meloidogyne incognita is one of the most economically damaging plant pathogens in agriculture and horticulture. Identifying and characterizing the effector proteins which M. incognita secretes into its host plants during infection is an important step toward finding new ways to manage this pest. In this study, we have identified the cDNAs for 18 putative effectors (i.e., proteins that have the potential to facilitate M. incognita parasitism of host plants). These putative effectors are secretory proteins that do not contain transmembrane domains and whose genes are specifically expressed in the secretory gland cells of the nematode, indicating that they are likely secreted from the nematode through its stylet. We have determined that, in the plant cells, these putative effectors are likely to localize to the cytoplasm. Furthermore, the transcripts of many of these novel effectors are specifically upregulated during different stages of the nematode's life cycle, indicating that they function at specific stages during M. incognita parasitism. The predicted proteins showed little to no homology to known proteins from free-living nematode species, suggesting that they evolved recently to support the parasitic lifestyle. On the other hand, several of the effectors are part of gene families within the M. incognita genome as well as that of M. hapla, which points to an important role that these putative effectors are playing in both parasites. With the discovery of these putative effectors, we have increased our knowledge of the effector repertoire utilized by root-knot nematodes to infect, feed on, and reproduce on their host plants. Future studies investigating the roles that these proteins play in planta will help mitigate the effects of this damaging pest.

The prevalence of species and strains in the human microbiome: a resource for experimental efforts.

Experimental efforts to characterize the human microbiota often use bacterial strains that were chosen for historical rather than biological reasons. Here, we report an analysis of 380 whole-genome shotgun samples from 100 subjects from the NIH Human Microbiome Project. By mapping their reads to 1,751 reference genome sequences and analyzing the resulting relative strain abundance in each sample we present metrics and visualizations that can help identify strains of interest for experimentalists. We also show that approximately 14 strains of 10 species account for 80% of the mapped reads from a typical stool sample, indicating that the function of a community may not be irreducibly complex. Some of these strains account for >20% of the sequence reads in a subset of samples but are absent in others, a dichotomy that could underlie biological differences among subjects. These data should serve as an important strain selection resource for the community of researchers who take experimental approaches to studying the human microbiota.

Identification and characterization of alternative splicing in parasitic nematode transcriptomes.

BACKGROUND: Alternative splicing (AS) of mRNA is a vital mechanism for enhancing genomic complexity in eukaryotes. Spliced isoforms of the same gene can have diverse molecular and biological functions and are often differentially expressed across various tissues, times, and conditions. Thus, AS has important implications in the study of parasitic nematodes with complex life cycles. Transcriptomic datasets are available from many species, but data must be revisited with splice-aware assembly protocols to facilitate the study of AS in helminthes. METHODS: We sequenced cDNA from the model worm Caenorhabditis elegans using 454/Roche technology for use as an experimental dataset. Reads were assembled with Newbler software, invoking the cDNA option. Several combinations of parameters were tested and assembled transcripts were verified by comparison with previously reported C. elegans genes and transcript isoforms and with Illumina RNAseq data. RESULTS: Thoughtful adjustment of program parameters increased the percentage of assembled transcripts that matched known C. elegans sequences, decreased mis-assembly rates (i.e., cis- and trans-chimeras), and improved the coverage of the geneset. The optimized protocol was used to update de novo transcriptome assemblies from nine parasitic nematode species, including important pathogens of humans and domestic animals. Our assemblies indicated AS rates in the range of 20-30%, typically with 2-3 transcripts per AS locus, depending on the species. Transcript isoforms from the nine species were translated and searched for similarity to known proteins and functional domains. Some 21 InterPro domains, including several involved in nucleotide and chromatin binding, were statistically correlated with AS genetic loci. In most cases, the Roche/454 data explored in this study are the only sequences available from the species in question; however, the recently published genome of the human hookworm Necator americanus provided an additional opportunity to validate our results. CONCLUSIONS: Our optimized assembly parameters facilitated the first survey of AS among parasitic nematodes. The nine transcriptome assemblies, their protein translations, and basic annotations are available from Nematode.net as a resource for the research community. These should be useful for studies of specific genes and gene families of interest as well as for curating draft genome assemblies as they become available.

Sepsis from the gut: the enteric habitat of bacteria that cause late-onset neonatal bloodstream infections.

BACKGROUND: Late-onset sepsis is a major problem in neonatology, but the habitat of the pathogens before bloodstream invasion occurs is not well established. METHODS: We examined prospectively collected stools from premature infants with sepsis to find pathogens that subsequently invaded their bloodstreams, and sought the same organisms in stools of infants without sepsis. Culture-based techniques were used to isolate stool bacteria that provisionally matched the bloodstream organisms, which were then genome sequenced to confirm or refute commonality. RESULTS: Of 11 children with late-onset neonatal bloodstream infections, 7 produced at least 1 stool that contained group B Streptococcus (GBS), Serratia marcescens, or Escherichia coli before their sepsis episode with provisionally matching organisms. Of 96 overlap comparison subjects without sepsis temporally associated with these cases, 4 were colonized with provisionally matching GBS or S. marcescens. Of 175 comparisons of stools from randomly selected infants without sepsis, 1 contained a GBS (this infant had also served as an overlap comparison subject and both specimens contained provisionally matching GBS). Genome sequencing confirmed common origin of provisionally matching fecal and blood isolates. The invasive E. coli were present in all presepticemic stools since birth, but gut colonization with GBS and S. marcescens occurred closer to time of bloodstream infection. CONCLUSIONS: The neonatal gut harbors sepsis-causing pathogens, but such organisms are not inevitable members of the normal microbiota. Surveillance microbiology, decolonization, and augmented hygiene might prevent dissemination of invasive bacteria between and within premature infants.

Genome of the human hookworm Necator americanus.

The hookworm Necator americanus is the predominant soil-transmitted human parasite. Adult worms feed on blood in the small intestine, causing iron-deficiency anemia, malnutrition, growth and development stunting in children, and severe morbidity and mortality during pregnancy in women. We report sequencing and assembly of the N. americanus genome (244 Mb, 19,151 genes). Characterization of this first hookworm genome sequence identified genes orchestrating the hookworm's invasion of the human host, genes involved in blood feeding and development, and genes encoding proteins that represent new potential drug targets against hookworms. N. americanus has undergone a considerable and unique expansion of immunomodulator proteins, some of which we highlight as potential treatments against inflammatory diseases. We also used a protein microarray to demonstrate a postgenomic application of the hookworm genome sequence. This genome provides an invaluable resource to boost ongoing efforts toward fundamental and applied postgenomic research, including the development of new methods to control hookworm and human immunological diseases.

Precise dissection of an Escherichia coli O157:H7 outbreak by single nucleotide polymorphism analysis.

The current pathogen-typing methods have suboptimal sensitivities and specificities. DNA sequencing offers an opportunity to type pathogens with greater degrees of discrimination using single nucleotide polymorphisms (SNPs) than with pulsed-field gel electrophoresis (PFGE) and other methodologies. In a recent cluster of Escherichia coli O157:H7 infections attributed to salad bar exposures and romaine lettuce, a subset of cases denied exposure to either source, although PFGE and multiple-locus variable-number tandem-repeat analysis (MLVA) suggested that all isolates had the same recent progenitor. Interrogation of a preselected set of 3,442,673 nucleotides in backbone open reading frames (ORFs) identified only 1 or 2 single nucleotide differences in 3 of 12 isolates from the cases who denied exposure. The backbone DNAs of 9 of 9 and 3 of 3 cases who reported or were unsure about exposure, respectively, were isogenic. Backbone ORF SNP set sequencing offers pathogen differentiation capabilities that exceed those of PFGE and MLVA.

Discovery of anthelmintic drug targets and drugs using chokepoints in nematode metabolic pathways.

Parasitic roundworm infections plague more than 2 billion people (1/3 of humanity) and cause drastic losses in crops and livestock. New anthelmintic drugs are urgently needed as new drug resistance and environmental concerns arise. A "chokepoint reaction" is defined as a reaction that either consumes a unique substrate or produces a unique product. A chokepoint analysis provides a systematic method of identifying novel potential drug targets. Chokepoint enzymes were identified in the genomes of 10 nematode species, and the intersection and union of all chokepoint enzymes were found. By studying and experimentally testing available compounds known to target proteins orthologous to nematode chokepoint proteins in public databases, this study uncovers features of chokepoints that make them successful drug targets. Chemogenomic screening was performed on drug-like compounds from public drug databases to find existing compounds that target homologs of nematode chokepoints. The compounds were prioritized based on chemical properties frequently found in successful drugs and were experimentally tested using Caenorhabditis elegans. Several drugs that are already known anthelmintic drugs and novel candidate targets were identified. Seven of the compounds were tested in Caenorhabditis elegans and three yielded a detrimental phenotype. One of these three drug-like compounds, Perhexiline, also yielded a deleterious effect in Haemonchus contortus and Onchocerca lienalis, two nematodes with divergent forms of parasitism. Perhexiline, known to affect the fatty acid oxidation pathway in mammals, caused a reduction in oxygen consumption rates in C. elegans and genome-wide gene expression profiles provided an additional confirmation of its mode of action. Computational modeling of Perhexiline and its target provided structural insights regarding its binding mode and specificity. Our lists of prioritized drug targets and drug-like compounds have potential to expedite the discovery of new anthelmintic drugs with broad-spectrum efficacy.

Gene expression analysis distinguishes tissue-specific and gender-related functions among adult Ascaris suum tissues.

Over a billion people are infected by Ascaris spp. intestinal parasites. To clarify functional differences among tissues of adult A. suum, we compared gene expression by various tissues of these worms by expression microarray methods. The A. suum genome was sequenced and assembled to allow generation of microarray elements. Expression of over 40,000 60-mer elements was investigated in a variety of tissues from both male and female adult worms. Nearly 50 percent of the elements for which signal was detected exhibited differential expression among different tissues. The unique profile of transcripts identified for each tissue clarified functional distinctions among tissues, such as chitin binding in the ovary and peptidase activity in the intestines. Interestingly, hundreds of gender-specific elements were characterized in multiple non-reproductive tissues of female or male worms, with most prominence of gender differences in intestinal tissue. A. suum genes from the same family were frequently expressed differently among tissues. Transcript abundance for genes specific to A. suum, by comparison to Caenorhabditis elegans, varied to a greater extent among tissues than for genes conserved between A. suum and C. elegans. Analysis using C. elegans protein interaction data identified functional modules conserved between these two nematodes, resulting in identification of functional predictions of essential subnetworks of protein interactions and how these networks may vary among nematode tissues. A notable finding was very high module similarity between adult reproductive tissues and intestine. Our results provide the most comprehensive assessment of gene expression among tissues of a parasitic nematode to date.

Using existing drugs as leads for broad spectrum anthelmintics targeting protein kinases.

As one of the largest protein families, protein kinases (PKs) regulate nearly all processes within the cell and are considered important drug targets. Much research has been conducted on inhibitors for PKs, leading to a wealth of compounds that target PKs that have potential to be lead anthelmintic drugs. Identifying compounds that have already been developed to treat neglected tropical diseases is an attractive way to obtain lead compounds inexpensively that can be developed into much needed drugs, especially for use in developing countries. In this study, PKs from nematodes, hosts, and DrugBank were identified and classified into kinase families and subfamilies. Nematode proteins were placed into orthologous groups that span the phylum Nematoda. A minimal kinome for the phylum Nematoda was identified, and properties of the minimal kinome were explored. Orthologous groups from the minimal kinome were prioritized for experimental testing based on RNAi phenotype of the Caenorhabditis elegans ortholog, transcript expression over the life-cycle and anatomic expression patterns. Compounds linked to targets in DrugBank belonging to the same kinase families and subfamilies in the minimal nematode kinome were extracted. Thirty-five compounds were tested in the non-parasitic C. elegans and active compounds progressed to testing against nematode species with different modes of parasitism, the blood-feeding Haemonchus contortus and the filarial Brugia malayi. Eighteen compounds showed efficacy in C. elegans, and six compounds also showed efficacy in at least one of the parasitic species. Hypotheses regarding the pathway the compounds may target and their molecular mechanism for activity are discussed.

Widespread colonization of the lung by Tropheryma whipplei in HIV infection.

RATIONALE: Lung infections caused by opportunistic or virulent pathogens are a principal cause of morbidity and mortality in HIV infection. It is unknown whether HIV infection leads to changes in basal lung microflora, which may contribute to chronic pulmonary complications that increasingly are being recognized in individuals infected with HIV. OBJECTIVES: To determine whether the immunodeficiency associated with HIV infection resulted in alteration of the lung microbiota. METHODS: We used 16S ribosomal RNA targeted pyrosequencing and shotgun metagenomic sequencing to analyze bacterial gene sequences in bronchoalveolar lavage (BAL) and mouths of 82 HIV-positive and 77 HIV-negative subjects. MEASUREMENTS AND MAIN RESULTS: Sequences representing Tropheryma whipplei, the etiologic agent of Whipple's disease, were significantly more frequent in BAL of HIV-positive compared with HIV-negative individuals. T. whipplei dominated the community (>50% of sequence reads) in 11 HIV-positive subjects, but only 1 HIV-negative individual (13.4 versus 1.3%; P = 0.0018). In 30 HIV-positive individuals sampled longitudinally, antiretroviral therapy resulted in a significantly reduced relative abundance of T. whipplei in the lung. Shotgun metagenomic sequencing was performed on eight BAL samples dominated by T. whipplei 16S ribosomal RNA. Whole genome assembly of pooled reads showed that uncultured lung-derived T. whipplei had similar gene content to two isolates obtained from subjects with Whipple's disease. CONCLUSIONS: Asymptomatic subjects with HIV infection have unexpected colonization of the lung by T. whipplei, which is reduced by effective antiretroviral therapy and merits further study for a potential pathogenic role in chronic pulmonary complications of HIV infection.

A nuclear single-nucleotide polymorphism (SNP) potentially useful for the separation of Rhodnius prolixus from members of the Rhodnius robustus cryptic species complex (Hemiptera: Reduviidae).

The design and application of rational strategies that rely on accurate species identification are pivotal for effective vector control. When morphological identification of the target vector species is impractical, the use of molecular markers is required. Here we describe a non-coding, single-copy nuclear DNA fragment that contains a single-nucleotide polymorphism (SNP) with the potential to distinguish the important domestic Chagas disease vector, Rhodnius prolixus, from members of the four sylvatic Rhodnius robustus cryptic species complex. A total of 96 primer pairs obtained from whole genome shotgun sequencing of the R. prolixus genome (12,626 random reads) were tested on 43 R. prolixus and R. robustus s.l. samples. One of the seven amplicons selected (AmpG) presented a SNP, potentially diagnostic for R. prolixus, on the 280th site. The diagnostic nature of this SNP was then confirmed based on the analysis of 154 R. prolixus and R. robustus s.l. samples representing the widest possible geographic coverage. The results of a 60% majority-rule Bayesian consensus tree and a median-joining network constructed based on the genetic variability observed reveal the paraphyletic nature of the R. robustus species complex, with respect to R. prolixus. The AmpG region is located in the fourth intron of the Transmembrane protein 165 gene, which seems to be in the R. prolixus X chromosome. Other possible chromosomal locations of the AmpG region in the R. prolixus genome are also presented and discussed.

Transcriptomic and proteomic analyses of a Wolbachia-free filarial parasite provide evidence of trans-kingdom horizontal gene transfer.

Most filarial parasites in the subfamilies Onchocercinae and Dirofilariinae depend on Wolbachia endobacteria to successfully carry out their life cycle. Recently published data indicate that the few Wolbachia-free species in these subfamilies were infected in the distant past and have subsequently shed their endosymbionts. We used an integrated transcriptomic and proteomic analysis of Onchocerca flexuosa to explore the molecular mechanisms that allow worms of this species to survive without a bacterial partner. Roche/454 sequencing of the adult transcriptome produced 16,814 isogroup and 47,252 singleton sequences that are estimated to represent approximately 41% of the complete gene set. Sequences similar to 97 Wolbachia genes were identified from the transcriptome, some of which appear on the same transcripts as sequences similar to nematode genes. Computationally predicted peptides, including those with similarity to Wolbachia proteins, were classified at the domain and pathway levels in order to assess the metabolic capabilities of O. flexuosa and compare against the Wolbachia-dependent model filaria, Brugia malayi. Transcript data further facilitated a shotgun proteomic analysis of O. flexuosa adult worm lysate, resulting in the identification of 1,803 proteins. Three of the peptides detected by mass spectroscopy map to two ABC transport-related proteins from Wolbachia. Antibodies raised to one of the Wolbachia-like peptides labeled a single 38 kDa band on Western blots of O. flexuosa lysate and stained specific worm tissues by immunohistology. Future studies will be required to determine the exact functions of Wolbachia-like peptides and proteins in O. flexuosa and to assess their roles in worm biology.

Exploring metazoan evolution through dynamic and holistic changes in protein families and domains.

BACKGROUND: Proteins convey the majority of biochemical and cellular activities in organisms. Over the course of evolution, proteins undergo normal sequence mutations as well as large scale mutations involving domain duplication and/or domain shuffling. These events result in the generation of new proteins and protein families. Processes that affect proteome evolution drive species diversity and adaptation. Herein, change over the course of metazoan evolution, as defined by birth/death and duplication/deletion events within protein families and domains, was examined using the proteomes of 9 metazoan and two outgroup species. RESULTS: In studying members of the three major metazoan groups, the vertebrates, arthropods, and nematodes, we found that the number of protein families increased at the majority of lineages over the course of metazoan evolution where the magnitude of these increases was greatest at the lineages leading to mammals. In contrast, the number of protein domains decreased at most lineages and at all terminal lineages. This resulted in a weak correlation between protein family birth and domain birth; however, the correlation between domain birth and domain member duplication was quite strong. These data suggest that domain birth and protein family birth occur via different mechanisms, and that domain shuffling plays a role in the formation of protein families. The ratio of protein family birth to protein domain birth (domain shuffling index) suggests that shuffling had a more demonstrable effect on protein families in nematodes and arthropods than in vertebrates. Through the contrast of high and low domain shuffling indices at the lineages of Trichinella spiralis and Gallus gallus, we propose a link between protein redundancy and evolutionary changes controlled by domain shuffling; however, the speed of adaptation among the different lineages was relatively invariant. Evaluating the functions of protein families that appeared or disappeared at the last common ancestors (LCAs) of the three metazoan clades supports a correlation with organism adaptation. Furthermore, bursts of new protein families and domains in the LCAs of metazoans and vertebrates are consistent with whole genome duplications. CONCLUSION: Metazoan speciation and adaptation were explored by birth/death and duplication/deletion events among protein families and domains. Our results provide insights into protein evolution and its bearing on metazoan evolution.

A case study for large-scale human microbiome analysis using JCVI's metagenomics reports (METAREP).

As metagenomic studies continue to increase in their number, sequence volume and complexity, the scalability of biological analysis frameworks has become a rate-limiting factor to meaningful data interpretation. To address this issue, we have developed JCVI Metagenomics Reports (METAREP) as an open source tool to query, browse, and compare extremely large volumes of metagenomic annotations. Here we present improvements to this software including the implementation of a dynamic weighting of taxonomic and functional annotation, support for distributed searches, advanced clustering routines, and integration of additional annotation input formats. The utility of these improvements to data interpretation are demonstrated through the application of multiple comparative analysis strategies to shotgun metagenomic data produced by the National Institutes of Health Roadmap for Biomedical Research Human Microbiome Project (HMP) (http://nihroadmap.nih.gov). Specifically, the scalability of the dynamic weighting feature is evaluated and established by its application to the analysis of over 400 million weighted gene annotations derived from 14 billion short reads as predicted by the HMP Unified Metabolic Analysis Network (HUMAnN) pipeline. Further, the capacity of METAREP to facilitate the identification and simultaneous comparison of taxonomic and functional annotations including biological pathway and individual enzyme abundances from hundreds of community samples is demonstrated by providing scenarios that describe how these data can be mined to answer biological questions related to the human microbiome. These strategies provide users with a reference of how to conduct similar large-scale metagenomic analyses using METAREP with their own sequence data, while in this study they reveal insights into the nature and extent of variation in taxonomic and functional profiles across body habitats and individuals. Over one thousand HMP WGS datasets and the latest open source code are available at http://www.jcvi.org/hmp-metarep.