The effect of library preparation protocol on the efficiency of heteroplasmy detection in mitochondrial DNA using two massively parallel sequencing Illumina systems.

Massively parallel sequencing (MPS) technology has become the gold standard in mitochondrial DNA research due to its high sensitivity in detecting mtDNA heteroplasmy, a prognostic marker in various medical applications. Various MPS technologies and platforms used for mtDNA analysis exist. Obtaining reliable and sensitive results requires deep and uniform coverage of the entire mtDNA sequence, which is heavily influenced by the choice of library preparation method and sequencing platform. Here, we present a comparison of the sequencing coverage and the ability to heteroplasmy detection using two library preparation protocols (Nextera XT DNA Library Preparation Kit and Nextera DNA Flex Library Preparation Kit) and two different (MiSeq FGx and ISeq 100) Illumina MPS platforms. Our study indicates that the Nextera DNA Flex Library protocol provides a more balanced coverage along the mitogenome and a reliable heteroplasmy detection with both MiSeq and iSeq Illumina MPS systems.

Capability of the iSeq 100 sequencing system from Illumina to detect low-level substitutions in the human mitochondrial genome.

The significance of mtDNA heteroplasmy in forensic and medical genetics has increased recently because massively parallel sequencing (MPS) technologies enable more accurate and precise detection of minority nucleotide variants. Recent reports have shown that detection of low-level substitutions may depend on library preparation or sequencing protocol, and can vary for different MPS platforms. The MiSeq (Illumina) and Ion S5 (Thermo Fisher Scientific) are mainly used for heteroplasmy detection, but no data are available regarding the iSeq 100, an Illumina platform of the smallest throughput. Notably, unlike the other systems, the machine utilizes sequencing by synthesis one-channel chemistry to determine DNA sequences. Thus, it is important to verify the capability of the iSeq 100 system to determine mitochondrial haplotypes and detect heteroplasmic substitutions. In this study, previously determined entire mitochondrial genomes were sequenced with the iSeq 100 system. Each mitogenome was sequenced twice, giving approximately 2000x and 10,000x coverage. All homoplasmic mutations and minority variants above the 19 % level detected with the iSeq 100 system were also observed after dideoxy sequencing. Moreover, all heteroplasmic substitutions above the 2 % level were consistently detected with SBS one-channel chemistry. However, detection of low-level mtDNA variants may require additional, confirmatory experiments. In summary, the iSeq 100 system enables reproducible and accurate sequencing of human mitochondrial genomes. Detection of mtDNA minority variants depends on the laboratory protocol and sequencing platform used, but homoplasmic mutations and heteroplasmy above the 2 % level can be correctly detected with the iSeq 100 system.

Low-level point heteroplasmy detection in human mitogenomes amplified with different polymerases and sequenced on MiSeq FGx platform.

Introduction: Massively parallel sequencing of mitogenomes usually requires prior amplification. The PCR step may influence the quality of the data obtained, especially when low-level heteroplasmy detection is applied. Aim: The aim of this study was to compare the reliability of two different DNA polymerases in detecting homoplasmic and heteroplasmic substitutions in human mitogenomes. Material and methods: Mitogenomes of five samples were amplified with Long PCR Enzyme Mix from Fermentas or TaKaRa LA Taq DNA Polymerase from TaKaRa. Then, NexteraTM XT DNA libraries were sequenced on MiSeq FGx platform (Illumina). mtDNA substitutions were called for alternative variants above the 1% level. Results: All homoplasmic substitutions detected in amplicons generated with polymerases studied here and sequenced on MiSeq FGx system were consistently identified as homoplasmies with alternative sequencing methods. TaKaRa LA Taq DNA Polymerase was found to be less accurate in low-level heteroplasmy detection than Long PCR Enzyme Mix enzyme as more false negative and false positive results were observed for minority variants called above the 1% level. Nevertheless, both PCR systems studied can be successfully used to detect authentic mtDNA substitutions, for which minority variants exceed the 3.61% level assuming at least 10,000x coverage and sequencing Nextera XT DNA libraries on MiSeq FGx machine. Conclusions: The accuracy and sensitivity of point heteroplasmy detection with the MiSeq FGx instrument varies on polymerase used for mtDNA amplification. Therefore, it is recommended to validate the laboratory protocols used for mtDNA substitution detection prior to their implementation for the forensic or medical genetics purposes.

