Generation Sequencing Reads Research Articles

First draft reference genome and annotation of the alternative oil species Physaria fendleri.

In the wake of increasing demand for renewable energy sources, plant-based sources including alternative oilseeds have come to the forefront of interest. Hydroxy fatty acids (HFAs), produced in a few oilseed species, are important chemical feed stocks for industrial applications. An integrated approach was taken to assemble the first draft genome of the alternative HFA producer Physaria fendleri (n = 6), an outcrossing species with high heterozygosity. Both de novo transcriptome assemblies and genome assemblies were produced with public and generated sequencing reads. Resulting intermediate assemblies were then scaffolded and patched with multiple data sources, followed by super-scaffolding onto a masked genome of Camelina laxa (n = 6). Despite a current lack of available resources for the physical mapping of genomic scaffolds of P. fendleri, topography of the genome with respect to repeat and gene content was preserved at the scaffold level and not significantly lost via super-scaffolding. Read representation, gene and genome completion statistics, and annotation results illustrated the creation of a functional draft genome and a tool for future research on alternative oil species.

Clinical characteristics and assessment of risk factors in patients with influenza A-induced severe pneumonia after the prevalence of SARS-CoV-2.

The aim of this study is to describe the novel epidemiological and clinical characteristics of influenza A-induced severe pneumonia occurring after the prevalence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and to further assess its potential risk factors for mortality. We retrospectively studied the consecutive case series of 30 patients with confirmed influenza A-induced severe pneumonia treated in the intensive care unit at Dazhou Central Hospital in Sichuan, China, from March 1 to April 30, 2023. Logistic regression was used to analyze the independent risk factors, and receiver operating characteristic (ROC) curves were applied to evaluate the predictive efficacy of associated risk factors for mortality. The mortality rate was 33.3% in this study. Independent risk factors for mortality of patients were acute respiratory distress syndrome (ARDS) (p = 0.044) and septic shock (p = 0.012). ROC statistics for ARDS and septic shock to predict mortality in patients with influenza A-induced severe pneumonia demonstrated an area under the curve of 0.800 (sensitivity 80.0%, specificity 80.0%) and 0.825 (sensitivity 70.0%, specificity 95.0%), respectively. ARDS and septic shock were the independent risk factors for mortality in patients with influenza A-induced severe pneumonia following the end of the SARS-CoV-2 pandemic. But high level of next generation sequencing reads Aspergillus coinfection, and comorbidities did not increase death risk of the study population.

Genome Resequencing of the Honeybee Apis mellifera jemenetica (Hymenoptera: Apidae): A Key Tool towards Characterization, Conservation, and Genomic Selection

We report the whole-genome sequence of the Arabian honeybee (Apis mellifera jemenetica). Seven A. m. jemenetica samples were sequenced representing three distinct subpopulations. Generated sequence reads were mapped to the reference honeybee Apis mellifera genome (Amel_HAv3.1). Data revealed genome-wide patterns of genetic variation which can be useful in the characterization and assessment of positive selection of the Arabian honeybee using different genetic markers. In total, 75.16 Gb of clean bases were generated, and the GC content of samples ranged between 31.9 and 35.3%. The effective reference genome size is 223,937,270 bp. The mapping rate of samples varied from 88.97% to 96.19%, and the effective mapping depth was between 41.80 and 48.84X. Single-nucleotide polymorphisms (SNPs) among sequenced individuals ranged between 2379499 and 2396116 with respect to the reference A. mellifera genome (Amel_HAv3.1), and 2% of the SNPs were nonsynonymous. Genome Analysis Toolkit (GATK) detected 1097962–1109829 InDels and 10090–11962 structural variations (SV) from which 22.1 to 33.8% were in the form of deletions. Copy number variation (CNV) ranged between 550 and 2824, and 45–91% of them were downregulated. These variations among interbreeding individuals or groups of the same species may reflect an adaptive environmental response and fitness among different subpopulations and can be very useful for subspecies characterization, conservation, and selection of the Arabian honeybee.

