Cell free extrachromosomal circular DNA is common in human urine

Cell free extrachromosomal circular DNA (eccDNA) is evolving as a potential biomarker in liquid biopsies for disease diagnosis. In this study, an optimized next generation sequencing-based Circle-Seq method was developed to investigate urinary cell free eccDNA (ucf-eccDNA) from 28 adult healthy volunteers (mean age = 28, 19 males/ 9 females). The genomic distributions and sequence compositions of ucf-eccDNAs were comprehensively characterized. Approximately 1.2 million unique ucf-eccDNAs are identified, covering 14.9% of the human genome. Comprehensive characterization of ucf-eccDNAs show that ucf-eccDNAs contain higher GC content than flanking genomic regions. Most eccDNAs are less than 1000 bp and present four pronounced peaks at 203, 361, 550 and 728 bp, indicating the association between eccDNAs and the numbers of intact nucleosomes. Analysis of genomic distribution of ucf-eccDNAs show that eccDNAs are found in all chromosomes but enriched in chromosomes i.e. chr.17, 19 and 20 with high density of protein-codding genes, CpG islands, SINE and simple repeat elements. Lastly, analysis of sequence motif signatures at eccDNA junction sites reveal that direct repeats (DRs) are commonly found, indicating a potential role of DRs in eccDNA biogenesis. This work underscores the deep sequencing analysis of ucf-eccDNAs and provides a valuable reference resource for exploring potential applications of ucf-eccDNA as diagnostic biomarkers of urogenital disorders in the future.Significance StatementExtrachromosomal circular DNA (eccDNA) is an important genetic element and a biomarker for disease diagnosis and treatment. In this study, we conduct a comprehensive characterization of urinary cell free eccDNA (ucf-eccDNA) in 28 heathy subjects. Over one million ucf-eccDNAs are identified. Ucf-eccDNAs are characterized as high GC content. The size of most ucf-eccDNAs is less than 1000 bp and enriched in four peaks resembling the size of single, double, triple, and quadruple nucleosomes. The genomic distribution of ucf-eccDNAs is enriched in generic regions, protein-coding genes, Alu, CpG islands, SINE and simple repeats. Sequence motif analysis of ucf-eccDNA junctions identified simple direct repeats (DRs) commonly presented in most eccDNAs, suggesting potential roles of DRs in eccDNA biogenesis.

The aadE*-sat4-aphA-3 Gene Cluster of Mycoplasma bovirhinis HAZ141_2 Undergoes Genomic Rearrangements Influencing the Primary Promoter Sequence

The 54 kb GC-rich prophage region of Mycoplasma bovirhinis HAZ141_2 contains three structural ‘compartments’, one of which is a highly transmittable cluster of three genes, aadE-like (aadE*), sat4, and aphA-3. In this study, we characterized recombination events and their consequences occurred within the aadE*-sat4-aphA-3 containing region. Analysis of this region revealed direct repeats (DRs) of 155 and invert repeats (IRs) of 197 base pairs (bps) each, flanking and overlapping with the primary promoter P* located upstream of the aadE*. Two recombination events, including inversions via both 197 and 155-bp IRs (the latter become inverted after the initial 197-bp IRs associated inversion) and the excision of the aadE*-sat4-aphA-3 cluster, were confirmed. Inversion via 155-IRs results in changes within the P* promoter region. Using Escherichia coli JM109 carrying plasmids containing derivatives of the aadE*-sat4-aphA-3 cluster, we validated the expression of those genes from different promoters. Our results showed no difference in the minimal inhibitory concentrations (MICs) to kanamycin and neomycin and only 2-fold decrease in MIC (from 512 to 256 μg/mL) to nourseothricin between the wild type and a P* derivative promoter. However, the MICs to kanamycin and neomycin were at least 4-fold lower in the construct where aphA-3 expressed under its P2 promoter (128 µg/mL) in comparison to the construct where aphA-3 expressed under P1” promoter located within the sat4 gene (512–1024 µg/mL). PCR confirmed the excision of the aadE*-sat4-aphA-3 cluster via 197- and 155-bp DRs, but no selection of antibiotic-sensitive M. bovirhinis were obtained after 100 passages in kanamycin-free medium.

APOBEC mutagenesis is low in most types of non-B DNA structures, unlike other types of cancer mutagenesis

