Search

Recent studies suggest that knowledge of differential expression of microRNAs (miRNA) in cancer may have substantial diagnostic and prognostic value. Here, we use a direct sequencing method to characterize the profiles of miRNAs and other small RNA segments for six human cervical carcinoma cell lines and five normal cervical samples. Of 166 miRNAs expressed in normal cervix and cancer cell lines, we observed significant expression variation of six miRNAs between the two groups. To further show the biological relevance of our findings, we examined the expression level of two significantly varying miRNAs in a panel of 29 matched pairs of human cervical cancer and normal cervical samples. Reduced expression of miR-143 and increased expression of miR-21 were reproducibly displayed in cancer samples, suggesting the potential value of these miRNAs as tumor markers. In addition to the known miRNAs, we found a number of novel miRNAs and an additional set of small RNAs that do not meet miRNA criteria. [Cancer Res 2007;67(13):6031–43]

Supplementary Table 6 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Figure 1 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 2 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Figure Legends 1-2 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 3 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 1 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 7 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Figure 2 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 4 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

Supplementary Table 5 from Patterns of Known and Novel Small RNAs in Human Cervical Cancer

The basis for acquired resistance to JAK inhibition in patients with JAK2-driven hematologic malignancies is not well understood. We report a patient with a myeloproliferative neoplasm (MPN) with a BCR activator of RhoGEF and GTPase (BCR)–JAK2 fusion with initial hematologic response to ruxolitinib who rapidly developed B-lymphoid blast transformation. We analyzed pre-ruxolitinib and blast transformation samples using genome sequencing, DNA mate-pair sequencing (MPseq), RNA sequencing (RNA-seq), and chromosomal microarray to characterize possible mechanisms of resistance. No resistance mutations in the BCR-JAK2 fusion gene or transcript were identified, and fusion transcript expression levels remained stable. However, at the time of blast transformation, MPseq detected a new IKZF1 copy-number loss, which is predicted to result in loss of normal IKZF1 protein translation. RNA-seq revealed significant upregulation of genes negatively regulated by IKZF1, including IL7R and CRLF2. Disease progression was also characterized by adaptation to an activated B-cell receptor (BCR)–like signaling phenotype, with marked upregulation of genes such as CD79A, CD79B, IGLL1, VPREB1, BLNK, ZAP70, RAG1, and RAG2. In summary, IKZF1 deletion and a switch from cytokine dependence to activated BCR-like signaling phenotype represent putative mechanisms of ruxolitinib resistance in this case, recapitulating preclinical data on resistance to JAK inhibition in CRLF2-rearranged Philadelphia chromosome-like acute lymphoblastic leukemia.

Background Single molecule measurements of DNA polymerization kinetics provide a sensitive means to detect both secondary structures in DNA and deviations from primary chemical structure as a result of modified bases. In one approach to such analysis, deviations can be inferred by monitoring the behavior of DNA polymerase using single-molecule, real-time sequencing with zero-mode waveguide. This approach uses a Single Molecule Real Time (SMRT)-sequencing measurement of time between fluorescence pulse signals from consecutive nucleosides incorporated during DNA replication, called the interpulse duration (IPD). Results In this paper we present an analysis of loci with high IPDs in two genomes, a bacterial genome (E. coli) and a eukaryotic genome (C. elegans). To distinguish the potential effects of DNA modification on DNA polymerization speed, we paired an analysis of native genomic DNA with whole-genome amplified (WGA) material in which DNA modifications were effectively removed. Adenine modification sites for E. coli are known and we observed the expected IPD shifts at these sites in the native but not WGA samples. For C. elegans, such differences were not observed. Instead, we found a number of novel sequence contexts where IPDs were raised relative to the average IPDs for each of the four nucleotides, but for which the raised IPD was present in both native and WGA samples. Conclusion The latter results argue strongly against DNA modification as the underlying driver for high IPD segments for C. elegans, and provide a framework for separating effects of DNA modification from context-dependent DNA polymerase kinetic patterns inherent in underlying DNA sequence for a complex eukaryotic genome.

