Genomics of Clostridium taeniosporum, an organism which forms endospores with ribbon-like appendages.

Clostridium taeniosporum, a non-pathogenic anaerobe closely related to the C. botulinum Group II members, was isolated from Crimean lake silt about 60 years ago. Its endospores are surrounded by an encasement layer which forms a trunk at one spore pole to which about 12-14 large, ribbon-like appendages are attached. The genome consists of one 3,264,813 bp, circular chromosome (with 26.6% GC) and three plasmids. The chromosome contains 2,892 potential protein coding sequences: 2,124 have specific functions, 147 have general functions, 228 are conserved but without known function and 393 are hypothetical based on the fact that no statistically significant orthologs were found. The chromosome also contains 101 genes for stable RNAs, including 7 rRNA clusters. Over 84% of the protein coding sequences and 96% of the stable RNA coding regions are oriented in the same direction as replication. The three known appendage genes are located within a single cluster with five other genes, the protein products of which are closely related, in terms of sequence, to the known appendage proteins. The relatedness of the deduced protein products suggests that all or some of the closely related genes might code for minor appendage proteins or assembly factors. The appendage genes might be unique among the known clostridia; no statistically significant orthologs were found within other clostridial genomes for which sequence data are available. The C. taeniosporum chromosome contains two functional prophages, one Siphoviridae and one Myoviridae, and one defective prophage. Three plasmids of 5.9, 69.7 and 163.1 Kbp are present. These data are expected to contribute to future studies of developmental, structural and evolutionary biology and to potential industrial applications of this organism.

RNA-seq of human reference RNA samples using a thermostable group II intron reverse transcriptase.

Next-generation RNA sequencing (RNA-seq) has revolutionized our ability to analyze transcriptomes. Current RNA-seq methods are highly reproducible, but each has biases resulting from different modes of RNA sample preparation, reverse transcription, and adapter addition, leading to variability between methods. Moreover, the transcriptome cannot be profiled comprehensively because highly structured RNAs, such as tRNAs and snoRNAs, are refractory to conventional RNA-seq methods. Recently, we developed a new method for strand-specific RNA-seq using thermostable group II intron reverse transcriptases (TGIRTs). TGIRT enzymes have higher processivity and fidelity than conventional retroviral reverse transcriptases plus a novel template-switching activity that enables RNA-seq adapter addition during cDNA synthesis without using RNA ligase. Here, we obtained TGIRT-seq data sets for well-characterized human RNA reference samples and compared them to previous data sets obtained for these RNAs by the Illumina TruSeq v2 and v3 methods. We find that TGIRT-seq recapitulates the relative abundance of human transcripts and RNA spike-ins in ribo-depleted, fragmented RNA samples comparably to non-strand-specific TruSeq v2 and better than strand-specific TruSeq v3. Moreover, TGIRT-seq is more strand specific than TruSeq v3 and eliminates sampling biases from random hexamer priming, which are inherent to TruSeq. The TGIRT-seq data sets also show more uniform 5' to 3' gene coverage and identify more splice junctions, particularly near the 5' ends of mRNAs, than do the TruSeq data sets. Finally, TGIRT-seq enables the simultaneous profiling of mRNAs and lncRNAs in the same RNA-seq experiment as structured small ncRNAs, including tRNAs, which are essentially absent with TruSeq.

High-throughput sequencing of human plasma RNA by using thermostable group II intron reverse transcriptases.

Next-generation RNA-sequencing (RNA-seq) has revolutionized transcriptome profiling, gene expression analysis, and RNA-based diagnostics. Here, we developed a new RNA-seq method that exploits thermostable group II intron reverse transcriptases (TGIRTs) and used it to profile human plasma RNAs. TGIRTs have higher thermostability, processivity, and fidelity than conventional reverse transcriptases, plus a novel template-switching activity that can efficiently attach RNA-seq adapters to target RNA sequences without RNA ligation. The new TGIRT-seq method enabled construction of RNA-seq libraries from <1 ng of plasma RNA in <5 h. TGIRT-seq of RNA in 1-mL plasma samples from a healthy individual revealed RNA fragments mapping to a diverse population of protein-coding gene and long ncRNAs, which are enriched in intron and antisense sequences, as well as nearly all known classes of small ncRNAs, some of which have never before been seen in plasma. Surprisingly, many of the small ncRNA species were present as full-length transcripts, suggesting that they are protected from plasma RNases in ribonucleoprotein (RNP) complexes and/or exosomes. This TGIRT-seq method is readily adaptable for profiling of whole-cell, exosomal, and miRNAs, and for related procedures, such as HITS-CLIP and ribosome profiling.

