NeSSM: A Next-Generation Sequencing Simulator for Metagenomics

Background

Metagenomics can reveal the vast majority of microbes that have been missed by traditional cultivation-based methods. Due to its extremely wide range of application areas, fast metagenome sequencing simulation systems with high fidelity are in great demand to facilitate the development and comparison of metagenomics analysis tools.

Results

We present here a customizable metagenome simulation system: NeSSM (Next-generation Sequencing Simulator for Metagenomics). Combining complete genomes currently available, a community composition table, and sequencing parameters, it can simulate metagenome sequencing better than existing systems. Sequencing error models based on the explicit distribution of errors at each base and sequencing coverage bias are incorporated in the simulation. In order to improve the fidelity of simulation, tools are provided by NeSSM to estimate the sequencing error models, sequencing coverage bias and the community composition directly from existing metagenome sequencing data. Currently, NeSSM supports single-end and pair-end sequencing for both 454 and Illumina platforms. In addition, a GPU (graphics processing units) version of NeSSM is also developed to accelerate the simulation. By comparing the simulated sequencing data from NeSSM with experimental metagenome sequencing data, we have demonstrated that NeSSM performs better in many aspects than existing popular metagenome simulators, such as MetaSim, GemSIM and Grinder. The GPU version of NeSSM is more than one-order of magnitude faster than MetaSim.

Conclusions

NeSSM is a fast simulation system for high-throughput metagenome sequencing. It can be helpful to develop tools and evaluate strategies for metagenomics analysis and it’s freely available for academic users at http://cbb.sjtu.edu.cn/~ccwei/pub/software/NeSSM.php.     

Human ISCA1 Interacts with IOP1/NARFL and Functions in Both Cytosolic and Mitochondrial Iron-Sulfur Protein Biogenesis

Iron-sulfur proteins play an essential role in many biologic processes. Hence, understanding their assembly is an important goal. In Escherichia coli, the protein IscA is a product of the isc (iron-sulfur cluster) operon and functions in the iron-sulfur cluster assembly pathway in this organism. IscA is conserved in evolution, but its function in mammalian cells is not known. Here, we provide evidence for a role for a human homologue of IscA, named IscA1, in iron-sulfur protein biogenesis. We observe that small interfering RNA knockdown of IscA1 in HeLa cells leads to decreased activity of two mitochondrial iron-sulfur enzymes, succinate dehydrogenase and mitochondrial aconitase, as well as a cytosolic iron-sulfur enzyme, cytosolic aconitase. IscA1 is observed both in cytosolic and mitochondrial fractions. We find that IscA1 interacts with IOP1 (iron-only hydrogenase-like protein 1)/NARFL (nuclear prelamin A recognition factor-like), a cytosolic protein that plays a role in the cytosolic iron-sulfur protein assembly pathway. We therefore propose that human IscA1 plays an important role in both mitochondrial and cytosolic iron-sulfur cluster biogenesis, and a notable component of the latter is the interaction between IscA1 and IOP1.     

Toward Doubled Haploid Production in the Fabaceae: Progress,Constraints, and Opportunities

The Fabaceae species have a major role to play in sustainable farming systems, but they have lagged behind other families in respect to the development of doubled haploid protocols for plant improvement. Currently, no plant improvement program uses doubled haploids on a routine basis for any member of the Fabaceae. There has recently been renewed interest in haploid research as the usefulness of doubled haploid material in molecular mapping has become clear. This review provides a comprehensive account of the current information regarding the development of haploid protocols in the Fabaceae. In the Fabaceae crop species there have been isolated reports of haploid plant induction in the phaseoloid clade; soybean, cowpea and pigeonpea, as well as promising progress towards haploidy in peanut and winged bean. As yet there have been no reports of haploid plant production in the galegoid clade, but early stage haploid embryogenesis has been achieved in chickpea, field pea, and lupin. Success in the production of haploid plants has also been reported within the pasture genera Lotus, Medicago, and Trifolium and the arboreal genera Cassia, Peltophorum, and Albizzia. A review of the literature has enabled us to identify some general similarities between the protocols developed for haploid plant induction across the various legumes. These are the culture of intact anthers; use of a cold pretreatment to induce sporophytic development; targeting of microspores at the uninucleate stage of development; and use of MS (Murashige and Skoog, 1962 Murashige, T. and Skoog, F. 1962. A revised medium for rapid growth and bioassays with tobacco tissue cultures. Physiol. Plant., 15: 473–497. [CSA][Crossref], [Web of Science ®] [Google Scholar]) based nutrient medium with plant growth regulators to encourage continued division following induction. These protocol commonalities will assist researchers to identify approaches suited to their target Fabaceae species. The paucity of research funding for haploid research in most Fabaceae species has highlighted the need for strong collaborative linkages between institutions and researchers.

     

