Using Drosophila melanogaster as a Model for Genotoxic Chemical Mutational Studies with a New Program, SnpSift

De novo mutations promote inflammation in children with STAT3 gain-of-function syndrome by affecting IL-1β expression

STAT3 gain-of-function syndrome, characterized by early-onset autoimmunity and primary immune regulatory disorder, remains poorly understood in terms of its immunological mechanisms. We employed whole-genome sequencing of familial trios to elucidate the pivotal role of de novo mutations in genetic diseases. We identified 37 high-risk pathogenic loci affecting 23 genes, including a novel STAT3 c.508G>A mutation. We also observed significant down-regulation of pathogenic genes in affected individuals, potentially associated with inflammatory responses regulated by PTPN14 via miR378c. These findings enhance our understanding of the pathogenesis of STAT3 gain-of-function syndrome and suggest potential therapeutic strategies. Notably, combined JAK inhibitors and IL-6R antagonists may offer promising treatment avenues for mitigating the severity of STAT3 gain-of-function syndrome.

The drought-induced plasticity of mineral nutrients contributes to drought tolerance discrimination in durum wheat

Drought is a major challenge for the cultivation of durum wheat, a crucial crop for global food security. Plants respond to drought by adjusting their mineral nutrient profiles to cope with water scarcity, showing the importance of nutrient plasticity for plant acclimation and adaptation to diverse environments. Therefore, it is essential to understand the genetic basis of mineral nutrient profile plasticity in durum wheat under drought stress to select drought-tolerant varieties. The research study investigated the responses of different durum wheat genotypes to severe drought stress at the seedling stage. The study employed an ionomic, molecular, biochemical and physiological approach to shed light on distinct behaviors among different genotypes. The drought tolerance of SVEMS16, SVEVO, and BULEL was related to their capacity of maintaining or increasing nutrient's accumulation, while the limited nutrient acquisition capability of CRESO and S.CAP likely resulted in their susceptibility to drought. The study highlighted the importance of macronutrients such as SO42-, NO3-, PO43-, and K+ in stress resilience and identified variant-containing genes potentially influencing nutritional variations under drought. These findings provide valuable insights for further field studies to assess the drought tolerance of durum wheat genotypes across various growth stages, ultimately ensuring food security and sustainable production in the face of changing environmental conditions.

QTL analysis to identify genes involved in the trade-off between silk protein synthesis and larva-pupa transition in silkworms

BACKGROUND

Insect-based food and feed are increasingly attracting attention. As a domesticated insect, the silkworm (Bombyx mori) has a highly nutritious pupa that can be easily raised in large quantities through large-scale farming, making it a highly promising source of food. The ratio of pupa to cocoon (RPC) refers to the proportion of the weight of the cocoon that is attributed to pupae, and is of significant value for edible utilization, as a higher RPC means a higher ratio of conversion of mulberry leaves to pupa. In silkworm production, there is a trade-off between RPC and cocoon shell ratiao(CSR), which refers the ratio of silk protein to the entire cocoon, during metamorphosis process. Understanding the genetic basis of this balance is crucial for breeding edible strains with a high RPC and further advancing its use as feed.

RESULTS

Using QTL-seq, we identified a quantitative trait locus (QTL) for the balance between RPC and CSR that is located on chromosome 11 and covers a 9,773,115-bp region. This locus is an artificial selection hot spot that contains ten non-overlapping genomic regions under selection that were involved in the domestication and genetic breeding processes. These regions include 17 genes, nine of which are highly expressed in the silk gland, which is a vital component in the trade-off between RPC and CSR. These genes are annotate with function related with epigenetic modifications and the regulation of DNA replication et al. We identified one and two single nucleotide polymorphisms (SNPs) in the exons of teh KWMTBOMO06541 and KWMTBOMO06485 genes that result in amino acid changes in the protein domains. These SNPs have been strongly selected for during the domestication process. The KWMTBOMO06485 gene encodes the Bombyx mori (Bm) tRNA methyltransferase (BmDnmt2) and its knockout results in a significant change in the trade-off between CSR and RPC in both sexes.

CONCLUSIONS

Taken together, our results contribute to a better understanding of the genetic basis of RPC and CSR. The identified QTL and genes that affect RPC can be used for marker-assisted and genomic selection of silkworm strains with a high RPC. This will further enhance the production efficiency of silkworms and of closely-related insects for edible and feed purposes.

Genomic and fitness consequences of a near-extinction event in the northern elephant seal

Understanding the genetic and fitness consequences of anthropogenic bottlenecks is crucial for biodiversity conservation. However, studies of bottlenecked populations combining genomic approaches with fitness data are rare. Theory predicts that severe bottlenecks deplete genetic diversity, exacerbate inbreeding depression and decrease population viability. However, actual outcomes are complex and depend on how a species' unique demography affects its genetic load. We used population genetic and veterinary pathology data, demographic modelling, whole-genome resequencing and forward genetic simulations to investigate the genomic and fitness consequences of a near-extinction event in the northern elephant seal. We found no evidence of inbreeding depression within the contemporary population for key fitness components, including body mass, blubber thickness and susceptibility to parasites and disease. However, we detected a genomic signature of a recent extreme bottleneck (effective population size = 6; 95% confidence interval = 5.0-7.5) that will have purged much of the genetic load, potentially leading to the lack of observed inbreeding depression in our study. Our results further suggest that deleterious genetic variation strongly impacted the post-bottleneck population dynamics of the northern elephant seal. Our study provides comprehensive empirical insights into the intricate dynamics underlying species-specific responses to anthropogenic bottlenecks.

Phage-antibiotic synergy suppresses resistance emergence of Klebsiella pneumoniae by altering the evolutionary fitness

Phage-antibiotic synergy (PAS) represents a superior treatment strategy for pathogen infections with less probability of resistance development. Here, we aim to understand the molecular mechanism by which PAS suppresses resistance in terms of population evolution. A novel hypervirulent Klebsiella pneumoniae (KP) phage H5 was genetically and structurally characterized. The combination of H5 and ceftazidime (CAZ) showed a robust synergistic effect in suppressing resistance emergence. Single-cell Raman analysis showed that the phage-CAZ combination suppressed bacterial metabolic activities, contrasting with the upregulation observed with phage alone. The altered population evolutionary trajectory was found to be responsible for the contrasting metabolic activities under different selective pressures, resulting in pleiotropic effects. A pre-existing wcaJ point mutation (wcaJG949A) was exclusively selected by H5, conferring a fitness advantage and up-regulated activity of carbohydrate metabolism, but also causing a trade-off between phage resistance and collateral sensitivity to CAZ. The wcaJ point mutation was counter-selected by H5-CAZ, inducing various mutations in galU that imposed evolutionary disadvantages with higher fitness costs, and suppressed carbohydrate metabolic activity. H5 and H5-CAZ treatments resulted in opposite effects on the transcriptional activity of the phosphotransferase system and the ascorbate and aldarate metabolism pathway, suggesting potential targets for phage resistance suppression. Our study reveals a novel mechanism of resistance suppression by PAS, highlighting how the complexity of bacterial adaptation to selective pressures drives treatment outcomes.

IMPORTANCE

