Pervasive migration of organellar DNA to the nucleus in plants

A surprisingly large number of plant nuclear DNA sequences inferred to be remnants of chloroplast and mitochondrial DNA migration events were detected through computer-assisted database searches. Nineteen independent organellar DNA insertions, with a median size of 117 by (range of 38 to >785 bp), occur in the proximity of 15 nuclear genes. One fragment appears to have been passed through a RNA intermediate, based on the presence of an edited version of the mitochondrial gene in the nucleus. Tandemly arranged fragments from disparate regions of organellar genomes and from different organellar genomes indicate that the fragments joined together from an intracellular pool of RNA and/or DNA before they integrated into the nuclear genome. Comparisons of integrated sequences to genes lacking the insertions, as well as the occurrence of coligated fragments, support a model of random integration by end joining. All transferred sequences were found in noncoding regions, but the positioning of organellar-derived DNA in introns, as well as regions 5 and 3 to nuclear genes, suggests that the random integration of organellar DNA has the potential to influence gene expression patterns. A semiquantitative estimate was performed on the amount of organellar DNA being transferred and assimilated into the nucleus. Based on this database survey, we estimate that 3–7% of the plant nuclear genomic sequence files contain organellar-derived DNA. The timing and the magnitude of genetic flux to the nuclear genome suggest that random integration is a substantial and ongoing process for creating sequence variation.Correspondence to: J.L. Blanchard     

Origin, evolution and genetic effects of nuclear insertions of organelle DNA

In eukaryotes, nuclear genomes are subject to an influx of DNA from mitochondria and plastids. The nuclear insertion of organellar sequences can occur during the illegitimate repair of double-stranded breaks. After integration, nuclear organelle DNA is modified by point mutations, and by deletions. Insertion of organelle DNA into nuclear genes is not rare and can potentially have harmful effects. In humans, some insertions of nuclear mitochondrial DNA are associated with heritable diseases. It remains to be determined whether nuclear organelle DNA can contribute beneficially to gene evolution.     

Plastid DNA in the nucleus: New genes for old

Nuclear genomes of eukaryotes are bombarded by a continuous deluge of organellar DNA which contributes significantly to eukaryote evolution. Here, we present a new PCR-based method that allows the specific amplification of nuclear integrants of organellar DNA (norgs) by exploiting recent deletions present in organellar genome sequences. We have used this method to amplify nuclear integrants of plastid DNA (nupts) from the nuclear genomes of several nicotiana species and to study the evolutionary forces acting upon these sequences. The role of nupts in endosymbiotic evolution and the different genetic factors influencing the time available for a chloroplastic gene to be functionally relocated in the nucleus are discussed.     

Detection of mitochondrial insertions in the nucleus (NuMts) of Pleistocene and modern muskoxen

Background  

Nuclear insertions of mitochondrial sequences (NuMts) have been identified in a wide variety of organisms. Trafficking of genetic material from the mitochondria to the nucleus has occurred frequently during mammalian evolution and can lead to the production of a large pool of sequences with varying degrees of homology to organellar mitochondrial DNA (mtDNA) sequences. This presents both opportunities and challenges for forensics, population genetics, evolutionary genetics, conservation biology and the study of DNA from ancient samples. Here we present a case in which difficulties in ascertaining the organellar mtDNA sequence from modern samples hindered their comparison to ancient DNA sequences.     

Assembly-free quantification of vagrant DNA inserts

Inserts of DNA from extranuclear sources, such as organelles and microbes, are common in eukaryote nuclear genomes. However, sequence similarity between the nuclear and extranuclear DNA, and a history of multiple insertions, make the assembly of these regions challenging. Consequently, the number, sequence and location of these vagrant DNAs cannot be reliably inferred from the genome assemblies of most organisms. We introduce two statistical methods to estimate the abundance of nuclear inserts even in the absence of a nuclear genome assembly. The first (intercept method) only requires low-coverage (<1×) sequencing data, as commonly generated for population studies of organellar and ribosomal DNAs. The second method additionally requires that a subset of the individuals carry extranuclear DNA with diverged genotypes. We validated our intercept method using simulations and by re-estimating the frequency of human NUMTs (nuclear mitochondrial inserts). We then applied it to the grasshopper Podisma pedestris, exceptional for both its large genome size and reports of numerous NUMT inserts, estimating that NUMTs make up 0.056% of the nuclear genome, equivalent to >500 times the mitochondrial genome size. We also re-analysed a museomics data set of the parrot Psephotellus varius, obtaining an estimate of only 0.0043%, in line with reports from other species of bird. Our study demonstrates the utility of low-coverage high-throughput sequencing data for the quantification of nuclear vagrant DNAs. Beyond quantifying organellar inserts, these methods could also be used on endosymbiont-derived sequences. We provide an R implementation of our methods called “vagrantDNA” and code to simulate test data sets.     

