FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires.

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation.

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci.

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Capturing pairwise and multi-way chromosomal conformations using chromosomal walks.

Chromosomes are folded into highly compacted structures to accommodate physical constraints within nuclei and to regulate access to genomic information. Recently, global mapping of pairwise contacts showed that loops anchoring topological domains (TADs) are highly conserved between cell types and species. Whether pairwise loops synergize to form higher-order structures is still unclear. Here we develop a conformation capture assay to study higher-order organization using chromosomal walks (C-walks) that link multiple genomic loci together into proximity chains in human and mouse cells. This approach captures chromosomal structure at varying scales. Inter-chromosomal contacts constitute only 7-10% of the pairs and are restricted by interfacing TADs. About half of the C-walks stay within one chromosome, and almost half of those are restricted to intra-TAD spaces. C-walks that couple 2-4 TADs indicate stochastic associations between transcriptionally active, early replicating loci. Targeted analysis of thousands of 3-walks anchored at highly expressed genes support pairwise, rather than hub-like, chromosomal topology at active loci. Polycomb-repressed Hox domains are shown by the same approach to enrich for synergistic hubs. Together, the data indicate that chromosomal territories, TADs, and intra-TAD loops are primarily driven by nested, possibly dynamic, pairwise contacts.

Mapping and characterizing N6-methyladenine in eukaryotic genomes using single-molecule real-time sequencing.

N6-Methyladenine (m6dA) has been discovered as a novel form of DNA methylation prevalent in eukaryotes; however, methods for high-resolution mapping of m6dA events are still lacking. Single-molecule real-time (SMRT) sequencing has enabled the detection of m6dA events at single-nucleotide resolution in prokaryotic genomes, but its application to detecting m6dA in eukaryotic genomes has not been rigorously examined. Herein, we identified unique characteristics of eukaryotic m6dA methylomes that fundamentally differ from those of prokaryotes. Based on these differences, we describe the first approach for mapping m6dA events using SMRT sequencing specifically designed for the study of eukaryotic genomes and provide appropriate strategies for designing experiments and carrying out sequencing in future studies. We apply the novel approach to study two eukaryotic genomes. For green algae, we construct the first complete genome-wide map of m6dA at single-nucleotide and single-molecule resolution. For human lymphoblastoid cells (hLCLs), it was necessary to integrate SMRT sequencing data with independent sequencing data. The joint analyses suggest putative m6dA events are enriched in the promoters of young full-length LINE-1 elements (L1s), but call for validation by additional methods. These analyses demonstrate a general method for rigorous mapping and characterization of m6dA events in eukaryotic genomes.

Elucidation of de novo small insertion/deletion biology with parent-of-origin phasing.

The mechanisms underlying de novo insertion/deletion (indel) genesis, such as polymerase slippage, have been hypothesized but not well characterized in the human genome. We implemented two methodological improvements, which were leveraged to dissect indel mutagenesis. We assigned de novo variants to parent-of-origin (i.e., phasing) with low-coverage long-read whole-genome sequencing, achieving better phasing compared to short-read sequencing (medians of 84% and 23%, respectively). We then wrote an application programming interface to classify indels into three subtypes according to sequence context. Across three cohorts with different phasing methods (Ntrios = 540, all cohorts), we observed that one de novo indel subtype, change in copy count (CCC), was significantly correlated with father's (p = 7.1 × 10-4 ) but not mother's (p = .45) age at conception. We replicated this effect in three cohorts without de novo phasing (ppaternal = 1.9 × 10-9 , pmaternal = .61; Ntrios = 3,391, all cohorts). Although this is consistent with polymerase slippage during spermatogenesis, the percentage of variance explained by paternal age was low, and we did not observe an association with replication timing. These results suggest that spermatogenesis-specific events have a minor role in CCC indel mutagenesis, one not observed for other indel subtypes nor for maternal age in general. These results have implications for indel modeling in evolution and disease.

