Continued colonization of the human genome by mitochondrial DNA

Integration of mitochondrial DNA fragments into nuclear chromosomes (giving rise to nuclear DNA sequences of mitochondrial origin, or NUMTs) is an ongoing process that shapes nuclear genomes. In yeast this process depends on double-strand-break repair. Since NUMTs lack amplification and specific integration mechanisms, they represent the prototype of exogenous insertions in the nucleus. From sequence analysis of the genome of Homo sapiens, followed by sampling humans from different ethnic backgrounds, and chimpanzees, we have identified 27 NUMTs that are specific to humans and must have colonized human chromosomes in the last 4–6 million years. Thus, we measured the fixation rate of NUMTs in the human genome. Six such NUMTs show insertion polymorphism and provide a useful set of DNA markers for human population genetics. We also found that during recent human evolution, Chromosomes 18 and Y have been more susceptible to colonization by NUMTs. Surprisingly, 23 out of 27 human-specific NUMTs are inserted in known or predicted genes, mainly in introns. Some individuals carry a NUMT insertion in a tumor-suppressor gene and in a putative angiogenesis inhibitor. Therefore in humans, but not in yeast, NUMT integrations preferentially target coding or regulatory sequences. This is indeed the case for novel insertions associated with human diseases and those driven by environmental insults. We thus propose a mutagenic phenomenon that may be responsible for a variety of genetic diseases in humans and suggest that genetic or environmental factors that increase the frequency of chromosome breaks provide the impetus for the continued colonization of the human genome by mitochondrial DNA.     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Allelic genome structural variations in maize detected by array comparative genome hybridization

DNA polymorphisms such as insertion/deletions and duplications affecting genome segments larger than 1 kb are known as copy-number variations (CNVs) or structural variations (SVs). They have been recently studied in animals and humans by using array-comparative genome hybridization (aCGH), and have been associated with several human diseases. Their presence and phenotypic effects in plants have not been investigated on a genomic scale, although individual structural variations affecting traits have been described. We used aCGH to investigate the presence of CNVs in maize by comparing the genome of 13 maize inbred lines to B73. Analysis of hybridization signal ratios of 60,472 60-mer oligonucleotide probes between inbreds in relation to their location in the reference genome (B73) allowed us to identify clusters of probes that deviated from the ratio expected for equal copy-numbers. We found CNVs distributed along the maize genome in all chromosome arms. They occur with appreciable frequency in different germplasm subgroups, suggesting ancient origin. Validation of several CNV regions showed both insertion/deletions and copy-number differences. The nature of CNVs detected suggests CNVs might have a considerable impact on plant phenotypes, including disease response and heterosis.     

Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies

Nuclear insertions of mitochondrial origin (NUMTs) can be useful tools in evolution and population studies. However, due to their similarity to mitochondrial DNA (mtDNA), NUMTs may also be a source of contamination in mtDNA studies. The main goal of this work is to present a database of NUMTs, based on the latest version of the human genome—GRCh37 draft. A total of 755 insertions were identified. There are 33 paralogous sequences with over 80% sequence similarity and of a greater length than 500 bp. The non-identical positions between paralogous sequences are listed for the first time. As an application example, the described database is used to evaluate the impact of NUMT contamination in cancer studies. The evaluation reveals that 220 positions from 256 with zero hits in the current mtDNA phylogeny could in fact be traced to one or more nuclear insertions of mtDNA. This is due to they are located in non-identical positions between mtDNA and nuclear DNA (nDNA). After in silico primer validation of each revised cancer study, risk of co-amplification between mtDNA and nDNA was detected in some cases, whereas in others no risk of amplification was identified. This approach to cancer studies clearly proves the potential of our NUMT database as a valuable new tool to validate mtDNA mutations described in different contexts. Moreover, due to the amount of information provided for each nuclear insertion, this database should play an important role in designing evolutionary, phylogenetic and epidemiological studies.     

Intraspecific variation of recombination rate in maize

Background

