Structural basis of adenine nucleotides regulation and neurodegenerative pathology in ClC-3 exchanger.

The ClC-3 chloride/proton exchanger is both physiologically and pathologically critical, as it is potentiated by ATP to detect metabolic energy level and point mutations in ClC-3 lead to severe neurodegenerative diseases in human. However, why this exchanger is differentially modulated by ATP, ADP or AMP and how mutations caused gain-of-function remains largely unknow. Here we determine the high-resolution structures of dimeric wildtype ClC-3 in the apo state and in complex with ATP, ADP and AMP, and the disease-causing I607T mutant in the apo and ATP-bounded state by cryo-electron microscopy. In combination with patch-clamp recordings and molecular dynamic simulations, we reveal how the adenine nucleotides binds to ClC-3 and changes in ion occupancy between apo and ATP-bounded state. We further observe I607T mutation induced conformational changes and augments in current. Therefore, our study not only lays the structural basis of adenine nucleotides regulation in ClC-3, but also clearly indicates the target region for drug discovery against ClC-3 mediated neurodegenerative diseases.

Adaptive gene loss in the common bean pan-genome during range expansion and domestication.

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Allosteric inhibition of CFTR gating by CFTRinh-172 binding in the pore.

Loss-of-function mutations of the CFTR gene cause the life-shortening genetic disease cystic fibrosis (CF), whereas overactivity of CFTR may lead to secretory diarrhea and polycystic kidney disease. While effective drugs targeting the CFTR protein have been developed for the treatment of CF, little progress has been made for diseases caused by hyper-activated CFTR. Here, we solve the cryo-EM structure of CFTR in complex with CFTRinh-172 (Inh-172), a CFTR gating inhibitor with promising potency and efficacy. We find that Inh-172 binds inside the pore of CFTR, interacting with amino acid residues from transmembrane segments (TMs) 1, 6, 8, 9, and 12 through mostly hydrophobic interactions and a salt bridge. Substitution of these residues lowers the apparent affinity of Inh-172. The inhibitor-bound structure reveals re-orientations of the extracellular segment of TMs 1, 8, and 12, supporting an allosteric modulation mechanism involving post-binding conformational changes. This allosteric inhibitory mechanism readily explains our observations that pig CFTR, which preserves all the amino acid residues involved in Inh-172 binding, exhibits a much-reduced sensitivity to Inh-172 and that the apparent affinity of Inh-172 is altered by the CF drug ivacaftor (i.e., VX-770) which enhances CFTR's activity through binding to a site also comprising TM8.

Variation in a Poaceae-conserved fatty acid metabolic gene cluster controls rice yield by regulating male fertility.

A wide variety of metabolic gene clusters exist in eukaryotic genomes, but fatty acid metabolic gene clusters have not been discovered. Here, combining with metabolic and phenotypic genome-wide association studies, we identify a major locus containing a six-gene fatty acid metabolic gene cluster on chromosome 3 (FGC3) that controls the cutin monomer hydroxymonoacylglycerols (HMGs) contents and rice yield, possibly through variation in the transcription of FGC3 members. We show that HMGs are sequentially synthesized in the endoplasmic reticulum by OsFAR2, OsKCS11, OsGPAT6, OsCYP704B2 and subsequently transported to the apoplast by OsABCG22 and OsLTPL82. Mutation of FGC3 members reduces HMGs, leading to defective male reproductive development and a significant decrease in yield. OsMADS6 and OsMADS17 directly regulate FGC3 and thus influence male reproduction and yield. FGC3 is conserved in Poaceae and likely formed prior to the divergence of Pharus latifolius. The eukaryotic fatty acid and plant primary metabolic gene cluster we identified show a significant impact on the origin and evolution of Poaceae and has potential for application in hybrid crop breeding.

Varied clinical presentations of RP1L1 variants in Chinese patients: a study of occult macular dystrophy and vitelliform macular dystrophy.

