def circular_assembly_fragments(f, overlap=35, maxlink=40):
fragments = assembly_fragments( f+f[0:1], overlap=overlap, maxlink=maxlink)
if hasattr( fragments[0], "template"):
fragments[0] = _pcr( (fragments[-1].forward_primer, fragments[0].reverse_primer), fragments[0])
return fragments[:-1]
def assembly_fragments(f, overlap=35, maxlink=40):
'''This function return a list of :mod:`pydna.amplicon.Amplicon` objects where
primers have been modified with tails so that the fragments can be fused in
the order they appear in the list by for example Gibson assembly or homologous
recombination.
Given that we have two linear :mod:`pydna.amplicon.Amplicon` objects a and b
we can modify the reverse primer of a and forward primer of b with tails to allow
fusion by fusion PCR, Gibson assembly or in-vivo homologous recombination.
The basic requirements for the primers for the three techniques are the same.
::
<-->
_________ a _________ __________ b ________
/ \\ / \\
agcctatcatcttggtctctgca TTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||| |||||||||||||||||||||||
<gacgt <AGAAA
agcct> TTTAT>
||||||||||||||||||||||| |||||||||||||||||||||||
tcggatagtagaaccagagacgt AAATATAGCGTACTGAGAAGAAA
agcctatcatcttggtctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATAGCGTACTGAGAAGAAA
\\___________________ c ______________________/
Design tailed primers incorporating a part of the next or previous fragment to be assembled.
::
agcctatcatcttggtctctgca
|||||||||||||||||||||||
gagacgtAAATATA
|||||||||||||||||||||||
tcggatagtagaaccagagacgt
TTTATATCGCATGACTCTTCTTT
|||||||||||||||||||||||
ctctgcaTTTATAT
|||||||||||||||||||||||
AAATATAGCGTACTGAGAAGAAA
PCR products with flanking sequences are formed in the PCR process.
::
agcctatcatcttggtctctgcaTTTATAT
||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATA
\\____________/
identical
sequences
____________
/ \\
ctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||
gagacgtAAATATAGCGTACTGAGAAGAAA
The fragments can be fused by any of the techniques mentioned earlier to form c:
::
agcctatcatcttggtctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATAGCGTACTGAGAAGAAA
The first argument of this function is a list of sequence objects containing
Amplicons and other similar objects.
**At least every second sequence object needs to be an Amplicon**
This rule exists because if a sequence object is that is not a PCR product
is to be fused with another fragment, that other fragment needs to be an Amplicon
so that the primer of the other object can be modified to include the whole stretch
of sequence homology needed for the fusion. See the example below where a is a
non-amplicon (a linear plasmid vector for instance)
::
_________ a _________ __________ b ________
/ \\ / \\
agcctatcatcttggtctctgca <--> TTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||| |||||||||||||||||||||||
tcggatagtagaaccagagacgt <AGAAA
TTTAT>
|||||||||||||||||||||||
<--> AAATATAGCGTACTGAGAAGAAA
agcctatcatcttggtctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATAGCGTACTGAGAAGAAA
\\___________________ c ______________________/
In this case only the forward primer of b is fitted with a tail with a part a:
::
agcctatcatcttggtctctgca
|||||||||||||||||||||||
tcggatagtagaaccagagacgt
TTTATATCGCATGACTCTTCTTT
|||||||||||||||||||||||
<AGAAA
tcttggtctctgcaTTTATAT
|||||||||||||||||||||||
AAATATAGCGTACTGAGAAGAAA
PCR products with flanking sequences are formed in the PCR process.
::
agcctatcatcttggtctctgcaTTTATAT
||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATA
\\____________/
identical
sequences
____________
/ \\
ctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||
gagacgtAAATATAGCGTACTGAGAAGAAA
The fragments can be fused by for example Gibson assembly:
::
agcctatcatcttggtctctgcaTTTATAT
||||||||||||||||||||||||||||||
tcggatagtagaacca
TCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||
gagacgtAAATATAGCGTACTGAGAAGAAA
to form c:
::
agcctatcatcttggtctctgcaTTTATATCGCATGACTCTTCTTT
||||||||||||||||||||||||||||||||||||||||||||||
tcggatagtagaaccagagacgtAAATATAGCGTACTGAGAAGAAA
The first argument of this function is a list of sequence objects containing
Amplicons and other similar objects.
The overlap argument controls how many base pairs of overlap required between
adjacent sequence fragments. In the junction between Amplicons, tails with the
length of about half of this value is added to the two primers
closest to the junction.
::
> <
Amplicon1
Amplicon2
> <
⇣
> <-
Amplicon1
Amplicon2
-> <
In the case of an Amplicon adjacent to a Dseqrecord object, the tail will
be twice as long (1*overlap) since the
recombining sequence is present entirely on this primer:
::
Dseqrecd1
Amplicon1
> <
⇣
Dseqrecd1
Amplicon1
--> <
Note that if the sequence of DNA fragments starts or stops with an Amplicon,
the very first and very last prinmer will not be modified i.e. assembles are
always assumed to be linear. There are simple tricks around that for circular
assemblies depicted in the last two examples below.
The maxlink arguments controls the cut off length for sequences that will be
synhtesized by adding them to primers for the adjacent fragment(s). The
argument list may contain short spacers (such as spacers between fusion proteins).
::
Example 1: Linear assembly of PCR products (pydna.amplicon.Amplicon class objects) ------
> < > <
Amplicon1 Amplicon3
Amplicon2 Amplicon4
> < > <
⇣
pydna.design.assembly_fragments
⇣
> <- -> <- pydna.assembly.Assembly
Amplicon1 Amplicon3
Amplicon2 Amplicon4 ➤ Amplicon1Amplicon2Amplicon3Amplicon4
-> <- -> <
Example 2: Linear assembly of alternating Amplicons and other fragments
> < > <
Amplicon1 Amplicon2
Dseqrecd1 Dseqrecd2
⇣
pydna.design.assembly_fragments
⇣
> <-- --> <-- pydna.assembly.Assembly
Amplicon1 Amplicon2
Dseqrecd1 Dseqrecd2 ➤ Amplicon1Dseqrecd1Amplicon2Dseqrecd2
Example 3: Linear assembly of alternating Amplicons and other fragments
Dseqrecd1 Dseqrecd2
Amplicon1 Amplicon2
> < --> <
⇣
pydna.design.assembly_fragments
⇣
pydna.assembly.Assembly
Dseqrecd1 Dseqrecd2
Amplicon1 Amplicon2 ➤ Dseqrecd1Amplicon1Dseqrecd2Amplicon2
--> <-- --> <
Example 4: Circular assembly of alternating Amplicons and other fragments
-> <==
Dseqrecd1 Amplicon2
Amplicon1 Dseqrecd1
--> <-
⇣
pydna.design.assembly_fragments
⇣
pydna.assembly.Assembly
-> <==
Dseqrecd1 Amplicon2 -Dseqrecd1Amplicon1Amplicon2-
Amplicon1 ➤ | |
--> <- -----------------------------
------ Example 5: Circular assembly of Amplicons
> < > <
Amplicon1 Amplicon3
Amplicon2 Amplicon1
> < > <
⇣
pydna.design.assembly_fragments
⇣
> <= -> <-
Amplicon1 Amplicon3
Amplicon2 Amplicon1
-> <- +> <
⇣
make new Amplicon using the Amplicon1.template and
the last fwd primer and the first rev primer.
⇣
pydna.assembly.Assembly
+> <= -> <-
Amplicon1 Amplicon3 -Amplicon1Amplicon2Amplicon3-
Amplicon2 ➤ | |
-> <- -----------------------------
Parameters
----------
f : list of :mod:`pydna.amplicon.Amplicon` and other Dseqrecord like objects
list Amplicon and Dseqrecord object for which fusion primers should be constructed.
overlap : int, optional
Length of required overlap between fragments.
maxlink : int, optional
Maximum length of spacer sequences that may be present in f. These will be included in tails for designed primers.
Returns
-------
seqs : list of :mod:`pydna.amplicon.Amplicon` and other Dseqrecord like objects :mod:`pydna.amplicon.Amplicon` objects
::
[Amplicon1,
Amplicon2, ...]
Examples
--------
>>> from pydna.dseqrecord import Dseqrecord
>>> from pydna.design import primer_design
>>> a=primer_design(Dseqrecord("atgactgctaacccttccttggtgttgaacaagatcgacgacatttcgttcgaaacttacgatg"))
>>> b=primer_design(Dseqrecord("ccaaacccaccaggtaccttatgtaagtacttcaagtcgccagaagacttcttggtcaagttgcc"))
>>> c=primer_design(Dseqrecord("tgtactggtgctgaaccttgtatcaagttgggtgttgacgccattgccccaggtggtcgtttcgtt"))
>>> from pydna.design import assembly_fragments
>>> # We would like a circular recombination, so the first sequence has to be repeated
>>> fa1,fb,fc,fa2 = assembly_fragments([a,b,c,a])
>>> # Since all fragments are Amplicons, we need to extract the rp of the 1st and fp of the last fragments.
>>> from pydna.amplify import pcr
>>> fa = pcr(fa2.forward_primer, fa1.reverse_primer, a)
>>> [fa,fb,fc]
[Amplicon(100), Amplicon(101), Amplicon(102)]
>>> from pydna.assembly import Assembly
>>> assemblyobj = Assembly([fa,fb,fc])
>>> assemblyobj
Assembly
fragments....: 100bp 101bp 102bp
limit(bp)....: 25
G.nodes......: 6
algorithm....: common_sub_strings
>>> assemblyobj.assemble_linear()
[Contig(-231), Contig(-166), Contig(-36)]
>>> assemblyobj.assemble_circular()[0].cseguid()
'V3Mi8zilejgyoH833UbjJOtDMbc'
>>> (a+b+c).looped().cseguid()
'V3Mi8zilejgyoH833UbjJOtDMbc'
>>> print(assemblyobj.assemble_circular()[0].figure())
-|100bp_PCR_prod|36
| \\/
| /\\
| 36|101bp_PCR_prod|36
| \\/
| /\\
| 36|102bp_PCR_prod|36
| \\/
| /\\
| 36-
| |
--------------------------------------------------------
>>>
'''
# sanity check for arguments
nf = [item for item in f if len(item)>maxlink]
if not all(hasattr(i[0],"template") or hasattr(i[1],"template") for i in zip(nf,nf[1:])):
raise ValueError("Every second fragment larger than maxlink has to be an Amplicon object")
_module_logger.debug("### assembly fragments ###")
_module_logger.debug("overlap = %s", overlap)
_module_logger.debug("max_link = %s", maxlink)
f = [_copy.copy(f) for f in f]
first_fragment_length = len(f[0])
last_fragment_length = len(f[-1])
if first_fragment_length<=maxlink:
# first fragment should be removed and added to second fragment (new first fragment) forward primer
f[1].forward_primer = f[0].seq._data + f[1].forward_primer
_module_logger.debug("first fragment removed since len(f[0]) = %s", first_fragment_length)
f=f[1:]
else:
_module_logger.debug("first fragment stays since len(f[0]) = %s", first_fragment_length)
if last_fragment_length<=maxlink:
f[-2].reverse_primer = f[-1].seq.reverse_complement()._data + f[-2].reverse_primer
f=f[:-1]
_module_logger.debug("last fragment removed since len(f[%s]) = %s", len(f), last_fragment_length)
else:
_module_logger.debug("last fragment stays since len(f[%s]) = %s", len(f),last_fragment_length)
empty = _Dseqrecord("")
_module_logger.debug(f)
_module_logger.debug("loop through fragments in groups of three:")
tail_length = _math.ceil(overlap/2)
for i in range(len(f)-1):
first = f[i]
secnd = f[i+1]
secnd_len = len(secnd)
_module_logger.debug( "first = %s", str(first.seq))
_module_logger.debug( "secnd = %s", str(secnd.seq))
if secnd_len <= maxlink:
_module_logger.debug("secnd is smaller or equal to maxlink; should be added to primer(s)")
third = f[i+2]
_module_logger.debug( "third = %s", str(third.seq))
if hasattr(f[i], "template") and hasattr(third, "template"):
_module_logger.debug("secnd is is flanked by amplicons, so half of secnd should be added each flanking primers")
first.reverse_primer = secnd.seq.reverse_complement()._data[secnd_len//2:] + first.reverse_primer
third.forward_primer = secnd.seq._data[secnd_len//2:] + third.forward_primer
lnk = (third.seq.reverse_complement()._data+secnd.reverse_complement().seq._data[:secnd_len//2])[-tail_length:]
_module_logger.debug("1 %s", lnk)
first.reverse_primer = lnk + first.reverse_primer
lnk = (first.seq._data + secnd.seq._data[:secnd_len//2])[-tail_length:]
_module_logger.debug("2 %s", lnk)
third.forward_primer = lnk + third.forward_primer
elif hasattr(first , "template"):
first.reverse_primer = secnd.seq.reverse_complement()._data + first.reverse_primer
lnk = str(third.seq[:overlap].reverse_complement())
first.reverse_primer = lnk + first.reverse_primer
elif hasattr(third , "template"):
third.forward_primer = secnd.seq._data + third.forward_primer
lnk = str(first.seq[-overlap:])
third.forward_primer = lnk + third.forward_primer
secnd=empty
f[i+2] = third
else: # secnd is larger than maxlink
if hasattr(first, "template") and hasattr(secnd, "template"):
lnk = str(first.seq[-tail_length:])
#_module_logger.debug("4 %s", lnk)
secnd.forward_primer = lnk + secnd.forward_primer
lnk = str(secnd.seq[:tail_length].reverse_complement())
#_module_logger.debug("5 %s", lnk)
first.reverse_primer = lnk + first.reverse_primer
elif hasattr(first , "template"):
lnk = str(secnd.seq[:overlap].reverse_complement())
#_module_logger.debug("6 %s", lnk)
first.reverse_primer = lnk + first.reverse_primer
elif hasattr(secnd , "template"):
lnk = str(first.seq[-overlap:])
#_module_logger.debug("7 %s", lnk)
secnd.forward_primer = lnk + secnd.forward_primer
f[i] = first
f[i+1] = secnd
_module_logger.debug("loop ended")
f = [item for item in f if len(item)]
return [_pcr(p.forward_primer, p.reverse_primer, p.template) if hasattr(p, "template") else p for p in f]