Circular Sequence Research Articles

Improving eukaryotic genome annotation using single molecule mRNA sequencing

BackgroundThe advantages of Pacific Biosciences (PacBio) single-molecule real-time (SMRT) technology include long reads, low systematic bias, and high consensus read accuracy. Here we use these attributes to improve on the genome annotation of the parasitic hookworm Ancylostoma ceylanicum using PacBio RNA-Seq.ResultsWe sequenced 192,888 circular consensus sequences (CCS) derived from cDNAs generated using the CloneTech SMARTer system. These SMARTer-SMRT libraries were normalized and size-selected providing a robust population of expressed structural genes for subsequent genome annotation. We demonstrate PacBio mRNA sequences based genome annotation improvement, compared to genome annotation using conventional sequencing-by-synthesis alone, by identifying 1609 (9.2%) new genes, extended the length of 3965 (26.7%) genes and increased the total genomic exon length by 1.9 Mb (12.4%). Non-coding sequence representation (primarily from UTRs based on dT reverse transcription priming) was particularly improved, increasing in total length by fifteen-fold, by increasing both the length and number of UTR exons. In addition, the UTR data provided by these CCS allowed for the identification of a novel SL2 splice leader sequence for A. ceylanicum and an increase in the number and proportion of functionally annotated genes. RNA-seq data also confirmed some of the newly annotated genes and gene features.ConclusionOverall, PacBio data has supported a significant improvement in gene annotation in this genome, and is an appealing alternative or complementary technique for genome annotation to the other transcript sequencing technologies.

Preparation of Circular RNA In Vitro.

This chapter describes a simple and straightforward way to obtain single-stranded circular RNA sequences in vitro. Linear RNA that is phosphorylated at the 5' end is first prepared by a chemical or enzymatic method, then circularized using ligase. The function of the prepared circular RNA molecule, such as an ability to induce translation, can then be investigated.

Pacific Biosciences Sequencing and IMGT/HighV-QUEST Analysis of Full-Length Single Chain Fragment Variable from an In Vivo Selected Phage-Display Combinatorial Library.

Phage-display selection of immunoglobulin (IG) or antibody single chain Fragment variable (scFv) from combinatorial libraries is widely used for identifying new antibodies for novel targets. Next-generation sequencing (NGS) has recently emerged as a new method for the high throughput characterization of IG and T cell receptor (TR) immune repertoires both in vivo and in vitro. However, challenges remain for the NGS sequencing of scFv from combinatorial libraries owing to the scFv length (>800 bp) and the presence of two variable domains [variable heavy (VH) and variable light (VL) for IG] associated by a peptide linker in a single chain. Here, we show that single-molecule real-time (SMRT) sequencing with the Pacific Biosciences RS II platform allows for the generation of full-length scFv reads obtained from an in vivo selection of scFv-phages in an animal model of atherosclerosis. We first amplified the DNA of the phagemid inserts from scFv-phages eluted from an aortic section at the third round of the in vivo selection. From this amplified DNA, 450,558 reads were obtained from 15 SMRT cells. Highly accurate circular consensus sequences from these reads were generated, filtered by quality and then analyzed by IMGT/HighV-QUEST with the functionality for scFv. Full-length scFv were identified and characterized in 348,659 reads. Full-length scFv sequencing is an absolute requirement for analyzing the associated VH and VL domains enriched during the in vivo panning rounds. In order to further validate the ability of SMRT sequencing to provide high quality, full-length scFv sequences, we tracked the reads of an scFv-phage clone P3 previously identified by biological assays and Sanger sequencing. Sixty P3 reads showed 100% identity with the full-length scFv of 767 bp, 53 of them covering the whole insert of 977 bp, which encompassed the primer sequences. The remaining seven reads were identical over a shortened length of 939 bp that excludes the vicinity of primers at both ends. Interestingly these reads were obtained from each of the 15 SMRT cells. Thus, the SMRT sequencing method and the IMGT/HighV-QUEST functionality for scFv provides a straightforward protocol for characterization of full-length scFv from combinatorial phage libraries.

