Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

RNA:RNA interaction in ternary complexes resolved by chemical probing

RNA regulation can be performed by a second targeting RNA molecule, such as in the microRNA regulation mechanism. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probes the structure of RNA molecules and can resolve RNA:protein interactions, but RNA:RNA interactions have not yet been addressed with this technique. Here, we apply SHAPE to investigate RNA-mediated binding processes in RNA:RNA and RNA:RNA-RBP complexes. We use RNA:RNA binding by SHAPE (RABS) to investigate microRNA-34a (miR-34a) binding its mRNA target, the silent information regulator 1 (mSIRT1), both with and without the Argonaute protein, constituting the RNA-induced silencing complex (RISC). We show that the seed of the mRNA target must be bound to the microRNA loaded into RISC to enable further binding of the compensatory region by RISC, while the naked miR-34a is able to bind the compensatory region without seed interaction. The method presented here provides complementary structural evidence for the commonly performed luciferase-assay-based evaluation of microRNA binding-site efficiency and specificity on the mRNA target site and could therefore be used in conjunction with it. The method can be applied to any nucleic acid-mediated RNA- or RBP-binding process, such as splicing, antisense RNA binding, or regulation by RISC, providing important insight into the targeted RNA structure.

Knockouts of TUT7 and 3'hExo show that they cooperate in histone mRNA maintenance and degradation

Metazoan histone mRNAs are the only cellular eukaryotic mRNAs that are not polyadenylated, ending instead in a conserved stem-loop. SLBP is bound to the 3' end of histone mRNAs and is required for translation of histone mRNA. The expression of histone mRNAs is tightly cell-cycle regulated. A major regulatory step is rapid degradation of histone mRNA at the end of S-phase or when DNA synthesis is inhibited in S-phase. 3'hExo, a 3' to 5' exonuclease, binds to the SLBP/SL complex and trims histone mRNA to 3 nt after the stem-loop. Together with a terminal uridyl transferase, 3'hExo maintains the length of the histone mRNA during S-phase. 3'hExo is essential for initiating histone mRNA degradation on polyribosomes, initiating degradation into the 3' side of the stem-loop. There is extensive uridylation of degradation intermediates in the 3' side of the stem when histone mRNA is degraded. Here, we knocked out TUT7 and 3'hExo and we show that both modification of histone mRNA during S-phase and degradation of histone mRNA involve the interaction of 3'hExo, and a specific TUTase, TENT3B (TUT7, ZCCHC6). Knockout of 3'hExo prevents the initiation of 3' to 5' degradation, stabilizing histone mRNA, whereas knockout of TUT7 prevents uridylation of the mRNA degradation intermediates, slowing the rate of degradation. In synchronized 3'hExo KO cells, histone mRNA degradation is delayed, but the histone mRNA is degraded prior to mitosis by a different pathway.

A-MYB/TCFL5 regulatory architecture ensures the production of pachytene piRNAs in placental mammals

In male mice, the transcription factor A MYB initiates the transcription of pachytene piRNA genes during meiosis. Here, we report that A MYB activates the transcription factor Tcfl5 produced in pachytene spermatocytes. Subsequently, A MYB and TCFL5 reciprocally reinforce their own transcription to establish a positive feedback circuit that triggers pachytene piRNA production. TCFL5 regulates the expression of genes required for piRNA maturation and promotes transcription of evolutionarily young pachytene piRNA genes, whereas A-MYB activates the transcription of older pachytene piRNA genes. Intriguingly, pachytene piRNAs from TCFL5-dependent young loci initiates the production of piRNAs from A-MYB-dependent older loci ensuring the self-propagation of pachytene piRNAs. A MYB and TCFL5 act via a set of incoherent feedforward loops that drive regulation of gene expression by pachytene piRNAs during spermatogenesis. This regulatory architecture is conserved in rhesus macaque, suggesting that it was present in the last common ancestor of placental mammals.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Widespread genetic heterogeneity of human ribosomal RNA genes

