A “Copernican” Reassessment of the Human Mitochondrial DNA Tree from its Root

Improving the reconstructed sapiens reference sequence of mitochondrial DNA

Comparing Phylogeny and the Predicted Pathogenicity of Protein Variations Reveals Equal Purifying Selection across the Global Human mtDNA Diversity

Objective To investigate the mutation of mitochondrial genome in lymphoma. Methods The peripheral blood or borrow fluid 2 ml from 14 lymphoma patients in the First Hospital of Qinhuangdao between May 2016 and July 2017 were collected. Polymerase chain reaction (PCR) was used to amplify and sequence mitochondrial DNA, and the results were compared with the revised Cambridge reference sequence (rCRS) and human mitochondrial genome database (mtDB), and then the mutation was also analyzed. Results There were 118 mutation genes, including 57.63% (68/118) in D-loop region, 18.64% (22/118) in NADH dehydrogenase 5 (ND5) region, 13.56% (16/118) in cytochrome b oxidase (CbO) region, 5.08% (6/118) in ND1 region, 3.39% (4/118) in cytochrome oxidase (COⅡ) region, 1.69% (2/118) in ND4 region. Conclusion Mitochondrial DNA mutation in lymphoma has a high mutation rate. Key words: Lymphoma; Genes, mitochondrial; Mutation

An Oligonucleotide Microarray for High-Throughput Sequencing of the Mitochondrial Genome

Objective To study the mitochondrial DNA mutation in leukemia. Methods Mitochondrial DNA of 16 leukemia patients in First Hospital of Qinhuangdao from February to June 2017 were amplified and sequenced by using polymerase chain reaction (PCR). The result was compared with revised Cambridge reference sequence (rCRS) and human mitochondrial genome database (mtDB), and the mutation was also analyzed. Results There were 106 mutation genes in total, including 47.17 % (50/106) in D-loop region, 2.83 % (3/106) in ND4 region, 17.92 % (19/106) in ND5 region, 22.64 % (24/106) in Cytb region, 7.55 % (8/106) in ND1 region, 1.89 % (2/106) in CoⅡregion. Conclusion There is a high mitochondrial DNA mutation rate in leukemia patients. Key words: Leukemia; Mitochondrial DNA mutation

Since the determination in 1981 of the sequence of the human mitochondrial DNA (mtDNA) genome, the Cambridge Reference Sequence (CRS), has been used as the reference sequence to annotate mtDNA in molecular anthropology, forensic science and medical genetics. The CRS was eventually upgraded to the revised version (rCRS) in 1999. This reference sequence is a convenient device for recording mtDNA variation, although it has often been misunderstood as a wild-type (WT) or consensus sequence by medical geneticists. Recently, there has been a proposal to replace the rCRS with the so-called Reconstructed Sapiens Reference Sequence (RSRS). Even if it had been estimated accurately, the RSRS would be a cumbersome substitute for the rCRS, as the new proposal fuses--and thus confuses--the two distinct concepts of ancestral lineage and reference point for human mtDNA. Instead, we prefer to maintain the rCRS and to report mtDNA profiles by employing the hitherto predominant circumfix style. Tree diagrams could display mutations by using either the profile notation (in conventional short forms where appropriate) or in a root-upwards way with two suffixes indicating ancestral and derived nucleotides. This would guard against misunderstandings about reporting mtDNA variation. It is therefore neither necessary nor sensible to change the present reference sequence, the rCRS, in any way. The proposed switch to RSRS would inevitably lead to notational chaos, mistakes and misinterpretations.

Canis MitoSNP is a tool allowing assignment of each mitochondrial genomic position a corresponding position in the mitochondrial gene and in the structure of tRNA, rRNA, and protein. The main aim of this bioinformatic tool was to use data from other bioinformatic tools (TMHMM, SOPMA, tRNA-SCAN, RNAfold, ConSurf) for dog and human mitochondrial genes in order to shorten the time necessary for the analysis of the whole genome single nucleotide polymorphism (SNP) as well as amino acid and protein analyses. Each position in the canine mitochondrial genome is assigned a position in genes, in codons, an amino acid position in proteins, or a position in tRNA or rRNA molecules. Therefore, a user analysing changes in the canine and human mitochondrial genome does not need to extract the sequences of individual genes from the mitochondrial genome for analysis and there is no need to rewrite them into amino acid sequences to assess whether the change is synonymous or nonsynonymous. Canis mitoSNP allows the comparison between the human and canine mitochondrial genomes as well. The Clustal W alignment of the dog and human mitochondrial DNA reference sequences for each gene obtained from GenBank (NC_002008.4 dog, NC_012920.1 human) was performed in order to determine which position in the canine mitochondrial genome corresponds to the position in the human mitochondrial genome. This function may be useful for the comparative analyses. The tool is available at: https://canismitosnp.pl .

Mutation analysis in the mitochondrial DNA (mtDNA) control region is widely used in population genetic studies as well as in forensic medicine. Among the difficulties linked to the mtDNA analysis, one can find the detection of heteroplasmy, which can be inherited or somatic. Recently, age-related point mutation A189G was described in mtDNA and shown to accumulate with age in muscles. We carried out the detection of this 189 heteroplasmic point mutation using three technologies: automated DNA sequencing, Southern blot hybridization using a digoxigenin-labeled oligonucleotide probe, and peptide nucleic acid (PNA)/real-time PCR combined method on different biological samples. Our results give additional information on the increase in mutation frequency with age in muscle tissue and revealed that the PNA/real-time PCR is a largely more sensitive method than DNA sequencing for heteroplasmy detection. These investigations could be of interest in the detection and interpretation of mtDNA heteroplasmy in anthropological and forensic studies.

