Genotype-phenotype Relationships Research Articles

Mechanisms and pathologies of human mitochondrial DNA replication and deletion formation.

Human mitochondria possess a multi-copy circular genome, mitochondrial DNA (mtDNA), that is essential for cellular energy metabolism. The number of copies of mtDNA per cell, and their integrity, are maintained by nuclear-encoded mtDNA replication and repair machineries. Aberrant mtDNA replication and mtDNA breakage are believed to cause deletions within mtDNA. The genomic location and breakpoint sequences of these deletions show similar patterns across various inherited and acquired diseases, and are also observed during normal ageing, suggesting a common mechanism of deletion formation. However, an ongoing debate over the mechanism by which mtDNA replicates has made it difficult to develop clear and testable models for how mtDNA rearrangements arise and propagate at a molecular and cellular level. These deletions may impair energy metabolism if present in a high proportion of the mtDNA copies within the cell, and can be seen in primary mitochondrial diseases, either in sporadic cases or caused by autosomal variants in nuclear-encoded mtDNA maintenance genes. These mitochondrial diseases have diverse genetic causes and multiple modes of inheritance, and show notoriously broad clinical heterogeneity with complex tissue specificities, which further makes establishing genotype-phenotype relationships challenging. In this review, we aim to cover our current understanding of how the human mitochondrial genome is replicated, the mechanisms by which mtDNA replication and repair can lead to mtDNA instability in the form of large-scale rearrangements, how rearranged mtDNAs subsequently accumulate within cells, and the pathological consequences when this occurs.

The distribution of fitness effects during adaptive walks using a simple genetic network.

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.

Distinct Network Patterns Emerge from Cartesian and XOR Epistasis Models: A Comparative Network Science Analysis.

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, can significantly contribute to the observed phenotypic variance of complex traits. To date, it has been generally assumed that genetic interactions can be detected using a Cartesian, or multiplicative, interaction model commonly utilized in standard regression approaches. However, a recent study investigating epistasis in obesity-related traits in rats and mice has identified potential limitations of the Cartesian model, revealing that it only detects some of the genetic interactions occurring in these systems. By applying an alternative approach, the exclusive-or (XOR) model, the researchers detected a greater number of epistatic interactions and identified more biologically relevant ontological terms associated with the interacting loci. This suggests that the XOR model may provide a more comprehensive understanding of epistasis in these species and phenotypes. To further explore these findings and determine if different interaction models also make up distinct epistatic networks, we leverage network science to provide a more comprehensive view into the genetic interactions underlying BMI in this system. Our comparative analysis of networks derived from Cartesian and XOR interaction models in rats (Rattus norvegicus) uncovers distinct topological characteristics for each model-derived network. Notably, we discover that networks based on the XOR model exhibit an enhanced sensitivity to epistatic interactions. This sensitivity enables the identification of network communities, revealing novel trait-related biological functions through enrichment analysis. Furthermore, we identify triangle network motifs in the XOR epistatic network, suggestive of higher-order epistasis, based on the topology of lower-order epistasis. These findings highlight the XOR model's ability to uncover meaningful biological associations as well as higher-order epistasis from lower-order epistatic networks. Additionally, our results demonstrate that network approaches not only enhance epistasis detection capabilities but also provide more nuanced understandings of genetic architectures underlying complex traits. The identification of community structures and motifs within these distinct networks, especially in XOR, points to the potential for network science to aid in the discovery of novel genetic pathways and regulatory networks. Such insights are important for advancing our understanding of phenotype-genotype relationships.

Massively parallel CRISPR-assisted homologous recombination enables saturation editing of full-length endogenous genes in yeast.

Performing saturation editing of chromosomal genes will enable the study of genetic variants in situ and facilitate protein and cell engineering. However, current in vivo editing of endogenous genes either lacks flexibility or is limited to discrete codons and short gene fragments, preventing a comprehensive exploration of genotype-phenotype relationships. To enable facile saturation editing of full-length genes, we used a protospacer adjacent motif-relaxed Cas9 variant and homology-directed repair to achieve above 60% user-defined codon replacement efficiencies in Saccharomyces cerevisiae genome. Coupled with massively parallel DNA design and synthesis, we developed a saturation gene editing method termed CRISPR-Cas9- and homology-directed repair-assisted saturation editing (CHASE) and achieved highly saturated codon swapping of long genomic regions. By applying CHASE to massively edit a well-studied global transcription factor gene, we found known and unreported genetic variants affecting an industrially relevant microbial trait. The user-defined codon editing capability and wide targeting windows of CHASE substantially expand the scope of saturation gene editing.