Li-Fraumeni Syndrome: The Two Faces of a Coin in Myeloid Malignancies

Background Li-Fraumeni Syndrome (LFS) is a rare, autosomal dominant genetic disorder caused by germline mutations in the TP53 tumor suppressor gene. It is characterized by a predisposition to a wide range of cancers, including hematological malignancies. Hematological malignancies account for approximately 3-4% of all LFS-related cancers of which are acute lymphoblastic leukemia (ALL), acute myeloid leukemia (AML), Myelodysplastic syndromes (MDS) and lymphomas. Method and Results We retrospectively analyzed our clinical data for patients with hematological malignances from the period of Jan 2016-Jan 2022 at the national center for cancer care and research (NCCCR), the only tertiary hematology center in the country. We report on two cases with Li Fraumeni syndrome with acute leukemia and their variable spectrum of the disease manifestations and prognosis. Case 1: 16 years old male known case of rheumatic heart disease presented with double malignancy; left mandible osteosarcoma and lymphoblastic leukemia/ lymphoma concurrently in 2020. For his Leukemia he received BFM*2009 protocol and after remission he underwent surgical resection and received radiotherapy for the osteosarcoma. Patient has family history of breast cancer in his maternal aunt and his paternal great aunt. Germline Genetic testing from buccal swab revealed a pathogenic variant in TP53 gene confirming Li Fraumeni syndrome which is most likely De novo as it was not detected in both parents and his 3 siblings. In 2022 patient presented with pancytopenia and therapy related AML/MDS with complex Karyotype. He was started on salvage therapy and was planned for allogenic stem cell transplant (SCT) from his HLA matched brother, however, unfortunately, patient didn't achieve complete remission thus could not receive the SCT and he passed away one year later in 2023. Case 2: 60-year-old female known case of hypertension who was diagnosed with right breast cancer in 2011 for which she underwent surgery, chemotherapy, and radiotherapy. A year later, patient was diagnosed with left upper lobe squamous cell lung carcinoma, and underwent left upper lobectomy, chemotherapy and radiotherapy. Patient was in remission for 6 years, until 2018 when she presented with pancytopenia and therapy related AML with complex monosomic Karyotype. Due to patient's age and comorbidities, she was unfit for intensive therapy, thus she received hypomethylation agent (Azacitidine) till disease progression and was referred to genetic counseling. Patient has a family history of unknown cancer in her paternal half-sister at the age of 77 and an unknown cancer in her paternal cousin. Genetic testing reveled a mosaic pathogenic variant in TP53 gene present in 17% of the 1137 next generation sequencing reads. Shortly later, patient passed away without confirmatory fibroblast testing. Conclusion: Genetic testing indications for Li-Fraumeni syndrome present challenges as patients with full mutationin TP53 often but not always present with the full spectrum of the disease, whereas patients with mosaic profilepresent with a later onset of the disease and might not be identified early. In both scenarios, patients can be miss-identified and exposed to radiotherapy thus testing criteria warrants further evaluation. Acute Leukemia in Li-Fraumeni syndrome is associated with poor prognosis and usually presents in the context of therapy related myeloid malignancies . Identification of patients with Li-Fraumeni syndrome early on is crucial to avoid exposures to radiotherapy whenever possible thus reducing the risk of subsequent malignancies and the tragic outcomes of the disease. There are unmet clinical needs in Li-Fraumeni syndrome diagnostic and therapeutic approaches which needs further attention.

Identification of odorant receptor protein genes in two Cephid stem borers (Hymenoptera: Cephidae) by high-throughput sequencing

