Oxford Nanopore MinION Sequencing and Genome Assembly

The revolution of genome sequencing is continuing after the successful second-generation sequencing (SGS) technology. The third-generation sequencing (TGS) technology, led by Pacific Biosciences (PacBio), is progressing rapidly, moving from a technology once only capable of providing data for small genome analysis, or for performing targeted screening, to one that pro-mises high quality de novo assembly and structural variation detection for human-sized genomes. In 2014, the MinION, the first commercial sequencer using nanopore technology, was released by Oxford Nanopore Technologies (ONT). MinION identifies DNA bases by measuring the changes in electrical conductivity generated as DNA strands pass through a biological pore. Its portability, affordability, and speed in data production makes it suitable for real-time applications, the release of the long read sequencer MinION has thus generated much excitement and interest in the geno-mics community. While de novo genome assemblies can be cheaply produced from SGS data, assem-bly continuity is often relatively poor, due to the limited ability of short reads to handle long repeats. Assembly quality can be greatly improved by using TGS long reads, since repetitive regions can be easily expanded into using longer sequencing lengths, despite having higher error rates at the base level. The potential of nanopore sequencing has been demonstrated by various studies in gen-ome surveillance at locations where rapid and reliable sequencing is needed, but where resources are limited.     

Whole Genome Amplification and De novo Assembly of Single Bacterial Cells

Background

Single-cell genome sequencing has the potential to allow the in-depth exploration of the vast genetic diversity found in uncultured microbes. We used the marine cyanobacterium Prochlorococcus as a model system for addressing important challenges facing high-throughput whole genome amplification (WGA) and complete genome sequencing of individual cells.

Methodology/Principal Findings

We describe a pipeline that enables single-cell WGA on hundreds of cells at a time while virtually eliminating non-target DNA from the reactions. We further developed a post-amplification normalization procedure that mitigates extreme variations in sequencing coverage associated with multiple displacement amplification (MDA), and demonstrated that the procedure increased sequencing efficiency and facilitated genome assembly. We report genome recovery as high as 99.6% with reference-guided assembly, and 95% with de novo assembly starting from a single cell. We also analyzed the impact of chimera formation during MDA on de novo assembly, and discuss strategies to minimize the presence of incorrectly joined regions in contigs.

Conclusions/Significance

The methods describe in this paper will be useful for sequencing genomes of individual cells from a variety of samples.     

真菌基因组de novo测序组装的方法与实践

真菌基因组较其他真核生物基因组结构简单,长度短,易于测序、组装与注释,因此真菌基因组是研究真核生物基因组的模型。为研究真菌基因组组装策略,本研究基于Illumina HiSeq测序平台对烟曲霉菌株An16007基因组测序,分别使用5种de novo组装软件ABySS、SOAP-denovo、Velvet、MaSuRCA和IDBA-UD组装基因组,然后通过Augustus软件进行基因预测,BUSCO软件评估组装结果。研究发现,5种组装软件对基因组组装结果不同,ABySS组装的基因组较其他4种组装软件具有更高的完整性和准确性,且预测的基因数量较高,因此,ABySS更适合本研究基因组的组装。本研究提供了真菌de novo测序、组装及组装质量评估的技术流程,为基因组<100 Mb的真菌或其他生物基因组的研究提供参考。     

Identification of Optimum Sequencing Depth Especially for De Novo Genome Assembly of Small Genomes Using Next Generation Sequencing Data

Next Generation Sequencing (NGS) is a disruptive technology that has found widespread acceptance in the life sciences research community. The high throughput and low cost of sequencing has encouraged researchers to undertake ambitious genomic projects, especially in de novo genome sequencing. Currently, NGS systems generate sequence data as short reads and de novo genome assembly using these short reads is computationally very intensive. Due to lower cost of sequencing and higher throughput, NGS systems now provide the ability to sequence genomes at high depth. However, currently no report is available highlighting the impact of high sequence depth on genome assembly using real data sets and multiple assembly algorithms. Recently, some studies have evaluated the impact of sequence coverage, error rate and average read length on genome assembly using multiple assembly algorithms, however, these evaluations were performed using simulated datasets. One limitation of using simulated datasets is that variables such as error rates, read length and coverage which are known to impact genome assembly are carefully controlled. Hence, this study was undertaken to identify the minimum depth of sequencing required for de novo assembly for different sized genomes using graph based assembly algorithms and real datasets. Illumina reads for E.coli (4.6 MB) S.kudriavzevii (11.18 MB) and C.elegans (100 MB) were assembled using SOAPdenovo, Velvet, ABySS, Meraculous and IDBA-UD. Our analysis shows that 50X is the optimum read depth for assembling these genomes using all assemblers except Meraculous which requires 100X read depth. Moreover, our analysis shows that de novo assembly from 50X read data requires only 6–40 GB RAM depending on the genome size and assembly algorithm used. We believe that this information can be extremely valuable for researchers in designing experiments and multiplexing which will enable optimum utilization of sequencing as well as analysis resources.     