Polymorphism drives survival under stress and provides adaptability. Genetic polymorphism of ribosomal RNA (rRNA) genes derives from internal repeat variation of this multicopy gene, and from interindividual variation. A considerable amount of rRNA sequence heterogeneity has been proposed but has been challenging to estimate given the scarcity of accurate reference sequences. We identified four rDNA copies on chromosome 21 (GRCh38) with 99% similarity to recently introduced reference sequence KY962518.1. We customized a GATK bioinformatics pipeline using the four rDNA loci, spanning a total 145 kb, for variant calling and used high-coverage whole-genome sequencing (WGS) data from the 1000 Genomes Project to analyze variants in 2504 individuals from 26 populations. We identified a total of 3791 variant positions. The variants positioned nonrandomly on the rRNA gene. Invariant regions included the promoter, early 5' ETS, most of 18S, 5.8S, ITS1, and large areas of the intragenic spacer. A total of 470 variant positions were observed on 28S rRNA. The majority of the 28S rRNA variants were located on highly flexible human-expanded rRNA helical folds ES7L and ES27L, suggesting that these represent positions of diversity and are potentially under continuous evolution. Several variants were validated based on RNA-seq analyses. Population analyses showed remarkable ancestry-linked genetic variance and the presence of both high penetrance and frequent variants in the 5' ETS, ITS2, and 28S regions segregating according to the continental populations. These findings provide a genetic view of rRNA gene array heterogeneity and raise the need to functionally assess how the 28S rRNA variants affect ribosome functions.

An association test of the spatial distribution of rare missense variants within protein structures identifies Alzheimer's disease-related patterns

More than 90% of genetic variants are rare in most modern sequencing studies, such as the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data. Furthermore, 54% of the rare variants in ADSP WES are singletons. However, both single variant and unit-based tests are limited in their statistical power to detect an association between rare variants and phenotypes. To best use missense rare variants and investigate their biological effect, we examine their association with phenotypes in the context of protein structures. We developed a protein structure-based approach, protein optimized kernel evaluation of missense nucleotides (POKEMON), which evaluates rare missense variants based on their spatial distribution within a protein rather than their allele frequency. The hypothesis behind this test is that the three-dimensional spatial distribution of variants within a protein structure provides functional context to power an association test. POKEMON identified three candidate genes (TREM2, SORL1, and EXOC3L4) and another suggestive gene from the ADSP WES data. For TREM2 and SORL1, two known Alzheimer's disease (AD) genes, the signal from the spatial cluster is stable even if we exclude known AD risk variants, indicating the presence of additional low-frequency risk variants within these genes. EXOC3L4 is a novel AD risk gene that has a cluster of variants primarily shared by case subjects around the Sec6 domain. This cluster is also validated in an independent replication data set and a validation data set with a larger sample size.

Genomic association with pathogen carriage in bighorn sheep (Ovis canadensis)

Genetic composition can influence host susceptibility to, and transmission of, pathogens, with potential population-level consequences. In bighorn sheep (Ovis canadensis), pneumonia epidemics caused by Mycoplasma ovipneumoniae have been associated with severe population declines and limited recovery across North America. Adult survivors either clear the infection or act as carriers that continually shed M. ovipneumoniae and expose their susceptible offspring, resulting in high rates of lamb mortality for years following the outbreak event. Here, we investigated the influence of genomic composition on persistent carriage of M. ovipneumoniae in a well-studied bighorn sheep herd in the Wallowa Mountains of Oregon, USA. Using 10,605 SNPs generated using RADseq technology for 25 female bighorn sheep, we assessed genomic diversity metrics and employed family-based genome-wide association methodologies to understand variant association and genetic architecture underlying chronic carriage. We observed no differences among genome-wide diversity metrics (heterozygosity and allelic richness) between groups. However, we identified two variant loci of interest and seven associated candidate genes, which may influence carriage status. Further, we found that the SNP panel explained ~55% of the phenotypic variance (SNP-based heritability) for M. ovipneumoniae carriage, though there was considerable uncertainty in these estimates. While small sample sizes limit conclusions drawn here, our study represents one of the first to assess the genomic factors influencing chronic carriage of a pathogen in a wild population and lays a foundation for understanding genomic influence on pathogen persistence in bighorn sheep and other wildlife populations. Future research should incorporate additional individuals as well as distinct herds to further explore the genomic basis of chronic carriage.