Variability of the rs333 in Polish patients with lupus erythematosus.

INTRODUCTION: Lupus erythematosus (LE) is an autoimmune disease with a strong influence of genetic and environmental factors. C-C motif chemokine receptor 5 (CCR5) gene expression may affect the development and intensity of LE. AIM: To evaluate the possible association between the 32bp deletion in rs333 locus located within the CCR5 gene and the development of LE or the occurrence of various clinical symptoms in the course of the disease. MATERIAL AND METHODS: One hundred and twenty patients with LE (77 with systemic lupus erythematosus (SLE) and 43 with discoid lupus erythematosus (DLE)) and 100 healthy controls from the Polish population were genotyped for deletion in rs333. RESULTS: 32 bp deletion in the rs333 was significantly more frequent among healthy individuals than DLE patients. Moreover, heterozygotes and homozygotes with deletion in rs333 were significantly more frequent within the control group than the group of patients with discoid lupus erythematosus. In contrast, any statistically significant differences in allele or genotype frequencies between healthy persons and SLE patients were observed. Furthermore, nucleotide sequence variability of rs333 was not associated with certain clinical symptoms of LE patients. CONCLUSIONS: Deletion in the rs333 might be a protective factor for DLE, but not SLE in the Polish population. Nevertheless further studies performed on larger populations are needed to confirm these observations.

Mitogenome germline mutations and colorectal cancer risk in Polish population.

INTRODUCTION: To date, several nuclear DNA variants have been shown to be associated with increased risk of developing colorectal cancer. Despite the fact that mitochondria play an important role in carcinogenesis, little is known about inherited mitochondrial DNA mutations that could be involved in this disease. Thus, potential associations between inherited mutations in the entire mitochondrial genomes and colorectal cancer were analysed in this study. MATERIAL AND METHODS: Two hundred mitogenome sequences determined for colorectal cancer patients and healthy individuals from Poland were used to investigate the association between mtDNA alleles or haplogroups and colorectal cancer. Additional mtDNA control region haplotypes determined for 1353 individuals from the general Polish population were used for comparison of haplogroup and certain allele frequencies between case and control groups. RESULTS: The non-R clades together with their diagnostic T alleles at positions 12705 and 16223 were observed with higher frequencies in healthy individuals than in colorectal cancer patients. Nevertheless, the differences of the R macrohaplogroup (as well as 12705 or 16223 alleles) frequencies between cases and controls were statistically insignificant after Bonferroni correction. Most of the non-R clades were of Asian and African origin, but none of them were prevalent in the control group. Moreover, neither mtDNA alleles nor haplogroups were associated with clinicopathological parameters of colorectal cancer patients. CONCLUSIONS: Contrary to some previous reports, the findings of this study do not support the hypothesis that mitochondrial DNA variants contribute to inherited predisposition to colorectal cancer.

TP53 somatic mutations are associated with somatic mitogenome substitutions but not indels in colorectal cancer cells.

BACKGROUND: p53 is a tumour suppressor protein that is involved in many cancer-related processes. Growing evidence suggests that p53 also plays an important role in mitochondrial (mtDNA) maintenance. Somatic mitogenome mutations are frequently observed in colorectal cancer (CC) cells. Thus, it was important to determine whether somatic mtDNA changes are associated with TP53 mutational status. METHODS: In the present study, we analysed the TP53 gene in 67 CC patients, for whom mitogenome haplotypes were previously described. In total, 134 TP53 sequences (of cancer and matched normal specimens) were determined using the dideoxy method. RESULTS: Nine hereditary polymorphisms in the TP53 gene were detected in normal colon cells. None of them (neither alleles, nor genotypes) was associated with somatic mitogenome mutations in CC cells. Moreover, 42 somatic TP53 mutations were found in approximately 36% of CC tissues. These somatic changes were significantly more frequent in CC cells with somatic mtDNA mutations (p = 0.0069). Furthermore, we show that only mitochondrial somatic substitutions (p = 0.0017), but not indels (p > 0.05), were associated with somatic TP53 mutations. CONCLUSIONS: The results of the present study suggest that changes in TP53 may modify p53 properties, which may result in the accumulation of somatic substitutions in CC mitogenomes.