BACKGROUND: Occult Macular Dystrophy (OMD), primarily caused by retinitis pigmentosa 1-like 1 (RP1L1) variants, is a complex retinal disease characterised by progressive vision loss and a normal fundus appearance. This study aims to investigate the diverse phenotypic expressions and genotypic correlations of OMD in Chinese patients, including a rare case of Vitelliform Macular Dystrophy (VMD) associated with RP1L1. METHODS: We analysed seven OMD patients and one VMD patient, all with heterozygous pathogenic RP1L1 variants. Clinical assessments included Best Corrected Visual Acuity (BCVA), visual field testing, Spectral Domain Optical Coherence Tomography (SD-OCT), multifocal Electroretinograms (mfERGs), and microperimetry. Next-generation sequencing was utilised for genetic analysis. RESULTS: The OMD patients displayed a range of phenotypic variability. Most (5 out of 7) had the RP1L1 variant c.133 C > T; p.R45W, associated with central vision loss and specific patterns in SD-OCT and mfERG. Two patients exhibited different RP1L1 variants (c.3599G > T; p.G1200V and c.2880G > C; p.W960C), presenting milder phenotypes. SD-OCT revealed photoreceptor layer changes, with most patients showing decreased mfERG responses in the central rings. Interestingly, a unique case of VMD linked to the RP1L1 variant was observed, distinct from traditional OMD presentations. CONCLUSIONS: This study highlights the phenotypic diversity within OMD and the broader spectrum of RP1L1-associated macular dystrophies, including a novel association with VMD. The findings emphasise the complexity of RP1L1 variants in determining clinical manifestations, underscoring the need for comprehensive genetic and clinical evaluations in macular dystrophies.

Development of prognostic model incorporating a ferroptosis/cuproptosis-related signature and mutational landscape analysis in muscle-invasive bladder cancer.

BACKGROUND: Muscle-invasive bladder cancer (MIBC) is a prevalent and aggressive malignancy. Ferroptosis and cuproptosis are recently discovered forms of programmed cell death (PCD) that have attracted much attention. However, their interactions and impacts on MIBC overall survival (OS) and treatment outcomes remain unclear. METHODS: Data from the TCGA-BLCA project (as the training set), cBioPortal database, and GEO datasets (GSE13507 and GSE32894, as the test sets) were utilized to identify hub ferroptosis/cuproptosis-related genes (FRGs and CRGs) and develop a prognostic signature. Differential expression analysis (DEA) was conducted, followed by univariate and multivariate Cox's regression analyses and multiple machine learning (ML) techniques to select genetic features. The performance of the ferroptosis/cuproptosis-related signature was evaluated using Kaplan-Meier (K-M) survival analysis and receiver-operating characteristics (ROC) curves. Mutational and tumour immune microenvironment landscapes were also explored. Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) experiments confirmed the expression patterns of the hub genes, and functional assays assessed the effects of SCD knockdown on cell viability, proliferation, and migration. RESULTS: DEA revealed dysregulated FRGs and CRGs in the TCGA MIBC cohort. SCD, DDR2, and MT1A were identified as hub genes. A prognostic signature based on the sum of the weighted expression of these genes demonstrated strong predictive efficacy in the training and test sets. Nomogram incorporating this signature accurately predicted 1-, 3-, and 5-year survival probabilities in the TCGA cohort and GSE13507 dataset. Copy number variation (CNV) and tumour immune microenvironment analysis revealed that high risk score level groups were associated with immunosuppression and lower tumour purity. The associations of risk scores with immunotherapy and chemical drugs were also explored, indicating their potential for guiding treatment for MIBC patients. The dysregulated expression patterns of three hub genes were validated by RT-qPCR experiments. CONCLUSIONS: Targeting hub FRGs and CRGs could be a promising therapeutic approach for MIBC. Our prognostic model offers a new framework for MIBC subtyping and can inform personalized therapeutic strategies.

