Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples.

BACKGROUND: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. RESULTS: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. CONCLUSIONS: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates.

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Human whole mitochondrial genome sequencing and analysis: optimization of the experimental workflow.

AIM: To evaluate critical steps in Illumina® Human mtDNA Genome assay: target enrichment, limited-cycle PCR, and library normalization, in order to optimize the protocol for analysis of whole mitochondrial genomes from human reference samples. METHODS: Three long-range high-fidelity DNA polymerases (PlatinumTM PCR SuperMix High Fidelity, LA Taq® Hot Start, and PrimeSTAR® GXL) were tested for their performance in the amplification of mtDNA fragments. Sequencing results of ten samples, as well as negative controls, which underwent library preparation with 12 and 15 cycles in limited-cycle PCR were compared. Additionally, two library normalization methods were compared: bead-based normalization vs quantification and individual normalization. RESULTS: PrimeSTAR® GXL performed best for mitochondrial DNA enrichment. Increment of amplification cycles to 15 in limited-cycle PCR step did not affect either the sequencing process or variant calling. Library quantification combined with individual library-by-library dilution outperformed bead-based normalization. CONCLUSION: Optimizations described herein provide beneficial insights for laboratories aiming at implementation and/or advancement of similar massively parallel sequencing workflows (eg, small genomes, PCR amplicons, and plasmids).

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines.

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Characterization of Nme5-Like Gene/Protein from the Red Alga Chondrus Crispus.

The Nme gene/protein family of nucleoside diphosphate kinases (NDPK) was originally named after its member Nm23-H1/Nme1, the first identified metastasis suppressor. Human Nme proteins are divided in two groups. They all possess nucleoside diphosphate kinase domain (NDK). Group I (Nme1-Nme4) display a single type NDK domain, whereas Group II (Nme5-Nme9) display a single or several different NDK domains, associated or not associated with extra-domains. Data strongly suggest that, unlike Group I, none of the members of Group II display measurable NDPK activity, although some of them autophosphorylate. The multimeric form is required for the NDPK activity. Group I proteins are known to multimerize, while there are no data on the multimerization of Group II proteins. The Group II ancestral type protein was shown to be conserved in several species from three eukaryotic supergroups. Here, we analysed the Nme protein from an early branching eukaryotic lineage, the red alga Chondrus crispus. We show that the ancestral type protein, unlike its human homologue, was fully functional multimeric NDPK with high affinity to various types of DNA and dispersed localization throughout the eukaryotic cell. Its overexpression inhibits both cell proliferation and the anchorage-independent growth of cells in soft agar but fails to deregulate cell apoptosis. We conclude that the ancestral gene has changed during eukaryotic evolution, possibly in correlation with the protein function.

Expanded Croatian 12 X-STR loci database with an overview of anomalous profiles.