Phage-antibiotic synergy (PAS) has been recently proposed as a superior strategy for the treatment of multidrug-resistant pathogens to effectively reduce bacterial load and slow down both phage and antibiotic resistance. However, the underlying mechanisms of resistance suppression by PAS have been poorly and rarely been studied. In this study, we tried to understand how PAS suppresses the emergence of resistance using a hypervirulent Klebsiella pneumoniae (KP) strain and a novel phage H5 in combination with ceftazidime (CAZ) as a model. Our study reveals a novel mechanism by which PAS drives altered evolutionary trajectory of bacterial populations, leading to suppressed emergence of resistance. The findings advance our understanding of how PAS suppresses the emergence of resistance, and are imperative for optimizing the efficacy of phage-antibiotic therapy to further improve clinical outcomes.

Efficacy of Pembrolizumab and Biomarker Analysis in Patients with WGS-Based Intermediate to High Tumor Mutational Load: Results from the Drug Rediscovery Protocol

PURPOSE

To evaluate the efficacy of pembrolizumab across multiple cancer types harboring different levels of whole-genome sequencing-based tumor mutational load (TML; total of nonsynonymous mutations across the genome) in patients included in the Drug Rediscovery Protocol (NCT02925234).

PATIENTS AND METHODS

Patients with solid, treatment-refractory, microsatellite-stable tumors were enrolled in cohort A: breast cancer cohort harboring a TML of 140 to 290, cohort B: tumor-agnostic cohort harboring a TML of 140 to 290, and cohort C: tumor-agnostic cohort harboring a TML >290. Patients received pembrolizumab 200 mg every 3 weeks. The primary endpoint was clinical benefit [CB; objective response or stable disease (SD) ≥16 weeks]. Pretreatment tumor biopsies were obtained for whole-genome sequencing and RNA sequencing.

RESULTS

Seventy-two evaluable patients with 26 different histotypes were enrolled. The CB rate was 13% in cohort A [3/24 with partial response (PR)], 21% in cohort B (3/24 with SD; 2/24 with PR), and 42% in cohort C (4/24 with SD; 6/24 with PR). In cohort C, neoantigen burden estimates and expression of inflammation and innate immune biomarkers were significantly associated with CB. Similar associations were not identified in cohorts A and B. In cohort A, CB was significantly associated with mutations in the chromatin remodeling gene PBRM1, whereas in cohort B, CB was significantly associated with expression of MICA/MICB and butyrophilins. CB and clonal TML were not significantly associated.

CONCLUSIONS

Although pembrolizumab lacked activity in cohort A, cohorts B and C met the study's primary endpoint. Further research is warranted to refine the selection of patients with tumors harboring lower TMLs and may benefit from a focus on innate immunity. See related commentary by Hsu and Yen, p. 3652.

Genetic drift drives faster-Z evolution in the salmon louse Lepeophtheirus salmonis

How sex chromosomes evolve compared to autosomes remains an unresolved question in population genetics. Most studies focus on only a handful of taxa, resulting in uncertainty over whether observed patterns reflect general processes or idiosyncrasies in particular clades. For example, in female heterogametic (ZW) systems, bird Z chromosomes tend to evolve quickly but not adaptively, while in Lepidopterans they evolve adaptively, but not always quickly. To understand how these observations fit into broader evolutionary patterns, we explore Z chromosome evolution outside of these two well-studied clades. We utilize a publicly available genome, gene expression, population, and outgroup data in the salmon louse Lepeophtheirus salmonis, an important agricultural pest copepod. We find that the Z chromosome is faster evolving than autosomes, but that this effect is driven by increased drift rather than adaptive evolution. Due to high rates of female reproductive failure, the Z chromosome exhibits a slightly lower effective population size than the autosomes which is nonetheless to decrease efficiency of hemizygous selection acting on the Z. These results highlight the usefulness of organismal life history in calibrating population genetic expectations and demonstrate the value of the ever-expanding wealth of publicly available data to help resolve outstanding evolutionary questions.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

Genetic heterogeneity in the Salmonella Typhi Vi capsule locus: a population genomic study from Fiji

Typhoid fever is endemic in many parts of the world and remains a major public health concern in tropical and sub-tropical developing nations, including Fiji. To address high rates of typhoid fever, the Northern Division of Fiji implemented a mass vaccination with typhoid conjugate vaccine (Vi-polysaccharide conjugated to tetanus toxoid) as a public health control measure in 2023. In this study we define the genomic epidemiology of Salmonella Typhi in the Northern Division prior to island-wide vaccination, sequencing 85% (n=419) of the total cases from the Northern and Central Divisions of Fiji that occurred in the period 2017-2019. We found elevated rates of nucleotide polymorphisms in the tviD and tviE genes (responsible for Vi-polysaccharide synthesis) relative to core genome levels within the Fiji endemic S. Typhi genotype 4.2. Expansion of these findings within a globally representative database of 12 382 S. Typhi (86 genotyphi clusters) showed evidence of convergent evolution of the same tviE mutations across the S. Typhi population, indicating that tvi selection has occurred both independently and globally. The functional impact of tvi mutations on the Vi-capsular structure and other phenotypic characteristics are not fully elucidated, yet commonly occurring tviE polymorphisms localize adjacent to predicted active site residues when overlayed against the predicted TviE protein structure. Given the central role of the Vi-polysaccharide in S. Typhi biology and vaccination, further integrated epidemiological, genomic and phenotypic surveillance is required to determine the spread and functional implications of these mutations.

Evolutionary dynamics in gut-colonizing Candida glabrata during caspofungin therapy: Emergence of clinically important mutations in sphingolipid biosynthesis

Invasive fungal infections are associated with high mortality, which is exacerbated by the limited antifungal drug armamentarium and increasing antifungal drug resistance. Echinocandins are a frontline antifungal drug class targeting β-glucan synthase (GS), a fungal cell wall biosynthetic enzyme. Echinocandin resistance is generally low but increasing in species like Candida glabrata, an opportunistic yeast pathogen colonizing human mucosal surfaces. Mutations in GS-encoding genes (FKS1 and FKS2 in C. glabrata) are strongly associated with clinical echinocandin failure, but epidemiological studies show that other, as yet unidentified factors also influence echinocandin susceptibility. Furthermore, although the gut is known to be an important reservoir for emergence of drug-resistant strains, the evolution of resistance is not well understood. Here, we studied the evolutionary dynamics of C. glabrata colonizing the gut of immunocompetent mice during treatment with caspofungin, a widely-used echinocandin. Whole genome and amplicon sequencing revealed rapid genetic diversification of this C. glabrata population during treatment and the emergence of both drug target (FKS2) and non-drug target mutations, the latter predominantly in the FEN1 gene encoding a fatty acid elongase functioning in sphingolipid biosynthesis. The fen1 mutants displayed high fitness in the gut specifically during caspofungin treatment and contained high levels of phytosphingosine, whereas genetic depletion of phytosphingosine by deletion of YPC1 gene hypersensitized the wild type strain to caspofungin and was epistatic to fen1Δ. Furthermore, high resolution imaging and mass spectrometry showed that reduced caspofungin susceptibility in fen1Δ cells was associated with reduced caspofungin binding to the plasma membrane. Finally, we identified several different fen1 mutations in clinical C. glabrata isolates, which phenocopied the fen1Δ mutant, causing reduced caspofungin susceptibility. These studies reveal new genetic and molecular determinants of clinical caspofungin susceptibility and illuminate the dynamic evolution of drug target and non-drug target mutations reducing echinocandin efficacy in patients colonized with C. glabrata.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Association of pathogenic determinants of Fusobacterium necrophorum with bacteremia, and Lemierre's syndrome