In sexually reproducing organisms, meiotic crossovers ensure the proper segregation of chromosomes and contribute to genetic diversity by shuffling allelic combinations. Such genetic reassortment is exploited in breeding to combine favorable alleles, and in genetic research to identify genetic factors underlying traits of interest via linkage or association-based approaches. Crossover numbers and distributions along chromosomes vary between species, but little is known about their intraspecies variation.

Results

Here, we report on the variation of recombination rates between 22 European maize inbred lines that belong to the Dent and Flint gene pools. We genotype 23 doubled-haploid populations derived from crosses between these lines with a 50 k-SNP array and construct high-density genetic maps, showing good correspondence with the maize B73 genome sequence assembly. By aligning each genetic map to the B73 sequence, we obtain the recombination rates along chromosomes specific to each population. We identify significant differences in recombination rates at the genome-wide, chromosome, and intrachromosomal levels between populations, as well as significant variation for genome-wide recombination rates among maize lines. Crossover interference analysis using a two-pathway modeling framework reveals a negative association between recombination rate and interference strength.

Conclusions

To our knowledge, the present work provides the most comprehensive study on intraspecific variation of recombination rates and crossover interference strength in eukaryotes. Differences found in recombination rates will allow for selection of high or low recombining lines in crossing programs. Our methodology should pave the way for precise identification of genes controlling recombination rates in maize and other organisms.     

Maize chromosomal knobs are located in gene-dense areas and suppress local recombination

Knobs are conspicuous heterochromatic regions found on the chromosomes of maize and its relatives. The number, locations, and sizes of knobs vary dramatically, with most lines containing between four and eight knobs in mid-arm positions. Prior data suggest that some knobs may reduce recombination. However, comprehensive tests have not been carried out, primarily because most knobs have not been placed on the genetic map. We used fluorescent in situ hybridization and two recombinant inbred populations to map seven knobs and to accurately place three knobs from the B73 inbred on the genomic sequence assembly. The data show that knobs lie in gene-dense regions of the maize genome. Comparisons to 23 other recombinant inbred populations segregating for knobs at the same sites confirm that large knobs can locally reduce crossing over by as much as twofold on a cM/Mb scale. These effects do not extend beyond regions ~10 cM to either side of knobs and do not appear to affect linkage disequilibrium among genes within and near knob repeat regions of the B73 RefGen_v2 assembly.     

Nuclear mitochondrial DNA activates replication in Saccharomyces cerevisiae

The nuclear genome of eukaryotes is colonized by DNA fragments of mitochondrial origin, called NUMTs. These insertions have been associated with a variety of germ-line diseases in humans. The significance of this uptake of potentially dangerous sequences into the nuclear genome is unclear. Here we provide functional evidence that sequences of mitochondrial origin promote nuclear DNA replication in Saccharomyces cerevisiae. We show that NUMTs are rich in key autonomously replicating sequence (ARS) consensus motifs, whose mutation results in the reduction or loss of DNA replication activity. Furthermore, 2D-gel analysis of the mrc1 mutant exposed to hydroxyurea shows that several NUMTs function as late chromosomal origins. We also show that NUMTs located close to or within ARS provide key sequence elements for replication. Thus NUMTs can act as independent origins, when inserted in an appropriate genomic context or affect the efficiency of pre-existing origins. These findings show that migratory mitochondrial DNAs can impact on the replication of the nuclear region they are inserted in.     

Genetic mapping shows intraspecific variation and transgressive segregation for caterpillar‐induced aphid resistance in maize