Heterozygous variants in transmembrane channel-like 1 gene cause autosomal recessive nonsyndromic hearing loss.

Autosomal recessive non-syndromic hearing loss (ARNSHL) can cause severe or very severe pre-speech hearing loss. Transmembrane channel-like 1 (TMC1) gene is the sixth deafness gene discovered, but the precise extent of its protein structure and function is unknown. First, history collection, audiology examination and imaging examination were performed on the proband and his family members. Peripheral blood of proband and family members was collected, genomic DNA was extracted, exon high-throughput sequencing technology was used to detect the deafness gene mutation of the proband, and Sanger sequencing was performed to verify the TMC1 gene of the proband's parents. The proband was born with hearing impairment, normal tympanic function, inability to induce acoustic reflex in both ears (acoustic reflex threshold is 100 dBHL), and severe sensorineural deafness. One of his sisters has severe sensorineural hearing loss, and neither his parents nor his other sister is hearing impaired. High-throughput sequencing of the proband identified mutations at c.741+3_741+6delAAGT (splicing) and c.884C>T (p.A295V) of the TMC1 gene, two of which were heterozygous mutations. Sanger sequencing confirmed that the c.884C > T mutation was inherited from the mother, while the c.741+3_741+6delAAGT mutation was derived from the father. Prediction of amino acid function suggested that both mutations were pathogenic mutations. In conclusion, we found a new pathogenic complex heterozygous mutation of the TMC1 gene, which enriched the mutation spectrum of the TMC1 gene and provided a basis for genetic counseling and prenatal diagnosis of ARNSHL.

Robustness quantification of a mutant library screen revealed key genetic markers in yeast.

BACKGROUND: Microbial robustness is crucial for developing cell factories that maintain consistent performance in a challenging environment such as large-scale bioreactors. Although tools exist to assess and understand robustness at a phenotypic level, the underlying metabolic and genetic mechanisms are not well defined, which limits our ability to engineer more strains with robust functions. RESULTS: This study encompassed four steps. (I) Fitness and robustness were analyzed from a published dataset of yeast mutants grown in multiple environments. (II) Genes and metabolic processes affecting robustness or fitness were identified, and 14 of these genes were deleted in Saccharomyces cerevisiae CEN.PK113-7D. (III) The mutants bearing gene deletions were cultivated in three perturbation spaces mimicking typical industrial processes. (IV) Fitness and robustness were determined for each mutant in each perturbation space. We report that robustness varied according to the perturbation space. We identified genes associated with increased robustness such as MET28, linked to sulfur metabolism; as well as genes associated with decreased robustness, including TIR3 and WWM1, both involved in stress response and apoptosis. CONCLUSION: The present study demonstrates how phenomics datasets can be analyzed to reveal the relationship between phenotypic response and associated genes. Specifically, robustness analysis makes it possible to study the influence of single genes and metabolic processes on stable microbial performance in different perturbation spaces. Ultimately, this information can be used to enhance robustness in targeted strains.

Two thiosemicarbazones derived from 1-indanone as potent non-nucleoside inhibitors of bovine viral diarrhea virus of different genotypes and biotypes.

Bovine viral diarrhea virus (BVDV) is a widespread pathogen of cattle and other mammals that causes major economic losses in the livestock industry. N4-TSC and 6NO2-TSC are two thiosemicarbazones derived from 1-indanone that exhibit anti-BVDV activity in vitro. These compounds selectively inhibit BVDV and are effective against both cytopathic and non-cytopathic BVDV-1 and BVDV-2 strains. We confirmed that N4-TSC acts at the onset of viral RNA synthesis, as previously reported for 6NO2-TSC. Moreover, resistance selection and characterization showed that N4-TSCR mutants were highly resistant to N4-TSC but remained susceptible to 6NO2-TSC. In contrast, 6NO2-TSCR mutants were resistant to both compounds. Additionally, mutations N264D and A392E were found in the viral RNA-dependent RNA polymerase (RdRp) of N4-TSCR mutants, whereas I261 M was found in 6NO2-TSCR mutants. These mutations lay in a hydrophobic pocket within the fingertips region of BVDV RdRp that has been described as a "hot spot" for BVDV non-nucleoside inhibitors.

