Complex pattern of coalescence and fast evolution of a mitochondrial rRNA pseudogene in a recent radiation of tiger beetles

Transposed copies of mitochondrial DNA into the nucleus (numts) are widespread, but to date they have not been described from the Coleoptera (beetles). Here we report the discovery of a numt derived from a mitochondrial ribosomal RNA gene in Australian tiger beetles (genus Rivacindela). The loss of function of the numt was confirmed by high proportion of transversions, numerous noncompensatory substitutions in stem regions, and large deletions in functionally important sequences. Phylogenetic analysis of orthologous numt sequences was performed together with the corresponding mtDNA lineage for a study of origination and establishment of the transposed copies in closely related populations and species. All numt sequences were strongly supported to be monophyletic, indicating a single origin of this element. However, populations were polymorphic for the presence of the numt, and phylogenetic trees based on the numt sequences showed inconsistencies with the corresponding mtDNA phylogeny, suggesting slower processes of fixation compared to the mtDNA sequences. In a side-by-side comparison with their mtDNA sister lineage, the nucleotide substitution rate of 1.66 x 10(-8) substitutions/site/year in the numts was approximately equal to the average rate of mtDNA in this group but substantially higher than previous estimates of neutral nuclear rates in vertebrates. The numt clade was affected by several deletions but no insertions, with estimates of nucleotide loss exceeding the rate of nucleotide substitutions by approximately five times. The young age of the Rivacindela numt clade, their absence in species outside of a narrow lineage of related individuals, and the high rate of deletions suggest that insertions do not persist in this group, which is consistent with the view that comparatively small genomes as those of Coleoptera harbor fewer mitochondrial and other nuclear pseudogenes.     

Evolution of the primate beta-globin gene region. High rate of variation in CpG dinucleotides and in short repeated sequences between man and chimpanzee

A 5500 base-pair fragment including the beta-globin gene downstream from codon 122 and about 4000 base-pairs of its 5' flanking sequence was cloned from chimpanzee DNA and thoroughly sequenced before being compared with the corresponding human sequence: 88 point differences (83 substitutions and 5 deletions or insertions of 1 base-pair) were detected as well as seven more important deletion/insertion events. These changes occur preferentially in two kinds of structure. First, 40% of the CpG dinucleotides present in either human or chimpanzee sequences are affected by nucleotide variations. This corresponds to a divergence level considerably higher than that expected. Second, most short repeated sequences found in the 5' extragenic sequence are involved in mutational events (amplification or contraction of the number of basic motifs as well as point substitutions or deletions/insertions of 1 base-pair). Considering the very low level of nucleotide sequence divergence between these two closely related species, our data provide direct evidence for CpG and tandem array instability.     

Rates of DNA Duplication and Mitochondrial DNA Insertion in the Human Genome

The hundreds of mitochondrial pseudogenes in the human nuclear genome sequence (numts) constitute an excellent system for studying and dating DNA duplications and insertions. These pseudogenes are associated with many complete mitochondrial genome sequences and through those with a good fossil record. By comparing individual numts with primate and other mammalian mitochondrial genome sequences, we estimate that these numts arose continuously over the last 58 million years. Our pairwise comparisons between numts suggest that most human numts arose from different mitochondrial insertion events and not by DNA duplication within the nuclear genome. The nuclear genome appears to accumulate mtDNA insertions at a rate high enough to predict within-population polymorphism for the presence/absence of many recent mtDNA insertions. Pairwise analysis of numts and their flanking DNA produces an estimate for the DNA duplication rate in humans of 2.2 × 10^–9 per numt per year. Thus, a nucleotide site is about as likely to be involved in a duplication event as it is to change by point substitution. This estimate of the rate of DNA duplication of noncoding DNA is based on sequences that are not in duplication hotspots, and is close to the rate reported for functional genes in other species.     

云南猕猴mtDNA控制区全序列测定

目的测定云南猕猴线粒体DNA控制区全序列,对其进行鉴定及进化分析。方法利用PCR技术扩增猕猴线粒体DNA控制区全序列,结合GenBank中下载的猕猴参考序列（AY612638）,采用多个生物学软件对序列碱基组成、同源性、转换／颠换比等遗传信息进行分析,并基于邻接法（NJ）和最小进化法（ME）构建系统进化树。结果云南猕猴线粒体DNA控制区全长为（1084-1089）bp,A、T、G和c四种碱基平均含量分别为29．9％、26．9％、12．3％和30．9％,A＋T含量（56．8％）高于G＋C含量（43．2％）。所分析序列间的同源性为91．5％-99．5％,平均核苷酸变异率为4．5％,变异类型包括转换、颠换、插入和缺失4种形式,转换／颠换比值平均为26．1。进化树显示云南猕猴存在两个平行进化的姐妹分支。结论本研究获得了云南猕猴mtDNA控制区全序列,为猕猴进化关系研究及mtDNA控制区功能研究奠定基础。     