Ankyrin Repeat Domain Protein 2 and Inhibitor of DNA Binding 3 Cooperatively Inhibit Myoblast Differentiation by Physical Interaction

Ankyrin repeat domain protein 2 (ANKRD2) translocates from the nucleus to the cytoplasm upon myogenic induction. Overexpression of ANKRD2 inhibits C2C12 myoblast differentiation. However, the mechanism by which ANKRD2 inhibits myoblast differentiation is unknown. We demonstrate that the primary myoblasts of mdm (muscular dystrophy with myositis) mice (pMB^mdm) overexpress ANKRD2 and ID3 (inhibitor of DNA binding 3) proteins and are unable to differentiate into myotubes upon myogenic induction. Although suppression of either ANKRD2 or ID3 induces myoblast differentiation in mdm mice, overexpression of ANKRD2 and inhibition of ID3 or vice versa is insufficient to inhibit myoblast differentiation in WT mice. We identified that ANKRD2 and ID3 cooperatively inhibit myoblast differentiation by physical interaction. Interestingly, although MyoD activates the Ankrd2 promoter in the skeletal muscles of wild-type mice, SREBP-1 (sterol regulatory element binding protein-1) activates the same promoter in the skeletal muscles of mdm mice, suggesting the differential regulation of Ankrd2. Overall, we uncovered a novel pathway in which SREBP-1/ANKRD2/ID3 activation inhibits myoblast differentiation, and we propose that this pathway acts as a critical determinant of the skeletal muscle developmental program.     

CoMeta: Classification of Metagenomes Using k-mers

Nowadays, the study of environmental samples has been developing rapidly. Characterization of the environment composition broadens the knowledge about the relationship between species composition and environmental conditions. An important element of extracting the knowledge of the sample composition is to compare the extracted fragments of DNA with sequences derived from known organisms. In the presented paper, we introduce an algorithm called CoMeta (Classification of metagenomes), which assigns a query read (a DNA fragment) into one of the groups previously prepared by the user. Typically, this is one of the taxonomic rank (e.g., phylum, genus), however prepared groups may contain sequences having various functions. In CoMeta, we used the exact method for read classification using short subsequences (k-mers) and fast program for indexing large set of k-mers. In contrast to the most popular methods based on BLAST, where the query is compared with each reference sequence, we begin the classification from the top of the taxonomy tree to reduce the number of comparisons. The presented experimental study confirms that CoMeta outperforms other programs used in this context. CoMeta is available at https://github.com/jkawulok/cometa under a free GNU GPL 2 license.     

Evolution of the B3 DNA Binding Superfamily: New Insights into REM Family Gene Diversification

Background

The B3 DNA binding domain includes five families: auxin response factor (ARF), abscisic acid-insensitive3 (ABI3), high level expression of sugar inducible (HSI), related to ABI3/VP1 (RAV) and reproductive meristem (REM). The release of the complete genomes of the angiosperm eudicots Arabidopsis thaliana and Populus trichocarpa, the monocot Orysa sativa, the bryophyte Physcomitrella patens,the green algae Chlamydomonas reinhardtii and Volvox carteri and the red algae Cyanidioschyzon melorae provided an exceptional opportunity to study the evolution of this superfamily.

