Genes

Environmental and dietary factors can exert multigenerational effects on health and development. In this study, we investigated whether early-life metabolic challenge affects the germline genome and epigenome across three generations. Using a murine model of early life obesity via litter size reduction (overnutrition group, ON) and a control group (CT), we followed the paternal lineage focusing on germline genomic and methylation changes employing Genotyping-by-Sequencing (GBS) coupled with methyl-immunoprecipitation (GBS-MeDIP). We found that unrelated ON families clustered together based on identified Single-Nucleotide Polymorphism (SNP), suggesting that the treatment may have genomic impact. Copy number variations (CNVs) events were identified in ON individuals, being enriched in Long Interspersed Nuclear Elements (LINEs) and Long Terminal Repeats (LTRs). While Principal Component Analysis (PCA) of the methylome showed no clear treatment effect, pathway enrichment and regional analyses revealed methylation changes associated with transposable elements and developmental genes. Notably, the ON group exhibited a disruption in the methylation of Repetitive Elements (RE), which was significant in the same type of RE that were also enriched in the observed CNVs. The ON also showed reduced emergence of novel SNPs in offspring compared to the CT group. These findings suggest that multigenerational metabolic challenge can constrain genetic variability and induce genome instability, potentially mediated by transposable element activity rather than widespread changes in DNA methylation. This work highlights the importance of studying both genome and epigenome dynamics under realistic, multigenerational exposure scenarios and suggests that early metabolic challenges can have long-lasting impacts on genomic architecture and evolutionary potential.

Mitochondrial diseases (MDs) are clinically heterogeneous rendering ascertainment challenging. Estimates of pathogenic mitochondrial DNA (mtDNA) variants in the population range from 1 in 200 to 1 in 4,000 individuals. Inclusion of mtDNA sequencing in genomic databases facilitates comprehensive estimation of mtDNA variation. However, interpretation of low heteroplasmy variation is complex, due in part to misalignment of nuclear mitochondrial DNA transcripts (NUMTs), whilst conservative heteroplasmy thresholds likely omit relevant variation. Cumulative burden of mtDNA variation contributes to aging and neurodegeneration, and recent characterisation of mitochondrial genome constraint allows quantitation of this burden. We analysed whole genome sequencing of blood DNA from 3,500 healthy older individuals in the Medical Genome Reference Bank using mity, considering pathogenic mtDNA variants [≥]1% heteroplasmy. We identified 34 distinct pathogenic mtDNA variants in 62 individuals, giving a combined population allele frequency of 1.77% (95% CI 1.36-2.27) or 1 in 56 individuals. We evaluated inclusion of false positive (FP) calls due to two common NUMTs, which accounted for up to 16% of variants. Increasing heteroplasmy thresholding to eliminate all NUMT-FPs also eliminated much of the total variation, including pathogenic variants. We propose a sample-specific, scaled heteroplasmy threshold to maximise variant retention and mitigate NUMT-FPs. Finally, we characterised measures of mitochondrial constraint in this healthy older cohort, observing an association between variant burden and summed constraint, whilst mean constraint was higher in pathogenic variant carriers. These findings suggest pathogenic mtDNA variation is more common in the population than is currently appreciated. Findings are comparable to larger genomic databases when heteroplasmy thresholding is adjusted, and support earlier population-based estimates. Incorporation of low heteroplasmy variation is relevant, but interpretation is nuanced, and optimising variant retention requires consideration of NUMT-FP rates.

