Sanger Sequencing Data Research Articles

Combination of Sanger and target-enrichment markers supports revised generic delimitation in the problematic ‘Urera clade’ of the nettle family (Urticaceae)

Urera Gaudich, s.l. is a pantropical genus comprising c. 35 species of trees, shrubs, and vines. It has a long history of taxonomic uncertainty, and is repeatedly recovered as polyphyletic within a poorly resolved complex of genera in the Urticeae tribe of the nettle family (Urticaceae). To provide generic delimitations concordant with evolutionary history, we use increased taxonomic and genomic sampling to investigate phylogenetic relationships among Urera and associated genera. A cost-effective two-tier genome-sampling approach provides good phylogenetic resolution by using (i) a taxon-dense sample of Sanger sequence data from two barcoding regions to recover clades of putative generic rank, and (ii) a genome-dense sample of target-enrichment data for a subset of representative species from each well-supported clade to resolve relationships among them. The results confirm the polyphyly of Urera s.l. with respect to the morphologically distinct genera Obetia, Poikilospermum and Touchardia. Afrotropic members of Urera s.l. are recovered in a clade sister to the xerophytic African shrubs Obetia; and Hawaiian ones with Touchardia, also from Hawaii. Combined with distinctive morphological differences between Neotropical and African members of Urera s.l., these results lead us to resurrect the previously synonymised name Scepocarpus Wedd. for the latter. The new species epiphet Touchardia oahuensis T.Wells & A.K. Monro is offered as a replacement name for Touchardia glabra non H.St.John, and subgenera are created within Urera s.s. to account for the two morphologically distinct Neotropical clades. This new classification minimises taxonomic and nomenclatural disruption, while more accurately reflecting evolutionary relationships within the group.

Cross-Site Evaluation of Commercial Sanger Sequencing Chemistries.

Sanger sequencing remains an essential tool utilized by researchers. Despite competition from commercial sequencing providers, many academic sequencing core facilities continue to offer these services based on a model of competitive pricing, knowledgeable technical support, and rapid turnaround time. In-house Sanger sequencing remains a viable core service and, until recently, Applied Biosystems BigDye Terminator chemistry was the only commercially available solution for Sanger DNA sequencing on Applied Biosystems (ABI) instruments; however, several new products employing novel dye chemistries and reaction configurations have entered the market. As a result, there is a need to benchmark the performance of these new chemistries on various DNA templates, including difficult-to-sequence templates, and their amenability to commonly employed cost-saving measures, such as dye dilution and reaction miniaturization. To evaluate these new reagents, a study was designed to compare the quality of Sanger sequencing data produced by ABI BigDye and commercially available kits from 2 other vendors using both control and difficult-to-sequence DNA templates under various reaction conditions. This study will serve as a valuable resource to core facilities conducting Sanger sequencing that wish to evaluate the use of an alternative chemistry in their sequencing core.

Molecular epidemiology of mumps viruses in the Netherlands, 2017-2019.

Mumps cases continue to occur, also in countries with a relatively high vaccination rate. The last major outbreaks of mumps in the Netherlands were in 2009–2012 and thereafter, only small clusters and single cases were reported. Molecular epidemiology can provide insights in the circulation of mumps viruses. The aims of the present study were to analyze the molecular epidemiology of mumps viruses in the Netherlands in 2017–2019 and to compare the phylogenetic trees built from sequence data of near complete mumps virus genomes or from the SH gene and non-coding regions (SH+NCRs). To this end, Sanger sequence data from SH+NCRs were analyzed from 82 mumps genotype G viruses. In addition, near complete genomes were obtained from 10 mumps virus isolates using next-generation sequencing. Analysis of SH+NCRs sequences of mumps genotype G viruses revealed the presence of two major genetic lineages in the Netherlands, which was confirmed by analysis of near complete genomes. Comparison of phylogenetic trees built with SH+NCRs or near complete genomes indicated that the topology was similar, while somewhat longer branches were present in the phylogenetic tree with near complete genomes. These results confirm that analysis of SH + NCRs sequence data is a useful approach for molecular surveillance. Furthermore, data from recent mumps genotype G viruses might indicate (intermittent) circulation of mumps genotype G viruses in the Netherlands in 2017–2019.