青蟹线粒体COI假基因的分离和特征分析

线粒体DNA标记在遗传结构和系统进化研究中得到广泛应用,然而核假基因的存在对此有很大威胁。本文以中国东南沿海的青蟹(Scylla paramamosain)为研究对象,利用线粒体COI基因的通用引物和特异性引物进行扩增,分别得到34个假基因（nuclear mitochondrial pseudogenes, Numts）和5个线粒体COI基因序列。在所获得的34个假基因中共定义了29种单倍型,根据序列的相似度,这些假基因可以分为2类,每类假基因都有各自保守的核苷酸序列。第Ⅰ类假基因存在2处插入序列和1处8 bp的缺失序列,这些位点导致了整个阅读框的移位;在第Ⅱ类假基因和5个线粒体COI序列中只有碱基替换,未发现插入和缺失序列。实验结果分析表明,这两类假基因分别代表了2次核整合事件,即核转移事件的最低值。研究结果提示了     

Probabilistic phylogenetic inference with insertions and deletions

A fundamental task in sequence analysis is to calculate the probability of a multiple alignment given a phylogenetic tree relating the sequences and an evolutionary model describing how sequences change over time. However, the most widely used phylogenetic models only account for residue substitution events. We describe a probabilistic model of a multiple sequence alignment that accounts for insertion and deletion events in addition to substitutions, given a phylogenetic tree, using a rate matrix augmented by the gap character. Starting from a continuous Markov process, we construct a non-reversible generative (birth-death) evolutionary model for insertions and deletions. The model assumes that insertion and deletion events occur one residue at a time. We apply this model to phylogenetic tree inference by extending the program dnaml in phylip. Using standard benchmarking methods on simulated data and a new "concordance test" benchmark on real ribosomal RNA alignments, we show that the extended program dnamlepsilon improves accuracy relative to the usual approach of ignoring gaps, while retaining the computational efficiency of the Felsenstein peeling algorithm.     

Insertion/deletion and nucleotide polymorphism data reveal constraints in Drosophila melanogaster introns and intergenic regions

Our study of nucleotide sequence and insertion/deletion polymorphism in Drosophila melanogaster noncoding DNA provides evidence for selective pressures in both intergenic regions and introns (of the large size class). Intronic and intergenic sequences show a similar polymorphic deletion bias. Insertions have smaller sizes and higher frequencies than deletions, supporting the hypothesis that insertions are selected to compensate for the loss of DNA caused by deletion bias. Analysis of a simple model of selective constraints suggests that the blocks of functional elements located in intergenic sequences are on average larger than those in introns, while the length distribution of relatively unconstrained sequences interspaced between these blocks is similar in intronic and intergenic regions.     

How intron splicing affects the deletion and insertion profile in Drosophila melanogaster

Studies of "dead-on-arrival" transposable elements in Drosophila melanogaster found that deletions outnumber insertions approximately 8:1 with a median size for deletions of approximately 10 bp. These results are consistent with the deletion and insertion profiles found in most other Drosophila pseudogenes. In contrast, a recent study of D. melanogaster introns found a deletion/insertion ratio of 1.35:1, with 84% of deletions being shorter than 10 bp. This discrepancy could be explained if deletions, especially long deletions, are more frequently strongly deleterious than insertions and are eliminated disproportionately from intron sequences. To test this possibility, we use analysis and simulations to examine how deletions and insertions of different lengths affect different components of splicing and determine the distribution of deletions and insertions that preserve the original exons. We find that, consistent with our predictions, longer deletions affect splicing at a much higher rate compared to insertions and short deletions. We also explore other potential constraints in introns and show that most of these also disproportionately affect large deletions. Altogether we demonstrate that constraints in introns may explain much of the difference in the pattern of deletions and insertions observed in Drosophila introns and pseudogenes.     

A DNA sequence evolution analysis generalized by simulation and the markov chain monte carlo method implicates strand slippage in a majority of insertions and deletions