Complex rearrangements are one of the rarest types of structural variants (SVs) and can be divided into two categories: complex chromosomal rearrangements (CCRs) and complex genomic rearrangements (CGRs). CCRs include structural rearrangements that present at least three breakpoints and show exchange of genetic material between more than two chromosomes and CGRs are rearrangements that present more than one junction and/or more than one SV in cis. They are usually formed by one of the chromoanagenesis mechanisms, where a massive disruptive cellular event leads to multiple structural rearrangements. Classical cytogenomic techniques have been commonly applied for their characterization, but methodologies that involve longer DNA molecules, namely optical genome mapping (OGM) and long-read genome sequencing (lrGS), present a considerably higher SV detection resolution, revealing more details about the rearrangements, including precise breakpoint location. Here, we describe six patients with complex rearrangements investigated through a combination of different techniques: karyotyping, chromosomal microarray, and OGM were performed to characterize the rearrangements. Subsequently, lrGS was used to further resolve the alterations, refine their breakpoints' location, and sequence their junction points. Three patients presented CCRs involving three, four, and six chromosomes, while three exhibited CGRs involving one different chromosome each, providing a variety of complex SVs to show the importance of each technique and their combination in rearrangement resolution. In total, the complex rearrangements presented 127 breakpoints, 66 junction points and involved 14 of the 24 chromosomes. Higher-resolution techniques revealed additional complexity in all cases. Despite the advances provided by OGM and lrGS, conventional karyotyping remained indispensable for complete rearrangement resolution. In two patients, the findings supported a novel mechanism combining features of the different chromoanagenesis processes. Furthermore, evidence of inherited alterations was identified, and the comprehensive characterization of the rearrangements enabled more accurate genotype-phenotype correlations. Our findings indicate that an integrated approach combining karyotyping, OGM, and lrGS can completely resolve SVs, including complex rearrangements.

The early onset of labour increases mortality and developmental risks for a human newborn. Key genes in human labour have been investigated using multiple modalities, but their regulation by non-coding RNA (e.g. lncRNA and miRNA) remains incomplete. This study explores the three-way relationship between labour-associated transcription factors (TFs), miRNA and lncRNA suggested by the competing endogenous RNA (ceRNA) hypothesis, to understand the underlying regulatory framework. Experimentally validated miRNA-lncRNA interactions are modelled using five distinct machine learning (ML) architectures to predict 20469 labour-linked miRNA-lncRNA interactions. Known mRNA-ncRNA interactions from databases were included to construct a tripartite network, and a subset of 9989 labour-linked network motifs containing TFs were isolated and analysed. Gene enrichment of nodes in TF-lncRNA-miRNA network, as well as validation from public myometrial datasets indicate high significance in contractile pathways including immune signalling. Experimentally unconfirmed tripartite network motifs have been found, and we elaborate on their potential regulation in labour using 8 TF-lncRNA-miRNA network motifs. A unified ncRNA-TF regulatory atlas in labour has been synthesized, and a complete summary of the tripartite network motifs can be accessed and visualised using the user-friendly, public database.

BackgroundHeterochromatic long arm of mouse Y chromosome harbors the multicopy species-specific sequences Ssty, Sly, Asty and Orly that are transcribed in testis and have known functions in male fertility. Of these Ssty and Sly encode proteins - yet all the transcripts are not translated. To investigate the roles of these Y-heterochromatic transcripts further, we analyzed them. MethodsMice with 2/3rd deletion of the Y-chromosome (XYRIIIqdel) and its wild type (XYRIII) were used in this study. Bioinformatic approaches, small RNA northern blots, Electrophoretic Mobility Shift Assays, Luciferase reporter assays, dPCR analysis, RT-qPCR assays and western blotting techniques were used to identify piRNAs that regulate autosomal genes. ResultsWe demonstrate that the multicopy gene families from mouse Y-long arm generate piRNAs predominantly in testis. We observed sequences homologous to these piRNAs in the UTRs of a few autosomal genes, which are differentially expressed in the sperms of XYRIIIqdel mice. Furthermore, the Endogenous Retrovirus Element (ERV) LTR, found in the Orly1 transcript identified piRNAs in the database, showed homology to UTRs and associated genomic regions of a few autosomal genes. Orly1 showed a reduction in genomic copy number by digital PCR in XYRIIIqdel mice. One of the four autosomal genes containing the ERV segment in their UTRs, showed a differential testicular protein expression in the mutant mice. ConclusionsThus, we further elucidate that different classes of repeats from Y-chromosome regulate autosomal gene expression via piRNAs. Besides, this study also identified novel roles for a Y-derived ERV in autosomal gene regulation in testis.

