Internucleotide correlations and nucleotide periodicity in Drosophila mtDNA: new evidence for panselective evolution

Analysis for the homogeneity of the distribution of the second base of dinucleotides in relation to the first, whose bases are separated by 0, 1, 2,... 21 nucleotide sites, was performed with the VIH-1 genome (cDNA), the Drosophila mtDNA, the Drosophila Torso gene and the human p-globin gene. These four DNA segments showed highly significant heterogeneities of base distributions that cannot be accounted for by neutral or nearly neutral evolution or by the "neighbor influence" of nucleotides on mutation rates. High correlations are found in the bases of dinucleotides separated by 0, 1 and more number of sites. A periodicity of three consecutive significance values (measured by the x29) was found only in Drosophila mtDNA. This periodicity may be due to an unknown structure or organization of mtDNA. This non-random distribution of the two bases of dinucleotides widespread throughout these DNA segments is rather compatible with panselective evolution and generalized internucleotide co-adaptation.     

The structure of selective dinucleotide interactions and periodicities in D melanogaster mtDNA

Background

We found a strong selective 3-sites periodicity of deviations from randomness of the dinucleotide (DN) distribution, where both bases of DN were separated by 1, 2, K sites in prokaryotes and mtDNA. Three main aspects are studied. I) the specific 3 K-sites periodic structure of the 16 DN. II) to discard the possibility that the periodicity was produced by the highly nonrandom interactive association of contiguous bases, by studying the interaction of non-contiguous bases, the first one chosen each I sites and the second chosen J sites downstream. III) the difference between this selective periodicity of association (distance to randomness) of the four bases with the described fixed periodicities of base sequences.

Results

I) The 16 pairs presented a consistent periodicity in the strength of association of both bases of the pairs; the most deviated pairs are those where G and C are involved and the least deviated ones are those where A and T are involved. II) we found significant non-random interactions when the first nucleotide is chosen every I sites and the second J sites downstream until I = J = 76. III) we showed conclusive differences between these internucleotide association periodicities and sequence periodicities.

Conclusions

This relational selective periodicity is different from sequence periodicities and indicates that any base strongly interacts with the bases of the residual genome; this interaction and periodicity is highly structured and systematic for every pair of bases. This interaction should be destroyed in few generations by recurrent mutation; it is only compatible with the Synthetic Theory of Evolution and agrees with the Wright’s adaptive landscape conception and evolution by shifting balanced adaptive peaks.

Electronic supplementary material

The online version of this article (doi:10.1186/0717-6287-47-18) contains supplementary material, which is available to authorized users.     

MitoAnalyzer, a computer program and interactive web site to determine the effects of single nucleotide polymorphisms and mutations in human mitochondrial DNA

MitoAnalyzer provides information about the effects of single nucleotide polymorphisms (SNPs) and mutations in human mitochondrial DNA (mtDNA). This program determines if a single base pair (bp) change is in the non-coding or coding region, in the first, second or third bp of the codon, in a rRNA, tRNA or a protein, causes an amino acid (aa) change, the nature of that change, the position of the aa change in the protein, and the new aa sequence of the changed protein. Mutations associated with published mitochondrial diseases are noted. This program, thus, facilitates rapid analysis and evaluation of SNPs and mutations in human mtDNA.     

Phylogeography of mitochondrial DNA in South Asian Dolly Varden char Salvelinus curilus Pallas, 1814 (Salmoniformes, Salmonidae): mediated gene introgression?