Fusobacterium necrophorum is a Gram-negative anaerobic bacterium responsible for localized infections of the oropharynx that can evolve into bacteremia and/or septic thrombophlebitis of the jugular vein or peritonsillar vein, called Lemierre's syndrome. To identify microbial genetic determinants associated with the severity of this life-threatening disease, 70 F. necrophorum strains were collected and grouped into two categories according to the clinical presentation: (i) localized infection, (ii) bacteremia with/without Lemierre's syndrome. Comparative genomic analyses revealed two clades with distinct genetic content, one clade being significantly enriched with isolates from subjects with bacteremia. To identify genetic determinants contributing to F. necrophorum pathogenicity, genomic islands and virulence factor orthogroups (OVFs) were predicted. The presence/absence profiles of OVFs did not group isolates according to their clinical category, but rather according to their phylogeny. However, a variant of lktA, a key virulence factor, with a frameshift deletion that results in two open reading frames, was associated with bacteremia. Moreover, a genome-wide association study identified three orthogroups associated with bacteremic strains: (i) cas8a1, (ii) a sodium/solute symporter, and (iii) a POP1 domain-containing protein. Further studies must be performed to assess the functional impact of lktA mutation and of these orthogroups on the physiopathological mechanisms of F. necrophorum infections.

SeqVerify: An accessible analysis tool for cell line genomic integrity, contamination, and gene editing outcomes

Over the last decade, advances in genome editing and pluripotent stem cell (PSC) culture have let researchers generate edited PSC lines to study a wide variety of biological questions. However, abnormalities in cell lines such as aneuploidy, mutations, on-target and off-target editing errors, and microbial contamination can arise during PSC culture or due to undesired editing outcomes. The ongoing decline of next-generation sequencing prices has made whole-genome sequencing (WGS) a promising option for detecting these abnormalities. However, this approach has been held back by a lack of easily usable data analysis software. Here, we present SeqVerify, a computational pipeline designed to take raw WGS data and a list of intended genome edits, and verify that the edits are present and that there are no abnormalities. We anticipate that SeqVerify will be a useful tool for researchers generating edited PSCs, and more broadly, for cell line quality control in general.

Clinical and Genomic Perspective of SARS CoV-2 Infection in Liver Disease Patients: A Single-Centre Retrospective Study

The limited literature on the clinical course of COVID-19 among patients with underlying liver disease (LD) is available from India. The present study aimed to evaluate the clinical and mutational profile of SARS-CoV-2 among LD cases. This was a retrospective study including admitted LD cases in whom SARS-CoV-2 RT-PCR testing was performed. Complete demographic and clinical details were retrieved from Hospital Information System. Detailed mutational analysis was performed by comparing LD COVID-19 positive study group, i.e. LD-CoV(+) with COVID-19 positive outpatients without any underlying LD as control, i.e. NLD-CoV(+). Out of 232 enrolled LD cases, 137 (59.1%) were LD-CoV(+). LD cases with existing co-morbidities were affected more (P = 0.002) and had 2.29 times (OR 2.29, CI 95%, 1.25-4.29) higher odds of succumbing to COVID-19 (P = 0.006). On multivariate regression analysis, ascites (P = 0.05), severe COVID-19 pneumonia (P = 0.046), and an increased levels of bilirubin (P = 0.005) and alkaline phosphatase (P = 0.003) were found to be associated with adverse outcome in LD-CoV(+).On mutational analysis, we found certain differences between LD- and NLD-CoV(+) infected with Delta [LD- and NLD-CoV (+ /D)] and Omicron [LD- and NLD-CoV(+/O)]. More mutations were shared between LD- and NLD-CoV(+/O) compared to LD- and NLD-CoV(+/D). There were differences in prevalence of indel mutations specific to LD-CoV ( +) for both Delta and Omicron. Moreover, we also reported an interesting genic bias between LD- and NLD-CoV( +) in harbouring deleterious/tolerated mutations. To conclude, LD cases with comorbidities were affected more and had higher odds of mortality due to COVID-19. The definite difference between LD- and NLD-CoV(+) groups with respect to frequency of harboured mutations and an inherent genic bias between them is of noteworthy importance.

A Novel Pathogenic TUBA1A Variant in a Croatian Infant Is Linked to a Severe Tubulinopathy with Walker-Warburg-like Features

Tubulinopathies are associated with malformations of cortical development but not Walker-Warburg Syndrome. Intensive monitoring of a Croatian infant presenting as Walker-Warburg Syndrome in utero began at 21 weeks due to increased growth of cerebral ventricles and foetal biparietal diameter. Monitoring continued until Caesarean delivery at 34 weeks where the infant was eutrophic. Clinical assessment of a progressive neurological disorder of unknown aetiology found a macrocephalic head and markedly hypoplastic genitalia with a micropenis. Neurological examination showed generalized hypotonia with very rare spontaneous movements, hypotonia-induced respiratory insufficiency and ventilator dependence, and generalized myoclonus intensifying during manipulation. With clinical features of hypotonia, lissencephaly, and brain malformations, Walker-Warburg Syndrome was suspected; however, eye anomalies were absent. Genetic trio analysis via whole-exome sequencing only identified a novel de novo mutation in the TUBA1A gene (NM_006009.4:c.848A>G; NP_006000.2:p.His283Arg) in the infant, who died at 2 months of age, as the likely cause. We report a previously unpublished, very rare heterozygous TUBA1A mutation with clinical features of macrocephaly and hypoplastic genitalia which have not previously been associated with the gene. The absence of eye phenotypes or mutations in Walker-Warburg-associated genes confirm this as not a new presentation of Walker-Warburg Syndrome but a novel TUBA1A tubulinopathy for neonatologists to be aware of.

Sequential Analysis of cfDNA Reveals Clonal Evolution in Patients with Neuroblastoma Receiving ALK-Targeted Therapy

PURPOSE

The study of cell-free DNA (cfDNA) enables sequential analysis of tumor cell-specific genetic alterations in patients with neuroblastoma.

EXPERIMENTAL DESIGN

Eighteen patients with relapsing neuroblastoma having received lorlatinib, a third-generation ALK inhibitor, were identified (SACHA national registry and/or in the institution). cfDNA was analyzed at relapse for nine patients and sequentially for five patients (blood/bone marrow plasma) by performing whole-genome sequencing library construction followed by ALK-targeted ddPCR of the hotspot mutations [F1174L, R1275Q, and I1170N; variant allele fraction (VAF) detection limit 0.1%] and whole-exome sequencing (WES) to evaluate disease burden and clonal evolution, following comparison with tumor/germline WES.

RESULTS

Overall response rate to lorlatinib was 33% (CI, 13%-59%), with response observed in 6/10 cases without versus 0/8 cases with MYCN amplification (MNA). ALK VAFs correlated with the overall clinical disease status, with a VAF < 0.1% in clinical remission, versus higher VAFs (>30%) at progression. Importantly, sequential ALK ddPCR detected relapse earlier than clinical imaging. cfDNA WES revealed new SNVs, not seen in the primary tumor, in all instances of disease progression after lorlatinib treatment, indicating clonal evolution, including alterations in genes linked to tumor aggressivity (TP53) or novel targets (EGFR). Gene pathway analysis revealed an enrichment for genes targeting cell differentiation in emerging clones, and cell adhesion in persistent clones. Evidence of clonal hematopoiesis could be observed in follow-up samples.

CONCLUSIONS

We demonstrate the clinical utility of combining ALK cfDNA ddPCR for disease monitoring and cfDNA WES for the study of clonal evolution and resistance mechanisms in patients with neuroblastoma receiving ALK-targeted therapy.