While somatic mutations are known to be enriched in genome regions with non-canonical DNA secondary structure, the impact of particular mutagens still needs to be elucidated. Here, we demonstrate that in human cancers, the APOBEC mutagenesis is not enriched in direct repeats, mirror repeats, short tandem repeats, and G-quadruplexes, and even decreased below its level in B-DNA for cancer samples with very high APOBEC activity. In contrast, we observe that the APOBEC-induced mutational density is positively associated with APOBEC activity in inverted repeats (cruciform structures), where the impact of cytosine at the 3'-end of the hairpin loop is substantial. Surprisingly, the APOBEC-signature mutation density per TC motif in the single-stranded DNA of a G-quadruplex (G4) is lower than in the four-stranded part of G4 and in B-DNA. The APOBEC mutagenesis, as well as the UV-mutagenesis in melanoma samples are absent in Z-DNA regions, due to depletion of their mutational signature motifs.

Identification of Three Clf-Sdr Subfamily Proteins in Staphylococcus warneri, and Comparative Genomics Analysis of a Locus Encoding CWA Proteins in Staphylococcus Species

Coagulase-negative Staphylococcus warneri is an opportunistic pathogen that is capable of causing several infections, especially in patients with indwelling medical devices. Here, we determined the complete genome sequence of a clinical S. warneri strain isolated from the blood culture of a 1-year-old nursling patient with acute upper respiratory infection. Genome-wide phylogenetic analysis confirmed the phylogenetic relationships between S. warneri and other Staphylococcus species. Using comparative genomics, we identified three cell wall-anchored (CWA) proteins at the same locus (sdr), named SdrJ, SdrK, and SdrL, on the chromosome sequences of different S. warneri strains. Structural predictions showed that SdrJ/K/L have structural features characteristic of Sdr proteins but exceptionally contained an unusual N-terminal repeat region. However, the C-terminal repetitive (R) region of SdrJ contains a significantly larger proportion of alanine (142/338, 42.01%) than the previously reported SdrI (37.00%). Investigation of the genetic organization revealed that the sdrJ/K/L genes were always followed by one or two glycosyltransferase genes, gtfA and gtfB and were present in an ∼56 kb region bordered by a pair of 8 bp identical direct repeats, named Sw-Sdr. This region was further found to be located on a 160-kb region subtended by a pair of 160-bp direct repeats along with other virulence genes and resistance genes. Sw-Sdr contained a putative integrase that was probably a remnant of a functional integrase. Evidence suggests that Sw-Sdr is improbably an efficient pathogenicity island. A large-scale investigation of Staphylococcus genomes showed that sdr loci were a potential hotspot of insertion sequences (ISs), which could lead to intraspecific diversity at these loci. Our work expanded the repository of Staphylococcus Sdr proteins, and for the first time, we established the connection between sdr loci and phylogenetic relationships and compared the sdr loci in different Staphylococcus species, which provided large insights into the genetic environment of CWA genes in Staphylococcus.

CRISPRbuilder-TB: “CRISPR-builder for tuberculosis”. Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA

Mycobacterium tuberculosis complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains. In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in M. tuberculosis. We first show that the process we set up, that we coin as “CRISPRbuilder-TB” (https://github.com/cguyeux/CRISPRbuilder-TB), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental in vitro approach and the basic in silico spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS6110 insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of M. tuberculosis CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.

Four large indels detected by cpTILLING in barley chloroplast mutator seedlings

In a previous work, a polymorphism detection strategy based on mismatch digestion was applied to the chloroplast genome of barley seedlings that carried the chloroplast mutator (cpm) genotype through many generations. Sixty-two different one- or two-nucleotide-polymorphisms were detected along with four large indels: an insertion of 15 bp in the intergenic region between tRNAHis and rps19 genes, a deletion of 620 bp in the psbA gene, a deletion of 79 bp in the intergenic region between rpl33 and rps18 genes and a deletion of 45 bp in the rps3 gene. In the present investigation, we analyzed direct repeats located at the borders of those four large indels. Furthermore, we investigated the consequences of protein expression of large indels located in coding regions. The deletion of 620 bp in the psbA gene was lethal at the second leaf stage when homoplastomic. The deletion of 45 bp in the rps3 gene, which eliminates 15 amino acids, did not affect the viability of the seedlings in homoplastomy. Interestingly, the deleted segment is also lacking in the wild type version of the rps3 gene of maize and sorghum. The presence of direct repeats at the borders of the four large indels suggests that they could have originated by illegitimate recombination. This would be in agreement with a previous hypothesis that the Cpm gene product would correspond to a mismatch repair (MMR) protein devoted to maintain plastome stability by playing fundamental roles in mismatch repair during replication and avoiding illegitimate recombination.

Mobile Plasmid Mediated Transition From Colistin-Sensitive to Resistant Phenotype in Klebsiella pneumoniae

Multidrug-resistant bacteria, including carbapenem-resistant Klebsiella pneumoniae (CRKP), are becoming an increasing health crisis worldwide. For CRKP, colistin is regarded as “the last treatment option.” In this study, we isolated a clinical CRKP strain named as K. pneumoniae R10-341. Phenotyping analysis showed that this strain could transit from a colistin-sensitive to a resistant phenotype by inserting an IS4 family ISKpn72 element into the colistin-resistance associated mgrB gene. To investigate the mechanism of this transition, we performed genome sequencing analysis of the colistin-sensitive parental strain and found that 12 copies of ISKpn72 containing direct repeats (DR) are located on the chromosome and 1 copy without DR is located on a multidrug-resistant plasmid pR10-341_2. Both types of ISKpn72 could be inserted into the mgrB gene to cause colistin-resistance, though the plasmid-derived ISKpn72 without DR was in higher efficiency. Importantly, we demonstrated that colistin-sensitive K. pneumoniae strain transferred with the ISKpn72 element also obtained the ability to switch from colistin-sensitive to colistin-resistant phenotype. Furthermore, we confirmed that the ISKpn72-containing pR10-341_2 plasmid was able to conjugate, suggesting that the ability of causing colistin-resistant transition is transferable through common conjugation. Our results point to new challenges for both colistin-resistance detection and CRKP treatment.

Deletion of the GAPDH gene contributes to genome stability in Saccharomyces cerevisiae

Cellular metabolism is directly or indirectly associated with various cellular processes by producing a variety of metabolites. Metabolic alterations may cause adverse effects on cell viability. However, some alterations potentiate the rescue of the malfunction of the cell system. Here, we found that the alteration of glucose metabolism suppressed genome instability caused by the impairment of chromatin structure. Deletion of the TDH2 gene, which encodes glyceraldehyde 3-phospho dehydrogenase and is essential for glycolysis/gluconeogenesis, partially suppressed DNA damage sensitivity due to chromatin structure, which was persistently acetylated histone H3 on lysine 56 in cells with deletions of both HST3 and HST4, encoding NAD+-dependent deacetylases. tdh2 deletion also restored the short replicative lifespan of cells with deletion of sir2, another NAD+-dependent deacetylase, by suppressing intrachromosomal recombination in rDNA repeats increased by the unacetylated histone H4 on lysine 16. tdh2 deletion also suppressed recombination between direct repeats in hst3∆ hst4∆ cells by suppressing the replication fork instability that leads to both DNA deletions among repeats. We focused on quinolinic acid (QUIN), a metabolic intermediate in the de novo nicotinamide adenine dinucleotide (NAD+) synthesis pathway, which accumulated in the tdh2 deletion cells and was a candidate metabolite to suppress DNA replication fork instability. Deletion of QPT1, quinolinate phosphoribosyl transferase, elevated intracellular QUIN levels and partially suppressed the DNA damage sensitivity of hst3∆ hst4∆ cells as well as tdh2∆ cells. qpt1 deletion restored the short replicative lifespan of sir2∆ cells by suppressing intrachromosomal recombination among rDNA repeats. In addition, qpt1 deletion could suppress replication fork slippage between direct repeats. These findings suggest a connection between glucose metabolism and genomic stability.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

Background Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes. Results We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats. Conclusion This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

A Possible Role of Insertion Sequence IS1216V in Dissemination of Multidrug-Resistant Elements MESPM1 and MES6272-2 between Enterococcus and ST59 Staphylococcus aureus

Sequence type 59 (ST59) is the dominant type of community-associated methicillin-resistant Staphylococcus aureus (MRSA) in Taiwan. Previously, we reported that ST59 MRSA harbors enterococcal IS1216V-mediated multidrug-resistant composite transposons MESPM1 or MES6272-2. The MES were found to have a mosaic structure, largely originating in enterococci and partly native to S. aureus. The current study aimed to track the origin of the MES and how they disseminated from enterococci to ST59 S. aureus. A total of 270 enterococcal isolates were analyzed, showing that two ST64 Enterococcus faecalis isolated in 1992 and 11 clonal complex 17 Enterococcus faecium harbored MESPM1-like and MES6272-2-like structures, respectively. Sequence analysis revealed that ST64 E. faecalis strain N48 acquired the MESPM1-like structure on the plasmid pEflis48. The pEflis48 harbored the enterococci-originated region (erythromycin, kanamycin, and streptomycin resistances) and the S.aureus-originated region (chloramphenicol resistance) of MESPM1 but was separated by the replication region of the plasmid. Homologous recombination between the two direct repeats of IS1216V resulted in excision of the replication region of the plasmid to regenerate MESPM1. The p4780-1 and pV19 of E. faecium carried MES6272-2-like structures with IS1216V, albeit with multiple insertions by other insertion sequences. The findings show that IS1216V plays important roles in bidirectional gene transfer of multidrug resistance between enterococci and S. aureus.