Whole mitochondrial genome diversity in two Hungarian populations.

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

X-linked TLR7 gene polymorphisms are associated with diverse immunological conditions but not with discoid lupus erythematosus in Polish patients.

INTRODUCTION: Toll-like receptor 7 (TLR7) is an important molecule involved in the development of autoimmunity and the response to different pathogens. Several polymorphisms within the TLR7 gene were previously found to be associated with systemic lupus erythematosus (SLE). However, none of those studies investigated the TLR7 promoter flanking variants rs1634318 and rs1616583. TLR7 gene diversity has not been analyzed with respect to discoid lupus erythematosus (DLE) development, while its role in the human immunological response to fungal infection is not fully known. AIM: To clarify the potential involvement of two novel single-nucleotide polymorphisms (SNPs) located in the TLR7 gene (rs1634318 and rs1616583) in a variety of immune-related conditions, we studied the variability of these loci in patients from a Polish population with SLE and DLE, as well as in immunocompromised patients who were affected by invasive aspergillosis (IA) and those who were not affected. MATERIAL AND METHODS: Real-time polymerase chain reaction was used to genotype SNPs. Statistically significant differences between case and control groups for both allele and genotype frequencies were assessed using the χ2 test with Yates' correction or two-tailed Fisher's exact test. The results were Bonferroni-corrected for multiple comparisons and odds ratios were calculated. RESULTS: Two polymorphisms located in TLR7 might be associated with the development of SLE but not DLE within the Polish population. Moreover, variation of the two investigated SNPs was found to be associated with IA in immunocompromised Polish patients. CONCLUSIONS: In Polish patients, TLR7 promoter flanking gene polymorphisms might be associated with IA and SLE but not DLE.

Mitogenomic differences between the normal and tumor cells of colorectal cancer patients.

So far, a reliable spectrum of mitochondrial DNA mutations in colorectal cancer cells is still unknown, and neither is their significance in carcinogenesis. Indeed, it remains debatable whether mtDNA mutations are "drivers" or "passengers" of colorectal carcinogenesis. Thus, we analyzed 200 mitogenomes from normal and cancer tissues of 100 colorectal cancer patients. Minority variant mutations were detected at the 1% level. We showed that somatic mutations frequently occur in colorectal cancer cells (75%) and are randomly distributed across the mitochondrial genome. Mutational signatures of somatic mitogenome mutations suggest that they might arise through nucleotide deamination due to oxidative stress. The majority of somatic mutations localized within the coding region (in positions not known from the human phylogeny) and was potentially pathogenic to cell metabolism. Further analysis suggested that the relaxation of negative selection in the mitogenomes of colorectal cancer cells may allow accumulation of somatic mutations. Thus, a shift in glucose metabolism from oxidative phosphorylation to glycolysis may create advantageous conditions for accumulation of mtDNA mutations. Considering the fact that the presence of somatic mtDNA mutations was not associated with any clinicopathological features, we suggested that mtDNA somatic mutations are "passengers" rather than the cause of colorectal carcinogenesis.

Genetic similarities and differences between discoid and systemic lupus erythematosus patients within the Polish population.