In 41 individuals of South Asian Dolly Varden char Salvelinus curilus, nucleotide sequences of tRNA-Pro gene fragment (27 bp) and mtDNA control region (483-484 bp) were analyzed. The fish were collected in 20 localities covering virtually the whole range of the species: Kuril Islands, Sakhalin Island. and Primorye. In addition, six individuals of three other char species (S. albus, S. malma, and S. leucomaenis), which are closely related to S. curilus and inhabit the Russian Far East, were examined. In all, we detected 12 different variants of mtDNA haplotypes that formed three distinct groups differing in 14--20 nucleotide positions. The first group consisted of six haplotypes found in S. curilus in Kuril Islands, Sakhalin, and Primorye (mtDNA phylogroup OKHOTSKIA). The second group comprised four haplotypes representing the mtDNA phylogroup BERING, which had been described earlier (Brunner et al, 2001); they were found in S. curilus in Kuril Islands and Sakhalin, as well as in S. albus and S. malma in Kamchatka and northern Kurils. The third group included two haplotypes detected in S. leucomaenis. The existence of two mtDNA lineages (OKHOTSKIA and BERING) in S. curilus from Kurils and Sakhalin was explained by hybridization and DNA transfer from S. malma to S. curilus. The absence of the BERING haplotypes in S. curilus from Primorye water reservoirs is related to the physical isolation of the Sea of Okhotsk and Sea of Japan basins in past epochs. On the basis of comparing phylogenetic trees, constructed from the data on allozyme and mtDNA variation, we suggest that in this case, an indirect transfer of mtDNA in Alpinoid chars--> S. malma-->S. curilus chain could occur.     

Genetic structure of the endangered species <Emphasis Type="Italic">Pinna nobilis</Emphasis> (Mollusca: Bivalvia) inferred from mtDNA sequences

This study examines the population genetic structure of the endangered bivalve Pinna nobilis (Mollusca: Bivalvia), based on novel mtDNA sequences (partial COI and 16S rDNA mtDNA genes). The analyzed nucleotide sequences of COI were 729 bp in size, coding for a 243 amino acid peptide, while the analyzed nucleotide sequences of 16S rDNA were 489 bp in size. These sequences of P. nobilis were the first DNA sequences of the species submitted to any Genetic Data Base. Population samples from four geographic regions from Greece, as well as a population sample of Atrina fragilis (as an outgroup) were used. High values of haplotypic diversity were found in the population samples of P. nobilis, based on the COI sequences. A single base in the analyzed 16S rDNA sequences was different in all analyzed individuals from a single population sample (Chios island) differentiating it from the other ones. These mtDNA sequences could be informative for further genetic analyses of the endangered species, contributing in conservation plans for its protection and/or aquaculture investigations.     

Internal symmetry in nucleotide sequences of genes encoding the dolichol cycle enzymes

In genes alg5, alg8 and swp1 of Saccharomyces cerevisiae, gpt of Schizosaccharomyces pombe and human gene alg6, encoding the dolichol cycle enzymes, a mirror type internal symmetry was found. The symmetry was detected in both complete nucleotide sequences and sequences of the first, second and third nucleotide bases of codons. In the encoding gene regions the density of single- and double-point centres of the internal symmetry for sequences of the second bases was higher in comparison with the sequences of the first and third bases of codons, whereas in the noncoding regions degrees of symmetry of the first, second and third bases sequences did not differ significantly. A clear positive correlation was revealed in the internal symmetry distribution in the second base sequences of codons in genes, on the one hand, and in the gene encoded amino acid sequences, on the other hand. The maximum internal symmetry of gene segments encoding the functionally important regions of proteins was found at the level of the second base sequences. The obtained results corroborate a hypothesis about the determining role of the second bases of codons in encoding amino acid residues. The investigation of internal symmetry in nucleotide sequences has first shown the existence of internal symmetry at the level of gene primary structure.     

Age-dependent 6kb deletion in human liver mitochondrial DNA.

Using PCR technique, restriction mapping and DNA sequencing, we analyzed liver mitochondrial DNA (mtDNA) of 2 stillborn babies and 62 Chinese subjects with non-liver disease from 27 to 86 years old. The results showed an age-dependent 6,063 bp deletion in the liver mtDNA of older subjects. We found a TAACAGAC sequence flanking the 5'-end breakpoint at 7,842 nucleotide position and an imperfect repeat sequence CAACATAC flanking the 3'-end breakpoint at 13,905 nucleotide position. The incidence of the deleted mtDNA was found to increase with age. The deleted mtDNA was not detected in the liver of the stillbirth or blood cells of all the subjects. This is the first account that an age-related 6,063 bp deletion occurs in the liver mtDNA of old humans. The occurrence of this and previously reported 4,977 bp deletions is consistent with our recent finding that liver mitochondrial respiratory functions decline with age and support the hypothesis that continuous accumulation of mtDNA mutations is an important contributor to ageing process in the human.     

