Division of labor among Mycobacterium smegmatis RNase H enzymes: RNase H1 activity of RnhA or RnhC is essential for growth whereas RnhB and RnhA guard against killing by hydrogen peroxide in stationary phase

A comparative analysis of mycobacterial ribonucleases: Towards a therapeutic novel drug target

Bacterial responses to continuously changing environments are addressed through modulation of gene expression at the level of transcription initiation, RNA processing and/or decay. Ribonucleases (RNases) are hydrolytic or phosphorolytic enzymes involved in a majority of RNA metabolism reactions. RNases play a crucial role in RNA degradation, either independently or in collaboration with various trans-acting regulatory factors. The genus Mycobacterium consists of five subgenera: Mycobacteroides, Mycolicibacterium, Mycobacterium, Mycolicibacter and Mycolicibacillus, which include 63 fully sequenced species (pathogenic/non-pathogenic) to date. These include 13 different RNases, among which 5 are exonucleases (RNase PH, PNPase, RNase D, nano-RNases and RNase AS) and 8 are endonucleases (RNase J, RNase H, RNase P, RNase III, RNase BN, RNase Z, RNase G and RNase E), although RNase J and RNase BN were later identified to have exoribonuclease functions also. Here, we provide a detailed comparative insight into the Escherichia coli and mycobacterial RNases with respect to their types, phylogeny, structure, function, regulation and mechanism of action, with the main emphasis on RNase E. Among these 13 different mycobacterial RNases, 10 are essential for cell survival and have diverse structures hence, they are promising drug targets. RNase E is also an essential endonuclease that is abundant in many bacteria, forms an RNA degradosome complex that controls central RNA processing/degradation and has a conserved 5' sensor domain/DNase-I like region in its RNase domain. The essential mycobacterial RNases especially RNase E provide a potential repertoire of drug targets that can be exploited for inhibitor/modulator screening against many deadly mycobacterial diseases.

Comparative genomic analysis identifies potential adaptive variation in Mycoplasma ovipneumoniae

Mycoplasma ovipneumoniae is associated with respiratory disease in wild and domestic Caprinae globally, with wide variation in disease outcomes within and between host species. To gain insight into phylogenetic structure and mechanisms of pathogenicity for this bacterial species, we compared M. ovipneumoniae genomes for 99 samples from 6 countries (Australia, Bosnia and Herzegovina, Brazil, China, France and USA) and 4 host species (domestic sheep, domestic goats, bighorn sheep and caribou). Core genome sequences of M. ovipneumoniae assemblies from domestic sheep and goats fell into two well-supported phylogenetic clades that are divergent enough to be considered different bacterial species, consistent with each of these two clades having an evolutionary origin in separate host species. Genome assemblies from bighorn sheep and caribou also fell within these two clades, indicating multiple spillover events, most commonly from domestic sheep. Pangenome analysis indicated a high percentage (91.4 %) of accessory genes (i.e. genes found only in a subset of assemblies) compared to core genes (i.e. genes found in all assemblies), potentially indicating a propensity for this pathogen to adapt to within-host conditions. In addition, many genes related to carbon metabolism, which is a virulence factor for Mycoplasmas, showed evidence for homologous recombination, a potential signature of adaptation. The presence or absence of annotated genes was very similar between sheep and goat clades, with only two annotated genes significantly clade-associated. However, three M. ovipneumoniae genome assemblies from asymptomatic caribou in Alaska formed a highly divergent subclade within the sheep clade that lacked 23 annotated genes compared to other assemblies, and many of these genes had functions related to carbon metabolism. Overall, our results suggest that adaptation of M. ovipneumoniae has involved evolution of carbon metabolism pathways and virulence mechanisms related to those pathways. The genes involved in these pathways, along with other genes identified as potentially involved in virulence in this study, are potential targets for future investigation into a possible genomic basis for the high variation observed in disease outcomes within and between wild and domestic host species.

Precise measurement of molecular phenotypes with barcode-based CRISPRi systems

Genome-wide CRISPR-Cas9 screens have untangled regulatory networks and revealed the genetic underpinnings of diverse biological processes. Their success relies on experimental designs that interrogate specific molecular phenotypes and distinguish key regulators from background effects. Here, we realize these goals with a generalizable platform for CRISPR interference with barcoded expression reporter sequencing (CiBER-seq) that dramatically improves the sensitivity and scope of genome-wide screens. We systematically address technical factors that distort phenotypic measurements by normalizing expression reporters against closely-matched control promoters, integrated together into the genome at single copy. To test our ability to capture post-transcriptional and post-translational regulation through sequencing, we screened for genes that affected nonsense-mediated mRNA decay and Doa10-mediated cytosolic protein decay. Our optimized CiBER-seq screens accurately capture the known components of well-studied RNA and protein quality control pathways with minimal background. These results demonstrate the precision and versatility of CiBER-seq for dissecting the genetic networks controlling cellular behaviors.

Spectroscopic approaches for studies of site-specific DNA base and backbone 'breathing' using exciton-coupled dimer-labeled DNA