Genome Resequencing, Improvement of Variant Calling, and Population Genomic Analyses Provide Insights into the Seedlessness in the Genus Vitis.

The seedlessness of grape derived from stenospermocarpy is one of the most prized traits of table or raisin grapes. It is controlled by a complex genetic system containing one dominant gene and multiple recessive genes. Here, we collected dense variation data from high-depth resequencing data of seeded, seedless, and wild relative grape genomes sequenced to > 37x mean depth. Variant calls were made using a modified variant calling pipeline that was suitable for highly diverse interspecific grape accessions. The modified pipeline enabled us to call several million more variants than the commonly recommended pipeline. The quality was validated by Sanger sequencing data and subsequently supported by the genetic population structure and the phylogenetic tree constructed using the obtained variation data, results of which were generally consistent with known pedigree and taxonomic classifications. Variation data enabled us to confirm a dominant gene and identify recessive loci for seedlessness. Incidentally, we found that grape cultivar Rizamat contains an ancestral chromosomal region of the dominant gene in Sultanina, a predominant seedlessness donor cultivar. Furthermore, we predicted new candidate causal genes including Vitvi01g00455, Vitvi08g01528, and Vitvi18g01237 associated with the recessive seedless-regulating loci, which showed high homology with genes that regulate seed development in Arabidopsis. This study provides fundamental insights relevant to variant calling from genome resequencing data of diverse interspecific hybrid germplasms such as grape and will accelerate future efforts aimed at crop improvement.

Selective Chemical Labeling and Sequencing of 5-Carboxylcytosine in DNA at Single-Base Resolution.

5-Carboxylcytosine (5caC) plays a vital role in the dynamics of DNA demethylation, and sequencing of its sites will help us dig out more biological functions of 5caC. Herein, we present a novel chemical method to efficiently label 5caC distinguished from other bases in DNA. Combined with bisulfite sequencing, 5caC sites can be located at single-base resolution, and the efficiency of 5caC labeling is 92% based on the Sanger sequencing data. Furthermore, dot blot assays have confirmed that 5caC-containing DNA isolated from HeLa cells was successfully labeled using our method. We expect that our strategy can be further applied to selectively tagging other carboxyl-modified bases and mapping their sites in RNA.

SOAPTyping: an open-source and cross-platform tool for sequence-based typing for HLA class I and II alleles

BackgroundThe human leukocyte antigen (HLA) gene family plays a key role in the immune response and thus is crucial in many biomedical and clinical settings. Utilizing Sanger sequencing, the golden standard technology for HLA typing enables accurate identification of HLA alleles in high-resolution. However, only the commercial software, such as uTYPE, SBT-Assign, and SBTEngine, and very few open-source tools could be applied to perform HLA typing based on Sanger sequencing.ResultsWe developed a user-friendly, cross-platform and open-source desktop application, known as SOAPTyping, for Sanger-based typing in HLA class I and II alleles. SOAPTyping can produce accurate results with a comprehensible protocol and featured functions. Moreover, SOAPTyping supports a more advanced group-specific sequencing primers (GSSP) module to solve the ambiguous typing results. We used SOAPTyping to analyze 36 samples with known HLA typing from the University of California Los Angeles (UCLA) International HLA DNA Exchange platform and 100 anonymous clinical samples, and the HLA typing results from SOAPTyping are identical to the golden results and 5.5 times faster than commercial software uTYPE, which shows the usability of SOAPTyping.ConclusionsWe introduce the SOAPTyping as the first open-source and cross-platform HLA typing software with the capability of producing high-resolution HLA typing predictions from Sanger sequence data.