Implementation of long-read sequencing for routine molecular diagnosis of familial mediterranean fever

Background

Long-read sequencing technology, widely used in research, is proving useful in clinical diagnosis, especially for infectious diseases. Despite recent advances, it hasn't been routinely applied to constitutional human diseases. Long-read sequencing detects intronic variants and phases variants, crucial for identifying recessive diseases.

Methods

We integrated long-read sequencing into the clinical diagnostic workflow for the MEFV gene, responsible for familial Mediterranean fever (FMF), using a Nanopore-based workflow. This involved long-range PCR amplification, native barcoding kit library preparation, GridION sequencing, and in-house bioinformatics. We compared this new workflow against our validated method using 39 patient samples and 3 samples from an external quality assessment scheme to ensure compliance with ISO15189 standards.

Results

Our evaluation demonstrated excellent performance, meeting ISO15189 requirements for reproducibility, repeatability, sensitivity, and specificity. Since October 2022, 150 patient samples were successfully analyzed with no failures. Among these samples, we identified 13 heterozygous carriers of likely pathogenic (LP) or pathogenic (P) variants, 1 patient with a homozygous LP/P variant in MEFV, and 4 patients with compound heterozygous variants.

Conclusion

This study represents the first integration of long-read sequencing for FMF clinical diagnosis, achieving 100 % sensitivity and specificity. Our findings highlight its potential to identify pathogenic variants without parental segregation analysis, offering faster, cost-effective, and accurate clinical diagnosis. This successful implementation lays the groundwork for future applications in other constitutional human diseases, advancing precision medicine.

Historical museum samples reveal signals of selection and drift in response to changing insecticide use in an agricultural pest moth

In response to environmental and human-imposed selective pressures, agroecosystem pests frequently undergo rapid evolution, with some species having a remarkable capacity to rapidly develop pesticide resistance. Temporal sampling of genomic data can comprehensively capture such adaptive changes over time, for example, by elucidating allele frequency shifts in pesticide resistance loci in response to different pesticides. Here, we leveraged museum specimens spanning over a century of collections to generate temporal contrasts between pre- and post-insecticide populations of an agricultural pest moth, Helicoverpa armigera. We used targeted exon sequencing of 254 samples collected across Australia from the pre-1950s (prior to insecticide introduction) to the 1990s, encompassing decades of changing insecticide use. Our sequencing approach focused on genes that are known to be involved in insecticide resistance, environmental sensation, and stress tolerance. We found an overall lack of spatial and temporal population structure change across Australia. In some decades (e.g., 1960s and 1970s), we found a moderate reduction of genetic diversity, implying stochasticity in evolutionary trajectories due to genetic drift. Temporal genome scans showed extensive evidence of selection following insecticide use, although the majority of selected variants were low impact. Finally, alternating trajectories of allele frequency change were suggestive of potential antagonistic pleiotropy. Our results provide new insights into recent evolutionary responses in an agricultural pest and show how temporal contrasts using museum specimens can improve mechanistic understanding of rapid evolution.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

SDHAF1 confers metabolic resilience to aging hematopoietic stem cells by promoting mitochondrial ATP production

Aging generally predisposes stem cells to functional decline, impairing tissue homeostasis. Here, we report that hematopoietic stem cells (HSCs) acquire metabolic resilience that promotes cell survival. High-resolution real-time ATP analysis with glucose tracing and metabolic flux analysis revealed that old HSCs reprogram their metabolism to activate the pentose phosphate pathway (PPP), becoming more resistant to oxidative stress and less dependent on glycolytic ATP production at steady state. As a result, old HSCs can survive without glycolysis, adapting to the physiological cytokine environment in bone marrow. Mechanistically, old HSCs enhance mitochondrial complex II metabolism during stress to promote ATP production. Furthermore, increased succinate dehydrogenase assembly factor 1 (SDHAF1) in old HSCs, induced by physiological low-concentration thrombopoietin (TPO) exposure, enables rapid mitochondrial ATP production upon metabolic stress, thereby improving survival. This study provides insight into the acquisition of resilience through metabolic reprogramming in old HSCs and its molecular basis to ameliorate age-related hematopoietic abnormalities.

A single base pair duplication in the SLC33A1 gene is associated with fetal losses and neonatal lethality in Manech Tête Rousse dairy sheep

We recently discovered that the Manech Tête Rousse (MTR) deficient homozygous haplotype 2 (MTRDHH2) probably carries a recessive lethal mutation in sheep. In this study, we fine-mapped this region through whole-genome sequencing of five MTRDHH2 heterozygous carriers and 95 non-carriers from various ovine breeds. We identified a single base pair duplication within the SLC33A1 gene, leading to a frameshift mutation and a premature stop codon (p.Arg246Alafs*3). SLC33A1 encodes a transmembrane transporter of acetyl-coenzyme A that is crucial for cellular metabolism. To investigate the lethality of this mutation in homozygous MTR sheep, we performed at-risk matings using artificial insemination (AI) between heterozygous SLC33A1 variant carriers (SLC33A1_dupG). Pregnancy was confirmed 15 days post-AI using a blood test measuring interferon Tau-stimulated MX1 gene expression. Ultrasonography between 45 and 60 days post-AI revealed a 12% reduction in AI success compared with safe matings, indicating embryonic/fetal loss. This was supported by the MX1 differential expression test suggesting fetal losses between 15 and 60 days of gestation. We also observed a 34.7% pre-weaning mortality rate in 49 lambs born from at-risk matings. Homozygous SLC33A1_dupG lambs accounted for 47% of this mortality, with deaths occurring mostly within the first 5 days without visible clinical signs. Therefore, appropriate management of SLC33A1_dupG with an allele frequency of 0.04 in the MTR selection scheme would help increase overall fertility and lamb survival.

A Proof-of-Concept Study of Sequential Treatment with the HDAC Inhibitor Vorinostat following BRAF and MEK Inhibitors in BRAFV600-Mutated Melanoma

PURPOSE

The development of resistance limits the clinical benefit of BRAF and MEK inhibitors (BRAFi/MEKi) in BRAFV600-mutated melanoma. It has been shown that short-term treatment (14 days) with vorinostat was able to initiate apoptosis of resistant tumor cells. We aimed to assess the antitumor activity of sequential treatment with vorinostat following BRAFi/MEKi in patients with BRAFV600-mutated melanoma who progressed after initial response to BRAFi/MEKi.

PATIENTS AND METHODS

Patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma were treated with vorinostat 360 mg once daily for 14 days followed by BRAFi/MEKi. The primary endpoint was an objective response rate of progressive lesions of at least 30% according to Response Evaluation Criteria in Solid Tumors 1.1. Secondary endpoints included progression-free survival, overall survival, safety, pharmacokinetics of vorinostat, and translational molecular analyses using ctDNA and tumor biopsies.

RESULTS

Of the 26 patients with progressive BRAFi/MEKi-resistant BRAFV600-mutated melanoma receiving treatment with vorinostat, 22 patients were evaluable for response. The objective response rate was 9%, with one complete response for 31.2 months and one partial response for 14.9 months. Median progression-free survival and overall survival were 1.4 and 5.4 months, respectively. Common adverse events were fatigue (23%) and nausea (19%). ctDNA analysis showed emerging secondary mutations in NRAS and MEK in eight patients at the time of BRAFi/MEKi resistance. Elimination of these mutations by vorinostat treatment was observed in three patients.

CONCLUSIONS