Abstract A010: Cancer associated fibroblast dependency in a novel immunocompetent mouse model of bladder cancer

Abstract Introduction: The relationship between genotype and phenotype in targeted therapy has proven to be crucial for the successful use of most kinase inhibitors. In bladder cancer, the PI3K/AKT/mTOR pathway is widely activated by mutation but untargeted. Our study focuses on developing an immunocompetent, CRISPR-edited bladder cancer mouse model that mirrors common PI3K/AKT/mTOR pathway mutations. We established effective genotype-kinase inhibitor pairings using CRISPR-edited organoids and cancer associated fibroblasts (CAFs) to create a pro-tumorigenic environment in vivo. Methods and Results: Urothelial cells were isolated from the bladder of Rosa26 Lox-Stop-Lox-Cas9/GFP transgenic mice and grown as organoids in matrigel. These organoids were co-edited for Trp53 and Rb1 KO, and further editing of the Pten gene resulted in a nine-fold increase in sensitivity to the kinase inhibitor ipatasertib. Despite successful initial engraftment into mouse bladder, these edited cells regressed before tumor formation, and subcutaneous injection did not yield subcutaneous tumors. An alternate model featuring a 9p21 deletion with Trp53 KO also displayed similar regression post-engraftment into bladder or subcutaneous injections. However, co-injection with a previously established CAF line subcutaneously yielded palpable subcutaneous tumors that enlarged for 12 weeks. We then established a BBN bladder-tumor-derived CAF line, leading to the formation of tumors with squamous differentiation upon intramural co-injection with edited organoids. These models displayed immune-excluded architecture and a predominance of M2 over M1 macrophages. We also aimed to enhance the model's adaptability by introducing common PI3K/AKT/mTOR mutations, specifically TSC1/TSC2 KO, and Pik3ca hotspots, into our 9p21/Trp53 KO model. While TSC1/TSC2 KO lines exhibited elevated P-S6K signaling, they displayed resistance to Torin 1, Voxtalisib, and Everolimus—kinase inhibitors targeting mTORC1/2. In parallel, we employed CRISPR-mediated homologous recombination to successfully introduce a Pik3ca H1047R hotspot into the 9p21/Trp53 KO genotype, as evidenced by increased Phosphorylation of AKT compared to wild-type 9p21/Trp53 KO cell lines. This genotype resulted in approximately 80-fold selective increased sensitivity to pictilisib, a PI3K inhibitor compared to Pik3ca WT control. Conclusions: In this model, CAFs were required for successful engraftment of tumors into the bladder wall or subcutaneously. BUTU is a customizable model of bladder cancer and is suitable for direct targeting by kinase inhibitors while providing a platform for understanding the genotype-phenotype relationship in bladder cancer. We hope to use it as a tool to develop preclinical evidence on which to base precision medicine clinical trials. Citation Format: Philip H. Abbosh, Henkel Valentine. Cancer associated fibroblast dependency in a novel immunocompetent mouse model of bladder cancer [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr A010.

Multimodal Learning for Mapping the Genotype-Phenotype Dynamics.

How complex phenotypes emerge from intricate gene expression patterns is a fundamental question in biology. Quantitative characterization of this relationship, however, is challenging due to the vast combinatorial possibilities and dynamic interplay between genotype and phenotype landscapes. Integrating high-content genotyping approaches such as single-cell RNA sequencing and advanced learning methods such as language models offers an opportunity for dissecting this complex relationship. Here, we present a computational integrated genetics framework designed to analyze and interpret the high-dimensional landscape of genotypes and their associated phenotypes simultaneously. We applied this approach to develop a multimodal foundation model to explore the genotype-phenotype relationship manifold for human transcriptomics at the cellular level. Analyzing this joint manifold showed a refined resolution of cellular heterogeneity, enhanced precision in phenotype annotating, and uncovered potential cross-tissue biomarkers that are undetectable through conventional gene expression analysis alone. Moreover, our results revealed that the gene networks are characterized by scale-free patterns and show context-dependent gene-gene interactions, both of which result in significant variations in the topology of the gene network, particularly evident during aging. Finally, utilizing contextualized embeddings, we investigated gene polyfunctionality which illustrates the multifaceted roles that genes play in different biological processes, and demonstrated that for VWF gene in endothelial cells. Overall, this study advances our understanding of the dynamic interplay between gene expression and phenotypic manifestation and demonstrates the potential of integrated genetics in uncovering new dimensions of cellular function and complexity.