Small Cell Carcinoma of the Ovary, Hypercalcemic Type (SCCOHT) beyond SMARCA4 Mutations: A Comprehensive Genomic Analysis.

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is an aggressive malignancy that occurs in young women, is characterized by recurrent loss-of-function mutations in the SMARCA4 gene, and for which effective treatments options are lacking. The aim of this study was to broaden the knowledge on this rare malignancy by reporting a comprehensive molecular analysis of an independent cohort of SCCOHT cases. We conducted Whole Exome Sequencing in six SCCOHT, and RNA-sequencing and array comparative genomic hybridization in eight SCCOHT. Additional immunohistochemical, Sanger sequencing and functional data are also provided. SCCOHTs showed remarkable genomic stability, with diploid profiles and low mutation load (mean, 5.43 mutations/Mb), including in the three chemotherapy-exposed tumors. All but one SCCOHT cases exhibited 19p13.2-3 copy-neutral LOH. SMARCA4 deleterious mutations were recurrent and accompanied by loss of expression of the SMARCA2 paralog. Variants in a few other genes located in 19p13.2-3 (e.g., PLK5) were detected. Putative therapeutic targets, including MAGEA4, AURKB and CLDN6, were found to be overexpressed in SCCOHT by RNA-seq as compared to benign ovarian tissue. Lastly, we provide additional evidence for sensitivity of SCCOHT to HDAC, DNMT and EZH2 inhibitors. Despite their aggressive clinical course, SCCOHT show remarkable inter-tumor homogeneity and display genomic stability, low mutation burden and few somatic copy number alterations. These findings and preliminary functional data support further exploration of epigenetic therapies in this lethal disease.

Copper-Mediated DNA-Compatible One-Pot Click Reactions of Alkynes with Aryl Borates and TMS-N3.

We report a DNA-compatible copper-mediated efficient synthesis of 1,2,3-triazoles via a one-pot reaction of aryl borates with TMS-N3 followed by a click cycloaddition reaction. Employing the binuclear macrocyclic nanocatalyst Cu(II)-β-cyclodextrin, the reactions were performed under mild conditions with high conversions and wide functional group tolerance. We also demonstrate the reaction application toward a one-pot DNA-compatible intramolecular macrocyclization. Our optimized reaction protocol results in no significant DNA damage as judged by qPCR analysis and Sanger sequencing data.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes.

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Complete Chromosomal Sequences of Two Borrelia miyamotoi Samples Obtained from Ixodes ricinus Eggs in Czechia.

Here, we present complete chromosome sequences of Borrelia miyamotoi samples CZ-F1E and CZ-F190E, which were obtained from Ixodes ricinus eggs from Czechia. The chromosome sequences, assembled from Illumina and Sanger sequencing data, had average coverage values of 647× and 3,216×, respectively. They belong to the European genotype, distinct from the Asian and American strains.

Tracy: basecalling, alignment, assembly and deconvolution of sanger chromatogram trace files

BackgroundDNA sequencing is at the core of many molecular biology laboratories. Despite its long history, there is a lack of user-friendly Sanger sequencing data analysis tools that can be run interactively as a web application or at large-scale in batch from the command-line.ResultsWe present Tracy, an efficient and versatile command-line application that enables basecalling, alignment, assembly and deconvolution of sequencing chromatogram files. Its companion web applications make all functionality of Tracy easily accessible using standard web browser technologies and interactive graphical user interfaces. Tracy can be easily integrated in large-scale pipelines and high-throughput settings, and it uses state-of-the-art file formats such as JSON and BCF for reporting chromatogram sequencing results and variant calls. The software is open-source and freely available at https://github.com/gear-genomics/tracy, the companion web applications are hosted at https://www.gear-genomics.com.ConclusionsTracy can be routinely applied in large-scale validation efforts conducted in clinical genomics studies as well as for high-throughput genome editing techniques that require a fast and rapid method to confirm discovered variants or engineered mutations. Molecular biologists benefit from the companion web applications that enable installation-free Sanger chromatogram analyses using intuitive, graphical user interfaces.

Investigating the utility of traditional and genomic multi-locus datasets to resolve relationships in Lipaugus and Tijuca (Cotingidae)