As genetic testing becomes increasingly integrated into Parkinson disease (PD) research, including targeted testing for variants in LRRK2 and GBA1, the return of individual research results is becoming more common. However, limited qualitative data exists regarding how research participants experience genetic results disclosure and post-test genetic counseling in PD research settings. We conducted semi-structured qualitative interviews with participants (n=13) enrolled in the Parkinson Precision Medicine Initiative (formerly Parkinson Progression Markers Initiative; PPMI) who had received PD-related genetic test results and post-test genetic counseling. Interviews were conducted 1 to 3 weeks following result disclosure and analyzed using thematic analysis with a primarily deductive coding approach informed by study aims and inductive identification of emergent themes. Four primary themes were identified: (1) personal connection and motivations for participation, (2) centrality of result disclosure and information preferences, (3) emotional experiences and support needs, and (4) communication quality and alignment with participant needs. Overall, our findings underscore the importance of person-centered genetic counseling within PD research. As return of genetic and biomarker results in research and clinical trial contexts expand, thoughtful integration of relational, informational, and communication-focused practices will be essential to support participant engagement and trust.

Genome sequencing of the heterogeneous primary mitochondrial disorders (PMD) frequently reveals variants of uncertain significance that require functional tests for diagnosis, and does not identify variants in all patients. We analyzed mitochondrial enzyme assays, blue native polyacrylamide gel electrophoresis (BN-PAGE) with in-gel activity staining, complex I assembly blot, and select protein abundances in fibroblasts of a case series of 204 PMD patients divided into functional classes, in comparison to 51 controls and 53 differential diagnostic conditions. Overall, sensitivity and specificity for respiratory chain enzyme assays were 46% and 93% respectively, for BN-PAGE 40% and 98%, for complex I assembly assay 49% and 99%. The overall sensitivity of all tests was 76%, specificity 93%, with positive predictive value 96% and negative predictive value 67%. Categories with high sensitivity were isolated complex deficiencies, nuclear DNA-encoded mitochondrial protein synthesis defects, co-factor defects, and mitochondrial amino-acyl-tRNA synthetase conditions when aided by protein abundance. Mitochondrial DNA mutations and maintenance disorders showed poor sensitivities. Secondary dysfunctions were rare. A complete battery of functional tests showed strong diagnostic clinical utility in fibroblasts.

The leading cause of mortality and morbidity in children under the age of 5 is preterm birth. The timing of birth is influenced by both genetic and environmental factors, but the underlying mechanisms remain poorly understood, making its prediction difficult. In this study, we investigated the potential of using machine learning models to predict preterm birth based on genetic data from the Norwegian Mother, Father and Child Cohort Study (MoBa). We trained and evaluated several classification algorithms on individual-level genetic data from over 15,000 mothers and children. Our results indicate that the predictive capacity of maternal gestational duration-associated loci for preterm birth is limited, with the highest AUC values around 0.57. Additionally, incorporating more SNPs within the associated loci did not improve prediction performance. As expected, the contribution of the maternal genome to preterm birth prediction was found to be larger than that of the fetal genome. Overall, our findings suggest that while genetic testing provides some information about an individual's risk for preterm birth, further research incorporating additional factors is necessary to enhance predictability.

Advances in high-throughput sequencing and genetic research have expanded the role of genetics in medicine and society. Population-based screening programs, including neonatal and preconception testing, are increasingly implemented globally, alongside the rise of direct-to-consumer (DTC) genetic testing. The "Genetics and the General Public" Ethics Working Group of the French Federation of Human Genetics (FFGH) assessed knowledge and awareness of genetics within the French population through a nationally representative survey (n=3,013) conducted by the polling firm Ipsos bva. Results indicated that 69% of respondents report an interest in genetics, although their level of knowledge remains limited. Most respondents expressed positive attitudes toward genetics, perceiving it as a major source of hope in healthcare. While a majority indicated willingness to undergo genetic testing for medical purposes, they also reported legitimate concerns regarding the potential results. Despite legal restrictions, 12% reported having ordered a DTC genetic test (5% for genealogical; 5% for medical and 2% for both purposes), and 45% of non-users expressed strong interest in this type of test. Notably, there is a substantial lack of awareness regarding the limitations of these tests and the French legal framework governing their use. These findings highlight critical gaps in public knowledge, emphasizing the need for improved genetic education, including incorporating genetics into school curricula and launching targeted awareness campaigns. These initiatives should help clarify the distinctions between clinically validated genetic tests and DTC genetic testing services, addressing both their benefits and their ethical, legal, and scientific limitations, in order to promote informed decision-making.

