Menu

Auto Tag: Read length

September 22, 2019

NextSV: a meta-caller for structural variants from low-coverage long-read sequencing data.

Structural variants (SVs) in human genomes are implicated in a variety of human diseases. Long-read sequencing delivers much longer read lengths than short-read sequencing and may greatly improve SV detection. However, due to the relatively high cost of long-read sequencing, it is unclear what coverage is needed and how to optimally use the aligners and SV callers.In this study, we developed NextSV, a meta-caller to perform SV calling from low coverage long-read sequencing data. NextSV integrates three aligners and three SV callers and generates two integrated call sets (sensitive/stringent) for different analysis purposes. We evaluated SV calling performance of NextSV under different PacBio coverages on two personal genomes, NA12878 and HX1. Our results showed that, compared with running any single SV caller, NextSV stringent call set had higher precision and balanced accuracy (F1 score) while NextSV sensitive call set had a higher recall. At 10X coverage, the recall of NextSV sensitive call set was 93.5 to 94.1% for deletions and 87.9 to 93.2% for insertions, indicating that ~10X coverage might be an optimal coverage to use in practice, considering the balance between the sequencing costs and the recall rates. We further evaluated the Mendelian errors on an Ashkenazi Jewish trio dataset.Our results provide useful guidelines for SV detection from low coverage whole-genome PacBio data and we expect that NextSV will facilitate the analysis of SVs on long-read sequencing data.

Read more
September 22, 2019

NPBSS: a new PacBio sequencing simulator for generating the continuous long reads with an empirical model.

PacBio sequencing platform offers longer read lengths than the second-generation sequencing technologies. It has revolutionized de novo genome assembly and enabled the automated reconstruction of reference-quality genomes. Due to its extremely wide range of application areas, fast sequencing simulation systems with high fidelity are in great demand to facilitate the development and comparison of subsequent analysis tools. Although there are several available simulators (e.g., PBSIM, SimLoRD and FASTQSim) that target the specific generation of PacBio libraries, the error rate of simulated sequences is not well matched to the quality value of raw PacBio datasets, especially for PacBio’s continuous long reads (CLR).By analyzing the characteristic features of CLR data from PacBio SMRT (single molecule real time) sequencing, we developed a new PacBio sequencing simulator (called NPBSS) for producing CLR reads. NPBSS simulator firstly samples the read sequences according to the read length logarithmic normal distribution, and choses different base quality values with different proportions. Then, NPBSS computes the overall error probability of each base in the read sequence with an empirical model, and calculates the deletion, substitution and insertion probabilities with the overall error probability to generate the PacBio CLR reads. Alignment results demonstrate that NPBSS fits the error rate of the PacBio CLR reads better than PBSIM and FASTQSim. In addition, the assembly results also show that simulated sequences of NPBSS are more like real PacBio CLR data.NPBSS simulator is convenient to use with efficient computation and flexible parameters setting. Its generating PacBio CLR reads are more like real PacBio datasets.

Read more
September 22, 2019

Genome-wide analysis of Mycoplasma bovirhinis GS01 reveals potential virulence factors and phylogenetic relationships.

Mycoplasma bovirhinis is a significant etiology in bovine pneumonia and mastitis, but our knowledge about the genetic and pathogenic mechanisms of M. bovirhinis is very limited. In this study, we sequenced the complete genome of M. bovirhinis strain GS01 isolated from the nasal swab of pneumonic calves in Gansu, China, and we found that its genome forms a 847,985 bp single circular chromosome with a GC content of 27.57% and with 707 protein-coding genes. The putative virulence determinants of M. bovirhinis were then analyzed. Results showed that three genomic islands and 16 putative virulence genes, including one adhesion gene enolase, seven surface lipoproteins, proteins involved in glycerol metabolism, and cation transporters, might be potential virulence factors. Glycerol and pyruvate metabolic pathways were defective. Comparative analysis revealed remarkable genome variations between GS01 and a recently reported HAZ141_2 strain, and extremely low homology with others mycoplasma species. Phylogenetic analysis demonstrated that M. bovirhinis was most genetically close to M. canis, distant from other bovine Mycoplasma species. Genomic dissection may provide useful information on the pathogenic mechanisms and genetics of M. bovirhinis. Copyright © 2018 Chen et al.

Read more
September 22, 2019

A whole genome assembly of the horn fly, Haematobia irritans, and prediction of genes with roles in metabolism and sex determination.