Long-term neural and physiological phenotyping of a single human.

Psychiatric disorders are characterized by major fluctuations in psychological function over the course of weeks and months, but the dynamic characteristics of brain function over this timescale in healthy individuals are unknown. Here, as a proof of concept to address this question, we present the MyConnectome project. An intensive phenome-wide assessment of a single human was performed over a period of 18 months, including functional and structural brain connectivity using magnetic resonance imaging, psychological function and physical health, gene expression and metabolomics. A reproducible analysis workflow is provided, along with open access to the data and an online browser for results. We demonstrate dynamic changes in brain connectivity over the timescales of days to months, and relations between brain connectivity, gene expression and metabolites. This resource can serve as a testbed to study the joint dynamics of human brain and metabolic function over time, an approach that is critical for the development of precision medicine strategies for brain disorders.

Next-generation sequencing fragment library construction.

All current next-generation sequencing (NGS) platforms and applications require the sequencing library to have specific characteristics: in particular, size, size distribution, and 5' and 3' flanking sequences. This unit presents a robust protocol for converting a wide variety of input DNA samples into appropriate NGS libraries and discusses important considerations in experimental design, failure modes, and typical results.

Next-Generation Sequencing RNA-Seq Library Construction.

This unit presents protocols for construction of next-generation sequencing (NGS) directional RNA sequencing libraries for the Illumina HiSeq and MiSeq from a wide variety of input RNA sources. The protocols are based on the New England Biolabs (NEB) small RNA library preparation set for Illumina, although similar kits exist from different vendors. The protocol preserves the orientation of the original RNA in the final sequencing library, enabling strand-specific analysis of the resulting data. These libraries have been used for differential gene expression analysis and small RNA discovery and are currently being tested for de novo transcriptome assembly. The protocol is robust and applicable to a broad range of RNA input types and RNA quality, making it ideal for high-throughput laboratories.

A capture-sequencing strategy identifies IRF8, EBF1, and APRIL as novel IGH fusion partners in B-cell lymphoma.

The characterization of immunoglobulin heavy chain (IGH) translocations provides information on the diagnosis and guides therapeutic decisions in mature B-cell malignancies while enhancing our understanding of normal and malignant B-cell biology. However, existing methodologies for the detection of IGH translocations are labor intensive, often require viable cells, and are biased toward known IGH fusions. To overcome these limitations, we developed a capture sequencing strategy for the identification of IGH rearrangements at nucleotide level resolution and tested its capabilities as a diagnostic and discovery tool in 78 primary diffuse large B-cell lymphomas (DLBCLs). We readily identified IGH-BCL2, IGH-BCL6, IGH-MYC, and IGH-CCND1 fusions and discovered IRF8, EBF1, and TNFSF13 (APRIL) as novel IGH partners in these tumors. IRF8 and TNFSF13 expression was significantly higher in lymphomas with IGH rearrangements targeting these loci. Modeling the deregulation of IRF8 and EBF1 in vitro defined a lymphomagenic profile characterized by up-regulation of AID and/or BCL6, down-regulation of PRMD1, and resistance to apoptosis. Using a capture sequencing strategy, we discovered the B-cell relevant genes IRF8, EBF1, and TNFSF13 as novel targets for IGH deregulation. This methodology is poised to change how IGH translocations are identified in clinical settings while remaining a powerful tool to uncover the pathogenesis of B-cell malignancies.

Engineered DNA ligases with improved activities in vitro.

The DNA ligase from bacteriophage T4 is one of the most widely used enzymes in molecular biology. It has evolved to seal single-stranded nicks in double-stranded DNA, but not to join double-stranded fragments with cohesive or blunt ends. Its poor activity in vitro, particularly with blunt-ended substrates, can lead to failed or sub-optimal experimental outcomes. We have fused T4 DNA ligase to seven different DNA-binding proteins, including eukaryotic transcription factors, bacterial DNA repair proteins and archaeal DNA-binding domains. Representatives from each of these classes improved the activity of T4 DNA ligase, by up to 7-fold, in agarose gel-based screens for cohesive- and blunt-ended fragment joining. Overall, the most active variants were p50-ligase (i.e. NF-κB p50 fused to T4 DNA ligase) and ligase-cTF (T4 DNA ligase fused to an artificial, chimeric transcription factor). Ligase-cTF out-performed T4 DNA ligase by ∼160% in blunt end 'vector + insert' cloning assays, and p50-ligase showed an improvement of a similar magnitude when it was used to construct a library for Illumina sequencing. The activity of the Escherichia coli DNA ligase was also enhanced by fusion to p50. Together, these results suggest that our protein design strategy is a generalizable one for engineering improved DNA ligases.