Insects are well adapted organisms to the terrestrial life on Earth. The evolution of the odorant receptor family is one of the causes underpinning this remarkable adaptation. Odorant receptors (ORs) sense aromas in the environment and cause the insect to respond. The ability of phytophagous insects to detect odor signals from their hosts is crucial for mating, oviposition, and feeding. The family of odorant receptor genes in Cephidae, pest on some economically important plants, is little understood. Bioinformatic tools were used to analyze the genomic data of the two pest species, Syrista parreyssii (Spinola, 1843) (Hymenoptera: Cephidae) (a rose pest) and Pachycephus smyrnensis J.P.E.F. Stein, 1876 (Hymenoptera: Cephidae), (a poppy pest), to determine their odorant receptors. The whole genome sequencing of P. smyrnensis collected in Sivas in 2020 was performed by next generation sequencing and short reads of S. parreyssii genome were obtained from previous studies. Following bioinformatic analyses, 67 and 82 putative odorant receptor genes were identified and annotated for P. smyrnensis and S. parreyssii, respectively. The ORs of these two species were found to be organized as repetitive genes in five separate clusters. No species-specific OR genes were identified in any of the investigated species. As a result, it was hypothesized that host specificity was acquired through the combined effect of multiple ORs.

In Silico Evaluation of Variant Calling Methods for Bacterial Whole-Genome Sequencing Assays.

Identification and analysis of clinically relevant strains of bacteria increasingly relies on whole-genome sequencing. The downstream bioinformatics steps necessary for calling variants from short-read sequences are well-established but seldom validated against haploid genomes. We devised an in silico workflow to introduce single nucleotide polymorphisms (SNP) and indels into bacterial reference genomes, and computationally generate sequencing reads based on the mutated genomes. We then applied the method to Mycobacterium tuberculosis H37Rv, Staphylococcus aureus NCTC 8325, and Klebsiella pneumoniae HS11286, and used the synthetic reads as truth sets for evaluating several popular variant callers. Insertions proved especially challenging for most variant callers to correctly identify, relative to deletions and single nucleotide polymorphisms. With adequate read depth, however, variant callers that use high quality soft-clipped reads and base mismatches to perform local realignment consistently had the highest precision and recall in identifying insertions and deletions ranging from1 to 50 bp. The remaining variant callers had lower recall values associated with identification of insertions greater than 20 bp.

Whole genome sequencing of simmental cattle for SNP and CNV discovery

BackgroudThe single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity.ResultsIn this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 − 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility.ConclusionsIn conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Virome diversity of Hyalomma dromedarii ticks collected from camels in the United Arab Emirates.

Viruses are important components of the microbiome of ticks. Ticks are capable of transmitting several serious viral diseases to humans and animals. Hitherto, the composition of viral communities in Hyalomma dromedarii ticks associated with camels in the United Arab Emirates (UAE) remains unexplored. This study aimed to characterize the RNA virome diversity in male and female H. dromedarii ticks collected from camels in Al Ain, UAE. We collected ticks, extracted, and sequenced RNA, using Illumina (NovaSeq 6000) and Oxford Nanopore (MinION). From the total generated sequencing reads, 180,559 (~0.35%) and 197,801 (~0.34%) reads were identified as virus-related reads in male and female tick samples, respectively. Taxonomic assignment of the viral sequencing reads was accomplished based on bioinformatic analyses. Further, viral reads were classified into 39 viral families. Poxiviridae, Phycodnaviridae, Phenuiviridae, Mimiviridae, and Polydnaviridae were the most abundant families in the tick viromes. Notably, we assembled the genomes of three RNA viruses, which were placed by phylogenetic analyses in clades that included the Bole tick virus. Overall, this study attempts to elucidate the RNA virome of ticks associated with camels in the UAE and the results obtained from this study improve the knowledge of the diversity of viruses in H. dromedarii ticks.

Benchmarking of Nanopore R10.4 and R9.4.1 flow cells in single-cell whole-genome amplification and whole-genome shotgun sequencing