HmtDB (http://www.hmtdb.uniba.it:8080/hmdb) is a open resource created to support population genetics and mitochondrial disease studies. The database hosts human mitochondrial genome sequences annotated with population and variability data, the latter being estimated through the application of the SiteVar software based on site-specific nucleotide and amino acid variability calculations. The annotations are manually curated thus adding value to the quality of the information provided to the end-user. Classifier tools implemented in HmtDB allow the prediction of the haplogroup for any human mitochondrial genome currently stored in HmtDB or externally submitted de novo by an end-user. Haplogroup definition is based on the Phylotree system. End-users accessing HmtDB are hence allowed to (i) browse the database through the use of a multi-criterion ‘query’ system; (ii) analyze their own human mitochondrial sequences via the ‘classify’ tool (for complete genomes) or by downloading the ‘fragment-classifier’ tool (for partial sequences); (iii) download multi-alignments with reference genomes as well as variability data.

The Diversity Present in 5140 Human Mitochondrial Genomes

Whole genome sequencing for generating SNP data is increasingly used in population genetic studies. However, obtaining genomes for massive numbers of samples is still not within the budgets of many researchers. It is thus imperative to select an appropriate reference genome and sequencing depth to ensure the accuracy of the results for a specific research question, while balancing cost and feasibility. To evaluate the effect of the choice of the reference genome and sequencing depth on downstream analyses, we used five confamilial reference genomes of variable relatedness and three levels of sequencing depth (3.5×, 7.5× and 12×) in a population genomic study on two caddisfly species: Himalopsyche digitata and H. tibetana. Using these 30 datasets (five reference genomes × three depths × two target species), we estimated population genetic indices (inbreeding coefficient, nucleotide diversity, pairwise FST, and genome‐wide distribution of FST) based on variants and population structure (PCA and admixture) based on genotype likelihood estimates. The results showed that both distantly related reference genomes and lower sequencing depth lead to degradation of resolution. In addition, choosing a more closely related reference genome may significantly remedy the defects caused by low depth. Therefore, we conclude that population genetic studies would benefit from closely related reference genomes, especially as the costs of obtaining a high‐quality reference genome continue to decrease. However, to determine a cost‐efficient strategy for a specific population genomic study, a trade‐off between reference genome relatedness and sequencing depth can be considered.

BackgroundThe complete and accurate human reference genome is important for functional genomics researches. Therefore, the incomplete reference genome and individual specific sequences have significant effects on various studies.Resultswe used two RNA-Seq datasets from human brain tissues and 10 mixed cell lines to investigate the completeness of human reference genome. First, we demonstrated that in previously identified ~5 Mb Asian and ~5 Mb African novel sequences that are absent from the human reference genome of NCBI build 36, ~211 kb and ~201 kb of them could be transcribed, respectively. Our results suggest that many of those transcribed regions are not specific to Asian and African, but also present in Caucasian. Then, we found that the expressions of 104 RefSeq genes that are unalignable to NCBI build 37 in brain and cell lines are higher than 0.1 RPKM. 55 of them are conserved across human, chimpanzee and macaque, suggesting that there are still a significant number of functional human genes absent from the human reference genome. Moreover, we identified hundreds of novel transcript contigs that cannot be aligned to NCBI build 37, RefSeq genes and EST sequences. Some of those novel transcript contigs are also conserved among human, chimpanzee and macaque. By positioning those contigs onto the human genome, we identified several large deletions in the reference genome. Several conserved novel transcript contigs were further validated by RT-PCR.ConclusionOur findings demonstrate that a significant number of genes are still absent from the incomplete human reference genome, highlighting the importance of further refining the human reference genome and curating those missing genes. Our study also shows the importance of de novo transcriptome assembly. The comparative approach between reference genome and other related human genomes based on the transcriptome provides an alternative way to refine the human reference genome.

Capturing the Perfect Reference Genome

The persistence of human DNA in soil following surface decomposition

BackgroundPopulation genetics studies based on the analysis of mtDNA and mitochondrial disease studies have produced a huge quantity of sequence data and related information. These data are at present worldwide distributed in differently organised databases and web sites not well integrated among them. Moreover it is not generally possible for the user to submit and contemporarily analyse its own data comparing them with the content of a given database, both for population genetics and mitochondrial disease data.ResultsHmtDB is a well-integrated web-based human mitochondrial bioinformatic resource aimed at supporting population genetics and mitochondrial disease studies, thanks to a new approach based on site-specific nucleotide and aminoacid variability estimation. HmtDB consists of a database of Human Mitochondrial Genomes, annotated with population data, and a set of bioinformatic tools, able to produce site-specific variability data and to automatically characterize newly sequenced human mitochondrial genomes. A query system for the retrieval of genomes and a web submission tool for the annotation of new genomes have been designed and will soon be implemented. The first release contains 1255 fully annotated human mitochondrial genomes. Nucleotide site-specific variability data and multialigned genomes can be downloaded. Intra-human and inter-species aminoacid variability data estimated on the 13 coding for proteins genes of the 1255 human genomes and 60 mammalian species are also available. HmtDB is freely available, upon registration, at .ConclusionThe HmtDB project will contribute towards completing and/or refining haplogroup classification and revealing the real pathogenic potential of mitochondrial mutations, on the basis of variability estimation.