Plants in nature have inducible defences that sometimes lead to targeted resistance against particular herbivores, but susceptibility to others. The metabolic diversity and genetic resources available for maize (Zea mays) make this a suitable system for a mechanistic study of within‐species variation in such plant‐mediated interactions between herbivores. Beet armyworms (Spodoptera exigua) and corn leaf aphids (Rhopalosiphum maidis) are two naturally occurring maize herbivores with different feeding habits. Whereas chewing herbivore‐induced methylation of 2,4‐dihydroxy‐7‐methoxy‐1,4‐benzoxazin‐3‐one glucoside (DIMBOA‐Glc) to form 2‐hydroxy‐4,7‐dimethoxy‐1,4‐benzoxazin‐3‐one glucoside (HDMBOA‐Glc) promotes caterpillar resistance, lower DIMBOA‐Glc levels favour aphid reproduction. Thus, caterpillar‐induced DIMBOA‐Glc methyltransferase activity in maize is predicted to promote aphid growth. To test this hypothesis, the impact of S. exigua feeding on R. maidis progeny production was assessed using seventeen genetically diverse maize inbred lines. Whereas aphid progeny production was increased by prior caterpillar feeding on lines B73, Ki11, Ki3 and Tx303, it decreased on lines Ky21, CML103, Mo18W and W22. Genetic mapping of this trait in a population of B73 × Ky21 recombinant inbred lines identified significant quantitative trait loci on maize chromosomes 1, 7 and 10. There is a transgressive segregation for aphid resistance, with the Ky21 alleles on chromosomes 1 and 7 and the B73 allele on chromosome 10 increasing aphid progeny production. The chromosome 1 QTL coincides with a cluster of three maize genes encoding benzoxazinoid O‐methyltransferases that convert DIMBOA‐Glc to HDMBOA‐Glc. Gene expression studies and benzoxazinoid measurements indicate that S. exigua ‐induced responses in this pathway differentially affect R. maidis resistance in B73 and Ky21.     

Multiple chromosomes in bacteria: structure and function of chromosome II of Rhodobacter sphaeroides 2.4.1T.

Although multiple chromosomes occur in bacteria, much remains to be learned about their structural and functional interrelationships. To study the structure-function relationships of chromosomes I and II of the facultative photosynthetic bacterium Rhodobacter sphaeroides 2.4.1T, auxotrophic mutants were isolated. Five strains having transposon insertions in chromosome II showed requirements for p-aminobenzoic acid (pABA)-dihydroxybenzoic acid (dHBA), serine, thymine, uracil, or histidine. The His, Thy, and pABA-dHBA mutants reverted to prototrophy at low frequency and concordantly lost their transposon insertions from the genome. The Ser, Ura, and pABA-dHBA mutants were complemented by cosmids that carried the region of chromosome II where the transposon insertions were located. The cosmids used for complementation analysis were selected, on the basis of map position, from a set of overlapping clones that had been ordered by a combination of hybridization and restriction endonuclease mapping. These experiments provide the basis for detailed studies of the structure, function, and interaction between each chromosome, and they demonstrate at this early stage of investigation that no fundamental differences exist between each chromosome.     

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.     

Population Polymorphism of Nuclear Mitochondrial DNA Insertions Reveals Widespread Diploidy Associated with Loss of Heterozygosity in Debaryomyces hansenii