Rapid diversification limits our ability to resolve evolutionary relationships and examine diversification history, as in the case of the Neotropical cotingas. Here we present an analysis with complete taxon sampling for the cotinga genera Lipaugus and Tijuca, which include some of the most range-restricted (e.g., T. condita) and also the most widespread and familiar (e.g., L. vociferans) forest birds in the Neotropics. We used two datasets: (1) Sanger sequencing data sampled from eight loci in 34 individuals across all described taxa and (2) sequence capture data linked to 1,079 ultraconserved elements and conserved exons sampled from one or two individuals per species. Phylogenies estimated from the Sanger sequencing data failed to resolve three nodes, but the sequence capture data produced a well-supported tree. Lipaugus and Tijuca formed a single, highly supported clade, but Tijuca species were not sister and were embedded within Lipaugus. A dated phylogeny confirmed Lipaugus and Tijuca diversified rapidly in the Miocene. Our study provides a detailed evolutionary hypothesis for Lipaugus and Tijuca and demonstrates that increasing genomic sampling can prove instrumental in resolving the evolutionary history of recent radiations.

I-Silence, Please! An Alternative for Gene Disruption via Adenine Base Editors.

i-Silence, Please! An Alternative for Gene Disruption via Adenine Base Editors.

Molecular Diagnosis of FLT3 Mutations in Acute Myeloid Leukemia Patients

Background. FLT3 gene is an important prognostic molecular marker in acute myeloid leukemia (AML). However, the detection of FLT3 mutations presents a challenge. Aim. To compare techniques used for the detection of FLT3 mutations, and to develop a test-system based on polymerase chain reaction (PCR) for quick and reliable determination of FLT3 mutation status. Materials & Methods. Bone marrow samples obtained from AML patients were subjected to examination. To detect FLT3-ITD and FLT3-TKD mutations PCR was performed with subsequent agarose gel electrophoresis visualization. The results were verified by Sanger sequencing. The data obtained using our test-system were compared with widely applied commercial kit ‘FLT3 Mutation Assay for Gel Detection’ by Invivoscribe. Results. To determine the FLT3 mutation status a PCR test was developed. This technique was validated on 22 bone marrow samples obtained from AML patients. FLT3-ITD mutation was detected in 4 patients, 3 patients showed FLT3-TKD mutation. In 1 patient both mutations were identified. These results fully corresponded to the molecular genetic analysis of FLT3, performed by ‘FLT3 Mutation Assay for Gel Detection’. The chosen technique was validated using Sanger sequencing data analysis. Conclusion. The article offers the review of all existing FLT3 mutation screening techniques and describes the experience of developing the PCR test for FLT3-ITD and FLT3-TKD mutation detection. The chosen technique is affordable and easy to use compared with the others. This study can be used as a guide for both doctors and researchers.

From metabarcoding to metaphylogeography: separating the wheat from the chaff.

Metabarcoding is by now a well‐established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α‐ and β‐diversity estimates, that is, interspecies (or inter‐MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra‐MOTU) information—mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra‐MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co‐occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two‐step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra‐ and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.

Long Range Sequencing Shows Improved Resolution in the Detection of Beta Globin Cluster Variants

