Structure and chromosomal distribution of human mitochondrial pseudogenes

Unconstrained precision mitochondrial genome editing with αDdCBEs

DddA-derived cytosine base editors (DdCBEs) enable the targeted introduction of C•G-to-T•A conversions in mitochondrial DNA (mtDNA). DdCBEs work in pairs, with each arm composed of a transcription activator-like effector (TALE), a split double-stranded DNA deaminase half, and a uracil glycosylase inhibitor. This pioneering technology has helped improve our understanding of cellular processes involving mtDNA and has paved the way for the development of models and therapies for genetic disorders caused by pathogenic mtDNA variants. Nonetheless, given the intrinsic properties of TALE proteins, several target sites in human mtDNA are predicted to remain out of reach to DdCBEs and other TALE-based technologies. Specifically, due to the conventional requirement for a thymine immediately upstream of the TALE target sequences (i.e., the 5'-T constraint), over 150 loci in the human mitochondrial genome are presumed to be inaccessible to DdCBEs. Previous attempts at circumventing this requirement, either by developing monomeric DdCBEs or utilizing DNA-binding domains alternative to TALEs, have resulted in suboptimal specificity profiles with reduced therapeutic potential. Here, aiming to challenge and elucidate the relevance of the 5'-T constraint in the context of DdCBE-mediated mtDNA editing, and to expand the range of motifs that are editable by this technology, we generated DdCBEs containing TALE proteins engineered to recognize all 5' bases. These modified DdCBEs are herein referred to as αDdCBEs. Notably, 5'-T-noncompliant canonical DdCBEs efficiently edited mtDNA at diverse loci. However, they were frequently outperformed by αDdCBEs, which exhibited significant improvements in activity and specificity, regardless of the most 5' bases of their TALE binding sites. Furthermore, we showed that αDdCBEs are compatible with the enhanced DddAtox variants DddA6 and DddA11, and we validated TALE shifting with αDdCBEs as an effective approach to optimize base editing outcomes. Overall, αDdCBEs enable efficient, specific, and unconstrained mitochondrial base editing.

Unconstrained Precision Mitochondrial Genome Editing with αDdCBEs

DddA-derived cytosine base editors (DdCBEs) enable the targeted introduction of C•G-to-T•A conversions in mitochondrial DNA (mtDNA). DdCBEs are often deployed as pairs, with each arm comprised of a transcription activator-like effector (TALE), a split double-stranded DNA deaminase half, and a uracil glycosylase inhibitor. This pioneering technology has helped improve our understanding of cellular processes involving mtDNA and has paved the way for the development of models and therapies for genetic disorders caused by pathogenic mtDNA variants. Nonetheless, given the intrinsic properties of TALE proteins, several target sites in human mtDNA remain out of reach to DdCBEs and other TALE-based technologies. Specifically, due to the conventional requirement for a thymine immediately upstream of the TALE target sequences (i.e., the 5'-T constraint), over 150 loci in the human mitochondrial genome are presumed to be inaccessible to DdCBEs. Previous attempts at circumventing this constraint, either by developing monomeric DdCBEs or utilizing DNA-binding domains alternative to TALEs, have resulted in suboptimal specificity profiles with reduced therapeutic potential. Here, aiming to challenge and elucidate the relevance of the 5'-T constraint in the context of DdCBE-mediated mtDNA editing, and to expand the range of motifs that are editable by this technology, we generated αDdCBEs that contain modified TALE proteins engineered to recognize all 5' bases. Notably, 5'-T-noncompliant, canonical DdCBEs efficiently edited mtDNA at diverse loci. However, DdCBEs were frequently outperformed by αDdCBEs, which consistently displayed significant improvements in activity and specificity, regardless of the 5'-most bases of their TALE binding sites. Furthermore, we showed that αDdCBEs are compatible with DddA tox and its derivatives DddA6, and DddA11, and we validated TALE shifting with αDdCBEs as an effective approach to optimize base editing outcomes at a single target site. Overall, αDdCBEs enable efficient, specific, and unconstrained mitochondrial base editing.

Interrogating 1000 insect genomes for NUMTs: A risk assessment for estimates of species richness

The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Applying genomic and transcriptomic advances to mitochondrial medicine

Next-generation sequencing (NGS) has increased our understanding of the molecular basis of many primary mitochondrial diseases (PMDs). Despite this progress, many patients with suspected PMD remain without a genetic diagnosis, which restricts their access to in-depth genetic counselling, reproductive options and clinical trials, in addition to hampering efforts to understand the underlying disease mechanisms. Although they represent a considerable improvement over their predecessors, current methods for sequencing the mitochondrial and nuclear genomes have important limitations, and molecular diagnostic techniques are often manual and time consuming. However, recent advances in genomics and transcriptomics offer realistic solutions to these challenges. In this Review, we discuss the current genetic testing approach for PMDs and the opportunities that exist for increased use of whole-genome NGS of nuclear and mitochondrial DNA (mtDNA) in the clinical environment. We consider the possible role for long-read approaches in sequencing of mtDNA and in the identification of novel nuclear genomic causes of PMDs. We examine the expanding applications of RNA sequencing, including the detection of cryptic variants that affect splicing and gene expression and the interpretation of rare and novel mitochondrial transfer RNA variants.

Mitochondrial DNA insert into CD40 ligand gene-associated X-linked hyper-IgM syndrome

BACKGROUND

X-linked hyper-IgM (X-HIGM), which results from mutations in the CD40LG gene located on chromosome Xq26.3, is the most common form of HIGM. To date, more than 130 variants of the CD40L gene have been reported. We described a patient with novel de novo nuclear mitochondrial DNA sequences (NUMTs) in the CD40LG gene that have resulted in X-HIGM.

METHODS

Whole-exome sequencing (WES) analysis was used to screen for causal variants in the genome, and the candidate breakpoint was confirmed by Sanger sequencing.

RESULTS

A new mutation of CD40LG, which deletes A at position 17 followed by a 147-nucleotide from mitochondrial DNA copies insertion in exon 1, was detected in a 20-month-old boy harbouring an X-HIGM combined with immunodeficiency syndrome.

CONCLUSION

This is one of the few cases of a human genetic disease caused by nuclear mitochondrial DNA sequences (NUMTs). The presented data serve to demonstrate that de novo NUMT transfer of nucleic acid is a novel mechanism of X-HIGM.

Mitochondria: their role in spermatozoa and in male infertility

BACKGROUND

The best-known role of spermatozoa is to fertilize the oocyte and to transmit the paternal genome to offspring. These highly specialized cells have a unique structure consisting of all the elements absolutely necessary to each stage of fertilization and to embryonic development. Mature spermatozoa are made up of a head with the nucleus, a neck, and a flagellum that allows motility and that contains a midpiece with a mitochondrial helix. Mitochondria are central to cellular energy production but they also have various other functions. Although mitochondria are recognized as essential to spermatozoa, their exact pathophysiological role and their functioning are complex. Available literature relative to mitochondria in spermatozoa is dense and contradictory in some cases. Furthermore, mitochondria are only indirectly involved in cytoplasmic heredity as their DNA, the paternal mitochondrial DNA, is not transmitted to descendants.

OBJECTIVE AND RATIONAL

This review aims to summarize available literature on mitochondria in spermatozoa, and, in particular, that with respect to humans, with the perspective of better understanding the anomalies that could be implicated in male infertility.

SEARCH METHODS

PubMed was used to search the MEDLINE database for peer-reviewed original articles and reviews pertaining to human spermatozoa and mitochondria. Searches were performed using keywords belonging to three groups: 'mitochondria' or 'mitochondrial DNA', 'spermatozoa' or 'sperm' and 'reactive oxygen species' or 'calcium' or 'apoptosis' or signaling pathways'. These keywords were combined with other relevant search phrases. References from these articles were used to obtain additional articles.

OUTCOMES

Mitochondria are central to the metabolism of spermatozoa and they are implicated in energy production, redox equilibrium and calcium regulation, as well as apoptotic pathways, all of which are necessary for flagellar motility, capacitation, acrosome reaction and gametic fusion. In numerous cases, alterations in one of the aforementioned functions could be linked to a decline in sperm quality and/or infertility. The link between the mitochondrial genome and the quality of spermatozoa appears to be more complex. Although the quantity of mtDNA, and the existence of large-scale deletions therein, are inversely correlated to sperm quality, the effects of mutations seem to be heterogeneous and particularly related to their pathogenicity.

WIDER IMPLICATIONS

The importance of the role of mitochondria in reproduction, and particularly in gamete quality, has recently emerged following numerous publications. Better understanding of male infertility is of great interest in the current context where a significant decline in sperm quality has been observed.

Intronic tRNAs of mitochondrial origin regulate constitutive and alternative splicing

BACKGROUND

The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive.

RESULTS

Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNALys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3'-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNATyr by an electrophoretic mobility shift assay.

CONCLUSIONS

We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

Characterization of nuclear mitochondrial insertions in the whole genomes of primates

The transfer and integration of whole and partial mitochondrial genomes into the nuclear genomes of eukaryotes is an ongoing process that has facilitated the transfer of genes and contributed to the evolution of various cellular pathways. Many previous studies have explored the impact of these insertions, referred to as NumtS, but have focused primarily on older events that have become fixed and are therefore present in all individual genomes for a given species. We previously developed an approach to identify novel Numt polymorphisms from next-generation sequence data and applied it to thousands of human genomes. Here, we extend this analysis to 79 individuals of other great ape species including chimpanzee, bonobo, gorilla, orang-utan and also an old world monkey, macaque. We show that recent Numt insertions are prevalent in each species though at different apparent rates, with chimpanzees exhibiting a significant increase in both polymorphic and fixed Numt sequences as compared to other great apes. We further assessed positional effects in each species in terms of evolutionary time and rate of insertion and identified putative hotspots on chromosome 5 for Numt integration, providing insight into both recent polymorphic and older fixed reference NumtS in great apes in comparison to human events.

Identification of Mitochondrial DNA (NUMTs) in the Nuclear Genome of Daphnia magna

This is the first study in which the Daphnia magna (D. magna) nuclear genome (nDNA) obtained from the GenBank database was analyzed for pseudogene sequences of mitochondrial origin. To date, there is no information about pseudogenes localized in D. magna genome. This study aimed to identify NUMTs, their length, homology, and location for potential use in evolutionary studies and to check whether their occurrence causes co-amplification during mitochondrial genome (mtDNA) analyses. Bioinformatic analysis showed 1909 fragments of the mtDNA of D. magna, of which 1630 were located in ten linkage groups (LG) of the nDNA. The best-matched NUMTs covering >90% of the gene sequence have been identified for two mt-tRNA genes, and they may be functional nuclear RNA molecules. Isolating the total DNA in mtDNA studies, co-amplification of nDNA fragments is unlikely in the case of amplification of the whole tRNA genes as well as fragments of other genes. It was observed that TRNA-MET fragments had the highest level of sequence homology, thus they could be evolutionarily the youngest. The lowest homology was found in the D-loop-derived pseudogene. It may probably be the oldest NUMT incorporated into the nDNA; however, further analysis is necessary.

Mitochondrial DNA Parameters in Blood of Infants Receiving Lopinavir/Ritonavir or Lamivudine Prophylaxis to Prevent Breastfeeding Transmission of HIV-1

Children who are human immunodeficiency virus (HIV)-exposed but uninfected (CHEU) accumulate maternal HIV and antiretroviral exposures through pregnancy, postnatal prophylaxis, and breastfeeding. Here, we compared the dynamics of mitochondrial DNA (mtDNA) parameters in African breastfed CHEU receiving lopinavir/ritonavir (LPV/r) or lamivudine (3TC) pre-exposure prophylaxis during the first year of life. The number of mtDNA copies per cell (MCN) and the proportion of deleted mtDNA (MDD) were assessed at day 7 and at week 50 post-delivery (PrEP group). mtDNA depletion was defined as a 50% or more decrease from the initial value, and mtDNA deletions was the detection of mtDNA molecules with large DNA fragment loss. We also performed a sub-analysis with CHEU who did not receive a prophylactic treatment in South Africa (control group). From day seven to week 50, MCN decreased with a median of 41.7% (interquartile range, IQR: 12.1; 64.4) in the PrEP group. The proportion of children with mtDNA depletion was not significantly different between the two prophylactic regimens. Poisson regressions showed that LPV/r and 3TC were associated with mtDNA depletion (reference: control group; LPV/r: PR = 1.75 (CI95%: 1.15-2.68), p < 0.01; 3TC: PR = 1.54 (CI95%: 1.00-2.37), p = 0.05). Moreover, the proportion of children with MDD was unexpectedly high before randomisation in both groups. Long-term health impacts of these mitochondrial DNA parameters should be investigated further for both CHEU and HIV-infected children receiving LPV/r- or 3TC- based regimens.

Survey of mitochondrial sequences integrated into the bovine nuclear genome

Nuclear copies of the mitochondrial DNA (NUMTs) have already been described in several species. In this context, we identified and analysed 166 bovine NUMT regions with a total length of 430 kbp, representing about 0.02% of the cattle nuclear genome. Copies of all mitochondrial regions were detected in the nuclear genome, with distinct degrees of sequence similarity to the mitogenome. Some NUMT regions include large mitogenome segments and show high similarity to the organelle genome sequence. NUMT regions are frequently modified by insertion of repetitive sequences and by sequence rearrangements. We confirmed the existence of 29 NUMT regions by PCR amplification using DNA from the cow (Dominette) which was used to generate the bovine genome reference sequence, ruling out the possibility that these NUMTs could be artifacts of the genome assembly. As there are NUMT regions with high similarity to the mitogenome, special care is needed when designing primers for mitochondrial DNA amplification. Our results can therefore be used to avoid co-amplification of bovine nuclear sequences similar to mitochondrial DNA.

A simple method for sequencing the whole human mitochondrial genome directly from samples and its application to genetic testing

Next-generation sequencing (NGS) is a revolutionary sequencing technology for analyzing genomes. However, preprocessing methods for mitochondrial DNA (mtDNA) sequencing remain complex, and it is required to develop an authenticated preprocessing method. Here, we developed a simple and easy preprocessing method based on isothermal rolling circle mtDNA amplification using commercially available reagents. Isothermal amplification of mtDNA was successfully performed using both nanoliter quantities of plasma directly and 25 ng of total DNA extracted from blood or tissue samples. Prior to mtDNA amplification, it was necessary to treat the extracted total DNA with Exonuclease V, but it was not required to treat plasma. The NGS libraries generated from the amplified mtDNA provided sequencing coverage of the entire human mitochondrial genome. Furthermore, the sequencing results successfully detected heteroplasmy in patient samples, with called mutations and variants matching those from previous, independent, Sanger sequencing analysis. Additionally, a novel single nucleotide variant was detected in a healthy volunteer. The successful analysis of mtDNA using very small samples from patients is likely to be valuable in clinical medicine, as it could reduce patient discomfort by reducing sampling-associated damage to tissues. Overall, the simple and convenient preprocessing method described herein may facilitate the future development of NGS-based clinical and forensic mtDNA tests.

Alignment-free approaches for predicting novel Nuclear Mitochondrial Segments (NUMTs) in the human genome

The nuclear human genome harbors sequences of mitochondrial origin, indicating an ancestral transfer of DNA from the mitogenome. Several Nuclear Mitochondrial Segments (NUMTs) have been detected by alignment-based sequence similarity search, as implemented in the Basic Local Alignment Search Tool (BLAST). Identifying NUMTs is important for the comprehensive annotation and understanding of the human genome. Here we explore the possibility of detecting NUMTs in the human genome by alignment-free sequence similarity search, such as k-mers (k-tuples, k-grams, oligos of length k) distributions. We find that when k=6 or larger, the k-mer approach and BLAST search produce almost identical results, e.g., detect the same set of NUMTs longer than 3 kb. However, when k=5 or k=4, certain signals are only detected by the alignment-free approach, and these may indicate yet unrecognized, and potentially more ancestral NUMTs. We introduce a "Manhattan plot" style representation of NUMT predictions across the genome, which are calculated based on the reciprocal of the Jensen-Shannon divergence between the nuclear and mitochondrial k-mer frequencies. The further inspection of the k-mer-based NUMT predictions however shows that most of them contain long-terminal-repeat (LTR) annotations, whereas BLAST-based NUMT predictions do not. Thus, similarity of the mitogenome to LTR sequences is recognized, which we validate by finding the mitochondrial k-mer distribution closer to those for transposable sequences and specifically, close to some types of LTR.

Pseudogenes as Functionally Significant Elements of the Genome

Pseudogene is a gene copy that has lost its original function. For a long time, pseudogenes have been considered as "junk DNA" that inevitably arises as a result of ongoing evolutionary process. However, experimental data obtained during recent years indicate this understanding of the nature of pseudogenes is not entirely correct, and many pseudogenes perform important genetic functions. In the review, we have addressed classification of pseudogenes, methods of their detection in the genome, and the problem of their evolutionary conservatism and prevalence among species belonging to different taxonomic groups in the light of modern data. The mechanisms of gene expression regulation by pseudogenes and the role of pseudogenes in pathogenesis of various human diseases are discussed.

Development of the MitoQ assay as a real-time quantification of mitochondrial DNA in degraded samples

Mitochondrial DNA is a reliable genetic material for estimating maternally related haplogroups and ancestries. Exploring maternal DNA inheritance is particularly useful when nuclear DNA is degraded or limited, as the copy number of mitochondrial DNA is far greater than the copy number of nuclear DNA. Normal mitochondrial DNA copy number has been estimated to 100 copies per buccal epithelial cell, 4000 copies in skeletal cells and 7000 copies in myocardial cells. This estimation is usually performed via extrapolation from the nuclear DNA quantitation. It is essential to reduce this variability and accurately quantify the exact number of copies of mitochondrial DNA, especially in compromised samples of a forensic or ancient nature. While useful, the testing of mitochondrial DNA is often long and costly and comes with limited success. The accurate quantification of mitochondrial DNA using specific quantitative PCR assays can be used to make better decisions on the downstream testing and success of amplification. As a result, this study develops a real-time assay for the quantification of mitochondrial DNA copy number and assesses its performance on a set of degraded DNA samples. The developed MitoQ assay has been shown to be highly specific to the human mitochondrial genome with no amplification of nuclear pseudogenes being observed and outperformed a previously published concordant assay. Additionally, a high sensitivity was measured to 280 copies of mitochondrial DNA. Minimal variation was observed between each replication cycle, indicating the assay to be robust and repeatable. Overall, this study presents a real-time assay that is sensitive and robust to quantifying mitochondrial DNA copy number in degraded samples. Furthermore, there is potential to incorporate the assay as an additional target in current qPCR assays which use a six-dye chemistry and provide a complete overview of a sample's quality and quantity.

Development of a triplex mtDNA qPCR assay to assess quantification, degradation, inhibition, and amplification target copy numbers

A hybrid absolute/relative qPCR assay which provides information regarding the condition of mitochondrial DNA (mtDNA) in a DNA sample is described. MtDNA concentration (copy number/μL) is determined via absolute quantification using a standard curve of a synthetic duplex DNA previously described (Kavlick et al., 2011). The state of mtDNA degradation is determined via the relative quantification of a mtDNA target found within the 16 s rRNA gene which is 3× longer than that of the short target in the former duplex assay, using the delta, delta Ct (ΔΔCt) method. The presence or absence of PCR inhibitors in the sample is qualitatively determined using a custom internal positive control (IPC) system which targets a unique and non-naturally occurring duplex DNA sequence. This IPC effectively detected inhibition by humic acid, tannic acid, melanin, and EDTA. All three assay components utilize sensitive and specific hydrolysis probes. The utility of ΔΔCt method was demonstrated in a series of experiments involving laboratory-fragmented DNA. Also described is a method for estimating copy number of any mtDNA target longer than the two targets amplified. The described triplex assay works well for intact and for fragmented or degraded mtDNA and therefore may be useful in forensic and ancient DNA disciplines as well as in biomedical research or practice.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

Numtogenesis as a mechanism for development of cancer

Transfer of genetic material from cytoplasmic organelles to the nucleus, an ongoing process, has implications in evolution, aging, and human pathologies such as cancer. The transferred mitochondrial DNA (mtDNA) fragments in the nuclear genome are called nuclear mtDNA or NUMTs. We have named the process numtogenesis, defining the term as the transfer of mtDNA into the nuclear genome, or, less specifically, the transfer of mitochondria or mitochondrial components into the nucleus. There is increasing evidence of the involvement of NUMTs in human biology and pathology. Although information pertaining to NUMTs and numtogenesis is sparse, the role of this aspect of mitochondrial biology to human cancers is apparent. In this review, we present available knowledge about the origin and mechanisms of numtogenesis, with special emphasis on the role of NUMTs in human malignancies. We describe studies undertaken in our laboratory and in others and discuss the influence of NUMTs in tumor initiation and progression and in survival of cancer patients. We describe suppressors of numtogenesis and evolutionary conserved mechanisms underlying numtogenesis in cancer. An understanding the emerging field of numtogenesis should allow comprehension of this process in various malignancies and other diseases and, more generally, in human health.

Unbiased Mitoproteome Analyses Confirm Non-canonical RNA, Expanded Codon Translations

Proteomic MS/MS mass spectrometry detections are usually biased towards peptides cleaved by experimentally added digestion enzyme(s). Hence peptides resulting from spontaneous degradation and natural proteolysis usually remain undetected. Previous analyses of tryptic human proteome data (cleavage after K, R) detected non-canonical tryptic peptides translated according to tetra- and pentacodons (codons expanded by silent mono- and dinucleotides), and from transcripts systematically (a) deleting mono-, dinucleotides after trinucleotides (delRNAs), (b) exchanging nucleotides according to 23 bijective transformations. Nine symmetric and fourteen asymmetric nucleotide exchanges (X ↔ Y, e.g. A ↔ C; and X → Y → Z → X, e.g. A → C → G → A) produce swinger RNAs. Here unbiased reanalyses of these proteomic data detect preferentially non-canonical tryptic peptides despite assuming random cleavage. Unbiased analyses couldn't reconstruct experimental tryptic digestion if most detected non-canonical peptides were false positives. Detected non-tryptic non-canonical peptides map preferentially on corresponding, previously described non-canonical transcripts, as for tryptic non-canonical peptides. Hence unbiased analyses independently confirm previous trypsin-biased analyses that showed translations of del- and swinger RNA and expanded codons. Accounting for natural proteolysis completes trypsin-biased mitopeptidome analyses, independently confirms non-canonical transcriptions and translations.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

In silico discovery of a nearly complete mitochondrial genome Numt in the dog (Canis lupus familiaris) nuclear genome

Through bacterial cloning, a non-specific product co-amplified in a previous whole mitochondrial genome study of Canis lupus familiaris was identified as part of a Numt on chromosome 29 of the dog. Even though further analysis confirmed the fidelity of the mitochondrial genome sequencing results, it still highlighted the risk of Numt contamination. A computer-based search of the dog's nuclear genome for segments homologous to the mtDNA sequence revealed the extent of this risk. Over 150 Numts of various sizes were observed throughout all but two chromosomes, covering all positions of the mtDNA. One of the Numts on chromosome 11 even covered over 95 % of the entire dog mtDNA sequence. This comprehensive list of Numts was provided to assist researchers with the evaluation of dog mtDNA sequencing protocols for Numt co-amplification.

Landscape of insertion polymorphisms in the human genome

Nucleotide substitutions, small (<50 bp) insertions or deletions (indels), and large (>50 bp) deletions are well-known causes of genetic variation within the human genome. We recently reported a previously unrecognized form of polymorphic insertions, termed templated sequence insertion polymorphism (TSIP), in which the inserted sequence was templated from a distant genomic region, and was inserted in the genome through reverse transcription of an RNA intermediate. TSIPs can be grouped into two classes based on nucleotide sequence features at the insertion junctions; class 1 TSIPs show target site duplication, polyadenylation, and preference for insertion at a 5′-TTTT/A-3′ sequence, suggesting a LINE-1 based insertion mechanism, whereas class 2 TSIPs show features consistent with repair of a DNA double strand break by nonhomologous end joining. To gain a more complete picture of TSIPs throughout the human population, we evaluated whole-genome sequence from 52 individuals, and identified 171 TSIPs. Most individuals had 25–30 TSIPs, and common (present in >20% of individuals) TSIPs were found in individuals throughout the world, whereas rare TSIPs tended to cluster in specific geographic regions. The number of rare TSIPs was greater than the number of common TSIPs, suggesting that TSIP generation is an ongoing process. Intriguingly, mitochondrial sequences were a frequent template for class 2 insertions, used more commonly than any nuclear chromosome. Similar to single nucleotide polymorphisms and indels, we suspect that these TSIPs may be important for the generation of human diversity and genetic diseases, and can be useful in tracking historical migration of populations.

Association of single nucleotide polymorphisms in the USF1, GTF2A1L and OR2W3 genes with non-obstructive azoospermia in the Chinese population

PURPOSE

To research the association between the single nucleotide polymorphisms (SNPs) of three spermatogenesis-related genes (USF1, GTF2A1L and OR2W3) and non-obstruction azoospermia (NOA).

METHODS

We investigated 361 NOA cases and 368 controls from the Chinese Han population, and we used Sequenom iplex technology to analyze the candidate 9 SNPs from the USF1, GTF2A1L and OR2W3 genes.

RESULTS

In this study, we found that the variant rs2516838 of USF1 was associated with NOA susceptibility (P = 0.020, OR = 1.436), and the haplotype TCG of the variants rs1556259, rs2516838, and rs2774276 of USF1 conferred an increased risk of NOA (P = 0.019, OR = 1.436). Furthermore, we found that the rs11204546 genotype of OR2W3 and the rs11677854 genotype of GTF2A1L were correlated with the FSH level in the patients (P = 0.004 and P = 0.018, respectively).

CONCLUSIONS

Our results provided a new insight into susceptibility of USF1 variant with male infertility. Clinically, the SNPs (rs11204546 of OR2W3 and rs11677854 of GTF2A1L ) might be additional valuable molecular predictive markers for assessing the treatment of NOA patients.

Transcription factors bind negatively selected sites within human mtDNA genes

Transcription of mitochondrial DNA (mtDNA)-encoded genes is thought to be regulated by a handful of dedicated transcription factors (TFs), suggesting that mtDNA genes are separately regulated from the nucleus. However, several TFs, with known nuclear activities, were found to bind mtDNA and regulate mitochondrial transcription. Additionally, mtDNA transcriptional regulatory elements, which were proved important in vitro, were harbored by a deletion that normally segregated among healthy individuals. Hence, mtDNA transcriptional regulation is more complex than once thought. Here, by analyzing ENCODE chromatin immunoprecipitation sequencing (ChIP-seq) data, we identified strong binding sites of three bona fide nuclear TFs (c-Jun, Jun-D, and CEBPb) within human mtDNA protein-coding genes. We validated the binding of two TFs by ChIP-quantitative polymerase chain reaction (c-Jun and Jun-D) and showed their mitochondrial localization by electron microscopy and subcellular fractionation. As a step toward investigating the functionality of these TF-binding sites (TFBS), we assessed signatures of selection. By analyzing 9,868 human mtDNA sequences encompassing all major global populations, we recorded genetic variants in tips and nodes of mtDNA phylogeny within the TFBS. We next calculated the effects of variants on binding motif prediction scores. Finally, the mtDNA variation pattern in predicted TFBS, occurring within ChIP-seq negative-binding sites, was compared with ChIP-seq positive-TFBS (CPR). Motifs within CPRs of c-Jun, Jun-D, and CEBPb harbored either only tip variants or their nodal variants retained high motif prediction scores. This reflects negative selection within mtDNA CPRs, thus supporting their functionality. Hence, human mtDNA-coding sequences may have dual roles, namely coding for genes yet possibly also possessing regulatory potential.

Interference of Co-amplified nuclear mitochondrial DNA sequences on the determination of human mtDNA heteroplasmy by Using the SURVEYOR nuclease and the WAVE HS system

High-sensitivity and high-throughput mutation detection techniques are useful for screening the homoplasmy or heteroplasmy status of mitochondrial DNA (mtDNA), but might be susceptible to interference from nuclear mitochondrial DNA sequences (NUMTs) co-amplified during polymerase chain reaction (PCR). In this study, we first evaluated the platform of SURVEYOR Nuclease digestion of heteroduplexed DNA followed by the detection of cleaved DNA by using the WAVE HS System (SN/WAVE-HS) for detecting human mtDNA variants and found that its performance was slightly better than that of denaturing high-performance liquid chromatography (DHPLC). The potential interference from co-amplified NUMTs on screening mtDNA heteroplasmy when using these 2 highly sensitive techniques was further examined by using 2 published primer sets containing a total of 65 primer pairs, which were originally designed to be used with one of the 2 techniques. We confirmed that 24 primer pairs could amplify NUMTs by conducting bioinformatic analysis and PCR with the DNA from 143B-ρ0 cells. Using mtDNA extracted from the mitochondria of human 143B cells and a cybrid line with the nuclear background of 143B-ρ0 cells, we demonstrated that NUMTs could affect the patterns of chromatograms for cell DNA during SN-WAVE/HS analysis of mtDNA, leading to incorrect judgment of mtDNA homoplasmy or heteroplasmy status. However, we observed such interference only in 2 of 24 primer pairs selected, and did not observe such effects during DHPLC analysis. These results indicate that NUMTs can affect the screening of low-level mtDNA variants, but it might not be predicted by bioinformatic analysis or the amplification of DNA from 143B-ρ0 cells. Therefore, using purified mtDNA from cultured cells with proven purity to evaluate the effects of NUMTs from a primer pair on mtDNA detection by using PCR-based high-sensitivity methods prior to the use of a primer pair in real studies would be a more practical strategy.

Circulating cell-free mitochondrial DNA as a novel cancer biomarker: opportunities and challenges

The unique characteristics of the mitochondrial genome, such as short length, simple molecular structure, and high copy number, have made monitoring aberrant changes of mitochondrial DNA (mtDNA) quantity an interesting molecular tool for early tumor detection with many advantages over the nuclear genome-based methods. Recently, circulating cell-free (ccf) mtDNA in blood has emerged on the platform as a non-invasive diagnostic and prognostic biomarker for many forms of solid tumors. Accumulating evidence demonstrate that plasma or serum ccf mtDNA levels are significantly different between cancer patients and healthy individuals. Furthermore, quantification of ccf mtDNA levels in circulation may assist in identifying patients from cancer-free healthy population. This minireview attempts to summarize our recent findings in this very promising field of cancer research. The potential technical challenges that we have encountered during the quantitative analysis of ccf mtDNA and mtDNA in general are also briefly discussed. Prospective studies with a larger cohort of patients in various cancer entities are beneficial to precisely define the clinical importance of assessing the ccf mtDNA amount for diagnosing and tracking malignant diseases and their progression.

Mammalian NUMT insertion is non-random

It is well known that remnants of partial or whole copies of mitochondrial DNA, known as Nuclear MiTochondrial sequences (NUMTs), are found in nuclear genomes. Since whole genome sequences have become available, many bioinformatics studies have identified putative NUMTs and from those attempted to infer the factors involved in NUMT creation. These studies conclude that NUMTs represent randomly chosen regions of the mitochondrial genome. There is less consensus regarding the nuclear insertion sites of NUMTs - previous studies have discussed the possible role of retrotransposons, but some recent ones have reported no correlation or even anti-correlation between NUMT sites and retrotransposons. These studies have generally defined NUMT sites using BLAST with default parameters. We analyze a redefined set of human NUMTs, computed with a carefully considered protocol. We discover that the inferred insertion points of NUMTs have a strong tendency to have high-predicted DNA curvature, occur in experimentally defined open chromatin regions and often occur immediately adjacent to A + T oligomers. We also show clear evidence that their flanking regions are indeed rich in retrotransposons. Finally we show that parts of the mitochondrial genome D-loop are under-represented as a source of NUMTs in primate evolution.

Mitochondrial DNA as a non-invasive biomarker: accurate quantification using real time quantitative PCR without co-amplification of pseudogenes and dilution bias

Circulating mitochondrial DNA (MtDNA) is a potential non-invasive biomarker of cellular mitochondrial dysfunction, the latter known to be central to a wide range of human diseases. Changes in MtDNA are usually determined by quantification of MtDNA relative to nuclear DNA (Mt/N) using real time quantitative PCR. We propose that the methodology for measuring Mt/N needs to be improved and we have identified that current methods have at least one of the following three problems: (1) As much of the mitochondrial genome is duplicated in the nuclear genome, many commonly used MtDNA primers co-amplify homologous pseudogenes found in the nuclear genome; (2) use of regions from genes such as β-actin and 18S rRNA which are repetitive and/or highly variable for qPCR of the nuclear genome leads to errors; and (3) the size difference of mitochondrial and nuclear genomes cause a "dilution bias" when template DNA is diluted. We describe a PCR-based method using unique regions in the human mitochondrial genome not duplicated in the nuclear genome; unique single copy region in the nuclear genome and template treatment to remove dilution bias, to accurately quantify MtDNA from human samples.

Nuclear copies of mitochondrial genes: another problem for ancient DNA

The application of ancient DNA techniques is subject to many problems caused primarily by low quality and by low quantity of DNA. For these reasons most studies employing ancient DNA rely on the characterization of mitochondrial DNA, which is present in many more copies per cell than nuclear DNA and hence more copies are likely to survive. We used universal and taxon specific mitochondrial primers to amplify DNA from museum specimens, and found many instances where the amplification of nuclear copies of the mitochondrial gene (numts) instead of the targeted mitochondrial fragment had occurred. Furthermore, the likelihood of amplifying numts increased dramatically when universal primers were utilized. Here we suggest that ancient DNA practitioners must consider the possibility that numts can be amplified at higher rates than previously thought. This is another complication for ancient DNA studies, but it also suggests that more extensive inclusion of nuclear markers in ancient DNA studies should be feasible.

The most redundant sequences in human CpG island library are derived from mitochondrial genome

An altered pattern of epigenetic modifications, such as DNA methylation and histone modification, is critical to many common human diseases, including cancer. Recently, mitochondrial DNA (mtDNA) was reported to be associated with tumorigenesis through epigenetic regulation of methylation patterns. One of the promising approaches to study DNA methylation and CpG islands (CGIs) is sequencing and analysis of clones derived from the physical library generated by methyl-CpG-binding domain proteins and restriction enzyme MseI. In this study, we observed that the most redundant sequences of 349 clones in a human CGI library were all generated from the human mitochondrial genome. Further analysis indicated that there was a 5,845-bp DNA transfer from mtDNA to chromosome 1, and all the clones should be the products of a 510-bp MseI fragment, which contained a putative CGI of 270 bp. The 510-bp fragment was annotated as part of cytochrome c oxidase subunit II (COXII), and phylogenetic analysis of homologous sequences containing COXII showed three DNA transfer events from mtDNA to nuclear genome, one of which underwent secondary transfer events between different chromosomes. These results may further our understanding of how the mtDNA regulates DNA methylation in the nucleus.

Genome-wide mapping of nuclear mitochondrial DNA sequences links DNA replication origins to chromosomal double-strand break formation in Schizosaccharomyces pombe

Chromosomal double-strand breaks (DSBs) threaten genome integrity and repair of these lesions is often mutagenic. How and where DSBs are formed is a major question conveniently addressed in simple model organisms like yeast. NUMTs, nuclear DNA sequences of mitochondrial origin, are present in most eukaryotic genomes and probably result from the capture of mitochondrial DNA (mtDNA) fragments into chromosomal breaks. NUMT formation is ongoing and was reported to cause de novo human genetic diseases. Study of NUMTs is likely to contribute to the understanding of naturally occurring chromosomal breaks. We show that Schizosaccharomyces pombe NUMTs are exclusively located in noncoding regions with no preference for gene promoters and, when located into promoters, do not affect gene transcription level. Strikingly, most noncoding regions comprising NUMTs are also associated with a DNA replication origin (ORI). Chromatin immunoprecipitation experiments revealed that chromosomal NUMTs are probably not acting as ORI on their own but that mtDNA insertions occurred directly next to ORIs, suggesting that these loci may be prone to DSB formation. Accordingly, induction of excessive DNA replication origin firing, a phenomenon often associated with human tumor formation, resulted in frequent nucleotide deletion events within ORI3001 subtelomeric chromosomal locus, illustrating a novel aspect of DNA replication-driven genomic instability. How mtDNA is fragmented is another important issue that we addressed by sequencing experimentally induced NUMTs. This highlighted regions of S. pombe mtDNA prone to breaking. Together with an analysis of human NUMTs, we propose that these fragile sites in mtDNA may correspond to replication pause sites.

Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes

The natural transfer of DNA from mitochondria to the nucleus generates nuclear copies of mitochondrial DNA (numts) and is an ongoing evolutionary process, as genome sequences attest. In humans, five different numts cause genetic disease and a dozen human loci are polymorphic for the presence of numts, underscoring the rapid rate at which mitochondrial sequences reach the nucleus over evolutionary time. In the laboratory and in nature, numts enter the nuclear DNA via non-homolgous end joining (NHEJ) at double-strand breaks (DSBs). The frequency of numt insertions among 85 sequenced eukaryotic genomes reveal that numt content is strongly correlated with genome size, suggesting that the numt insertion rate might be limited by DSB frequency. Polymorphic numts in humans link maternally inherited mitochondrial genotypes to nuclear DNA haplotypes during the past, offering new opportunities to associate nuclear markers with mitochondrial markers back in time.

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.