In order to implement X-chromosome short tandem repeat (X-STR) typing into routine forensic practice, reference database of a given population should be established. Therefore we extended already published data with additional 397 blood samples from unrelated Croatian citizens, and analyzed the total of 995 samples (549 male and 446 female) typed by Investigator® Argus X-12 Kit. To test genetic homogeneity of consecutively processed five historic-cultural regions covering the entire national territory, we calculated pairwise Fst genetic distances between regions based on allele and full haplotype frequencies. Since the comparison did not yield any statistically significant difference, we integrated STR profile information from all regions and used the whole data set to calculate forensic parameters. The most informative marker is DXS10135 (polymorphism information content (PIC = 0.929) and the most informative linkage group (LG) is LG1 (PIC = 0.996). We confirmed linkage disequilibrium (LD) for seven marker pairs belonging to LG2, LG3 and LG4. By including LD information, we calculated cumulative power of discrimination that amounted to 0.999999999997 in females and 0.999999005 in males. We also compared Croatia with 13 European populations based on haplotype frequencies and detected no statistically significant Fst values after Bonferroni correction in any LG. Multi-dimensional scaling plot revealed tight grouping of four Croatian regions amongst populations of southern, central and northern Europe, with the exception of northern Croatia. In this study we gave the first extensive overview of aberrant profiles encountered during Investigator® Argus X-12 typing. We found ten profiles consistent with single locus duplication followed by tetranucleotide tract length polymorphism. Locus DXS10079 is by far the most frequently affected one, presumably mutated in eight samples. We also found four profiles consistent with X-chromosome aneuploidy (three profiles with XXX pattern and one profile with XXY pattern). In conclusion, we established integral forensic Croatian X-chromosome database, proved forensic pertinence of Investigator® Argus X-12 Kit for the entire Croatian population and identified locus DXS10079 as a potential duplication hotspot.

Investigator Argus X-12 study on the population of northern Croatia

Abstract X chromosome STR typing has emerged recently as a powerful tool, complementary to autosomal STR typing, in solving complex forensic and missing person cases. Investigator® Argus X-12 is a commercial product that allows co-amplification of 12 X chromosomal markers belonging to four linkage groups (LGs). In this study, we analyzed by capillary electrophoresis blood samples from 100 females and 102 males from a population of northern Croatia. Statistical analysis included calculation of allele and haplotype frequencies, as well as forensic parameters. The most informative marker for the northern Croatia population was DXS10135 with PIC=0.9211 and a total of 27 alleles. The least polymorphic marker was DXS8378 with 6 alleles. The proportion of observed haplotypes from the number of possible haplotypes varied from 2.74-8.57% across all LGs, with LG1 being the most informative. Of the 11 tested world populations compared to the population of northern Croatia, significant differences in genetic distance (FST) were found for Greenlandic and all non-European populations. We found that all tested markers are in HWE and can thus be used for match probability calculation. Because of high combined power of discrimination in both men and women, Investigator® Argus X-12 is applicable for the northern Croatia population in routine forensic casework.

Analysis of 12 X-STR loci in the population of south Croatia.

The aim of the study was to assess forensic pertinence of 12 short tandem repeats (STRs) on X-chromosome in south Croatia population. Investigator® Argus X-12 kit was used to co-amplify 12 STR loci belonging to four linkage groups (LGs) on X-chromosome in 99 male and 98 female DNA samples of unrelated donors. PCR products were analyzed by capillary electrophoresis. Population genetic and forensic parameters were calculated by the Arlequin and POPTREE2 software, and an on-line tool available at ChrX-STR.org. Hardy-Weinberg equilibrium was confirmed for all X-STR markers in female samples. Biallelic patterns at DXS10079 locus were detected in four male samples. Polymorphism information content for the most (DXS10135) and the least (DXS8378) informative markers was 0.9212 and 0.6347, respectively. In both male and female samples, combined power of discrimination exceeded 0.999999999. As confirmed by linkage disequilibrium test, significant association of marker pair DXS10074-DXS10079 (P = 0.0004) within LG2 and marker pair DXS10101-DXS10103 (P = 0.0003) within LG3 was found only in male samples. Number of observed haplotypes in our sample pool amounted 3.01, 7.53, 5 and 3.25% of the number of possible haplotypes for LG1, LG2, LG3 and LG4, respectively. According to haplotype diversity value of 0.9981, LG1 was the most informative. In comparison of south Croatia with 26 world populations, pair-wise [Formula: see text] values increase in parallel with geographical distance. Overall statistical assessment confirmed suitability of Investigator® Argus X-12 kit for forensic casework in both identification and familial testing in the population of south Croatia.

Investigator Argus X-12 study on the population of northern Croatia.

X chromosome STR typing has emerged recently as a powerful tool, complementary to autosomal STR typing, in solving complex forensic and missing person cases. Investigator® Argus X-12 is a commercial product that allows co-amplification of 12 X chromosomal markers belonging to four linkage groups (LGs). In this study, we analyzed by capillary electrophoresis blood samples from 100 females and 102 males from a population of northern Croatia. Statistical analysis included calculation of allele and haplotype frequencies, as well as forensic parameters. The most informative marker for the northern Croatia population was DXS10135 with PIC=0.9211 and a total of 27 alleles. The least polymorphic marker was DXS8378 with 6 alleles. The proportion of observed haplotypes from the number of possible haplotypes varied from 2.74-8.57% across all LGs, with LG1 being the most informative. Of the 11 tested world populations compared to the population of northern Croatia, significant differences in genetic distance (FST) were found for Greenlandic and all non-European populations. We found that all tested markers are in HWE and can thus be used for match probability calculation. Because of high combined power of discrimination in both men and women, Investigator® Argus X-12 is applicable for the northern Croatia population in routine forensic casework.

Analysis of 12 X-chromosomal markers in the population of central Croatia.

Investigator® Argus X-12 Kit is a commercially available set that allows simultaneous PCR amplification of 12 X-STR markers belonging to four linkage groups (LG). To assess the forensic efficiency of these markers for the population of central Croatia and consequent applicability in routine forensic casework, DNA from 200 blood samples of unrelated donors (100 female and 100 male) was amplified by Investigator® Argus X-12 Kit and analyzed by capillary electrophoresis. Statistical computations based on allele and haplotype frequencies for LG1 - LG4 were performed using Arlequin 3.5 software and on-line tool available at ChrX-STR.org. In female samples, all X-STR markers were in Hardy-Weinberg equilibrium (HWE). The most informative marker for central Croatia population was DXS10135 with polymorphism information content (PIC) 0.9296. The least polymorphic locus was DXS8378 (PIC=0.6363). Power of discrimination (PD) varied from 0.6968 to 0.9336 in male and from 0.8476 to 0.9916 in female samples. Combined PD exceeded 0.999999999 in both men and women. In male samples, linkage disequilibrium (LD) test revealed significant association (P=0.0000) of one marker pair in LG4 and two marker pairs in LG3. Portion of observed haplotypes in the number of possible haplotypes varied from 2.86% to 7.47% across all LGs. LG1 was the most informative with haplotype diversity (H) 0.9972. High PD of all analyzed markers exhibited for central Croatia population confirms suitability of Investigator® Argus X-12 for forensic pertinence. Moreover, results of this study will be included in establishing a national reference X-STR database based on 12 X-STR loci, which is necessary for the correct interpretation of the forensic casework results.

Response to hyperoxia is associated with similar ho-1 gene expression level in lungs of aging CBA mice of both sexes.

Beside classical antioxidative enzymes, the response to hyperoxia might be mediated via regulation of other systems, such as heme oxygenase (HO). Ho-1 gene expression is found to be upregulated by hyperoxia in all groups of mice, while HO-1 protein isoform was increased only in 4 months old male mice. In steady-state conditions ho-1 and ho-2 gene expression remained unchanged irrespective of sex or age, which was not the case with protein level of both isoforms. This study suggests that in lungs of CBA mice the response to oxidative stress may be mediated through the interaction of other systems such as heme oxygenase, primarily via upregulation of ho-1 gene expression in both sexes. Contrary to our previous study in liver of hyperoxia treated mice, current results might imply that at conventional oxygen conditions lungs of female mice with the emphasis on aging females, are better prepared for oxidative stress conditions through the increase of HO-activity.

Functional and Structural Characterization of FAU Gene/Protein from Marine Sponge Suberites domuncula.

Finkel-Biskis-Reilly murine sarcoma virus (FBR-MuSV) ubiquitously expressed (FAU) gene is down-regulated in human prostate, breast and ovarian cancers. Moreover, its dysregulation is associated with poor prognosis in breast cancer. Sponges (Porifera) are animals without tissues which branched off first from the common ancestor of all metazoans. A large majority of genes implicated in human cancers have their homologues in the sponge genome. Our study suggests that FAU gene from the sponge Suberites domuncula reflects characteristics of the FAU gene from the metazoan ancestor, which have changed only slightly during the course of animal evolution. We found pro-apoptotic activity of sponge FAU protein. The same as its human homologue, sponge FAU increases apoptosis in human HEK293T cells. This indicates that the biological functions of FAU, usually associated with "higher" metazoans, particularly in cancer etiology, possess a biochemical background established early in metazoan evolution. The ancestor of all animals possibly possessed FAU protein with the structure and function similar to evolutionarily more recent versions of the protein, even before the appearance of true tissues and the origin of tumors and metastasis. It provides an opportunity to use pre-bilaterian animals as a simpler model for studying complex interactions in human cancerogenesis.

Female headstart in resistance to hyperoxia-induced oxidative stress in mice.

Increased oxygen concentration (hyperoxia) induces oxidative damage of tissues and organs. Oxygen toxicity in hyperoxia is controlled by factors such as sex, age, tissue, strain and hormones. In most species females show lower incidence of some age-related pathologies linked with oxidative stress, which has been attributed to a beneficial effect of ovarian hormones. In this study we found that hyperoxia induced hepatic oxidative damage exclusively in male CBA/H mice, followed by their decreased survival. Histopathological examination revealed that the observed differences in survival were not the consequence of acute lung injury induced by hyperoxia. Next, we observed that an increased Sirt1 protein level in hyperoxia-exposed female CBA/H mice correlated with their lower PPAR-γ and higher eNOS and Sod2 protein levels. In males, higher PPAR-γ and lower Sod2 protein levels were associated with unchanged Sirt1 expression. Although these results are of a correlative nature only, they clearly show that females show better survival, increased resistance to hyperoxia and have generally more efficient defense systems, which suggests that their headstart in resistance to hyperoxia could be a consequence of the beneficial effect of ovarian hormones.

Polymorphisms in the IL-18 and IL-12B genes and their association with the clinical outcome in Croatian patients with Type 1 diabetes.

Genetic variants of IL-18 and IL-12B may be important in immunoregulatory abnormalities, observed in the patients with Type 1 diabetes mellitus (T1DM), that contribute to individual differences in response to a treatment. Therefore, we examined the significance of IL-18-137G/C, IL-18-607C/A, and IL-12B A/C polymorphisms in Croatians (187 patients, 236 controls), not only as factors that contribute to susceptibility to T1DM, but also as determinants of the clinical presentation of disease. The polymorphism screening has been performed using PCR sequence-specific primers (IL-18) or PCR-RFLP (IL-12B) approach. Results were evaluated by GraphPad Prism and Sigma Stat 3.5, Arlequin software and calculator for Hardy-Weinberg equilibrium. The genotype, allele and haplotype distribution were not statistically different between the patients and control subjects. The clinical parameter analysis revealed that patients with minor alleles at each locus, IL-18-137C/-607A, were significantly younger at T1DM onset than carriers of major alleles, IL-18-137G/-607C (20 vs 23.5 years). Moreover, the concomitant presence of minor alleles not only of IL-18 but also of IL-12B, is associated with the risk of disease progression even at younger age. These patients developed diabetes at 16 years of age, what is significantly earlier (p=0.044) compared to 25.5 years of age in patients with common alleles IL-18-137G/-607C/IL-12B A. Furthermore, combined genotype analysis of IL-18 and IL-12B has pointed out that patients with CC/AA/AA genotype have the worst glucose control based on HbA1c (8.7%, range 6.8-13.1%). In conclusion, susceptibility to T1DM in Croatians is not strongly associated with IL-18-137/-607 and IL-12B polymorphisms. These SNPs are associated with the higher risk of earlier disease development and might be implicated in the effectiveness of glycemic control.

Structural and functional characterization of ribosomal protein gene introns in sponges.

Ribosomal protein genes (RPGs) are a powerful tool for studying intron evolution. They exist in all three domains of life and are much conserved. Accumulating genomic data suggest that RPG introns in many organisms abound with non-protein-coding-RNAs (ncRNAs). These ancient ncRNAs are small nucleolar RNAs (snoRNAs) essential for ribosome assembly. They are also mobile genetic elements and therefore probably important in diversification and enrichment of transcriptomes through various mechanisms such as intron/exon gain/loss. snoRNAs in basal metazoans are poorly characterized. We examined 449 RPG introns, in total, from four demosponges: Amphimedon queenslandica, Suberites domuncula, Suberites ficus and Suberites pagurorum and showed that RPG introns from A. queenslandica share position conservancy and some structural similarity with "higher" metazoans. Moreover, our study indicates that mobile element insertions play an important role in the evolution of their size. In four sponges 51 snoRNAs were identified. The analysis showed discrepancies between the snoRNA pools of orthologous RPG introns between S. domuncula and A. queenslandica. Furthermore, these two sponges show as much conservancy of RPG intron positions between each other as between themselves and human. Sponges from the Suberites genus show consistency in RPG intron position conservation. However, significant differences in some of the orthologous RPG introns of closely related sponges were observed. This indicates that RPG introns are dynamic even on these shorter evolutionary time scales.

Alzheimer's disease and type 2 diabetes: the association study of polymorphisms in tumor necrosis factor-alpha and apolipoprotein E genes.

Type 2 diabetes (T2D) and Alzheimer's disease (AD) are two progressive disorders with high prevalence worldwide. Polymorphisms in tumor necrosis factor-alpha (TNF-α) and apolipoprotein E (ApoE) genes might be associated with both T2D and AD, representing possible genetic markers for the development of the AD in subjects with T2D. The aim was to determine ApoE and G-308A TNF-α gene polymorphisms in unrelated Croatian Caucasians: 207 patients with sporadic AD, 196 T2D patients and 456 healthy controls. Patients with AD had higher frequency of ApoE4 allele compared to T2D patients and controls. The significant association, observed between ApoE2 allele and T2D, disappeared after the data were adjusted for age and sex. The genotype or allele frequencies of G-308A TNF-α gene polymorphism were similar among the patients with AD, T2D and healthy controls. In conclusion, these results do not support the hypothesis that the A allele of G-308A TNF-α gene polymorphism is associated either with AD or T2D. Our data confirm the association between the ApoE4 allele and AD, and point out the E2 allele of ApoE gene as the possible risk factor for T2D.

Rational variety mapping for contrast-enhanced nonlinear unsupervised segmentation of multispectral images of unstained specimen.

A methodology is proposed for nonlinear contrast-enhanced unsupervised segmentation of multispectral (color) microscopy images of principally unstained specimens. The methodology exploits spectral diversity and spatial sparseness to find anatomical differences between materials (cells, nuclei, and background) present in the image. It consists of rth-order rational variety mapping (RVM) followed by matrix/tensor factorization. Sparseness constraint implies duality between nonlinear unsupervised segmentation and multiclass pattern assignment problems. Classes not linearly separable in the original input space become separable with high probability in the higher-dimensional mapped space. Hence, RVM mapping has two advantages: it takes implicitly into account nonlinearities present in the image (ie, they are not required to be known) and it increases spectral diversity (ie, contrast) between materials, due to increased dimensionality of the mapped space. This is expected to improve performance of systems for automated classification and analysis of microscopic histopathological images. The methodology was validated using RVM of the second and third orders of the experimental multispectral microscopy images of unstained sciatic nerve fibers (nervus ischiadicus) and of unstained white pulp in the spleen tissue, compared with a manually defined ground truth labeled by two trained pathophysiologists. The methodology can also be useful for additional contrast enhancement of images of stained specimens.

Over-represented localized sequence motifs in ribosomal protein gene promoters of basal metazoans.

Equimolecular presence of ribosomal proteins (RPs) in the cell is needed for ribosome assembly and is achieved by synchronized expression of ribosomal protein genes (RPGs) with promoters of similar strengths. Over-represented motifs of RPG promoter regions are identified as targets for specific transcription factors. Unlike RPs, those motifs are not conserved between mammals, drosophila, and yeast. We analyzed RPGs proximal promoter regions of three basal metazoans with sequenced genomes: sponge, cnidarian, and placozoan and found common features, such as 5'-terminal oligopyrimidine tracts and TATA-boxes. Furthermore, we identified over-represented motifs, some of which displayed the highest similarity to motifs abundant in human RPG promoters and not present in Drosophila or yeast. Our results indicate that humans over-represented motifs, as well as corresponding domains of transcription factors, were established very early in metazoan evolution. The fast evolving nature of RPGs regulatory network leads to formation of other, lineage specific, over-represented motifs.