家马核基因组中线粒体核内插入序列分析

在各种真核生物核基因组中,存在一些由线粒体基因组转移进入核基因组中的DNA片段,这些被认为是分子化石的片段叫做线粒体核内插入序列(Numt)。由于Numt与真实的线粒体序列高度相似,因此它的存在必然会成为PCR扩增线粒体DNA的不利因素。利用已经公布的家马(Equus caballus)基因组序列(2007年9月公布,GenBank登录号为NC_009144-NC_009175)对家马Numt进行了深入分析,共发现200个可能的Numt,长度范围为29到3727bp,其中有10个的长度大于800bp。分析结果显示由于不存在线粒体控制区域的疑似Numt,因此对基于此区域的群体遗传学研究不会产生影响。本研究还发现在家马进化过程中,第1号和27号染色体更倾向于接受线粒体序列的转移。以上结果将为今后马科动物的研究提供重要的参考信息,有助于避免在线粒体DNA研究中由于Numt污染的存在而得出错误的实验结果。     

Chloroplast DNA insertions into the nuclear genome of rice: the genes,sites and ages of insertion involved

Rice (Oryza sativa) is one of three predominant grain crops, and its nuclear and organelle genomes have been sequenced. Following genome analysis revealed many exchanges of DNA sequences between the nuclear and organelle genomes. In this study, a total of 45 chloroplast DNA insertions more than 2 kb in length were detected in rice nuclear genome. A homologous recombination mechanism is expected for those chloroplast insertions with high similarity between their flanking sequences. Only five chloroplast insertions with high sequence similarity between two flanking sequences from an insertion were found in the 45 insertions, suggesting that rice might follow the non-homologous end-joining (NHEJ) repair of double-stranded breaks mechanism, which is suggested to be common to all eukaryotes. Our studies indicate that the most chloroplast insertions occurred at a nuclear region characterized by a sharp change of repetitive sequence density. One potential explanation is that regions such as this might be susceptible target sites or “hotspots” of DNA damage. Our results also suggest that the insertion of retrotransposon elements or non-chloroplast DNA into chloroplast DNA insertions may contribute significantly to their fragmentation process. Moreover, based on chloroplast insertions in nuclear genomes of two subspecies (indica and japonica) of cultivated rice, our results strongly suggest that they diverged during 0.06–0.22 million years ago. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Odintifier - A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.     

NUPTs in sequenced eukaryotes and their genomic organization in relation to NUMTs

NUPTs (nuclear plastid DNA) derive from plastid-to-nucleus DNA transfer and exist in various plant species. Experimental data imply that the DNA transfer is an ongoing, highly frequent process, but for the interspecific diversity of NUPTs, no clear explanation exists. Here, an inventory of NUPTs in the four sequenced plastid-bearing species and their genomic organization is presented. Large genomes with a predicted low gene density contain more NUPTs. In Chlamydomonas and Plasmodium, DNA transfer occurred but was limited, probably because of the presence of only one plastid per cell. In Arabidopsis and rice, NUPTs are frequently organized as clusters. Tight clusters can contain both NUPTs and NUMTs (nuclear mitochondrial DNA), indicating that preNUPTs and preNUMTs might have concatamerized before integration. The composition of such a hypothetical preNUPT-preNUMT pool seems to be variable, as implied by substantially different NUPTs:NUMTs ratios in different species. Loose clusters can span several dozens of kbps of nuclear DNA, and they contain markedly more NUPTs or NUMTs than expected from a random genomic distribution of nuclear organellar DNA. The level of sequence similarity between NUPTs/NUMTs and plastid/mitochondrial DNA correlates with the size of the integrant. This implies that original insertions are large and decay over evolutionary time into smaller fragments with diverging sequences. We suggest that tight and loose clusters represent intermediates of this decay process.     

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.     

Dual-domain, dual-targeting organellar protein presequences in Arabidopsis can use non-AUG start codons