Methodology

In order to better understand the origin and the diversification of B3 domains in plants, we combined comparative phylogenetic analysis with exon/intron structure and duplication events. In addition, we investigated the conservation and divergence of the B3 domain during the origin and evolution of each family.

Conclusions

Our data indicate that showed that the B3 containing genes have undergone extensive duplication events, and that the REM family B3 domain has a highly diverged DNA binding. Our results also indicate that the founding member of the B3 gene family is likely to be similar to the ABI3/HSI genes found in C. reinhardtii and V. carteri. Among the B3 families, ABI3, HSI, RAV and ARF are most structurally conserved, whereas the REM family has experienced a rapid divergence. These results are discussed in light of their functional and evolutionary roles in plant development.     

HeurAA: Accurate and Fast Detection of Genetic Variations with a Novel Heuristic Amplicon Aligner Program for Next Generation Sequencing

Next generation sequencing (NGS) of PCR amplicons is a standard approach to detect genetic variations in personalized medicine such as cancer diagnostics. Computer programs used in the NGS community often miss insertions and deletions (indels) that constitute a large part of known human mutations. We have developed HeurAA, an open source, heuristic amplicon aligner program. We tested the program on simulated datasets as well as experimental data from multiplex sequencing of 40 amplicons in 12 oncogenes collected on a 454 Genome Sequencer from lung cancer cell lines. We found that HeurAA can accurately detect all indels, and is more than an order of magnitude faster than previous programs. HeurAA can compare reads and reference sequences up to several thousand base pairs in length, and it can evaluate data from complex mixtures containing reads of different gene-segments from different samples. HeurAA is written in C and Perl for Linux operating systems, the code and the documentation are available for research applications at http://sourceforge.net/projects/heuraa/     

Cysteine Is Not the Sulfur Source for Iron-Sulfur Cluster and Methionine Biosynthesis in the Methanogenic Archaeon Methanococcus maripaludis

Three multiprotein systems are known for iron-sulfur (Fe-S) cluster biogenesis in prokaryotes and eukaryotes as follows: the NIF (nitrogen fixation), the ISC (iron-sulfur cluster), and the SUF (mobilization of sulfur) systems. In all three, cysteine is the physiological sulfur source, and the sulfur is transferred from cysteine desulfurase through a persulfidic intermediate to a scaffold protein. However, the biochemical nature of the sulfur source for Fe-S cluster assembly in archaea is unknown, and many archaea lack homologs of cysteine desulfurases. Methanococcus maripaludis is a methanogenic archaeon that contains a high amount of protein-bound Fe-S clusters (45 nmol/mg protein). Cysteine in this archaeon is synthesized primarily via the tRNA-dependent SepRS/SepCysS pathway. When a ΔsepS mutant (a cysteine auxotroph) was grown with ³⁴S-labeled sulfide and unlabeled cysteine, <8% of the cysteine, >92% of the methionine, and >87% of the sulfur in the Fe-S clusters in proteins were labeled, suggesting that the sulfur in methionine and Fe-S clusters was derived predominantly from exogenous sulfide instead of cysteine. Therefore, this investigation challenges the concept that cysteine is always the sulfur source for Fe-S cluster biosynthesis in vivo and suggests that Fe-S clusters are derived from sulfide in those organisms, which live in sulfide-rich habitats.     

CEP192 interacts physically and functionally with the K63-deubiquitinase CYLD to promote mitotic spindle assembly