DNA regulation and repair processes require direct interactions between proteins and DNA at specific sites. Local fluctuations of the sugar-phosphate backbones and bases of DNA (a form of DNA 'breathing') play a central role in such processes. Here we review the development and application of novel spectroscopic methods and analyses - both at the ensemble and single-molecule levels - to study structural and dynamic properties of exciton-coupled cyanine and fluorescent nucleobase analogue dimer-labeled DNA constructs at key positions involved in protein-DNA complex assembly and function. The exciton-coupled dimer probes act as 'sensors' of the local conformations adopted by the sugar-phosphate backbones and bases immediately surrounding the dimer probes. These methods can be used to study the mechanisms of protein binding and function at these sites.

Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants

Many membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for the mammalian gonadotropin-releasing hormone receptor GPCRs (GnRHR). We recently demonstrated that evolutionary GnRHR modifications appear to have coincided with adaptive changes in cotranslational folding efficiency. Though protein stability is known to shape evolution, it is unclear how cotranslational folding constraints modulate the synergistic, epistatic interactions between mutations. We therefore compared the pairwise interactions formed by mutations that disrupt the membrane topology (V276T) or tertiary structure (W107A) of GnRHR. Using deep mutational scanning, we evaluated how the plasma membrane expression of these variants is modified by hundreds of secondary mutations. An analysis of 251 mutants in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the severity and the mechanism of destabilization. V276T forms predominantly negative epistatic interactions with destabilizing mutations in soluble loops. In contrast, W107A forms positive interactions with mutations in both loops and transmembrane domains that reflect the diminishing impacts of the destabilizing mutations in variants that are already unstable. These findings reveal how epistasis is remodeled by conformational defects in membrane proteins and in unstable proteins more generally.

Structure of the M. tuberculosis DnaK-GrpE complex reveals how key DnaK roles are controlled

The molecular chaperone DnaK is essential for viability of Mycobacterium tuberculosis (Mtb). DnaK hydrolyzes ATP to fold substrates, and the resulting ADP is exchanged for ATP by the nucleotide exchange factor GrpE. It has been unclear how GrpE couples DnaK's nucleotide exchange with substrate release. Here we report a cryo-EM analysis of GrpE bound to an intact Mtb DnaK, revealing an asymmetric 1:2 DnaK-GrpE complex. The GrpE dimer ratchets to modulate both DnaK nucleotide-binding domain and the substrate-binding domain. We further show that the disordered GrpE N-terminus is critical for substrate release, and that the DnaK-GrpE interface is essential for protein folding activity both in vitro and in vivo. Therefore, the Mtb GrpE dimer allosterically regulates DnaK to concomitantly release ADP in the nucleotide-binding domain and substrate peptide in the substrate-binding domain.

Differences in white blood cell proportions between schizophrenia cases and controls are influenced by medication and variations in time of day

Cases with schizophrenia (SCZ) and healthy controls show differences in white blood cell (WBC) counts and blood inflammation markers. Here, we investigate whether time of blood draw and treatment with psychiatric medications are related to differences in estimated WBC proportions between SCZ cases and controls. DNA methylation data from whole blood was used to estimate proportions of six subtypes of WBCs in SCZ patients (n = 333) and healthy controls (n = 396). We tested the association of case-control status with estimated cell-type proportions and the neutrophil-to-lymphocyte ratio (NLR) in 4 models: with/without adjusting for time of blood draw, and then compared results from blood samples drawn during a 12-h (07:00-19:00) or 7-h (07:00-14:00) period. We also investigated WBC proportions in a subgroup of medication-free patients (n = 51). Neutrophil proportions were significantly higher in SCZ cases (mean=54.1%) vs. controls (mean=51.1%; p = <0.001), and CD8⁺T lymphocyte proportions were lower in SCZ cases (mean=12.1%) vs. controls (mean=13.2%; p = 0.001). The effect sizes in the 12-h sample (07:00-19:00) showed a significant difference between SCZ vs. controls for neutrophils, CD4⁺T, CD8⁺T, and B-cells, which remained significant after adjusting for time of blood draw. In the samples matched for time of blood draw during 07.00-14.00, we also observed an association with neutrophils, CD4⁺T, CD8⁺T, and B-cells that was unaffected by further adjustment for time of blood draw. In the medication-free patients, we observed differences that remained significant in neutrophils (p = 0.01) and CD4⁺T (p = 0.01) after adjusting for time of day. The association of SCZ with NLR was significant in all models (range: p < 0.001 to p = 0.03) in both medicated and unmedicated patients. In conclusion, controlling for pharmacological treatment and circadian cycling of WBC is necessary for unbiased estimates in case-control studies. Nevertheless, the association of WBC with SCZ remains, even after adjusting for the time of day.

The C-Terminal Acid Phosphatase Module of the RNase HI Enzyme RnhC Controls Rifampin Sensitivity and Light-Dependent Colony Pigmentation of Mycobacterium smegmatis