BRCA1 secondary splice-site mutations drive exon-skipping and PARP inhibitor resistance.

PARP inhibitor (PARPi) therapy has transformed outcomes for patients with homologous recombination DNA repair (HRR) deficient ovarian cancers, for example those with BRCA1 or BRCA2 gene defects. Unfortunately, PARPi resistance is common. Multiple resistance mechanisms have been described, including secondary mutations that restore the HR gene reading frame. BRCA1 splice isoforms △11 and △11q can contribute to PARPi resistance by splicing out the mutation-containing exon, producing truncated, partially functional proteins. However, the clinical impacts and underlying drivers of BRCA1 exon skipping are not fully understood.We analyzed nine ovarian and breast cancer patient derived xenografts (PDX) with BRCA1 exon 11 frameshift mutations for exon skipping and therapy response, including a matched PDX pair derived from a patient pre- and post-chemotherapy/PARPi. BRCA1 exon 11 skipping was elevated in PARPi resistant PDX tumors. Two independent PDX models acquired secondary BRCA1 splice site mutations (SSMs) that drive exon skipping, confirmed using qRT-PCR, RNA sequencing, immunoblotting and minigene modelling. CRISPR/Cas9-mediated disruption of splicing functionally validated exon skipping as a mechanism of PARPi resistance. SSMs were also enriched in post-PARPi ovarian cancer patient cohorts from the ARIEL2 and ARIEL4 clinical trials.Few PARPi resistance mechanisms have been confirmed in the clinical setting. While secondary/reversion mutations typically restore a gene's reading frame, we have identified secondary mutations in patient cohorts that hijack splice sites to enhance mutation-containing exon skipping, resulting in the overexpression of BRCA1 hypomorphs, which in turn promote PARPi resistance. Thus, BRCA1 SSMs can and should be clinically monitored, along with frame-restoring secondary mutations.

Clinical and Molecular Analysis in Patients with Peutz-Jeghers Syndrome.

Peutz-Jeghers syndrome (PJS) is a rare hereditary disorder linked to increased cancer risk due to specific genetic variants in the STK11 gene. This study aimed to assess disease manifestations, genetic profiles, and genotype-phenotype correlations in PJS patients. Twenty patients from 14 families with PJS who were followed up at our clinic between 2011 and 2021 were included. Genetic susceptibility to hereditary cancers was assess-ed using targeted next-generation sequencing (NGS) and multiplex ligation-dependent probe amplification (MLPA) of the STK11 gene. Clinical data were also collected and analyzed in conjunction with the genetic findings. Initial symptoms appeared around 18.9 years, predominantly abdominal pain and intussusception. Mucocutaneous lesions were found in 85%, and hamartomatous polyps in 90%. Dysplastic polyps were found in 4 patients, with 3 cases of malignancy. Nextgeneration sequencing identified 11 pathogenic and 3 likely pathogenic mutations, including 3 novel STK11 variants (LRG_319: c.598- 8_601del, LRG_319: c.708_718del, and LRG_319: c.146_147del). Next-generation sequencing diagnostic rate was 78.5% (11/14), and the overall diagnostic rate with NGS and MLPA studies was 85.7% (12/14). Patients without STK11 mutations had later symptom onset and potentially lower cancer risk. Truncated mutations are associated with earlier symptoms and elevated cancer risk. This is the first PJS case series in Turkey using the NGS and MLPA methods. It reports 3 novel mutations and emphasizes the genotype-phenotype relationship of PJS. With further studies, the genotype-phenotype relationship of STK11 variants will be better understood.