CEP192 is a centrosome protein that plays a critical role in centrosome biogenesis and function in mammals, Drosophila and C. elegans.¹ Dix CI, Raff JW. Drosophila Spd-2 recruits PCM to the sperm centriole, but is dispensable for centriole duplication. Curr Biol 2007; 17:1759 - 64; http://dx.doi.org/10.1016/j.cub.2007.08.065; PMID: 17919907 [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]^-6 Gomez-Ferreria MA, Rath U, Buster DW, Chanda SK, Caldwell JS, Rines DR, et al. Human Cep192 is required for mitotic centrosome and spindle assembly. Curr Biol 2007; 17:1960 - 6; http://dx.doi.org/10.1016/j.cub.2007.10.019; PMID: 17980596 [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] Moreover, CEP192-depleted cells arrest in mitosis with disorganized microtubules, suggesting that CEP192’s function in spindle assembly goes beyond its role in centrosome activity and pointing to a potentially more direct role in the regulation of the mitotic microtubule landscape.⁷ Gomez-Ferreria MA, Sharp DJ. Cep192 and the generation of the mitotic spindle. Cell Cycle 2008; 7:1507 - 10; http://dx.doi.org/10.4161/cc.7.11.5957; PMID: 18469523 [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] To better understand CEP192 function in mitosis, we used mass spectrometry to identify CEP192-interacting proteins. We previously reported that CEP192 interacts with NEDD1, a protein that associates with the γ-tubulin ring complex (γ-TuRC) and regulates its phosphorylation status during mitosis.⁸ Gomez-Ferreria M, Bashkurov M, Helbig A, Larsen B, Pawson T, Gingras AC, et al. Novel NEDD1 phosphorylation sites regulate γ-tubulin bindingand mitotic spindle assembly. J Cell Sci 2012; http://dx.doi.org/10.1242/jcs.105130; PMID: 22595525 [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] Additionally, within the array of proteins that interact with CEP192, we identified the microtubule binding K63-deubiquitinase CYLD. Further analyses show that co-depletion of CYLD alleviates the bipolar spindle assembly defects observed in CEP192-depleted cells. This functional relationship exposes an intriguing role for CYLD in spindle formation and raises the tantalizing possibility that CEP192 promotes robust mitotic spindle assembly by regulating K63-polyubiquitin-mediated signaling through CYLD.     

Systems analysis of immune responses to attenuated P. falciparum malaria sporozoite vaccination reveals excessive inflammatory signatures correlating with impaired immunity

Immunization with radiation-attenuated sporozoites (RAS) can confer sterilizing protection against malaria, although the mechanisms behind this protection are incompletely understood. We performed a systems biology analysis of samples from the Immunization by Mosquito with Radiation Attenuated Sporozoites (IMRAS) trial, which comprised P. falciparum RAS-immunized (PfRAS), malaria-naive participants whose protection from malaria infection was subsequently assessed by controlled human malaria infection (CHMI). Blood samples collected after initial PfRAS immunization were analyzed to compare immune responses between protected and non-protected volunteers leveraging integrative analysis of whole blood RNA-seq, high parameter flow cytometry, and single cell CITEseq of PBMCs. This analysis revealed differences in early innate immune responses indicating divergent paths associated with protection. In particular, elevated levels of inflammatory responses early after the initial immunization were detrimental for the development of protective adaptive immunity. Specifically, non-classical monocytes and early type I interferon responses induced within 1 day of PfRAS vaccination correlated with impaired immunity. Non-protected individuals also showed an increase in Th2 polarized T cell responses whereas we observed a trend towards increased Th1 and T-bet+ CD8 T cell responses in protected individuals. Temporal differences in genes associated with natural killer cells suggest an important role in immune regulation by these cells. These findings give insight into the immune responses that confer protection against malaria and may guide further malaria vaccine development.Trial registration: ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT01994525","term_id":"NCT01994525"}}NCT01994525.     

Chemically Modified DNA Aptamers Bind Interleukin-6 with High Affinity and Inhibit Signaling by Blocking Its Interaction with Interleukin-6 Receptor

