Mitochondrial DNA enrichment reduced NUMT contamination in porcine NGS analyses

Biological and translational attributes of mitochondrial DNA copy number: Laboratory perspective to clinical relevance

The mitochondrial DNA copy number (mtDNAcn) plays a vital role in cellular energy metabolism and mitochondrial health. As mitochondria are responsible for adenosine triphosphate production through oxidative phosphorylation, maintaining an appropriate mtDNAcn level is vital for the overall cellular function. Alterations in mtDNAcn have been linked to various diseases, including neurodegenerative disorders, metabolic conditions, and cancers, making it an important biomarker for understanding the disease pathogenesis. The accurate estimation of mtDNAcn is essential for clinical applications. Quantitative polymerase chain reaction and next-generation sequencing are commonly employed techniques with distinct advantages and limitations. Clinically, mtDNAcn serves as a valuable indicator for early diagnosis, disease progression, and treatment response. For instance, in oncology, elevated mtDNAcn levels in blood samples are associated with tumor aggressiveness and can aid in monitoring treatment efficacy. In neurodegenerative diseases such as Alzheimer’s and Parkinson’s, altered mtDNAcn patterns provide insights into disease mechanisms and progression. Understanding and estimating mtDNAcn are critical for advancing diagnostic and therapeutic strategies in various medical fields. As research continues to uncover the implications of mtDNAcn alterations, its potential as a clinical biomarker is likely to expand, thereby enhancing our ability to diagnose and manage complex diseases.

Non-Coding RNAs of Mitochondrial Origin: Roles in Cell Division and Implications in Cancer

Non-coding RNAs (ncRNAs) are a heterogeneous group, in terms of structure and sequence length, consisting of RNA molecules that do not code for proteins. These ncRNAs have a central role in the regulation of gene expression and are virtually involved in every process analyzed, ensuring cellular homeostasis. Although, over the years, much research has focused on the characterization of non-coding transcripts of nuclear origin, improved bioinformatic tools and next-generation sequencing (NGS) platforms have allowed the identification of hundreds of ncRNAs transcribed from the mitochondrial genome (mt-ncRNA), including long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miR). Mt-ncRNAs have been described in diverse cellular processes such as mitochondrial proteome homeostasis and retrograde signaling; however, the function of the majority of mt-ncRNAs remains unknown. This review focuses on a subgroup of human mt-ncRNAs whose dysfunction is associated with both failures in cell cycle regulation, leading to defects in cell growth, cell proliferation, and apoptosis, and the development of tumor hallmarks, such as cell migration and metastasis formation, thus contributing to carcinogenesis and tumor development. Here we provide an overview of the mt-ncRNAs/cancer relationship that could help the future development of new biomedical applications in the field of oncology.

Experimental evaluation of genetic variability based on DNA metabarcoding from the aquatic environment: Insights from the Leray COI fragment

Intraspecific genetic variation is important for the assessment of organisms' resistance to changing environments and anthropogenic pressures. Aquatic DNA metabarcoding provides a non-invasive method in biodiversity research, including investigations at the within-species level. Through the analysis of eDNA samples collected from the Peter the Great Gulf of the Japan Sea, in this study, we aimed to evaluate the identification of Amplicon Sequence Variants (ASVs) in marine eDNA among abundant species of the Zostera sp. community: Hexagrammos octogrammus, Pholidapus dybowskii (Teleostei: Perciformes), and Pandalus latirostris (Arthropoda: Decapoda). These species were collected from two distant locations to produce mock communities and gather aquatic eDNA both on the community and individual level. Our approach highlights the efficacy of eDNA metabarcoding in capturing haplotypic diversity and the potential for this methodology to track genetic diversity accurately, contributing to conservation efforts and ecosystem management. Additionally, our results elucidate the impact of nuclear mitochondrial DNA segments (NUMTs) on the reliability of metabarcoding data, indicating the necessity for cautious interpretation of such data in ecological studies. Moreover, we analyzed 83 publicly available COI sequence datasets from common groups of multicellular organisms (Mollusca, Echinodermata, Crustacea, Polychaeta, and Actinopterygii). The results reflect the decrease in population diversity that arises from using the metabarcode compared to the COI barcode.