BackgroundGray is one of the relatively rare coat colors in donkeys. The Hetian Gray donkey is a distinctive indigenous breed from the Xinjiang Uygur Autonomous Region of Northwestern China, characterized by progressive hair depigmentation with aging while retaining dark skin pigmentation. However, the genetic basis underlying this unique gray coat color phenotype remains unclear. ResultsTo elucidate the genetic basis, we conducted whole-genome resequencing of Gray and non-Gray donkeys. Genome-wide selection signature analyses identified a candidate region on chromosome 15. Subsequent fine-mapping using mass spectrometry-based genotyping of 42 loci refined the candidate interval and revealed a SNP within intron 2 of the ASIP gene, located in a genomic fragment with highly similar sequences, showing complete association with the gray coat color. Association analysis in an expanded population further confirmed a strong correlation between this variant and the gray phenotype. Gene expression analyses also supported the role of ASIP in regulating pigmentation in donkeys. ConclusionsThese findings identify a genetic determinant of gray coat color in donkeys and provide new insights into the molecular mechanisms underlying age-related depigmentation in domestic animals.

Long COVID (LC) is a complex condition characterized by persistent, chronic multisystem manifestations, with a significant proportion of patients exhibiting neurological symptoms. Human ion channels (HICs), particularly potassium channels, are abundantly expressed in the nervous system and linked to key metabolic processes, making them potential candidates for understanding LC pathophysiology and drug repurposing. Meta-analysis of RNA-Seq datasets from COVID-19 recovered and LC patients was performed to identify altered HICs in LC. Differential gene expression analysis, functional enrichment analysis, and weighted gene co-expression network analysis (WGCNA) were performed to uncover key genes, pathways, and co-expression modules consisting of HICs, lipid metabolism-, and immune signaling-related genes. Drug-gene interaction analysis was performed to identify approved drugs targeting potential HICs. A total of 715 dysregulated genes, including eighteen HICs were identified, among which seven were potassium channels. Three significant modules containing HICs, lipid metabolism-, and immune signaling-related genes were identified and found to be associated with antigen processing and presentation, complement and coagulation cascades, and cytokine-related pathways. Approved drugs targeting KCNA6, KCNJ10, KCNN3, and KCNH4 were identified. With further experimental validation, these dysregulated potassium channels, supported by their co-expression networks and pathway associations, may act as potential candidates for drug repurposing in LC patients.

Purpose: To identify transcript-level variants in crystallin genes in paediatric patients with unilateral cataracts. Methods: Anterior capsulorhexis (n=12) from patients underwent surgical management of congenital unilateral cataracts was collected. Total RNA was isolated from lens epithelial cells, and complementary DNA (cDNA) was synthesized. Full-length RNA transcripts of 10 lens-specific crystallin genes were PCR-amplified and analysed via Sanger sequencing. Identified transcript variants were further validated using genomic DNA (gDNA) through Sanger sequencing. In addition, the full-length (~7,535 bp) CRYBA1 genomic region was sequenced using Oxford Nanopore Technology. Results: Aberrant low molecular weight (LMW) amplicons (~370 bp) of the CRYBA1 transcript were identified in three patients presented with unilateral cataract. Of 3 patients, 2 had persistent fetal vasculature (PFV) and 1 had pre-existing posterior capsular defect (PPCD). Sanger sequencing revealed a precise loss of exons 2 to 4 in the CRYBA1 RNA transcript. No coding, splice-site, or large deletion variants were detected in the genomic DNA of the patients or their parents. In silico analysis predicted two possible truncated proteins arising from these alternatively spliced transcripts: one comprising the first 11 amino acids of the N-terminal region with a loss of all Greek key motifs, and another comprising 90 amino acids encoded by exons 5 and 6, initiated from an alternative start codon in exon 5, and loss of Greek key motifs 1 & 2. Conclusion: The precise skipping of exons 2 to 4, consistent with canonical splicing signals (5-prime-GU...AG-3-prime), in the absence of genomic alterations, suggests the presence of alternatively spliced (AS) CRYBA1 transcripts in human lenses. This is the first report documenting AS-CRYBA1 transcripts in association with childhood cataracts with PFV and PPCD.

