Circular DNA in the human germline and its association with recombination.

Extrachromosomal circular DNA (eccDNA) is common in somatic tissue, but its existence and effects in the human germline are unexplored. We used microscopy, long-read DNA sequencing, and new analytic methods to document thousands of eccDNAs from human sperm. EccDNAs derived from all genomic regions and mostly contained a single DNA fragment, although some consisted of multiple fragments. The generation of eccDNA inversely correlates with the meiotic recombination rate, and chromosomes with high coding-gene density and Alu element abundance form the least eccDNA. Analysis of insertions in human genomes further indicates that eccDNA can persist in the human germline when the circular molecules reinsert themselves into the chromosomes. Our results suggest that eccDNA has transient and permanent effects on the germline. They explain how differences in the physical and genetic map might arise and offer an explanation of how Alu elements coevolved with genes to protect genome integrity against deleterious mutations producing eccDNA.

Isolation and Genomic Characterization of Schaalia turicensis from a Patient with Gonococcal Urethritis, Thailand.

Schaalia turicensis is facultative anaerobic Gram-positive bacillus that commonly inhabits the oropharynx, gastrointestinal, and genitourinary tract of healthy individuals. This organism has been co-isolated with Neisseria gonorrhoeae from 15-year-old Thai male patient with gonococcal urethritis in Bangkok, Thailand. In this study, we characterized the class 1 integron in S. turicensis isolate using whole-genome sequencing and bioinformatics analysis. Sequencing analysis confirmed the presence of an imperfect class 1 integron located on chromosome and a novel 24.5-kb-long composite transposon, named Tn7083. The transposon Tn7083 carried genes encoding chloramphenicol resistance (cmx), sulfonamide resistance (sul1), and aminoglycoside resistance [aph(6)-Id (strB), aph(3'')-Ib (strA), aph(3')-Ia].

Genomic analysis of Enterococcus faecium strain RAOG174 associated with acute chorioamnionitis carried antibiotic resistance gene: is it time for precise microbiological identification for appropriate antibiotic use?

BACKGROUND: Preterm labor syndrome is associated with high perinatal morbidity and mortality, and intra-amniotic infection is a cause of preterm labor. The standard identification of causative microorganisms is based on the use of biochemical phenotypes, together with broth dilution-based antibiotic susceptibility from organisms grown in culture. However, such methods could not provide an accurate epidemiological aspect and a genetic basis of antimicrobial resistance leading to an inappropriate antibiotic administration. Hybrid genome assembly is a combination of short- and long-read sequencing, which provides better genomic resolution and completeness for genotypic identification and characterization. Herein, we performed a hybrid whole genome assembly sequencing of a pathogen associated with acute histologic chorioamnionitis in women presenting with PPROM. RESULTS: We identified Enterococcus faecium, namely E. faecium strain RAOG174, with several antibiotic resistance genes, including vancomycin and aminoglycoside. Virulence-associated genes and potential bacteriophage were also identified in this genome. CONCLUSION: We report herein the first study demonstrating the use of hybrid genome assembly and genomic analysis to identify E. faecium ST17 as a pathogen associated with acute histologic chorioamnionitis. The analysis provided several antibiotic resistance-associated genes/mutations and mobile genetic elements. The occurrence of E. faecium ST17 raised the awareness of the colonization of clinically relevant E. faecium and the carrying of antibiotic resistance. This finding has brought the advantages of genomic approach in the identification of the bacterial species and antibiotic resistance gene for E. faecium for appropriate antibiotic use to improve maternal and neonatal care.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae.

BACKGROUND: Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed. RESULTS: The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis. CONCLUSIONS: In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

The occurrence and genomic characteristics of the blaIMI-1 carbapenemase-producing Enterobacter cloacae complex retrieved from natural water sources in central Thailand.

AIM: Carbapenem resistance among Enterobacteriaceae is a serious threat to humans worldwide. This study aims to evaluate the phenotypic and genotypic characterization of carbapenemase-producing Enterobacter cloacae complex (ECC) retrieved from water sources in the central part of Thailand. METHODS AND RESULTS: Samples were collected from water bodies surrounding farms and communities in central Thailand. The species were identified by using MALDI-TOF MS. The minimum inhibitory concentration (MIC) and antibiotic susceptibility were determined. The carbapenemase-producing genes were detected by PCR and whole genome sequencing (WGS). ECC with chromosome-encoded blaIMI-1 carbapenemase were detected. These isolates were resistant to last-resort antibiotics such as carbapenems and colistin as well as penicillin. In addition, all blaIMI-1 genes isolated from this study were found to be associated with chromosomally integrated Xer-dependent integrative mobile elements (IMEXs). CONCLUSION: These findings highlight the diversity and dissemination of carbapenemases-producing Enterobacterales in environmental sources. With the increasing detection of carbapenemase genes worldwide, we should be aware of the blaIMI-producing E. cloacae complex with a high resistance profile and the ability to mobilize within the environment.

Decoding the epitranscriptional landscape from native RNA sequences.

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

Rapid diagnosis of intra-amniotic infection using nanopore-based sequencing.

OBJECTIVES: Early diagnosis and treatment of intra-amniotic infection is crucial. Rapid pathogen identification allows for a definite diagnosis and enables proper management. We determined whether the 16S amplicon sequencing performed by a nanopore sequencing technique make possible rapid bacterial identification at the species level in intra-amniotic infection. METHODS: Five cases of confirmed intra-amniotic infection, determined by either cultivation or 16S rDNA polymerase chain reaction (PCR) Sanger sequencing, and 10 cases of women who underwent mid-trimester genetic amniocentesis were included. DNA was extracted from amniotic fluid and PCR was performed on the full-length 16S rDNA. Nanopore sequencing was performed. The results derived from nanopore sequencing were compared with those derived from cultivation and Sanger sequencing methods. RESULTS: Bacteria were successfully detected from amniotic fluid using nanopore sequencing in all cases of intra-amniotic infection. Nanopore sequencing identified additional bacterial species and polymicrobial infections. All patients who underwent a mid-trimester amniocentesis had negative cultures, negative 16S PCR Sanger sequencing and nanopore sequencing. Identification of the microorganisms using nanopore sequencing technique at the bacterial species level was achieved within 5-9 h from DNA extraction. CONCLUSIONS: This is the first study demonstrating that the nanopore sequencing technique is capable of rapid diagnosis of intra-amniotic infection using fresh amniotic fluid samples.

Oligosarcomas, IDH-mutant are distinct and aggressive.

Oligodendrogliomas are defined at the molecular level by the presence of an IDH mutation and codeletion of chromosomal arms 1p and 19q. In the past, case reports and small studies described gliomas with sarcomatous features arising from oligodendrogliomas, so called oligosarcomas. Here, we report a series of 24 IDH-mutant oligosarcomas from 23 patients forming a distinct methylation class. The tumors were recurrences from prior oligodendrogliomas or developed de novo. Precursor tumors of 12 oligosarcomas were histologically and molecularly indistinguishable from conventional oligodendrogliomas. Oligosarcoma tumor cells were embedded in a dense network of reticulin fibers, frequently showing p53 accumulation, positivity for SMA and CALD1, loss of OLIG2 and gain of H3K27 trimethylation (H3K27me3) as compared to primary lesions. In 5 oligosarcomas no 1p/19q codeletion was detectable, although it was present in the primary lesions. Copy number neutral LOH was determined as underlying mechanism. Oligosarcomas harbored an increased chromosomal copy number variation load with frequent CDKN2A/B deletions. Proteomic profiling demonstrated oligosarcomas to be highly distinct from conventional CNS WHO grade 3 oligodendrogliomas with consistent evidence for a smooth muscle differentiation. Expression of several tumor suppressors was reduced with NF1 being lost frequently. In contrast, oncogenic YAP1 was aberrantly overexpressed in oligosarcomas. Panel sequencing revealed mutations in NF1 and TP53 along with IDH1/2 and TERT promoter mutations. Survival of patients was significantly poorer for oligosarcomas as first recurrence than for grade 3 oligodendrogliomas as first recurrence. These results establish oligosarcomas as a distinct group of IDH-mutant gliomas differing from conventional oligodendrogliomas on the histologic, epigenetic, proteomic, molecular and clinical level. The diagnosis can be based on the combined presence of (a) sarcomatous histology, (b) IDH-mutation and (c) TERT promoter mutation and/or 1p/19q codeletion, or, in unresolved cases, on its characteristic DNA methylation profile.

Detection and Discrimination of DNA Adducts Differing in Size, Regiochemistry, and Functional Group by Nanopore Sequencing.

Chemically induced DNA adducts can lead to mutations and cancer. Unfortunately, because common analytical methods (e.g., liquid chromatography-mass spectrometry) require adducts to be digested or liberated from DNA before quantification, information about their positions within the DNA sequence is lost. Advances in nanopore sequencing technologies allow individual DNA molecules to be analyzed at single-nucleobase resolution, enabling us to study the dynamic of epigenetic modifications and exposure-induced DNA adducts in their native forms on the DNA strand. We applied and evaluated the commercially available Oxford Nanopore Technology (ONT) sequencing platform for site-specific detection of DNA adducts and for distinguishing individual alkylated DNA adducts. Using ONT and the publicly available ELIGOS software, we analyzed a library of 15 plasmids containing site-specifically inserted O6- or N2-alkyl-2'-deoxyguanosine lesions differing in sizes and regiochemistries. Positions of DNA adducts were correctly located, and individual DNA adducts were clearly distinguished from each other.

Toward predictive R-loop computational biology: genome-scale prediction of R-loops reveals their association with complex promoter structures, G-quadruplexes and transcriptionally active enhancers.

R-loops are three-stranded RNA:DNA hybrid structures essential for many normal and pathobiological processes. Previously, we generated a quantitative R-loop forming sequence (RLFS) model, quantitative model of R-loop-forming sequences (QmRLFS) and predicted ∼660 000 RLFSs; most of them located in genes and gene-flanking regions, G-rich regions and disease-associated genomic loci in the human genome. Here, we conducted a comprehensive comparative analysis of these RLFSs using experimental data and demonstrated the high performance of QmRLFS predictions on the nucleotide and genome scales. The preferential co-localization of RLFS with promoters, U1 splice sites, gene ends, enhancers and non-B DNA structures, such as G-quadruplexes, provides evidence for the mechanical linkage between DNA tertiary structures, transcription initiation and R-loops in critical regulatory genome regions. We introduced and characterized an abundant class of reverse-forward RLFS clusters highly enriched in non-B DNA structures, which localized to promoters, gene ends and enhancers. The RLFS co-localization with promoters and transcriptionally active enhancers suggested new models for in cis and in trans regulation by RNA:DNA hybrids of transcription initiation and formation of 3D-chromatin loops. Overall, this study provides a rationale for the discovery and characterization of the non-B DNA regulatory structures involved in the formation of the RNA:DNA interactome as the basis for an emerging quantitative R-loop biology and pathobiology.

Complete genomic and transcriptional landscape analysis using third-generation sequencing: a case study of Saccharomyces cerevisiae CEN.PK113-7D.

Completion of eukaryal genomes can be difficult task with the highly repetitive sequences along the chromosomes and short read lengths of second-generation sequencing. Saccharomyces cerevisiae strain CEN.PK113-7D, widely used as a model organism and a cell factory, was selected for this study to demonstrate the superior capability of very long sequence reads for de novo genome assembly. We generated long reads using two common third-generation sequencing technologies (Oxford Nanopore Technology (ONT) and Pacific Biosciences (PacBio)) and used short reads obtained using Illumina sequencing for error correction. Assembly of the reads derived from all three technologies resulted in complete sequences for all 16 yeast chromosomes, as well as the mitochondrial chromosome, in one step. Further, we identified three types of DNA methylation (5mC, 4mC and 6mA). Comparison between the reference strain S288C and strain CEN.PK113-7D identified chromosomal rearrangements against a background of similar gene content between the two strains. We identified full-length transcripts through ONT direct RNA sequencing technology. This allows for the identification of transcriptional landscapes, including untranslated regions (UTRs) (5' UTR and 3' UTR) as well as differential gene expression quantification. About 91% of the predicted transcripts could be consistently detected across biological replicates grown either on glucose or ethanol. Direct RNA sequencing identified many polyadenylated non-coding RNAs, rRNAs, telomere-RNA, long non-coding RNA and antisense RNA. This work demonstrates a strategy to obtain complete genome sequences and transcriptional landscapes that can be applied to other eukaryal organisms.

R-loopDB: a database for R-loop forming sequences (RLFS) and R-loops.

R-loopDB (http://rloop.bii.a-star.edu.sg) was originally constructed as a collection of computationally predicted R-loop forming sequences (RLFSs) in the human genic regions. The renewed R-loopDB provides updates, improvements and new options, including access to recent experimental data. It includes genome-scale prediction of RLFSs for humans, six other animals and yeast. Using the extended quantitative model of RLFSs (QmRLFS), we significantly increased the number of RLFSs predicted in the human genes and identified RLFSs in other organism genomes. R-loopDB allows searching of RLFSs in the genes and in the 2 kb upstream and downstream flanking sequences of any gene. R-loopDB exploits the Ensembl gene annotation system, providing users with chromosome coordinates, sequences, gene and genomic data of the 1 565 795 RLFSs distributed in 121 056 genic or proximal gene regions of the covered organisms. It provides a comprehensive annotation of Ensembl RLFS-positive genes including 93 454 protein coding genes, 12 480 long non-coding RNA and 7 568 small non-coding RNA genes and 7 554 pseudogenes. Using new interface and genome viewers of R-loopDB, users can search the gene(s) in multiple species with keywords in a single query. R-loopDB provides tools to carry out comparative evolution and genome-scale analyses in R-loop biology.

The interaction between 30b-5p miRNA and MBNL1 mRNA is involved in vascular smooth muscle cell differentiation in patients with coronary atherosclerosis.

Vascular smooth muscle cells (VSMCs) in the arterial wall have diverse functions. In pathological states, the interplay between transcripts and microRNAs (miRNAs) leads to phenotypic changes. Understanding the regulatory role of miRNAs and their target genes may reveal how VSMCs modulate the pathogenesis of coronary artery disease. Laser capture microdissection was performed on aortic wall tissues obtained from coronary artery bypass graft patients with and without recent acute myocardial infarction (MI). The mSMRT-qPCR miRNA assay platform (MiRXES, Singapore) was used to profile miRNA. The miRNA data were co-analyzed with significant mRNA transcripts. TargetScan 7.1 was applied to evaluate miRNA-mRNA interactions. The miRNA profiles of 29 patients (16 MI and 13 non-MI) were evaluated. Thirteen VSMC-related miRNAs were differentially expressed between the MI and non-MI groups. Analysis revealed seven miRNA-targeted mRNAs related to muscular tissue differentiation and proliferation. TargetScan revealed that among the VSMC-related transcripts, MBNL1 had a recognition site that matched the hsa-miR-30b-5p target seed sequence. In addition to predicted analysis, our experiment in vitro with human VSMC culture confirmed that hsa-miR-30b-5p negatively correlated with MBNL1. Our data showed that overexpression of hsa-miR-30b-5p led to downregulation of MBNL1 in VSMCs. This process influences VSMC proliferation and might be involved in VSMC differentiation.

Gene expression profile analysis of aortic vascular smooth muscle cells reveals upregulation of cadherin genes in myocardial infarction patients.

Myocardial infarction (MI) induced by acute coronary arterial occlusion is usually secondary to atherosclerotic plaque rupture. Dysregulated response of vascular smooth muscle cells (VSMCs) in atherosclerotic plaques may promote plaque rupture. Cadherins (CDHs) form adherens junctions and are known stabilizers of atherosclerotic plaques. To date, the expression patterns of cadherin have not been well investigated in MI aortic VSMCs. We aimed to investigate the expression of cadherin genes in the aortic wall of patients with and without MI. Laser capture microdissected VSMCs were obtained from aortic tissue samples of patients undergoing coronary artery bypass graft surgery. Integrative bioinformatic analysis of the microarray profiles of the VSMCs revealed that MI is discriminated at the whole transcriptome level by hundreds of differentially expressed genes, including genes involved in cell adhesion, of which the cadherin superfamily genes were among the top structural category. Eleven significantly deregulated candidates of the cadherin superfamily were chosen and formed a new classifier that collectively discriminated MI vs. non-MI with ~95% accuracy. Significance validation was performed with an independent cohort by quantitative RT-quantitative PCR, confirming overexpression of CDH2, CDH12, PCDH17, and PCDH18 in MI VSMCs. The dysregulation of these cadherin superfamily genes might be related to an MI-induced remote effect on aortic wall VSMCs and to imbalances in signaling pathways and myocardial repair mechanisms. Although pathophysiological significance of our findings requires functional studies, mRNA upregulation of the identified cadherin superfamily members in VSMCs might be associated with the progression of atherosclerosis and angiogenesis activation in MI.

Case of Microcephaly after Congenital Infection with Asian Lineage Zika Virus, Thailand.

We sequenced the virus genomes from 3 pregnant women in Thailand with Zika virus diagnoses. All had infections with the Asian lineage. The woman infected at gestational week 9, and not those infected at weeks 20 and 24, had a fetus with microcephaly. Asian lineage Zika viruses can cause microcephaly.

Genome-Based Comparison of Clostridioides difficile: Average Amino Acid Identity Analysis of Core Genomes.

Infections due to Clostridioides difficile (previously known as Clostridium difficile) are a major problem in hospitals, where cases can be caused by community-acquired strains as well as by nosocomial spread. Whole genome sequences from clinical samples contain a lot of information but that needs to be analyzed and compared in such a way that the outcome is useful for clinicians or epidemiologists. Here, we compare 663 public available complete genome sequences of C. difficile using average amino acid identity (AAI) scores. This analysis revealed that most of these genomes (640, 96.5%) clearly belong to the same species, while the remaining 23 genomes produce four distinct clusters within the Clostridioides genus. The main C. difficile cluster can be further divided into sub-clusters, depending on the chosen cutoff. We demonstrate that MLST, either based on partial or full gene-length, results in biased estimates of genetic differences and does not capture the true degree of similarity or differences of complete genomes. Presence of genes coding for C. difficile toxins A and B (ToxA/B), as well as the binary C. difficile toxin (CDT), was deduced from their unique PfamA domain architectures. Out of the 663 C. difficile genomes, 535 (80.7%) contained at least one copy of ToxA or ToxB, while these genes were missing from 128 genomes. Although some clusters were enriched for toxin presence, these genes are variably present in a given genetic background. The CDT genes were found in 191 genomes, which were restricted to a few clusters only, and only one cluster lacked the toxin A/B genes consistently. A total of 310 genomes contained ToxA/B without CDT (47%). Further, published metagenomic data from stools were used to assess the presence of C. difficile sequences in blinded cases of C. difficile infection (CDI) and controls, to test if metagenomic analysis is sensitive enough to detect the pathogen, and to establish strain relationships between cases from the same hospital. We conclude that metagenomics can contribute to the identification of CDI and can assist in characterization of the most probable causative strain in CDI patients.

QmRLFS-finder: a model, web server and stand-alone tool for prediction and analysis of R-loop forming sequences.

The possible formation of three-stranded RNA and DNA hybrid structures (R-loops) in thousands of functionally important guanine-rich genic and inter-genic regions could suggest their involvement in transcriptional regulation and even development of diseases. Here, we introduce the first freely available R-loop prediction program called Quantitative Model of R-loop Forming Sequence (RLFS) finder (QmRLFS-finder), which predicts RLFSs in nucleic acid sequences based on experimentally supported structural models of RLFSs. QmRLFS-finder operates via a web server or a stand-alone command line tool. This tool identifies and visualizes RLFS coordinates from any natural or artificial DNA or RNA input sequences and creates standards-compliant output files for further annotation and analysis. QmRLFS-finder demonstrates highly accurate predictions of the detected RLFSs, proposing new perspective to further discoveries in R-loop biology, biotechnology and molecular therapy. QmRLFS-finder is freely available at http://rloop.bii.a-star.edu.sg/?pg=qmrlfs-finder.

Quantitative model of R-loop forming structures reveals a novel level of RNA-DNA interactome complexity.

R-loop is the structure co-transcriptionally formed between nascent RNA transcript and DNA template, leaving the non-transcribed DNA strand unpaired. This structure can be involved in the hyper-mutation and dsDNA breaks in mammalian immunoglobulin (Ig) genes, oncogenes and neurodegenerative disease related genes. R-loops have not been studied at the genome scale yet. To identify the R-loops, we developed a computational algorithm and mapped R-loop forming sequences (RLFS) onto 66,803 sequences defined by UCSC as 'known' genes. We found that ∼59% of these transcribed sequences contain at least one RLFS. We created R-loopDB (http://rloop.bii.a-star.edu.sg/), the database that collects all RLFS identified within over half of the human genes and links to the UCSC Genome Browser for information integration and visualisation across a variety of bioinformatics sources. We found that many oncogenes and tumour suppressors (e.g. Tp53, BRCA1, BRCA2, Kras and Ptprd) and neurodegenerative diseases related genes (e.g. ATM, Park2, Ptprd and GLDC) could be prone to significant R-loop formation. Our findings suggest that R-loops provide a novel level of RNA-DNA interactome complexity, playing key roles in gene expression controls, mutagenesis, recombination process, chromosomal rearrangement, alternative splicing, DNA-editing and epigenetic modifications. RLFSs could be used as a novel source of prospective therapeutic targets.