Haematobia irritans, commonly known as the horn fly, is a globally distributed blood-feeding pest of cattle that is responsible for significant economic losses to cattle producers. Chemical insecticides are the primary means for controlling this pest but problems with insecticide resistance have become common in the horn fly. To provide a foundation for identification of genomic loci for insecticide resistance and for discovery of new control technology, we report the sequencing, assembly, and annotation of the horn fly genome. The assembled genome is 1.14 Gb, comprising 76,616 scaffolds with N50 scaffold length of 23 Kb. Using RNA-Seq data, we have predicted 34,413 gene models of which 19,185 have been assigned functional annotations. Comparative genomics analysis with the Dipteran flies Musca domestica L., Drosophila melanogaster, and Lucilia cuprina, show that the horn fly is most closely related to M. domestica, sharing 8,748 orthologous clusters followed by D. melanogaster and L. cuprina, sharing 7,582 and 7,490 orthologous clusters respectively. We also identified a gene locus for the sodium channel protein in which mutations have been previously reported that confers target site resistance to the most common class of pesticides used in fly control. Additionally, we identified 276 genomic loci encoding members of metabolic enzyme gene families such as cytochrome P450s, esterases and glutathione S-transferases, and several genes orthologous to sex determination pathway genes in other Dipteran species. Copyright © 2018 Konganti et al.

Read more
September 22, 2019

Genome sequencing and comparative analysis of Stenotrophomonas acidaminiphila reveal evolutionary insights into sulfamethoxazole resistance.

Stenotrophomonas acidaminiphila is an aerobic, glucose non-fermentative, Gram-negative bacterium that been isolated from various environmental sources, particularly aquatic ecosystems. Although resistance to multiple antimicrobial agents has been reported in S. acidaminiphila, the mechanisms are largely unknown. Here, for the first time, we report the complete genome and antimicrobial resistome analysis of a clinical isolate S. acidaminiphila SUNEO which is resistant to sulfamethoxazole. Comparative analysis among closely related strains identified common and strain-specific genes. In particular, comparison with a sulfamethoxazole-sensitive strain identified a mutation within the sulfonamide-binding site of folP in SUNEO, which may reduce the binding affinity of sulfamethoxazole. Selection pressure analysis indicated folP in SUNEO is under purifying selection, which may be owing to long-term administration of sulfonamide against Stenotrophomonas.

Read more
September 22, 2019

The complete chloroplast genome sequence of Actinidia arguta using the PacBio RS II platform.

Actinidia arguta is the most basal species in a phylogenetically and economically important genus in the family Actinidiaceae. To better understand the molecular basis of the Actinidia arguta chloroplast (cp), we sequenced the complete cp genome from A. arguta using Illumina and PacBio RS II sequencing technologies. The cp genome from A. arguta was 157,611 bp in length and composed of a pair of 24,232 bp inverted repeats (IRs) separated by a 20,463 bp small single copy region (SSC) and an 88,684 bp large single copy region (LSC). Overall, the cp genome contained 113 unique genes. The cp genomes from A. arguta and three other Actinidia species from GenBank were subjected to a comparative analysis. Indel mutation events and high frequencies of base substitution were identified, and the accD and ycf2 genes showed a high degree of variation within Actinidia. Forty-seven simple sequence repeats (SSRs) and 155 repetitive structures were identified, further demonstrating the rapid evolution in Actinidia. The cp genome analysis and the identification of variable loci provide vital information for understanding the evolution and function of the chloroplast and for characterizing Actinidia population genetics.

Read more
September 22, 2019

Whole genome sequence and comparative analysis of Borrelia burgdorferi MM1.

Lyme disease is caused by spirochaetes of the Borrelia burgdorferi sensu lato genospecies. Complete genome assemblies are available for fewer than ten strains of Borrelia burgdorferi sensu stricto, the primary cause of Lyme disease in North America. MM1 is a sensu stricto strain originally isolated in the midwestern United States. Aside from a small number of genes, the complete genome sequence of this strain has not been reported. Here we present the complete genome sequence of MM1 in relation to other sensu stricto strains and in terms of its Multi Locus Sequence Typing. Our results indicate that MM1 is a new sequence type which contains a conserved main chromosome and 15 plasmids. Our results include the first contiguous 28.5 kb assembly of lp28-8, a linear plasmid carrying the vls antigenic variation system, from a Borrelia burgdorferi sensu stricto strain.

Read more
September 22, 2019

Insect symbionts as valuable grist for the biotechnological mill: an alkaliphilic silkworm gut bacterium for efficient lactic acid production.