A novel Tn1696-like composite transposon (Tn6404) harboring blaIMP-4 in a Klebsiella pneumoniae isolate carrying a rare ESBL gene blaSFO-1

Genetic determinants of a clinical Klebsiella pneumoniae isolate (KP1814) coproducing IMP-4 and a rare ESBL gene SFO-1 was investigated. KP1814 belongs to a novel sequence type (ST) assigned to ST2270. WGS identified four circular DNA sequences in KP1814, including two multidrug-resistance (MDR) plasmids, one virulence plasmid, and one circular form. The MDR plasmid pKP1814-1 (299.9 Kb) is untypeable, and carries two large mosaic multiresistance regions (MRRs). blaSFO-1 and blaIMP-4 co-exists on MRR1, and blaSFO-1 is associated with an IS/Tn-independent genetic context. blaIMP-4 is carried by a novel In804-like integron (intlI-blaIMP-4-Kl.pn.I3-qacG2-aacA4-catB3∆) associated with a novel Tn1696-like transposon (designed Tn6404) flanked by IS5075. The other MDR plasmid pKP1814-3 is a 95,701-bp IncFII plasmid, and is a hybrid of a Shigella flexneri plasmid pSF07201 and an E. coli plasmid pCA08. All resistance genes of pKP1814-3 were detected in a ~16-kb IS26-flanked composite transposon carried by a Tn5396 transposon. The circular form (18.3 Kb) was composed of two parts belonging to pKP1814-1 and pKP1814-3, respectively. The plasmid pKP1814-2, carrying multiple virulence factors, encodes IncFIBK and IncFIIK replicons with a size of 187,349 bp. The coexistence of MDR and virulence plasmids largely enhances the bacterial fitness in the host and environment.

Genome Sequence of Bacillus velezensis S141, a New Strain of Plant Growth-Promoting Rhizobacterium Isolated from Soybean Rhizosphere.

ABSTRACTBacillus velezensis strain S141 is a plant growth-promoting rhizobacterium isolated from soybean (Glycine max) rhizosphere that enhances soybean growth, nodulation, and N2 fixation efficiency by coinoculation with Bradyrhizobium diazoefficiens USDA110. The S141 genome was identified to comprise a 3,974,582-bp-long circular DNA sequence encoding at least 3,817 proteins.

Profust reliability of linear and circular type F and G systems having two failure criteria under Markov dependency

ABSTRACTIn this study, the profust reliabilities of (n, f, k): F(G) and < n, f, k > : F(G) systems for Markov dependent components are investigated. Having two failure criteria are the common features of these four systems. The usage of both fuzzy approach and two failure criteria in the same system provides us more realistic approach to evaluate the reliability of more complex systems. The component configurations are examined for both linear and circular sequences and the working principle of systems are studied for both F and G systems. Under all these assumptions, the profust reliabilities of (n, f, k): F(G) and < n, f, k > : F(G) systems are obtained using the distribution of run statistics. Also a new membership function different from the linear membership function which is generally used in the literature is proposed. Some numerical results which allow the comparison of the systems from various perspectives and various figures for both linear and circular type systems are presented. Some special cases (between Markov – iid assumption, conventional – profust reliability) are also considered.

Parallel approach on sorting of genes in search of optimal solution

An important tool for comparing genome analysis is the rearrangement event that can transform one given genome into other. For finding minimum sequence of fission and fusion, we have proposed here an algorithm and have shown a transformation example for converting the source genome into the target genome. The proposed algorithm comprises of circular sequence i.e. “cycle graph” in place of mapping. The main concept of algorithm is based on optimal result of permutation. These sorting processes are performed in constant running time by showing permutation in the form of cycle. In biological instances it has been observed that transposition occurs half of the frequency as that of reversal. In this paper we are not dealing with reversal instead commencing with the rearrangement of fission, fusion as well as transposition.

Complete Circular Genome Sequence and Temperature Independent Adaptation to Anaerobiosis of Listeria weihenstephanensis DSM 24698.