Background: It is well known that genetic variants contribute to cellular sensitivity to chemotherapeutic agents and ionizing radiation (IR). The aim of this study was to identify single nucleotide polymorphisms (SNPs) and genes associated with the spectrum of normal cellular sensitivity of lymphoblastoid cell lines (LCLs) towards ionizing radiation and mitomycin C (MMC). Methods: In a first step, we determined the viability of LCLs established from male participants of the Berlin Aging Study II (BASE-II) aged >=62 years following treatments with increasing doses of IR (n=137 cell lines) or MMC (n=140 cell lines) using the alamarBlue assay. Results from intra-experimental triplicates and three independent experiments for each cell line and treatment were used to calculate the area under the curves (AUCs) representing the specific sensitivity to IR and MMC of each LCL. The data from these experiments were subsequently used as outcomes in genome-wide association studies (GWASs). In addition, we calculated polygenic risk scores (PGS) from UK Biobank GWAS results for four cancer-related phenotypes and assessed the extent to which the variance in the IR and MMC sensitivity is explained by these PGS. Results: The GWAS analyses revealed one variant, rs74728080, located in CDH13 on chromosome 16, to show genome-wide significant (p < 5 x 10-8, beta = 2.81) association with cellular viability after treatment with IR. In the GWAS on MMC sensitivity the most interesting signal was elicited by SNP rs113978558 in an intron of the PLD5 gene on chromosome 1 (p = 9.232 x 10-8; beta = 1.44). Several other SNPs with statistically suggestive (i.e., p < 1 x 10-5) evidence of association with IR or MMC sensitivity were identified. PGSs calculations from GWAS of four cancer-related traits in UKB explained ~5% and ~3% of phenotypic variance in IR- and MMC-induced cell viability, respectively. Conclusion: The genome-wide significant association of rs74728080 with IR sensitivity and the location of this variant in CDH13 is interesting and functionally highly plausible given its known involvement in oxidative-stress response and function as tumor suppressor. Taken together, our novel data suggest that CDH13 may be genuinely involved in regulating cellular IR sensitivity.

The Gene Version Iteration Hypothesis (GVIH) proposes that mutant genes may originate from the Y chromosome, traverse through the X chromosome to autosomes, undergo interchromosomal transfer, and potentially return to the Y chromosome via the X chromosome. This hypothetical closed transmission loop may facilitate the storage, screening, and elimination of different versions of mutant genes. The hypothesis comprises five core propositions: (1) Mutation reservoir: The Y chromosome may serve as a specialized carrier for generating mutant genes, characterized by elevated mutation rates, reduced gene density, and accelerated evolutionary dynamics; (2) Closed-loop transmission: Mutant genes may follow a unidirectional pathway Y[-&gt;]X[-&gt;]autosomes[-&gt;]X[-&gt;]Y, forming a complete transmission circuit; (3) Coexistence of multiple versions: A single functional gene may exist in multiple versions across different chromosomes, constituting a dynamic gene version library; (4) Reproductive screening: Environmentally adaptive gene versions may persist across generations and potentially migrate to upstream chromosomes, while maladaptive versions may be eliminated; (5) Terminal elimination: Gene versions reaching the Y chromosome may undergo elimination processes, potentially preventing version monopolization and maintaining evolutionary dynamics. This hypothesis provides a novel framework for understanding adaptive evolution at the genetic level. If empirically validated, it may offer new insights into the molecular mechanisms underlying certain genetic phenomena and evolutionary processes.