Interleukin-6 (IL-6) is a pleiotropic cytokine that regulates immune and inflammatory responses, and its overproduction is a hallmark of inflammatory diseases. Inhibition of IL-6 signaling with the anti-IL-6 receptor antibody tocilizumab has provided some clinical benefit to patients; however, direct cytokine inhibition may be a more effective option. We used the systematic evolution of ligands by exponential enrichment (SELEX) process to discover slow off-rate modified aptamers (SOMAmers) with hydrophobic base modifications that inhibit IL-6 signaling in vitro. Two classes of IL-6 SOMAmers were isolated from modified DNA libraries containing 40 random positions and either 5-(N-benzylcarboxamide)-2′-deoxyuridine (Bn-dU) or 5-[N-(1-naphthylmethyl)carboxamide]-2′-deoxyuridine (Nap-dU) replacing dT. These modifications facilitate the high affinity binding interaction with IL-6 and provide resistance against degradation by serum endonucleases. Post-SELEX optimization of one Bn-dU and one Nap-dU SOMAmer led to improvements in IL-6 binding (10-fold) and inhibition activity (greater than 20-fold), resulting in lead SOMAmers with sub-nanomolar affinity (K_d = 0.2 nm) and potency (IC₅₀ = 0.2 nm). Although similar in inhibition properties, the two SOMAmers have unique sequences and different ortholog specificities. Furthermore, these SOMAmers were stable in human serum in vitro for more than 48 h. Both SOMAmers prevented IL-6 signaling by blocking the interaction of IL-6 with its receptor and inhibited the proliferation of tumor cells in vitro as effectively as tocilizumab. This new class of IL-6 inhibitor may be an effective therapeutic alternative for patients suffering from inflammatory diseases.     

Mutants in the Candida glabrata Glycerol Channels Are Sensitized to Cell Wall Stress

Many fungal species use glycerol as a compatible solute with which to maintain osmotic homeostasis in response to changes in external osmolarity. In Saccharomyces cerevisiae, intracellular glycerol concentrations are regulated largely by the high osmolarity glycerol (HOG) response pathway, both through induction of glycerol biosynthesis and control of its flux through the plasma membrane Fps1 glycerol channel. The channel activity of Fps1 is also controlled by a pair of positive regulators, Rgc1 and Rgc2. In this study, we demonstrate that Candida glabrata, a fungal pathogen that possesses two Fps1 orthologs and two Rgc1/-2 orthologs, accumulates glycerol in response to hyperosmotic stress. We present an initial characterization of mutants with deletions in the C. glabrata FPS1 (CAGL0C03267 [www.candidagenome.org]) and FPS2 (CAGL0E03894) genes and find that a double mutant accumulates glycerol, experiences constitutive cell wall stress, and is hypersensitive to treatment by caspofungin, an antifungal agent that targets the cell wall. This mutant is cleared more efficiently in mouse infections than is wild-type C. glabrata by caspofungin treatment. Finally, we demonstrate that one of the C. glabrata RGC orthologs complements an S. cerevisiae rgc1 rgc2 null mutant, supporting the conclusion that this regulatory assembly is conserved between these species.     

MetaMetaDB: A Database and Analytic System for Investigating Microbial Habitability

MetaMetaDB (http://mmdb.aori.u-tokyo.ac.jp/) is a database and analytic system for investigating microbial habitability, i.e., how a prokaryotic group can inhabit different environments. The interaction between prokaryotes and the environment is a key issue in microbiology because distinct prokaryotic communities maintain distinct ecosystems. Because 16S ribosomal RNA (rRNA) sequences play pivotal roles in identifying prokaryotic species, a system that comprehensively links diverse environments to 16S rRNA sequences of the inhabitant prokaryotes is necessary for the systematic understanding of the microbial habitability. However, existing databases are biased to culturable prokaryotes and exhibit limitations in the comprehensiveness of the data because most prokaryotes are unculturable. Recently, metagenomic and 16S rRNA amplicon sequencing approaches have generated abundant 16S rRNA sequence data that encompass unculturable prokaryotes across diverse environments; however, these data are usually buried in large databases and are difficult to access. In this study, we developed MetaMetaDB (Meta-Metagenomic DataBase), which comprehensively and compactly covers 16S rRNA sequences retrieved from public datasets. Using MetaMetaDB, users can quickly generate hypotheses regarding the types of environments a prokaryotic group may be adapted to. We anticipate that MetaMetaDB will improve our understanding of the diversity and evolution of prokaryotes.