Third-generation sequencing can be used in human cancer genomics and epigenomic research. Oxford Nanopore Technologies (ONT) recently released R10.4 flow cell, which claimed an improved read accuracy compared to R9.4.1 flow cell. To evaluate the benefits and defects of R10.4 flow cell for cancer cell profiling on MinION devices, we used the human non-small-cell lung-carcinoma cell line HCC78 to construct libraries for both single-cell whole-genome amplification (scWGA) and whole-genome shotgun sequencing. The R10.4 and R9.4.1 reads were benchmarked in terms of read accuracy, variant detection, modification calling, genome recovery rate and compared with the next generation sequencing (NGS) reads. The results highlighted that the R10.4 outperforms R9.4.1 reads, achieving a higher modal read accuracy of over 99.1%, superior variation detection, lower false-discovery rate (FDR) in methylation calling, and comparable genome recovery rate. To achieve high yields scWGA sequencing in the ONT platform as NGS, we recommended multiple displacement amplification with a modified T7 endonuclease Ⅰ cutting procedure as a promising method. In addition, we provided a possible solution to filter the likely false positive sites among the whole genome region with R10.4 by using scWGA sequencing result as a negative control. Our study is the first benchmark of whole genome single-cell sequencing using ONT R10.4 and R9.4.1 MinION flow cells by clarifying the capacity of genomic and epigenomic profiling within a single flow cell. A promising method for scWGA sequencing together with the methylation calling results can benefit researchers who work on cancer cell genomic and epigenomic profiling using third-generation sequencing.

Benchmarking of ATAC Sequencing Data From BGI's Low-Cost DNBSEQ-G400 Instrument for Identification of Open and Occupied Chromatin Regions.

Background: Chromatin falls into one of two major subtypes: closed heterochromatin and euchromatin which is accessible, transcriptionally active, and occupied by transcription factors (TFs). The most widely used approach to interrogate differences in the chromatin state landscape is the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). While library generation is relatively inexpensive, sequencing depth requirements can make this assay cost-prohibitive for some laboratories. Findings: Here, we benchmark data from Beijing Genomics Institute’s (BGI) DNBSEQ-G400 low-cost sequencer against data from a standard Illumina instrument (HiSeqX10). For comparisons, the same bulk ATAC-seq libraries generated from pluripotent stem cells (PSCs) and fibroblasts were sequenced on both platforms. Both instruments generate sequencing reads with comparable mapping rates and genomic context. However, DNBSEQ-G400 data contained a significantly higher number of small, sub-nucleosomal reads (>30% increase) and a reduced number of bi-nucleosomal reads (>75% decrease), which resulted in narrower peak bases and improved peak calling, enabling the identification of 4% more differentially accessible regions between PSCs and fibroblasts. The ability to identify master TFs that underpin the PSC state relative to fibroblasts (via HOMER, HINT-ATAC, TOBIAS), namely, foot-printing capacity, were highly similar between data generated on both platforms. Integrative analysis with transcriptional data equally enabled direct recovery of three published 3-factor combinations that have been shown to induce pluripotency. Conclusion: Other than a small increase in peak calling sensitivity for DNBSEQ-G400 data (BGI), both platforms enable comparable levels of open chromatin identification for ATAC-seq library sequencing, yielding similar analytical outcomes, albeit at low-data generation costs in the case of the BGI instrument.

Metagenomic analysis of viral community in the Yangtze River expands known eukaryotic and prokaryotic virus diversity in freshwater.

Viruses in aquatic ecosystems are characterized by extraordinary abundance and diversity. Thus far, there have been limited studies focused on viral communities in river water systems. Here, we investigated the virome of the Yangtze River Delta using viral metagenomic analysis. The compositions of viral communities from six sampling sites were analyzed and compared. By using library construction and next generation sequencing, contigs and singlet reads similar to viral sequences were classified into 17 viral families, including nine dsDNA viral families, four ssDNA viral families and four RNA viral families. Statistical analysis using Friedman test suggested that there was no significant difference among the six sampling sites (P ​> ​0.05). The viromes in this study were all dominated by the order Caudovirales, and a group of Freshwater phage uvFW species were particularly prevalent among all the samples. The virome from Nanjing presented a unique pattern of viral community composition with a relatively high abundance of family Parvoviridae. Phylogenetic analyses based on virus hallmark genes showed that the Caudovirales order and CRESS-DNA viruses presented high genetic diversity, while viruses in the Microviridae and Parvoviridae families and the Riboviria realm were relatively conservative. Our study provides the first insight into viral community composition in large river ecosystem, revealing the diversity and stability of river water virome, contributing to the proper utilization of freshwater resource.