Intermittent treatment with vorinostat in patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma is well tolerated. Although the primary endpoint of this study was not met, durable antitumor responses were observed in a minority of patients (9%).

Genomic analysis reveals the association of KIT and MITF variants with the white spotting in swamp buffaloes

BACKGROUND

Swamp-type buffaloes with varying degrees of white spotting are found exclusively in Tana Toraja, South Sulawesi, Indonesia, where spotted buffalo bulls are highly valued in accordance with the Torajan customs. The white spotting depigmentation is caused by the absence of melanocytes. However, the genetic variants that cause this phenotype have not been fully characterized. The objective of this study was to identify the genomic regions and variants responsible for this unique coat-color pattern.

RESULTS

Genome-wide association study (GWAS) and selection signature analysis identified MITF as a key gene based on the whole-genome sequencing data of 28 solid and 39 spotted buffaloes, while KIT was also found to be involved in the development of this phenotype by a candidate gene approach. Alternative candidate mutations included, in addition to the previously reported nonsense mutation c.649 C > T (p.Arg217*) and splice donor mutation c.1179 + 2T > A in MITF, a nonsense mutation c.2028T > A (p.Tyr676*) in KIT. All these three mutations were located in the genomic regions that were highly conserved exclusively in Indonesian swamp buffaloes and they accounted largely (95%) for the manifestation of white spotting. Last but not the least, ADAMTS20 and TWIST2 may also contribute to the diversification of this coat-color pattern.

CONCLUSIONS

The alternative mutations identified in this study affect, at least partially and independently, the development of melanocytes. The presence and persistence of such mutations may be explained by significant financial and social value of spotted buffaloes used in historical Rambu Solo ceremony in Tana Toraja, Indonesia. Several de novo spontaneous mutations have therefore been favored by traditional breeding for the spotted buffaloes.

Co-occurrence of ST412 Klebsiella pneumoniae isolates with hypermucoviscous and non-mucoviscous phenotypes in a short-term hospitalized patient

Hypermucoviscosity (HMV) is a phenotype that is commonly associated with hypervirulence in Klebsiella pneumoniae. The factors that contribute to the emergence of HMV subpopulations remain unclear. In this study, eight K. pneumoniae strains were recovered from an inpatient who had been hospitalized for 20 days. Three of the isolates exhibited a non-HMV phenotype, which was concomitant with higher biofilm formation than the other five HMV isolates. All eight isolates were highly susceptible to serum killing, albeit HMV strains were remarkably more infective than non-HMV counterparts in a mouse model of infection. Whole genome sequencing (WGS) showed that the eight isolates belonged to the K57-ST412 lineage. Average nucleotide identity (FastANIb) analysis indicated that eight isolates share 99.96% to 99.99% similarity and were confirmed to be the same clone. Through comparative genomics analysis, 12 non-synonymous mutations were found among these isolates, eight of which in the non-HMV variants, including rmpA (c.285delG) and wbaP (c.1305T > A), which are assumed to be associated with the non-HMV phenotype. Mutations in manB (c.1318G > A), dmsB (c.577C > T) and tkt (c.1928C > A) occurred in HMV isolates only. RNA-Seq revealed transcripts of genes involved in energy metabolism, carbohydrate metabolism and membrane transport, including cysP, cydA, narK, tktA, pduQ, aceB, metN, and lsrA, to be significantly dysregulated in the non-HMV strains, suggesting a contribution to HMV phenotype development. This study suggests that co-occurrence of HMV and non-HMV phenotypes in the same clonal population may be mediated by mutational mechanisms as well as by certain genes involved in membrane transport and central metabolism.

IMPORTANCE

K. pneumoniae with a hypermucoviscosity (HMV) phenotype is a community-acquired pathogen that is associated with increased invasiveness and pathogenicity, and underlying diseases are the most common comorbid risk factors inducing metastatic complications. HMV was earlier attributed to the overproduction of capsular polysaccharide, and more data point to the possibility of several causes contributing to this bacterial phenotype. Here, we describe a unique event in which the same clonal population showed both HMV and non-HMV characteristics. Studies have demonstrated that this process is influenced by mutational processes and genes related to transport and central metabolism. These findings provide fresh insight into the mechanisms behind co-occurrence of HMV and non-HMV phenotypes in monoclonal populations as well as potentially being critical in developing strategies to control the further spread of HMV K. pneumoniae.

The rate of W chromosome degeneration across multiple avian neo-sex chromosomes

When sex chromosomes evolve recombination suppression, the sex-limited chromosome (Y/W) commonly degenerate by losing functional genes. The rate of Y/W degeneration is believed to slow down over time as the most essential genes are maintained by purifying selection, but supporting data are scarce especially for ZW systems. Here, we study W degeneration in Sylvioidea songbirds where multiple autosomal translocations to the sex chromosomes, and multiple recombination suppression events causing separate evolutionary strata, have occurred during the last ~ 28.1-4.5 million years (Myr). We show that the translocated regions have maintained 68.3-97.7% of their original gene content, compared to only 4.2% on the much older ancestral W chromosome. By mapping W gene losses onto a dated phylogeny, we estimate an average gene loss rate of 1.0% per Myr, with only moderate variation between four independent lineages. Consistent with previous studies, evolutionarily constrained and haploinsufficient genes were preferentially maintained on W. However, the gene loss rate did not show any consistent association with strata age or with the number of W genes at strata formation. Our study provides a unique account on the pace of W gene loss and reinforces the significance of purifying selection in maintaining essential genes on sex chromosomes.

Sex Influences Genetic Susceptibility to Depression-Like Behaviors in Chronic Unpredictable Mild Stress-Exposed Wistar Rats

Depression is one of the most common mood disorders among psychiatric diseases. It affects about 10% of the adult population. However, its etiopathogenesis remains poorly understood. Exploring the dynamics of stress-susceptibility and resilience will help in understanding the molecular and biological mechanisms underlying the etiopathogenesis of depression. This study aimed to determine the differences and/or similarities in factors responsible for susceptibility to depression-like behaviors in male and female Wistar rats subjected to chronic unpredictable mild stress (CUMS). Sixty Wistar rats (30 male and 30 female) weighing between 120 and 150 g were used for this study. The rats were divided into two sub-groups: control (10) and test (20) groups. Rats in the test groups were subjected to CUMS. Depression-like behaviors were assessed using light-dark box, sucrose preference, and tail suspension tests. Rats that showed depression-like behaviors following the behavioral tests (CUMS-susceptible group) were sacrificed, and their hippocampi were excised. Genomic deoxyribonucleic acid (gDNA) was purified from the hippocampal samples. Purified gDNA was subjected to whole genome sequencing (WGS). Base-calling of sequence reads from raw sequencing signal (FAST5) files was carried out, and variants were called from alignment BAM files. The corresponding VCF files generated from the variant calling experiment were filtered. Genes were identified, their impacts estimated, and variants annotated. Functional enrichment analysis was then carried out. Approximately 41% of the male and 49% of the female rats subjected to CUMS showed significant (p < 0.05) depression-like behaviors following assessment on behavioral tests. WGS of the hippocampal DNA revealed 289,839 single nucleotide polymorphisms variant types, 7002 insertions, and 34,459 deletions in males, and 1,570,186 single nucleotide polymorphisms variant types, 109,860 insertions, and 597,241 deletions in female Wistar rats. Three genes with high-impact variants were identified in male and 22 in female Wistar rats, respectively. In conclusion, female Wistar rats are more susceptible to depression-like behaviors after exposure to CUMS than males. They also have more gene variants (especially high-impact variants) than male Wistar rats.