The aim of this study was to analyze the adaptation of the environmental Listeria weihenstephanensis DSM 24698 to anaerobiosis. The complete circular genome sequence of this species is reported and the adaptation of L. weihenstephanensis DSM 24698 to oxygen availability was investigated by global transcriptional analyses via RNAseq at 18 and 34°C. A list of operons was created based on the transcriptional data. Forty-two genes were upregulated anaerobically and 62 genes were downregulated anaerobically. The oxygen dependent gene expression of selected genes was further validated via qPCR. Many of the differentially regulated genes encode metabolic enzymes indicating broad metabolic adaptations with respect to oxygen availability. Genes showing the strongest oxygen-dependent adaption encoded nitrate (narGHJI) and nitrite (nirBD) reductases. Together with the observation that nitrate supported anaerobic growth, these data indicate that L. weihenstephanensis DSM 24698 performs anaerobic nitrate respiration. The wide overlap between the oxygen-dependent transcriptional regulation at 18 and 34°C suggest that temperature does not play a key role in the oxygen-dependent transcriptional regulation of L. weihenstephanensis DSM 24698.

Analyses of Mitogenome Sequences Revealed that Asian Citrus Psyllids (Diaphorina citri) from California Were Related to Those from Florida

Asian citrus psyllid (ACP, Diaphorina citri Kuwayama) transmits “Candidatus Liberibacter asiaticus” (CLas), an unculturable alpha-proteobacterium associated with citrus Huanglongbing (HLB). CLas has recently been found in California. Understanding ACP population diversity is necessary for HLB regulatory practices aimed at reducing CLas spread. In this study, two circular ACP mitogenome sequences from California (mt-CApsy, ~15,027 bp) and Florida (mt-FLpsy, ~15,012 bp), USA, were acquired. Each mitogenome contained 13 protein coding genes, 2 ribosomal RNA and 22 transfer RNA genes, and a control region varying in sizes. The Californian mt-CApsy was identical to the Floridian mt-FLpsy, but different from the mitogenome (mt-GDpsy) of Guangdong, China, in 50 single nucleotide polymorphisms (SNPs). Further analyses were performed on sequences in cox1 and trnAsn regions with 100 ACPs, SNPs in nad1-nad4-nad5 locus through PCR with 252 ACP samples. All results showed the presence of a Chinese ACP cluster (CAC) and an American ACP cluster (AAC). We proposed that ACP in California was likely not introduced from China based on our current ACP collection but somewhere in America. However, more studies with ACP samples from around the world are needed. ACP mitogenome sequence analyses will facilitate ACP population research.

On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data.

To benchmark algorithms for automated plasmid sequence reconstruction from short-read sequencing data, we selected 42 publicly available complete bacterial genome sequences spanning 12 genera, containing 148 plasmids. We predicted plasmids from short-read data with four programs (PlasmidSPAdes, Recycler, cBar and PlasmidFinder) and compared the outcome to the reference sequences. PlasmidSPAdes reconstructs plasmids based on coverage differences in the assembly graph. It reconstructed most of the reference plasmids (recall=0.82), but approximately a quarter of the predicted plasmid contigs were false positives (precision=0.75). PlasmidSPAdes merged 84 % of the predictions from genomes with multiple plasmids into a single bin. Recycler searches the assembly graph for sub-graphs corresponding to circular sequences and correctly predicted small plasmids, but failed with long plasmids (recall=0.12, precision=0.30). cBar, which applies pentamer frequency analysis to detect plasmid-derived contigs, showed a recall and precision of 0.76 and 0.62, respectively. However, cBar categorizes contigs as plasmid-derived and does not bin the different plasmids. PlasmidFinder, which searches for replicons, had the highest precision (1.0), but was restricted by the contents of its database and the contig length obtained from de novo assembly (recall=0.36). PlasmidSPAdes and Recycler detected putative small plasmids (<10 kbp), which were also predicted as plasmids by cBar, but were absent in the original assembly. This study shows that it is possible to automatically predict small plasmids. Prediction of large plasmids (>50 kbp) containing repeated sequences remains challenging and limits the high-throughput analysis of plasmids from short-read whole-genome sequencing data.