The type X collagen gene. Intron sequences split the 5'-untranslated region and separate the coding regions for the non-collagenous amino-terminal and triple-helical domains

Type X collagen, expressed by hypertrophic chondrocytes, consists of homotrimeric molecules with subunits that are only about one-half the size of the polypeptides of fibrillar collagens. In this report we describe for the first time the complete primary structure of type X collagen, based on cloning and sequencing of cDNA and genomic DNA. A comparison between the nucleotide sequences of the cDNA and genomic DNA clones has also allowed determination of the complete exon structure of the type X collagen gene. Our results demonstrate that the primary translation product of the chicken type X collagen mRNA is 682 amino acid residues long with a calculated molecular mass of 67,317 Da for the nonhydroxylated form. This calculated molecular mass is in excellent agreement with the observed electrophoretic mobility of cell-free translation products with both poly(A)+ RNA isolated from chondrocytes as well as RNA transcribed in vitro from a full length cDNA construct. It is also in agreement with the observed size of type X collagen polypeptides isolated from the media of cultured hypertrophic chondrocytes. Thus, our data exclude the possibility of a high molecular weight precursor form of type X collagen. Our results also confirm that the chicken type X gene has a most unusual exon structure when compared to other vertebrate collagen genes. The gene has only three exons. One exon (97 base pairs (bp)), codes for most of the 5'-untranslated region of the mRNA, a second exon (159 bp) codes for the signal peptide and a short non-triple-helical domain, while the third exon (2136 bp) contains the coding region for the entire triple-helix and a large non-triple-helical carboxyl domain.     

Effect of the structure of oligonucleotide substrates on interaction with methylase Eco dam

The interaction of Eco dam methylase with various synthetic oligonucleotide substrates was investigated. The "imperfect" duplexes contained a normal GATC recognition sequence in one chain of the enzyme recognition site and had some defects in the complementary chain, i.e., the absence of one or several nucleotide residues or the presence of S-methyl thiophosphate groups at the 3'-termini. The 3'-S-methyl thiophosphate residue has the same effect on the methylation of oligonucleotide complexes as does the absence of internucleotide phosphate in the analogous complexes. The presence of both GA dinucleotides in the recognition site is necessary for a productive enzyme-substrate interaction. The experimental data suggest that Eco dam methylase does form a symmetrical enzyme-substrate complex which is similar to that formed by type II restriction enzymes.     

Non-random pre-transcriptional evolution in HIV-1. A refutation of the foundational conditions for neutral evolution

The complete base sequence of HIV-1 virus and GP120 ENV gene were analyzed to establish their distance to the expected neutral random sequence. An especial methodology was devised to achieve this aim. Analyses included: a) proportion of dinucleotides (signatures); b) homogeneity in the distribution of dinucleotides and bases (isochores) by dividing both segments in ten and three sub-segments, respectively; c) probability of runs of bases and No-bases according to the Bose-Einstein distribution. The analyses showed a huge deviation from the random distribution expected from neutral evolution and neutral-neighbor influence of nucleotide sites. The most significant result is the tremendous lack of CG dinucleotides (p < 10(-50) ), a selective trait of eukaryote and not of single stranded RNA virus genomes. Results not only refute neutral evolution and neutral neighbor influence, but also strongly indicate that any base at any nucleotide site correlates with all the viral genome or sub-segments. These results suggest that evolution of HIV-1 is pan-selective rather than neutral or nearly neutral.     

On the structure of poly d(A-S4T).