A multifaceted kinase axis regulates plant organ abscission through conserved signaling mechanisms

Plants have evolved mechanisms to abscise organs as they develop or when exposed to unfavorable conditions.1 Uncontrolled abscission of petals, fruits, or leaves can impair agricultural productivity.2,3,4,5 Despite its importance for abscission progression, our understanding of the IDA signaling pathway and its regulation remains incomplete. IDA is secreted to the apoplast, where it is perceived by the receptors HAESA (HAE) and HAESA-LIKE2 (HSL2) and somatic embryogenesis receptor kinase (SERK) co-receptors.6,7,8,9 These plasma membrane receptors activate an intracellular cascade of mitogen-activated protein kinases (MAPKs) by an unknown mechanism.10,11,12 Here, we characterize brassinosteroid signaling kinases (BSKs) as regulators of floral organ abscission in Arabidopsis. BSK1 localizes to the plasma membrane of abscission zone cells, where it interacts with HAESA receptors to regulate abscission. Furthermore, we demonstrate that YODA (YDA) has a leading role among other MAPKKKs in controlling abscission downstream of the HAESA/BSK complex. This kinase axis, comprising a leucine-rich repeat receptor kinase, a BSK, and an MAPKKK, is known to regulate stomatal patterning, early embryo development, and immunity.10,13,14,15,16 How specific cellular responses are obtained despite signaling through common effectors is not well understood. We show that the identified abscission-promoting allele of BSK1 also enhances receptor signaling in other BSK-mediated pathways, suggesting conservation of signaling mechanisms. Furthermore, we provide genetic evidence supporting independence of BSK1 function from its kinase activity in several developmental processes. Together, our findings suggest that BSK1 facilitates signaling between plasma membrane receptor kinases and MAPKKKs via conserved mechanisms across multiple facets of plant development.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Introduction

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.

Methods

We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.

Results

Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Identification of a pathogenic mutation in ARPP21 in patients with amyotrophic lateral sclerosis

BACKGROUND AND OBJECTIVE

Between 5% and 10% of amyotrophic lateral sclerosis (ALS) cases have a family history of the disease, 30% of which do not have an identifiable underlying genetic cause after a comprehensive study of the known ALS-related genes. Based on a significantly increased incidence of ALS in a small geographical region from Spain, the aim of this work was to identify novel ALS-related genes in ALS cases with negative genetic testing.

METHODS

We detected an increased incidence of both sporadic and, especially, familial ALS cases in a small region from Spain compared with available demographic and epidemiological data. We performed whole genome sequencing in a group of 12 patients with ALS (5 of them familial) from this unique area. We expanded the study to include affected family members and additional cases from a wider surrounding region.

RESULTS

We identified a shared missense mutation (c.1586C>T; p.Pro529Leu) in the cyclic AMP regulated phosphoprotein 21 (ARPP21) gene that encodes an RNA-binding protein, in a total of 10 patients with ALS from 7 unrelated families. No mutations were found in other ALS-causing genes.

CONCLUSIONS

While previous studies have dismissed a causal role of ARPP21 in ALS, our results strongly support ARPP21 as a novel ALS-causing gene.

IDeRare: a lightweight and extensible open-source phenotype and exome analysis pipeline for germline rare disease diagnosis

Objectives

Diagnosing rare diseases is an arduous and challenging process in clinical settings, resulting in the late discovery of novel variants and referral loops. To help clinicians, we built IDeRare pipelines to accelerate phenotype-genotype analysis for patients with suspected rare diseases.

Materials and Methods

IDeRare pipeline is separated into phenotype and genotype parts. The phenotype utilizes our handmade Python library, while the genotype part utilizes command line (bash) and Python script to combine bioinformatics executable and Docker image.

Results

We described various implementations of IDeRare phenotype and genotype parts with real-world clinical and exome data using IDeRare, accelerating the terminology conversion process and giving insight on the diagnostic pathway based on disease linkage analysis until exome analysis and HTML-based reporting for clinicians.

Conclusion

IDeRare is freely available under the BSD-3 license, obtainable via GitHub. The portability of IDeRare pipeline could be easily implemented for semi-technical users and extensible for advanced users.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

The RHCE gene encodes the chicken blood system I

BACKGROUND

There are 13 known chicken blood systems, which were originally detected by agglutination of red blood cells by specific alloantisera. The genomic region or specific gene responsible has been identified for four of these systems (A, B, D and E). We determined the identity of the gene responsible for the chicken blood system I, using DNA from multiple birds with known chicken I blood system serology, 600K and 54K single nucleotide polymorphism (SNP) data, and lowpass sequence information.

RESULTS

The gene responsible for the chicken I blood system was identified as RHCE, which is also one of the genes responsible for the highly polymorphic human Rh blood group locus, for which maternal/fetal antigenic differences can result in fetal hemolytic anemia with fetal mortality. We identified 17 unique RHCE haplotypes in the chicken, with six haplotypes corresponding to known I system serological alleles. We also detected deletions in the RHCE gene that encompass more than 6000 bp and that are predicted to remove its last seven exons.

CONCLUSIONS

RHCE is the gene responsible for the chicken I blood system. This is the fifth chicken blood system for which the responsible gene and gene variants are known. With rapid DNA-based testing now available, the impact of I blood system variation on response against disease, general immune function, and animal production can be investigated in greater detail.

Known pathogenic gene variants and new candidates detected in sudden unexpected infant death using whole genome sequencing

The purpose of this study is to gain insights into potential genetic factors contributing to the infant's vulnerability to Sudden Unexpected Infant Death (SUID). Whole Genome Sequencing (WGS) was performed on 144 infants that succumbed to SUID, and 573 healthy adults. Variants were filtered by gnomAD allele frequencies and predictions of functional consequences. Variants of interest were identified in 88 genes, in 64.6% of our cohort. Seventy-three of these have been previously associated with SIDS/SUID/SUDP. Forty-three can be characterized as cardiac genes and are related to cardiomyopathies, arrhythmias, and other conditions. Variants in 22 genes were associated with neurologic functions. Variants were also found in 13 genes reported to be pathogenic for various systemic disorders and in two genes associated with immunological function. Variants in eight genes are implicated in the response to hypoxia and the regulation of reactive oxygen species (ROS) and have not been previously described in SIDS/SUID/SUDP. Seventy-two infants met the triple risk hypothesis criteria. Our study confirms and further expands the list of genetic variants associated with SUID. The abundance of genes associated with heart disease and the discovery of variants associated with the redox metabolism have important mechanistic implications for the pathophysiology of SUID.

FLI1 is associated with regulation of DNA methylation and megakaryocytic differentiation in FPDMM caused by a RUNX1 transactivation domain mutation