RNase H enzymes participate in various processes that require processing of RNA-DNA hybrids, including DNA replication, transcription, and ribonucleotide excision from DNA. Mycobacteria encode multiple RNase H enzymes, and prior data indicate that RNase HI activity is essential for mycobacterial viability. However, the additional roles of mycobacterial RNase Hs are unknown, including whether RNase HII (RnhB and RnhD) excises chromosomal ribonucleotides misincorporated during DNA replication and whether individual RNase HI enzymes (RnhA and RnhC) mediate additional phenotypes. We find that loss of RNase HII activity in Mycobacterium smegmatis (through combined deletion of rnhB/rnhD) or individual RNase HI enzymes does not affect growth, hydroxyurea sensitivity, or mutagenesis, whereas overexpression (OE) of either RNase HII severely compromises bacterial viability. We also show that deletion of rnhC, which encodes a protein with an N-terminal RNase HI domain and a C-terminal acid phosphatase domain, confers sensitivity to rifampin and oxidative stress as well as loss of light-induced carotenoid pigmentation. These phenotypes are due to loss of the activity of the C-terminal acid phosphatase domain rather than the RNase HI activity, suggesting that the acid phosphatase activity may confer rifampin resistance through the antioxidant properties of carotenoid pigment production. IMPORTANCE Mycobacteria encode multiple RNase H enzymes, with RNase HI being essential for viability. Here, we examine additional functions of RNase H enzymes in mycobacteria. We find that RNase HII is not involved in mutagenesis but is highly toxic when overexpressed. The RNase HI enzyme RnhC is required for tolerance to rifampin, but this role is surprisingly independent of its RNase H activity and is instead mediated by an autonomous C-terminal acid phosphatase domain. This study provides new insights into the functions of the multiple RNase H enzymes of mycobacteria.

PolyAMiner-Bulk: A Machine Learning Based Bioinformatics Algorithm to Infer and Decode Alternative Polyadenylation Dynamics from bulk RNA-seq data

More than half of human genes exercise alternative polyadenylation (APA) and generate mRNA transcripts with varying 3' untranslated regions (UTR). However, current computational approaches for identifying cleavage and polyadenylation sites (C/PASs) and quantifying 3'UTR length changes from bulk RNA-seq data fail to unravel tissue- and disease-specific APA dynamics. Here, we developed a next-generation bioinformatics algorithm and application, PolyAMiner-Bulk, that utilizes an attention-based machine learning architecture and an improved vector projection-based engine to infer differential APA dynamics accurately. When applied to earlier studies, PolyAMiner-Bulk accurately identified more than twice the number of APA changes in an RBM17 knockdown bulk RNA-seq dataset compared to current generation tools. Moreover, on a separate dataset, PolyAMiner-Bulk revealed novel APA dynamics and pathways in scleroderma pathology and identified differential APA in a gene that was identified as being involved in scleroderma pathogenesis in an independent study. Lastly, we used PolyAMiner-Bulk to analyze the RNA-seq data of post-mortem prefrontal cortexes from the ROSMAP data consortium and unraveled novel APA dynamics in Alzheimer's Disease. Our method, PolyAMiner-Bulk, creates a paradigm for future alternative polyadenylation analysis from bulk RNA-seq data.

A comprehensive method to study the DNA's association with lamin and chromatin compaction in intact cell nuclei at super resolution

Super-resolution fluorescence microscopy has revolutionized multicolor imaging of nuclear structures due to the combination of high labeling specificity and high resolution. Here we expanded the recently developed fBALM (DNA structure fluctuation-assisted binding activated localization microscopy) method by developing a stable methodological sequence that enables dual-color imaging of high-resolution genomic DNA together with an immunofluorescently labeled intranuclear protein. Our measurements of the nuclear periphery, imaging DNA and LaminB1 in biologically relevant samples, show that this novel dual-color imaging method is feasible for further quantitative evaluations. We were able to study the relative spatial signal organization between DNA and LaminB1 by means of highly specific colocalization measurements at nanometer resolution. Measurements were performed with and without the antifade embedding medium ProLong Gold, which proved to be essential for imaging of LaminB1, but not for imaging of SytoxOrange labeled DNA. The localization precision was used to differentiate between localizations with higher and lower amounts of emitting photons. We interpret high intensity localizations to be renatured DNA sections in which a high amount of Sytox Orange molecules were bound. This could give insight into the denaturation kinetics of DNA during fBALM. These results were further complemented by measurements of γH2AX and H3K9me3 signal organization to demonstrate differences within the chromatin landscape, which were quantified with image processing methods such as Voronoi segmentation.

The catalytic mechanism, metal dependence, substrate specificity, and biodiversity of ribonuclease H

Ribonucleoside monophosphates are inevitably misincorporated into the DNA genome inside cells, and they need to be excised to avoid chromosome instability. Ribonucleases H (RNases H) are enzymes that specifically hydrolyze the RNA strand of RNA/DNA hybrids or the RNA moiety from DNA containing a stretch of RNA, they therefore are required for DNA integrity. Extensive studies have drawn a mostly clear picture of the mechanisms of RNase H catalysis, but some questions are still lacking definitive answers. This review summarizes three alternative models of RNase H catalysis. The two-metal model is prevalent, but a three-metal model suggests the involvement of a third cation in catalysis. Apparently, the mechanisms underlying metal-dependent hydrolyzation are more complicated than initially thought. We also discuss the metal choices of RNases H and analyze how chemically similar cations function differently. Substrate and cleavage-site specificities vary among RNases H, and this is explicated in detail. An intriguing phenomenon is that organisms have diverse RNase H combinations, which may provide important hints to how rnh genes were transferred during evolution. Whether RNase H is essential for cellular growth, a key question in the study of in vivo functions, is also discussed. This article may aid in understanding the mechanisms underlying RNase H and in developing potentially promising applications of it.