Planarians have long been studied for their regenerative abilities. Moving forward, tools for ectopic expression of non-native proteins will be of substantial value. Using a luminescent reporter to overcome the strong autofluorescence background of planarian tissues, we demonstrate heterologous protein expression in planarian cells and live animals. Our approach is based on the introduction of mRNA through several nanotechnological and chemical transfection methods. We improve reporter expression by altering untranslated region (UTR) sequences and codon bias, facilitating measurement of expression kinetics both in isolated cells and in whole planarians using luminescence imaging. We also examine protein expression as a function of variations in the UTRs of delivered mRNA, demonstrating a framework to investigate gene regulation at the post-transcriptional level. Together, these advances expand the toolbox for the mechanistic analysis of planarian biology and establish a strong foundation for the development and expansion of transgenic techniques in this unique model system.MotivationThe study of planarians has contributed to advances in our understanding of regeneration, stem cell dynamics, and many other fundamental biological processes. However, the persistent challenge of expressing transgenes in planarians has led to the speculation that they may be resistant to transfection. In this work, we develop methods to express exogenous mRNAs in both isolated planarian cells and whole animals by optimizing delivery techniques, genetic constructs, and detection methods. These methods allow us to study transfection kinetics and post-transcriptional regulation of gene expression in a quantitative manner. Beyond planarian research, this work should also provide a broadly applicable strategy to develop similar tools for animals that are also challenging to modify genetically.

Piwi-interacting RNAs (piRNAs) are RNA effectors with key roles in maintaining genome integrity and promoting fertility in metazoans. In Caenorhabditis elegans loss of piRNAs leads to a transgenerational sterility phenotype. The plethora of piRNAs and their ability to silence transcripts with imperfect complementarity have raised several (non-exclusive) models for the underlying drivers of sterility. Here, we report the extranuclear and transferable nature of the sterility driver, its suppression via mutations disrupting the endogenous RNAi and poly-uridylation machinery, and copy-number amplification at the ribosomal DNA locus. In piRNA-deficient animals, several small interfering RNA (siRNA) populations become increasingly overabundant in the generations preceding loss of germline function, including ribosomal siRNAs (risiRNAs). A concomitant increase in uridylated sense rRNA fragments suggests that poly-uridylation may potentiate RNAi-mediated gene silencing of rRNAs. We conclude that loss of the piRNA machinery allows for unchecked amplification of siRNA populations, originating from abundant highly structured RNAs, to deleterious levels.

RNA viruses in aquatic environments remain poorly studied. Here, we analysed the RNA virome from approximately 10 l water from Yangshan Deep-Water Harbour near the Yangtze River estuary in China and identified more than 4,500 distinct RNA viruses, doubling the previously known set of viruses. Phylogenomic analysis identified several major lineages, roughly, at the taxonomic ranks of class, order and family. The 719-member-strong Yangshan virus assemblage is the sister clade to the expansive class Alsuviricetes and consists of viruses with simple genomes that typically encode only RNA-dependent RNA polymerase (RdRP), capping enzyme and capsid protein. Several clades within the Yangshan assemblage independently evolved domain permutation in the RdRP. Another previously unknown clade shares ancestry with Potyviridae, the largest known plant virus family. The ‘Aquatic picorna-like viruses/Marnaviridae’ clade was greatly expanded, with more than 800 added viruses. Several RdRP-linked protein domains not previously detected in any RNA viruses were identified, such as the small ubiquitin-like modifier (SUMO) domain, phospholipase A2 and PrsW-family protease domain. Multiple viruses utilize alternative genetic codes implying protist (especially ciliate) hosts. The results reveal a vast RNA virome that includes many previously unknown groups. However, phylogenetic analysis of the RdRPs supports the previously established five-branch structure of the RNA virus evolutionary tree, with no additional phyla., Metagenomic analysis of a single RNA virome from the Yangshan Deep-Water Harbour in China enabled the recovery of more than 4,500 distinct RNA viruses, doubling the known set of RNA viruses to date, and provided insights into their biology.

The germline genome is guarded against invading foreign genetic elements by small RNA-dependent gene-silencing pathways. Components of these pathways localize to, or form distinct aggregates in the vicinity of, germ granules. These components and their dynamics in and out of granules are currently being intensively studied. Here, we report the identification of PLP-1, a Caenorhabditis elegans protein related to the human single-stranded nucleic acid-binding protein Pur-alpha, as a component of germ granules in C. elegans. We show that PLP-1 is essential for silencing different types of transgenes in the germ line and for suppressing the expression of several endogenous genes controlled by the germline gene-silencing pathways. Our results reveal that PLP-1 functions downstream of small RNA biogenesis during initiation of gene silencing. Based on these results and the earlier findings that Pur-alpha proteins interact with both RNA and protein, we propose that PLP-1 couples certain RNAs with their protein partners in the silencing complex. PLP-1 orthologs localized on RNA granules may similarly contribute to germline gene silencing in other organisms.

