I would have thought DNA composition would have played some role (for example, if 10% of the genome came from different organisms, with a different codon bias, that might result in a non-50:50 ratio even if that 10% had a random insertion direction).
I like this idea, because it would seem to follow that if codon bias affected the 50:50 ratio, then it would explain why organisms that do horizontal gene transfer break Chargraff's second parity rule, and why organelles - genomes highly dependant and interacting with an "external" genome - are the worst offenders. In short, my theory would be that the more your genome interacts with others, the less optimised your codon/nucleoside ratios are, the further you deviate from 50:50 even in the presence of random insertion direction. So is matted's random-insertion-direction sufficient to give a 50:50 ratio?
The code below tests random 50:50 reverse-complimenting of DNA, of user-editable codon usage:
import string
import random
def rc(DNA):
return DNA.translate(string.maketrans('ACGT','TGCA'))[::-1]
genome = [
'AAA' for _ in xrange(1000) ]+[
'AAT' for _ in xrange(10000) ]+[
'ATT' for _ in xrange(100000) ]+[
'TTT' for _ in xrange(1000000) ]+[
'CCC' for _ in xrange(1000000) ]+[
'CCG' for _ in xrange(100000) ]+[
'CGG' for _ in xrange(10000) ]+[
'GGG' for _ in xrange(1000) ]
before = ''.join(genome)
A,C,G,T = before.count('A'),before.count('C'),before.count('G'),before.count('T')
print 'Before random insertion:'
print float(A)/(A+T), float(C)/(C+G)
for position,codon in enumerate(genome):
if random.randint(1,100) > 50: # Flip codon 50%
genome[position] = rc(codon) # of the time
after = ''.join(genome)
A,C,G,T = after.count('A'),after.count('C'),after.count('G'),after.count('T')
print 'After random insertion:'
print float(A)/(A+T), float(C)/(C+G)
Sample output:
Before random insertion:
0.036903690369 0.963096309631
After random insertion:
0.500247890819 0.499687968797
So as you can see, even if you start with a highly 'deviant' genome, you still end up with near 50:50 ratios if every codon has a 50% chance of being reverse-complimented. Therefore, I would say my idea of codon usage is false, and matted's answer is correct and sufficient to explain Chargraff's second rule :)
It also stops working when the codon number get small - obviously, because if you only had a single, highly biased gene, it doesnt matter how you orientate it, it's going to result in a highly-biased nucleoside composition. In this test i'm flipping every codon, but in nature you can only flip whole genes. So this suggests that the genomes that break the second rule are probably due to the number and size of the genes, or non-random insertion direction, and bad luck.
hi,
look here - https://en.wikipedia.org/wiki/Chargaff%27s_rules
Chargaff experiment has been done on whole genom , which contains both strand + and -. In this case, complementary explains why A/T and C/G are constant..
This is not the same for my exemple. I working on a single strand !
Sacha, you're referring to his first rule
I really don't see the second rule as relevant or important. These rules were thought up 60 years ago, back when we didn't know how cells stored their genetic material. Its like aliens coming down to earth and seeing that books can be transfered between humans, and is somehow used to transfer information.
The first rule is like saying that there is a relationship between the amount of black ink on the page, and the number of little black dots "." - this gives some insight into how our language works.
The second rule is like saying different o'Reilly books which convey different information, must have a different number of pages - (except for the small books which occasionally have the same number of pages).
This isn't interesting or relevant, and is likely due to chance. Perhaps this was interesting 50 years ago when we didnt know how books worked, but these days we can read and write our own books...
Run your script on duomers (AA, AC, AT, AG, CA .... ) and get back to us on the next really interesting question of why one of them is totally weird.
If you don't have the answer, please don't answer. Thanks!
Not good etiquette. The forum encourages discussion, and the current relevance of a question is as important as, if not more important than, the answer itself. One of the primary skills you learn as you move up the ladder is when to solve a problem ab initio and when to use an existing solution. You also learn when a problem no longer needs solving because the question it answers is of relatively low significance.
In any case, we do not encourage comments that discourage participation. You may want to stick around, contribute more and earn people's trust before trying to change how stuff works on an online discussion forum.
Seems I m not the only one who ask this "unrelevant" question ...
http://www.sciencedirect.com.gate2.inist.fr/science/article/pii/S0006291X05027130
This link cannot be accessed by the general public. It seems to route through an authentication gateway.
This is related to Chargaff's second rule. Have a look on wiki
You are right! Sorry! I can explain the first rule, but cannot explain the second one.. I m reading it :
Second parity rule
The second rule holds that both %A ~ %T and %G ~ %C are valid for each of the two DNA strands. [3]This describes only a global feature of the base composition in a single DNA strand.[4]
If you read further, the page says there is evidence from 2006 that the second rule was proven for major kingdoms. I guess you could dig a bit deeper to examine the actual evidence, but the statement lends itself to easy understanding, so there should be no problem comprehending it. "%A~%T and %G~%C" is the logical representation of "In most
eukaryoticsingle DNA strands, the percentage of A is approximately equal to the percentage of T, and the percentage of G is approximately equal to the percentage of C"