Introduction: Globin gene sequencing has long been impeded by repetitive regions, homology and structural variation. Currently, long-read sequencing technology has been limited by high cost and error rate. These challenges can be mitigated by sample multiplexing of targeted capture sequencing and the use of consensus reads, respectively. This ongoing study analyzes the utility of long range sequencing technology in identifying historically problematic beta globin cluster mutations. Materials and methods: Long range sequencing was performed on complex beta globin gene cluster cases utilizing the Sequel platform (Pacific Biosciences, Menlo Park, California, United States). Sample input was 1 microgram of DNA extracted from whole blood. Circular Consensus Sequence (CCS) FASTQ sequences files were mapped to reference genome (GRCh37/hg19) using NGMLR, followed by annotation of variants using a custom bioinformatics workflow. Comparison to multiplex ligation-dependent probe amplification (MLPA), Sanger sequencing, array CGH, short-range NGS (300 bp), hemoglobin protein, and clinical data was performed. Results: CCS read length varied from 1 to 10 kb, average 3.5. Analysis of the beta globin cluster region showed impressive superiority and comprehensiveness to Sanger, MLPA, array CGH and short range sequencing technology: 1) single nucleotide variants (SNVs) were identified to a sensitivity and specificity of 99.5%; 2) large structural variants (SVs) such as large deletions, duplications, insertions, crossovers and fusions spanning many kilobases were characterized in the heterozygous, homozygous and compound heterozygous state to a precise genomic coordinate (table 1); 3) phasing SNVs identified Hb S haplotypes (Central African region (CAR), Benin, Senegal, Arab-Indian and Cameroon). Conclusion: 1) Long range sequencing holds significant promise in genotyping hemoglobin disorders. This method improves efficiency with potential for a single test offering with increased resolution, particularly in large deletion analysis. Complex cases greatly benefited from improved range and resolution. 2) The decreased requirement of blood volume for extensive testing is of particular interest for the pediatric population. 3) Genotyping of Hb S patients detected by positive newborn screening will help predict need for prophylaxis or therapy by assessing variant zygosity, phasing and haplotypes. In conclusion, we believe that long range sequencing has immense potential to contribute extensively to thalassemia and hemoglobinopathy diagnostics and is far superior to the currently available technologies for beta globin cluster analysis. Disclosures Jen: Celgene: Employment.

A novel FLCN mutation in family members diagnosed with primary spontaneous pneumothorax.

BackgroundPrimary spontaneous pneumothorax (PSP) is a disease characterized by the accumulation of air in the pleural space between the lung and thoracic wall. It is more common in young, tall, thin, and asthenic men. A family history was reported for approximately 11.5% of individuals admitted with PSP. The literature has reported cases diagnosed with familial PSP, who have no manifestations of Birt–Hogg–Dubé (BHD) syndrome but mutations in different exons of the Folliculin (FLCN) gene. The aim of this study is to present a Turkish family in which 13 members from three generations of the same family developed recurrent isolated spontaneous pneumothorax with a novel mutation in the FLCN.MethodsA male proband was diagnosed with spontaneous pneumothorax in the emergency department of the University of Health Sciences Haydarpasa Numune Training and Research Center, Istanbul, Turkey. His 12 relatives from three generations diagnosed with PSP, as revealed by his family history, were invited to the hospital to give blood samples for mutation analysis. The Sanger sequence data of FLCN were analyzed on the ENSEMBL website using SeqScape 3 and Codon Aligner software.ResultsA novel heterozygous mutation c. 1273C>T (p.Gln425Ter) was detected in exon 11 of the FLCN, which caused PSP in the proband and his 12 relatives tested using Sanger sequencing.ConclusionWe found that a heterozygous mutation in exon 11 of FLCN c. 1273C>T (p.Gln425Ter), which was identified for the first time in our study, might cause isolated familial spontaneous pneumothorax.

Old World and New World Phasmatodea: Phylogenomics Resolve the Evolutionary History of Stick and Leaf Insects