A PCR-independent approach for mtDNA enrichment and next-generation sequencing: comprehensive evaluation and clinical application

BACKGROUND

Sequencing the mitochondrial genome has been increasingly important for the investigation of primary mitochondrial diseases (PMD) and mitochondrial genetics. To overcome the limitations originating from PCR-based mtDNA enrichment, we set out to develop and evaluate a PCR-independent approach in this study, named Pime-Seq (PCR-independent mtDNA enrichment and next generation Sequencing).

RESULTS

By using the optimized mtDNA enrichment procedure, the mtDNA reads ratio reached 88.0 ± 7.9% in the sequencing library when applied on human PBMC samples. We found the variants called by Pime-Seq were highly consistent among technical repeats. To evaluate the accuracy and reliability of this method, we compared Pime-Seq with lrPCR based NGS by performing both methods simultaneously on 45 samples, yielding 1677 concordant variants, as well as 146 discordant variants with low-level heteroplasmic fraction, in which Pime-Seq showed higher reliability. Furthermore, we applied Pime-Seq on 4 samples of PMD patients retrospectively, and successfully detected all the pathogenic mtDNA variants. In addition, we performed a prospective study on 192 apparently healthy pregnant women during prenatal screening, in which Pime-Seq identified pathogenic mtDNA variants in 4 samples, providing extra information for better health monitoring in these cases.

CONCLUSIONS

Pime-Seq can obtain highly enriched mtDNA in a PCR-independent manner for high quality and reliable mtDNA deep-sequencing, which provides us an effective and promising tool for detecting mtDNA variants for both clinical and research purposes.

ATF4 renders human T-cell acute lymphoblastic leukemia cell resistance to FGFR1 inhibitors through amino acid metabolic reprogramming

Abnormalities of FGFR1 have been reported in multiple malignancies, suggesting FGFR1 as a potential target for precision treatment, but drug resistance remains a formidable obstacle. In this study, we explored whether FGFR1 acted a therapeutic target in human T-cell acute lymphoblastic leukemia (T-ALL) and the molecular mechanisms underlying T-ALL cell resistance to FGFR1 inhibitors. We showed that FGFR1 was significantly upregulated in human T-ALL and inversely correlated with the prognosis of patients. Knockdown of FGFR1 suppressed T-ALL growth and progression both in vitro and in vivo. However, the T-ALL cells were resistant to FGFR1 inhibitors AZD4547 and PD-166866 even though FGFR1 signaling was specifically inhibited in the early stage. Mechanistically, we found that FGFR1 inhibitors markedly increased the expression of ATF4, which was a major initiator for T-ALL resistance to FGFR1 inhibitors. We further revealed that FGFR1 inhibitors induced expression of ATF4 through enhancing chromatin accessibility combined with translational activation via the GCN2-eIF2α pathway. Subsequently, ATF4 remodeled the amino acid metabolism by stimulating the expression of multiple metabolic genes ASNS, ASS1, PHGDH and SLC1A5, maintaining the activation of mTORC1, which contributed to the drug resistance in T-ALL cells. Targeting FGFR1 and mTOR exhibited synergistically anti-leukemic efficacy. These results reveal that FGFR1 is a potential therapeutic target in human T-ALL, and ATF4-mediated amino acid metabolic reprogramming contributes to the FGFR1 inhibitor resistance. Synergistically inhibiting FGFR1 and mTOR can overcome this obstacle in T-ALL therapy.

Mito-SiPE is a sequence-independent and PCR-free mtDNA enrichment method for accurate ultra-deep mitochondrial sequencing