To study the mechanisms for local evolutionary changes in DNA sequences involving slippage-type insertions and deletions, an alignment approach is explored that can consider the posterior probabilities of alignment models. Various patterns of insertion and deletion that can link the ancestor and descendant sequences are proposed and evaluated by simulation and compared by the Markov chain Monte Carlo (MCMC) method. Analyses of pseudogenes reveal that the introduction of the parameters that control the probability of slippage-type events markedly augments the probability of the observed sequence evolution, arguing that a cryptic involvement of slippage occurrences is manifested as insertions and deletions of short nucleotide segments. Strikingly, approximately 80% of insertions in human pseudogenes and approximately 50% of insertions in murids pseudogenes are likely to be caused by the slippage-mediated process, as represented by BC in ABCD --> ABCBCD. We suggest that, in both human and murids, even very short repetitive motifs, such as CAGCAG, CACACA, and CCCC, have approximately 10- to 15-fold susceptibility to insertions and deletions, compared to nonrepetitive sequences. Our protocol, namely, indel-MCMC, thus seems to be a reasonable approach for statistical analyses of the early phase of microsatellite evolution.     

Models of sequence evolution for DNA sequences containing gaps

Most evolutionary tree estimation methods for DNA sequences ignore or inefficiently use the phylogenetic information contained within shared patterns of gaps. This is largely due to the computational difficulties in implementing models for insertions and deletions. A simple way to incorporate this information is to treat a gap as a fifth character (with the four nucleotides being the other four) and to incorporate it within a Markov model of nucleotide substitution. This idea has been dismissed in the past, since it treats a multiple-site insertion or deletion as a sequence of independent events rather than a single event. While this is true, we have found that under many circumstances it is better to incorporate gap information inadequately than to ignore it, at least for topology estimation. We propose an extension to a class of nucleotide substitution models to incorporate the gap character and show that, for data sets (both real and simulated) with short and medium gaps, these models do lead to effective use of the information contained within insertions and deletions. We also implement an ad hoc method in which the likelihood at columns containing multiple-site gaps is downweighted in order to avoid giving them undue influence. The precision of the estimated tree, assessed using Markov chain Monte Carlo techniques to find the posterior distribution over tree space, improves under these five-state models compared with standard methods which effectively ignore gaps.     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

In vitro and in vivo analysis of transcription within the replication region of plasmid pIP501

Summary The tobacco (Nicotiana tabacum) nuclear genome contains long tracts of DNA (i.e. in excess of 18 kb) with high sequence homology to the tobacco plastid genome. Five lambda clones containing these nuclear DNA sequences encompass more than one-third of the tobacco plastid genome. The absolute size of these five integrants is unknown but potentially includes uninterrupted sequences that are as large as the plastid genome itself. An additional sequence was cloned consisting of both nuclear and plastid-derived DNA sequences. The nuclear component of the clone is part of a family of repeats, which are present in about 400 locations in the nuclear genome. The homologous sequences present in chromosomal DNA were very similar to those of the corresponding sequences in the plastid genome. However significant sequence divergence, including base substitutions, insertions and deletions of up to 41 bp, was observed between these nuclear sequences and the plastid genome. Associated with the larger deletions were sequence motifs suggesting that processes such as DNA replication slippage and excision of hairpin loops may have been involved in deletion formation.     

Genetic Diversity of Moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian, and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species origin or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of -Casein and exon 2 of MhcAlal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075–0.15 Myr ago.     

A Comparative Approach Shows Differences in Patterns of Numt Insertion During Hominoid Evolution

Nuclear integrations of mitochondrial DNA (numts) are widespread among eukaryotes, although their prevalence differs greatly among taxa. Most knowledge of numt evolution comes from analyses of whole-genome sequences of single species or, more recently, from genomic comparisons across vast phylogenetic distances. Here we employ a comparative approach using human and chimpanzee genome sequence data to infer differences in the patterns and processes underlying numt integrations. We identified 66 numts that have integrated into the chimpanzee nuclear genome since the human–chimp divergence, which is significantly greater than the 37 numts observed in humans. By comparing these closely related species, we accurately reconstructed the preintegration target site sequence and deduced nucleotide changes associated with numt integration. From >100 species-specific numts, we quantified the frequency of small insertions, deletions, duplications, and instances of microhomology. Most human and chimpanzee numt integrations were accompanied by microhomology and short indels of the kind typically observed in the nonhomologous end-joining pathway of DNA double-strand break repair. Human-specific numts have integrated into regions with a significant deficit of transposable elements; however, the same was not seen in chimpanzees. From a separate data set, we also found evidence for an apparent increase in the rate of numt insertions in the last common ancestor of humans and the great apes using a polymerase chain reaction–based screen. Last, phylogenetic analyses indicate that mitochondrial-numt alignments must be at least 500 bp, and preferably >1 kb in length, to accurately reconstruct hominoid phylogeny and recover the correct point of numt insertion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic diversity of moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of kappa-Casein and exon 2 of MhcALal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075-0.15 Myr ago.     