Accurate profiling of forensic autosomal STRs using the Oxford Nanopore Technologies MinION device

The high variability characteristic of short tandem repeat (STR) markers is harnessed for human identification in forensic genetic analyses. Despite the power and reliability of current typing techniques, sequence-level information both within and around STRs are masked in the length-based profiles generated. Forensic STR typing using next generation sequencing (NGS) has therefore gained attention as an alternative to traditional capillary electrophoresis (CE) approaches. In this proof-of-principle study, we evaluate the forensic applicability of the newest and smallest NGS platform available – the Oxford Nanopore Technologies (ONT) MinION device. Although nanopore sequencing on the handheld MinION offers numerous advantages, including low startup cost and on-site sample processing, the relatively high error rate and lack of forensic-specific analysis software has prevented accurate profiling across STR panels in previous studies. Here we present STRspy, a streamlined method capable of producing length- and sequence-based STR allele designations from noisy, error-prone third generation sequencing reads. To assess the capabilities of STRspy, seven reference samples (female: n = 2; male: n = 5) were amplified at 15 and 30 PCR cycles using the Promega PowerSeq 46GY System and sequenced on the ONT MinION device in triplicate. Basecalled reads were then processed with STRspy using a custom database containing alleles reported in the STRSeq BioProject NIST 1036 dataset. Resultant STR allele designations and flanking region single nucleotide polymorphism (SNP) calls were compared to the manufacturer-validated genotypes for each sample. STRspy generated robust and reliable genotypes across all autosomal STR loci amplified with 30 PCR cycles, achieving 100% concordance based on both length and sequence. Furthermore, we were able to identify flanking region SNPs in the 15-cycle dataset with > 90% accuracy. These results demonstrate that when analyzed with STRspy ONT reads can reveal additional variation in and around STR loci depending on read coverage. As the first and only third generation sequencing platform-specific method to successfully profile the entire panel of autosomal STRs amplified by a commercially available multiplex, STRspy significantly increases the feasibility of nanopore sequencing in forensic applications.

SWALO: scaffolding with assembly likelihood optimization.

Scaffolding, i.e. ordering and orienting contigs is an important step in genome assembly. We present a method for scaffolding using second generation sequencing reads based on likelihoods of genome assemblies. A generative model for sequencing is used to obtain maximum likelihood estimates of gaps between contigs and to estimate whether linking contigs into scaffolds would lead to an increase in the likelihood of the assembly. We then link contigs if they can be unambiguously joined or if the corresponding increase in likelihood is substantially greater than that of other possible joins of those contigs. The method is implemented in a tool called Swalo with approximations to make it efficient and applicable to large datasets. Analysis on real and simulated datasets reveals that it consistently makes more or similar number of correct joins as other scaffolders while linking very few contigs incorrectly, thus outperforming other scaffolders and demonstrating that substantial improvement in genome assembly may be achieved through the use of statistical models. Swalo is freely available for download at https://atifrahman.github.io/SWALO/.

ZWA: Viral genome assembly and characterization hindrances from virus-host chimeric reads; a refining approach.