Epigenetic preconditioning with decitabine sensitizes glioblastoma to temozolomide via induction of MLH1.

INTRODUCTION: To improve the standard treatment paradigm for glioblastoma (GBM), efforts have been made to explore the efficacy of epigenetic agents as chemosensitizers. Recent data suggest possible synergy between decitabine (DAC), a DNA hypomethylating agent, and temozolomide (TMZ) in GBM, but the mechanism remains unclear. The objective of this study was to determine the effects of DAC on TMZ sensitization in a consecutively derived set of primary GBM cultures, with a focus on mismatch repair (MMR) proteins. METHODS: Half maximal inhibitory concentrations (IC50) of TMZ were calculated in eleven consecutive patient-derived GBM cell lines before and after preconditioning with DAC. MMR protein expression changes were determined by quantitative immunoblots and qPCR arrays. Single-molecule real-time (SMRT) sequencing of bisulfite (BS)-converted PCR amplicons of the MLH1 promoter was performed to determine methylation status. RESULTS: TMZ IC50 significantly changed in 6 of 11 GBM lines of varying MGMT promoter methylation status in response to DAC preconditioning. Knockdown of MLH1 after preconditioning reversed TMZ sensitization. SMRT-BS sequencing of the MLH1 promoter region revealed higher levels of baseline methylation at proximal CpGs in desensitized lines compared to sensitized lines. CONCLUSIONS: DAC enhances TMZ cytotoxicity in a subset of GBM cell lines, comprising lines both MGMT methylated and unmethylated tumors. This effect may be driven by levels of MLH1 via E2F1 transcription factor binding. Using unbiased long-range next-generation bisulfite-sequencing, we identified a region of the proximal MLH1 promoter with differential methylation patterns that has potential utility as a clinical biomarker for TMZ sensitization.

A comparison of immunoglobulin IGHV, IGHD and IGHJ genes in wild-derived and classical inbred mouse strains.

The genomes of classical inbred mouse strains include genes derived from all three major subspecies of the house mouse, Mus musculus. We recently posited that genetic diversity in the immunoglobulin heavy chain (IGH) gene loci of C57BL/6 and BALB/c mice reflects differences in subspecies origin. To investigate this hypothesis, we conducted high-throughput sequencing of IGH gene rearrangements to document IGH variable (IGHV), joining (IGHJ) and diversity (IGHD) genes in four inbred wild-derived mouse strains (CAST/EiJ, LEWES/EiJ, MSM/MsJ and PWD/PhJ) and a single disease model strain (NOD/ShiLtJ), collectively representing genetic backgrounds of several major mouse subspecies. A total of 341 germline IGHV sequences were inferred in the wild-derived strains, including 247 not curated in the international ImMunoGeneTics information system. By contrast, 83/84 inferred NOD IGHV genes had previously been observed in C57BL/6 mice. Variability among the strains examined was observed for only a single IGHJ gene, involving a description of a novel allele. By contrast, unexpected variation was found in the IGHD gene loci, with four previously unreported IGHD gene sequences being documented. Very few IGHV sequences of C57BL/6 and BALB/c mice were shared with strains representing major subspecies, suggesting that their IGH loci may be complex mosaics of genes of disparate origins. This suggests a similar level of diversity is likely present in the IGH loci of other classical inbred strains. This must now be documented if we are to properly understand interstrain variation in models of antibody-mediated disease.

A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak.

Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices. In this study, we established a long-read technology-based WGS screening program of all first-episode methicillin-resistant Staphylococcus aureus (MRSA) blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads enabled detailed characterization of an outbreak lasting several months of a CC5/ST105/USA100 clone among 18 infants in a neonatal intensive care unit (NICU). Available hospital-wide genome surveillance data traced the origins of the outbreak to three patients admitted to adult wards during a 4-month period preceding the NICU outbreak. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its progression, which was characterized by multiple subtransmissions and likely precipitated by equipment sharing between adults and infants. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance, and persistence. This included a DNA recognition domain recombination in the hsdS gene of a type I restriction modification system that altered DNA methylation. Transcriptome sequencing (RNA-Seq) profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall, our findings demonstrate the utility of long-read sequencing for hospital surveillance and for characterizing accessory genomic elements that may impact MRSA virulence and persistence.