Evolution of the terminal regions of the Streptomyces linear chromosome

Comparative analysis of the Streptomyces chromosome sequences, between Streptomyces coelicolor, Streptomyces avermitilis, and Streptomyces ambofaciens ATCC23877 (whose partial sequence is released in this study), revealed a highly compartmentalized genetic organization of their genome. Indeed, despite the presence of specific genomic islands, the central part of the chromosome appears highly syntenic. In contrast, the chromosome of each species exhibits large species-specific terminal regions (from 753 to 1,393 kb), even when considering closely related species (S. ambofaciens and S. coelicolor). Interestingly, the size of the central conserved region between species decreases as the phylogenetic distance between them increases, whereas the specific terminal fraction reciprocally increases in size. Between highly syntenic central regions and species-specific chromosomal parts, there is a notable degeneration of synteny due to frequent insertions/deletions. This reveals a massive and constant genomic flux (from lateral gene transfer and DNA rearrangements) affecting the terminal contingency regions. We speculate that a gradient of recombination rate (i.e., insertion/deletion events) toward the extremities is the force driving the exclusion of essential genes from the terminal regions (i.e., chromosome compartmentalization) and generating a fast gene turnover for strong adaptation capabilities.     

Development of a system useful for studying the formation of unstable alleles of IS2

Summary IS2-induced deletions of the gal control region were isolated in a plasmid carrying gal OP-308 :: IS2-7. This contains a 54 basepair long, unstable mini insertion within IS2, thus allowing constitutive expression of the gal structural genes. Deletion PPI is 11.9 kilobasepairs (kb) long and is Gal⁺ because it has retained the mini insertion. In PP4 7.2 kb DNA material including markers gal OP, chlD and pgl are deleted. PP4 has lost the mini insertion and is therefore Gal negative. DNA sequencing of the newly formed junction in PP4 reveals that the deletion terminates precisely at nucleotide 1 of IS2 and that no DNA sequence homology is involved in this IS2-mediated deletion formation. PPI segregates Gal^- clones due to the loss of the mini insertion. One such segregant PPIS and PP4 both give only constitutive Gal⁺ revertants, which consist of the previously known mini insertions and also a new class of supermini inserts within IS2 of about 10 to 20 basepairs long. Therefore, PPIS and PP4 can be used to study various parameters involved in the formation of mini insertions.     

Molecular evolutionary processes and conflicting gene trees: The hominoid case

Molecular evolutionary processes modify DNA over time, creating both newly derived substitutions shared by related descendant lineages (phylogenetic signal) and “false” similarities which confound phylogenetic reconstruction (homoplasy). However, some types of DNA regions, for example those containing tandem duplicate repeats, are preferentially subject to homoplasy-inducing processes such as sporadically occurring concerted evolution and DNA insertion/deletion. This added level of homoplasic “noise” can make DNA regions with repeats less reliable in phylogenetic reconstruction than those without repeats. Most molecular datasets which distinguish among African hominoids support a human-chimpanzee clade; the most notable exception is from the involucrin gene. However, phylogenetic resolution supporting a chimpanzee-gorilla clade is based entirely on involucrin DNA repeat regions. This is problematic because (1) involucrin repeats are difficult to align, and published alignments are contradictory; (2) involucrin repeats are subject to DNA insertion/deletion; (3) gorillas are polymorphic in that some do not have repeats reported to be synapomorphies linking chimpanzees and gorillas. Gene tree/species tree conflicts can occur due to the sorting of ancestrally polymorphic alleles during speciation. Because hominoid females transfer between groups, mitochondrial and nuclear gene flow occur to the same extent, and the probability of conflict between mitochondrial and nuclear gene trees is theoretically low. When hominoid intraspecific mitochondrial variability is taken into account [based on cytochrome oxidase subunit II (COII) gene sequences], humans and chimpanzees are most closely related, showing the same relative degree of separation from gorillas as when single individuals representing species are analyzed. Conflicting molecular phylogenies can be explained in terms of molecular evolutionary processes and sorting of ancient polymorphisms. This perspective can enhance our understanding of hominoid molecular phylogenies. © 1994 Wiley-Liss, Inc.