Holocentric chromosomes allow rapid genome changes through chromosomal rearrangements such as fissions, fusions, inversions or translocations. The plant genus Carex shows one of the highest rates of karyotypic evolution among holocentric organisms. We studied the genomic patterns underlying chromosomal rearrangements in the karyotypic radiation of the narrow endemic species Carex helodes (2n = 68-75). Comparing genome assemblies of C. helodes from the two karyologically distinct extremes of its European distribution, revealed a striking number of eight chromosomal rearrangements including fusions, translocations and inversions. Genomic breakpoints are gene-poor and TE-rich, corroborating findings in other species and suggesting common genomic characteristics that facilitate the evolution and establishment of chromosomal rearrangements. We identified a chromosomal inversion exhibiting patterns of purifying selection and enrichment in functional genes that potentially mediate rearrangement tolerance. Conversely, another inversion displayed elevated sequence divergence and enrichment in response to temperature stress and phosphate limitation, matching key environmental variables that differ between the study localities. The establishment of chromosomal rearrangements along Carex helodes European populations was likely driven by demographic bottlenecks and distinct genomic features at breakpoints. Our findings provide preliminary evidence on the rearrangement role in population differentiation either as reproductive barriers or as genomic islands of differentiation.

Hurthle cell carcinoma (HCC) is an aggressive form of thyroid cancer. While mitochondrial DNA mutations and chromosomal losses have been identified in HCC, isoform switching, and its functional consequences remain uncharacterized. This study reanalyzed NCBI GEO dataset GSE228870 (n = 32), using Salmon and IsoformSwitchAnalyzeR() to identify isoform switching. The analysis resulted in 371 switches across 335 genes showing functional consequences including loss of protein domains, shorter open reading frames (ORFs), loss of signal peptides and novel sub-cellular localizations. Most significant isoform switches (q-value < 0.05, |dIF| > 0.1) were observed in LAMA2, LSP1, MAD2L2, FBLN2 and CXCL12, implicating extracellular matrix dysregulation, DNA damage response, immune signaling and cytoskeleton regulation. These genes are expressed in normal thyroid (median TPM 20.69, 11.66, 14.79, 134.1 & 80.76). However, specific isoforms of LAMA2 and MAD2L2 are not expressed in normal thyroid, explaining tumor-specific expression in HCC. Alternative transcription termination site (ATTS) gain was significant, suggesting altered 3 end in HCC transcripts. TCGA SpliceSeq showed LSP1, FBLN2 and CXCL12 undergo alternative promoter (LSP1 exon1 PSI=94.5%, FBLN2 exon2 PSI=99.0%) and alternative termination (CXCL12 exon3.3 PSI=53.9%) in thyroid cancer, suggesting ATTS and alternative transcription start site (ATSS) as shared splicing dysregulation mechanisms. This is the first systematic characterization of isoform-level dysregulation in HCC.

The AMY2B gene encodes pancreatic amylase, a critical enzyme for starch digestion. While previous studies have examined AMY2B copy number variation (CNV) in domestic and some wild animals, less is known about wild carnivores inhabiting regions with limited anthropogenic starch exposure. We analyzed blood samples for serum amylase activity and copy number variation in AMY2B gene from 8 wolves (Canis lupus), 11 brown bears (Ursus arctos), and 3 red foxes (Vulpes vulpes) from Bosnia and Herzegovina. AMY2B gene copy number was assessed using droplet digital PCR (ddPCR), and serum amylase activity and glucose levels were quantified. Although the number of fox samples was limited, foxes and wolves consistently harbored two copies of AMY2B, while brown bears exhibited higher CNV (3.67-8.40, mean 5.88). Serum amylase activity was highest in foxes, moderate in wolves, and variable but lower in bears. Despite differences in AMY2B copy number and serum amylase activity, circulating glucose concentrations did not differ significantly among species. Our findings suggest that variation in AMY2B copy number among wild carnivores may be associated with species-specific evolutionary histories and dietary adaptations, providing insight into genomic mechanisms underlying carbohydrate utilization in natural populations.