Draft Genome Sequence of Thiohalobacter thiocyanaticus Strain FOKN1, a Neutrophilic Halophile Capable of Thiocyanate Degradation.

ABSTRACTA draft genome sequence of a neutrophilic halophile capable of thiocyanate degradation, Thiohalobacter thiocyanaticus FOKN1, was determined using a PacBio RSII sequencer. A 3.23-Mb circular genome sequence was assembled, in which 3,026 gene-coding sequences, 45 tRNAs, and 1 rrn operon were annotated.

Retraction Note to: The complete mitochondrial genomes of two freshwater snails provide new protein-coding gene rearrangement models and phylogenetic implications

Mitochondrial (mt) genome sequences are widely used for species identification and to study the phylogenetic relationships among Gastropoda. However, to date, limited data are available as taxon sampling is narrow. In this study we sequenced the complete mt genomes of the freshwater gastropods Radix swinhoei (Lymnaeidae) and Planorbarius corneus (Planorbidae). Based on these sequences, we investigated the gene rearrangement in these two species and the relationships with respect to the ancestral gene order and assessed their phylogenetic relationships. The complete mt genomes of R. swinhoei and P. corneus were sequenced using Illumina-based paired-end sequencing and annotated by comparing the sequence information with that of related gastropod species. Putative models of mitochondrial gene rearrangements were predicted for both R. swinhoei and P. corneus, using Reishia clavigera mtDNA structure as the ancestral gene order. The phylogenetic relationships were inferred using thirteen protein sequences based on Maximum likelihood and Bayesian inference analyses. The complete circular mt genome sequences of R. swinhoei and P. corneus were 14,241 bp and 13,687 bp in length, respectively. Comparison of the gene order demonstrated complex rearrangement events in Gastropoda, both for tRNA genes and protein-coding genes. The phylogenetic analyses showed that the family Lymnaeidae was more closely related to the family Planorbidae, consistent with previous classification. Nevertheless, due to the position recovered for R. swinhoei, the family Lymnaeidae was not monophyletic. This study provides the complete mt genomes of two freshwater snails, which will aid the development of useful molecular markers for epidemiological, ecological and phylogenetic studies. Additionally, the predicted models for mt gene rearrangement might provide novel insights into mt genome evolution in gastropods.

First report of two complete Clostridium chauvoei genome sequences and detailed in silico genome analysis

Clostridium (C.) chauvoei is a Gram-positive, spore forming, anaerobic bacterium. It causes black leg in ruminants, a typically fatal histotoxic myonecrosis. High quality circular genome sequences were generated for the C. chauvoei type strain DSM 7528T (ATCC 10092T) and a field strain 12S0467 isolated in Germany. The origin of replication (oriC) was comparable to that of Bacillus subtilis in structure with two regions containing DnaA boxes. Similar prophages were identified in the genomes of both C. chauvoei strains which also harbored hemolysin and bacterial spore formation genes. A CRISPR type I-B system with limited variations in the repeat number was identified. Sporulation and germination process related genes were homologous to that of the Clostridia cluster I group but novel variations for regulatory genes were identified indicative for strain specific control of regulatory events. Phylogenomics showed a higher relatedness to C. septicum than to other so far sequenced genomes of species belonging to the genus Clostridium. Comparative genome analysis of three C. chauvoei circular genome sequences revealed the presence of few inversions and translocations in locally collinear blocks (LCBs). The species genome also shows a large number of genes involved in proteolysis, genes for glycosyl hydrolases and metal iron transportation genes which are presumably involved in virulence and survival in the host. Three conserved flagellar genes (fliC) were identified in each of the circular genomes. In conclusion this is the first comparative analysis of circular genomes for the species C. chauvoei, enabling insights into genome composition and virulence factor variation.

Molecular and phylogenetic identification of unique isolates of hammerhead viroid-like RNA from ‘Pacific Gala’ apple (Malus domestica) in Canada