Use of chlorisondamine to assess the neurogenic contribution to blood pressure in mice: An evaluation of method

Chlorisondamine (CSD) has been used to assess the neurogenic contribution to blood pressure (BP) and vasomotor sympathetic tone in animal models. It is assumed that the reduction in BP following CSD administration is associated to decreases in cardiac output (CO) and peripheral resistance, reflecting cardiac and vasomotor sympathetic tone, respectively. Surprisingly, this has not been characterized experimentally in mice, despite the extensive use of this animal model in cardiovascular research. We hypothesize that a specific dose of CSD can selectively block the sympathetic vasomotor tone. To test this hypothesis, we evaluated the effects of different doses of CSD (intraperitoneal) on BP and heart rate (HR) using telemetry, and on CO using echocardiography. BP and HR in normotensive C57Bl/6J mice reduced to a similar extent by all CSD doses tested (1-6 mg/kg). CSD at 6 mg/kg also reduced CO without affecting left ventricular stroke volume or fractional shortening. On the other hand, lower doses of CSD (1 and 2 mg/kg) produced significantly larger BP and HR reductions in DOCA-salt-induced hypertensive mice, indicating a greater neurogenic BP response. In addition, all doses of CSD reduced CO in hypertensive mice. Our data suggest that the BP response to CSD in mice likely reflects reduced CO and vasomotor sympathetic tone. We conclude that CSD can be used to assess the neurogenic contribution to BP in mice but may not be appropriate for specifically estimating vasomotor sympathetic tone.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Molecular docking-assisted screening reveals tannic acid as a natural protein disulphide isomerase inhibitor with antiplatelet and antithrombotic activities

Protein disulphide isomerase (PDI) promotes platelet activation and constitutes a novel antithrombotic target. In this study, we reported that a PDI‐binding plant polyphenol, tannic acid (TA), inhibits PDI activity, platelet activation and thrombus formation. Molecular docking using plant polyphenols from dietary sources with cardiovascular benefits revealed TA as the most potent binding molecule with PDI active centre. Surface plasmon resonance demonstrated that TA bound PDI with high affinity. Using Di‐eosin‐glutathione disulphide fluorescence assay and PDI assay kit, we showed that TA inhibited PDI activity. In isolated platelets, TA inhibited platelet aggregation stimulated by either GPVI or ITAM pathway agonists. Flow cytometry showed that TA inhibited thrombin‐ or CRP‐stimulated platelet activation, as reflected by reduced granule secretion and integrin activation. TA also reduced platelet spreading on immobilized fibrinogen and platelet adhesion under flow conditions. In a laser‐induced vascular injury mouse model, intraperitoneal injection of TA significantly decreased the size of cremaster arteriole thrombi. No prolongation of mouse jugular vein and tail‐bleeding time was observed after TA administration. Therefore, we identified TA from natural polyphenols as a novel inhibitor of PDI function. TA inhibits platelet activation and thrombus formation, suggesting it as a potential antithrombotic agent.

Inactivation of fission yeast Erh1 de-represses pho1 expression: evidence that Erh1 is a negative regulator of prt lncRNA termination