The processes accompanying endosymbiosis have led to a complex network of interorganellar protein traffic that originates from nuclear genes encoding mitochondrial and plastid proteins. A significant proportion of nucleus-encoded organellar proteins are dual targeted, and the process by which a protein acquires the capacity for both mitochondrial and plastid targeting may involve intergenic DNA exchange coupled with the incorporation of sequences residing upstream of the gene. We evaluated targeting and sequence alignment features of two organellar DNA polymerase genes from Arabidopsis thaliana. Within one of these two loci, protein targeting appeared to be plastidic when the 5' untranslated leader region (UTR) was deleted and translation could only initiate at the annotated ATG start codon but dual targeted when the 5' UTR was included. Introduction of stop codons at various sites within the putative UTR demonstrated that this region is translated and influences protein targeting capacity. However, no ATG start codon was found within this upstream, translated region, suggesting that translation initiates at a non-ATG start. We identified a CTG codon that likely accounts for much of this initiation. Investigation of the 5' region of other nucleus-encoded organellar genes suggests that several genes may incorporate upstream sequences to influence targeting capacity. We postulate that a combination of intergenic recombination and some relaxation of constraints on translation initiation has acted in the evolution of protein targeting specificity for those proteins capable of functioning in both plastids and mitochondria.     

Evolutionary rates of insertion and deletion in noncoding nucleotide sequences of primates

Insertions and deletions are responsible for gaps in aligned nucleotide sequences, but they have been usually ignored when the number of nucleotide substitutions was estimated. We compared six sets of nuclear and mitochondrial noncoding DNA sequences of primates and obtained the estimates of the evolutionary rate of insertion and deletion. The maximum-parsimony principle was applied to locate insertions and deletions on a given phylogenetic tree. Deletions were about twice as frequent as insertions for nuclear DNA, and single-nucleotide insertions and deletions were the most frequent in all events. The rate of insertion and deletion was found to be rather constant among branches of the phylogenetic tree, and the rate (approximately 2.0/kb/Myr) for mitochondrial DNA was found to be much higher than that (approximately 0.2/kb/Myr) for nuclear DNA. The rates of nucleotide substitution were about 10 times higher than the rate of insertion and deletion for both nuclear and mitochondrial DNA.      

Aphidicolin and ethidium bromide as tools to verify dependence of early post-germinative growth on nuclear and organellar DNA synthesis

Embryos of Haplopappus gracilis (Nutt.) Gray were sown in aphidicolin (inhibitor of nuclear DNA synthesis) used alone or in combination with ethidium bromide (inhibitor of organellar DNA synthesis). In both cases germination and the early processes of post-germinative growth, including enlargement and greening of cotyledons, some elongation of radicle, and formation of root hairs, occurred at the same rate as in the controls. Comparable results were obtained using fluorodeoxyuridine. The effect of aphidicolin is, unlike fluorodeoxyuridine, completely reversible and the side effect on protein and RNA synthesis are negligible. Total non-dependence on nuclear and organellar DNA synthesis during early post-germinative growth may not be a generalized situation in all seeds, as was shown by some observations on lettuce and watermelon.     

Modulation of DNA end joining by nuclear proteins

DNA double strand breaks in mammalian cells are primarily repaired by homologous recombination and non-homologous end joining (NHEJ). NHEJ may either be error-free or mutagenic with deletions or insertions at the joint. Recent studies showed that DNA ends can also be joined via microhomologous sequences flanking the break point especially when proteins responsible for NHEJ, such as Ku, are absent. Microhomology-mediated end joining (MHEJ) is always accompanied by a deletion that spans one of the two homologous sequences and the intervening sequence, if any. In this study we evaluated several factors affecting the relative contribution of MHEJ to DNA end joining using nuclear extracts and DNA substrates containing 10-bp repeats at the ends. We found that the occurrence of MHEJ is determined by the relative abundance of nuclear proteins. At low DNA/protein ratios, an error-free end-joining mechanism predominated over MHEJ. As the DNA/protein ratio increased, MHEJ became predominant. We show that the nuclear proteins that contribute to the inhibition of the error-prone MHEJ include Ku and histone H1. Treatment of extracts with flap endonuclease 1 antiserum significantly reduced MHEJ. Addition of a 17-bp intervening sequence between the microhomologous sequences significantly reduced the efficiency of MHEJ. Thus, this cell-free assay provides a platform for evaluating factors modulating end joining.     