Thermostable group II intron reverse transcriptase fusion proteins and their use in cDNA synthesis and next-generation RNA sequencing.

Mobile group II introns encode reverse transcriptases (RTs) that function in intron mobility ("retrohoming") by a process that requires reverse transcription of a highly structured, 2-2.5-kb intron RNA with high processivity and fidelity. Although the latter properties are potentially useful for applications in cDNA synthesis and next-generation RNA sequencing (RNA-seq), group II intron RTs have been difficult to purify free of the intron RNA, and their utility as research tools has not been investigated systematically. Here, we developed general methods for the high-level expression and purification of group II intron-encoded RTs as fusion proteins with a rigidly linked, noncleavable solubility tag, and we applied them to group II intron RTs from bacterial thermophiles. We thus obtained thermostable group II intron RT fusion proteins that have higher processivity, fidelity, and thermostability than retroviral RTs, synthesize cDNAs at temperatures up to 81°C, and have significant advantages for qRT-PCR, capillary electrophoresis for RNA-structure mapping, and next-generation RNA sequencing. Further, we find that group II intron RTs differ from the retroviral enzymes in template switching with minimal base-pairing to the 3' ends of new RNA templates, making it possible to efficiently and seamlessly link adaptors containing PCR-primer binding sites to cDNA ends without an RNA ligase step. This novel template-switching activity enables facile and less biased cloning of nonpolyadenylated RNAs, such as miRNAs or protein-bound RNA fragments. Our findings demonstrate novel biochemical activities and inherent advantages of group II intron RTs for research, biotechnological, and diagnostic methods, with potentially wide applications.

High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire.

Each B-cell receptor consists of a pair of heavy and light chains. High-throughput sequencing can identify large numbers of heavy- and light-chain variable regions (V(H) and V(L)) in a given B-cell repertoire, but information about endogenous pairing of heavy and light chains is lost after bulk lysis of B-cell populations. Here we describe a way to retain this pairing information. In our approach, single B cells (>5 × 10(4) capacity per experiment) are deposited in a high-density microwell plate (125 pl/well) and lysed in situ. mRNA is then captured on magnetic beads, reverse transcribed and amplified by emulsion V(H):V(L) linkage PCR. The linked transcripts are analyzed by Illumina high-throughput sequencing. We validated the fidelity of V(H):V(L) pairs identified by this approach and used the method to sequence the repertoire of three human cell subsets-peripheral blood IgG(+) B cells, peripheral plasmablasts isolated after tetanus toxoid immunization and memory B cells isolated after seasonal influenza vaccination.

The caBIG® Life Science Business Architecture Model.

MOTIVATION: Business Architecture Models (BAMs) describe what a business does, who performs the activities, where and when activities are performed, how activities are accomplished and which data are present. The purpose of a BAM is to provide a common resource for understanding business functions and requirements and to guide software development. The cancer Biomedical Informatics Grid (caBIG®) Life Science BAM (LS BAM) provides a shared understanding of the vocabulary, goals and processes that are common in the business of LS research. RESULTS: LS BAM 1.1 includes 90 goals and 61 people and groups within Use Case and Activity Unified Modeling Language (UML) Diagrams. Here we report on the model's current release, LS BAM 1.1, its utility and usage, and plans for future use and continuing development for future releases. AVAILABILITY AND IMPLEMENTATION: The LS BAM is freely available as UML, PDF and HTML (https://wiki.nci.nih.gov/x/OFNyAQ).

Monoclonal antibodies isolated without screening by analyzing the variable-gene repertoire of plasma cells.

Isolation of antigen-specific monoclonal antibodies (mAbs) and antibody fragments relies on high-throughput screening of immortalized B cells or recombinant antibody libraries. We bypassed the screening step by using high-throughput DNA sequencing and bioinformatic analysis to mine antibody variable region (V)-gene repertoires from bone marrow plasma cells (BMPC) of immunized mice. BMPCs, which cannot be immortalized, produce the vast majority of circulating antibodies. We found that the V-gene repertoire of BMPCs becomes highly polarized after immunization, with the most abundant sequences represented at frequencies between approximately 1% and >10% of the total repertoire. We paired the most abundant variable heavy (V(H)) and variable light (V(L)) genes based on their relative frequencies, reconstructed them using automated gene synthesis, and expressed recombinant antibodies in bacteria or mammalian cells. Antibodies generated in this manner from six mice, each immunized with one of three antigens were overwhelmingly antigen specific (21/27 or 78%). Those generated from a mouse with high serum titers had nanomolar binding affinities.