The helical twist of poly d(A-s4T) was determined from the periodicity of the cleavage patterns of the double stranded polydeoxynucleotide adsorbed on calcium phosphate and found to be 14 bp per turn. Both cleavage patterns and 31P NMR spectra indicate a mononucleotide structure rather than an alternating B DNA like poly d(A-T). The failure of nucleosome formation excludes a B type structure. The discrepancy of the mononucleotide structure found in 31P NMR spectra and the dinucleotide structure given by X ray fiber diffraction is explained by an alternating tilt of the planes of the base pairs (base roll) as a consequence of a strong propeller twist. The importance of interstrand stacking interactions of adjacent 4-thiothymidines for the helical stability is discussed.     

The mitochondrial genome structure of the clouded leopard (Neofelis nebulosa).

The complete 16 844 bp mitochondrial genome of Neofelis nebulosa has been sequenced and compared with the complete mitochondrial genomes of Felis catus and the Acinonyx jubatus. The base composition of the mitochondrial genome of N. nebulosa is as follows: A, 5343 bp (31.7%); C, 4441 bp (26.4%); G, 2491 bp (14.8%); T, 4569 bp (27.1%). The genome complement and the gene order of this mitochondrial genome was found to be typical of those reported for other mammals. Several unusual features of this genome, however, were found. First, in protein-coding regions, AT bias in the genome was not prevalent in the third position of codons, as it is in most other mammals, but was found in the second position of codons. Second, in tRNA regions, tRNASer (AGY), which lacked the "DHU" arm, could not be folded into the typical cloverleaf-shaped structure. Third, in the control region, no repetitive sequences (RS)-2 were found. However, RS-2 repetitive motifs usually occurr in the control regions of most great cats. In addition, 4 variable sites were found in CSB-3 of the control region. Fourth, AT content in the control region of the mtDNA from the clouded leopard was lower than it is in other regions.     

Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (∼10.3 bp) at the two ends, with larger (∼11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames’ start and end sites were occupied by nucleosomes.     

Differential accumulations of 4,977 bp deletion in mitochondrial DNA of various tissues in human ageing

Several types of deletions in mitochondrial DNA (mtDNA) have been recetly identified in various tissues of old humans. In order to determine whether there are differences in the incidence and proportion of deleted mtDNAs in different tissues during human ageing, we examined tha 4,977 bp deletion in mtDNA of various tissues from subjects of different ages. Total DNA was extracted from each of the biopsied tissues and was serially diluted by two-fold with distilled water. A 533 bp DNA fragment was amplified by PCR from total mtDNA using a pair of primers L3304-3323 and H3817-3836, and another 524 bp PCR product was amplified from 4,977 bp deleted mtDNA by identical conditions using another pair of primers L8150-8166 and H13631-13650. The maximum dilution fold of each sample that still allowed the ethidium bromide-stained PCR product (533 bp or 524 bp) in the agarose gel to be visible under UV light illumination was taken as the relative abundance of the mtDNA (wild-type or mutant) in the original sample. By this method, we were able to determine the proportion of deleted mtDNA in human tissues. We found that the 4,977 bp deletion started to appear in the second and third decades of life in human muscle and liver tissues. But the deletion was not detectable in the testis until the age of 60 years. Moreover, the proportion of deleted mtDNA varied greatly in different tissues. Among the tissues examined, muscle was found to harbor higher proportin of deleted mtDNA than the other tissues. The average proportion of the 4,977 bp depleted mtDNA of the muscle from subjects over 70 years old was approximately 0.06%, and that of the liver and the testis was 0.0076% and 0.05%, respectively. These findings suggest that the frequency and proportion of the deleted mtDNA in human tissues increase with age and that the mtDNA deletions occur more frequently and abundantly in high energy-demanding tissues during the ageing process of the human.     

Comparison of the sequence specificity of bleomycin cleavage in two slightly different DNA sequences.

The sequence specificity of bleomycin damage was investigated utilising 340 bp alpha-DNA (a middle repetitive sequence in the human genome) as a target sequence. The following significant facts were found:- i) The dinucleotides GT and GC were cleaved on all occasions, GA most of the time, and AT, AC, GG and AA cleaved some of the time; ii) The base immediately 5' to the purine-pyrimidine dinucleotides was found to be statistically highly significant in determining the degree of damage caused by bleomycin, while other nearest neighbour bases had no significant effect; iii) The sequence specificity of bleomycin damage was determined on both strands and it was found that damage on either strand follows the above dinucleotide preference and is independent of the extent of damage on the opposite strand; iv) Bleomycin damage was compared between genomic 340 bp alpha-DNA and a cloned alpha-DNA with eleven base substitutions relative to the "consensus" sequence. There were forty-nine detectable differences in intensity of damage between these two DNA molecules. Although four of the differences can be directly attributed to changes in base sequence, the remaining differences were not at the base substitution sites. Some of the differences were over fifty base pairs from the nearest base substitution. We propose that the majority of these differences are due to microvariation in the structure of DNA with a slightly different DNA sequence.     