Fission yeast Erh1 exists in a complex with RNA-binding protein Mmi1. Deletion of erh1 up-regulates the phosphate homeostasis gene pho1, which is normally repressed by transcription in cis of a 5' flanking prt lncRNA. Here we present evidence that de-repression of pho1 by erh1Δ is achieved through precocious 3'-processing/termination of prt lncRNA synthesis, to wit: (i) erh1Δ does not affect the activity of the prt or pho1 promoters per se; (ii) de-repression by erh1Δ depends on CPF (cleavage and polyadenylation factor) subunits Ctf1, Dis2, Ssu72, Swd22, and Ppn1 and on termination factor Rhn1; (iii) de-repression requires synthesis by the Asp1 IPP kinase of inositol 1-pyrophosphates (1-IPPs); (iv) de-repression is effaced by mutating Thr4 of the RNA polymerase II CTD to alanine; and (v) erh1Δ exerts an additive effect on pho1 de-repression in combination with mutating CTD Ser7 to alanine and with deletion of the IPP pyrophosphatase Aps1. These findings point to Erh1 as an antagonist of lncRNA termination in the prt-pho1 axis. In contrast, in mmi1Δ cells there is a reduction in pho1 mRNA and increase in the formation of a prt-pho1 read-through transcript, consistent with Mmi1 being an agonist of prt termination. We envision that Erh1 acts as a brake on Mmi1's ability to promote CPF-dependent termination during prt lncRNA synthesis. Consistent with this idea, erh1Δ de-repression of pho1 was eliminated by mutating the Mmi1-binding sites in the prt lncRNA.

HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads

Complete and accurate genome assemblies form the basis of most downstream genomic analyses and are of critical importance. Recent genome assembly projects have relied on a combination of noisy long-read sequencing and accurate short-read sequencing, with the former offering greater assembly continuity and the latter providing higher consensus accuracy. The recently introduced Pacific Biosciences (PacBio) HiFi sequencing technology bridges this divide by delivering long reads (>10 kbp) with high per-base accuracy (>99.9%). Here we present HiCanu, a modification of the Canu assembler designed to leverage the full potential of HiFi reads via homopolymer compression, overlap-based error correction, and aggressive false overlap filtering. We benchmark HiCanu with a focus on the recovery of haplotype diversity, major histocompatibility complex (MHC) variants, satellite DNAs, and segmental duplications. For diploid human genomes sequenced to 30× HiFi coverage, HiCanu achieved superior accuracy and allele recovery compared to the current state of the art. On the effectively haploid CHM13 human cell line, HiCanu achieved an NG50 contig size of 77 Mbp with a per-base consensus accuracy of 99.999% (QV50), surpassing recent assemblies of high-coverage, ultralong Oxford Nanopore Technologies (ONT) reads in terms of both accuracy and continuity. This HiCanu assembly correctly resolves 337 out of 341 validation BACs sampled from known segmental duplications and provides the first preliminary assemblies of nine complete human centromeric regions. Although gaps and errors still remain within the most challenging regions of the genome, these results represent a significant advance toward the complete assembly of human genomes.

Nitric oxide and oxidative stress pathways do not contribute to sex differences in renal injury and function in Dahl SS/Jr rats