Neurodegenerative diseases are complex disorders characterised by progressive neuronal loss and widespread transcriptomic dysregulation; however, the coordinated interactions among coding and non-coding RNAs that contribute to disease progression remain incompletely understood. In this study, RNA-seq datasets from disease-relevant neuronal populations and brain regions representing Alzheimer's disease (AD), Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS) were analysed using an integrative network-based framework. Differential expression analysis coupled with weighted gene co-expression network analysis identified modules significantly correlated with disease and prioritised highly connected hub genes. Integration of these hub genes with curated RNA interaction database enabled the construction of candidate lncRNA-miRNA-mRNA regulatory networks. Functional enrichment analysis revealed Gene Ontology biological processes associated with synaptic signalling, mitochondrial function, RNA metabolism and neuroinflammatory responses across neurodegenerative conditions. The inferred regulatory networks suggested both disease-specific and shared post-transcriptional regulatory modules involving key hub genes and non-coding RNAs. Additionally, putative sequence variants were identified within untranslated regions of selected hub genes, suggesting potential alterations in miRNA-mediated regulations. Therefore, this study provides a systems-level view of transcriptomic dysregulation across major neurodegenerative diseases and identifies candidate regulatory interactions and molecular targets for future functional investigation

AnnotationMitochondrial heteroplasmy --the coexistence of both wild-type and mutant copies of mitochondrial DNA (mtDNA) within a cell--is a key factor in the pathogenesis of mitochondrial diseases. Classical approaches, which rely solely on the scalar fraction of mutant DNA, fail to fully account for threshold effects, the stochastic nature of heteroplasmy dynamics, and tissue specificity. The aim of the work is to construct a complex stochastic model of heteroplasmy dynamics, which for the first time combines the effects of selection, genetic drift, migration of mitochondrial genomes between tissues and threshold mechanisms of pathology development, for a quantitative assessment of the risk of mitochondrial diseases. In this paper, we propose a complex-phase formalism in which the state of a cells mitochondrial genome is described by a complex number Z = a + ib, where a and b are the absolute numbers of normal and mutant mtDNA copies, respectively. This approach naturally combines information on copy number and heteroplasmy level, and the argument{phi} = arctan (b / a) is interpreted as a phase characterizing the mutant load. Based on this formalism, we developed a stochastic model of tissue dynamics that includes the processes of selection, genetic drift, and intertissue migration of mitochondrial genomes. Using Monte Carlo methods (1000 simulations), we demonstrated that neuronal tissues are characterized by high heteroplasmy variability and a significant probability of reaching a pathological threshold even with a relatively low systemic mutant load. Kaplan-Meier survival analysis demonstrates that the development of pathology is probabilistic and can be described as a time -to-event process . The proposed approach enables quantitative assessment of the individual risk of developing mitochondrial diseases and opens the door to personalized prognosis.

Purpose: Elective genomic sequencing (EGS) returns monogenic disease findings in multiple genes, including potentially novel variants, and may also provide participants with carrier status, pharmacogenomic and other health-related information. The PeopleSeq Study assessed participants' motivations for and concerns about EGS and the associated clinical and psychosocial outcomes across diverse EGS providers. Methods: We administered a shared questionnaire to participants who chose to undergo EGS via 18 academic, clinical, or commercial EGS platforms. Results: We enrolled 1575 participants, of whom 1147 (72.8%) completed a questionnaire after receiving their EGS results. A majority (60.3%) of the participants who completed a post-result questionnaire self-reported receiving results they assessed as important, including negative findings, and 75.9% reported a form of health-related utility. Among a subset (19.4%) who shared their EGS reports, 16.6% (37 of n=223) received a monogenic finding and self-reported results deemed "important" were consistent with EGS reports. Most participants (74.1%) discussed their results with their family, but fewer discussed their results with a healthcare provider other than the site team (41.7%) or had one or more medical visits as a direct result of their EGS testing (23.1%). Participants expressed diverse motivations for EGS, with 91.4% expressing interest in their personal disease risk and 54% who expressed quasi-indication-based motivations related to family medical history. Individuals motivated by family history reported important results at a significantly higher rate. Conclusions: Early adopters of EGS are motivated by general interest in their health as well as quasi-indication-based considerations such as family history. A majority of participants learned results they considered medically important, but a much smaller segment engaged healthcare providers with their results.