Debaryomyces hansenii, a yeast that participates in the elaboration of foodstuff, displays important genetic diversity. Our recent phylogenetic classification of this species led to the subdivision of the species into three distinct clades. D. hansenii harbors the highest number of nuclear mitochondrial DNA (NUMT) insertions known so far for hemiascomycetous yeasts. Here we assessed the intraspecific variability of the NUMTs in this species by testing their presence/absence first in 28 strains, with 21 loci previously detected in the completely sequenced strain CBS 767^T, and second in a larger panel of 77 strains, with 8 most informative loci. We were able for the first time to structure populations in D. hansenii, although we observed little NUMT insertion variability within the clades. We determined the chronology of the NUMT insertions, which turned out to correlate with the previously defined taxonomy and provided additional evidence that colonization of nuclear genomes by mitochondrial DNA is a dynamic process in yeast. In combination with flow cytometry experiments, the NUMT analysis revealed the existence of both haploid and diploid strains, the latter being heterozygous and resulting from at least four crosses among strains from the various clades. As in the diploid pathogen Candida albicans, to which D. hansenii is phylogenetically related, we observed a differential loss of heterozygosity in the diploid strains, which can explain some of the large genetic diversity found in D. hansenii over the years.Debaryomyces hansenii is a ubiquist, hemiascomycetous yeast that can be found in soil, fruits, and various manufactured foodstuff in which it participates by contributing to the maturation or as a contaminant. Its ability to grow at low temperatures and in high salinity environments makes it the most common yeast in cheeses, to which it brings a number of proteolytic and lipolytic activities and aromas in the course of maturation. D. hansenii has also been implicated as an emerging pathogen, sometimes under the name of Candida famata var. famata (see reference 17). Taxonomic classification of the species related to D. hansenii has always been subject to debate. Recent analyses have reinstated D. hansenii (previously D. hansenii var. hansenii), Debaryomyces fabryi (previously D. hansenii var. fabryi), and Debaryomyces subglobosus (previously Candida famata var. flareri) (13, 25). Phylogenetic analysis using conserved spliceosomal intron sequence comparison has shown that D. hansenii is a complex of species, which comprises at least four members: D. hansenii, Debaryomyces tyrocola, D. fabryi, and Candida flareri (previously Candida famata var. flareri) (18). In addition, our study has revealed the existence of at least three populations (clades 1 to 3) in D. hansenii, with the first one containing the strain CBS 767^T, which has been entirely sequenced (8), and the last one containing Candida famata var. famata CBS 1795.Most eukaryotic nuclear genomes contain pieces of mitochondrial sequences (designated NUMT [nuclear mitochondrial DNA] for nuclear sequences of mitochondrial origin) that result from the transfer of fragments of mitochondrial DNA (mtDNA) to the chromosomes. The number and size of the NUMTs varies greatly between eukaryotic genomes (33). A recent investigation of six hemiascomycetous yeasts has shown that even within this monophyletic group, the number of NUMTs varies greatly, from 1 in Kluyveromyces thermotolerans CBS 6340^T to 145 in D. hansenii CBS 767^T (36). The mtDNA is thought to invade nuclear genomes during the repair of chromosomal DNA double-strand breaks (DSB) by nonhomologous end joining (NHEJ), as shown experimentally in the yeast Saccharomyces cerevisiae (31, 44). The colonization of nuclear genomes by mtDNA is a dynamic evolutionary process, as observed in yeast and humans (3, 32).D. hansenii harbors the highest number of NUMTs known so far for hemiascomycetous yeasts, making it of particular interest for NUMT studies. Conversely, NUMTs are potentially interesting markers to differentiate strains of this species. The 145 NUMTs of type strain CBS 767^T are distributed in 86 loci (61 single NUMTs and 25 clusters). Most clusters (23, 25) are mosaics of NUMTs formed from noncontiguous mtDNA fragments inserted in random orientation at the same chromosomal locus. In the other two clusters, the NUMTs are all in the same orientation and order, as in the mitochondrial genome. These clusters (designated “processions”) correspond to a single ancient mtDNA insertion, followed by mutational decay, leaving recognizable mtDNA segments separated by more diverged sequences (36).Few studies have attempted to evaluate the variability of NUMTs within the same species (2, 23, 32). Here, we have studied natural isolates to assess the intraspecific variability of the NUMT insertions in the nuclear genome of the yeast species D. hansenii. We were able to structure populations in this species, to determine the chronology of the NUMT insertions, and to correlate this chronology to the taxonomy of the D. hansenii complex species. Moreover, NUMT analysis revealed the existence of both haploid and diploid strains, the latter resulting from crosses between different D. hansenii clades.     

Microsatellite megatracts in the maize (Zea mays L.) genome.

Long tracts (megatracts) of (CAG)n, (TAG)n, and (GAA)n microsatellite sequences capable of forming composite DNA segments were found in the maize (Zea mays L.) genome. Some of the (CAG)n and (TAG)n megatracts were organized in clusters of up to 1 Mb on several chromosomes, as detected by fluorescence in situ hybridization (FISH), as well as on extended DNA fibers. Extensive polymorphism was found among different maize inbred lines with respect to the number and size of microsatellite megatract clusters on the A chromosomes. Polymorphism was also common among B chromosomes of different nuclei in the inbred line Zapalote Chico. Different retrotransposable elements were often inserted into the microsatellite tracts. Size variation in some (TAG)n and (GAA)n megatracts was observed in consecutive generations among siblings of the inbred lines, indicating that these loci are highly unstable and predisposed to dynamic mutations similar to those described in mammalian systems.     

Analysis of a 120-Kilobase Mitochondrial Chromosome in Maize