HIV-1 Infection of Primary CD4+ T Cells Regulates the Expression of Specific Human Endogenous Retrovirus HERV-K (HML-2) Elements.

Endogenous retroviruses (ERVs) occupy extensive regions of the human genome. Although many of these retroviral elements have lost their ability to replicate, those whose insertion took place more recently, such as the HML-2 group of HERV-K elements, still retain intact open reading frames and the capacity to produce certain viral RNA and/or proteins. Transcription of these ERVs is, however, tightly regulated by dedicated epigenetic control mechanisms. Nonetheless, it has been reported that some pathological states, such as viral infections and certain cancers, coincide with ERV expression, suggesting that transcriptional reawakening is possible. HML-2 elements are reportedly induced during HIV-1 infection, but the conserved nature of these elements has, until recently, rendered their expression profiling problematic. Here, we provide comprehensive HERV-K HML-2 expression profiles specific for productively HIV-1-infected primary human CD4+ T cells. We combined enrichment of HIV-1 infected cells using a reporter virus expressing a surface reporter for gentle and efficient purification with long-read single-molecule real-time sequencing. We show that three HML-2 proviruses-6q25.1, 8q24.3, and 19q13.42-are upregulated on average between 3- and 5-fold in HIV-1-infected CD4+ T cells. One provirus, HML-2 12q24.33, in contrast, was repressed in the presence of active HIV replication. In conclusion, this report identifies the HERV-K HML-2 loci whose expression profiles differ upon HIV-1 infection in primary human CD4+ T cells. These data will help pave the way for further studies on the influence of endogenous retroviruses on HIV-1 replication.IMPORTANCE Endogenous retroviruses inhabit big portions of our genome. Moreover, although they are mainly inert, some of the evolutionarily younger members maintain the ability to express both RNA and proteins. We have developed an approach using long-read single-molecule real-time (SMRT) sequencing that produces long reads that allow us to obtain detailed and accurate HERV-K HML-2 expression profiles. We applied this approach to study HERV-K expression in the presence or absence of productive HIV-1 infection of primary human CD4+ T cells. In addition to using SMRT sequencing, our strategy also includes the magnetic selection of the infected cells so that levels of background expression due to uninfected cells are kept at a minimum. The results presented here provide a blueprint for in-depth studies of the interactions of the authentic upregulated HERV-K HML-2 elements and HIV-1.

Genome Plasticity of agr-Defective Staphylococcus aureus during Clinical Infection.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes.

Assembly and diploid architecture of an individual human genome via single-molecule technologies.

We present the first comprehensive analysis of a diploid human genome that combines single-molecule sequencing with single-molecule genome maps. Our hybrid assembly markedly improves upon the contiguity observed from traditional shotgun sequencing approaches, with scaffold N50 values approaching 30 Mb, and we identified complex structural variants (SVs) missed by other high-throughput approaches. Furthermore, by combining Illumina short-read data with long reads, we phased both single-nucleotide variants and SVs, generating haplotypes with over 99% consistency with previous trio-based studies. Our work shows that it is now possible to integrate single-molecule and high-throughput sequence data to generate de novo assembled genomes that approach reference quality.

Parallel Epidemics of Community-Associated Methicillin-Resistant Staphylococcus aureus USA300 Infection in North and South America.