The analysis of somatic variation in the mitochondrial genome requires deep sequencing of mitochondrial DNA. This is ordinarily achieved by selective enrichment methods, such as PCR amplification or probe hybridization. These methods can introduce bias and are prone to contamination by nuclear-mitochondrial sequences (NUMTs), elements that can introduce artefacts into heteroplasmy analysis. We isolated intact mitochondria using differential centrifugation and alkaline lysis and subjected purified mitochondrial DNA to a sequence-independent and PCR-free method to obtain ultra-deep (>80,000X) sequencing coverage of the mitochondrial genome. This methodology avoids false-heteroplasmy calls that occur when long-range PCR amplification is used for mitochondrial DNA enrichment. Previously published methods employing mitochondrial DNA purification did not measure mitochondrial DNA enrichment or utilise high coverage short-read sequencing. Here, we describe a protocol that yields mitochondrial DNA and have quantified the increased level of mitochondrial DNA post-enrichment in 7 different mouse tissues. This method will enable researchers to identify changes in low frequency heteroplasmy without introducing PCR biases or NUMT contamination that are incorrectly identified as heteroplasmy when long-range PCR is used.

Do pseudogenes pose a problem for metabarcoding marine animal communities?

Because DNA metabarcoding typically employs sequence diversity among mitochondrial amplicons to estimate species composition, nuclear mitochondrial pseudogenes (NUMTs) can inflate diversity. This study quantifies the incidence and attributes of NUMTs derived from the 658 bp barcode region of cytochrome c oxidase I (COI) in 156 marine animal genomes. NUMTs were examined to ascertain if they could be recognized by their possession of indels or stop codons. In total, 309 NUMTs ≥ 150 bp were detected, with an average of 1.98 per species (range = 0-33) and a mean length of 391 bp ± 200 bp. Among this total, 75 (24.3 %) lacked indels or stop codons. NUMTs appear to pose the greatest interpretational risk when short (< 313 bp) amplicons are used, such as in eDNA studies, dietary analyses, or processed fish identification. Employing the standard amplicon length (313 bp) for marine metabarcoding, NUMTs could potentially inflate the OTU count by 21% above the true species count while also raising intraspecific variation at COI by 15%. However, when both amplicon length and position are considered, inflation in OTU counts and in barcode variation were just 9% and 10%, respectively, suggesting NUMTs will not seriously distort biodiversity assessments. There was a weak positive correlation between genome size and NUMT count but no variation among phyla or trophic groups. Until bioinformatic advances improve NUMT detection, the best defense involves targeting long amplicons and developing reference databases that include both mitochondrial sequences and their NUMT derivatives.

Breed-specific reference sequence optimized mapping accuracy of NGS analyses for pigs

Background

Reference sequences play a vital role in next-generation sequencing (NGS), impacting mapping quality during genome analyses. However, reference genomes usually do not represent the full range of genetic diversity of a species as a result of geographical divergence and independent demographic events of different populations. For the mitochondrial genome (mitogenome), which occurs in high copy numbers in cells and is strictly maternally inherited, an optimal reference sequence has the potential to make mitogenome alignment both more accurate and more efficient. In this study, we used three different types of reference sequences for mitogenome mapping, i.e., the commonly used reference sequence (CU-ref), the breed-specific reference sequence (BS-ref) and the sample-specific reference sequence (SS-ref), respectively, and compared the accuracy of mitogenome alignment and SNP calling among them, for the purpose of proposing the optimal reference sequence for mitochondrial DNA (mtDNA) analyses of specific populations

Results

Four pigs, representing three different breeds, were high-throughput sequenced, subsequently mapping reads to the reference sequences mentioned above, resulting in a largest mapping ratio and a deepest coverage without increased running time when aligning reads to a BS-ref. Next, single nucleotide polymorphism (SNP) calling was carried out by 18 detection strategies with the three tools SAMtools, VarScan and GATK with different parameters, using the bam results mapping to BS-ref. The results showed that all eighteen strategies achieved the same high specificity and sensitivity, which suggested a high accuracy of mitogenome alignment by the BS-ref because of a low requirement for SNP calling tools and parameter choices.

Conclusions

This study showed that different reference sequences representing different genetic relationships to sample reads influenced mitogenome alignment, with the breed-specific reference sequences being optimal for mitogenome analyses, which provides a refined processing perspective for NGS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08030-1.