the mammalian order primate that includes humans , apes , and monkeys in addition to several other organisms can be traced to the late cretaceous period .
the rhesus macaque is in many ways an ideal model organism , being closely related to humans ( sharing a common ancestor about 25 million years ago ) and also sharing similar physiology , neurobiology , and susceptibility to infectious and metabolic diseases .
since both the m. mulatta and h. sapiens genomes have been sequenced , it is known that the evolutionary distance between them is small , with local fluctuations and low divergence , particularly in chromosome x. on average , orthologs have about 97% identity between the genomes both at the nucleotide and amino acid sequence levels .
approximately 50% of the rhesus macaque genome consists of various repetitive sequences , similar to the human genome .
the phenomenon of retrotransposition that occurs in eukaryotic genomes of diverse taxonomic groups is implicated in various human genetic diseases .
insertion sites of many non - long - terminal repeat ( ltr ) retrotransposons play an important role in genome evolution and are distributed throughout the genome .
the phenomenon behind the selection of the insertion sites of these elements has been shown to be correlated with patterns found at pre - insertion loci .
it is well known that mobile elements insertions are capable of altering gene expression , generating genomic deletions , and they can even create new genes and gene families .
existing repetitive elements can also cause ectopic recombinations . despite the overall similarity in retrotransposon mobilization activity in the old world monkeys and hominid lineages ,
the retroelements that are presumed to have had the most dramatic impact in shaping primate genomes are the l1 family of line elements and the alu elements , their partner sines . besides contributing large amounts of dna to many genomes ( including at least 40% of the human genome ) they have also provided new genes , exons and other motifs involved in the physical and sequence structure of chromosomes .
there are instances in which a previous element is now a part of the machinery that regulates gene expression .
repetitive elements account for about 50% of the genome among all of the presently sequenced primate species ( table 1 ) . in m. mulatta ,
two classes of mobile elements are present , class i dna transposons and class ii retrotransposons .
the transposons can also be categorized into different families and subfamilies , based on the relationships between their sequences .
the rhesus family consists of about 320,000 copies of many families of dna transposons and about half a million copies of endogenous retroviruses .
the l1s and alu elements account for most of the lineage specific insertions and these have been playing an important role in shaping the complete genome .
elan was developed earlier as a suite of tools for genome - wide retrotransposon element analysis .
the application of modules of this bioinformatics pipeline is described as follows : ( 1 ) elefinder performs a whole genome distribution analysis of the mges through a blastn search by making the use of perl / bioperl scripts by which the output files are parsed .
it also extracts sequence 100 bp up and downstream at each mge site identified as a preinsertion locus .
( 2 ) dnascanner scans dna sequences such as preinsertion loci and analyses insertion hotspots of elements in detail so as to provide a set of signals or characteristics that are potentially recognized by an element for its insertion .
the m. mulatta genome has a total of 22 chromosomes including x and y ( although the sequencing project did not sequence y ) .
the only y chromosome that has been completely sequenced is of humans , and sequencing of the chimpanzee and mouse y chromosomes is in progress .
the sequencing of the mammalian y chromosomes , of the organisms like rhesus macaque ( macaca mulatta ) , the white - tufted - ear marmoset ( callithrix jacchus ) , the rat ( rattus novergicus ) , the bull ( bos taurus ) and the opossum ( monodelphis domestica ) , is proposed and still under process .
we present here a study of the primate macaca mulatta genome to identify and characterize insertion sites of the two representative retroelements present , and further , comparison with similar features of the human genome ( excluding the y chromosome ) .
the structural and thermodynamic features as well as protein interaction measures are computed in preinsertion loci using the tool dna scanner .
in the human genome , a full - length l1 element is around 6 kb long , and is reported to be the most successful tes in human genome by mass while alu elements , typically ~300 bp long are most successful in terms of copy number . currently , there are three macaque consensus sequences for alu : alumacya3 , alumacyb2 , and alumacyb4 and five of that for l1 : l1p4a , l1p4b , l1p4c , l1p4d , l1p4e in repbase ( version 13.5 ) .
we first compute the nucleotide sequence divergence of l1 ( human ) with respect to the corresponding l1 ( macaca genome ) ( table 2(a ) ) .
as can be seen , pairwise alignment shows a low percentage identity , with the major areas of dissimilarity in the 3 region .
the macaca lineages show higher percentage similarity with each other ; with regard to the human genome , the l1 nearest to its human analog is l1p4d , while that which has diverged the most is l1p4e .
macaca lineages show high similarity with each other and human alu seems to be most closely related to the macaca specific alumcaya3 ( table 2(b ) ) .
tables 3 and 4 give details of the numbers of alu and l1 elements in each chromosome .
the four groups constructed for the pre insertion loci were : intact on both ends , intact on 5 , intact on 3 and intact on neither end .
about 1077410 copies of alus were found uniformly distributed on each chromosome . of these approximately 13.49% ( 145428 ) are truncated on the 5 end , 456752 copies ( 42.39% ) are truncated at both the ends , 30.29% are the truncated at 3 end ( 326447 ) .
a negligible fraction ( 1.6% ) of l1 elements was intact on both sides : of the 94616 elements identified 22% ( 21044 ) were truncated at the 5 end and 1.38% ( 1314 ) were truncated at the 3 end and about 75% at both ends . as in the human genome , there are few functional alu copies in m. mulatta , and their distribution on the various chromosomes follows similar patterns ( fig . 
the x chromosome is known to have an exceptionally high number of the l1 elements ( fig . 
1b ) , which also contained the largest number of truncated elements ( at both the ends ) .
the y - axis represents the frequency of elements found on the different chromosomes ( marked along the x - axis ) .
the y - axis represents the frequency of elements found on the different chromosomes ( marked along the x - axis ) .
the pre - insertion loci were extracted and evaluated for various physicochemical properties as described in our earlier paper .
the positions of the extrema were similar to those seen in other cases in the physicochemical profiles generated by dnascanner .
similar tables were calculated for each m. mulatta chromosome and for all the 14 characteristics extrema were seen in the range of -9 to -11 bp for alu element and that for the l1 element was between -2 to -19 for the majority of the cases . for l1 elements , results are similar and table 5 gives the complete information about the values and the extrema for each of the properties .
figures  2 and 3 show the graphs obtained for 4 physicochemical properties for alu and l1 elements insertion sites in macaca genome .
control sequences were generated by scrambling the positive data set of pre - insertion sequences ; all these above properties gave a featureless distribution ( namely no extrema ) .
another independent set of control sequences that were obtained by randomly selecting genomic sequences of 100 bp also gave similar featureless results .
various signals upstream of the insertion sites of alu in chromosome 1 , for ( a ) stacking energy ,
( b ) at content , ( c ) propeller twist and ( d ) protein induced deformability .
the y axis represents value of the property and the x - axis gives the relative position with respect to the insertion site ( taken to be 0 ) .
various signals upstream of the insertion sites of l1 in chromosome 1 , for ( a ) a rule , ( b ) nucleosomal bending , ( c ) nucleosomal positioning,(d ) protein induced deformability .
the y axis represents value of the property and the x - axis gives the relative position with respect to the insertion site ( taken to be 0 ) . our rationale for the choice of the various parameters
1 ) regions with alternating purines / pyrimidines steps and at rich regions melt more readily .
we found regions of low gc and high at content , indicating that a relatively less energy is required to melt dna near insertion sites , which in turn favors retrotransposition ( figs . 
( 2 ) propeller twist is a property , involved in the distortion of the hydrogen bonds that hold two bases together . regions with specific dinucleotides with large propeller twist , followed by a lower propeller twist were obtained ( fig .  2c ) , which shows that latter regions may be easily distorted and are suitable for insertion .
( 3 ) nucleosomes are involved in dna compacting and providing transcription factors access to the respective regulatory regions .
two different nucleosomal related features , the bending energy / persistence length and the position profiles of the nucleosomes , were studied .
regions with comparatively low energies were obtained , within the upstream areas of the insertion sites .
( 4 ) stacking energy profiles showed a maximum near -10 , indicating that this region is unstable , leading to easy de - stacking of dna sequence , which would thereby enable an easy insertion of alu .
( 5 ) duplex stability is a measure of the relative stability of the dna - duplex structure , which is directly dependent on sequence .
we obtained the region around the -10 position for alu and at -19 for l1 , with a peak , representing a region which would de - stack or melt easily .
( 6 ) dna deformability is an important property , dependent on the sequence and required for interaction with proteins .
the dna deformability was calculated and a region ( at -10 for alu and -18 for l1 ) of low deformability was seen ; this facilitates retrotransposon insertion .
the results of elefinder as well as the information curated in the inside database were used to find the distribution of truncation sites in the whole m. mulatta genome .
the x - axis represents the length of element divided in bins of 25 , i.e.,125 , 2550 , and so on .
the criteria used for the plotting was the occurrence of 5truncated and 3 intact ends and 3 truncated and 5 intact ends in the macaca genome .
the basic aim was to identify the positions in the element sequence , where most of the times truncation occurs for the macaca genome as a whole .
the graph plotted for the 3 truncated ends , the maximum number of hits were obtained in the last bin , i.e. , the 276300 bin has the maximum number of truncations in the macaca genome ( fig . 
similarly , for the 5 truncated category , the maximum number of truncations was found to be in the first bin , i.e. , the 125 bin ( fig . 
( a ) the truncation distribution of alu element in macaca mulatta for bin size 25 and the graph is plotted for 3 end of alu truncated and its 5 end being intact .
( b ) the truncation distribution of alu element in macaca mulatta for bin size 25 and the graph is plotted for 5 end of alu truncated and its 3 end being intact .
previously , we analyzed several genomes for distribution of mges as well as their insertion sites along with the signals facilitating their insertion . in the present study of alus ( sine ) and l1s ( line ) in the macaca genome we found that the insertion sites had physicochemical characteristics that were similar to those observed in other organisms , suggesting that these are generally important .
the present study confirmed that the region 100 bp upstream from alu and l1 insertion sites show statistically significant distinctive properties both in the physical and structural characteristics , as well as in the energetics .
these properties seem to play an important role in the insertion of mge . during insertion ,
a mge causes the target site to distort in a number of ways and requires the co - operative action of a number of proteins to break bonds , unwind the dna , and to nick the target site strand .
it is due to this series of requirements , that the insertion sites for all the chromosomes , show a characteristic set of physicochemical properties , signified by the extremum peaks in each case .
each of the peaks for the several properties carry biological significance and may be used further for the identification of potential new insertion sites . in each case , that the nature of the extremum that has a role in defining the trend of each property was identical for both alu and l1 ( table 5 ) , explaining that the signals that are needed for the insertion of the element is same , and are probably necessary for the insertion to actually occur , although also probably are not sufficient .
detection of the most probable truncation sites for alu elements within the complete macaca genome revealed that during insertion the truncation distribution of alu peaks toward the starting positions in the case of the trailing edge truncation and is reversed for leading edge truncations .
we have found that lines are present in large number in all the primate genomes which includes the recently studied gorilla genome as well as callithrix jaqqus , pan troglodytes , orangutan ( fig . 
there could be an evolutionary link in primate genomes through the spread of lines and sines ; alternately , there may have been a  master  l1 or retrotransposon copy in the genome of last common ancestor of all primates .
the present work adds insight into primate genome architecture , showing the common structural features that promote mge insertion and genome expansion . in future work
, we aim to compare the insertion sites across species- a task that will be facilitated as the diversity of sequenced genomes increases further .
the genome sequence of macaca mulatta was retrieved from the ncbi ( ftp server : ftp://ftp.ncbi.nih.gov/genomes/ ) . the element sequence ( for alu and l1 )
were obtained from repbase and pairwise alignment was performed for the l1 for the human and macaca specific l1s and its lineages using standard procedures .
elefinder was used to find the insertion site of alu and l1 in the macaca genome .
this tool finds the nature , distribution , genomic location , and the site of truncation for each of the insertion sites and performs comparative genome analysis and also generates several set of sequences .
since we have previously shown that intact copies show the presence of signal as compared with the truncated groups , the analysis of the full - length elements , i.e. , of those capable of transposition was also performed .
the tool dna scanner analyses the dna for many physico - chemical properties by using various thermodynamic , protein interactions and sequence - based features , which are beyond the t density and at density . in accordance to the choice of input parameters , the program evaluates a number of properties , using windows that move along the length of the query dna sequence .
since downstream sequences tend to show signals that are not sufficiently statistically significant only upstream sequences were investigated for both alu and l1 elements .
the controls were selected by scrambling sequences , randomly picking sequences from genome as well as gene sequences .
 mobile genetic elements ( mges ) are fragments of dna that can move around within the genome through retrotransposition . 
 these are responsible for various important events such as gene inactivation , transduction , regulation of gene expression and genome expansion . 
 the present work involves the identification and study of the distribution of alu and l1 retrotransposons in the genome of macaca mulatta , an extensively used organism in biomedical studies . 
 we also make comparisons with mge distributions in other primate genomes and study the physicochemical properties of the local dna structure around the transposon insertion site using elan . 
 the present work also includes computational testing of the pre - insertion loci in order to detect unique features based on dna structure , thermodynamic considerations and protein interaction measures . 
 although there is significant sequence divergence between the elements of m. mulatta and h. sapiens , their genome wide distribution is very similar ; comparing the distribution of l1 s in all available x chromosome sequences suggests a common mechanism behind the spread of mge s in primate genomes . 
Introduction
Results
Discussion
Methods