Viral metagenomics, also known as virome studies, have yielded an unprecedented number of novel sequences, essential in recognizing and characterizing the etiological agent and the origin of emerging infectious diseases. Several tools and pipelines have been developed, to date, for the identification and assembly of viral genomes. Assembly pipelines often result in viral genomes contaminated with host genetic material, some of which are currently deposited into public databases. In the current report, we present a group of deposited sequences that encompass ribosomal RNA (rRNA) contamination. We highlight the detrimental role of chimeric next generation sequencing reads, between host rRNA sequences and viral sequences, in virus genome assembly and we present the hindrances these reads may pose to current methodologies. We have further developed a refining pipeline, the Zero Waste Algorithm (ZWA) that assists in the assembly of low abundance viral genomes. ZWA performs context-depended trimming of chimeric reads, precisely removing their rRNA moiety. These, otherwise discarded, reads were fed to the assembly pipeline and assisted in the construction of larger and cleaner contigs making a substantial impact on current assembly methodologies. ZWA pipeline may significantly enhance virus genome assembly from low abundance samples and virus metagenomics approaches in which a small number of reads determine genome quality and integrity.

Clinical and Microbiological Characteristics of a Community-Acquired Carbapenem-Resistant Escherichia coli ST410 Isolate Harbouring blaNDM-5-Encoding IncX3-Type Plasmid From Blood.

Objectives: The aim of this research was to investigate the clinical and microbiological characteristics of a case of community-acquired carbapenem-resistant Escherichia coli isolated from a patient with a bloodstream infection in China.Methods: Escherichia coli Huamei202001 was recovered from the first blood culture from a patient hospitalised in China. An antimicrobial susceptibility test was performed, and the genome was sequenced on an Illumina HiSeq X 10 platform with a 150-bp paired-end approach. The generated sequence reads were assembled using Unicycler, and the whole genome sequence data were analysed using bioinformatics tools. Moreover, the patient and her main family members obtained a faecal sample screening test for CRE, the positive strain was further isolated and the identification and antimicrobial susceptibility testing was performed.Results: Escherichia coli Huamei202001 belonged to sequence type 410. In addition, a blaNDM-5-encoding IncX3-type plasmid was responsible for the spreading of carbapenem resistance. Only the patient was detected as having a positive faecal sample screening test for CRE. Strain Fec01 was identified as E. coli, and the antibiotic susceptibility profile was the same as that of E. coli Huamei202001.Conclusions: Escherichia coli Huamei202001 is defined as community-acquired carbapenem-resistant Enterobacteriaceae. The clone ST410 that harbours the blaNDM-5-encoding IncX3-type plasmid is causing new high-risk clones globally. Thus, infection control measures should be strengthened to curb the dissemination of IncX3.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes.

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

IsoDetect: Detection of Splice Isoforms from Third Generation Long Reads Based on Short Feature Sequences

Background: Transcriptome annotation is the basis for understanding gene structures and analysing gene expression. The transcriptome annotation of many organisms such as humans is far from incomplete, due partly to the challenge in the identification of isoforms that are produced from the same gene through alternative splicing. Third generation sequencing (TGS) reads provide unprecedented opportunity for detecting isoforms due to their long length that exceeds the length of most isoforms. One limitation of current TGS reads-based isoform detection methods is that they are exclusively based on sequence reads, without incorporating the sequence information of annotated isoforms. Objective: We aim to develop a method to detect isoforms by incorporating annotated isoforms. Methods: Based on annotated isoforms, we propose a splice isoform detection method called IsoDetect. First, the sequence at exon-exon junctions is extracted from annotated isoforms as “short feature sequences”, which is used to distinguish splice isoforms. Second, we align these feature sequences to long reads and partition long reads into groups that contain the same set of feature sequences, thereby avoiding the pair-wise comparison among the large number of long reads. Third, clustering and consensus generation are carried out based on sequence similarity. For the long reads that do not contain any short feature sequence, clustering analysis based on sequence similarity is performed to identify isoforms. Therefore, our method can detect not only known but also novel isoforms. Result: Tested on two datasets from Calypte anna and Zebra Finch, IsoDetect shows higher speed and good accuracies compared with four existing methods. Conclusion: IsoDetect may become a promising method for isoform detection.

A strategy for complete telomere-to-telomere assembly of ciliate macronuclear genome using ultra-high coverage Nanopore data