Congenital Long QT Syndrome in Children and Adolescents: A General Overview.

Congenital long QT syndrome (LQTS) represents a disorder of myocardial repolarization characterized by a prolongation of QTc interval on ECG, which can degenerate into fast polymorphic ventricular arrhythmias. The typical symptoms of LQTS are syncope and palpitations, mainly triggered by adrenergic stimuli, but it can also manifest with cardiac arrest. At least 17 genotypes have been associated with LQTS, with a specific genotype-phenotype relationship described for the three most common subtypes (LQTS1, -2, and -3). β-Blockers are the first-line therapy for LQTS, even if the choice of the appropriate patients needing to be treated may be challenging. In specific cases, interventional measures, such as an implantable cardioverter-defibrillator (ICD) or left cardiac sympathetic denervation (LCSD), are useful. The aim of this review is to highlight the current state-of-the-art knowledge on LQTS, providing an updated picture of possible diagnostic algorithms and therapeutic management.

Clinical, genetic, and cognitive correlates of seizure occurrences in Phelan-McDermid syndrome

BackgroundPhelan-McDermid syndrome (PMS) is a genetic neurodevelopmental disorder caused by SHANK3 haploinsufficiency and is associated with an increased risk for seizures. Previous literature indicates that around one third of individuals with PMS also have epilepsy or seizures, with a wide range of types and ages of onset. Investigating the impact of seizures on intellectual and adaptive functioning for PMS is a primary concern for caregivers and is important to understanding the natural history of this syndrome.MethodsWe report on results from 98 individuals enrolled in a prospective, longitudinal study. We detailed seizure frequency, type, and age of onset, and we analyzed seizure occurrence with best estimate IQ, adaptive functioning, clinical features, and genotype. We conducted multiple linear regression analyses to assess the relationship between the presence of seizures and the Vineland Adaptive Behavior Scale, Second Edition (VABS-II) Adaptive Behavior Composite score and the best estimate full-scale IQ. We also performed Chi-square tests to explore associations between seizure prevalence and genetic groupings. Finally, we performed Chi-square tests and t-tests to explore the relationship between seizures and demographic features, features that manifest in infancy, and medical features.ResultsSeizures were present in 41% of the cohort, and age of onset was widely variable. The presence of seizures was associated with significantly lower adaptive and intellectual functioning. Genotype–phenotype analyses were discrepant, with no differences in seizure prevalence across genetic classes, but with more genes included in deletions of participants with 22q13 deletions and seizures compared to those with 22q13 deletions and no seizures. No clinical associations were found between the presence of seizures and sex, history of pre- or neonatal complications, early infancy, or medical features. In this cohort, generalized seizures were associated with developmental regression, which is a top concern for PMS caregivers.ConclusionsThese results begin to eludicate correlates of seizures in individuals with PMS and highlight the importance of early seizure management. Importantly, presence of seizures was associated with adaptive and cognitive functioning. A larger cohort might be able to identify additional associations with medical features. Genetic findings suggest an increased capability to realize genotype–phenotype relationships when deletion size is taken into account.

Identification of bacterial determinants of tuberculosis infection and treatment outcomes: a phenogenomic analysis of clinical strains