RNase HI Depletion Strongly Potentiates Cell Killing by Rifampicin in Mycobacteria

Multidrug-resistant (MDR) tuberculosis (TB) is defined by the resistance of Mycobacterium tuberculosis, the causative organism, to the first-line antibiotics rifampicin and isoniazid. Mitigating or reversing resistance to these drugs offers a means of preserving and extending their use in TB treatment. R-loops are RNA/DNA hybrids that are formed in the genome during transcription, and they can be lethal to the cell if not resolved. RNase HI is an enzyme that removes R-loops, and this activity is essential in M. tuberculosis: knockouts of rnhC, the gene encoding RNase HI, are nonviable. This essentiality makes it a candidate target for the development of new antibiotics. In the model organism Mycolicibacterium smegmatis, RNase HI activity is provided by two enzymes, RnhA and RnhC. We show that the partial depletion of RNase HI activity in M. smegmatis, by knocking out either of the genes encoding RnhA or RnhC, led to the accumulation of R-loops. The sensitivity of the knockout strains to the antibiotics moxifloxacin, streptomycin, and rifampicin was increased, the latter by a striking near 100-fold. We also show that R-loop accumulation accompanies partial transcriptional inhibition, suggesting a mechanistic basis for the synergy between RNase HI depletion and rifampicin. A model of how transcriptional inhibition can potentiate R-loop accumulation is presented. Finally, we identified four small molecules that inhibit recombinant RnhC activity and that also potentiated rifampicin activity in whole-cell assays against M. tuberculosis, supporting an on-target mode of action and providing the first step in developing a new class of antimycobacterial drug.

Cabin1 domain-containing gene picd-1 interacts with pry-1/Axin to regulate multiple processes in Caenorhabditis elegans

The Axin family of scaffolding proteins control diverse processes, such as facilitating the interactions between cellular components and providing specificity to signaling pathways. While several Axin family members have been discovered in metazoans and shown to play crucial roles, their mechanism of action are not well understood. The Caenorhabditis elegans Axin homolog, pry-1, is a powerful tool for identifying interacting genes and downstream effectors that function in a conserved manner to regulate Axin-mediated signaling. Our lab and others have established pry-1's essential role in developmental processes that affect the reproductive system, seam cells, and a posterior P lineage cell, P11.p. Additionally, pry-1 is crucial for lipid metabolism, stress responses, and aging. In this study, we expanded on our previous work on pry-1 by reporting a novel interacting gene named picd-1 (pry-1-interacting and Cabin1 domain-containing). PICD-1 protein shares sequence conservation with CABIN1, a component of the HUCA complex. Our findings have revealed that PICD-1 is involved in several pry-1-mediated processes, including stress response and lifespan maintenance. picd-1's expression overlapped with that of pry-1 in multiple tissues throughout the lifespan. Furthermore, PRY-1 and PICD-1 inhibited CREB-regulated transcriptional coactivator homolog CRTC-1, which promotes longevity in a calcineurin-dependent manner. Overall, our study has demonstrated that picd-1 is necessary for mediating pry-1 function and provides the basis to investigate whether Cabin-1 domain-containing protein plays a similar role in Axin signaling in other systems.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

An AraC/XylS Family Transcriptional Regulator Modulates the Oxidative Stress Response of Francisella tularensis

Francisella tularensis is a Gram-negative bacterium that causes a fatal human disease known as tularemia. The Centers for Disease Control and Prevention have classified F. tularensis as a category A tier 1 select agent. The virulence mechanisms of Francisella are not entirely understood. Francisella possesses very few transcription regulators, and most of these regulate the expression of genes involved in intracellular survival and virulence. The F. tularensis genome sequence analysis reveals an AraC (FTL_0689) transcriptional regulator homologous to the AraC/XylS family of transcriptional regulators. In Gram-negative bacteria, AraC activates genes required for l-arabinose utilization and catabolism. The role of the FTL_0689 regulator in F. tularensis is not known. In this study, we characterized the role of FTL_0689 in the gene regulation of F. tularensis and investigated its contribution to intracellular survival and virulence. The results demonstrate that FTL_0689 in Francisella is not required for l-arabinose utilization. Instead, FTL_0689 specifically regulates the expression of the oxidative and global stress response, virulence, metabolism, and other key pathways genes required by Francisella when exposed to oxidative stress. The FTL_0689 mutant is attenuated for intramacrophage growth and virulence in mice. Based on the deletion mutant phenotype, FTL_0689 was termed osrR (oxidative stress response regulator). Altogether, this study elucidates the role of the osrR transcriptional regulator in tularemia pathogenesis. IMPORTANCE The virulence mechanisms of category A select agent Francisella tularensis, the causative agent of a fatal human disease known as tularemia, remain largely undefined. The present study investigated the role of a transcriptional regulator and its overall contribution to the oxidative stress resistance of F. tularensis. The results provide an insight into a novel gene regulatory mechanism, especially when Francisella is exposed to oxidative stress conditions. Understanding such Francisella- specific regulatory mechanisms will help identify potential targets for developing effective therapies and vaccines to prevent tularemia.