Using partial genomic fosmid libraries for sequencing complete organellar genomes

Organellar genome sequences provide numerous phylogenetic markers and yield insight into organellar function and molecular evolution. These genomes are much smaller in size than their nuclear counterparts; thus, their complete sequencing is much less expensive than total nuclear genome sequencing, making broader phylogenetic sampling feasible. However; for some organisms, it is challenging to isolate plastid DNA for sequencing using standard methods. To overcome these difficulties, we constructed partial genomic libraries from total DNA preparations of two heterotrophic and two autotrophic angiosperm species using fosmid vectors. We then used macroarray screening to isolate clones containing large fragments of plastid DNA. A minimum tiling path of clones comprising the entire genome sequence of each plastid was selected, and these clones were shotgun-sequenced and assembled into complete genomes. Although this method worked well for both heterotrophic and autotrophic plants, nuclear genome size had a dramatic effect on the proportion of screened clones containing plastid DNA and, consequently, the overall number of clones that must be screened to ensure full plastid genome coverage. This technique makes it possible to determine complete plastid genome sequences for organisms that defy other available organellar genome sequencing methods, especially those for which limited amounts of tissue are available.     

Organellar genome,nuclear ribosomal DNA repeat unit,and microsatellites isolated from a small-scale of 454 GS FLX sequencing on two mosses

Recent innovations in high-throughput DNA sequencing methodology (next generation sequencing technologies [NGS]) allow for the generation of large amounts of high quality data that may be particularly critical for resolving ambiguous relationships such as those resulting from rapid radiations. Application of NGS technology to bryology is limited to assembling entire nuclear or organellar genomes of selected exemplars of major lineages (e.g., classes). Here we outline how organellar genomes and the entire nuclear ribosomal DNA repeat can be obtained from minimal amounts of moss tissue via small-scale 454 GS FLX sequencing. We sampled two Funariaceae species, Funaria hygrometrica and Entosthodon obtusus, and assembled nearly complete organellar genomes and the whole nuclear ribosomal DNA repeat unit (18S-ITS1-5.8S-ITS2-26S-IGS1-5S-IGS2) for both taxa. Sequence data from these species were compared to sequences from another Funariaceae species, Physcomitrella patens, revealing low overall degrees of divergence of the organellar genomes and nrDNA genes with substitutions spread rather evenly across their length, and high divergence within the external spacers of the nrDNA repeat. Furthermore, we detected numerous microsatellites among the 454 assemblies. This study demonstrates that NGS methodology can be applied to mosses to target large genomic regions and identify microsatellites.     

Isolation and analysis of recombinant DNA molecules containing yeast DNA.

2500 recombinant plasmids containing insertions of yeast nuclear DNA have been cloned in Escherichia coli. It can be calculated that about 85% of the yeast genome is represented in this collection. The clones have been characterized by hybridization to purified RNA species. Of the 2000 clones examined, 75 contain insertions of yeast ribosomal DNA, 201 contain insertions of yeast tRNA genes, and 26 contain DNA sequences that are complementary to abundant mRNA species.     

Analysis of nuclear copies of mitochondrial sequences in honeybee (Apis mellifera) genome

At least 0.08% of the Apis mellifera nuclear genome contains sequences that originated from mitochondria. These nuclear copies of mitochondrial sequences (numts) are scattered all over the honeybee chromosomes and have originated by multiple independent insertions of mitochondrial DNA (mtDNA) as evident by phylogenetic analysis. Apart from original insertions, moderate duplications of numts also contributed to the present pattern and distribution of mitochondrial sequences in honeybee chromosomes. Assimilation of mitochondrial genes in the nuclear genome is mediated by extensive fragmentations of the original inserts. Replication slippage seems to be a major mechanism by which small sequences are inserted or deleted from mtDNA destined to nucleus. Most of the honeybee numts (84%) are located in the nongenic regions. The majority (94%) of the numts that are located in predicted nuclear genes have originated from mitochondrial genes coding for cytochrome oxidase and NADH dehydrogenase subunits. On the other hand, the mitochondrial rRNA or tRNA gene sequences are predominantly (88%) located in nongenic regions of the genome. Evidences also support for exertion of purifying selection on numts located in specific genes. Comparative analysis of numts of European, African, and Africanized honeybees suggests that numt evolution in A. mellifera is probably not demarked by speciation time frame but may be a continuous and dynamic process.