Small RNA profiling reveals antisense transcription throughout the KSHV genome and novel small RNAs.

Kaposi's sarcoma-associated herpesvirus (KSHV) is a human tumor virus that encodes 12 precursor microRNAs (pre-miRNAs) that give rise to 17 different known approximately 22-nucleotide (nt) effector miRNAs. Like all herpesviruses, KSHV has two modes of infection: (1) a latent mode whereby only a subset of viral genes are expressed and (2) a lytic mode during which the full remaining viral genes are expressed. To date, KSHV miRNAs have been mostly identified via analysis of cells that are undergoing latent infection. Here, we developed a method to profile small RNAs ( approximately 18-75 nt) from populations of cells undergoing predominantly lytic infection. Using two different next-generation sequencing platforms, we cloned and sequenced both pre-miRNAs and derivative miRNAs. Our analysis shows that the vast majority of viral and host 5p miRNAs are co-terminal with the 5' end of the cloned pre-miRNAs, consistent with both being defined by microprocessor cleavage. We report the complete repertoire (25 total) of 5p and 3p derivative miRNAs from all 12 previously described KSHV pre-miRNAs. Two KSHV pre-miRNAs, pre-miR-K12-8 and pre-miR-K12-12, encode abundant derivative miRNAs from the previously unreported strands of the pre-miRNA. We identify several novel small RNAs of low abundance, including viral miRNA-offset-RNAs (moRNAs), and antisense viral miRNAs (miRNA-AS) that are encoded antisense to previously reported KSHV pre-miRNAs. Finally, we observe widespread antisense transcription relative to known coding sequences during lytic replication. Despite the enormous potential to form double-stranded RNA in KSHV-infected cells, we observe no evidence for the existence of abundant viral-derived small interfering RNAs (siRNAs).

Rapid construction and characterization of synthetic antibody libraries without DNA amplification.

We report on a simple method to rapidly generate very large libraries of genes encoding mutant proteins without the use of DNA amplification, and the application of this methodology in the construction of synthetic immunoglobulin variable heavy (V(H)) and light (V(kappa)) libraries. Four high quality, chemically synthesized polynucleotides (90-140 bases) were annealed and extended using T4 DNA polymerase. Following electroporation, >10(9) transformants could be synthesized within 1 day. Fusion to beta-lactamase and selection on ampicillin resulted in 3.7 x 10(8) V(H) and 6.9 x 10(8) V(kappa) clones highly enriched for full-length, in-frame genes. High-throughput 454 DNA sequencing of >250,000 V(H) and V(kappa) genes from the pre- and post-selection libraries revealed that, in addition to the expected reduction in reading-frame shifts and stop codons, selection for functional expression also resulted in a statistical decrease in the cysteine content. Apart from these differences, there was a good agreement between the expected and actual diversity, indicating that neither oligonucleotide synthesis nor biological constrains due to protein synthesis of V(H)/V(kappa)-beta-lactamase fusions introduce biases in the amino acid composition of the randomized regions. This methodology can be employed for the rapid construction of highly diverse libraries with the near elimination of PCR errors in invariant regions.

High-throughput RNAi screening in vitro: from cell lines to primary cells.

Small interfering RNAs (siRNAs) are being used to induce sequence-specific gene silencing in cultured cells to study mammalian gene function. Libraries of siRNAs targeting entire human gene classes can be used to identify genes with specific cellular functions. Here we describe high-throughput siRNA delivery methods to facilitate siRNA library screening experiments with both immortalized and primary cells. We adapted chemical reverse transfection for immortalized adherent cell lines in a 96-well format. The method is fast, robust, and exceptionally effective for many cell types. For primary cells and immortalized cells that are recalcitrant to lipofection-based methods, we developed electropermeabilization (electroporation) conditions that facilitate siRNA delivery to a broad range of cell types, including primary human T-cells, hMSC, NHA, NDHF-Neo, HUVEC, DI TNC1, RPTEC, PC12, and K562 cells. To enable high-throughput electropermeabilization of primary cells, we developed a novel 96-well electroporation device that provides highly efficient and reproducible delivery of siRNAs. The combination of high-throughput chemical reverse transfection and electroporation makes it possible to deliver libraries of siRNAs to virtually any cell type, enabling gene function analysis and discovery on a genome scale.