Bacillus subtilis PcrA Helicase Removes Trafficking Barriers

Bacillus subtilis PcrA interacts with the RNA polymerase and might contribute to mitigate replication-transcription conflicts (RTCs). We show that PcrA depletion lethality is partially suppressed by rnhB inactivation, but cell viability is significantly reduced by rnhC or dinG inactivation. Following PcrA depletion, cells lacking RnhC or DinG are extremely sensitive to DNA damage. Chromosome segregation is not further impaired by rnhB or dinG inactivation but is blocked by rnhC or recA inactivation upon PcrA depletion. Despite our efforts, we could not construct a ΔrnhC ΔrecA strain. These observations support the idea that PcrA dismantles RTCs. Purified PcrA, which binds single-stranded (ss) DNA over RNA, is a ssDNA-dependent ATPase and preferentially unwinds DNA in a 3'→5'direction. PcrA unwinds a 3'-tailed RNA of an RNA-DNA hybrid significantly faster than that of a DNA substrate. Our results suggest that a replicative stress, caused by mis-incorporated rNMPs, indirectly increases cell viability upon PcrA depletion. We propose that PcrA, in concert with RnhC or DinG, contributes to removing spontaneous or enzyme-driven R-loops, to counteract deleterious trafficking conflicts and preserve to genomic integrity.

DNA Breaks-Mediated Fitness Cost Reveals RNase HI as a New Target for Selectively Eliminating Antibiotic-Resistant Bacteria

Antibiotic resistance often generates defects in bacterial growth called fitness cost. Understanding the causes of this cost is of paramount importance, as it is one of the main determinants of the prevalence of resistances upon reducing antibiotics use. Here we show that the fitness costs of antibiotic resistance mutations that affect transcription and translation in Escherichia coli strongly correlate with DNA breaks, which are generated via transcription–translation uncoupling, increased formation of RNA–DNA hybrids (R-loops), and elevated replication–transcription conflicts. We also demonstrated that the mechanisms generating DNA breaks are repeatedly targeted by compensatory evolution, and that DNA breaks and the cost of resistance can be increased by targeting the RNase HI, which specifically degrades R-loops. We further show that the DNA damage and thus the fitness cost caused by lack of RNase HI function drive resistant clones to extinction in populations with high initial frequency of resistance, both in laboratory conditions and in a mouse model of gut colonization. Thus, RNase HI provides a target specific against resistant bacteria, which we validate using a repurposed drug. In summary, we revealed key mechanisms underlying the fitness cost of antibiotic resistance mutations that can be exploited to specifically eliminate resistant bacteria.

Bisulfite-free epigenomics and genomics of single cells through methylation-sensitive restriction

Single-cell multi-omics are powerful means to study cell-to-cell heterogeneity. Here, we present a single-tube, bisulfite-free method for the simultaneous, genome-wide analysis of DNA methylation and genetic variants in single cells: epigenomics and genomics of single cells analyzed by restriction (epi-gSCAR). By applying this method, we obtained DNA methylation measurements of up to 506,063 CpGs and up to 1,244,188 single-nucleotide variants from single acute myeloid leukemia-derived cells. We demonstrate that epi-gSCAR generates accurate and reproducible measurements of DNA methylation and allows to differentiate between cell lines based on the DNA methylation and genetic profiles.

RnhP is a plasmid-borne RNase HI that contributes to genome maintenance in the ancestral strain Bacillus subtilis NCIB 3610

RNA-DNA hybrids form throughout the chromosome during normal growth and under stress conditions. When left unresolved, RNA-DNA hybrids can slow replication fork progression, cause DNA breaks, and increase mutagenesis. To remove hybrids, all organisms use ribonuclease H (RNase H) to specifically degrade the RNA portion. Here we show that, in addition to chromosomally encoded RNase HII and RNase HIII, Bacillus subtilis NCIB 3610 encodes a previously uncharacterized RNase HI protein, RnhP, on the endogenous plasmid pBS32. Like other RNase HI enzymes, RnhP incises Okazaki fragments, ribopatches, and a complementary RNA-DNA hybrid. We show that while chromosomally encoded RNase HIII is required for pBS32 hyper-replication, RnhP compensates for the loss of RNase HIII activity on the chromosome. Consequently, loss of RnhP and RNase HIII impairs bacterial growth. We show that the decreased growth rate can be explained by laggard replication fork progression near the terminus region of the right replichore, resulting in SOS induction and inhibition of cell division. We conclude that all three functional RNase H enzymes are present in B. subtilis NCIB 3610 and that the plasmid-encoded RNase HI contributes to chromosome stability, while the chromosomally encoded RNase HIII is important for chromosome stability and plasmid hyper-replication.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