Bacterial diversity could contribute to the diversity of tuberculosis infection and treatment outcomes observed clinically, but the biological basis of this association is poorly understood. The aim of this study was to identify associations between phenogenomic variation in Mycobacterium tuberculosis and tuberculosis clinical features. We developed a high-throughput platform to define phenotype-genotype relationships in M tuberculosis clinical isolates, which we tested on a set of 158 drug-sensitive M tuberculosis strains sampled from a large tuberculosis clinical study in Ho Chi Minh City, Viet Nam. We tagged the strains with unique genetic barcodes in multiplicate, allowing us to pool the strains for in-vitro competitive fitness assays across 16 host-relevant antibiotic and metabolic conditions. Relative fitness was quantified by deep sequencing, enumerating output barcode read counts relative to input normalised values. We performed a genome-wide association study to identify phylogenetically linked and monogenic mutations associated with the in-vitro fitness phenotypes. These genetic determinants were further associated with relevant clinical outcomes (cavitary disease and treatment failure) by calculating odds ratios (ORs) with binomial logistic regressions. We also assessed the population-level transmission of strains associated with cavitary disease and treatment failure using terminal branch length analysis of the phylogenetic data. M tuberculosis clinical strains had diverse growth characteristics in host-like metabolic and drug conditions. These fitness phenotypes were highly heritable, and we identified monogenic and phylogenetically linked variants associated with the fitness phenotypes. These data enabled us to define two genetic features that were associated with clinical outcomes. First, mutations in Rv1339, a phosphodiesterase, which were associated with slow growth in glycerol, were further associated with treatment failure (OR 5·34, 95% CI 1·21-23·58, p=0·027). Second, we identified a phenotypically distinct slow-growing subclade of lineage 1 strains (L1.1.1.1) that was associated with cavitary disease (OR 2·49, 1·11-5·59, p=0·027) and treatment failure (OR 4·76, 1·53-14·78, p=0·0069), and which had shorter terminal branch lengths on the phylogenetic tree, suggesting increased transmission. Slow growth under various antibiotic and metabolic conditions served as in-vitro intermediate phenotypes underlying the association between M tuberculosis monogenic and phylogenetically linked mutations and outcomes such as cavitary disease, treatment failure, and transmission potential. These data suggest that Mtuberculosis growth regulation is an adaptive advantage for bacterial success in human populations, at least in some circumstances. These data further suggest markers for the underlying bacterial processes that contribute to these clinical outcomes. National Health and Medical Research Council/A∗STAR, National Institutes of Allergy and Infectious Diseases, National Institute of Child Health and Human Development, and the Wellcome Trust Fellowship in Public Health and Tropical Medicine.

Genotype-phenotype association and functional analysis of hnRNPA1 mutations in amyotrophic lateral sclerosis

Background Pathogenic variants in hnRNPA1 have been reported in amyotrophic lateral sclerosis (ALS) patients. However, studies on hnRNPA1 mutant spectrum and pathogenicity of variants were rare. Methods We performed whole exome sequencing of ALS-associated genes and subsequent verification of rare variants in hnRNPA1 in our ALS patients. The hnRNPA1 mutations reported in literature were reviewed and combined with our results to determine the genotype-phenotype relationship. Functional analysis of the novel variant p.G195A was performed in vitro by transfection of mutant hnRNPA1 into 293T cell. Results Among 207 ALS patients recruited, 3 rare hnRNPA1 variants were identified (mutant frequency 1.45%), including two recurrent mutations (p.P340S and p.G283R), and a novel rare variant p.G195A. In combination with previous reports, there are 27 ALS patients with 15 hnRNPA1 mutations identified. Disease onset age was 47.90 ± 1.52 years with predominant limb onset. The p.P340S mutation caused flail arm syndrome (FAS) in two independent families with extended life expectancy. The newly identified p.G195A mutation, lying at the start of the PrLD (“prion-like” domain)/LCD (low-complexity domain), causes local structural changes in 3D protein prediction. Upon sodium arsenite exposure, mutant hnRNPA1 retained in the nucleus but deficit of cytoplasmic G3BP1-positive stress granule clearance was observed. This is different from the p.P340S mutation which caused both cytoplasmic translocation and stress granule formation. No cytoplasmic TDP-43 translocation was observed. Conclusion Mutations in hnRNPA1 are overall minor in ALS patients. The p.P340S mutation is associated with manifestation of FAS. Mutations in LCD of hnRNPA1 cause stress granule misprocessing.

Genotype and clinical phenotype of children with Marfan syndrome in Southeastern Anatolia