Familial platelet disorder with associated myeloid malignancies (FPDMM) is an autosomal dominant disease caused by heterozygous germline mutations in RUNX1. It is characterized by thrombocytopenia, platelet dysfunction, and a predisposition to hematological malignancies. Although FPDMM is a precursor for diseases involving abnormal DNA methylation, the DNA methylation status in FPDMM remains unknown, largely due to a lack of animal models and challenges in obtaining patient-derived samples. Here, using genome editing techniques, we established two lines of human induced pluripotent stem cells (iPSCs) with different FPDMM-mimicking heterozygous RUNX1 mutations. These iPSCs showed defective differentiation of hematopoietic progenitor cells (HPCs) and megakaryocytes (Mks), consistent with FPDMM. The FPDMM-mimicking HPCs showed DNA methylation patterns distinct from those of wild-type HPCs, with hypermethylated regions showing the enrichment of ETS transcription factor (TF) motifs. We found that the expression of FLI1, an ETS family member, was significantly downregulated in FPDMM-mimicking HPCs with a RUNX1 transactivation domain (TAD) mutation. We demonstrated that FLI1 promoted binding-site-directed DNA demethylation, and that overexpression of FLI1 restored their megakaryocytic differentiation efficiency and hypermethylation status. These findings suggest that FLI1 plays a crucial role in regulating DNA methylation and correcting defective megakaryocytic differentiation in FPDMM-mimicking HPCs with a RUNX1 TAD mutation.

Long-Read Structural and Epigenetic Profiling of a Kidney Tumor-Matched Sample with Nanopore Sequencing and Optical Genome Mapping

Carcinogenesis often involves significant alterations in the cancer genome architecture, marked by large structural and copy number variations (SVs and CNVs) that are difficult to capture with short-read sequencing. Traditionally, cytogenetic techniques are applied to detect such aberrations, but they are limited in resolution and do not cover features smaller than several hundred kilobases. Optical genome mapping and nanopore sequencing are attractive technologies that bridge this resolution gap and offer enhanced performance for cytogenetic applications. These methods profile native, individual DNA molecules, thus capturing epigenetic information. We applied both techniques to characterize a clear cell renal cell carcinoma (ccRCC) tumor's structural and copy number landscape, highlighting the relative strengths of each method in the context of variant size and average read length. Additionally, we assessed their utility for methylome and hydroxymethylome profiling, emphasizing differences in epigenetic analysis applicability.

Immune Checkpoint Blockade Delays Cancer and Extends Survival in Murine DNA Polymerase Mutator Syndromes

Mutations in polymerases Pold1 and Pole exonuclease domains in humans are associated with increased cancer incidence, elevated tumor mutation burden (TMB) and response to immune checkpoint blockade (ICB). Although ICB is approved for treatment of several cancers, not all tumors with elevated TMB respond. Here we generated Pold1 and Pole proofreading mutator mice and show that ICB treatment of mice with high TMB tumors did not improve survival as only a subset of tumors responded. Similarly, introducing the mutator alleles into mice with Kras/p53 lung cancer did not improve survival, however, passaging mutator tumor cells in vitro without immune editing caused rejection in immune-competent hosts, demonstrating the efficiency by which cells with antigenic mutations are eliminated. Finally, ICB treatment of mutator mice earlier, before observable tumors delayed cancer onset, improved survival, and selected for tumors without aneuploidy, suggesting the use of ICB in individuals at high risk for cancer prevention.

Highlights

Germline somatic and conditional Pold1 and Pole exonuclease domain mutations in mice produce a mutator phenotype. Spontaneous cancers arise in mutator mice that have genomic features comparable to human tumors with these mutations.ICB treatment of mutator mice with tumors did not improve survival as only a subset of tumors respond. Introduction of the mutator alleles into an autochthonous mouse lung cancer model also did not produce immunogenic tumors, whereas passaging mutator tumor cells in vitro caused immune rejection indicating efficient selection against antigenic mutations in vivo . Prophylactic ICB treatment delayed cancer onset, improved survival, and selected for tumors with no aneuploidy.

Critical assessment of missense variant effect predictors on disease-relevant variant data

Regular, systematic, and independent assessment of computational tools used to predict the pathogenicity of missense variants is necessary to evaluate their clinical and research utility and suggest directions for future improvement. Here, as part of the sixth edition of the Critical Assessment of Genome Interpretation (CAGI) challenge, we assess missense variant effect predictors (or variant impact predictors) on an evaluation dataset of rare missense variants from disease-relevant databases. Our assessment evaluates predictors submitted to the CAGI6 Annotate-All-Missense challenge, predictors commonly used by the clinical genetics community, and recently developed deep learning methods for variant effect prediction. To explore a variety of settings that are relevant for different clinical and research applications, we assess performance within different subsets of the evaluation data and within high-specificity and high-sensitivity regimes. We find strong performance of many predictors across multiple settings. Meta-predictors tend to outperform their constituent individual predictors; however, several individual predictors have performance similar to that of commonly used meta-predictors. The relative performance of predictors differs in high-specificity and high-sensitivity regimes, suggesting that different methods may be best suited to different use cases. We also characterize two potential sources of bias. Predictors that incorporate allele frequency as a predictive feature tend to have reduced performance when distinguishing pathogenic variants from very rare benign variants, and predictors supervised on pathogenicity labels from curated variant databases often learn label imbalances within genes. Overall, we find notable advances over the oldest and most cited missense variant effect predictors and continued improvements among the most recently developed tools, and the CAGI Annotate-All-Missense challenge (also termed the Missense Marathon) will continue to assess state-of-the-art methods as the field progresses. Together, our results help illuminate the current clinical and research utility of missense variant effect predictors and identify potential areas for future development.

A de novo mutation in CACNA1A is associated with autosomal dominant bovine familial convulsions and ataxia in Angus cattle

Bovine familial convulsions and ataxia (BFCA) is considered an autosomal dominant syndrome with incomplete penetrance. Nine Angus calves from the same herd were diagnosed with BFCA within days of birth. Necropsy revealed cerebellar and spinal cord lesions associated with the condition. Parentage testing confirmed that all affected calves had a common sire. The sire was then bred to 36 cows across two herds using artificial insemination, producing an additional 14 affected calves. The objective of this investigation was to identify hypothesized dominant genetic variation underlying the condition. Whole-genome sequencing was performed on the sire, six affected and seven unaffected paternal half-sibling calves and combined with data from 135 unrelated controls. The sire and five of the six affected calves were heterozygous for a nonsense variant (Chr7 g.12367906C>T, c.5073C>T, p.Arg1681*) in CACNA1A. The other affected calves (N = 8) were heterozygous for the variant but it was absent in the other unaffected calves (N = 7) and parents of the sire. This variant was also absent in sequence data from over 6500 other cattle obtained via public repositories and collaborator projects. The variant in CACNA1A is expressed in the cerebellum of the ataxic calves as detected in the transcriptome and was not differentially expressed compared with controls. The CACNA1A protein is part of a highly expressed cerebellar calcium voltage gated channel. The nonsense variant is proposed to cause haploinsufficiency, preventing proper transmission of neuronal signals through the channel and resulting in BFCA.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

RAS/RAF mutations and microsatellite instability status in primary colorectal cancers according to HER2 amplification

HER2 amplification-associated molecular alterations and clinicopathologic features in colorectal cancers (CRCs) have not been well established. In this study, we assessed the prevalence of HER2 amplification and microsatellite instability (MSI) status of 992 patients with primary CRC. In addition, molecular alterations of HER2 amplified and unamplified CRCs were examined and compared by next-generation sequencing. HER2 amplifications were found in 41 (4.1%) of 992 primary CRCs. HER2 amplification was identified in 1.0% of the right colonic tumors, 5.1% of the left colonic tumors, and 4.8% of the rectal tumors. Approximately 95% of HER2 amplification was observed in the left colon and rectum. Seven (87.5%) of eight metastatic tumors showed HER2 amplification. Most clinicopathologic features were unrelated to HER2 amplification except tumor size and MSI status. All 41 HER2 amplified CRCs were microsatellite stable. In a molecular analysis of frequently identified somatic mutations in CRCs, HER2 amplified CRCs showed a lower rate of KRAS mutations (24.4%) but a higher rate of TP53 mutations (83%) than unamplified CRCs. No BRAF and NRAS mutations were identified in HER2 amplified CRCs. Our study suggests that HER2 amplified CRCs are mutually exclusive of MSI and harbor less frequent KRAS/NRAS/BRAF mutations but frequent T53 mutations.