Revisiting Dominance in Population Genetics.

Dominance refers to the effect of a heterozygous genotype relative to that of the two homozygous genotypes. The degree of dominance of mutations for fitness can have a profound impact on how deleterious and beneficial mutations change in frequency over time as well as on the patterns of linked neutral genetic variation surrounding such selected alleles. Since dominance is such a fundamental concept, it has received immense attention throughout the history of population genetics. Early work from Fisher, Wright, and Haldane focused on understanding the conceptual basis for why dominance exists. More recent work has attempted to test these theories and conceptual models by estimating dominance effects of mutations. However, estimating dominance coefficients has been notoriously challenging and has only been done in a few species in a limited number of studies. In this review, we first describe some of the early theoretical and conceptual models for understanding the mechanisms for the existence of dominance. Second, we discuss several approaches used to estimate dominance coefficients and summarize estimates of dominance coefficients. We note trends that have been observed across species, types of mutations, and functional categories of genes. By comparing estimates of dominance coefficients for different types of genes, we test several hypotheses for the existence of dominance. Lastly, we discuss how dominance influences the dynamics of beneficial and deleterious mutations in populations and how the degree of dominance of deleterious mutations influences the impact of inbreeding on fitness.

Impact of reference design on estimating SARS-CoV-2 lineage abundances from wastewater sequencing data.

BACKGROUND: Sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA from wastewater samples has emerged as a valuable tool for detecting the presence and relative abundances of SARS-CoV-2 variants in a community. By analyzing the viral genetic material present in wastewater, researchers and public health authorities can gain early insights into the spread of virus lineages and emerging mutations. Constructing reference datasets from known SARS-CoV-2 lineages and their mutation profiles has become state-of-the-art for assigning viral lineages and their relative abundances from wastewater sequencing data. However, selecting reference sequences or mutations directly affects the predictive power. RESULTS: Here, we show the impact of a mutation- and sequence-based reference reconstruction for SARS-CoV-2 abundance estimation. We benchmark 3 datasets: (i) synthetic "spike-in"' mixtures; (ii) German wastewater samples from early 2021, mainly comprising Alpha; and (iii) samples obtained from wastewater at an international airport in Germany from the end of 2021, including first signals of Omicron. The 2 approaches differ in sublineage detection, with the marker mutation-based method, in particular, being challenged by the increasing number of mutations and lineages. However, the estimations of both approaches depend on selecting representative references and optimized parameter settings. By performing parameter escalation experiments, we demonstrate the effects of reference size and alternative allele frequency cutoffs for abundance estimation. We show how different parameter settings can lead to different results for our test datasets and illustrate the effects of virus lineage composition of wastewater samples and references. CONCLUSIONS: Our study highlights current computational challenges, focusing on the general reference design, which directly impacts abundance allocations. We illustrate advantages and disadvantages that may be relevant for further developments in the wastewater community and in the context of defining robust quality metrics.

INSIGHTS INTO THE ORIGINS OF CYSTIC FIBROSIS LUNG DISEASE.

In this paper, I will discuss recent studies using a cystic fibrosis pig model to better understand the origins of cystic fibrosis lung disease. Specifically, I will review our work investigating how loss of the cystic fibrosis transmembrane conductance regulator function (CFTR) impairs mucociliary transport in the cystic fibrosis airway. These studies reveal new insights into the early, underlying mechanisms of cystic fibrosis lung disease and could lead to novel therapeutic interventions.

Cellular transitions during cranial suture establishment in zebrafish.