The cardinal phenotypic hallmarks of Marfan syndrome (MFS) include cardiac, ocular, and skeletal abnormalities. Since the clinical phenotype of MFS is highly heterogeneous, with certain symptoms appearing as children age, the diagnostic process and establishing a genotype-phenotype association in childhood MFS can be challenging. The lack of sufficient childhood studies also makes it difficult to interpret the subject. This study aims to evaluate the relationship between clinical symptoms used as diagnostic criteria and FBN1 variations in children with MFS. This study investigated the relationships between genotypes and phenotypes in 131 children suspected of having Marfan syndrome (MFS). Diagnosis of MFS was made according to the revised Ghent nosology. FBN1 variants were categorized based on exon regions, type of variant, and pathogenicity classes. These FBN1 variants were then correlated with the clinical manifestations including cardiovascular, ocular, facial, and skeletal abnormalities. Out of the children, 43 were diagnosed with MFS. FBN1 variant was identified in 32 (74.4%) of the MFS children. MFS diagnosis could not be made in five (15.6%) FBN1 variant-positive children. The most common cardinal finding is cardiac anomalies n = 38 (88.3%). The most common FBN1 pathogenic variant was c.1786 T > C/p.Cys596Arg n = 4 (12.5%). The distribution of pathogenic variants was as follows: 29 (90.6%) missense, 2 (6.3%) frameshift, and 1 (3.1%) nonsense. The numbers of AD and EL of the variant-positive children were 16 (50%) and 14 (43.7%), respectively. Ocular abnormalities were more common in children with FBN1-positive MFS (p = 0.009). There was no difference in the number of cardiac abnormalities between FBN1-positive and FBN1-negative MFS patients (p = 0.139). Conclusion: This study examines the relationship between FBN1 variants and clinical features used as diagnostic criteria in MFS children. The findings emphasize the importance of long-term monitoring of heterogeneous clinical phenotypes and bioinformatic reanalysis in determining the genotype-phenotype relationship in children, as MFS symptoms can vary with age.What is Known:• Marfan syndrome has highly variable phenotypic heterogeneity.• The genotype-phenotype relationship in childhood Marfan syndrome is not clear enough due to the variation in the time of onset of the findings.What is New:• This article provides regional data for the field of research on genotype-phenotype relationships in childhood Marfan syndrome.• Long-term follow-up of clinical findings and bioinformatics reanalysis is an important requirement for a well-established genotype-phenotype relationship in childhood Marfan syndrome.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Interplay between α-thalassemia and β-hemoglobinopathies: Translating genotype-phenotype relationships into therapies.

α-Thalassemia represents one of the most important genetic modulators of β-hemoglobinopathies. During this last decade, the ongoing interest in characterizing genotype-phenotype relationships has yielded incredible insights into α-globin gene regulationand its impact on β-hemoglobinopathies. In this review, we provide a holistic update on α-globin gene expression stemming from DNA to RNA to protein, as well as epigenetic mechanisms that can impact gene expression and potentially influence phenotypic outcomes. Here, we highlight defined α-globin targeted strategies and rationalize the use of distinct molecular targets based on the restoration of balanced α/β-like globin chain synthesis. Considering the therapies that either increase β-globin synthesis or reactivate γ-globin gene expression, the modulation of α-globin chains as a disease modifier for β-hemoglobinopathies still remains largely uncharted in clinical studies.

Expanding the clinical phenotype and genetic spectrum of GEMIN5 disorders: Early-infantile developmental and epileptic encephalopathies.

Several biallelic truncating and missense variants of the gem nuclear organelle-associated protein 5 (GEMIN5) gene have been reported to cause neurodevelopmental disorders characterized by cerebellar atrophy, intellectual disability, and motor dysfunction. However, the association between biallelic GEMIN5 variants and early-infantile developmental and epileptic encephalopathies (EIDEEs) has not been reported. This study aimed to expand the phenotypic spectrum of GEMIN5 and explore the correlations between epilepsy and molecular sub-regional locations. We performed whole-exome sequencing in two patients with EIDEE with unexplained etiologies. The damaging effects of variants were predicted using multiple in silico tools and modeling. All reported patients with GEMIN5 pathogenic variants and detailed neurological phenotypes were analyzed to evaluate the genotype-phenotype relationship. Novel biallelic GEMIN5 variants were identified in two unrelated female patients with EIDEE, including a frameshift variant (Hg19, chr5:154284147-154284148delCT: NM_015465: c.2551_c.2552delCT: p.(Leu851fs*30)), a nonsense mutation (Hg19, chr5:154299603-154299603delTinsAGA: NM_015465: c.1523delTinsAGA: p.(Leu508*)), and two missense variants (Hg19, chr5:154282663T>A: NM_015465: c.2705T>A: p.(Leu902Gln) and Hg19, chr5:154281002C>G: NM_015465: c.2911C>G: p.(Gln971Glu)), which were inherited from asymptomatic parents and predicted to be damaging or probably damaging using in silico tools. Except p.Leu508*, all these mutations are located in tetratricopeptide repeat (TPR) domain. Our two female patients presented with seizures less than 1 month after birth, followed by clusters of spasms. Brain magnetic resonance imaging suggests dysgenesis of the corpus callosum and cerebellar hypoplasia. Video electroencephalogram showed suppression-bursts. Through a literature review, we found 5 published papers reporting 48 patients with biallelic variants in GEMIN5. Eight of 48 patients have epilepsy, and 5 patients started before 1 year old, which reminds us of the relevance between GEMIN5 variants and EIDEE. Further analysis of the 49 GEMIN5 variants in those 50 patients demonstrated that variants in TPR-like domain or RBS domain were more likely to be associated with epilepsy. We found novel biallelic variants of GEMIN5 in two individuals with EIDEE and expanded the clinical phenotypes of GEMIN5 variants. It is suggested that the GEMIN5 gene should be added to the EIDEE gene panel to aid in the clinical diagnosis of EIDEE and to help determine patient prognosis.