暗纹东方鲀线粒体基因组核苷酸全序列测定与分析

以暗纹东方鲀(Takifugu fasciatus)肝的线粒体DNA为模板,参照红鳍东方鲀(T.rubripes)等近源鱼类的线粒体基因组DNA序列,设计合成14对特异引物,进行PCR扩增并测序,首次获得了暗纹东方鲀线粒体基因组全序列。结果表明,暗纹东方鲀线粒体基因组序列全长16 444 bp(GenBank登录号为GQ409967),A+T含量为55.8%,其mtDNA结构与其他脊椎动物相似,由22个tRNA基因、2个rRNA基因、13个蛋白质编码基因和1段819 bp非编码的控制区(D-loop)所组成。蛋白质基因除COⅠ和ND6的起始密码子为GTG、CCT以外,均为典型的起始密码子ATG。ND1、ATPase8、COⅢ、ND4L、ND5、Cyt b使用典型的终止密码子TAA,其他的使用不完全终止密码子。除ND6和tRNAGln、tRNAAla、tRNAAsn、tRNACys、tRNATyr、tRNASer、tRNAGlu、tRNAPro在L-链上编码之外,其余基因均在H-链编码。基因排列顺序与已测定的鲀类一致,这显示了鲀类线粒体基因排列顺序上的保守性。tRNA基因核苷酸长度为64～73nt,预测了22个tRNA基因的二级结构,均呈较为典型的三叶草状。基于19种鲀类mtDNA全序列构建的进化树表明,暗纹东方鲀与红鳍东方鲀、中华东方鲀(T.chinensis)聚成一个姊妹群。结果还支持东方鲀属鱼类为一单系类群。     

GENETIC VARIATION OF KOGIA SPP. WITH PRELIMINARY EVIDENCE FOR TWO SPECIES OF KOGIA SIMA

Concordance between mitochondrial DNA (mtDNA) markers and morphologically based species identifications was examined for the two currently recognized Kogia species. We sequenced 406 base pairs of the control region and 398 base pairs of the cytochrome b gene from 108 Kogia breviceps and 47 K. sima samples. As expecred, the two sister species were reciprocally monophyletic to each other in phylogenetic reconstructions, but within K. sima , we unexpectedly observed another reciprocally monophyletic relationship. The two K. sima clades resolved were phylogeographically concordant with all of the haplotypes in one clade observed solely among specimens sampled from the Atlantic Ocean and with those in the other clade observed solely among specimens sampled from the Indo-Pacific Ocean. These apparently allopatric clades were observed in all phylogenetic reconstructions using the maximum parsimony, maximum likelihood, and neighborjoining algorithms, with the mtDNA gene sequences analyzed separately and combined. The nucleotide diversity for the combined gene sequence haplotypes of the two K. sima clades resolved in our analyses was 0.58% and 1.03% for the Atlantic and Indo-Pacific, respectively, whereas for the two recognized sister species, nucleotide diversity was 1.65% and 4.02% for K. breviceps and K. sima , respectively. The combined gene sequence haplotypes have accumulated 44 fixed base pair differences between the two K. sima clades compared to 20 fixed base pair differences between the two recognized sister species. Although our results are consistent with species-level differences between the two K. sima clades, recognition of a third Kogia species awaits supporting evidence that these two apparently allopatric clades represent reproductively isolated groups of animals.