The organization of the mitochondrial genome in plants is not well understood. In maize mitochondrial DNA (mtDNA) several subgenomic circular molecules as well as an abundant fraction of linear molecules have been seen by electron microscopy. It has been hypothesized that the circular molecules are the genetic entities of the mitochondrial genome while the linear molecules correspond to randomly sheared mtDNA. A model has been proposed that explains the mechanism of generation of subgenomic circles (of a predictable size) by homologous recombination between pairs of large direct repeats found on a large (approximately 570 kb for the fertile (N) cytoplasm) master circle. So far the physical entities of the mitochondrial genome, as they exist in vivo, and the genes they carry, have not been identified. For this purpose, we used two gel systems (pulsed field gel electrophoresis and Eckhardt gels) designed to resolve large DNA. Large DNA was prepared from the Black Mexican Sweet (BMS) cultivar. We resolved several size classes of mtDNA circles and designate these as chromosomes. A 120 kb chromosome was mapped in detail. It is shown to contain the three ribosomal genes (rrn26, rrn18 and rrn5) plus two genes encoding subunits of cytochrome oxidase (Cox1 and Cox3); it appears to be colinear with the 570-kb master circle map of another fertile cytoplasm (B37N) except at the "breakpoints" required to form the 120-kb circle. The presence of the 120-kb chromosome could not have been predicted by homologous recombination through any of the known repetitive sequences nor is it a universal feature of normal maize mitochondria. It is present in mitochondria of BMS suspension cultures and seedlings, but is not detectable in seedlings of B37N. No master genome was detected in BMS.     

Promiscuous DNA in the nuclear genomes of hemiascomycetous yeasts

Transfer of fragments of mtDNA to the nuclear genome is a general phenomenon that gives rise to NUMTs (NUclear sequences of MiTochondrial origin). We present here the first comparative analysis of the NUMT content of entirely sequenced species belonging to a monophyletic group, the hemiascomycetous yeasts ( Candida glabrata, Kluyveromyces lactis, Kluyveromyces thermotolerans, Debaryomyces hansenii and Yarrowia lipolytica , along with the updated NUMT content of Saccharomyces cerevisiae ). This study revealed a huge diversity in NUMT number and organization across the six species. Debaryomyces hansenii harbors the highest number of NUMTs (145), half of which are distributed in numerous large mosaics of up to eight NUMTs arising from multiple noncontiguous mtDNA fragments inserted at the same chromosomal locus. Most NUMTs, in all species, are found within intergenic regions including seven NUMTs in pseudogenes. However, five NUMTs overlap a gene, suggesting a positive impact of NUMTs on protein evolution. Contrary to the other species, K. lactis and K. thermotolerans harbor only a few diverged NUMTs, suggesting that mitochondrial transfer to the nuclear genome has decreased or ceased in these phylogenetic branches. The dynamics of NUMT acquisition and loss are illustrated here by their species-specific distribution.     

Analysis of plastid and mitochondrial DNA insertions in the nucleus (NUPTs and NUMTs) of six plant species: size,relative age and chromosomal localization

We analysed the size, relative age and chromosomal localization of nuclear sequences of plastid and mitochondrial origin (NUPTs-nuclear plastid DNA and NUMTs-nuclear mitochondrial DNA) in six completely sequenced plant species. We found that the largest insertions showed lower divergence from organelle DNA than shorter insertions in all species, indicating their recent origin. The largest NUPT and NUMT insertions were localized in the vicinity of the centromeres in the small genomes of Arabidopsis and rice. They were also present in other chromosomal regions in the large genomes of soybean and maize. Localization of NUPTs and NUMTs correlated positively with distribution of transposable elements (TEs) in Arabidopsis and sorghum, negatively in grapevine and soybean, and did not correlate in rice or maize. We propose a model where new plastid and mitochondrial DNA sequences are inserted close to centromeres and are later fragmented by TE insertions and reshuffled away from the centromere or removed by ectopic recombination. The mode and tempo of TE dynamism determines the turnover of NUPTs and NUMTs resulting in their species-specific chromosomal distributions.     

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.