INTRODUCTION: Many studies have shown that some SNPs might be a risk factor for systemic lupus erythematosus (SLE), but little is known about potential susceptibility loci of the skin types of the disease. Discoid lupus erythematosus (DLE) is the most common form of the cutaneous lupus erythematosus. Nevertheless, a genetic contribution to DLE is not fully recognized. AIM: We aimed to analyze three SNPs located in the STAT4 (rs7574865), ITGAM (rs1143679) and TNXB (rs1150754) genes in both DLE and SLE patients from Poland. MATERIAL AND METHODS: SNPs were genotyped using real-time polymerase chain reaction (PCR). Statistical significance of the differences between patient and control groups in both allele and genotype frequencies were calculated using two tailed Fisher's exact test. The correction for multiple testing by the Bonferroni adjustment and odds ratio were also calculated. RESULTS: For the first time, we have shown that the polymorphisms located in the STAT4 (rs7574865), but not in the ITGAM (rs1143679) nor the TNXB (rs1150754) genes, might be associated with the development of DLE within the Polish population. The variation of the three investigated SNPs was found to be associated with SLE in our dataset. CONCLUSIONS: The results of our study suggest differences in the molecular background between DLE and SLE within the Polish population.

Mitogenomic diversity in Russians and Poles.

Complete mtDNA genome sequencing improves molecular resolution for distinguishing variation between individuals and populations, but there is still deficiency of mitogenomic population data. To overcome this limitation, we used Sanger-based protocol to generate complete mtDNA sequences of 376 Russian individuals from six populations of European part of Russia and 100 Polish individuals from northern Poland. Nearly complete resolution of mtDNA haplotypes was achieved - about 97% of haplotypes were unique both in Russians and Poles, and no haplotypes overlapped between them when indels were considered. While European populations showed a low, but statistically significant level of between-population differentiation (Fst=0.66%, p=0), Russians demonstrate lack of between-population differences (Fst=0.22%, p=0.15). Results of the Bayesian skyline analysis of Russian mitogenomes demonstrate not only post-Last Glacial Maximum expansion, but also rapid population growth starting from about 4.3kya (95% CI: 2.9-5.8kya), i.e. in the Bronze Age. This expansion strongly correlates with the Kurgan model established by archaeologists and confirmed by paleogeneticists.

Mitochondrial super-haplogroup U diversity in Serbians.

BACKGROUND: Available mitochondrial (mtDNA) data demonstrate genetic differentiation among South Slavs inhabiting the Balkan Peninsula. However, their resolution is insufficient to elucidate the female-specific aspects of the genetic history of South Slavs, including the genetic impact of various migrations which were rather common within the Balkans, a region having a turbulent demographic history. AIM: The aim was to thoroughly study complete mitogenomes of Serbians, a population linking westward and eastward South Slavs. SUBJECTS AND METHODS: Forty-six predominantly Serbian super-haplogroup U complete mitogenomes were analysed phylogenetically against ∼4000 available complete mtDNAs of modern and ancient Western Eurasians. RESULTS: Serbians share a number of U mtDNA lineages with Southern, Eastern-Central and North-Western Europeans. Putative Balkan-specific lineages (e.g. U1a1c2, U4c1b1, U5b3j, K1a4l and K1a13a1) and lineages shared among Serbians (South Slavs) and West and East Slavs were detected (e.g. U2e1b1, U2e2a1d, U4a2a, U4a2c, U4a2g1, U4d2b and U5b1a1). CONCLUSION: The exceptional diversity of maternal lineages found in Serbians may be associated with the genetic impact of both autochthonous pre-Slavic Balkan populations whose mtDNA gene pool was affected by migrations of various populations over time (e.g. Bronze Age pastoralists) and Slavic and Germanic newcomers in the early Middle Ages.

Massively parallel targeted resequencing reveals novel genetic variants associated with aspergillosis in paediatric patients with haematological malignancies.