Cranial sutures separate neighboring skull bones and are sites of bone growth. A key question is how osteogenic activity is controlled to promote bone growth while preventing aberrant bone fusions during skull expansion. Using single-cell transcriptomics, lineage tracing, and mutant analysis in zebrafish, we uncover key developmental transitions regulating bone formation at sutures during skull expansion. In particular, we identify a subpopulation of mesenchyme cells in the mid-suture region that upregulate a suite of genes including BMP antagonists (e.g. grem1a) and pro-angiogenic factors. Lineage tracing with grem1a:nlsEOS reveals that this mid-suture subpopulation is largely non-osteogenic. Moreover, combinatorial mutation of BMP antagonists enriched in this mid-suture subpopulation results in increased BMP signaling in the suture, misregulated bone formation, and abnormal suture morphology. These data reveal establishment of a non-osteogenic mesenchyme population in the mid-suture region that restricts bone formation through local BMP antagonism, thus ensuring proper suture morphology.

Inherited mitochondrial dysfunction triggered by OPA1 mutation impacts the sensory innervation fibre identity, functionality and regenerative potential in the cornea.

Mitochondrial dysfunctions are detrimental to organ metabolism. The cornea, transparent outmost layer of the eye, is prone to environmental aggressions, such as UV light, and therefore dependent on adequate mitochondrial function. While several reports have linked corneal defects to mitochondrial dysfunction, the impact of OPA1 mutation, known to induce such dysfunction, has never been studied in this context. We used the mouse line carrying OPA1delTTAG mutation to investigate its impact on corneal biology. To our surprise, neither the tear film composition nor the corneal epithelial transcriptomic signature were altered upon OPA1 mutation. However, when analyzing the corneal innervation, we discovered an undersensitivity of the cornea upon the mutation, but an increased innervation volume at 3 months. Furthermore, the fibre identity changed with a decrease of the SP + axons. Finally, we demonstrated that the innervation regeneration was less efficient and less functional in OPA1+/- corneas. Altogether, our study describes the resilience of the corneal epithelial biology, reflecting the mitohormesis induced by the OPA1 mutation, and the adaptation of the corneal innervation to maintain its functionality despite its morphogenesis defects. These findings will participate to a better understanding of the mitochondrial dysfunction on peripheral innervation.

Effect of mutation of phaC on carbon supply, extracellular polysaccharide production, and pathogenicity of Xanthomonas oryzae pv. oryzae.

Xanthomonas oryzae pv. oryzae (Xoo) causes bacterial blight in rice. Polyhydroxyalkanoates (PHAs) consitute a diverse group of biopolyesters synthesized by bacteria under nutrient-limited conditions. The phaC gene is important for PHA polymerization. We investigated the effects of phaC gene mutagensis in Xoo strain PXO99A. The phaC gene knock-out mutant exhibited reduced swarming ability relative to that of the wild-type. Under conditions where glucose was the sole sugar source, extracellular polysaccharide (EPS) production by ΔphaC declined by 44.8%. ΔphaC showed weak hypersensitive response (HR) induction in the leaves of non-host Nicotiana tabacum, concomitant with downregulation of hpa1 gene expression. When inoculated in rice leaves by the leaf-clipping method, ΔphaC displayed reduced virulence in terms of lesion length compared with the wild-type strain. The complemented strain showed no significant difference from the wild-type strain, suggesting that the deletion of phaC in Xoo induces significant alterations in various physiological and biological processes. These include bacterial swarming ability, EPS production, transcription of hrp genes, and glucose metabolism. These changes are intricately linked to the energy utilization and virulence of Xoo during plant infection. These findings revealed involvement of phaC in Xoo is in the maintaining carbon metabolism by functioning in the PHA metabolic pathway.

Integrated genomic/epigenomic analysis stratifies subtypes of clear cell ovarian carcinoma, highlighting their cellular origin.