BACKGROUND: The community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) epidemic in the United States is attributed to the spread of the USA300 clone. An epidemic of CA-MRSA closely related to USA300 has occurred in northern South America (USA300 Latin-American variant, USA300-LV). Using phylogenomic analysis, we aimed to understand the relationships between these 2 epidemics. METHODS: We sequenced the genomes of 51 MRSA clinical isolates collected between 1999 and 2012 from the United States, Colombia, Venezuela, and Ecuador. Phylogenetic analysis was used to infer the relationships and times since the divergence of the major clades. RESULTS: Phylogenetic analyses revealed 2 dominant clades that segregated by geographical region, had a putative common ancestor in 1975, and originated in 1989, in North America, and in 1985, in South America. Emergence of these parallel epidemics coincides with the independent acquisition of the arginine catabolic mobile element (ACME) in North American isolates and a novel copper and mercury resistance (COMER) mobile element in South American isolates. CONCLUSIONS: Our results reveal the existence of 2 parallel USA300 epidemics that shared a recent common ancestor. The simultaneous rapid dissemination of these 2 epidemic clades suggests the presence of shared, potentially convergent adaptations that enhance fitness and ability to spread.

Methicillin-Susceptible, Vancomycin-Resistant Staphylococcus aureus, Brazil.

We report characterization of a methicillin-susceptible, vancomycin-resistant bloodstream isolate of Staphylococcus aureus recovered from a patient in Brazil. Emergence of vancomycin resistance in methicillin-susceptible S. aureus would indicate that this resistance trait might be poised to disseminate more rapidly among S. aureus and represents a major public health threat.

Global analysis of mRNA decay intermediates in Bacillus subtilis wild-type and polynucleotide phosphorylase-deletion strains.

Messenger RNA decay in Bacillus subtilis is accomplished by a combination of exoribonucleases and endoribonucleases. Intermediates in the decay process have not been readily detectable, and previous studies on mRNA decay have used a handful of highly expressed transcripts as models. Here, we use RNA-Seq analysis to probe mRNA turnover globally. A significant fraction of messages showed differential accumulation of RNA fragments that mapped near the 5' or 3' end of the coding sequence, consistent with initiation of decay from either the 5' end or from an internal cleavage site. Patterns of mRNA decay in the wild type were compared with patterns in a mutant strain lacking polynucleotide phosphorylase (PNPase), which is considered the major 3' exonuclease activity in mRNA decay and which is one of four known 3' exonucleases in B. subtilis. The results showed a striking dependence on PNPase for mRNA turnover in many cases, suggesting specificity in the ability of 3' exonucleases to degrade from 3'-hydroxyl termini. RNA-Seq data demonstrated a sharp decrease in expression of Sigma D in the PNPase-deletion strain. Reduction in sigD regulon expression explained the chain growth phenotype of the PNPase mutant and also predicted a defect in swarming motility.

Evolution of hypervirulence by a MRSA clone through acquisition of a transposable element.

Staphylococcus aureus has evolved as a pathogen that causes a range of diseases in humans. There are two dominant modes of evolution thought to explain most of the virulence differences between strains. First, virulence genes may be acquired from other organisms. Second, mutations may cause changes in the regulation and expression of genes. Here we describe an evolutionary event in which transposition of an IS element has a direct impact on virulence gene regulation resulting in hypervirulence. Whole-genome analysis of a methicillin-resistant S. aureus (MRSA) strain USA500 revealed acquisition of a transposable element (IS256) that is absent from close relatives of this strain. Of the multiple copies of IS256 found in the USA500 genome, one was inserted in the promoter sequence of repressor of toxins (Rot), a master transcriptional regulator responsible for the expression of virulence factors in S. aureus. We show that insertion into the rot promoter by IS256 results in the derepression of cytotoxin expression and increased virulence. Taken together, this work provides new insight into evolutionary strategies by which S. aureus is able to modify its virulence properties and demonstrates a novel mechanism by which horizontal gene transfer directly impacts virulence through altering toxin regulation.

Whole-genome sequencing identifies emergence of a quinolone resistance mutation in a case of Stenotrophomonas maltophilia bacteremia.