Ten isolates of apple hammerhead viroid-like RNA (AHVd RNA) were detected by next-generation sequencing (NGS) and/or by RT-PCR analysis of apple (Malus domestica) ‘Pacific Gala’ (PG) from Canada. All consisted of 433 nucleotides in size, and their circular structure and sequence were supported by reverse transcription amplification with two separate pairs of abutted primers, cloning and sequencing. Two regions were defined for the computer-predicted secondary structure of AHVd RNA-PG isolates – a hammerhead region (HHR) and a loop-rich region (LRR). The HHR contained the core and flanking nucleotides typical of natural hammerhead ribozymes. The secondary structure associated with the LRR was distinct from that predicted for the corresponding region of AHVd RNA isolates previously described in ‘Fuji’ apple in China. However, both structures share conserved stem-loops in the LRR, identified as S1, S2, S3 and S4.

A faster and more accurate heuristic for cyclic edit distance computation

Sequence comparison is the core computation of many applications involving textual representations of data. Edit distance is the most widely used measure to quantify the similarity of two sequences. Edit distance can be defined as the minimal total cost of a sequence of edit operations to transform one sequence into the other; for a sequence x of length m and a sequence y of length n, it can be computed in time O(mn). In many applications, it is common to consider sequences with circular structure: for instance, the orientation of two images or the leftmost position of two linearised circular DNA sequences may be irrelevant. To this end, an algorithm to compute the cyclic edit distance in time O(mnlogm) was proposed (Maes, 2003 [18]) and several heuristics have been proposed to speed up this computation. Recently, a new algorithm based on q-grams was proposed for circular sequence comparison (Grossi et al., 2016 [13]). We extend this algorithm for cyclic edit distance computation and show that this new heuristic is faster and more accurate than the state of the art. The aim of this letter is to give visibility to this idea in the pattern recognition community.

MARS: improving multiple circular sequence alignment using refined sequences

BackgroundA fundamental assumption of all widely-used multiple sequence alignment techniques is that the left- and right-most positions of the input sequences are relevant to the alignment. However, the position where a sequence starts or ends can be totally arbitrary due to a number of reasons: arbitrariness in the linearisation (sequencing) of a circular molecular structure; or inconsistencies introduced into sequence databases due to different linearisation standards. These scenarios are relevant, for instance, in the process of multiple sequence alignment of mitochondrial DNA, viroid, viral or other genomes, which have a circular molecular structure. A solution for these inconsistencies would be to identify a suitable rotation (cyclic shift) for each sequence; these refined sequences may in turn lead to improved multiple sequence alignments using the preferred multiple sequence alignment program.ResultsWe present MARS, a new heuristic method for improving Multiple circular sequence Alignment using Refined Sequences. MARS was implemented in the C++ programming language as a program to compute the rotations (cyclic shifts) required to best align a set of input sequences. Experimental results, using real and synthetic data, show that MARS improves the alignments, with respect to standard genetic measures and the inferred maximum-likelihood-based phylogenies, and outperforms state-of-the-art methods both in terms of accuracy and efficiency. Our results show, among others, that the average pairwise distance in the multiple sequence alignment of a dataset of widely-studied mitochondrial DNA sequences is reduced by around 5% when MARS is applied before a multiple sequence alignment is performed.ConclusionsAnalysing multiple sequences simultaneously is fundamental in biological research and multiple sequence alignment has been found to be a popular method for this task. Conventional alignment techniques cannot be used effectively when the position where sequences start is arbitrary. We present here a method, which can be used in conjunction with any multiple sequence alignment program, to address this problem effectively and efficiently.

RETRACTED ARTICLE: The complete mitochondrial genomes of two freshwater snails provide new protein-coding gene rearrangement models and phylogenetic implications