Pleistocene population dynamics can inform the consequences of current climate change. This phylogeography of 35 complete American mastodon mitochondrial genomes suggests distinct lineages in this species repeatedly expanded northwards and then went locally extinct in response to glacial cycles.

Functional Disassociation Between the Protein Domains of MSMEG_4305 of Mycolicibacterium smegmatis (Mycobacterium smegmatis) in vivo

MSMEG_4305 is a two-domain protein of Mycolicibacterium smegmatis (Mycobacterium smegmatis) (Mycolicibacterium smegmatis). The N-terminal domain of MSMEG_4305 encodes an RNase H type I. The C-terminal domain is a presumed CobC, predicted to be involved in the aerobic synthesis of vitamin B12. Both domains reach their maximum at distinct pH, approximately 8.5 and 4.5, respectively. The presence of the CobC domain influenced RNase activity in vitro in homolog Rv2228c. Here, we analyzed the role of MSMEG_4305 in vitamin B12 synthesis and the functional association between both domains in vivo in M. smegmatis. We used knock-out mutant of M. smegmatis, deficient in MSMEG_4305. Whole-cell lysates of the mutants strain contained a lower concentration of vitamin B12, as it determined with immunoenzimatic assay. We observed growth deficits, related to vitamin B12 production, on media containing sulfamethazine and propionate. Removal of the CobC domain of MSMEG_4305 in ΔrnhA background hardly affected the growth rate of M. smegmatis in vivo. The strain carrying truncation showed no fitness deficit in the competitive assay and it did not show increased level of RNA/DNA hybrids in its genome. We show that homologs of MSMEG_4305 are present only in the Actinomycetales phylogenetic branch (according to the old classification system). The domains of MSMEG_4305 homologs accumulate mutations at a different rate, while the linker region is highly variable. We conclude that MSMEG_4305 is a multidomain protein that most probably was fixed in the phylogenetic tree of life due to genetic drift.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

Mycobacterial DNA polymerase I: activities and crystal structures of the POL domain as apoenzyme and in complex with a DNA primer-template and of the full-length FEN/EXO-POL enzyme

Mycobacterial Pol1 is a bifunctional enzyme composed of an N-terminal DNA flap endonuclease/5' exonuclease domain (FEN/EXO) and a C-terminal DNA polymerase domain (POL). Here we document additional functions of Pol1: FEN activity on the flap RNA strand of an RNA:DNA hybrid and reverse transcriptase activity on a DNA-primed RNA template. We report crystal structures of the POL domain, as apoenzyme and as ternary complex with 3'-dideoxy-terminated DNA primer-template and dNTP. The thumb, palm, and fingers subdomains of POL form an extensive interface with the primer-template and the triphosphate of the incoming dNTP. Progression from an open conformation of the apoenzyme to a nearly closed conformation of the ternary complex entails a disordered-to-ordered transition of several segments of the thumb and fingers modules and an inward motion of the fingers subdomain-especially the O helix-to engage the primer-template and dNTP triphosphate. Distinctive structural features of mycobacterial Pol1 POL include a manganese binding site in the vestigial 3' exonuclease subdomain and a non-catalytic water-bridged magnesium complex at the protein-DNA interface. We report a crystal structure of the bifunctional FEN/EXO-POL apoenzyme that reveals the positions of two active site metals in the FEN/EXO domain.

Multiplexed CRISPR technologies for gene editing and transcriptional regulation

Multiplexed CRISPR technologies, in which numerous gRNAs or Cas enzymes are expressed at once, have facilitated powerful biological engineering applications, vastly enhancing the scope and efficiencies of genetic editing and transcriptional regulation. In this review, we discuss multiplexed CRISPR technologies and describe methods for the assembly, expression and processing of synthetic guide RNA arrays in vivo. Applications that benefit from multiplexed CRISPR technologies, including cellular recorders, genetic circuits, biosensors, combinatorial genetic perturbations, large-scale genome engineering and the rewiring of metabolic pathways, are highlighted. We also offer a glimpse of emerging challenges and emphasize experimental considerations for future studies.

The Rhinella arenarum transcriptome: de novo assembly, annotation and gene prediction