The burden of hypertension in the United States is increasing and yields significant morbidity and mortality, and sex differences in hypertension are widely recognized. Reduced nitric oxide (NO) bioavailability and increased oxidative stress are known to contribute to the pathogenesis of hypertensive renal injury, and but their contributions to sex differences in injury progression of are undefined. Our purpose was to test the hypothesis that male hypertensive rats have accelerated renal injury compared to females and to determine the contributions of the nitric oxide pathway and oxidative stress in these differences. Male and female Dahl SS/Jr rats, a model that spontaneously develops hypertension with age, were allowed to age on a 0.3% NaCl diet until 3 or 6 months of age, at which points blood pressure was measured and plasma, tissue, and urine were collected. While no significant sex differences in blood pressure were present at either time point, renal injury measured by urine protein excretion was more severe (male = 44.9 ± 6; female = 15±3 mg/day/100 g bw, p = .0001), and renal function was reduced (male = 0.48 ± 0.02; female = 0.7 ± 0.03 ml min-1  g-1 kw, p = .001) in males compared to females with age. Both male and female rats exhibited reduced nitric oxide metabolites (3 months: male = 0.65 ± 0.1; female = 0.74 ± 0.3; 6 months: male = 0.16 ± 0.1; female = 0.41 ± 0.1 ml min-1  g-1 kw, p, age = 0.02, p, sex = 0.3). Levels of urinary TBARS were similar (3 months: male = 20±1.5; female = 23±1.8; 6 months: male = 26±4.8; female = 23±4.7µM day g-1 kw, p, age = 0.4, p, sex = 0.9), extracellular superoxide dismutase (EC SOD) mRNA was greater in females (3 months: male = 0.35 ± 0.03; female = 1.4 ± 0.2; 6 months: male = 0.4 ± 0.05; female = 1.3 ± 0.1 normalized counts, p, age = 0.7, p, sex < 0.0001), but EC SOD protein expression was not different (3 months: male = 0.01 ± 0.002; female = 0.01 ± 0.002; 6 months: male = 0.02 ± 0.004; female = 0.01 ± 0.002 relative density, p, age = 0.2, p, sex = 0.8). These data support the presence of significant sex differences in renal injury and function in the Dahl S rat and identify a need for further study into the mechanisms involved.

Measuring the neutral zone of spinal motion segments: Comparison of multiple analysis methods to quantify spinal instability

PURPOSE

Neutral zone (NZ) parameters in spinal biomechanics studies are sensitive to spinal instability, disc degeneration, and repair. Multiple methods in the literature quantify NZ, yet no consensus exists on applicability and comparability of methods. This study compares five different NZ quantification methods using two different load-deflection profiles.

METHODS

Rat caudal and lumbar motion segments were tested in axial rotation to generate load-deflection curves with profiles exhibiting prominent distinction between elastic and NZ regions (ie, triphasic) and profiles that did not (ie, viscoelastic). NZ was quantified using five methods: trilinear, double sigmoid (DS), zero load, stiffness threshold (ST), and extrapolated elastic zone. Absolute agreement and consistency of NZ parameters were assessed using intraclass correlation (ICC), Bland-Altman analyses, and analysis of variance.

RESULTS

For triphasic profiles, NZ magnitude exhibited high consistency (methods correlate but differ in absolute values), and only some methods exhibited agreement. For viscoelastic profiles, NZ magnitude showed limited consistency and no absolute agreement. NZ stiffness had high agreement and consistency across most methods and profiles. For triphasic profiles, the linear NZ regions for all methods were not well-described by a linear fit yet for viscoelastic profiles all methods characterized a linear NZ region.

CONCLUSION

This NZ comparison study showed surprisingly limited agreement and consistency among NZ parameters with approximately 5% to 100% difference depending on the method and load-deflection profile. Nevertheless, the DS and ST methods appeared to be most comparable. We conclude that most NZ quantification methods cannot be applied interchangeably, highlighting a need to clearly state NZ calculation methods. Future studies are required to identify which methods are most sensitive to disc degeneration and repair in order to identify a "best" method.

Reproducibility of CRISPR-Cas9 methods for generation of conditional mouse alleles: a multi-center evaluation

BACKGROUND

CRISPR-Cas9 gene-editing technology has facilitated the generation of knockout mice, providing an alternative to cumbersome and time-consuming traditional embryonic stem cell-based methods. An earlier study reported up to 16% efficiency in generating conditional knockout (cKO or floxed) alleles by microinjection of 2 single guide RNAs (sgRNA) and 2 single-stranded oligonucleotides as donors (referred herein as "two-donor floxing" method).

RESULTS

