Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Name_less_
    Member
    • Oct 2014
    • 12

    #1

    Questions about small RNA-seq (single-end and strand-specific)

    Hello everyone!
    I'm new analyzing transcriptomic.
    I have small RNA-seq data, single-end from Rat.
    Initially, I mapped the data with Bowtie with next command:

    bowtie -p 50 --mm -n 0 --best --strata -m 1 /media/stdata/bowtie1_index/rn6/rn6_bw1 -q $each_file\.fastq -S $each_file\_n0_m1.sam

    ** -n 0 and -m 1, because we are probing parameters.

    After, I noticed that the kit to make library was "TruSeq® Small RNA Library".
    Reading in the web, I found that:

    "For example, the Illumina TruSeq small RNA sample preparation kit produces strand-specific libraries. The kit specifically selects the RNA species that have a monophosphate group at its 5′-end and a hydroxyl group at 3′-end, a typical structure for miRNAs."
    LINK: https://www.liebertpub.com/doi/full/.../nat.2012.0367

    Then, I think that bowtie also must have option for strand-specific.
    Considering that the data is single-end and default is --fr, Must I need to add option for strand-specific?

    If so, what options are they?
    I found:

    --fr/--rf/--ff -1, -2 mates align fw/rev, rev/fw, fw/fw (default: --fr)
    --nofw/--norc do not align to forward/reverse-complement reference strand

    But I have single end, so I think is necessary to use "--nofw/--norc", but I'm confused which one to choose.
  • Opplatek
    Junior Member
    • Nov 2010
    • 2

    #2
    Hi Name_less_,

    You don't say if you do a genomic alignment to you map directly to a smallRNA database(s).
    In case you do genomic alignment you should not use --nofw/--norc for the Bowtie alignment as you would disallow alignment to the - strand of the genome. Instead, you should use strand-specific counting afterward (for example in FeatureCounts option -s 1).
    If you do database (miRBase, piRBase, ...) you should use the --nofw/--norc option because you assume the sequences in the database are already in the same direction as they could be found in the organism.

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      How Immunogenomics Decodes Immunity’s Genetic Blueprint
      by SEQadmin2




      The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

      This convergence of genetics, immunology, and computation...
      Today, 05:41 AM
    • SEQadmin2
      Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
      by SEQadmin2



      CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

      Despite this, “CRISPR helped turn genome editing from a specialized technique into
      ...
      07-31-2026, 11:01 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 08-24-2026, 10:32 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-20-2026, 11:17 AM
    0 responses
    48 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-18-2026, 10:05 AM
    0 responses
    55 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-13-2026, 12:22 PM
    0 responses
    50 views
    0 reactions
    Last Post SEQadmin2  
    Working...