Phasmatodea comprises over 3,000 extant species and stands out as one of the last remaining insect orders for which a robust, higher-level phylogenetic hypothesis is lacking. New research suggests that the extant diversity is the result of a surprisingly recent and rapid radiation that has been difficult to resolve with standard Sanger sequence data. In order to resolve the early branching events of stick and leaf insects, we analyzed transcriptomes from 61 species, including 38 Phasmatodea species comprising all major clades and 23 outgroup taxa, including all other Polyneoptera orders. Using a custom-made ortholog set based on reference genomes from four species, we identified on average 2,274 orthologous genes in the sequenced transcriptomes. We generated various sub-alignments and performed maximum-likelihood analyses on several representative datasets to evaluate the effect of missing data and matrix composition on our phylogenetic estimates. Based on our new data, we are able to reliably resolve the deeper nodes between the principal lineages of extant Phasmatodea. Among Euphasmatodea, we provide strong evidence for a basal dichotomy of Aschiphasmatodea and all remaining euphasmatodeans, the Neophasmatodea. Within the latter clade, we recovered a previously unrecognized major New World and Old World lineage, for which we introduce the new names Oriophasmata tax. nov. (“Eastern phasmids”) and Occidophasmata tax. nov. (“Western phasmids”). Occidophasmata comprise Diapheromerinae, Pseudophasmatinae, and Agathemera, whereas all remaining lineages form the Oriophasmata, including Heteropterygidae, Phylliinae, Bacillus, Lonchodidae (Necrosciinae + Lonchodinae), Clitumninae, Cladomorphinae, and Lanceocercata. We furthermore performed a divergence time analysis and reconstructed the historical biogeography for stick and leaf insects. Phasmatodea either originated in Southeast Asia or in the New World. Our results suggest that the extant distribution of Phasmatodea is largely the result of dispersal events in a recently and rapidly diversified insect lineage rather than the result of vicariant processes.

De Novo Genome Sequence Assemblies of Gossypium raimondii and Gossypium turneri.

Cotton is an agriculturally important crop. Because of its importance, a genome sequence of a diploid cotton species (Gossypium raimondii, D-genome) was first assembled using Sanger sequencing data in 2012. Improvements to DNA sequencing technology have improved accuracy and correctness of assembled genome sequences. Here we report a new de novo genome assembly of G. raimondii and its close relative G. turneri. The two genomes were assembled to a chromosome level using PacBio long-read technology, HiC, and Bionano optical mapping. This report corrects some minor assembly errors found in the Sanger assembly of G. raimondii. We also compare the genome sequences of these two species for gene composition, repetitive element composition, and collinearity. Most of the identified structural rearrangements between these two species are due to intra-chromosomal inversions. More inversions were found in the G. turneri genome sequence than the G. raimondii genome sequence. These findings and updates to the D-genome sequence will improve accuracy and translation of genomics to cotton breeding and genetics.

Identification and characterization of novel type of RNAs, circRNAs in crucian carp Carassius auratus gibelio

Circular RNAs (circRNAs) with regulatory potency activity was identified from varieties of species. Crucian carp (Carassius auratus gibelio) is one of the most freshwater aquaculture species in China. Every year, huge economic damage to the farming was caused by the virus and bacterial infection. Until now, there is any information about circRNA reported from the Crucian carp. In this study, the expression pattern of circRNA in Crucian carp was investigated with transcriptomic analysis. The results showed that only 37 circRNAs were identified from the Crucian carp, and these circRNAs biogenesis was formed with canonical GU-AG splicing mechanism with unevenly distributed on the chromosomes. Wherein, most of the circRNAs were derived from the sense overlapping strategy. Reverse transcript PCR and Sanger sequencing data indicated that these circRNAs were existed authenticity in Crucian carp. The bioinformatics analysis indicated that circRNAs identified from the Crucian carp with potential miRNA sponge regulate the expression level of mRNAs. GO annotation and KEGG pathway analysis of these circRNAs showed that more than 20% circRNAs were related with catalytic activity and binding in the category of molecular function, and these circRNAs were enriched in 9 signaling pathways, such as, Wnt signaling pathway, MAPK signaling pathway, Ubiquitin mediated proteolysis et al. 220 mRNAs would be regulated by the circRNAs via miRNAs mediation. These target mRNAs were further analyzed with functional annotation and KEGG analysis. GO annotation analysis showed that several genes were related with function of nucleotide binding, transcription regulatory activity. KEGG pathway analysis showed that 5 genes were enriched in the pathway of Endocytosis. The circRNA-miRNA-mRNA regulation network indicated that one miRNA can link one or more circRNA and one or more mRNA. Overall, these results will not only help us to further understand the novel RNA transcripts in Crucian carp, but also provide the novel clue to investigate the interaction between host and pathogens from this novel circRNA molecule.