BackgroundMitochondrial (mt) genome sequences are widely used for species identification and to study the phylogenetic relationships among Gastropoda. However, to date, limited data are available as taxon sampling is narrow. In this study we sequenced the complete mt genomes of the freshwater gastropods Radix swinhoei (Lymnaeidae) and Planorbarius corneus (Planorbidae). Based on these sequences, we investigated the gene rearrangement in these two species and the relationships with respect to the ancestral gene order and assessed their phylogenetic relationships.MethodsThe complete mt genomes of R. swinhoei and P. corneus were sequenced using Illumina-based paired-end sequencing and annotated by comparing the sequence information with that of related gastropod species. Putative models of mitochondrial gene rearrangements were predicted for both R. swinhoei and P. corneus, using Reishia clavigera mtDNA structure as the ancestral gene order. The phylogenetic relationships were inferred using thirteen protein sequences based on Maximum likelihood and Bayesian inference analyses.ResultsThe complete circular mt genome sequences of R. swinhoei and P. corneus were 14,241 bp and 13,687 bp in length, respectively. Comparison of the gene order demonstrated complex rearrangement events in Gastropoda, both for tRNA genes and protein-coding genes. The phylogenetic analyses showed that the family Lymnaeidae was more closely related to the family Planorbidae, consistent with previous classification. Nevertheless, due to the position recovered for R. swinhoei, the family Lymnaeidae was not monophyletic.ConclusionThis study provides the complete mt genomes of two freshwater snails, which will aid the development of useful molecular markers for epidemiological, ecological and phylogenetic studies. Additionally, the predicted models for mt gene rearrangement might provide novel insights into mt genome evolution in gastropods.

Radar Communication Integrated Waveform Design Based on OFDM and Circular Shift Sequence

Based on the orthogonal frequency division multiplexing (OFDM) technique, an intelligent waveform is designed, which is suitable for simultaneously performing data transmission and radar sensing. In view of the inherent high peak-to-mean envelope power ratio (PMERP) and poor peak-to-side-lobe ratio (PSLR) problems in the OFDM based radar and communication (RadCom) waveform design, we propose two technologies to deal with that. To be specific, we adopt Gray code technology to reduce the PMERP and simultaneously choose an optimal cyclic sequence to improve the PSLR of RadCom waveform. In our method, the optimal cyclic sequence is dynamically generated to continuously provide the best waveform according to the change of communication data. In addition, to meet the requirements of different radar detection tasks, two simple methods are utilized to adjust the bandwidth of RadCom waveform. To verify the advantages of the designed waveform, we conduct several simulation experiments to compare with some existing RadCom waveforms. The results show that our designed RadCom waveform can simultaneously achieve lower PMERP and higher PSLR. In addition, our designed RadCom waveform has a thumbtack type fuzzy function and shows the good ability to do multitarget detection.

The complete mitochondrial genome of Paramesotriton chinensis (Caudata: Salamandridae) and phylogenetic studies of Paramesotriton

The complete mitochondrial genome of the Asian warty newt Paramesotriton chinensis was sequenced. The complete mitogenome of P. chinensis is a circular double-stranded DNA sequence that is 16,361 bp long and was biased toward A + T content at 61.3% (33.0% A, 28.3% T, 23.9% C, and 14.7% G). The complete mitogenome of P. chinensis consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 1 ribosomal RNAs (16S rRNA), and 1 putative control region. This study presented the complete mitogenome of P. chinensis and provided essential and important DNA molecular data for further phylogenetic and evolutionary analysis for genus Paramesotriton.

Inverse PCR for Point Mutation Introduction.

Inverse PCR is a powerful tool for the rapid introduction of desired mutations at desired positions in a circular double-stranded DNA sequence. Here, custom-designed mutant primers oriented in the inverse direction are used to amplify the entire circular template with incorporation of the required mutation(s). By careful primer design it can be used to perform such diverse modifications as the introduction of point mutations and multiple mutations, the insertion of new sequences, and even sequence deletions. Three primer formats are commonly used; nonoverlapping, partially overlapping and fully overlapping primers, and here we describe the use of nonoverlapping primers for introduction of a point mutation. Use of such a primer setup in the PCR reaction, with one of the primers containing the desired mismatch mutation, results in the amplification of a linear, double-stranded, mutated product. Methylated template DNA is removed from the nonmethylated PCR product by DpnI digestion and the PCR product is then phosphorylated by polynucleotide kinase treatment before being recircularized by ligation, and transformed to E. coli. This relatively simple site-directed mutagenesis procedure is of major importance in biology and biotechnology today where it is commonly employed for the study and engineering of DNA, RNA, and proteins.