Insects constitute the most abundant and diverse animal class and act as hosts to an extraordinary variety of symbiotic microorganisms. These microbes living inside the insects play critical roles in host biology and are also valuable bioresources. Enterococcus mundtii EMB156, isolated from the larval gut (gut pH >10) of the model organism Bombyx mori (Lepidoptera: Bombycidae), efficiently produces lactic acid, an important metabolite for industrial production of bioplastic materials. E. mundtii EMB156 grows well under alkaline conditions and stably converts various carbon sources into lactic acid, offering advantages in downstream fermentative processes. High-yield lactic acid production can be achieved by the strain EMB156 from renewable biomass substrates under alkaline pretreatments. Single-molecule real-time (SMRT) sequencing technology revealed its 3.01 Mbp whole genome sequence. A total of 2956 protein-coding sequences, 65 tRNA genes, and 6 rRNA operons were predicted in the EMB156 chromosome. Remarkable genomic features responsible for lactic acid fermentation included key enzymes involved in the pentose phosphate (PP)/glycolytic pathway, and an alpha amylase and xylose isomerase were characterized in EMB156. This genomic information coincides with the phenotype of E. mundtii EMB156, reflecting its metabolic flexibility in efficient lactate fermentation, and established a foundation for future biotechnological application. Interestingly, enzyme activities of amylase were quite stable in high-pH broths, indicating a possible mechanism for strong EMB156 growth in an alkaline environment, thereby facilitating lactic acid production. Together, these findings implied that valuable lactic acid-producing bacteria can be discovered efficiently by screening under the extremely alkaline conditions, as exemplified by gut microbial symbionts of Lepidoptera insects.

Read more
September 22, 2019

The genome of Artemisia annua provides insight into the evolution of Asteraceae family and artemisinin biosynthesis.

Artemisia annua, commonly known as sweet wormwood or Qinghao, is a shrub native to China and has long been used for medicinal purposes. A. annua is now cultivated globally as the only natural source of a potent anti-malarial compound, artemisinin. Here, we report a high-quality draft assembly of the 1.74-gigabase genome of A. annua, which is highly heterozygous, rich in repetitive sequences, and contains 63 226 protein-coding genes, one of the largest numbers among the sequenced plant species. We found that, as one of a few sequenced genomes in the Asteraceae, the A. annua genome contains a large number of genes specific to this large angiosperm clade. Notably, the expansion and functional diversification of genes encoding enzymes involved in terpene biosynthesis are consistent with the evolution of the artemisinin biosynthetic pathway. We further revealed by transcriptome profiling that A. annua has evolved the sophisticated transcriptional regulatory networks underlying artemisinin biosynthesis. Based on comprehensive genomic and transcriptomic analyses we generated transgenic A. annua lines producing high levels of artemisinin, which are now ready for large-scale production and thereby will help meet the challenge of increasing global demand of artemisinin. Copyright © 2018 The Author. Published by Elsevier Inc. All rights reserved.

Read more
September 22, 2019

Mycobacterial biomaterials and resources for researchers.

There are many resources available to mycobacterial researchers, including culture collections around the world that distribute biomaterials to the general scientific community, genomic and clinical databases, and powerful bioinformatics tools. However, many of these resources may be unknown to the research community. This review article aims to summarize and publicize many of these resources, thus strengthening the quality and reproducibility of mycobacterial research by providing the scientific community access to authenticated and quality-controlled biomaterials and a wealth of information, analytical tools and research opportunities.

Read more
September 22, 2019

Parallels between experimental and natural evolution of legume symbionts.

The emergence of symbiotic interactions has been studied using population genomics in nature and experimental evolution in the laboratory, but the parallels between these processes remain unknown. Here we compare the emergence of rhizobia after the horizontal transfer of a symbiotic plasmid in natural populations of Cupriavidus taiwanensis, over 10 MY ago, with the experimental evolution of symbiotic Ralstonia solanacearum for a few hundred generations. In spite of major differences in terms of time span, environment, genetic background, and phenotypic achievement, both processes resulted in rapid genetic diversification dominated by purifying selection. We observe no adaptation in the plasmid carrying the genes responsible for the ecological transition. Instead, adaptation was associated with positive selection in a set of genes that led to the co-option of the same quorum-sensing system in both processes. Our results provide evidence for similarities in experimental and natural evolutionary transitions and highlight the potential of comparisons between both processes to understand symbiogenesis.

Read more
September 22, 2019

The impact of Staphylococcus aureus genomic variation on clinical phenotype of children with acute hematogenous osteomyelitis.