Artificial selection footprints in indigenous and commercial chicken genomes

BackgroundAlthough many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China.ResultsWe found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds.ConclusionWe identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

Incomplete-penetrant hypertrophic cardiomyopathy MYH7 G256E mutation causes hypercontractility and elevated mitochondrial respiration

Determining the pathogenicity of hypertrophic cardiomyopathy-associated mutations in the β-myosin heavy chain (MYH7) can be challenging due to its variable penetrance and clinical severity. This study investigates the early pathogenic effects of the incomplete-penetrant MYH7 G256E mutation on myosin function that may trigger pathogenic adaptations and hypertrophy. We hypothesized that the G256E mutation would alter myosin biomechanical function, leading to changes in cellular functions. We developed a collaborative pipeline to characterize myosin function across protein, myofibril, cell, and tissue levels to determine the multiscale effects on structure-function of the contractile apparatus and its implications for gene regulation and metabolic state. The G256E mutation disrupts the transducer region of the S1 head and reduces the fraction of myosin in the folded-back state by 33%, resulting in more myosin heads available for contraction. Myofibrils from gene-edited MYH7WT/G256E human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) exhibited greater and faster tension development. This hypercontractile phenotype persisted in single-cell hiPSC-CMs and engineered heart tissues. We demonstrated consistent hypercontractile myosin function as a primary consequence of the MYH7 G256E mutation across scales, highlighting the pathogenicity of this gene variant. Single-cell transcriptomic and metabolic profiling demonstrated upregulated mitochondrial genes and increased mitochondrial respiration, indicating early bioenergetic alterations. This work highlights the benefit of our multiscale platform to systematically evaluate the pathogenicity of gene variants at the protein and contractile organelle level and their early consequences on cellular and tissue function. We believe this platform can help elucidate the genotype-phenotype relationships underlying other genetic cardiovascular diseases.

Syndromic craniosynostosis caused by a novel missense variant in MAP4K4: Expanding the genotype-phenotype relationship in RASopathies.

RASopathies represent a distinct class of neurodevelopmental syndromes caused by germline variants in the Ras/MAPK pathways. Recently, a novel disease-gene association was implicated in MAPK kinase kinase kinase 4 (MAP4K4), which regulates the upstream signals of the MAPK pathways. However, to our knowledge, only two studies have reported the genotype-phenotype relationships in the MAP4K4-related disorder. This study reports on a Korean boy harboring a novel de novo missense variant in MAP4K4 (NM_001242559:c.569G>T, p.Gly190Val), revealed by trio exome sequencing, and located in the hotspot of the protein kinase domain. The patient exhibited various clinical features, including craniofacial dysmorphism, language delay, congenital heart defects, genitourinary anomalies, and sagittal craniosynostosis. Our study expands the phenotypic association of the MAP4K4-related disorder to include syndromic craniosynostosis, thereby providing further insights into the role of the RAS/MAPK pathways in the development of premature fusion of calvarial sutures.

Clinical and Laboratory Characteristics of MODY Cases, Genetic Mutation Spectrum and Phenotype-genotype Relationship