De novo Genome Assembly of the Fungal Plant Pathogen Pyrenophora semeniperda

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.     

Whole Genome Sequencing Reveals a De Novo SHANK3 Mutation in Familial Autism Spectrum Disorder

IntroductionClinical genomics promise to be especially suitable for the study of etiologically heterogeneous conditions such as Autism Spectrum Disorder (ASD). Here we present three siblings with ASD where we evaluated the usefulness of Whole Genome Sequencing (WGS) for the diagnostic approach to ASD.MethodsWe identified a family segregating ASD in three siblings with an unidentified cause. We performed WGS in the three probands and used a state-of-the-art comprehensive bioinformatic analysis pipeline and prioritized the identified variants located in genes likely to be related to ASD. We validated the finding by Sanger sequencing in the probands and their parents.ResultsThree male siblings presented a syndrome characterized by severe intellectual disability, absence of language, autism spectrum symptoms and epilepsy with negative family history for mental retardation, language disorders, ASD or other psychiatric disorders. We found germline mosaicism for a heterozygous deletion of a cytosine in the exon 21 of the SHANK3 gene, resulting in a missense sequence of 5 codons followed by a premature stop codon ({"type":"entrez-nucleotide","attrs":{"text":"NM_033517","term_id":"380748962","term_text":"NM_033517"}}NM_033517:c.3259_3259delC, p.Ser1088Profs*6).ConclusionsWe reported an infrequent form of familial ASD where WGS proved useful in the clinic. We identified a mutation in SHANK3 that underscores its relevance in Autism Spectrum Disorder.     

A Probabilistic Model to Predict Clinical Phenotypic Traits from Genome Sequencing

Genetic screening is becoming possible on an unprecedented scale. However, its utility remains controversial. Although most variant genotypes cannot be easily interpreted, many individuals nevertheless attempt to interpret their genetic information. Initiatives such as the Personal Genome Project (PGP) and Illumina''s Understand Your Genome are sequencing thousands of adults, collecting phenotypic information and developing computational pipelines to identify the most important variant genotypes harbored by each individual. These pipelines consider database and allele frequency annotations and bioinformatics classifications. We propose that the next step will be to integrate these different sources of information to estimate the probability that a given individual has specific phenotypes of clinical interest. To this end, we have designed a Bayesian probabilistic model to predict the probability of dichotomous phenotypes. When applied to a cohort from PGP, predictions of Gilbert syndrome, Graves'' disease, non-Hodgkin lymphoma, and various blood groups were accurate, as individuals manifesting the phenotype in question exhibited the highest, or among the highest, predicted probabilities. Thirty-eight PGP phenotypes (26%) were predicted with area-under-the-ROC curve (AUC)>0.7, and 23 (15.8%) of these were statistically significant, based on permutation tests. Moreover, in a Critical Assessment of Genome Interpretation (CAGI) blinded prediction experiment, the models were used to match 77 PGP genomes to phenotypic profiles, generating the most accurate prediction of 16 submissions, according to an independent assessor. Although the models are currently insufficiently accurate for diagnostic utility, we expect their performance to improve with growth of publicly available genomics data and model refinement by domain experts.     

De novo construction of cell-to-cell channels

Summary Nexus (gap junctions), which are considered to contain cell-to-cell channels, are newly formed in uterine smooth muscle during parturition or in response to estrogen treatment of virginal animals. A mRNA preparation was isolated from estrogen-dominated rat myometria and was encapsulated into liposomes. Subsequently the liposomes were fused with cultured cells of a mouse cell line CL-1D. It is established that these tumor cells normally are neither electrically coupled nor do they contain nexus. The cells, however, become electrically coupled a few hours after being loaded with the mRNA preparation. This de novo expression of cell coupling persisted for a little more than 24 hr after a single loading procedure. Freeze-fracture electron microscopy revealed small nexus-like particle aggregates at the time coupling was present. In control experiments the cells remained noncoupling when the RNA preparation was pretreated with ribonuclease, when cycloheximide was applied to the cells, or when liposomes filled with buffer solution only were used. These data suggest that the de novo expression of cell-to-cell coupling is accomplished by mRNA-induced protein biosynthesis resulting in the formation of cell-to-cell channels. Presented in the symposium on Molecular Morphological Aspects of Cell-Cell Communication at the 31 st Annual Meeting of the Tissue Culture Association, St. Louis, Missouri, June 1–5, 1980. This symposium was supported in part by Contract 263-MD-025754 from the National Cancer Institute and the Fogarty International Center.     

Dna Sequencing After the Human Genome Project

Abstract

In future DNA sequencing, gel electrophoresis, which is particularly effective for de novo sequencing, is likely to be replaced by sequencing by hybridization, mass spectrometry, or combinations of these two methods, which are particularly effective for comparative or diagnostic sequencing.     

Complete Genome Sequence of Amycolatopsis mediterranei S699 Based on De Novo Assembly via a Combinatorial Sequencing Strategy