We re-evaluate the two-donor method from a consortium of 20 laboratories across the world. The dataset constitutes 56 genetic loci, 17,887 zygotes, and 1718 live-born mice, of which only 15 (0.87%) mice contain cKO alleles. We subject the dataset to statistical analyses and a machine learning algorithm, which reveals that none of the factors analyzed was predictive for the success of this method. We test some of the newer methods that use one-donor DNA on 18 loci for which the two-donor approach failed to produce cKO alleles. We find that the one-donor methods are 10- to 20-fold more efficient than the two-donor approach.

CONCLUSION

We propose that the two-donor method lacks efficiency because it relies on two simultaneous recombination events in cis, an outcome that is dwarfed by pervasive accompanying undesired editing events. The methods that use one-donor DNA are fairly efficient as they rely on only one recombination event, and the probability of correct insertion of the donor cassette without unanticipated mutational events is much higher. Therefore, one-donor methods offer higher efficiencies for the routine generation of cKO animal models.

A piggyBac-based toolkit for inducible genome editing in mammalian cells

We describe the development and application of a novel series of vectors that facilitate CRISPR-Cas9-mediated genome editing in mammalian cells, which we call CRISPR-Bac. CRISPR-Bac leverages the piggyBac transposon to randomly insert CRISPR-Cas9 components into mammalian genomes. In CRISPR-Bac, a single piggyBac cargo vector containing a doxycycline-inducible Cas9 or catalytically dead Cas9 (dCas9) variant and a gene conferring resistance to Hygromycin B is cotransfected with a plasmid expressing the piggyBac transposase. A second cargo vector, expressing a single-guide RNA (sgRNA) of interest, the reverse-tetracycline TransActivator (rtTA), and a gene conferring resistance to G418, is also cotransfected. Subsequent selection on Hygromycin B and G418 generates polyclonal cell populations that stably express Cas9, rtTA, and the sgRNA(s) of interest. We show that CRISPR-Bac can be used to knock down proteins of interest, to create targeted genetic deletions with high efficiency, and to activate or repress transcription of protein-coding genes and an imprinted long noncoding RNA. The ratio of sgRNA-to-Cas9-to-transposase can be adjusted in transfections to alter the average number of cargo insertions into the genome. sgRNAs targeting multiple genes can be inserted in a single transfection. CRISPR-Bac is a versatile platform for genome editing that simplifies the generation of mammalian cells that stably express the CRISPR-Cas9 machinery.

The accessible chromatin landscape of the murine hippocampus at single-cell resolution

Here we present a comprehensive map of the accessible chromatin landscape of the mouse hippocampus at single-cell resolution. Substantial advances of this work include the optimization of a single-cell combinatorial indexing assay for transposase accessible chromatin (sci-ATAC-seq); a software suite, scitools, for the rapid processing and visualization of single-cell combinatorial indexing data sets; and a valuable resource of hippocampal regulatory networks at single-cell resolution. We used sci-ATAC-seq to produce 2346 high-quality single-cell chromatin accessibility maps with a mean unique read count per cell of 29,201 from both fresh and frozen hippocampi, observing little difference in accessibility patterns between the preparations. By using this data set, we identified eight distinct major clusters of cells representing both neuronal and nonneuronal cell types and characterized the driving regulatory factors and differentially accessible loci that define each cluster. Within pyramidal neurons, we identified four major clusters, including CA1 and CA3 neurons, and three additional subclusters. We then applied a recently described coaccessibility framework, Cicero, which identified 146,818 links between promoters and putative distal regulatory DNA. Identified coaccessibility networks showed cell-type specificity, shedding light on key dynamic loci that reconfigure to specify hippocampal cell lineages. Lastly, we performed an additional sci-ATAC-seq preparation from cultured hippocampal neurons (899 high-quality cells, 43,532 mean unique reads) that revealed substantial alterations in their epigenetic landscape compared with nuclei from hippocampal tissue. This data set and accompanying analysis tools provide a new resource that can guide subsequent studies of the hippocampus.

Rare variants of the FMN riboswitch class in Clostridium difficile and other bacteria exhibit altered ligand specificity