Maturity onset diabetes of the young (MODY) occurs due to mutations in genes involved in pancreatic beta cell function and insulin secretion, has heterogeneous clinical and laboratory features, and account for 1-5% of all diabetes cases. The prevalence and distribution of MODY subtypes vary between countries. The aim of this study was to evaluate the clinical and laboratory characteristics, mutation distribution, and phenotype-genotype relationship in a large case series of pediatric Turkish patients genetically diagnosed with MODY. MODY cases from 14 different pediatric endocrinology departments were included. Diagnosis, treatment, follow-up data, and results of genetic analysis were evaluated. A total of 224 patients were included, of whom 101 (45%) were female, and the mean age at diagnosis was 9.4±4.1 years. Gene variant distribution was: 146 (65%) GCK; 43 (19%) HNF1A; 8 (3.6%) HNF4A, 8 (3.6%) KLF11 and 7 (3.1%) HNF1B. The remaining 12 variants were: PDX (n=1), NEUROD1 (n=3), CEL (n=1), INS (n=3), ABCC8 (n= 3) and KJNC11 (n=1). Of the cases, 197 (87.9%) were diagnosed with incidental hyperglycemia, 16 with ketosis (7%) and 7 (3%) with diabetic ketoacidosis (DKA), while 30% presented with classical symptoms of diabetes. Two-hundred (89%) had a family history of diabetes. Anti-GAD antibody was detected in 13 cases, anti-islet antibody in eight and anti-insulin antibody in four. Obesity was present in 16. Distribution of therapy was: 158 (71%) diet only; 23 (11%) intensive insulin treatment; 17 (7.6%) sulfonylureas; 10 (4.5%) metformin; and 6 (2.7%) insulin and oral anti-diabetic treatment. This was the largest genetically diagnosed series from Turkey. The most common gene variants were GCK and HNF1A with much lower proportions for other MODY types. Hyperglycemia was the most common presenting symptom while 11% of patients had diabetes-associated autoantibodies and 7% were obese. The majority of patients received dietary management only.

Early embryonic failure caused by a novel mutation in the TUBB8 gene: A case report

BACKGROUND This study aimed to explore the relationship between gene mutations and early embryonic development arrest and to provide more possibilities for the diagnosis and treatment of repeated implantation failure. CASE SUMMARY Here, we collected and described the clinical data of a patient with early embryonic development stagnation after repeated in vitro fertilization attempts for primary infertility at the Department Reproductive Center of Zaozhuang Maternal and Child Healthcare Hospital. We also detected the whole-exon gene of the patient's spouse and parents, and conducted bioinformatics analysis to determine the pathogenesis of the gene. CONCLUSION A novel mutant of the TUBB8 gene [c.602G&gt;T(p.C201F)] was identified, and this mutant provided new data on the genotype-phenotype relationships of related diseases.

Variable Penetrance and Expressivity of a Rare Pore Loss-of-Function Mutation (p.L889V) of Nav1.5 Channels in Three Spanish Families.

A novel rare mutation in the pore region of Nav1.5 channels (p.L889V) has been found in three unrelated Spanish families that produces quite diverse phenotypic manifestations (Brugada syndrome, conduction disease, dilated cardiomyopathy, sinus node dysfunction, etc.) with variable penetrance among families. We clinically characterized the carriers and recorded the Na+ current (INa) generated by p.L889V and native (WT) Nav1.5 channels, alone or in combination, to obtain further insight into the genotypic-phenotypic relationships in patients carrying SCN5A mutations and in the molecular determinants of the Nav1.5 channel function. The variant produced a strong dominant negative effect (DNE) since the peak INa generated by p.L889V channels expressed in Chinese hamster ovary cells, either alone (-69.4 ± 9.0 pA/pF) or in combination with WT (-62.2 ± 14.6 pA/pF), was significantly (n ≥ 17, p < 0.05) reduced compared to that generated by WT channels alone (-199.1 ± 44.1 pA/pF). The mutation shifted the voltage dependence of channel activation and inactivation to depolarized potentials, did not modify the density of the late component of INa, slightly decreased the peak window current, accelerated the recovery from fast and slow inactivation, and slowed the induction kinetics of slow inactivation, decreasing the fraction of channels entering this inactivated state. The membrane expression of p.L889V channels was low, and in silico molecular experiments demonstrated profound alterations in the disposition of the pore region of the mutated channels. Despite the mutation producing a marked DNE and reduction in the INa and being located in a critical domain of the channel, its penetrance and expressivity are quite variable among the carriers. Our results reinforce the argument that the incomplete penetrance and phenotypic variability of SCN5A loss-of-function mutations are the result of a combination of multiple factors, making it difficult to predict their expressivity in the carriers despite the combination of clinical, genetic, and functional studies.