Whole-genome sequences for Stenotrophomonas maltophilia serial isolates from a bacteremic patient before and after development of levofloxacin resistance were assembled de novo and differed by one single-nucleotide variant in smeT, a repressor for multidrug efflux operon smeDEF. Along with sequenced isolates from five contemporaneous cases, they displayed considerable diversity compared against all published complete genomes. Whole-genome sequencing and complete assembly can conclusively identify resistance mechanisms emerging in S. maltophilia strains during clinical therapy.

HIV-1 interacts with human endogenous retrovirus K (HML-2) envelopes derived from human primary lymphocytes.

UNLABELLED: Human endogenous retroviruses (HERVs) are viruses that have colonized the germ line and spread through vertical passage. Only the more recently acquired HERVs, such as the HERV-K (HML-2) group, maintain coding open reading frames. Expression of HERV-Ks has been linked to different pathological conditions, including HIV infection, but our knowledge on which specific HERV-Ks are expressed in primary lymphocytes currently is very limited. To identify the most expressed HERV-Ks in an unbiased manner, we analyzed their expression patterns in peripheral blood lymphocytes using Pacific Biosciences (PacBio) single-molecule real-time (SMRT) sequencing. We observe that three HERV-Ks (KII, K102, and K18) constitute over 90% of the total HERV-K expression in primary human lymphocytes of five different donors. We also show experimentally that two of these HERV-K env sequences (K18 and K102) retain their ability to produce full-length and posttranslationally processed envelope proteins in cell culture. We show that HERV-K18 Env can be incorporated into HIV-1 but not simian immunodeficiency virus (SIV) particles. Moreover, HERV-K18 Env incorporation into HIV-1 virions is dependent on HIV-1 matrix. Taken together, we generated high-resolution HERV-K expression profiles specific for activated human lymphocytes. We found that one of the most abundantly expressed HERV-K envelopes not only makes a full-length protein but also specifically interacts with HIV-1. Our findings raise the possibility that these endogenous retroviral Env proteins could directly influence HIV-1 replication. IMPORTANCE: Here, we report the HERV-K expression profile of primary lymphocytes from 5 different healthy donors. We used a novel deep-sequencing technology (PacBio SMRT) that produces the long reads necessary to discriminate the complexity of HERV-K expression. We find that primary lymphocytes express up to 32 different HERV-K envelopes, and that at least two of the most expressed Env proteins retain their ability to make a protein. Importantly, one of them, the envelope glycoprotein of HERV-K18, is incorporated into HIV-1 in an HIV matrix-specific fashion. The ramifications of such interactions are discussed, as the possibility of HIV-1 target tissue broadening and immune evasion are considered.

Isoform-resolved transcriptome of the human preimplantation embryo.

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

5' End-independent RNase J1 endonuclease cleavage of Bacillus subtilis model RNA.

Bacillus subtilis trp leader RNA is a small (140-nucleotide) RNA that results from attenuation of trp operon transcription upon binding of the regulatory TRAP complex. Previously, endonucleolytic cleavage by ribonuclease RNase J1 in a 3'-proximal, single-stranded region was shown to be critical for initiation of trp leader RNA decay. RNase J1 is a dual-specificity enzyme, with both 5' exonucleolytic and endonucleolytic activities. Here, we provide in vivo and in vitro evidence that RNase J1 accesses its internal target site on trp leader RNA in a 5' end-independent manner. This has important implications for the role of RNase J1 in RNA decay. We also tested the involvement in trp leader RNA decay of the more recently discovered endonuclease RNase Y. Half-lives of several trp leader RNA constructs, which were designed to probe pathways of endonucleolytic versus exonucleolytic decay, were measured in an RNase Y-deficient mutant. Remarkably, the half-lives of these constructs were indistinguishable from their half-lives in an RNase J1-deficient mutant. These results suggest that lowering RNase Y concentration may affect RNA decay indirectly via an effect on RNase J1, which is thought to exist with RNase Y in a degradosome complex. To generalize our findings with trp leader RNA to other RNAs, we show that the mechanism of trp leader RNA decay is not dependent on TRAP binding.