The cellular origin of clear cell ovarian carcinoma (CCOC), a major histological subtype of ovarian carcinoma remains elusive. Here, we explored the candidate cellular origin and identify molecular subtypes using integrated genomic/epigenomic analysis. We performed whole exome-sequencing, microarray, and DNA methylation array in 78 CCOC samples according to the original diagnosis. The findings revealed that ARID1A and/or PIK3CA mutations were mutually exclusive with DNA repair related genes, including TP53, BRCA1, and ATM. Clustering of CCOC and other ovarian carcinomas (n = 270) with normal tissues from the fallopian tube, ovarian surface epithelium, endometrial epithelium, and pelvic peritoneum mesothelium (PPM) in a methylation array showed that major CCOC subtypes (with ARID1A and/or PIK3CA mutations) were associated with the PPM-lile cluster (n = 64). This cluster was sub-divided into three clusters: (1) mismatch repair (MMR) deficient with tumor mutational burden-high (n = 2), (2) alteration of ARID1A (n = 51), and (3) ARID1A wild-type (n = 11). The remaining samples (n = 14) were subdivided into (4) ovarian surface epithelium-like (n = 11) and (5) fallopian tube-like (considered as high-grade serous histotype; n = 3). Among these, subtypes (1-3) and others (4 and 5) were found to be associated with immunoreactive signatures and epithelial-mesenchymal transition, respectively. These results contribute to the stratification of CCOC into biological subtypes.

The genetic differences between types 1 and 2 in von Hippel-Lindau syndrome: comprehensive meta-analysis.

BACKGROUND: Patients with von Hippel-Lindau (VHL) disease are at risk of developing tumors in the eye, brain, kidney, adrenal gland, and other organs based on their gene mutations. The VHL tumor suppressor gene contains pathogenic variants responsible for these events. This meta-analysis aims to investigate the genetic differences among the various types of VHL syndrome and their correlation with the location of mutations (exons and domains) in the VHL gene. METHOD: Papers eligible for publication until September 2023 were identified using the electronic databases of PubMed, Google Scholar, Scopus, and EMBASE. The Random Effect model was utilized to evaluate the genetic differences between type 1 and type 2 VHL syndromes. RESULTS: The prevalence of missense mutations (MSs) was found to be 58.9% in type 1, while it was 88.1% in type 2. Interestingly, the probability of observing MSs in type 1 was 0.42 times lower compared to type 2. The mutation hotspots of the VHL gene were R167Q/W, Y98H, R238W, and S65L, respectively. Although type 2 had a high presentation of Y98H and R238W, it did not have a higher S65L than type 1. The analysis demonstrated a statistically significant higher prevalence of truncated mutations (PTMs) in type 1. Among type 1, large/complete deletions (L/C DELs) were found in 16.9% of cases, whereas in type 2 only 3.7%. This difference was statistically significant with a p-value < 0.001. Overall, the probability of identifying mutations in domain 2 compared to domain 1 was found to be 2.13 times higher in type 1 (p-value < 0.001). Furthermore, the probability of detecting exon 1 in comparison with observing exon 2 in type 1 was 2.11 times higher than type 2 and revealed a statistically significant result (p-value < 0.001). The detection of exon 2 was 2.18 times higher in type 1 (p-value < 0.001). In addition, the likelihood of discovering exon 2 compared with others was significantly lower in type 1 compared with type 2 VHL (OR = 0.63, p-value = 0.015). CONCLUSIONS: We have revealed a comprehensive genetic difference between types 1 and 2 of VHL syndrome. The significant differences in MS, PTMs, L/C DELs, and the location of the mutations between type 1 and type 2 VHL patients in the Asian, European, and American populations emphasize the genetic heterogeneity of the syndrome. These findings may pave the way for the diagnosis, treatment, and further investigation of the mechanisms behind this complex genetic disorder.

The GC-content at the 5' ends of human protein-coding genes is undergoing mutational decay.

BACKGROUND: In vertebrates, most protein-coding genes have a peak of GC-content near their 5' transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations. RESULTS: Our data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5' end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5' end of protein-coding is increasing. We show that these patterns extend into the 5' end of the open reading frame, thus impacting synonymous codon position choices. CONCLUSIONS: Our results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species.