This study aimed to find novel genetic variants of susceptibility to aspaergillosis in paediatric patients with haematological malignancies. Complete sequences of fifteen genes of human innate immunity (CCL2, CCR2, CD209, CLEC6A, CLEC7A and ten TLR genes) were studied in 40 patients diagnosed with haematological disorders (20 unaffected and 20 affected by aspergillosis). All samples were sequenced with MiSeq (Illumina) and 454 (Roche Diagnostics) technologies. Statistical significance of the differences between studied groups was determined using the two-tailed Fisher's exact test. Sixty variants of potential importance were identified, the vast majority of which are located in non-coding parts of the targeted genes. At the threshold of p < 0.000005, one intergenic (TLR2 rs4585282) and one intronic variant (CLEC6A rs12099687) were found significant between the case and control groups for genotype and allele frequencies, respectively. Rs12099687 in CLEC6A was predicted to constitute an alternative isoform or cryptic splice site, which potentially changes activity of the Dectin-2 protein. Overall, we assume that the two strongest associations reported in this study are expected to be reproducible even in the absence of other evidence, while another twelve associations may be strong enough to justify additional research in larger cohorts.

RNA isolation from bloodstains collected on FTA cards - application in clinical and forensic genetics.

Aim of the study: In recent years, RNA analysis has been increasingly used in clinical and forensic genetics. Nevertheless, a major limitation of RNA-based applications is very low RNA stability in biological material, due to the RNAse activity. This highlights the need for improving the methods of RNA collection and storage. Technological approaches such as FTA Classic Cards (Whatman) could provide a solution for the problem of RNA degradation. However, different methods of RNA isolation from FTA cards could have diverse effects on RNA quantity and quality. The purpose of this research was to analyze the utility of three different methods of RNA isolation from peripheral blood collected on FTA Classic Cards (Whatman). The study also aimed at assessing RNA stability in bloodstains deposited on FTA cards. Material and methods: The study was performed on peripheral bloodstains collected from 59 individuals on FTA Classic Cards (Whatman). RNA was isolated with High Pure RNA Isolation Kit (Roche Diagnostics), Universal RNA/miRNA Purification (EURx) and TRIzol Reagent (Life Technologies). RNA was subjected to quantitative analysis followed by reverse transcription and Real - Time PCR reaction. Results: The study has shown that FTA Classic Cards (Whatman) are useful tools for storing bloodstains at room temperature for RNA analysis. Moreover, the method of RNA extraction employing TRIzol Reagent (Life Technologies) provides the highest efficiency and reproducibility for samples stored for no more than 2 years. Conclusions: The FTA cards are suitable for collecting and storing bloodstains for RNA analysis in clinical and forensic genetics.

Searching for association of the CAG repeat polymorphism in the mitochondrial DNA polymerase gamma gene (POLG) with colorectal cancer.

Mitochondrial DNA polymerase gamma (POLG) is the only DNA polymerase involved in maintaining the mitochondrial genome. Recent studies demonstrated an association of CAG repeat polymorphism in the second exon of POLG gene with the risk of cancer. We investigated the CAG repeat variability in the POLG gene in tumor and non-tumor tissues from colorectal cancer patients and in DNA samples isolated from blood obtained from age-matched healthy persons. Somatically occuring CAG-repeat alterations in cancer tissues have been observed in 10% of patients, but no association has been found between the CAG repeat variants in the POLG gene and colorectal cancer risk.

Complete mitochondrial genome database and standardized classification system for Canis lupus familiaris.

To contribute to the complete mitogenome database of the species Canis lupus familiaris and shed more light on its origin, we have sequenced mitochondrial genomes of 120 modern dogs from worldwide populations. Together with all the previously published mitogenome sequences of acceptable quality, we have reconstructed a global phylogenetic tree of 555 C. l. familiaris mitogenomes and standardized haplogroup nomenclature. The phylogenetic tree presented here and available online at http://clf.mtdna.tree.cm.umk.pl/ could be further used by forensic and evolutionary geneticists as well cynologists, for data quality control and unambiguous haplogroup classification. Our in-depth phylogeographic analysis of all C. l. familiaris mitogenomes confirmed that domestic dogs may have originated in East Asia during the Mesolithic and Upper Paleolithic time periods and started to expand to other parts of the world during Neolithic times.