Meta-metagenomic searches allow for high-speed, low-cost identification of potentially significant biological niches for sequences of interest., In numerous instances, tracking the biological significance of a nucleic acid sequence can be augmented through the identification of environmental niches in which the sequence of interest is present. Many metagenomic data sets are now available, with deep sequencing of samples from diverse biological niches. While any individual metagenomic data set can be readily queried using web-based tools, meta-searches through all such data sets are less accessible. In this brief communication, we demonstrate such a meta-metagenomic approach, examining close matches to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in all high-throughput sequencing data sets in the NCBI Sequence Read Archive accessible with the “virome” keyword. In addition to the homology to bat coronaviruses observed in descriptions of the SARS-CoV-2 sequence (F. Wu, S. Zhao, B. Yu, Y. M. Chen, et al., Nature 579:265–269, 2020, https://doi.org/10.1038/s41586-020-2008-3; P. Zhou, X. L. Yang, X. G. Wang, B. Hu, et al., Nature 579:270–273, 2020, https://doi.org/10.1038/s41586-020-2012-7), we note a strong homology to numerous sequence reads in metavirome data sets generated from the lungs of deceased pangolins reported by Liu et al. (P. Liu, W. Chen, and J. P. Chen, Viruses 11:979, 2019, https://doi.org/10.3390/v11110979). While analysis of these reads indicates the presence of a similar viral sequence in pangolin lung, the similarity is not sufficient to either confirm or rule out a role for pangolins as an intermediate host in the recent emergence of SARS-CoV-2. In addition to the implications for SARS-CoV-2 emergence, this study illustrates the utility and limitations of meta-metagenomic search tools in effective and rapid characterization of potentially significant nucleic acid sequences. IMPORTANCE Meta-metagenomic searches allow for high-speed, low-cost identification of potentially significant biological niches for sequences of interest.

Revisiting replicating RNAs DNA-dependent RNA polymerases are well known for their ability to produce RNA from DNA templates. Much less is known about their noncanonical activity: the generation and replication of RNA from RNA templates. Deeper insights into this process could shed light on questions relating to the origin of life, molecular evolution, and the replication of certain RNA pathogens such as hepatitis delta virus and plant viroids. Jain et al. explore in detail how the bacteriophage T7 RNA polymerase generates and amplifies diverse RNA sequences in vitro. Using sequencing, microfluidics, and bioinformatics, they chart the emergence and evolution of replicating RNA motifs and identify mechanisms that explain their selection and structures. Science , this issue p. eaay0688

A goal of HIV vaccine development is to elicit antibodies with neutralizing breadth. Broadly neutralizing antibodies (bNAbs) to HIV often have unusual sequences with long heavy-chain complementarity-determining region loops, high somatic mutation rates and polyreactivity. A subset of HIV-infected individuals develops such antibodies, but it is unclear whether this reflects systematic differences in their antibody repertoires or is a consequence of rare stochastic events involving individual clones. We sequenced antibody heavy-chain repertoires in a large cohort of HIV-infected individuals with bNAb responses or no neutralization breadth and uninfected controls, identifying consistent features of bNAb repertoires, encompassing thousands of B cell clones per individual, with correlated T cell phenotypes. These repertoire features were not observed during chronic cytomegalovirus infection in an independent cohort. Our data indicate that the development of numerous B cell lineages with antibody features associated with autoreactivity may be a key aspect in the development of HIV neutralizing antibody breadth., A subset of HIV-infected individuals can develop broadly neutralizing antibodies. Boyd and colleagues studied such HIV-infected individuals and found significant perturbations in their antibody repertoires, including increased frequencies of B cells expressing antibodies with features associated with autoreactivity.

In numerous instances, tracking the biological significance of a nucleic acid sequence can be augmented through the identification of environmental niches in which the sequence of interest is present. Many metagenomic datasets are now available, with deep sequencing of samples from diverse biological niches. While any individual metagenomic dataset can be readily queried using web-based tools, meta-searches through all such datasets are less accessible. In this brief communication, we demonstrate such a meta-meta-genomic approach, examining close matches to the Wuhan coronavirus 2019-nCoV in all high-throughput sequencing datasets in the NCBI Sequence Read Archive accessible with the keyword “virome”. In addition to the homology to bat coronaviruses observed in descriptions of the 2019-nCoV sequence (F. Wu et al. 2020, Nature, doi.org/10.1038/s41586-020-2008-3; P. Zhou et al. 2020, Nature, doi.org/10.1038/s41586-020-2012-7), we note a strong homology to numerous sequence reads in a metavirome dataset generated from the lungs of deceased Pangolins reported by Liu et al. (Viruses 11:11, 2019, http://doi.org/10.3390/v11110979). Our observations are relevant to discussions of the derivation of 2019-nCoV and illustrate the utility and limitations of meta-metagenomic search tools in effective and rapid characterization of potentially significant nucleic acid sequences.ImportanceMeta-metagenomic searches allow for high-speed, low-cost identification of potentially significant biological niches for sequences of interest.

The germline genome is guarded against invading foreign genetic elements by small RNA-dependent gene-silencing pathways. Components of these pathways localize to, or form distinct aggregates in the vicinity of, germ granules. These components and their dynamics in and out of granules are currently being intensively studied. Here, we report the identification of PLP-1, a Caenorhabditiselegans protein related to the human single-stranded nucleic acid-binding protein Pur-alpha, as a component of germ granules in C. elegans. We show that PLP-1 is essential for silencing different types of transgenes in the germ line and for suppressing the expression of several endogenous genes controlled by the germline gene-silencing pathways. Our results reveal that PLP-1 functions downstream of small RNA biogenesis during initiation of gene silencing. Based on these results and the earlier findings that Pur-alpha proteins interact with both RNA and protein, we propose that PLP-1 couples certain RNAs with their protein partners in the silencing complex. PLP-1 orthologs localized on RNA granules may similarly contribute to germline gene silencing in other organisms.

Investigations aimed at defining the 3D configuration of eukaryotic chromosomes have consistently encountered an endogenous population of chromosome-derived circular genomic DNA, referred to as extrachromosomal circular DNA (eccDNA). While the production, distribution, and activities of eccDNAs remain understudied, eccDNA formation from specific regions of the linear genome has profound consequences on the regulatory and coding capabilities for these regions. Here, we define eccDNA distributions in Caenorhabditis elegans and in three human cell types, utilizing a set of DNA topology-dependent approaches for enrichment and characterization. The use of parallel biophysical, enzymatic, and informatic approaches provides a comprehensive profiling of eccDNA robust to isolation and analysis methodology. Results in human and nematode systems provide quantitative analysis of the eccDNA loci at both unique and repetitive regions. Our studies converge on and support a consistent picture, in which endogenous genomic DNA circles are present in normal physiological states, and in which the circles come from both coding and noncoding genomic regions. Prominent among the coding regions generating DNA circles are several genes known to produce a diversity of protein isoforms, with mucin proteins and titin as specific examples.

Interactions between Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) RNAs and CRISPR-associated (Cas) proteins form an RNA-guided adaptive immune system in prokaryotes. The adaptive immune system utilizes segments of the genetic material of invasive foreign elements in the CRISPR locus. The loci are transcribed and processed to produce small CRISPR RNAs (crRNAs), with degradation of invading genetic material directed by a combination of complementarity between RNA and DNA and in some cases recognition of adjacent motifs called PAMs (Protospacer Adjacent Motifs). Here we describe a general, high-throughput procedure to test the efficacy of thousands of targets, applying this to the Escherichia coli type I-E Cascade (CRISPR-associated complex for antiviral defense) system. These studies were followed with reciprocal experiments in which the consequence of CRISPR activity was survival in the presence of a lytic phage. From the combined analysis of the Cascade system, we found that (i) type I-E Cascade PAM recognition is more expansive than previously reported, with at least 22 distinct PAMs, with many of the noncanonical PAMs having CRISPR-interference abilities similar to the canonical PAMs; (ii) PAM positioning appears precise, with no evidence for tolerance to PAM slippage in interference; and (iii) while increased guanine-cytosine (GC) content in the spacer is associated with higher CRISPR-interference efficiency, high GC content (>62.5%) decreases CRISPR-interference efficiency. Our findings provide a comprehensive functional profile of Cascade type I-E interference requirements and a method to assay spacer efficacy that can be applied to other CRISPR-Cas systems.

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted ≥53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology.

Next-generation DNA-sequencing (NGS) technologies, which are designed to streamline the acquisition of massive amounts of sequencing data, are nonetheless dependent on various preparative steps to generate DNA fragments of required concentration, purity, and average size (molecular weight). Current automated electrophoresis systems for DNA- and RNA-sample quality control, such as Agilent’s Bioanalyzer ®and TapeStation ® products, are costly to acquire and use; they also provide limited information for samples having broad size distributions. Here we describe a software tool that helps determine the size distribution of DNA fragments in an NGS library, or other DNA sample, based on gel-electrophoretic line profiles. The software, developed as an ImageJ plug-in, allows for straightforward processing of gel images, including lane selection and fitting of univariate functions to intensity distributions. The user selects the option of fitting either discrete profiles in cases where discrete gel bands are visible, or continuous profiles, having multiple bands buried under a single broad peak. The method requires only modest imaging capabilities and is a cost-effective, rigorous alternative characterization method to augment existing techniques for library quality control.

Assessments of biological complexity for populations that are of mixed species are central in many biological contexts, including microbiomes, tumor cell population structure, and immune cell populations. Here we address the problem of quantifying the population diversity in experiments where high throughput DNA sequencing is used to distinguish a large number of cell subpopulations. Our model assumes a list of clonal species and their observed frequencies in each of several replicate sequencing libraries. Though the underlying distribution of frequencies cannot be estimated well from data coming from only a small fraction of the total cell population, one can estimate well the population-level clonality, defined as the sum of squared underlying fractions of the respective clones, the complement of the Gini–Simpson index. Specifically, we proposed to adaptively combine multiple unbiased estimators of clonality derived from pairs of replicates to construct a single estimator without relying on the commonly used but restrictive multinomial assumption. The new estimator performs particularly well for replicates of unequal size. We further illustrate the proposed methods with extensive simulations and a small real data example.

Caenorhabditis elegans provides an amenable system to explore whether newly composed ribosomes are required to progress through development. Despite the complex pattern of tissues that are formed during embryonic development, we found that null homozygotes lacking any of the five different ribosomal proteins (RPs) can produce fully functional first-stage larvae, with similar developmental competence seen upon complete deletion of the multi-copy ribosomal RNA locus. These animals, relying on maternal but not zygotic contribution of ribosomal components, are capable of completing embryogenesis. In the absence of new ribosomal components, the resulting animals are arrested before progression from the first larval stage and fail in two assays for postembryonic plasticity of neuronal structure. Mosaic analyses of larvae that are a mixture of ribosome-competent and non-competent cells suggest a global regulatory mechanism in which ribosomal insufficiency in a subset of cells triggers organism-wide growth arrest.

Cas9, a CRISPR-associated RNA-guided nuclease, has been rapidly adopted as a tool for biochemical and genetic manipulation of DNA. Although complexes between Cas9 and guide RNAs (gRNAs) offer remarkable specificity and versatility for genome manipulation, mis-targeted events occur. To extend the understanding of gRNA::target homology requirements, we compared mutational tolerance for a set of Cas9::gRNA complexes in vitro and in vivo (in Saccharomyces cerevisiae). A variety of gRNAs were tested with variant libraries based on four different targets (with varying GC content and sequence features). In each case, we challenged a mixture of matched and mismatched targets, evaluating cleavage activity on a wide variety of potential target sequences in parallel through high-throughput sequencing of the products retained after cleavage. These experiments evidenced notable and consistent differences between in vitro and S. cerevisiae (in vivo) Cas9 cleavage specificity profiles including (i) a greater tolerance for mismatches in vitro and (ii) a greater specificity increase in vivo with truncation of the gRNA homology regions.

Background Specific immunotherapy (SIT) is the only treatment with proved long-term curative potential in patients with allergic disease. Allergen-specific IgE is the causative agent of allergic disease, and antibodies contribute to SIT, but the effects of SIT on aeroallergen-specific B-cell repertoires are not well understood. Objective We sought to characterize the IgE sequences expressed by allergen-specific B cells and track the fate of these B-cell clones during SIT. Methods We used high-throughput antibody gene sequencing and identification of allergen-specific IgE with combinatorial antibody fragment library technology to analyze immunoglobulin repertoires of blood and the nasal mucosa from aeroallergen-sensitized subjects before and during the first year of subcutaneous SIT. Results Of 52 distinct allergen-specific IgE heavy chains from 8 allergic donors, 37 were also detected by using high-throughput antibody gene sequencing of blood samples, nasal mucosal samples, or both. The allergen-specific clones had increased persistence, higher likelihood of belonging to clones expressing other switched isotypes, and possibly larger clone size than the rest of the IgE repertoire. Clone members in nasal tissue showed close mutational relationships. Conclusion In the future, combining functional binding studies, deep antibody repertoire sequencing, and information on clinical outcomes in larger studies might aid assessment of SIT mechanisms and efficacy.

Chikungunya is caused by the mosquito-borne arthrogenic alphavirus, chikungunya virus (CHIKV). Chikungunya was introduced into the Americas in late 2013 and Nicaragua in mid-2014. Here, we sequenced five imported and 30 autochthonous Nicaraguan CHIKV from cases identified in the first epidemic in the country between August 2014 and April 2015. One full-length and two partial genomic sequences were obtained by deep sequencing; Sanger methodology yielded 33 E1 sequences from five imported and 28 autochthonous cases. Phylogenetic analysis indicates that Nicaraguan CHIKV all belonged to the Asian genotype, Caribbean clade. Moreover, E1 gene sequences revealed accumulation of mutations in later months of the epidemic, including four silent mutations in 11 autochthonous cases and three non-synonymous mutations in three autochthonous cases. No mutations contributing to increased transmissibility by Aedes albopictus were identified in the E1 gene. This represents the most comprehensive set of CHIKV sequences available from the Americas to date.

In the course of identifying and cleaving RNA, the RNAi machinery must encounter and contend with the megadalton-sized ribosomes that carry out translation. We investigated this interface by examining the fate of actively translated mRNAs subjected to RNAi in C. elegans. Quantifying RNA levels (RNA-seq) and ongoing translation (Ribo-seq), we found there is a greater fold repression of ongoing translation than expected from loss of RNA alone, observing stronger translation repression relative to RNA repression for multiple, independent double-stranded RNA triggers, and for multiple genes. In animals that lack the RNA helicase SKI complex and the ribosome rescue factor PELOTA, ribosomes stall on the 3′ edges of mRNAs at and upstream of the RNAi trigger. One model to explain these observations is that ribosomes are actively cleared from mRNAs by SKI and PELO during or following mRNA cleavage. Our results expand prior studies that show a role for the SKI RNA helicase complex in removing RNA targets following RNAi in flies and plants, illuminating the widespread role of the nonstop translation surveillance in RNA silencing during RNAi. Our results are also consistent with proposals that RNAi can attack messages during active translation.

Background and objectives A kidney disease of unknown cause is common in Sri Lanka’s lowland (dry) region. Detailed clinical characterizations of patients with biopsy-proven disease are limited, and there is no current consensus on criteria for a noninvasive diagnosis. Design, setting, participants, & measurements We designed a prospective study in a major Sri Lankan hospital servicing endemic areas to ascertain pathologic and clinical characteristics of and assess risk factors for primary tubulointerstitial kidney disease. We used logistic regression to determine whether common clinical characteristics could be used to predict the presence of primary tubulointerstitial kidney disease on kidney biopsy. Results From 600 new patients presenting to a tertiary nephrology clinic over the course of 1 year, 87 underwent kidney biopsy, and 43 (49%) had a biopsy diagnosis of primary tubulointerstitial kidney disease. On detailed biopsy review, 13 (30%) had evidence of moderate to severe active kidney disease, and six (15%) had evidence of moderate to severe chronic tubulointerstitial kidney disease. Patients with tubulointerstitial kidney disease were exclusively born in endemic provinces; 91% spent a majority of their lifespan there. They were more likely men and farmers (risk ratio, 2.0; 95% confidence interval, 1.2 to 2.9), and they were more likely to have used tobacco (risk ratio, 1.7; 95% confidence interval, 1.0 to 2.3) and well water (risk ratio, 1.5; 95% confidence interval, 1.1 to 2.0). Three clinical characteristics—age, urine dipstick for protein, and serum albumin—could predict likelihood of tubulointerstitial kidney disease on biopsy (model sensitivity of 79% and specificity of 84%). Patients referred for kidney biopsy despite comorbid diabetes or hypertension did not experience lower odds of tubulointerstitial kidney disease. Conclusions A primary tubulointerstitial kidney disease occurs commonly in specific regions of Sri Lanka with characteristic environmental and lifestyle exposures.

Next generation high-throughput sequencing has enabled sensitive and unambiguous analysis of RNA populations in cells. Here, we describe a method for isolation and strand-specific sequencing of small RNA pools from bacteria that can be multiplexed to accommodate multiple biological samples in a single experiment. Small RNAs are isolated by polyacrylamide gel electrophoresis and treated with T4 polynucleotide kinase. This allows for 3' adapter ligation to CRISPR RNAs, which don't have pre-existing 3'-OH ends. Pre-adenylated adapters are then ligated using T4 RNA ligase 1 in the absence of ATP and with a high concentration of polyethylene glycol (PEG). The 3' capture step enables precise determination of the 3' ends of diverse RNA molecules. Additionally, a random hexamer in the ligated adapter helps control for potential downstream amplification bias. Following reverse-transcription, the cDNA product is circularized and libraries are prepared by PCR. We show that the amplified library need not be visible by gel electrophoresis for efficient sequencing of the desired product. Using this method, we routinely prepare RNA sequencing libraries from minute amounts of purified small RNA. This protocol is tailored to assay for CRISPR RNA biogenesis in bacteria through sequencing of mature CRISPR RNAs, but can be used to sequence diverse classes of small RNAs. We also provide a fully worked example of our data processing pipeline, with instructions for running the provided scripts.

CRISPR-Cas-mediated defense utilizes information stored as spacers in CRISPR arrays to defend against genetic invaders. We define the mode of target interference and role in antiviral defense for two CRISPR-Cas systems in Marinomonas mediterranea. One system (type I-F) targets DNA. A second system (type III-B) is broadly capable of acquiring spacers in either orientation from RNA and DNA, and exhibits transcription-dependent DNA interference. Examining resistance to phages isolated from Mediterranean seagrass meadows, we found that the type III-B machinery co-opts type I-F CRISPR-RNAs. Sequencing and infectivity assessments of related bacterial and phage strains suggests an ‘arms race’ in which phage escape from the type I-F system can be overcome through use of type I-F spacers by a horizontally-acquired type III-B system. We propose that the phage-host arms race can drive selection for horizontal uptake and maintenance of promiscuous type III interference modules that supplement existing host type I CRISPR-Cas systems.

Cas1 integrase is the key enzyme of the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas adaptation module that mediates acquisition of spacers derived from foreign DNA by CRISPR arrays. In diverse bacteria, the cas1 gene is fused (or adjacent) to a gene encoding a reverse transcriptase (RT) related to group II intron RTs. An RT-Cas1 fusion protein has been recently shown to enable acquisition of CRISPR spacers from RNA. Phylogenetic analysis of the CRISPR-associated RTs demonstrates monophyly of the RT-Cas1 fusion, and coevolution of the RT and Cas1 domains. Nearly all such RTs are present within type III CRISPR-Cas loci, but their phylogeny does not parallel the CRISPR-Cas type classification, indicating that RT-Cas1 is an autonomous functional module that is disseminated by horizontal gene transfer and can function with diverse type III systems. To compare the sequence pools sampled by RT-Cas1-associated and RT-lacking CRISPR-Cas systems, we obtained samples of a commercially grown cyanobacterium—Arthrospira platensis. Sequencing of the CRISPR arrays uncovered a highly diverse population of spacers. Spacer diversity was particularly striking for the RT-Cas1-containing type III-B system, where no saturation was evident even with millions of sequences analyzed. In contrast, analysis of the RT-lacking type III-D system yielded a highly diverse pool but reached a point where fewer novel spacers were recovered as sequencing depth was increased. Matches could be identified for a small fraction of the non-RT-Cas1-associated spacers, and for only a single RT-Cas1-associated spacer. Thus, the principal source(s) of the spacers, particularly the hypervariable spacer repertoire of the RT-associated arrays, remains unknown., IMPORTANCE While the majority of CRISPR-Cas immune systems adapt to foreign genetic elements by capturing segments of invasive DNA, some systems carry reverse transcriptases (RTs) that enable adaptation to RNA molecules. From analysis of available bacterial sequence data, we find evidence that RT-based RNA adaptation machinery has been able to join with CRISPR-Cas immune systems in many, diverse bacterial species. To investigate whether the abilities to adapt to DNA and RNA molecules are utilized for defense against distinct classes of invaders in nature, we sequenced CRISPR arrays from samples of commercial-scale open-air cultures of Arthrospira platensis, a cyanobacterium that contains both RT-lacking and RT-containing CRISPR-Cas systems. We uncovered a diverse pool of naturally occurring immune memories, with the RT-lacking locus acquiring a number of segments matching known viral or bacterial genes, while the RT-containing locus has acquired spacers from a distinct sequence pool for which the source remains enigmatic.

Investigations aimed at defining the 3-D configuration of eukaryotic chromosomes have consistently encountered an endogenous population of chromosome-derived circular genomic DNA, referred to as extrachromosomal circular DNA (eccDNA). While the production, distribution, and activities of eccDNAs remain understudied, eccDNA formation from specific regions of the linear genome has profound consequences on the regulatory and coding capabilities for these regions. High-throughput sequencing has only recently made extensive genomic mapping of eccDNA sequences possible and had yet to be applied using a rigorous approach that distinguishes ascertainment bias from true enrichment. Here, we define eccDNA distribution, utilizing a set of unbiased topology-dependent approaches for enrichment and characterization. We use parallel biophysical, enzymatic, and informatic approaches to obtain a comprehensive profiling of eccDNA inC. elegansand in three human cell types, where eccDNAs were previously uncharacterized. We also provide quantitative analysis of the eccDNA loci at both unique and repetitive regions. Our studies converge on and support a consistent picture in which endogenous genomic DNA circles are present in normal physiological DNA metabolism, and in which the circles come from both coding and noncoding genomic regions. Prominent among the coding regions generating DNA circles are several genes known to produce a diversity of protein isoforms, with mucin proteins and titin as specific examples.

Conventional plasmid vectors are incapable of achieving sustained levels of transgene expression in vivo even in quiescent mammalian tissues because the transgene expression cassette is silenced. Transcriptional silencing results from the presence of the bacterial plasmid backbone or virtually any DNA sequence of >1 kb in length placed outside of the expression cassette. Here, we show that transcriptional silencing can be substantially forestalled by increasing the An/Tn sequence composition in the plasmid bacterial backbone. Increasing numbers of An/Tn sequences increased sustained transcription of both backbone sequences and adjacent expression cassettes. In order to recapitulate these expression profiles in compact and portable plasmid DNA backbones, we engineered the standard kanamycin or ampicillin antibiotic resistance genes, optimizing the number of An/Tn sequence without altering the encoded amino acids. The resulting vector backbones yield sustained transgene expression from mouse liver, providing generic DNA vectors capable of sustained transgene expression without additional genes or mammalian regulatory elements.

FLT3 fusions are associated with myeloid and lymphoid neoplasms with eosinophilia. We describe a patient presenting with clinicopathologic features of both chronic eosinophilic leukemia, not otherwise specified (CEL, NOS) and systemic mastocytosis (SM). The bone marrow demonstrated a myeloproliferative neoplasm with eosinophilia and aggregates of atypical mast cells. Cytogenetic analysis revealed a t(13;14)(q12;q32), which was subsequently molecularly characterized as a novel TRIP11-FLT3 rearrangement. A KIT D816V mutation was also identified. The patient rapidly transformed to T-lymphoblastic leukemia/lymphoma and expired shortly after diagnosis. This is the fifth FLT3 fusion gene described in the literature; the presence of both myeloid and lymphoid neoplasms implicates involvement of an early hematopoietic progenitor by rearranged FLT3. We suggest that leukemias and lymphomas with FLT3 fusion genes exhibit similar clinicopathologic features to, and should be included in, the WHO category of "Myeloid and lymphoid neoplasms with eosinophilia and abnormalities of PDGFRA, PDGFRB, or FGFR1, or with PCM1-JAK2."