The genome of Amycolatopsis mediterranei S699 was resequenced and assembled de novo. By comparing the sequences of S699 previously released and that of A. mediterranei U32, about 10 kb of major indels was found to differ between the two S699 genomes, and the differences are likely attributable to their different assembly strategies.     

De novo genesis of enhancers in vertebrates

Evolutionary innovation relies partially on changes in gene regulation. While a growing body of evidence demonstrates that such innovation is generated by functional changes or translocation of regulatory elements via mobile genetic elements, the de novo generation of enhancers from non-regulatory/non-mobile sequences has, to our knowledge, not previously been demonstrated. Here we show evidence for the de novo genesis of enhancers in vertebrates. For this, we took advantage of the massive gene loss following the last whole genome duplication in teleosts to systematically identify regions that have lost their coding capacity but retain sequence conservation with mammals. We found that these regions show enhancer activity while the orthologous coding regions have no regulatory activity. These results demonstrate that these enhancers have been de novo generated in fish. By revealing that minor changes in non-regulatory sequences are sufficient to generate new enhancers, our study highlights an important playground for creating new regulatory variability and evolutionary innovation.     

植物病原生物基因组测序进展

随着基因组学的迅猛发展,越来越多的生物全基因组序列已经测定完成或正在进行之中。植物病原生物基因组序列的测定为理解植物与病原物互作分子机制有重要意义,并为植物病理学的发展作出了重要贡献。目前已有9种病原细菌和1种病原真菌的基因组序列彻底完成,另外还有更多的基因组草图正在组装或测序工作正在进行之中。对NCBI上主要的植物病原真菌和细菌全基因组测序进展作了整理和概述。     

De novo malignant solitary fibrous tumor of the kidney

Background

To study the expression of MK-1 and RegIV and to detect their pathological significances in benign and malignant lesions of gallbladder.

Methods

The expression of MK-1 and RegIV was detected by immunohistochemical method in paraffin-embedded sections of surgical resected specimens from gallbladder adenocarcinoma (n = 108), peritumoral tissues (n = 46), adenomatous polyp (n = 15), and chronic cholecystitis (n = 35).

Results

The positive rate of MK-1 or RegIV expression was significantly higher in gallbladder adenocarcinoma than that in peritumoral tissues (X² _MK-1 = 18.76, P < 0.01; X² _RegIV = 9.92, P < 0.01), denomatous polyp (X² _MK-1 = 9.49, P < 0.01; X² _RegIV = 8.59, P < 0.01) and chronic cholecystitis (X² _MK-1 = 24.11, P < 0.01; X² _RegIV = 19.24, P < 0.01). The positive cases of MK-1 and/or RegIV in the benign lesions showed moderately- or severe-atypical hyperplasia of gallbladder epitheli. The positive rates of MK-1 were significantly higher in the cases of well-differentiated adenocarcinoma, no-metastasis of lymph node, and no-invasiveness of regional tissues than those in the ones of differentiated adenocarcinoma, metastasis of lymph node, and invasiveness of regional tissues in gallbladder adenocarcinoma (P < 0.05 or P < 0.01). On the contrary, the positive rates of RegIV were significantly lower in the cases of well-differentiated adenocarcinoma, no-metastasis of lymph node, and no-invasiveness of regional tissues than those in the ones of differentiated adenocarcinoma, metastasis of lymph node, and invasiveness of regional tissues in gallbladder adenocarcinoma (P < 0.05 or P < 0.01). Univariate Kaplan-Meier analysis showed that decreased expression of MK-1 (P = 0.09) or increased expression of RegIV (P = 0.003) was associated with decreased overall survival. Multivariate Cox regression analysis showed that decreased expression of MK-1 (P = 0.033) and increased expression of RegIV (P = 0.008) was an independent prognostic predictor in gallbladder adenocarcinoma.

Conclusions

The expression of MK-1 and/or RegIV might be closely related to the carcinogenesis, clinical biological behaviors, and prognosis of gallbladder adenocarcinoma.     

Personalized Oncogenomics: Clinical Experience with Malignant Peritoneal Mesothelioma Using Whole Genome Sequencing

Peritoneal mesothelioma is a rare and sometimes lethal malignancy that presents a clinical challenge for both diagnosis and management. Recent studies have led to a better understanding of the molecular biology of peritoneal mesothelioma. Translation of the emerging data into better treatments and outcome is needed. From two patients with peritoneal mesothelioma, we derived whole genome sequences, RNA expression profiles, and targeted deep sequencing data. Molecular data were made available for translation into a clinical treatment plan. Treatment responses and outcomes were later examined in the context of molecular findings. Molecular studies presented here provide the first reported whole genome sequences of peritoneal mesothelioma. Mutations in known mesothelioma-related genes NF2, CDKN2A, LATS2, amongst others, were identified. Activation of MET-related signaling pathways was demonstrated in both cases. A hypermutated phenotype was observed in one case (434 vs. 18 single nucleotide variants) and was associated with a favourable outcome despite sarcomatoid histology and multifocal disease. This study represents the first report of whole genome analyses of peritoneal mesothelioma, a key step in the understanding and treatment of this disease.