Re: How to find scattered (non-tandem) repeats in DNA sequences?

Article View: bionet.software

Article #2428

#2428

From: Mikhail Fursov
Date: Wed, 07 Apr 2010 06:41

77 lines
7148 bytes

Hi Sunny,

Your results are not related to the nested repeats filter. The nested
repeat in UGENE is the repeat that has all its parts inside of some
other repeat. For example the repeat join(20..30, 120..130) can be
cosidered as nested for the repeat join(10..40, 110, 140)

In your example we see that the region 3+ looks very similar to 4
other regions: 1625+, 1995+, 1806+ and 1488+
Actually this can be considered as a single 5-parts repeat, but UGENE
doesn't have repeats grouping function today and reports 5 pairs. The
grouping is not a simple task for repeats with <100% identity,  but
should be easy to implement for 100% option. We will add this feature
to our TODO list for the future releases.

Now you can also try to use other types of analysis available in UGENE
together with repeats finder. For example you can find ORFs, run
HMMER, BLAST, CDD, TFBS searches for the region between your repeats
or just align them in alignment editor. This can be like a constructor
of complex elements. Note that markup results for the given sequence
are always stored as annotations in Genbank file, so there is no need
to repeat computations every time - just load the annotations for the
sequence you found previous time.

Best regards,
Mikhail.

On Apr 7, 4:19�am, RY Bao <ry_...@yahoo.com> wrote:
> Hi Mike,
>
> Thank you very much! It's working great. I obtained the list of scattered repeats with their positions, length, and distance between the repeats (minimum repeat length = 6bp, repeats identity = 100%, minimum distance between repeats �= 1bp) . I have one quick question. For example, in case of
>
> 1. repeat unit, join(2..9, 1625..1632 )
> 2. repeat unit, join(3..10, 1995..2002)
> 3. repeat unit, join(3..9, 1806..1812)
> 4. repeat unit, join(3..8, 1488..1493)
>
> Those 4 query repeats overlap, however, their target repeats do not overlap at all. Is this because in case of overlapping/nested query repeats, only non-overlapping target repeats are counted? I read the repeat finder plugin section on the instruction page but was not sure. Thanks a lot!
>
> Greetings,
> Sunny
>
> ________________________________
> From: Mikhail Fursov <mike.fur...@gmail.com>
> To: bionet-softw...@moderators.isc.org
> Sent: Tue, April 6, 2010 8:38:06 AM
> Subject: [Bio-software] Re: How to find scattered (non-tandem) repeats in DNA sequences?
>
> On Apr 6, 2:54 am, RY Bao <ry_...@yahoo.com> wrote:
>
> > Hello Everybody,
>
> > I have a 2kb DNA sequence and would like to analyze the scattered repeats in it, which can be allowed with certain mismatches and may be separated by 1 to hundreds of base pairs. Is there a computer program or Web service available to perform this task ?
>
> > Greetings,
> > Sunny Bao (Wayne State University, US)
>
> Hi Sunny Bao,
>
> Here is a description of one of the possible solutions for your:
>
> 1) Open your sequence in UGENE(http://ugene.unipro.ru). Just run
> UGENE and use "Open" button or menu.
>
> 2) Select "Actions->Analyze->Find repeats" from the main menu
>
> 3) Select the repeats parameters: size/identity/distance etc...
>
> 4) Push start button and see the repeated regions as annotations
> (highlighted regions) to the sequence.
>
> Hope it helps,
> Mike.
>
> _______________________________________________
> Bio-soft mailing list
> Bio-s...@net.bio.nethttp://www.bio.net/biomail/listinfo/bio-soft

Technical Details

Message-ID: <mailman.419.1270623833.25217.bio-soft@net.bio.net>

Path:

arc1-ca.pugleaf.net!g2news1.google.com!news2.google.com!Xl.tags.giganews.com!border1.nntp.dca.giganews.com!nntp.giganews.com!local2.nntp.dca.giganews.com!news.giganews.com.POSTED!not-for-mail

References:

<mailman.400.1270497598.25217.bio-soft@net.bio.net> <8c18c392-f0f4-4ee9-825f-b5804eff4bb8@r18g2000yqd.googlegroups.com> <mailman.418.1270589432.25217.bio-soft@net.bio.net>

Show Full Headers

Complete Article Headers

