CRISPR-mediated adaptive immune systems in bacteria and archaea.

Effective clearance of an infection requires that the immune system rapidly detects and neutralizes invading parasites while strictly avoiding self-antigens that would result in autoimmunity. The cellular machinery and complex signaling pathways that coordinate an effective immune response have generally been considered properties of the eukaryotic immune system. However, a surprisingly sophisticated adaptive immune system that relies on small RNAs for sequence-specific targeting of foreign nucleic acids was recently discovered in bacteria and archaea. Molecular vaccination in prokaryotes is achieved by integrating short fragments of foreign nucleic acids into a repetitive locus in the host chromosome known as a CRISPR (clustered regularly interspaced short palindromic repeat). Here we review the mechanisms of CRISPR-mediated immunity and discuss the ecological and evolutionary implications of these adaptive defense systems.

SAFB2 Enables the Processing of Suboptimal Stem-Loop Structures in Clustered Primary miRNA Transcripts.

Many microRNAs (miRNAs) are generated from primary transcripts containing multiple clustered stem-loop structures that are thought to be recognized and cleaved by the Microprocessor complex as independent units. Here, we uncover an unexpected mode of processing of the bicistronic miR-15a-16-1 cluster. We find that the primary miR-15a stem-loop is not processed on its own but that the presence of the neighboring primary miR-16-1 stem-loop on the same transcript can compensate for this deficiency in cis. Using a CRISPR/Cas9 screen, we identify SAFB2 (scaffold attachment factor B2) as an essential co-factor in this miR-16-1-assisted pri-miR-15 cleavage and describe SAFB2 as an accessory protein of the Microprocessor. Notably, SAFB2-mediated cleavage expands to other clustered pri-miRNAs, indicating a general mechanism. Together, our study reveals an unrecognized function of SAFB2 in miRNA processing and suggests a scenario in which SAFB2 enables the binding and processing of suboptimal Microprocessor substrates in clustered primary miRNA transcripts.

Epigenetic therapy induces transcription of inverted SINEs and ADAR1 dependency.

Cancer therapies that target epigenetic repressors can mediate their effects by activating retroelements within the human genome. Retroelement transcripts can form double-stranded RNA (dsRNA) that activates the MDA5 pattern recognition receptor1-6. This state of viral mimicry leads to loss of cancer cell fitness and stimulates innate and adaptive immune responses7,8. However, the clinical efficacy of epigenetic therapies has been limited. To find targets that would synergize with the viral mimicry response, we sought to identify the immunogenic retroelements that are activated by epigenetic therapies. Here we show that intronic and intergenic SINE elements, specifically inverted-repeat Alus, are the major source of drug-induced immunogenic dsRNA. These inverted-repeat Alus are frequently located downstream of 'orphan' CpG islands9. In mammals, the ADAR1 enzyme targets and destabilizes inverted-repeat Alu dsRNA10, which prevents activation of the MDA5 receptor11. We found that ADAR1 establishes a negative-feedback loop, restricting the viral mimicry response to epigenetic therapy. Depletion of ADAR1 in patient-derived cancer cells potentiates the efficacy of epigenetic therapy, restraining tumour growth and reducing cancer initiation. Therefore, epigenetic therapies trigger viral mimicry by inducing a subset of inverted-repeats Alus, leading to an ADAR1 dependency. Our findings suggest that combining epigenetic therapies with ADAR1 inhibitors represents a promising strategy for cancer treatment.

A mutually exclusive stem-loop arrangement in roX2 RNA is essential for X-chromosome regulation in Drosophila.

The X chromosome provides an ideal model system to study the contribution of RNA-protein interactions in epigenetic regulation. In male flies, roX long noncoding RNAs (lncRNAs) harbor several redundant domains to interact with the ubiquitin ligase male-specific lethal 2 (MSL2) and the RNA helicase Maleless (MLE) for X-chromosomal regulation. However, how these interactions provide the mechanics of spreading remains unknown. By using the uvCLAP (UV cross-linking and affinity purification) methodology, which provides unprecedented information about RNA secondary structures in vivo, we identified the minimal functional unit of roX2 RNA. By using wild-type and various MLE mutant derivatives, including a catalytically inactive MLE derivative, MLEGET, we show that the minimal roX RNA contains two mutually exclusive stem-loops that exist in a peculiar structural arrangement: When one stem-loop is unwound by MLE, an alternate structure can form, likely trapping MLE in this perpetually structured region. We show that this functional unit is necessary for dosage compensation, as mutations that disrupt this formation lead to male lethality. Thus, we propose that roX2 lncRNA contains an MLE-dependent affinity switch to enable reversible interactions of the MSL complex to allow dosage compensation of the X chromosome.