The common toad Rhinella arenarum is widely distributed in Argentina, where it is utilised as an autochthonous model in ecotoxicological research and environmental toxicology. However, the lack of a reference genome makes molecular assays and gene expression studies difficult to carry out on this non-model species. To address this issue, we performed a genome-wide transcriptome analysis on R. arenarum larvae through massive RNA sequencing, followed by de novo assembly, annotation, and gene prediction. We obtained 57,407 well-annotated transcripts representing 99.4% of transcriptome completeness (available at http://rhinella.uncoma.edu.ar). We also defined a set of 52,800 high-confidence lncRNA transcripts and demonstrated the reliability of the transcriptome data to perform phylogenetic analysis. Our comprehensive transcriptome analysis of R. arenarum represents a valuable resource to perform functional genomic studies and to identify potential molecular biomarkers in ecotoxicological research.

Complete mitochondrial genome of Sesia siningensis (Lepidoptera: Sesiidae)

Sesia siningensis is an important trunk borer of poplar and is widely distributed in China. Here, the complete mitochondrial genome of S. siningensis was sequenced. The circle genome of the clearwing moth is 15,454 bp in length. There are 38 sequence elements including 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, and a control region. The order of most elements was consistent with that of Chilo suppressalis, with the exception of one tRNA gene. As the first reported mitochondrial genome in the Sesiidae family, it will provide useful information to the development and application of better markers and primers in the molecular taxonomy of this family.

RNases H: Structure and mechanism

RNases H are a family of endonucleases that hydrolyze RNA residues in various nucleic acids. These enzymes are present in all branches of life, and their counterpart domains are also found in reverse transcriptases (RTs) from retroviruses and retroelements. RNases H are divided into two main classes (RNases H1 and H2 or type 1 and type 2 enzymes) with common structural features of the catalytic domain but different range of substrates for enzymatic cleavage. Additionally, a third class is found in some Archaea and bacteria. Besides distinct cellular functions specific for each type of RNases H, this family of proteins is generally involved in the maintenance of genome stability with overlapping and cooperative role in removal of R-loops thus preventing their accumulation. Extensive biochemical and structural studies of RNases H provided not only a comprehensive and complete picture of their mechanism but also revealed key basic principles of nucleic acid recognition and processing. RNase H1 is present in prokaryotes and eukaryotes and cleaves RNA in RNA/DNA hybrids. Its main function is hybrid removal, notably in the context of R-loops. RNase H2, which is also present in all branches of life, can play a similar role but it also has a specialized function in the cleavage of single ribonucleotides embedded in the DNA. RNase H3 is present in Archaea and bacteria and is closely related to RNase H2 in sequence and structure but has RNase H1-like biochemical properties. This review summarizes the mechanisms of substrate recognition and enzymatic cleavage by different classes of RNases H with particular insights into structural features of nucleic acid binding, specificity towards RNA and/or DNA strands and catalysis.

The interaction of lncRNA EZR-AS1 with SMYD3 maintains overexpression of EZR in ESCC cells

EZR, a member of the ezrin-radixin-moesin (ERM) family, is involved in multiple aspects of cell migration and cancer. SMYD3, a histone H3–lysine 4 (H3–K4)-specific methyltransferase, regulates EZR gene transcription, but the molecular mechanisms of epigenetic regulation remain ill-defined. Here, we show that antisense lncRNA EZR-AS1 was positively correlated with EZR expression in both human esophageal squamous cell carcinoma (ESCC) tissues and cell lines. Both in vivo and in vitro studies revealed that EZR-AS1 promoted cell migration through up-regulation of EZR expression. Mechanistically, antisense lncRNA EZR-AS1 formed a complex with RNA polymerase II to activate the transcription of EZR. Moreover, EZR-AS1 could recruit SMYD3 to a binding site, present in a GC-rich region downstream of the EZR promoter, causing the binding of SMYD3 and local enrichment of H3K4me3. Finally, the interaction of EZR-AS1 with SMYD3 further enhanced EZR transcription and expression. Our findings suggest that antisense lncRNA EZR-AS1, as a member of an RNA polymerase complex and through enhanced SMYD3-dependent H3K4 methylation, plays an important role in enhancing transcription of the EZR gene to promote the mobility and invasiveness of human cancer cells.

Precisely measured protein lifetimes in the mouse brain reveal differences across tissues and subcellular fractions

The turnover of brain proteins is critical for organism survival, and its perturbations are linked to pathology. Nevertheless, protein lifetimes have been difficult to obtain in vivo. They are readily measured in vitro by feeding cells with isotopically labeled amino acids, followed by mass spectrometry analyses. In vivo proteins are generated from at least two sources: labeled amino acids from the diet, and non-labeled amino acids from the degradation of pre-existing proteins. This renders measurements difficult. Here we solved this problem rigorously with a workflow that combines mouse in vivo isotopic labeling, mass spectrometry, and mathematical modeling. We also established several independent approaches to test and validate the results. This enabled us to measure the accurate lifetimes of ~3500 brain proteins. The high precision of our data provided a large set of biologically significant observations, including pathway-, organelle-, organ-, or cell-specific effects, along with a comprehensive catalog of extremely long-lived proteins (ELLPs).

Measuring precise protein turnover rates in animals is technically challenging at the proteomic level. Here, Fornasiero and colleagues use isotopic labeling with mass spectrometry and mathematical modeling to accurately determine protein lifetimes in the mouse brain

A tRNA-based multiplex sgRNA expression system in zebrafish and its application to generation of transgenic albino fish

The CRISPR/Cas9 system can be introduced into zebrafish as transgenes. Namely, expression of single-guide RNA (sgRNA) and controlled expression of Cas9 in transgenic zebrafish enables the study of gene functions in specific cell types. This transgenic CRISPR/Cas9 approach would be more useful if multiple sgRNAs could be expressed simultaneously since we could knock-out a gene more efficiently or disrupt multiple genes simultaneously. Here we describe a novel system to express multiple sgRNAs efficiently in zebrafish, that relies on the endogenous tRNA processing machinery. We cloned nine endogenous zebrafish tRNA genes, fused them to sgRNAs, and demonstrated that an active sgRNA can be produced from a precursor transcript containing either of these tRNAs. To show a proof of principle, we constructed transgenic fish expressing Cas9 under the control of the ubiquitin promoter and a single transcript containing three distinct sgRNAs, that targeted the slc45a2 (albino) gene, fused to tRNAs under the control of the U6 promoter. We found that the Tg(ubb:SpCas9,u6c:3xslc45a2-sgRNA) harbored mutations in all of the target sites in the albino gene and showed nearly complete albino phenotypes, which were amenable to imaging experiments. Thus, the tRNA-based multiplex sgRNA expression system should facilitate gene knock-out studies in transgenic zebrafish.

Reactivity and Specificity of RNase T1, RNase A, and RNase H toward Oligonucleotides of RNA Containing 8-Oxo-7,8-dihydroguanosine

Understanding how oxidatively damaged RNA interacts with ribonucleases is important because of its proposed role in the development and progression of disease. Thus, understanding structural aspects of RNA containing lesions generated under oxidative stress, as well as its interactions with other biopolymers, is fundamental. We explored the reactivity of RNase A, RNase T₁, and RNase H toward oligonucleotides of RNA containing 8-oxo-7,8-dihydroguanosine (8oxoG). This is the first example that addresses this relationship and will be useful for understanding (1) how these RNases can be used to characterize the structural impact that this lesion has on RNA and (2) how oxidatively modified RNA may be handled intracellularly. 8-OxoG was incorporated into 10-16-mers of RNA, and its reactivity with each ribonuclease was assessed via electrophoretic analyses, circular dichroism, and the use of other C8-purine-modified analogues (8-bromoguanosine, 8-methoxyguanosine, and 8-oxoadenosine). RNase T₁ does not recognize sites containing 8-oxoG, while RNase A recognizes and cleaves RNA at positions containing this lesion while differentiating if it is involved in H-bonding. The selectivity of RNase A followed the order C > 8-oxoG ≈ U. In addition, isothermal titration calorimetry showed that an 8-oxoG-C3'-methylphosphate derivative can inhibit RNase A activity. Cleavage patterns obtained from RNase H displayed changes in reactivity in a sequence- and concentration-dependent manner and displayed recognition at sites containing the modification in some cases. These data will aid in understanding how this modification affects reactivity with ribonucleases and will enable the characterization of global and local structural changes in oxidatively damaged RNA.

Guardians of the mycobacterial genome: A review on DNA repair systems in Mycobacterium tuberculosis

The genomic integrity of Mycobacterium tuberculosis is continuously threatened by the harsh survival conditions inside host macrophages, due to immune and antibiotic stresses. Faithful genome maintenance and repair must be accomplished under stress for the bacillus to survive in the host, necessitating a robust DNA repair system. The importance of DNA repair systems in pathogenesis is well established. Previous examination of the M. tuberculosis genome revealed homologues of almost all the major DNA repair systems, i.e. nucleotide excision repair (NER), base excision repair (BER), homologous recombination (HR) and non-homologous end joining (NHEJ). However, recent developments in the field have pointed to the presence of novel proteins and pathways in mycobacteria. Homologues of archeal mismatch repair proteins were recently reported in mycobacteria, a pathway previously thought to be absent. RecBCD, the major nuclease-helicase enzymes involved in HR in E. coli, were implicated in the single-strand annealing (SSA) pathway. Novel roles of archeo-eukaryotic primase (AEP) polymerases, previously thought to be exclusive to NHEJ, have been reported in BER. Many new proteins with a probable role in DNA repair have also been discovered. It is now realized that the DNA repair systems in M. tuberculosis are highly evolved and have redundant backup mechanisms to mend the damage. This review is an attempt to summarize our current understanding of the DNA repair systems in M. tuberculosis.

Comparative transcriptional profiling of tildipirosin-resistant and sensitive Haemophilus parasuis

Numerous studies have been conducted to examine the molecular mechanism of Haemophilus parasuis resistance to antibiotic, but rarely to tildipirosin. In the current study, transcriptional profiling was applied to analyse the variation in gene expression of JS0135 and tildipirosin-resistant JS32. The growth curves showed that JS32 had a higher growth rate but fewer bacteria than JS0135. The cell membranes of JS32 and a resistant clinical isolate (HB32) were observed to be smoother than those of JS0135. From the comparative gene expression profile 349 up- and 113 downregulated genes were observed, covering 37 GO and 63 KEGG pathways which are involved in biological processes (11), cellular components (17), molecular function (9), cellular processes (1), environmental information processing (4), genetic information processing (9) and metabolism (49) affected in JS32. In addition, the relative overexpression of genes of the metabolism pathway (HAPS_RS09315, HAPS_RS09320), ribosomes (HAPS_RS07815) and ABC transporters (HAPS_RS10945) was detected, particularly the metabolism pathway, and verified with RT-qPCR. Collectively, the gene expression profile in connection with tildipirosin resistance factors revealed unique and highly resistant determinants of H. parasuis to macrolides that warrant further attention due to the significant threat of bacterial resistance.