Genomic Context of SARS-CoV-2 Outbreaks in Farmed Mink in Spain during Pandemic: Unveiling Host Adaptation Mechanisms

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects various mammalian species, with farmed minks experiencing the highest number of outbreaks. In Spain, we analyzed 67 whole genome sequences and eight spike sequences from 18 outbreaks, identifying four distinct lineages: B.1, B.1.177, B.1.1.7, and AY.98.1. The potential risk of transmission to humans raises crucial questions about mutation accumulation and its impact on viral fitness. Sequencing revealed numerous not-lineage-defining mutations, suggesting a cumulative mutation process during the outbreaks. We observed that the outbreaks were predominantly associated with different groups of mutations rather than specific lineages. This clustering pattern by the outbreaks could be attributed to the rapid accumulation of mutations, particularly in the ORF1a polyprotein and in the spike protein. Notably, the mutations G37E in NSP9, a potential host marker, and S486L in NSP13 were detected. Spike protein mutations may enhance SARS-CoV-2 adaptability by influencing trimer stability and binding to mink receptors. These findings provide valuable insights into mink coronavirus genetics, highlighting both host markers and viral transmission dynamics within communities.

The contribution of an X chromosome QTL to non-Mendelian inheritance and unequal chromosomal segregation in Auanema freiburgense

Auanema freiburgense is a nematode with males, females, and selfing hermaphrodites. When XO males mate with XX females, they typically produce a low proportion of XO offspring because they eliminate nullo-X spermatids. This process ensures that most sperm carry an X chromosome, increasing the likelihood of X chromosome transmission compared to random segregation. This occurs because of an unequal distribution of essential cellular organelles during sperm formation, likely dependent on the X chromosome. Some sperm components are selectively segregated into the X chromosome's daughter cell, while others are discarded with the nullo-X daughter cell. Intriguingly, the interbreeding of 2 A. freiburgense strains results in hybrid males capable of producing viable nullo-X sperm. Consequently, when these hybrid males mate with females, they yield a high percentage of male offspring. To uncover the genetic basis of nullo-spermatid elimination and X chromosome drive, we generated a genome assembly for A. freiburgense and genotyped the intercrossed lines. This analysis identified a quantitative trait locus spanning several X chromosome genes linked to the non-Mendelian inheritance patterns observed in A. freiburgense. This finding provides valuable clues to the underlying factors involved in asymmetric organelle partitioning during male meiotic division and thus non-Mendelian transmission of the X chromosome and sex ratios.

A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

A recessive CLN3 variant is responsible for delayed-onset retinal degeneration in Hereford cattle

Thirteen American Hereford cattle were reported blind with presumed onset when ~12-mo-old. All blind cattle shared a common ancestor through both the maternal and paternal pedigrees, suggesting a recessive genetic origin. Given the pedigree relationships and novel phenotype, we characterized the ophthalmo-pathologic changes associated with blindness and identified the responsible gene variant. Ophthalmologic examinations of 5 blind cattle revealed retinal degeneration. Histologically, 2 blind cattle had loss of the retinal photoreceptor layer. Whole-genome sequencing (WGS) of 7 blind cattle and 9 unaffected relatives revealed a 1-bp frameshift deletion in ceroid lipofuscinosis neuronal 3 (CLN3; chr25 g.26043843del) for which the blind cattle were homozygous and their parents heterozygous. The identified variant in exon 16 of 17 is predicted to truncate the encoded protein (p. Pro369Argfs*8) battenin, which is involved in lysosomal function necessary for photoreceptor layer maintenance. Of 462 cattle genotyped, only blind cattle were homozygous for the deletion. A query of WGS data of > 5,800 animals further revealed that the variant was only observed in related Hereford cattle. Mutations in CLN3 are associated with human juvenile neuronal ceroid lipofuscinosis (JNCL), or Batten disease, which results in early-onset retinal degeneration and lesions similar to those observed in our cases. Our data support the frameshift variant of CLN3 as causative of blindness in these Hereford cattle, and provide additional evidence of the role of this gene in retinal lesions, possibly as a model for human non-syndromic JNCL.

Monkeypox virus genomic accordion strategies

The 2023 monkeypox (mpox) epidemic was caused by a subclade IIb descendant of a monkeypox virus (MPXV) lineage traced back to Nigeria in 1971. Person-to-person transmission appears higher than for clade I or subclade IIa MPXV, possibly caused by genomic changes in subclade IIb MPXV. Key genomic changes could occur in the genome's low-complexity regions (LCRs), which are challenging to sequence and are often dismissed as uninformative. Here, using a combination of highly sensitive techniques, we determine a high-quality MPXV genome sequence of a representative of the current epidemic with LCRs resolved at unprecedented accuracy. This reveals significant variation in short tandem repeats within LCRs. We demonstrate that LCR entropy in the MPXV genome is significantly higher than that of single-nucleotide polymorphisms (SNPs) and that LCRs are not randomly distributed. In silico analyses indicate that expression, translation, stability, or function of MPXV orthologous poxvirus genes (OPGs), including OPG153, OPG204, and OPG208, could be affected in a manner consistent with the established "genomic accordion" evolutionary strategies of orthopoxviruses. We posit that genomic studies focusing on phenotypic MPXV differences should consider LCR variability.

From buds to shoots: insights into grapevine development from the Witch's Broom bud sport

BACKGROUND

Bud sports occur spontaneously in plants when new growth exhibits a distinct phenotype from the rest of the parent plant. The Witch's Broom bud sport occurs occasionally in various grapevine (Vitis vinifera) varieties and displays a suite of developmental defects, including dwarf features and reduced fertility. While it is highly detrimental for grapevine growers, it also serves as a useful tool for studying grapevine development. We used the Witch's Broom bud sport in grapevine to understand the developmental trajectories of the bud sports, as well as the potential genetic basis. We analyzed the phenotypes of two independent cases of the Witch's Broom bud sport, in the Dakapo and Merlot varieties of grapevine, alongside wild type counterparts. To do so, we quantified various shoot traits, performed 3D X-ray Computed Tomography on dormant buds, and landmarked leaves from the samples. We also performed Illumina and Oxford Nanopore sequencing on the samples and called genetic variants using these sequencing datasets.

RESULTS

The Dakapo and Merlot cases of Witch's Broom displayed severe developmental defects, with no fruit/clusters formed and dwarf vegetative features. However, the Dakapo and Merlot cases of Witch's Broom studied were also phenotypically different from one another, with distinct differences in bud and leaf development. We identified 968-974 unique genetic mutations in our two Witch's Broom cases that are potential causal variants of the bud sports. Examining gene function and validating these genetic candidates through PCR and Sanger-sequencing revealed one strong candidate mutation in Merlot Witch's Broom impacting the gene GSVIVG01008260001.

CONCLUSIONS

The Witch's Broom bud sports in both varieties studied had dwarf phenotypes, but the two instances studied were also vastly different from one another and likely have distinct genetic bases. Future work on Witch's Broom bud sports in grapevine could provide more insight into development and the genetic pathways involved in grapevine.