LIRBase: a comprehensive database of long inverted repeats in eukaryotic genomes.

Small RNAs (sRNAs) constitute a large portion of functional elements in eukaryotic genomes. Long inverted repeats (LIRs) can be transcribed into long hairpin RNAs (hpRNAs), which can further be processed into small interfering RNAs (siRNAs) with vital biological roles. In this study, we systematically identified a total of 6 619 473 LIRs in 424 eukaryotic genomes and developed LIRBase (https://venyao.xyz/lirbase/), a specialized database of LIRs across different eukaryotic genomes aiming to facilitate the annotation and identification of LIRs encoding long hpRNAs and siRNAs. LIRBase houses a comprehensive collection of LIRs identified in a wide range of eukaryotic genomes. In addition, LIRBase not only allows users to browse and search the identified LIRs in any eukaryotic genome(s) of interest available in GenBank, but also provides friendly web functionalities to facilitate users to identify LIRs in user-uploaded sequences, align sRNA sequencing data to LIRs, perform differential expression analysis of LIRs, predict mRNA targets for LIR-derived siRNAs, and visualize the secondary structure of candidate long hpRNAs encoded by LIRs. As demonstrated by two case studies, collectively, LIRBase bears the great utility for systematic investigation and characterization of LIRs and functional exploration of potential roles of LIRs and their derived siRNAs in diverse species.

Opposing roles for Egalitarian and Staufen in transport, anchoring and localization of oskar mRNA in the Drosophila oocyte.

Localization of oskar mRNA includes two distinct phases: transport from nurse cells to the oocyte, a process typically accompanied by cortical anchoring in the oocyte, followed by posterior localization within the oocyte. Signals within the oskar 3' UTR directing transport are individually weak, a feature previously hypothesized to facilitate exchange between the different localization machineries. We show that alteration of the SL2a stem-loop structure containing the oskar transport and anchoring signal (TAS) removes an inhibitory effect such that in vitro binding by the RNA transport factor, Egalitarian, is elevated as is in vivo transport from the nurse cells into the oocyte. Cortical anchoring within the oocyte is also enhanced, interfering with posterior localization. We also show that mutation of Staufen recognized structures (SRSs), predicted binding sites for Staufen, disrupts posterior localization of oskar mRNA just as in staufen mutants. Two SRSs in SL2a, one overlapping the Egalitarian binding site, are inferred to mediate Staufen-dependent inhibition of TAS anchoring activity, thereby promoting posterior localization. The other three SRSs in the oskar 3' UTR are also required for posterior localization, including two located distant from any known transport signal. Staufen, thus, plays multiple roles in localization of oskar mRNA.

Soybean DICER-LIKE2 Regulates Seed Coat Color via Production of Primary 22-Nucleotide Small Interfering RNAs from Long Inverted Repeats.

In plants, 22-nucleotide small RNAs trigger the production of secondary small interfering RNAs (siRNAs) and enhance silencing. DICER-LIKE2 (DCL2)-dependent 22-nucleotide siRNAs are rare in Arabidopsis (Arabidopsis thaliana) and are thought to function mainly during viral infection; by contrast, these siRNAs are abundant in many crops such as soybean (Glycine max) and maize (Zea mays). Here, we studied soybean 22-nucleotide siRNAs by applying CRISPR-Cas9 to simultaneously knock out the two copies of soybean DCL2, GmDCL2a and GmDCL2b, in the Tianlong1 cultivar. Small RNA sequencing revealed that most 22-nucleotide siRNAs are derived from long inverted repeats (LIRs) and disappeared in the Gmdcl2a/2b double mutant. De novo assembly of a Tianlong1 reference genome and transcriptome profiling identified an intronic LIR formed by the chalcone synthase (CHS) genes CHS1 and CHS3 This LIR is the source of primary 22-nucleotide siRNAs that target other CHS genes and trigger the production of secondary 21-nucleotide siRNAs. Disruption of this process in Gmdcl2a/2b mutants substantially increased CHS mRNA levels in the seed coat, thus changing the coat color from yellow to brown. Our results demonstrated that endogenous LIR-derived transcripts in soybean are predominantly processed by GmDCL2 into 22-nucleotide siRNAs and uncovered a role for DCL2 in regulating natural traits.

DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome.

Transposable elements are viewed as 'selfish genetic elements', yet they contribute to gene regulation and genome evolution in diverse ways. More than half of the human genome consists of transposable elements. Alu elements belong to the short interspersed nuclear element (SINE) family of repetitive elements, and with over 1 million insertions they make up more than 10% of the human genome. Despite their abundance and the potential evolutionary advantages they confer, Alu elements can be mutagenic to the host as they can act as splice acceptors, inhibit translation of mRNAs and cause genomic instability. Alu elements are the main targets of the RNA-editing enzyme ADAR and the formation of Alu exons is suppressed by the nuclear ribonucleoprotein HNRNPC, but the broad effect of massive secondary structures formed by inverted-repeat Alu elements on RNA processing in the nucleus remains unknown. Here we show that DHX9, an abundant nuclear RNA helicase, binds specifically to inverted-repeat Alu elements that are transcribed as parts of genes. Loss of DHX9 leads to an increase in the number of circular-RNA-producing genes and amount of circular RNAs, translational repression of reporters containing inverted-repeat Alu elements, and transcriptional rewiring (the creation of mostly nonsensical novel connections between exons) of susceptible loci. Biochemical purifications of DHX9 identify the interferon-inducible isoform of ADAR (p150), but not the constitutively expressed ADAR isoform (p110), as an RNA-independent interaction partner. Co-depletion of ADAR and DHX9 augments the double-stranded RNA accumulation defects, leading to increased circular RNA production, revealing a functional link between these two enzymes. Our work uncovers an evolutionarily conserved function of DHX9. We propose that it acts as a nuclear RNA resolvase that neutralizes the immediate threat posed by transposon insertions and allows these elements to evolve as tools for the post-transcriptional regulation of gene expression.

YB-1 unwinds mRNA secondary structures in vitro and negatively regulates stress granule assembly in HeLa cells.

In the absence of the scanning ribosomes that unwind mRNA coding sequences and 5'UTRs, mRNAs are likely to form secondary structures and intermolecular bridges. Intermolecular base pairing of non polysomal mRNAs is involved in stress granule (SG) assembly when the pool of mRNAs freed from ribosomes increases during cellular stress. Here, we unravel the structural mechanisms by which a major partner of dormant mRNAs, YB-1 (YBX1), unwinds mRNA secondary structures without ATP consumption by using its conserved cold-shock domain to destabilize RNA stem/loops and its unstructured C-terminal domain to secure RNA unwinding. At endogenous levels, YB-1 facilitates SG disassembly during arsenite stress recovery. In addition, overexpression of wild-type YB-1 and to a lesser extent unwinding-defective mutants inhibit SG assembly in HeLa cells. Through its mRNA-unwinding activity, YB-1 may thus inhibit SG assembly in cancer cells and package dormant mRNA in an unfolded state, thus preparing mRNAs for translation initiation.

Integrity of a heterochromatic domain ensured by its boundary elements.

In fission yeast, the inverted repeats IR-L and IR-R function as boundary elements at the edges of a 20-kb silent heterochromatic domain where nucleosomes are methylated at histone H3K9. Each repeat contains a series of B-box motifs physically associated with the architectural TFIIIC complex and with other factors including the replication regulator Sap1 and the Rix1 complex (RIXC). We demonstrate here the activity of these repeats in heterochromatin formation and maintenance. Deletion of the entire IR-R repeat or, to a lesser degree, deletion of just the B boxes impaired the de novo establishment of the heterochromatic domain. Nucleation proceeded normally at the RNA interference (RNAi)-dependent element cenH but subsequent propagation to the rest of the region occurred at reduced rates in the mutants. Once established, heterochromatin was unstable in the mutants. These defects resulted in bistable populations of cells occupying alternate "on" and "off" epigenetic states. Deleting IR-L in combination with IR-R synergistically tipped the balance toward the derepressed state, revealing a concerted action of the two boundaries at a distance. The nuclear rim protein Amo1 has been proposed to tether the mating-type region and its boundaries to the nuclear envelope, where Amo1 mutants displayed milder phenotypes than boundary mutants. Thus, the boundaries might facilitate heterochromatin propagation and maintenance in ways other than just through Amo1, perhaps by constraining a looped domain through pairing.

Protein arginine methyltransferase CARM1 attenuates the paraspeckle-mediated nuclear retention of mRNAs containing IRAlus.

In many cells, mRNAs containing inverted repeated Alu elements (IRAlus) in their 3' untranslated regions (UTRs) are inefficiently exported to the cytoplasm. Such nuclear retention correlates with paraspeckle-associated protein complexes containing p54(nrb). However, nuclear retention of mRNAs containing IRAlus is variable, and how regulation of retention and export is achieved is poorly understood. Here we show one mechanism of such regulation via the arginine methyltransferase CARM1 (coactivator-associated arginine methyltransferase 1). We demonstrate that disruption of CARM1 enhances the nuclear retention of mRNAs containing IRAlus. CARM1 regulates this nuclear retention pathway at two levels: CARM1 methylates the coiled-coil domain of p54(nrb), resulting in reduced binding of p54(nrb) to mRNAs containing IRAlus, and also acts as a transcription regulator to suppress NEAT1 transcription, leading to reduced paraspeckle formation. These actions of CARM1 work together synergistically to regulate the export of transcripts containing IRAlus from paraspeckles under certain cellular stresses, such as poly(I:C) treatment. This work demonstrates how a post-translational modification of an RNA-binding protein affects protein-RNA interaction and also uncovers a mechanism of transcriptional regulation of the long noncoding RNA NEAT1.

Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability.

Transposable elements (TEs) are a rich source of genetic variability. Among TEs, miniature inverted-repeat TEs (MITEs) are of particular interest as they are present in high copy numbers in plant genomes and are closely associated with genes. MITEs are deletion derivatives of class II transposons, and can be mobilized by the transposases encoded by the latter through a typical cut-and-paste mechanism. However, MITEs are typically present at much higher copy numbers than class II transposons. We present here an analysis of 103 109 transposon insertion polymorphisms (TIPs) in 738 Oryza sativa genomes representing the main rice population groups. We show that an important fraction of MITE insertions has been fixed in rice concomitantly with its domestication. However, another fraction of MITE insertions is present at low frequencies. We performed MITE TIP-genome-wide association studies (TIP-GWAS) to study the impact of these elements on agronomically important traits and found that these elements uncover more trait associations than single nucleotide polymorphisms (SNPs) on important phenotypes such as grain width. Finally, using SNP-GWAS and TIP-GWAS we provide evidence of the replicative amplification of MITEs.

Hairpin Transfer-Independent Parvovirus DNA Replication Produces Infectious Virus.

Parvoviruses package a linear single-stranded DNA genome with hairpin structures at both ends. It has been thought that terminal hairpin sequences are indispensable for viral DNA replication. Here, we provide evidence that the hairpin-deleted duplex genomes of human bocavirus 1 (HBoV1) replicate in human embryonic kidney 293 (HEK293) cells. We propose an alternative model for HBoV1 DNA replication in which the leading strand can initiate strand displacement without hairpin transfer. The transfection of the HBoV1 duplex genomes that retain a minimal replication origin at the right end (OriR) but with extensive deletions in the right-end hairpin (REH) generated viruses in HEK293 cells at a level 10 to 20 times lower than that of the wild-type (WT) duplex genome. Importantly, these viruses that have a genome with various deletions after the OriR but not the one retaining only the OriR replicated in polarized human airway epithelia. We discovered that the 18-nucleotide (nt) sequence (nt 5403 to 5420) beyond the OriR was sufficient to confer virus replication in polarized human airway epithelia, although its progeny virus production was ∼5 times lower than that of the WT virus. Thus, our study demonstrates that hairpin transfer-independent productive parvovirus DNA replication can occur. IMPORTANCE Hairpin transfer-independent parvovirus replication was modeled with human bocavirus 1 (HBoV1) duplex genomes whose 5' hairpin structure was ablated by various deletions. In HEK293 cells, these duplex viral genomes with ablated 5' hairpin sequence replicated efficiently and generated viruses that productively infected polarized human airway epithelium. Thus, for the first time, we reveal a previously unknown phenomenon that productive parvovirus DNA replication does not depend on the hairpin sequence at REH to initiate rolling-hairpin DNA replication. Notably, the intermediates of viral DNA replication, as revealed by two-dimensional electrophoresis, from transfections of hairpin sequence-deleted duplex genome and full-length genome in HEK293 cells as well as from virus infection of polarized human airway epithelia are similar. Thus, the establishment of the hairpin transfer-independent parvoviral DNA replication deepens our understanding of viral DNA replication and may have implications in the development of parvovirus-based viral vectors with alternative properties.

CtIP maintains stability at common fragile sites and inverted repeats by end resection-independent endonuclease activity.

Chromosomal rearrangements often occur at genomic loci with DNA secondary structures, such as common fragile sites (CFSs) and palindromic repeats. We developed assays in mammalian cells that revealed CFS-derived AT-rich sequences and inverted Alu repeats (Alu-IRs) are mitotic recombination hotspots, requiring the repair functions of carboxy-terminal binding protein (CtBP)-interacting protein (CtIP) and the Mre11/Rad50/Nbs1 complex (MRN). We also identified an endonuclease activity of CtIP that is dispensable for end resection and homologous recombination (HR) at I-SceI-generated "clean" double-strand breaks (DSBs) but is required for repair of DSBs occurring at CFS-derived AT-rich sequences. In addition, CtIP nuclease-defective mutants are impaired in Alu-IRs-induced mitotic recombination. These studies suggest that an end resection-independent CtIP function is important for processing DSB ends with secondary structures to promote HR. Furthermore, our studies uncover an important role of MRN, CtIP, and their associated nuclease activities in protecting CFSs in mammalian cells.

NF90 modulates processing of a subset of human pri-miRNAs.

MicroRNAs (miRNAs) are predicted to regulate the expression of >60% of mammalian genes and play fundamental roles in most biological processes. Deregulation of miRNA expression is a hallmark of most cancers and further investigation of mechanisms controlling miRNA biogenesis is needed. The double stranded RNA-binding protein, NF90 has been shown to act as a competitor of Microprocessor for a limited number of primary miRNAs (pri-miRNAs). Here, we show that NF90 has a more widespread effect on pri-miRNA biogenesis than previously thought. Genome-wide approaches revealed that NF90 is associated with the stem region of 38 pri-miRNAs, in a manner that is largely exclusive of Microprocessor. Following loss of NF90, 22 NF90-bound pri-miRNAs showed increased abundance of mature miRNA products. NF90-targeted pri-miRNAs are highly stable, having a lower free energy and fewer mismatches compared to all pri-miRNAs. Mutations leading to less stable structures reduced NF90 binding while increasing pri-miRNA stability led to acquisition of NF90 association, as determined by RNA electrophoretic mobility shift assay (EMSA). NF90-bound and downregulated pri-miRNAs are embedded in introns of host genes and expression of several host genes is concomitantly reduced. These data suggest that NF90 controls the processing of a subset of highly stable, intronic miRNAs.

Do the P1 and P2 hairpins of the Guanidine-II riboswitch interact?

Riboswitches regulate genes by adopting different structures in responds to metabolite binding. The guanidine-II riboswitch is the smallest representative of the ykkC class with the mechanism of its function being centred on the idea that its two stem loops P1 and P2 form a kissing hairpin interaction upon binding of guanidinium (Gdm+). This mechanism is based on in-line probing experiments with the full-length riboswitch and crystal structures of the truncated stem loops P1 and P2. However, the crystal structures reveal only the formation of the homodimers P1 | P1 and P2 | P2 but not of the proposed heterodimer P1 | P2. Here, site-directed spin labeling (SDSL) in combination with Pulsed Electron-Electron Double Resonance (PELDOR or DEER) is used to study their structures in solution and how they change upon binding of Gdm+. It is found that both hairpins adopt different structures in solution and that binding of Gdm+ does indeed lead to the formation of the heterodimer but alongside the homodimers in a statistical 1:2:1 fashion. These results do thus support the proposed switching mechanism.

Stem-loop 4 of U1 snRNA is essential for splicing and interacts with the U2 snRNP-specific SF3A1 protein during spliceosome assembly.

The pairing of 5' and 3' splice sites across an intron is a critical step in spliceosome formation and its regulation. Interactions that bring the two splice sites together during spliceosome assembly must occur with a high degree of specificity and fidelity to allow expression of functional mRNAs and make particular alternative splicing choices. Here, we report a new interaction between stem-loop 4 (SL4) of the U1 snRNA, which recognizes the 5' splice site, and a component of the U2 small nuclear ribonucleoprotein particle (snRNP) complex, which assembles across the intron at the 3' splice site. Using a U1 snRNP complementation assay, we found that SL4 is essential for splicing in vivo. The addition of free U1-SL4 to a splicing reaction in vitro inhibits splicing and blocks complex assembly prior to formation of the prespliceosomal A complex, indicating a requirement for a SL4 contact in spliceosome assembly. To characterize the interactions of this RNA structure, we used a combination of stable isotope labeling by amino acids in cell culture (SILAC), biotin/Neutravidin affinity pull-down, and mass spectrometry. We show that U1-SL4 interacts with the SF3A1 protein of the U2 snRNP. We found that this interaction between the U1 snRNA and SF3A1 occurs within prespliceosomal complexes assembled on the pre-mRNA. Thus, SL4 of the U1 snRNA is important for splicing, and its interaction with SF3A1 mediates contact between the 5' and 3' splice site complexes within the assembling spliceosome.

Mutation of the 5'-untranslated region stem-loop mRNA structure reduces type I collagen deposition and arterial stiffness in male obese mice.

Arterial stiffening, a characteristic feature of obesity and type 2 diabetes, contributes to the development and progression of cardiovascular diseases (CVD). Currently, no effective prophylaxis or therapeutics is available to prevent or treat arterial stiffening. A better understanding of the molecular mechanisms underlying arterial stiffening is vital to identify newer targets and strategies to reduce CVD burden. A major contributor to arterial stiffening is increased collagen deposition. In the 5'-untranslated regions of mRNAs encoding for type I collagen, an evolutionally conserved stem-loop (SL) structure plays an essential role in its stability and post-transcriptional regulation. Here, we show that feeding a high-fat/high-sucrose (HFHS) diet for 28 wk increases adiposity, insulin resistance, and blood pressure in male wild-type littermates. Moreover, arterial stiffness, assessed in vivo via aortic pulse wave velocity, and ex vivo using atomic force microscopy in aortic explants or pressure myography in isolated femoral and mesenteric arteries, was also increased in those mice. Notably, all these indices of arterial stiffness, along with collagen type I levels in the vasculature, were reduced in HFHS-fed mice harboring a mutation in the 5'SL structure, relative to wild-type littermates. This protective vascular phenotype in 5'SL-mutant mice did not associate with a reduction in insulin resistance or blood pressure. These findings implicate the 5'SL structure as a putative therapeutic target to prevent or reverse arterial stiffening and CVD associated with obesity and type 2 diabetes.NEW & NOTEWORTHY In the 5'-untranslated (UTR) regions of mRNAs encoding for type I collagen, an evolutionally conserved SL structure plays an essential role in its stability and posttranscriptional regulation. We demonstrate that a mutation of the SL mRNA structure in the 5'-UTR decreases collagen type I deposition and arterial stiffness in obese mice. Targeting this evolutionarily conserved SL structure may hold promise in the management of arterial stiffening and CVD associated with obesity and type 2 diabetes.

Processing-independent CRISPR RNAs limit natural transformation in Neisseria meningitidis.

CRISPR interference confers adaptive, sequence-based immunity against viruses and plasmids and is specified by CRISPR RNAs (crRNAs) that are transcribed and processed from spacer-repeat units. Pre-crRNA processing is essential for CRISPR interference in all systems studied thus far. Here, our studies of crRNA biogenesis and CRISPR interference in naturally competent Neisseria spp. reveal a unique crRNA maturation pathway in which crRNAs are transcribed from promoters that are embedded within each repeat, yielding crRNA 5' ends formed by transcription and not by processing. Although crRNA 3' end formation involves RNase III and trans-encoded tracrRNA, as in other type II CRISPR systems, this processing is dispensable for interference. The meningococcal pathway is the most streamlined CRISPR/Cas system characterized to date. Endogenous CRISPR spacers limit natural transformation, which is the primary source of genetic variation that contributes to immune evasion, antibiotic resistance, and virulence in the human pathogen N. meningitidis.

Coding palindromes in mitochondrial genes of Nematomorpha.

Inverted repeats are common DNA elements, but they rarely overlap with protein-coding sequences due to the ensuing conflict with the structure and function of the encoded protein. We discovered numerous perfect inverted repeats of considerable length (up to 284 bp) embedded within the protein-coding genes in mitochondrial genomes of four Nematomorpha species. Strikingly, both arms of the inverted repeats encode conserved regions of the amino acid sequence. We confirmed enzymatic activity of the respiratory complex I encoded by inverted repeat-containing genes. The nucleotide composition of inverted repeats suggests strong selection at the amino acid level in these regions. We conclude that the inverted repeat-containing genes are transcribed and translated into functional proteins. The survey of available mitochondrial genomes reveals that several other organisms possess similar albeit shorter embedded repeats. Mitochondrial genomes of Nematomorpha demonstrate an extraordinary evolutionary compromise where protein function and stringent secondary structure elements within the coding regions are preserved simultaneously.