Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Suggestion for mitochondrial genome assembly

    1.I have denovo assembled my illumina paired-end reads (2x101 bp with abyss genome assembler and got contigs.
    2.I have blasted my denovo assembled contig assembly against my closest reference mitochondria genome (from NCBI organelle genome)
    3.Extracted contigs which have greater than 80% identity with my reference mitochondria genome
    4.Through blast alignment, I extracted 40 putative mitochondrial contigs from my denovo assembled genome


    Now I want my extracted contigs as single mitochondrial genome assembly as in NCBI organelle genome (for e.g.NC_012119.1: http://www.ncbi.nlm.nih.gov/nuccore/...9?report=fasta)


    I tried reference mapping with my close reference, iterative mapping using genious, codon code assembler and also MITObim. But could not able to get
    as single complete genome. Please let me know how to get single complete genome?

  • #2
    you can reoder your denovo assembled contigs into a single molecule (with some gaps) using ABACAS (it needs MUMmer) and the closest reference mitochondria genome, after that you can use IMAGE to fill some gaps with your illumina reads.

    ABACAS ----> http://abacas.sourceforge.net/
    IMAGE ----> https://www.sanger.ac.uk/resources/s...e/pagit/#IMAGE

    Comment


    • #3
      Since you used ABySS, you can use ABySS-Explorer to investigate connections between contigs. It will probably be difficult, since most mitochondrial genomes contain at least some repeated sequences, and sequence shared with plastid genome, and nuclear genome, which may make the assembly ambiguous. And even if your reference is the same species, there could be major rearrangements. If it is a different species, there most certainly are rearrangements. Good luck!

      Comment


      • #4
        Your coverage may be too high for good assembly. Approximately how much coverage do you have? Once you have isolated the mitochondrial contigs, you can map to them and collect all the reads that hit. Then subsample or normalize those reads to an appropriate depth (maybe 100x) and assemble that. Error-correction may help as well.

        Comment


        • #5
          Brian's suggestion sounds great although I'd tweak it to map all your reads to a combination of your contigs and any related reference mitochondria genomes you can find. And when collecting the reads that hit, make sure you include the unmapped mates of reads where only one from the pair hit. Then de novo assemble that data set, subsampling at an appropriate depth as suggested by Brian.

          An alternative is that once you know what an appropriate subsample fraction for mitochondria is from the above step, de novo assemble the same fraction of your entire original data set.

          Comment


          • #6
            Thanks all for your suggestion. When I map my reads to denovo assembled genome (nuclear+mito) i get around 5x coverage.

            Also, I tried by mapping my reads to reference mitochondria genome and I got mapped reads. Then mapped reads where assembled with sub-sampled WGS illumina data (10% random reduction) with MITObim. However, I got lot of gaps in the assembly.

            So, I tried second method, by denovo assembling the WGS illumia paired-end reads and collected putative contigs by blasting against my closest reference mitochondiral genome.

            1.I get around 40 contigs which got hit above 80% similarity and are above 1Kbp.
            2.I get around 3 contigs which are above 3Kbp and 1 contig around 5kbp. Does that mean this one contig 5Kbp will be mitochondrial contig, can I select this for iterative mapping and assembling?

            Also I would like to try Brian's suggestion, to sub-sample with apporpriate depth (100x), how to get the reads which have greater depth (I read some of papers that they reads which have greater depth and denovo assemble). Does your BBMap can get reads for appropriate depth? If not, can you suggest some tools?

            Comment


            • #7
              The BBTools package comes with BBNorm, which can bin reads by depth. Sometimes we use that to separate mitochondrial reads from nuclear reads. e.g:

              bbnorm.sh in=reads.fq passes=1 keepall lowbindepth=20 highbindepth=40 outlow=nuclear.fq outmid=unknown.fq outhigh=mito.fq

              But for this to work you need to know the coverage specifically over the mitochondrial part to figure out what the thresholds should be. You can determine that with my pileup program if you have a sam file (or generate it directly from BBMap):

              pileup.sh in=mapped.sam covstats=stats.txt

              That will tell you the average coverage for every contig, so if you know which are the mito contigs and which are the genomic contigs, it will tell you what limits to use.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Genetic Variation in Immunogenetics and Antibody Diversity
                by seqadmin



                The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                11-06-2024, 07:24 PM
              • seqadmin
                Choosing Between NGS and qPCR
                by seqadmin



                Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                10-18-2024, 07:11 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 11-08-2024, 11:09 AM
              0 responses
              53 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-08-2024, 06:13 AM
              0 responses
              36 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-01-2024, 06:09 AM
              0 responses
              34 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-30-2024, 05:31 AM
              0 responses
              23 views
              0 likes
              Last Post seqadmin  
              Working...
              X