1,383 research outputs found

    Methods for Transcriptome Assembly in the Allopolyploid Brassica napus

    Get PDF
    Canada is the world’s largest producer of canola and the trend of production is ever increasing with an annual growth rate of 9.38% according to FAOSTAT. In 2017, canola acreage surpassed wheat in Saskatchewan, the highest producer of both crops in Canada. Country-wide, the total farming area of canola increased by 9.9% to 22.4 million acres while wheat area saw a slight decline to 23.3 million acres. While Canada is the highest producer of the crop, yields are lower than other countries. To maximize the benefit of this market, canola cultivation could be made more efficient with further characterization of the organism’s genes and their involvement in plant robustness. Such studies using transcriptome analysis have been successful in organisms with relatively small and simple genomes. However, such analyses in B. napus are complicated by the allopolyploid genome structure resulting from ancestral whole genome duplications in the species’ evolutionary history. Homeologous gene pairs originating from the orthology between the two B. napus progenitor species complicate the process of transcriptome assembly. Modern assemblers: Trinity, Oases and SOAPdenovo-Trans were used to generate several de novo transcriptome assemblies for B. napus. A variety of metrics were used to determine the impact that the complex genome structure has on transcriptome studies. In particular, the most important questions for transcriptome assembly in B. napus were how does varying the k-mer parameter effect assembly quality, and to what extent do similar genes resulting from homeology within B. napus complicate the process of assembly. These metrics used for evaluating the assemblies include basic assembly statistics such as the number of contigs and contig lengths (via N25, N50 and N75 statistics); as well as more involved investigation via comparison to annotated coding DNA sequences; evaluation softwares scores for de novo transcriptome assemblies and finally; quantification of homeolog differentiation by alignment to previously identified pairs of homeologous genes. These metrics provided a picture of the trade-offs between assembly softwares and the k-parameter determining the length of subsequences used to build de Bruijn graphs for de novo transcriptome assembly. It was shown that shorter k-mer lengths produce fewer, and more complete contigs due to the shorter required overlap between read sequences; while longer k-mer lengths increase the sensitivity of an assembler to sequence variation between similar gene sequences. The Trinity assembler outperformed Oases and SOAPdenovo-Trans when considering the total breadth of evaluation metrics, generating longer transcripts with fewer chimers between homeologous gene pairs

    Efficient Algorithms for Prokaryotic Whole Genome Assembly and Finishing

    Get PDF
    De-novo genome assembly from DNA fragments is primarily based on sequence overlap information. In addition, mate-pair reads or paired-end reads provide linking information for joining gaps and bridging repeat regions. Genome assemblers in general assemble long contiguous sequences (contigs) using both overlapping reads and linked reads until the assembly runs into an ambiguous repeat region. These contigs are further bridged into scaffolds using linked read information. However, errors can be made in both phases of assembly due to high error threshold of overlap acceptance and linking based on too few mate reads. Identical as well as similar repeat regions can often cause errors in overlap and mate-pair evidence. In addition, the problem of setting the correct threshold to minimize errors and optimize assembly of reads is not trivial and often requires a time-consuming trial and error process to obtain optimal results. The typical trial-and-error with multiple assembler, which can be computationally intensive, and is very inefficient, especially when users must learn how to use a wide variety of assemblers, many of which may be serial requiring long execution time and will not return usable or accurate results. Further, we show that the comparison of assembly results may not provide the users with a clear winner under all circumstances. Therefore, we propose a novel scaffolding tool, Correlative Algorithm for Repeat Placement (CARP), capable of joining short low error contigs using mate pair reads, computationally resolved repeat structures and synteny with one or more reference organisms. The CARP tool requires a set of repeat sequences such as insertion sequences (IS) that can be found computationally found without assembling the genome. Development of methods to identify such repeating regions directly from raw sequence reads or draft genomes led to the development of the ISQuest software package. ISQuest identifies bacterial ISs and their sequence elements—inverted and direct repeats—in raw read data or contigs using flexible search parameters. ISQuest is capable of finding ISs in hundreds of partially assembled genomes within hours; making it a valuable high-throughput tool for a global search of IS and repeat elements. The CARP tool matches very low error contigs with strong overlap using the ambiguous partial repeat sequence at the ends of the contig annotated using the repeat sequences discovered using ISQuest. These matches are verified by synteny with genomes of one or more reference organisms. We show that the CARP tool can be used to verify low mate pair evidence regions, independently find new joins and significantly reduce the number of scaffolds. Finally, we are demonstrate a novel viewer that presents to the user the computationally derived joins along with the evidence used to make the joins. The viewer allows the user to independently assess their confidence in the joins made by the finishing tools and make an informed decision of whether to invest the resources necessary to confirm a particular portion of the assembly. Further, we allow users to manually record join evidence, re-order contigs, and track the assembly finishing process

    Transcriptomic Profiling Using Next Generation Sequencing - Advances, Advantages, and Challenges

    Get PDF
    Transcriptome, the functional element of the genome, is comprised of different kinds of RNA molecules such as mRNA, miRNA, ncRNA, rRNA, and tRNA to name a few. Each of these RNA molecules plays a vital role in the physiological response, and understanding the regulation of these molecules is extremely critical for the better understanding of the functional genome. RNA Sequencing (RNASeq) is one of the latest techniques applied to study genome-wide transcriptome characterization and profiling using high-throughput sequenced data. As compared to array-based methods, RNASeq provides in-depth and more precise information on transcriptome characterization and quantification. Based upon availability of reference genome, transcriptome assembly can be reference-guided or de novo. Once transcripts are assembled, downstream analysis such as expression profiling, gene ontology, and pathway enrichment analyses can give more insight into gene regulation. This chapter describes the significance of RNASeq study over array-based traditional methods, approach to analyze RNASeq data, available methods and tools, challenges associated with the data analysis, application areas, some of the recent advancement made in the area of transcriptome study and its application

    Diverse syntrophic partnerships from deep-sea methane vents revealed by direct cell capture and metagenomics

    Get PDF
    Microorganisms play a fundamental role in the cycling of nutrients and energy on our planet. A common strategy for many microorganisms mediating biogeochemical cycles in anoxic environments is syntrophy, frequently necessitating close spatial proximity between microbial partners. We are only now beginning to fully appreciate the diversity and pervasiveness of microbial partnerships in nature, the majority of which cannot be replicated in the laboratory. One notable example of such cooperation is the interspecies association between anaerobic methane oxidizing archaea (ANME) and sulfate-reducing bacteria. These consortia are globally distributed in the environment and provide a significant sink for methane by substantially reducing the export of this potent greenhouse gas into the atmosphere. The interdependence of these currently uncultured microbes renders them difficult to study, and our knowledge of their physiological capabilities in nature is limited. Here, we have developed a method to capture select microorganisms directly from the environment, using combined fluorescence in situ hybridization and immunomagnetic cell capture. We used this method to purify syntrophic anaerobic methane oxidizing ANME-2c archaea and physically associated microorganisms directly from deep-sea marine sediment. Metagenomics, PCR, and microscopy of these purified consortia revealed unexpected diversity of associated bacteria, including Betaproteobacteria and a second sulfate-reducing Deltaproteobacterial partner. The detection of nitrogenase genes within the metagenome and subsequent demonstration of 15N2 incorporation in the biomass of these methane-oxidizing consortia suggest a possible role in new nitrogen inputs by these syntrophic assemblages
    • …
    corecore