X-Google-Thread: fa6b6,96ca273231ca4a2b
X-Google-NewGroupId: yes
X-Google-Attributes: gid6710630482,domainid0,public,usenet
X-Google-Language: ENGLISH,ASCII
Path: arc1-ca.pugleaf.net!g2news1.google.com!news2.google.com!Xl.tags.giganews.com!border1.nntp.dca.giganews.com!nntp.giganews.com!local2.nntp.dca.giganews.com!news.giganews.com.POSTED!not-for-mail
NNTP-Posting-Date: Wed, 07 Apr 2010 02:04:17 -0500
X-Envelope-From: news@google.com
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma;
 h=domainkey-signature:received:received:received
 :x-authentication-warning:to:path:from:newsgroups:subject:date
 :organization:lines:message-id:references:nntp-posting-host
 :mime-version:content-type:content-transfer-encoding:x-trace
 :x-complaints-to:nntp-posting-date:complaints-to:injection-info
 :user-agent:x-http-useragent;
 bh=YYKu2ScfHqYW7dYum5ji91NjT3Ku1uCtd4M2Z16cnKc=;
 b=p4R5sby1pgvnS33MJyFBB19yA83RO/JKUwu9+K1JOKe1eKeNSUHFskBNXDVUbFoILK
 cmHG0PhFSSNZXwtjQMW3ieke0HxPjJB7wrfMDvgodcNtRDzGfZW6tz8frkSwGKezwAJq
 jr8Hb3Y7C//+UO6z1tPcVwm4xO9XY3jYhfJCY=
DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma;
 h=x-authentication-warning:to:path:from:newsgroups:subject:date
 :organization:lines:message-id:references:nntp-posting-host
 :mime-version:content-type:content-transfer-encoding:x-trace
 :x-complaints-to:nntp-posting-date:complaints-to:injection-info
 :user-agent:x-http-useragent;
 b=JgN5ByVwr0AiEKBvBYW0KWBWkZRj0a1eFNwkjYiuQFtnQFaipw6ICaF1dlzA0hx9mw
 X0A0ZCOZtMi5HiZAvlqxZRfEXIU1JP3miWZvxiX5BorwCZ8Rdb8g4sm3JxHhFNgD1yVT
 ojQasEtrQhrUq4RXkxDJ858Oqw1nFoBX3jRuM=
X-Authentication-Warning: ywd9.prod.google.com: news set sender to
 news@google.com using -f
To: bionet-software@moderators.isc.org
From: Mikhail Fursov <mike.fursov@gmail.com>
Date: Tue, 6 Apr 2010 23:41:21 -0700 (PDT)
Organization: http://groups.google.com
References: <mailman.400.1270497598.25217.bio-soft@net.bio.net> 
 <8c18c392-f0f4-4ee9-825f-b5804eff4bb8@r18g2000yqd.googlegroups.com>
 <mailman.418.1270589432.25217.bio-soft@net.bio.net>
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Complaints-To: groups-abuse@google.com
Injection-Info: h27g2000yqm.googlegroups.com; posting-host=89.31.118.250; 
 posting-account=GUh41QoAAAC2bcZAU6oyL7TynOgTtusR
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/5.0 (Windows; U; Windows NT 5.1; en-GB; rv:1.9.2.3) 
 Gecko/20100401 Firefox/3.6.3,gzip(gfe)
Content-Transfer-Encoding: 8bit
X-MIME-Autoconverted: from quoted-printable to 8bit by net.bio.net id
 o376fTd08621
X-Mailman-Approved-At: Wed, 07 Apr 2010 02:03:51 -0500
X-BeenThere: bio-soft@net.bio.net
X-Mailman-Version: 2.1.4
Precedence: list
subject: Re: How to find scattered (non-tandem) repeats in DNA sequences?
Newsgroups: bionet.software
Message-ID: <mailman.419.1270623833.25217.bio-soft@net.bio.net>
Approved: BIOSCI Approval Key <biosci-approval-key@net.bio.net>
X-Auth: PGPMoose V1.1 PGP bionet.software
 iQCVAwUBS7wuWqCIqsWvyF5FAQHNqwP+KECif20ufy6rgfXYY3SpSmi8vIpANSnr
 jWbNn7ieYvDxpgxm36HhEtkahfC9uOOzNkLuRI0pGgMtk+J49kHShrmcLEV9qCYs
 GQfCtx6Zd/ZYqEy+h5S49zn0bTzZAW/OglSm7KmuLyQSNgpkqllwlntsjOAyLeTN
 VFWxLLLc0WA= =LcSn
Lines: 79
X-Usenet-Provider: http://www.giganews.com
X-Trace: sv3-TVJzkdiy+zcANkI7m42RETdnjb4LNLLy0VkpRgNlFcj0iqN7Mj26ropqhW/rZyoHTkFBfjQnYmtzPQ7!P7u42PQ99S4KTsmZ55UxaeSeXlFJ6oRRmoBdTUjUHgN4ePzb0xkIN2XEj07PGxnVwcxvIE/d69DK!2+Ny
X-Complaints-To: abuse@giganews.com
X-DMCA-Notifications: http://www.giganews.com/info/dmca.html
X-Abuse-and-DMCA-Info: Please be sure to forward a copy of ALL headers
X-Abuse-and-DMCA-Info: Otherwise we will be unable to process your complaint properly
X-Postfilter: 1.3.40

Reply to Parent Message