Children with acute hematogenous osteomyelitis (AHO) have a broad spectrum of illness ranging from mild to severe. The purpose of this study is to evaluate the impact of genomic variation of Staphylococcus aureus on clinical phenotype of affected children and determine which virulence genes correlate with severity of illness.De novo whole genome sequencing was conducted for a strain of Community Acquired Methicillin Resistant Staphylococcus aureus (CA-MRSA), using PacBio Hierarchical Genome Assembly Process (HGAP) from 6 Single Molecule Real Time (SMRT) Cells, as a reference for DNA library assembly of 71 Staphylococcus aureus isolates from children with AHO. Virulence gene annotation was based on exhaustive literature review and genomic data in NCBI for Staphylococcus aureus. Clinical phenotype was assessed using a validated severity score. Kruskal-Wallis rank sum test determined association between clinical severity and virulence gene presence using False Discovery Rate (FDR), significance <0.01.PacBio produced an assembled genome of 2,898,306 bp and 2054 Open Reading Frames (ORFs). Annotation confirmed 201 virulence genes. Statistical analysis of gene presence by clinical severity found 40 genes significantly associated with severity of illness (FDR =0.009). MRSA isolates encoded a significantly greater number of virulence genes than did MSSA (p < 0.0001). Phylogenetic analysis by maximum likelihood (PAML) demonstrated the relatedness of genomic distance to clinical phenotype.The Staphylococcus aureus genome contains virulence genes which are significantly associated with severity of illness in children with osteomyelitis. This study introduces a novel reference strain and detailed annotation of Staphylococcus aureus virulence genes. While this study does not address bacterial gene expression, a platform is created for future transcriptome investigations to elucidate the complex mechanisms involved in childhood osteomyelitis.

Read more
September 22, 2019

Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris.

A parasitic lifestyle, where plants procure some or all of their nutrients from other living plants, has evolved independently in many dicotyledonous plant families and is a major threat for agriculture globally. Nevertheless, no genome sequence of a parasitic plant has been reported to date. Here we describe the genome sequence of the parasitic field dodder, Cuscuta campestris. The genome contains signatures of a fairly recent whole-genome duplication and lacks genes for pathways superfluous to a parasitic lifestyle. Specifically, genes needed for high photosynthetic activity are lost, explaining the low photosynthesis rates displayed by the parasite. Moreover, several genes involved in nutrient uptake processes from the soil are lost. On the other hand, evidence for horizontal gene transfer by way of genomic DNA integration from the parasite’s hosts is found. We conclude that the parasitic lifestyle has left characteristic footprints in the C. campestris genome.

Read more
September 22, 2019

Landscape of the genome and host cell response of Mycobacterium shigaense reveals pathogenic features.

A systems approach was used to explore the genome and transcriptome of Mycobacterium shigaense, a new opportunistic pathogen isolated from a patient with a skin infection, and the host response transcriptome was assessed using a macrophage infection model. The M. shigaense genome comprises 5,207,883?bp, with 67.2% G+C content and 5098 predicted coding genes. Evolutionarily, the bacterium belongs to a cluster in the phylogenetic tree along with three target opportunistic pathogenic strains, namely, M. avium, M. triplex and M. simiae. Potential virulence genes are indeed expressed by M. shigaense under culture conditions. Phenotypically, M. shigaense had similar infection and replication capacities in a macrophage model as the opportunistic species compared to M. tuberculosis. M. shigaense activated NF-?B, TNF, cytokines and chemokines in the host innate immune-related signaling pathways and elicited an early response shared with pathogenic bacilli except M. tuberculosis. M. shigaense upregulated specific host response genes such as TLR7, CCL4 and CXCL5. We performed an integrated and comparative analysis of M. shigaense. Multigroup comparison indicated certain differences with typical pathogenic bacilli in terms of gene features and the macrophage response.

Read more
September 22, 2019

Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation.

Echinoderms exhibit several fascinating evolutionary innovations that are rarely seen in the animal kingdom, but how these animals attained such features is not well understood. Here we report the sequencing and analysis of the genome and extensive transcriptomes of the sea cucumber Apostichopus japonicus, a species from a special echinoderm group with extraordinary potential for saponin synthesis, aestivation and organ regeneration. The sea cucumber does not possess a reorganized Hox cluster as previously assumed for all echinoderms, and the spatial expression of Hox7 and Hox11/13b potentially guides the embryo-to-larva axial transformation. Contrary to the typical production of lanosterol in animal cholesterol synthesis, the oxidosqualene cyclase of sea cucumber produces parkeol for saponin synthesis and has “plant-like” motifs suggestive of convergent evolution. The transcriptional factors Klf2 and Egr1 are identified as key regulators of aestivation, probably exerting their effects through a clock gene-controlled process. Intestinal hypometabolism during aestivation is driven by the DNA hypermethylation of various metabolic gene pathways, whereas the transcriptional network of intestine regeneration involves diverse signaling pathways, including Wnt, Hippo and FGF. Decoding the sea cucumber genome provides a new avenue for an in-depth understanding of the extraordinary features of sea cucumbers and other echinoderms.

Read more

Subscribe for blog updates:

Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.