Ciliates contain two kinds of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC) in a single cell. The MAC usually have fragmented chromosomes. These fragmented chromosomes, capped with telomeres at both ends, could be gene size to several megabases in length among different ciliate species. So far, no telomere-to-telomere assembly of entire MAC genome in ciliate species has been finished. Development of the third generation sequencing technologies allows to generate sequencing reads up to megabases in length that could possibly span an entire MAC chromosome. Taking advantage of the ultra-long Nanopore reads, we established a simple strategy for the complete assembly of ciliate MAC genomes. Using this strategy, we assembled the complete MAC genomes of two ciliate species Tetrahymena thermophila and Tetrahymena shanghaiensis, composed of 181 and 214 chromosomes telomere-to-telomere respectively. The established strategy as well as the high-quality genome data will provide a useful approach for ciliate genome assembly, and a valuable community resource for further biological, evolutionary and population genomic studies.

Draft genome sequence of a multidrug-resistant Stenotrophomonas sp. B1-1 strain isolated from radiation-polluted soil and its pathogenic potential

ObjectivesStenotrophomonas is a genus of Gram-negative bacteria with several potential industrial uses as well as an increasingly relevant pathogen that may cause dangerous nosocomial infections. Here we present the draft genome sequence of a multidrug-resistant Stenotrophomonas sp. B1-1 isolated from radiation-polluted soil in Xinjiang Uyghur Autonomous Region, China. MethodsThe genome of Stenotrophomonas sp. B1-1 was sequenced using a BGISEQ-500 platform. The generated sequencing reads were de novo assembled using SOAPdenovo and the resulting sequences were predicted and annotated to identify antimicrobial resistance genes and virulence factors using the ARDB and VFDB databases, respectively. ResultsThe Stenotrophomonas sp. B1-1 genome assembly resulted in a total genome size of 4,723,769 bp with a GC content of 67.47%. There were 4280 predicted genes with 68 tRNAs, 2 rRNAs and 163 sRNAs. A number of antimicrobial resistance genes were identified conferring resistance to various antibiotics as well as numerous virulence genes. ConclusionThe genome sequence of Stenotrophomonas sp. B1-1 will provide timely information for comparison of the Stenotrophomonas genus and to help further understand the pathogenesis and antimicrobial resistance of this genus.

Complete Draft Genome Sequence of Cutibacterium (Propionibacterium) acnes Type Strain ATCC6919

Abstract Objectives: The wide use of antibiotics for acne is the main cause of drug-resistant Cutibacterium acnes. However, the complete genetic sequence reference of C. acnes has not yet been reported. Herein, we report the complete genetic sequence of Cutibacterium acnes type strain ATCC6919. Methods: Whole-genome sequencing of Cutibacterium acnes ATCC6919 was performed using both the Illumina Hiseq4000 and PacBio RSII single-molecule real-time sequencing platforms. The generated sequence reads were de novo assembled using the HGAP3 protocol of single-molecule real-time analysis v.2.3.0 software. The draft genome was annotated and further analyzed by various bioinformatics tools. Results: The total genome length of Cutibacterium acnes ATCC6919 was 2,495,001 bp (60.02% G + C content), with 2,450 coding sequences including 2,358 proteins, 45 tRNAs, 9 rRNAs, and 1 sRNA. The Comprehensive Antibiotic Resistance Database showed that the gene (thy A) found in ATCC6919 was associated with resistance to aminosalicylic acid. Antimicrobial susceptibility testing showed that this isolate was not sensitive to metronidazole and sulfamethoxazole, but was sensitive to all other tested antimicrobials (erythromycin, azithromycin, clarithromycin, tetracycline, minocycline, trimethoprim, levofloxacin, ceftriaxone, and doxycycline). Conclusion: This study reports the complete genome sequence of Cutibacterium acnes ATCC6919, which provides a valuable reference for comparative analysis based on the genomic characteristics and resistance gene diversity of Cutibacterium acnes.