Many bacteria use flavin mononucleotide (FMN) riboswitches to control the expression of genes responsible for the biosynthesis and transport of this enzyme cofactor or its precursor, riboflavin. Rare variants of FMN riboswitches found in strains of Clostridium difficile and some other bacteria typically control the expression of proteins annotated as transporters, including multidrug efflux pumps. These RNAs no longer recognize FMN, and differ from the original riboswitch consensus sequence at nucleotide positions normally involved in binding of the ribityl and phosphate moieties of the cofactor. Representatives of one of the two variant subtypes were found to bind the FMN precursor riboflavin and the FMN degradation products lumiflavin and lumichrome. Although the biologically relevant ligand sensed by these variant FMN riboswitches remains uncertain, our findings suggest that many strains of C. difficile might use rare riboswitches to sense flavin degradation products and activate transporters for their detoxification.

Poly(A) polymerase is required for RyhB sRNA stability and function in Escherichia coli

Small regulatory RNAs (sRNAs) are an important class of bacterial post-transcriptional regulators that control numerous physiological processes, including stress responses. In Gram-negative bacteria including Escherichia coli, the RNA chaperone Hfq binds many sRNAs and facilitates pairing to target transcripts, resulting in changes in mRNA transcription, translation, or stability. Here, we report that poly(A) polymerase (PAP I), which promotes RNA degradation by exoribonucleases through the addition of poly(A) tails, has a crucial role in the regulation of gene expression by Hfq-dependent sRNAs. Specifically, we show that deletion of pcnB, encoding PAP I, paradoxically resulted in an increased turnover of certain Hfq-dependent sRNAs, including RyhB. RyhB instability in the pcnB deletion strain was suppressed by mutations in hfq or ryhB that disrupt pairing of RyhB with target RNAs, by mutations in the 3' external transcribed spacer of the glyW-cysT-leuZ transcript (3'ETS^LeuZ) involved in pairing with RyhB, or an internal deletion in rne, which encodes the endoribonuclease RNase E. Finally, the reduced stability of RyhB in the pcnB deletion strain resulted in impaired regulation of some of its target mRNAs, specifically sodB and sdhCDAB. Altogether our data support a model where PAP I plays a critical role in ensuring the efficient decay of the 3'ETS^LeuZ In the absence of PAP I, the 3'ETS^LeuZ transcripts accumulate, bind Hfq, and pair with RyhB, resulting in its depletion via RNase E-mediated decay. This ultimately leads to a defect in RyhB function in a PAP I deficient strain.

Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Alu elements, the short interspersed element numbering more than 1 million copies per human genome, can mediate the formation of copy number variants (CNVs) between substrate pairs. These Alu/Alu-mediated rearrangements (AAMRs) can result in pathogenic variants that cause diseases. To investigate the impact of AAMR on gene variation and human health, we first characterized Alus that are involved in mediating CNVs (CNV-Alus) and observed that these Alus tend to be evolutionarily younger. We then computationally generated, with the assistance of a supercomputer, a test data set consisting of 78 million Alu pairs and predicted ∼18% of them are potentially susceptible to AAMR. We further determined the relative risk of AAMR in 12,074 OMIM genes using the count of predicted CNV-Alu pairs and experimentally validated the predictions with 89 samples selected by correlating predicted hotspots with a database of CNVs identified by clinical chromosomal microarrays (CMAs) on the genomes of approximately 54,000 subjects. We fine-mapped 47 duplications, 40 deletions, and two complex rearrangements and examined a total of 52 breakpoint junctions of simple CNVs. Overall, 94% of the candidate breakpoints were at least partially Alu mediated. We successfully predicted all (100%) of Alu pairs that mediated deletions (n = 21) and achieved an 87% positive predictive value overall when including AAMR-generated deletions and duplications. We provided a tool, AluAluCNVpredictor, for assessing AAMR hotspots and their role in human disease. These results demonstrate the utility of our predictive model and provide insights into the genomic features and molecular mechanisms underlying AAMR.