DELLY: structural variant discovery by integrated paired-end and split-read analysis

Colibactin Exerts Androgen-dependent and -independent Effects on Prostate Cancer

BACKGROUND AND OBJECTIVE

The etiology of prostate cancer (PC) is multifactorial and poorly understood. It has been suggested that colibactin-producing Escherichia coli positive for the pathogenicity island pks (pks+) initiate cancers via induction of genomic instability. In PC, androgens promote oncogenic translocations. Our aim was to investigate the association of pks+E. coli with PC diagnosis and molecular architecture, and its relationship with androgens.

METHODS

We quantified the association of pks+E. coli with PC diagnosis in a volunteer-sampled 235-person cohort from two institutional practices (UT San Antonio). We then used colibactin 742 and DNA/RNA sequencing to evaluate the effects of colibactin 742, dihydrotestosterone (DHT), and their combination in vitro.

KEY FINDINGS AND LIMITATIONS

Colibactin exposure was positively associated with PC diagnosis (p = 0.04) in our clinical cohort, and significantly increased replication fork stalling and fusions in vitro (p < 0.01). Combined in vitro exposure to colibactin 742 and DHT induced more somatic mutations of all types than exposure to either alone. The combination also elicited kataegis, with a higher density of somatic point mutations. Laboratory analyses were conducted using a single cell line, which limited our ability to fully recapitulate the complexity of PC etiology.

CONCLUSIONS AND CLINICAL IMPLICATIONS

Our findings are consistent with synergistic induction of genome instability and kataegis by colibactin 742 and DHT in cell culture. Colibactin-producing pks+ E. coli may plausibly contribute to PC etiology.

PATIENT SUMMARY

We investigated whether a bacterial toxin that is linked to colon cancer can also cause prostate cancer. Our results support this idea by showing a link between the toxin and prostate cancer diagnosis in a large patient population. We also found that this toxin causes genetic dysfunction in prostate cancer cells when combined with testosterone.

Transposition of a non-autonomous element into the Gβ gene of Schizophyllum commune causes the streak mutation

Streak mutants of Schizophyllum commune are characterized by ropy, hyperbranching hyphae, suppressed aerial hyphae formation, and the production of pigments. Additionally, these mutants dikaryotize unilaterally, with the mutant fertilizing its compatible mating partner, but not accepting its nucleus. Here we show that a 512 bp non-autonomous transposable element had integrated in the Gβ protein of a streak mutant of S. commune. This element has the same 50 bp inverted repeat as an autonomous element, dubbed Bike transposon. Its transposase has homologues in various Agaricomycetes. Introducing the Gβ gene in the streak mutant restored the wild-type phenotype showing that the integration of the 512 bp element in the Gβ gene is responsible for the streak phenotype.

Deciphering the complex molecular architecture of the genetically modified soybean FG72 through paired-end whole genome sequencing

The clear molecular characterization of genetically modified (GM) plants and animals is a prerequisite for obtaining regulatory approval and safety certification for commercial cultivation. This characterization includes the identification of the transferred DNA (T-DNA) insertion site, its flanking sequences, the copy number of inserted genes, and the detection of any unintended genomic alterations accompanying the transformation process. In this study, we performed a comprehensive molecular characterization of the well-known GM soybean event FG72 using paired-end whole-genome sequencing (PE-WGS). We examined the T-DNA insertion site, flanking sequences, the entire structure and copy number of the T-DNA integration, the presence of plasmid backbone sequences, and genome-wide structural variations (SVs). Our analysis revealed that the T-DNA integrated into two distinct sites on chromosome 15 of the host genome, accompanied by a translocation of host genomic sequences. One site harbored a partial T-DNA integration, while the other site contained two tandem repeats of the full T-DNA. Importantly, no plasmid backbone sequences were detected in the host genome, indicating a clean T-DNA integration during the biolistic transformation process. Furthermore, we identified numerous genome-wide SVs, with chromosome 15 ranking second among all 20 chromosomes in terms of SV frequency, and most of these variations occurring within gene-coding regions. These results provide a refined and comprehensive molecular characterization of the FG72 soybean event, which could further support its commercial approval and cultivation. Our work highlights the utility of the PE-WGS approach as a sensitive and labor-efficient alternative to conventional molecular characterization techniques, generating comprehensive data to facilitate the safety assessment of GM crops during research and commercialization pipelines.

Identification of Genomic Structural Variations in Xinjiang Brown Cattle by Deep Sequencing and Their Association with Body Conformation Traits

Xinjiang Brown cattle is an elite dual-purpose breed (raised for dairy and beef) developed in China. To elucidate its genomic architecture, we conducted whole-genome resequencing of 169 Xinjiang Brown cattle, followed by structural variation (SV) detection and a genome-wide association study (GWAS). We identified 71,668 SVs, among which deletions were the most prevalent, followed by translocations, inversions, duplications, and insertions. We further identified 1286 high-frequency SVs involving 2016 protein-coding genes. Through functional enrichment analysis of these genes, we revealed associations of genetic variation at genomic positions near genes implicated in immune response and disease resistance (NFKBIZ and PTPRT), growth and development (HDAC4 and MEF2A), and milk production (TP63, FABP4, and MEF2A). GWAS analysis of 31 body conformation traits revealed 58 SVs significantly associated with five traits (chest width, rear udder width, udder depth, rump width, and heel depth) at the genome-wide level. Additionally, nine candidate genes (CLINT1, EBF1, PAM16, GRIP1, CFAP54, SLC22A16, DOK5, ETAA1, and IPMK) were identified as potentially involved in the genetic regulation of body conformation traits. These findings provide novel insights for genetic improvement strategies and indicate that precision breeding could further enhance the production performance of this breed in the future.

Whole-genome sequencing analyses suggest novel genetic factors associated with Alzheimer's disease and a cumulative effects model for risk liability

Genome-wide association studies (GWAS) on Alzheimer's disease (AD) have predominantly focused on identifying common variants in Europeans. Here, we performed whole-genome sequencing (WGS) of 1,559 individuals from a Korean AD cohort to identify various genetic variants and biomarkers associated with AD. Our GWAS analysis identified a previously unreported locus for common variants (APCDD1) associated with AD. Our WGS analysis was extended to explore the less-characterized genetic factors contributing to AD risk. We identified rare noncoding variants located in cis-regulatory elements specific to excitatory neurons associated with cognitive impairment. Moreover, structural variation analysis showed that short tandem repeat expansion was associated with an increased risk of AD, and copy number variant at the HPSE2 locus showed borderline statistical significance. APOE ε4 carriers with high polygenic burden or structural variants exhibited severe cognitive impairment and increased amyloid beta levels, suggesting a cumulative effects model of AD risk.

Clinical Implication of Sequential Circulating Tumor DNA Assessments for the Treatment of Diffuse Large B-Cell Lymphoma

BACKGROUND/OBJECTIVES

Circulating tumor DNA (ctDNA) has emerged as a promising biomarker for non-invasive tumor monitoring in diffuse large B-cell lymphoma (DLBCL).

METHODS

In this study, 52 patients with newly diagnosed advanced-stage DLBCL treated with R-CHOP underwent serial ctDNA analysis at baseline, interim (after three cycles), and end of treatment. The prognostic significance of ctDNA dynamics was evaluated, and its predictive value was compared with the PET/CT response.

RESULTS

Targeted next-generation sequencing revealed baseline ctDNA in 98.1% of patients, with 74.7% concordance to tumor tissue genotyping. Higher baseline ctDNA levels correlated with elevated LDH, older age, and high IPI scores. A ≥2-log reduction in ctDNA at interim was significantly associated with improved overall survival (p = 0.004), though not with progression-free survival. Notably, combining interim ctDNA dynamics with PET/CT results enhanced the predictive accuracy for treatment outcomes, particularly among patients with partial metabolic responses.

CONCLUSIONS

These findings support the clinical utility of ctDNA for dynamic risk assessment in DLBCL, and suggest that integrating ctDNA with imaging biomarkers may guide more personalized therapeutic strategies. Further validation using highly sensitive ctDNA assays is warranted to optimize its role in routine clinical practice.

A bovine model of rhizomelic chondrodysplasia punctata caused by a deep intronic splicing variant in the GNPAT gene

BACKGROUND

Genetic defects that occur naturally in livestock species provide valuable models for investigating the molecular mechanisms underlying rare human diseases. Livestock breeds are subject to the regular emergence of recessive genetic defects due to genetic drift and recent inbreeding. At the same time, their large population sizes provide easy access to case and control individuals and to massive amounts of pedigree, genomic and phenotypic information recorded for management and selection purposes. In this study, we investigated a lethal form of recessive chondrodysplasia observed in 21 stillborn calves of the Aubrac beef cattle breed.

RESULTS

Detailed examinations of three affected calves revealed proximal limb shortening, epiphyseal calcific deposits, and other pathological signs consistent with human rhizomelic chondrodysplasia punctata, a rare peroxisomal disorder caused by recessive variants in one of five genes (AGPS, FAR1, GNPAT, PEX5, and PEX7). Using homozygosity mapping, whole genome sequencing of two affected individuals, and filtering for variants found in 1867 control genomes, we reduced the list of candidate variants to a single deep intronic substitution in GNPAT (NC_037355.1:g.4039268G > A on chromosome 28 of the ARS-UCD1.2 bovine genome assembly). For verification, we performed large-scale genotyping of this variant using a custom SNP array and found a perfect genotype-phenotype correlation in 21 cases and 26 of their parents, and a complete absence of homozygotes in 1195 unaffected Aubrac controls. The g.4039268A allele segregated at a frequency of 2.6% in this population and was absent in 375,535 additional individuals from 17 breeds. Then, using in vivo and in vitro analyses, we demonstrated that the derived allele activates cryptic splice sites within intron 11 resulting in abnormal transcripts. Finally, by mining the wealth of records available in the French bovine database, we also reported suggestive effects on juvenile mortality (and not just stillbirth) in homozygotes and on muscle development in heterozygotes, which merit further investigation.

CONCLUSIONS

We report the first spontaneous large animal model of rhizomelic chondrodysplasia punctata and provide a diagnostic test to select against this defect in cattle. Our work also brings interesting insights into the molecular consequences of complete or partial GNPAT insufficiency in mammals.

Understanding the epidemiology and pathogenesis of Mycobacterium tuberculosis with non-redundant pangenome of epidemic strains in China

Tuberculosis is a major public health threat resulting in more than one million lives lost every year. Many challenges exist to defeat this deadly infectious disease which address the importance of a thorough understanding of the biology of the causative agent Mycobacterium tuberculosis (MTB). We generated a non-redundant pangenome of 420 epidemic MTB strains from China including 344 Lineage 2 strains, 69 Lineage 4 strains, six Lineage 3 strains, and one Lineage 1 strain. We estimate that MTB strains have a pangenome of 4,278 genes encoding 4,183 proteins, of which 3,438 are core genes. However, due to 99,694 interruptions in 2,447 coding genes, we can only confidently confirm 1,651 of these genes are translated in all samples. Of these interruptions, 67,315 (67.52%) could be classified by various genetic variations detected by currently available tools, and more than half of them are due to structural variations, mostly small indels. Assuming a proportion of these interruptions are artifacts, the number of active core genes would still be much lower than 3,438. We further described differential evolutionary patterns of genes under the influences of selective pressure, population structure and purifying selection. While selective pressure is ubiquitous among these coding genes, evolutionary adaptations are concentrated in 1,310 genes. Genes involved in cell wall biogenesis are under the strongest selective pressure, while the biological process of disruption of host organelles indicates the direction of the most intensive positive selection. This study provides a comprehensive view on the genetic diversity and evolutionary patterns of coding genes in MTB which may deepen our understanding of its epidemiology and pathogenicity.

New insights into the cold tolerance of upland switchgrass by integrating a haplotype-resolved genome and multi-omics analysis

BACKGROUND

Switchgrass (Panicum virgatum L.) is a bioenergy and forage crop. Upland switchgrass exhibits superior cold tolerance compared to the lowland ecotype, but the underlying molecular mechanisms remain unclear.

RESULTS

Here, we present a high-quality haplotype-resolved genome of the upland ecotype "Jingji31." We then conduct multi-omics analysis to explore the mechanism underlying its cold tolerance. By comparative transcriptome analysis of the upland and lowland ecotypes, we identify many genes with ecotype-specific differential expression, particularly members of the cold-responsive (COR) gene family, under cold stress. Notably, AFB1, ATL80, HOS10, and STRS2 gene families show opposite expression changes between the two ecotypes. Based on the haplotype-resolved genome of "Jingji31," we detect more cold-induced allele-specific expression genes in the upland ecotype than in the lowland ecotype, and these genes are significantly enriched in the COR gene family. By genome-wide association study, we detect an association signal related to the overwintering rate, which overlaps with a selective sweep region containing a cytochrome P450 gene highly expressed under cold stress. Heterologous overexpression of this gene in rice alleviates leaf chlorosis and wilting under cold stress. We also verify that expression of this gene is suppressed by a structural variation in the promoter region.

CONCLUSIONS

Based on the high-quality haplotype-resolved genome and multi-omics analysis of upland switchgrass, we characterize candidate genes responsible for cold tolerance. This study advances our understanding of plant cold tolerance, which provides crop breeding for improved cold tolerance.

Genome wide population genetics and molecular surveillance of insecticide resistance in Anopheles stephensi mosquitoes from Awash Sebat Kilo in Ethiopia

Since the detection of the Asian mosquito Anopheles stephensi in Dijbouti in 2012, it has spread throughout the Horn of Africa. This invasive vector continues to expand across the continent and is a significant threat to malaria control programs. Vector control methods, including insecticide-treated nets and indoor residual spraying, have substantially reduced the malaria burden. However, the increasing prevalence of mosquitoes resistant to insecticides, including An. stephensi populations, undermines ongoing malaria elimination efforts. Understanding population structure, gene flow between populations, and the distribution of insecticide resistance mutations is essential for guiding effective malaria control strategies. Here, we generated whole genome sequencing data for An. stephensi sourced from Awash Sebat Kilo, Ethiopia (n = 27) and compared with South Asian populations (n = 45; India and Pakistan) to assess genomic diversity, population structure, and uncovering insecticide resistance mutations. Population structure analysis using genome-wide single nucleotide polymorphisms (n = 15,533,476) revealed Ethiopian isolates clustering as a distinct ancestral group, separate from South Asian isolates. Three insecticide resistance-associated SNPs (gaba gene: A296S and V327I; vgsc L1014F) were detected. Evidence of ongoing selection was found in several loci, including genes previously associated with neonicotinoids, ivermectin, DDT, and pyrethroid resistance. This study represents the first whole genome population genetics study of invasive An. stephensi, revealing genomic differences from South Asian populations, which can be used for future assessments of vector population dispersal and detection of insecticide resistance mechanisms.

Mapping the genetic landscape establishing a tumor immune microenvironment favorable for anti-PD-1 response

Identifying host genetic factors modulating immune checkpoint inhibitor (ICI) efficacy is experimentally challenging. Our approach, utilizing the Collaborative Cross mouse genetic resource, fixes the tumor genomic configuration while varying host genetics. We find that response to anti-PD-1 (aPD1) immunotherapy is significantly heritable in four distinct murine tumor models (H2: 0.18-0.40). For the MC38 colorectal carcinoma system, we map four significant ICI response quantitative trait loci (QTLs) with significant epistatic interactions. The differentially expressed genes within these QTLs that define responder genetics are highly enriched for processes involving antigen processing and presentation, allograft rejection, and graft vs. host disease (all p < 1 × 10-10). Functional blockade of two top candidate immune targets, GM-CSF and IL-2RB, completely abrogates the MC38 transcriptional response to aPD1 therapy. Thus, our in vivo experimental platform is a powerful approach for discovery of host genetic factors that establish the tumor immune microenvironment propitious for ICI response.

Genomic landscape of diffuse glioma revealed by whole genome sequencing

Diffuse gliomas are the commonest malignant primary brain tumour in adults. Herein, we present analysis of the genomic landscape of adult glioma, by whole genome sequencing of 403 tumours (256 glioblastoma, 89 astrocytoma, 58 oligodendroglioma; 338 primary, 65 recurrence). We identify an extended catalogue of recurrent coding and non-coding genetic mutations that represents a source for future studies and provides a high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and extrachromosomal DNA. Finally, we relate these to clinical outcome. As well as identifying drug targets for treatment of glioma our findings offer the prospect of improving treatment allocation with established targeted therapies.

A Systematic Review of the Advances and New Insights into Copy Number Variations in Plant Genomes

Copy number variations (CNVs), as an important structural variant in genomes, are widely present in plants, affecting their phenotype and adaptability. In recent years, CNV research has not only focused on changes in gene copy numbers but has also been linked to complex mechanisms such as genome rearrangements, transposon activity, and environmental adaptation. The advancement in sequencing technologies has made the detection and analysis of CNVs more efficient, not only revealing their crucial roles in plant disease resistance, adaptability, and growth development, but also demonstrating broad application potential in crop improvement, particularly in selective breeding and genomic selection. By studying CNV changes during the domestication process, researchers have gradually recognized the important role of CNVs in plant domestication and evolution. This article reviews the formation mechanisms of CNVs in plants, methods for their detection, their relationship with plant traits, and their applications in crop improvement. It emphasizes future research directions involving the integration of multi-omics to provide new perspectives on the structure and function of plant genomes.

RPE65 variant p.(E519K) causes a novel dominant adult-onset maculopathy in 83 affected individuals

Recessive RPE65-related retinopathy is an inherited retinal disease (IRD) that is a well-established target for gene therapy. Dominant RPE65-related retinopathy, however, due to Irish founder variant p.(D477G), is extremely rare. Here, we report the discovery, replication and characterization of a novel dominant retinopathy caused by RPE65 variant p.(E519K), identified in 83 individuals of European ancestry across IRD registries (Belgian discovery cohort, n=2,873; replication cohort, n=18,796). Long-read sequencing-based haplotyping revealed a shared region of 464 kb, supporting a founder effect. Genotype-phenotype data support dominant inheritance and phenotypic variability respectively, characterized by late-onset macular dystrophy with two main subtypes, a pathognomonic mottled subtype and a pattern dystrophy subtype. Functional studies showed that the p.(E519K) variant affects RPE65 enzymatic activity, correlating with lower protein expression. Protein modelling and cellular thermal shift assays further supported a destabilizing effect on protein structure. Overall, our work provides strong genetic, clinical, molecular and functional evidence for a novel dominant RPE65 retinopathy in multiple families in Europe and North America due to a Belgian founder variant. This discovery reduces the diagnostic gap in dominant IRD, particularly in individuals of European ancestry. Finally, it lays the foundation for developing therapeutic strategies targeting dominant RPE65 retinopathy.

The Germline and Somatic Origins of Prostate Cancer Heterogeneity

SIGNIFICANCE

This study uncovered 223 recurrently mutated driver regions using the largest cohort of prostate tumors to date. It reveals associations between germline SNPs, somatic drivers, and tumor aggression, offering significant insights into how prostate tumor evolution is shaped by germline factors and the timing of somatic mutations.

The contribution of de novo coding mutations to meningomyelocele

Meningomyelocele (also known as spina bifida) is considered to be a genetically complex disease resulting from a failure of the neural tube to close. Individuals with meningomyelocele display neuromotor disability and frequent hydrocephalus, requiring ventricular shunting. A few genes have been proposed to contribute to disease susceptibility, but beyond that it remains unexplained1. We postulated that de novo mutations under purifying selection contribute to the risk of developing meningomyelocele2. Here we recruited a cohort of 851 meningomyelocele trios who required shunting at birth and 732 control trios, and found that de novo likely gene disruption or damaging missense mutations occurred in approximately 22.3% of subjects, with 28% of such variants estimated to contribute to disease risk. The 187 genes with damaging de novo mutations collectively define networks including actin cytoskeleton and microtubule-based processes, Netrin-1 signalling and chromatin-modifying enzymes. Gene validation demonstrated partial or complete loss of function, impaired signalling and defective closure of the neural tube in Xenopus embryos. Our results indicate that de novo mutations make key contributions to meningomyelocele risk, and highlight critical pathways required for neural tube closure in human embryogenesis.

What Are We Learning from Plant Pangenomes?

A single reference genome does not fully capture species diversity. By contrast, a pangenome incorporates multiple genomes to capture the entire set of nonredundant genes in a given species, along with its genome diversity. New sequencing technologies enable researchers to produce multiple high-quality genome sequences and catalog diverse genetic variations with better precision. Pangenomic studies have detected structural variants in plant genomes, dissected the genetic architecture of agronomic traits, and helped unravel molecular underpinnings and evolutionary origins of plant phenotypes. The pangenome concept has further evolved into a so-called super-pangenome that includes wild relatives within a genus or clade and shifted to graph-based reference systems. Nevertheless, building pangenomes and representing complex structural variants remain challenging in many crops. Standardized computing pipelines and common data structures are needed to compare and interpret pangenomes. The growing body of plant pangenomics data requires new algorithms, huge data storage capacity, and training to help researchers and breeders take advantage of newly discovered genes and genetic variants.

Characterization of Extrachromosomal Circular DNA in Primary and Cisplatin-Resistant High-Grade Serous Ovarian Cancer

BACKGROUND

Cisplatin resistance is a major cause of tumor recurrence and mortality in high-grade serous ovarian cancer (HGSOC). Extrachromosomal circular DNA (eccDNA) has emerged as a critical factor in tumor evolution and drug resistance. However, the specific contribution of eccDNA to cisplatin resistance in HGSOC remains unclear.

METHODS

We performed whole-genome sequencing, Circle-Seq, and RNA-Seq in four pairs of primary and cisplatin-resistant (cisR) HGSOC cell lines to characterize genome-wide eccDNA distribution and features. Functional enrichment analyses were subsequently conducted on differentially expressed eccDNA-related genes.

RESULTS

In the SKOV3 cisR cell line, we identified a large extrachromosomal circular DNA (ecDNA) carrying the HIF1A gene, which regulates DNA repair, drug efflux, and epithelial-mesenchymal transition, contributing to cisplatin resistance. Using Circle-Seq, we detected a total of 161,062 eccDNAs, most of which were less than 1000 bp and distributed across all chromosomes. Notably, the number of eccDNAs on chromosome 21 differed significantly between the primary and cisR cell lines. Additionally, eccDNAs were predominantly located in non-coding repetitive elements. Functional analysis of eccDNA-related differentially expressed genes revealed that, compared to primary cell lines, cisR cell lines were associated with mitotic spindle assembly, regulation of vascular permeability, and cell differentiation. eccDNA-related genes involved in these pathways include MISP, WIPF1, RHOD, KRT80, and PLVAP.

CONCLUSIONS

Our findings suggest that eccDNAs, particularly ecDNA amplifications like HIF1A, contribute significantly to cisplatin resistance mechanisms in HGSOC. These insights highlight eccDNA as a potential target for overcoming therapeutic resistance and improving treatment outcomes in ovarian cancer.

CanSeer: a translational methodology for developing personalized cancer models and therapeutics

Computational modeling and analysis of biomolecular network models annotated with omics data are emerging as a versatile tool for designing personalized therapies. Current endeavors aimed at employing in silico models towards personalized cancer therapeutics remain limited in providing all-in-one approach that ascertains actionable targets, re-positions FDA (Food and Drug Administration) approved drugs, furnishes quantitative cues on therapy responses such as efficacy and cytotoxic effect, and identifies novel drug combinations. Here we propose "CanSeer"-a methodology for developing personalized therapeutics. CanSeer employs patient-specific genetic alterations and RNA-seq data to annotate in silico models followed by dynamical network analyses towards assessment of treatment responses. To exemplify, three use cases involving paired samples, unpaired samples, and cancer samples only, of lung squamous cell carcinoma (LUSC) patients are provided. CanSeer reveals the effectiveness of repositioned drugs along with the identification of several novel LUSC treatment combinations including Afuresertib + Palbociclib, Dinaciclib + Trametinib, Afatinib + Oxaliplatin, Ulixertinib + Olaparib, etc.

Genetic Etiology of Epilepsy: A Retrospective Study From a Single-Center Cohort

Next generation sequencing (NGS) technology has made significant progress in the genetic diagnosis and treatment of epilepsy. However, genetic studies on epilepsy with different etiologies remain relatively limited. In this study, whole-genome or whole-exome sequencing was performed on 158 unrelated patients with epilepsy of various etiologies, and the identified variants were analyzed for their association with 1356 seizure-related genes in the database. Additionally, the pathogenicity or likely pathogenicity of those variants associated with known epilepsy genes was evaluated. The results showed that pathogenic or likely pathogenic variants were detected in 31.65% (50/158) of the patients in our cohort study. Further analysis revealed significant differences in the diagnostic rates among different epilepsy categories: 29.60% (37/125) for idiopathic epilepsy and 39.39% (13/33) for symptomatic epilepsy. Moreover, the genes PRRT2, KMT2C, PRKRA, NOTCH3, NAGLU, and SCN1A were identified as potentially important for epilepsy, suggesting they could become key targets for clinical diagnosis and treatment. In conclusion, NGS technology demonstrates high diagnostic efficiency for epilepsy of different etiologies and highlights significant differences among various types. This provides novel genetic insights for the diagnosis and treatment of epilepsy.

Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline

Summary

The price, quality and throughout of DNA sequencing continue to improve. Algorithmic innovations have allowed inference of a growing range of features from DNA sequencing data, quantifying nuclear, mitochondrial and evolutionary aspects of both germline and somatic genomes. To automate analyses of the full range of genomic characteristics, we created an extensible Nextflow meta-pipeline called metapipeline-DNA. Metapipeline-DNA analyzes targeted and whole-genome sequencing data from raw reads through pre-processing, feature detection by multiple algorithms, quality-control and data- visualization. Each step can be run independently and is supported robust software engineering including automated failure-recovery, robust testing and consistent verifications of inputs, outputs and parameters. Metapipeline-DNA is cloud-compatible and highly configurable, with options to subset and optimize each analysis. Metapipeline-DNA facilitates high-scale, comprehensive analysis of DNA sequencing data.

Availability

Metapipeline-DNA is an open-source Nextflow pipeline under the GPLv2 license and is available at https://github.com/uclahs-cds/metapipeline-DNA .

DeBasher: a flow-based programming bash extension for the implementation of complex and interactive workflows with stateful processes

BACKGROUND

Bioinformatics data analysis faces significant challenges. As data analysis often takes the form of pipelines or workflows, workflow managers (WfMs) have become essential. Data flow programming constitutes the preferred approach in WfMs, enabling parallel processes activated reactively based on input availability. While this paradigm typically follows a linear, acyclic progression, cyclic workflows are sometimes necessary in bioinformatics analyses. These cyclic workflows also present an opportunity to explore workflow interactivity, a feature not widely implemented in existing WfMs.

RESULTS

We propose DeBasher, a tool that adopts the flow-based programming (FBP) paradigm, in which the workflow components are in control of their life cycle and can store state information, allowing the execution of complex workflows that include cycles. DeBasher also incorporates a powerful model of interactivity, where the user can alter the behavior of a running workflow. Additionally, DeBasher allows the user to define triggers so as to initiate the execution of a complete workflow or a part of it. The ability to execute processes with state and in control of their life cycle also has applications in dynamic scheduling tasks. Furthermore, DeBasher presents a series of extra features, including the combination of multiple workflows at runtime through a feature we have called runtime piping, switching to static scheduling to increase scalability, or implementing processes in multiple languages. DeBasher has been successfully used to process 131.7 TB of genomic data by means of a variant calling pipeline.

CONCLUSIONS

DeBasher is an FBP Bash extension that can be useful in a wide range of situations and in particular when implementing complex workflows, workflows with interactivity or triggers, or when a high scalability is required.

A Hitchhiker's Guide to long-read genomic analysis

Over the past decade, long-read sequencing has evolved into a pivotal technology for uncovering the hidden and complex regions of the genome. Significant cost efficiency, scalability, and accuracy advancements have driven this evolution. Concurrently, novel analytical methods have emerged to harness the full potential of long reads. These advancements have enabled milestones such as the first fully completed human genome, enhanced identification and understanding of complex genomic variants, and deeper insights into the interplay between epigenetics and genomic variation. This mini-review provides a comprehensive overview of the latest developments in long-read DNA sequencing analysis, encompassing reference-based and de novo assembly approaches. We explore the entire workflow, from initial data processing to variant calling and annotation, focusing on how these methods improve our ability to interpret a wide array of genomic variants. Additionally, we discuss the current challenges, limitations, and future directions in the field, offering a detailed examination of the state-of-the-art bioinformatics methods for long-read sequencing.

Candida albicans isolates contain frequent heterozygous structural variants and transposable elements within genes and centromeres

The human fungal pathogen Candida albicans poses a significant burden on global health, causing high rates of mortality and antifungal drug resistance. C. albicans is a heterozygous diploid organism that reproduces asexually. Structural variants (SVs) are an important source of genomic rearrangement, particularly in species that lack sexual recombination. To comprehensively investigate SVs across clinical isolates of C. albicans, we conducted long-read sequencing and genome-wide SV analysis in three distantly related clinical isolates. Our work includes a new, comprehensive analysis of transposable element (TE) composition, location, and diversity. SVs and TEs are frequently close to coding sequences and many SVs are heterozygous, suggesting that SVs might impact gene and allele-specific expression. Most SVs are uniquely present in only one clinical isolate, indicating that SVs represent a significant source of intraspecies genetic variation. We identify multiple, distinct SVs at the centromeres of Chromosome 4 and Chromosome 5, including inversions and transposon polymorphisms. These two chromosomes are often aneuploid in drug-resistant clinical isolates and can form isochromosome structures with breakpoints near the centromere. Further screening of 100 clinical isolates confirms the widespread presence of centromeric SVs in C. albicans, often appearing in a heterozygous state, indicating that SVs are contributing to centromere evolution in C. albicans Together, these findings highlight that SVs and TEs are common across diverse clinical isolates of C. albicans and that the centromeres of this organism are important sites of genome rearrangement.

Benchmarking of germline copy number variant callers from whole genome sequencing data for clinical applications

Motivation

Whole-genome sequencing (WGS) is increasingly preferred for clinical applications due to its comprehensive coverage, effectiveness in detecting copy number variants (CNVs), and declining costs. However, systematic evaluations of WGS CNV callers tailored to germline clinical testing-where high sensitivity and confirmation of reported CNVs are essential-remain necessary. Clinical reporting typically emphasizes CNVs affecting coding regions over precise breakpoint detection. This study benchmarks several short-read WGS CNV detection tools using reference cell lines to inform their clinical use.

Results

While tools vary in sensitivity (7%-83%) and precision (1%-76%), few meet the sensitivity needed for clinical testing. Callers generally perform better for deletions (up to 88% sensitivity) than duplications (up to 47% sensitivity), with poor detection of duplications under 5 kb. Notably, for CNVs in genes commonly included in clinical panels, significantly improved sensitivity and precision were observed when benchmarking against 25 cell lines with known CNVs. DRAGEN v4.2 high-sensitivity CNV calls, post-processed with custom filters, achieved 100% sensitivity and 77% precision on the optimized gene panel after excluding recurring artifacts. This level of performance may support clinical use with orthogonal confirmation of reportable CNVs, pending validation on laboratory-specific samples.

Availability and implementation

The data underlying this article are available in the European Nucleo-tide Archive under project accession PRJEB87628.

The continued promise of genomic technologies and software in neurogenetics

The continued evolution of genomic technologies over the past few decades has revolutionized the field of neurogenetics, offering profound insights into the genetic underpinnings of neurological disorders. Identification of causal genes for numerous monogenic neurological conditions has informed key aspects of disease mechanisms and facilitated research into critical proteins and molecular pathways, laying the groundwork for therapeutic interventions. However, the question remains: has this transformative trend reached its zenith? In this review, we suggest that despite significant strides in genome sequencing and advanced computational analyses, there is still ample room for methodological refinement. We anticipate further major genetic breakthroughs corresponding with the increased use of long-read genomes, variant calling software, AI tools, and data aggregation databases. Genetic progress has historically been driven by technological advancements from the commercial sector, which are developed in response to academic research needs, creating a continuous cycle of innovation and discovery. This review explores the potential of genomic technologies to address the challenges of neurogenetic disorders. By outlining both established and modern resources, we aim to emphasize the importance of genetic technologies as we enter an era poised for discoveries.

FGFR2 fusion/rearrangement analysis in intrahepatic cholangiocarcinoma using DNA/RNA-based NGS and FISH

Patients with intrahepatic cholangiocarcinoma (iCCA) harboring FGFR2 fusion/rearrangement benefit from targeted therapies, highlighting the need for reliable testing strategies to identify FGFR2 alterations. We assessed 226 iCCA cases using RNA-based NGS, DNA-based NGS, and break-apart FISH to evaluate the effectiveness of these methods in detecting FGFR2 fusion/rearrangement. The detection rates for FGFR2 fusion/rearrangement were 9.7% (22/226) for RNA-based NGS, 7.1% (16/226) for DNA-based NGS, and 10.2% (23/226) for FISH. Among the 26 FGFR2 fusion/rearrangement-positive cases identified by any method, only 15 (57.7%) were positive by all three techniques, yielding a concordance rate of 95.1% (215/226). RNA-based NGS confirmed oncogenic FGFR2 fusion in 81% (21/26) of positive cases and identified five novel oncogenic fusions. Thirty-five percent (6/17) of the partner genes were located on chromosome 10, with BICC1 being the most common fusion partner, while the rest were distributed across the other 9 chromosomes. FISH demonstrated a sensitivity of 95.2% and specificity of 98.5%, compared to oncogenic FGFR2 fusions confirmed by RNA-based NGS, while DNA-based NGS exhibited a sensitivity of 71.4% and specificity of 99.5%, identifying FGFR2 mutations in 4 cases. FGFR2-FISH positive cases displayed no significant heterogeneity in positive cell distribution. Oncogenic FGFR2 fusion/rearrangement was associated with small duct type iCCA, especially in cases with positive serum HBsAg and absent cholangiolocarcinoma components and peripheral liver steatosis. This study provides a comprehensive comparison of three assays for detecting FGFR2 fusion/rearrangement, along with clinicopathologic characterization of oncogenic FGFR2 fusion in iCCA.

Mars: simplifying bioinformatics workflows through a containerized approach to tool integration and management

Summary

Bioinformatics is a rapidly evolving field with numerous specialized tools developed for essential genomic analysis tasks, such as read simulation, mapping, and variant calling. However, managing these tools presents significant challenges due to varied dependencies, execution steps, and output formats, complicating the installation and configuration processes. To address these issues, we introduce "Mars" a bioinformatics solution encapsulated within a singularity container that preloads a comprehensive suite of widely used genomic tools. Mars not only simplifies the installation of these tools but also automates critical workflow functions, including sequence sample preparation, read simulation, read mapping, variant calling, and result comparison. By streamlining the execution of these workflows, Mars enables users to easily manage input-output formats and compare results across different tools, thereby enhancing reproducibility and efficiency. Furthermore, by providing a cohesive environment that integrates tool management with a flexible workflow interface, Mars empowers researchers to focus on their analyses rather than the complexities of tool configuration. This integrated solution facilitates the testing of various combinations of tools and algorithms, enabling users to evaluate performance based on different metrics and identify the optimal tools for their specific genomic analysis needs. Through Mars, we aim to enhance the accessibility and usability of bioinformatics tools, ultimately advancing research in genomic analysis.

Availability and implementation

Mars is freely available at https://github.com/GenomicAI/mars. It is implemented within a Singularity container environment and supports modular extension for additional genomic tools and custom workflows.

The tandem duplicator phenotype may be a novel targetable subgroup in pancreatic cancer

Tandem duplicator phenotype (TDP) consists of distinct genomic rearrangements where tandem duplications are randomly distributed. In this study, we characterized the prevalence and outcomes of TDP in a large series of prospectively sequenced tumors from patients with pancreatic ductal adenocarcinomas (PDAC). Whole-genome sequencing (WGS) was performed in 530 PDAC cases from the PanCuRx Initiative, COMPASS and PanGen/POG trials in Canada. Of 530 cases, 52 were identified as TDP (9.8%; 13 resected, 39 advanced). Etiological subgroups of TDP included BRCA1 (n = 9), CCNE1 (n = 4), and unknown (n = 39). Presence of TDP was not prognostic in resected specimens (p = 0.77) compared with non-HRD and non-TDP cases, described as typicals. In advanced cases, when stratified for only classical subtype cases, platinum therapy was correlated with longer response in non-BRCA1 TDP vs. typicals (p = 0.0036). There was no difference in overall survival between TDP and typicals (p = 0.5).TDP represents a potential novel targetable subgroup for chemotherapy selection in PDAC.

Severus detects somatic structural variation and complex rearrangements in cancer genomes using long-read sequencing

For the detection of somatic structural variation (SV) in cancer genomes, long-read sequencing is advantageous over short-read sequencing with respect to mappability and variant phasing. However, most current long-read SV detection methods are not developed for the analysis of tumor genomes characterized by complex rearrangements and heterogeneity. Here, we present Severus, a breakpoint graph-based algorithm for somatic SV calling from long-read cancer sequencing. Severus works with matching normal samples, supports unbalanced cancer karyotypes, can characterize complex multibreak SV patterns and produces haplotype-specific calls. On a comprehensive multitechnology cell line panel, Severus consistently outperforms other long-read and short-read methods in terms of SV detection F1 score (harmonic mean of the precision and recall). We also illustrate that compared to long-read methods, short-read sequencing systematically misses certain classes of somatic SVs, such as insertions or clustered rearrangements. We apply Severus to several clinical cases of pediatric leukemia/lymphoma, revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Genomic insights and the conservation potential of captive breeding: The case of Chinese alligator

Despite 40 years of conservation of the critically endangered Chinese alligator (Alligator sinensis), the genomic underpinnings of its status remained uncharted. Genome sequencing data of 244 individuals uncovered relatively low overall genomic diversity/heterozygosity and long runs of homozygosity, with captive populations exhibiting higher heterozygosity and smaller inbreeding coefficients compared to wild individuals. The decreased level of inbreeding in the captive population demonstrates the contribution of the large captive breeding population. The estimated recent effective population size was around a few dozen. To combat challenges of inbreeding depression and reduced adaptability, we used genome-wide SNP-based kinship analysis on captive populations to enable a genome-informed breeding program that minimizes inbreeding. Long-term field monitoring revealed that the Chinese government greatly advanced the conservation of A. sinensis through conservation measures and reintroduction programs. Our research enriches the understanding of the Chinese alligator's genetic landscape, offering invaluable genomic resources for breeding and conservation strategies.

Unveiling the intriguing relationship: oncogenic KRAS, morphological shifts, and mutational complexity in pancreatic mucinous cystic neoplasms

Pancreatic ductal adenocarcinoma (PDAC) often arises from preexisting cystic lesions such as intraductal papillary mucinous neoplasms (IPMN) and mucinous cystic neoplasms (MCN). This study investigated the molecular heterogeneity and mutational landscape of MCN in relation to PDAC, highlighting the significance of KRAS mutations in tumor progression. Utilizing targeted next-generation sequencing on low-grade MCN and invasive PDAC samples, we identified a substantial overlap in mutational profiles, particularly mutations in KRAS, TP53, and FBXW7. Specifically, 69.2% of MCN exhibited somatic mutations, with KRAS mutations being a predominant oncogenic driver. The characterization of mutant versus wildtype KRAS variant allele frequencies (VAF) indicated higher mutation levels in PDAC compared to MCN, suggesting an evolutionary trajectory toward malignancy. Further histological analysis of 12 additional MCN cases revealed significant intratumor heterogeneity, with variant KRAS mutation distributions correlating with distinct cellular morphologies and dysplastic features. Additionally, we explored the potential of liquid biopsies, demonstrating a concordance rate of 71.4% for KRAS mutation detection in circulating tumor DNA (ctDNA) relative to tissue biopsies across cohorts. Our findings underscore the relevance of evaluating KRAS mutations-herein referred to as VAF per microdissected region-as they relate to histopathological markers of dysplasia, contributing to improved stratification of pancreatic lesions and facilitating personalized treatment strategies. In conclusion, this comprehensive analysis of MCN highlights the importance of KRAS as a crucial biomarker for both malignant progression and therapeutic decision-making in pancreatic pathology. Ultimately, our study suggests that characterizing the mutational landscape and histological features of MCN can enhance early detection and intervention strategies for at-risk patients. © 2025 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Multiomic profiling of T cell lymphoma after therapy with anti-BCMA CAR T cells and GPRC5D-directed bispecific antibody

Chimeric antigen receptor (CAR) T cells and bispecific T cell engagers have become integral components in the treatment of relapsed/refractory multiple myeloma. We report a 63-year-old male who received ciltacabtagene autoleucel CAR T cells and the GPRC5D × CD3 bispecific talquetamab for early relapse of his multiple myeloma. Nine months after CAR T therapy, he developed a symptomatic leukemic peripheral T cell lymphoma with cutaneous and intestinal involvement. Longitudinal single-cell RNA and T cell receptor sequencing of peripheral blood and bone marrow revealed two hyperexpanded CAR-carrying T cell clones. These expanded clones exhibited an exhausted effector-memory T cell transcriptional signature, and the neoplasm itself was sensitive to dexamethasone treatment. The immunophenotypic and transcriptional alterations of these abnormal T cells resembled those of T-large granular lymphocytic leukemia. Spatial transcriptomes of skin lesions confirmed the aberrant CAR-expressing T cells. Whole-genome sequencing revealed three distinct integration sites, within the introns of ZGPAT, KPNA4 and polycomb-associated noncoding RNAs. Before and after CAR T whole-genome analyses implicated clonal outgrowth of a TET2-mutated precursor propelled by additional subclone-specific loss of heterozygosity and other secondary mechanisms. This case highlights the evolution of a CAR-carrying peripheral T cell lymphoma following CAR T cell and bispecific T cell engager therapy, offering critical insights into the clonal evolution from a predisposed hematopoietic precursor to a mature neoplasm.

Evolutionary fingerprints of epithelial-to-mesenchymal transition

Mesenchymal plasticity has been extensively described in advanced epithelial cancers; however, its functional role in malignant progression is controversial1-5. The function of epithelial-to-mesenchymal transition (EMT) and cell plasticity in tumour heterogeneity and clonal evolution is poorly understood. Here we clarify the contribution of EMT to malignant progression in pancreatic cancer. We used somatic mosaic genome engineering technologies to trace and ablate malignant mesenchymal lineages along the EMT continuum. The experimental evidence clarifies the essential contribution of mesenchymal lineages to pancreatic cancer evolution. Spatial genomic analysis, single-cell transcriptomic and epigenomic profiling of EMT clarifies its contribution to the emergence of genomic instability, including events of chromothripsis. Genetic ablation of mesenchymal lineages robustly abolished these mutational processes and evolutionary patterns, as confirmed by cross-species analysis of pancreatic and other human solid tumours. Mechanistically, we identified that malignant cells with mesenchymal features display increased chromatin accessibility, particularly in the pericentromeric and centromeric regions, in turn resulting in delayed mitosis and catastrophic cell division. Thus, EMT favours the emergence of genomic-unstable, highly fit tumour cells, which strongly supports the concept of cell-state-restricted patterns of evolution, whereby cancer cell speciation is propagated to progeny within restricted functional compartments. Restraining the evolutionary routes through ablation of clones capable of mesenchymal plasticity, and extinction of the derived lineages, halts the malignant potential of one of the most aggressive forms of human cancer.

Gastric cancer genomics study using reference human pangenomes

A pangenome is the sum of the genetic information of all individuals in a species or a population. Genomics research has been gradually shifted to a paradigm using a pangenome as the reference. However, in disease genomics study, pangenome-based analysis is still in its infancy. In this study, we introduced a graph-based pangenome GGCPan from 185 patients with gastric cancer. We then systematically compared the cancer genomics study results using GGCPan, a linear pangenome GCPan, and the human reference genome as the reference. For small variant detection and microsatellite instability status identification, there is little difference in using three different genomes. Using GGCPan as the reference had a significant advantage in structural variant identification. A total of 24 candidate gastric cancer driver genes were detected using three different reference genomes, of which eight were common and five were detected only based on pangenomes. Our results showed that disease-specific pangenome as a reference is promising and a whole set of tools are still to be developed or improved for disease genomics study in the pangenome era.

Population Genetic Structure of Three-Spined Sticklebacks in the St. Lawrence: A Gradient of Change

Understanding how environmental gradients shape population genetic structure is critical for elucidating evolutionary dynamics in heterogeneous landscapes. The St. Lawrence Estuary, spanning fluvial, middle, and marine zones, presents a steep salinity gradient that serves as an ideal setting to study such a question. Three-spined sticklebacks (Gasterosteus aculeatus) thrive across these zones, offering an ideal model system to investigate the interplay of gene flow and natural selection in shaping population structure. Using whole-genome resequencing of sticklebacks from 12 sites, this study aimed to resolve fine-scale population structure and investigate how genetic diversity and differentiation are influenced by selection and gene flow. By integrating single nucleotide polymorphisms (SNPs) and structural variants (SVs), we assessed differentiation patterns, examined clinal variation, and evaluated the relative roles of gene flow and selection in shaping population dynamics. Our findings reveal clear genetic differentiation between fluvial and saltwater populations, with Baie-Saint-Paul forming a potential third group. Salinity emerged as a key driver of genetic structure, with clinal variation in allele frequencies suggesting ongoing adaptation along the gradient. Demographic modeling indicated a history of secondary contact with recent and weak gene flow. Structural variants, particularly indels, complemented SNP-based analyses, underscoring their importance in detecting fine-scale population structure. These results highlight the complex interplay of evolutionary forces shaping biodiversity in transitional environments, providing a basis for exploring local adaptation in connected populations and contributing to broader efforts in conservation genomics.

Role of genomic analysis in the classification of well differentiated hepatocellular lesions

BACKGROUND

The distinction of focal nodular hyperplasia (FNH) and hepatocellular adenoma (HCA) from well-differentiated hepatocellular carcinoma (WD-HCC) in noncirrhotic liver can be challenging. High-grade dysplastic nodule (HGDN) in cirrhosis can have overlapping features with WD-HCC. In some cases, HCA diagnosis is evident but glutamine synthetase (GS) staining is indeterminate for β-catenin activation, which does not allow reliable risk assessment. This study examines the role of genomic analysis in better categorization of WD hepatocellular lesions (WDHL).

DESIGN

Genomic analysis using capture-based NGS assay was done in 23 WDHLs that could not be definitely classified based on morphology, reticulin stain and IHC, and were designated as 'atypical hepatocellular neoplasms' (AHNs). GS staining was classified as diffuse homogeneous (moderate to strong staining in >90 % of tumor cells), diffuse heterogeneous (50-90 %), not diffuse (<50 %) and borderline (not clear if more or less than 50 %).

RESULTS

The genomic profile provided additional information for the diagnosis and/or risk assessment enabling a benign diagnosis in 15/23 cases (66 %) and HCC in 4/23 cases (17 %), while the diagnosis remained as atypical in the remaining 4 cases. Of the 4 cases with final HCC diagnosis, findings were suspicious but not diagnostic based on morphology/IHC; additional changes like TERT promoter mutation (n = 2), AXIN mutation (n = 1), CDKN2A loss (n = 2) and copy number alterations (n = 3) helped to support HCC. Of the 15 cases with a final benign diagnosis, the status of β-catenin activation was unclear based on GS stain in 8 cases, 2 of which showed CTNNB1 exon 7 mutation, 1 showed CTNNB1 exon 8 mutation, while genomic changes in 5 cases did not show any evidence of Wnt activation. FNH-like features were seen in 2 cases, but the genomic changes excluded FNH (CTNNB1 and ARID1A mutation). The final diagnosis was unchanged from the initial diagnosis of AHN in 4/23 cases (17 %) as the molecular findings did not favor HCC.

CONCLUSION

Genomic changes were helpful in characterization of WDHLs, supporting HCC in 17 % of cases and clarifying β-catenin activation status in all 7 cases with borderline GS staining. Genomic changes are not specific but can provide diagnostic clues in selected challenging cases that cannot be classified on morphology and IHC. Given the significant treatment implications of distinguishing between HCC and benign/premalignant entities, routine use of genomic analysis in diagnostically challenging settings should be considered.

Genome-wide detection of runs of homozygosity in Ding'an pigs revealed candidate genes relating to meat quality traits

BACKGROUND

Ding'an (DA) pig, a native Chinese breed, is renowned for its excellent meat quality, disease resistance, high reproductive performance, and adaptability. Its meat quality traits hold significant economic value. However, its conservation population has been declining due to the impact of commercialized breeds and African swine fever, which is not conducive to its development and utilization.

RESULTS

This study utilized whole-genome resequencing data from 15 DA pigs to reveal their genetic characteristics and current resource status. We analyzed the length, number, and distribution patterns of Runs of Homozygosity (ROH) in DA pigs, as well as high-frequency ROH regions. The results identified 23,208,098 single nucleotide polymorphisms (SNPs), 4,497,242 insertion and deletion (InDels), 13,622 copy number variation (CNVs), and 399,934 structure variation (SVs). Further analysis revealed relatively high genetic diversity and low inbreeding levels in DA pigs. Through functional gene enrichment analysis of high-frequency ROH regions, we identified multiple candidate genes associated with specific traits in DA pigs, including meat quality (ANKRD1, CPNE5, MYOM1), fat deposition (OBSCN, MAPK4, PNPLA1, PACSIN1, GRM4), and skeletal muscle development (LRPPRC, WNT9A).

CONCLUSIONS

This study conducted whole-genome sequencing and ROH analysis on DA pigs, revealing high genetic diversity and low inbreeding levels within the population. Through functional gene enrichment analysis of high-frequency ROH regions, we identified multiple candidate genes associated with meat quality, fat deposition, and skeletal muscle development. These findings not only enhance our understanding of the genetic mechanisms underlying the unique traits of DA pigs but also provide valuable insights for practical applications. Specifically, the identified candidate genes and genomic regions can guide conservation efforts to maintain genetic diversity and mitigate inbreeding risks. Meanwhile, these genetic insights can be integrated into breeding programs to improve meat quality and other economically important traits, thereby supporting the sustainable development and utilization of DA pigs.

OctopusV and TentacleSV: a one-stop toolkit for multi-sample, cross-platform structural variant comparison and analysis

Structural variants (SVs) significantly influence genomic variability and disease, but their accurate analysis across multiple samples and sequencing platforms remains challenging. We developed OctopusV, a tool that standardizes ambiguous breakend (BND) annotations into canonical SV types (inversions, duplications, translocations) and integrates variant calls using flexible set operations, such as union, intersection, difference, and complement, enabling cohort-specific variant identification. Together with TentacleSV, an automated pipeline, OctopusV provides an end-to-end solution from raw data to final callsets. Evaluations show improved precision, recall, and consistency, highlighting its value in cancer genomics and rare disease diagnostics. Both tools are available at https://github.com/ylab-hi/OctopusV and https://github.com/ylab-hi/TentacleSV.

Multi-tool copy number detection highlights common body size-associated variants in miniature pig breeds from different geographical regions

BACKGROUND

Copy number variations (CNVs) represent a common and highly specific type of variation in the genome, potentially influencing genetic diversity and mammalian phenotypic development. Structural variants, such as deletions, duplications, and insertions, have frequently been highlighted as key factors influencing traits in high-production pigs. However, comprehensive CNV analyses in miniature pig breeds are limited despite their value in biomedical research.

RESULTS

This study performed whole-genome sequencing in 36 miniature pigs from nine breeds from America, Asia and Oceania, and Europe. By employing a multi-tool approach (CNVpytor, Delly, GATK gCNV, Smoove), the accuracy of CNV identification was improved. In total, 34 homozygous CNVs overlapped with exonic regions in all samples, suggesting a role in expressing specific phenotypes such as uniform growth patterns, fertility, or metabolic function. In addition, 386 copy number variation regions (CNVRs) shared by all breeds were detected, covering 33.6 Mb (1.48% of the autosomal genome). Further, 132 exclusive CNVRs were identified for American breeds, 47 for Asian and Oceanian breeds, and 114 for European breeds. Functional enrichment analysis revealed genes within the common CNVRs involved in body height determination and other growth-related parameters. Exclusive CNVRs were located in the region of genes enriched for lipid metabolism in American minipigs, reproductive traits in Asian and Oceanian breeds, and cardiovascular features and body height in European breeds. In the selected groups, quantitative trait loci associated with body size, meat quality, reproduction, and disease susceptibility were highlighted.

CONCLUSION

This investigation of the CNV landscape of minipigs underlines the impact of selective breeding on structural variants and its role in the development of specific breed phenotypes across geographical areas. The multi-tool approach provides a valuable resource for future studies on the effects of artificial selection on livestock genomes.

Genetic parameter estimation and molecular foundation of Double-yolk eggs trait in white leghorn

Double-yolk (DY) eggs refer to the presence of two yolks in an egg, and they are often present in fowl flocks. As abnormal eggs, DY eggs occur frequently in the early stages of egg-laying in hens, as well as in hens with early sexual maturity. In order to understand the inheritant pattern of DY eggs and explore candidate genes associated with DY eggs, we selected over 10,000 white Leghorn (WL) chickens from 4 generations and recorded the data during the egg laying period, including total egg number and the rate of DY egg number during the first 2 months and the period of 18-58 weeks of age (EN2month, EN38, EN58, DY2month, DY38, and DY58), age at first egg (AFE), first egg weight (FEW), and body weight at first egg (BWA). The results of genetic parameter analysis showed that the DY egg rate was a trait with low to medium heritabilities with the values from 0.15 to 0.29. And there were strong positive phenotypic and genetic correlations between DY egg rate and egg production at different age stages, and they were all strongly negatively correlated with AFE. However, the DY egg rate and egg production at different stages had strong positive and negative genetic correlations with BWA and FEW, respectively. We also found that significant differences in these trait values between different generations and cage layers, indicating that generations and cage layers had a certain influence on these traits. Furtherly, we used whole genome-wide association (GWA) analysis to identify genes underlying DY, and 5 candidate genes (EZH2, CNTNAP2, TMEM163, GPC1, and ACMSD) associated with DY2month in WL. Our study improved the understanding of DY eggs in hens, and the genetic parameters of DY eggs, and also provided insights into reducing the production of DY eggs by various selection strategies.

Independent mechanisms of benzimidazole resistance across Caenorhabditis nematodes

Benzimidazoles (BZs), a widely used class of anthelmintic drugs, target beta-tubulin proteins, disrupt microtubule formation, and cause nematode death. In parasitic nematode species, mutations in beta-tubulin genes (e.g., isotype-1 beta-tubulin) are predicted to inhibit BZ binding and are associated with BZ resistance. Similarly, in the free-living nematode Caenorhabditis elegans, mutations in an isotype-1 beta-tubulin ortholog, ben-1, are the primary drivers of BZ resistance. The recurrent association of BZ resistance with beta-tubulins suggests that BZ resistance is repeatedly caused by mutations in beta-tubulin genes, an example of repeated evolution of drug resistance across nematode species. To evaluate the hypothesis of repeated evolution of BZ resistance mediated by beta-tubulin, we identified predicted resistance alleles in beta-tubulin genes across wild strains from three Caenorhabditis species: C. elegans, Caenorhabditis briggsae, and Caenorhabditis tropicalis. We hypothesized that, if these species experienced similar selective pressures, they would evolve resistance to BZs by mutations in any of three beta-tubulin genes (ben-1, tbb-1, and tbb-2). Using high-throughput development assays, we tested the association of predicted beta-tubulin alleles with BZ resistance. We found that a heterogeneous set of variants identified in C. elegans ben-1 were associated with BZ resistance. In C. briggsae, only two variants in ben-1, predicted to encode a premature stop codon (W21stop) and a missense substitution (Q134H), were associated with BZ resistance. In C. tropicalis, two missense variants were identified in ben-1, but neither was associated with BZ resistance. C. briggsae and C. tropicalis might have evolved BZ resistance by mutations in other beta-tubulin genes, but we found that variants in tbb-1 or tbb-2 in these species were not associated with BZ resistance. Our findings reveal a lack of repeated evolution of BZ resistance across the three Caenorhabditis species and highlight the importance of defining BZ resistance mechanisms outside of beta-tubulins.

FAN1-mediated translesion synthesis and POLQ/HELQ-mediated end joining generate interstrand crosslink-induced mutations

To counteract the damaging effects of DNA interstrand crosslinks (ICLs), cells have evolved various specialized ICL repair pathways. However, how ICL repair impacts genetic integrity remains incompletely understood. Here, we determined the mutagenic consequences of psoralen ICL repair in the animal model C. elegans and identify two mutagenic repair mechanisms: (i) translesion synthesis through POLH and REV1/3-mediated bypass, leading to single nucleotide polymorphisms (SNVs), and (ii) end joining via POLQ or HELQ action resulting in deletions. While we found no role for the Fanconi anemia genes FANCD2 and FANCI, disruption of TRAIP, which triggers unloading of the CMG helicase at sites of blocked replication, led to a strikingly altered repair profile, suggesting a role for DNA replication in the etiology of ICL-induced deletions. TRAIP deficiency did not affect SNV formation; instead, we found these SNVs to depend on the functionality of the Fanconi anemia-associated nuclease FAN1.

Benchmarking long-read structural variant calling tools and combinations for detecting somatic variants in cancer genomes

Cancer genomes have a complicated landscape of mutations, including large-scale rearrangements known as structural variants (SVs). These SVs can disrupt genes or regulatory elements, playing a critical role in cancer development and progression. Despite their importance, accurate identification of somatic structural variants (SVs) remains a significant bottleneck in cancer genomics. Long-read sequencing technologies hold great promise in SV discovery, and there is an increasing number of efforts to develop new tools to detect them. In this study, we employ eight widely used SV callers on paired tumor and matched normal samples from both the NCI-H2009 lung cancer cell line and the COLO829 melanoma cell line, the latter of which has a well-established somatic SV truth set. Following separate variation detection in both tumor and normal DNA, the VCF merging procedure and a subtraction method were used to identify candidate somatic SVs. Additionally, we explored different combinations of the tools to enhance the accuracy of true somatic SV detection. Our analysis adopts a comprehensive approach, evaluating the performance of each SV caller across a spectrum of variant types and numbers in finding cancer-related somatic SVs. This study, by comparing eight different tools and their combinations, not only reveals the benefits and limitations of various techniques but also establishes a framework for developing more robust SV calling pipelines. Our findings highlight the strengths and weaknesses of current SV calling tools and suggest that combining multiple tools and testing different combinations can significantly enhance the validation of somatic alterations.

Paracentric inversion disrupting the SHANK2 gene

In this study, we employed a multifaceted approach combining short-read whole genome sequencing (WGS) analyzed using Delly, cytogenomics using Bionano technology, and Sanger sequencing to identify the breakpoints of a balanced de novo paracentric inversion on chromosome 11, spanning approximately 64 Mb (inv11q13.3; q25). This inversion was discovered in a girl who presented with mild intellectual disability (ID), speech and language delays, a delay in motor development and attention deficit hyperactivity disorder (ADHD). Detailed analysis of the breakpoints revealed the disruption of two genes; SHANK2, which is critical for encoding a postsynaptic scaffolding protein at glutamatergic synapses in the brain, and LINC02714, a long non-coding RNA (lncRNA). Although SHANK2 is not listed in the OMIM database as a causative gene to this date, literature reports at least 21 cases where (likely) pathogenic variants in SHANK2 have been identified in patients with neurodevelopmental disorders (NDDs). A loss of function variant of the SHANK2 gene is in line with the clinical presentation of this patient. No additional genetic variants that could explain her phenotype were identified. In conclusion, by combining WGS, cytogenomics and Sanger sequencing techniques, we identified the exact breakpoints of a large inversion providing a likely molecular diagnosis for our patient.

Benchmarking strategies for CNV calling from whole genome bisulfite data in humans

It's important to dissect the relationship between copy number variations (CNVs) and DNA methylation, because both greatly change the dosages of genes and are responsible for diverse human cancers. Although whole genome bisulfite sequencing (WGBS) informs CNVs and DNA methylation, no study has provided a systematic benchmark for detecting CNVs from WGBS data. Herein, based on simulated and real WGBS datasets of 84.62 billion reads, we undertook 714 CNV detections to comprehensively benchmark the performance of 35 strategies, 5 alignment algorithms (bismarkbt2, bsbolt, bsmap, bwameth, and walt) wrapping with 7 CNV detection applications (BreakDancer, cn.mops, CNVkit, CNVnator, DELLY, GASV and Pindel). The results highlighted a subset of strategies that accurately called CNVs depending on numbers, lengths, precision, recall, and F1 scores of CNV detections. We found that bwameth-DELLY and bwameth-BreakDancer were the best strategies for calling deletions, and walt-CNVnator and bismarkbt2-CNVnator were the best strategies for calling duplications. These works provided investigators with useful information to accurately explore CNVs from WGBS data in humans.

Genetic and Anatomical Determinants of Olfaction in Dogs and Wild Canids

Understanding the anatomical and genetic basis of complex phenotypic traits has long been a challenge for biological research. Domestic dogs offer a compelling model as they demonstrate more phenotypic variation than any other vertebrate species. Dogs have been intensely selected for specific traits and abilities, directly or indirectly, over the past 15,000 years since their initial domestication from the gray wolf. Because olfaction plays a central role in critical tasks, such as the detection of drugs, diseases, and explosives, as well as human rescue, we compared relative olfactory capacity across dog breeds and assessed changes to the canine olfactory system to their direct ancestors, wolves, and coyotes. We conducted a cross-disciplinary survey of olfactory anatomy, olfactory receptor (OR) gene variation, and OR gene expression in domestic dogs. Through comparisons to their closest wild canid relatives, the gray wolf and coyote, we show that domestic dogs might have lost functional OR genes commensurate with a documented reduction in nasal morphology as an outcome of the domestication process prior to breed formation. Critically, within domestic dogs alone, we found no genetic or morphological profile shared among functional or genealogical breed groupings, such as scent hounds, that might indicate evidence of any human-directed selection for enhanced olfaction. Instead, our results suggest that superior scent detection dogs likely owe their success to advantageous behavioral traits and training rather than an "olfactory edge" provided by morphology or genes.

Structural variants linked to Alzheimer's disease and other common age-related clinical and neuropathologic traits

BACKGROUND

Alzheimer's disease (AD) is a complex neurodegenerative disorder with substantial genetic influence. While genome-wide association studies (GWAS) have identified numerous risk loci for late-onset AD (LOAD), the functional mechanisms underlying most of these associations remain unresolved. Large genomic rearrangements, known as structural variants (SVs), represent a promising avenue for elucidating such mechanisms within some of these loci.

METHODS

By leveraging data from two ongoing cohort studies of aging and dementia, the Religious Orders Study and Rush Memory and Aging Project (ROS/MAP), we performed genome-wide association analysis testing 20,205 common SVs from 1088 participants with whole genome sequencing (WGS) data. A range of Alzheimer's disease and other common age-related clinical and neuropathologic traits were examined.

RESULTS

First, we mapped SVs across 81 AD risk loci and discovered 22 SVs in linkage disequilibrium (LD) with GWAS lead variants and directly associated with the phenotypes tested. The strongest association was a deletion of an Alu element in the 3'UTR of the TMEM106B gene, in high LD with the respective AD GWAS locus and associated with multiple AD and AD-related disorders (ADRD) phenotypes, including tangles density, TDP-43, and cognitive resilience. The deletion of this element was also linked to lower TMEM106B protein abundance. We also found a 22-kb deletion associated with depression in ROS/MAP and bearing similar association patterns as GWAS SNPs at the IQCK locus. In addition, we leveraged our catalog of SV-GWAS to replicate and characterize independent findings in SV-based GWAS for AD and five other neurodegenerative diseases. Among these findings, we highlight the replication of genome-wide significant SVs for progressive supranuclear palsy (PSP), including markers for the 17q21.31 MAPT locus inversion and a 1483-bp deletion at the CYP2A13 locus, along with other suggestive associations, such as a 994-bp duplication in the LMNTD1 locus, suggestively linked to AD and a 3958-bp deletion at the DOCK5 locus linked to Lewy body disease (LBD) (P = 3.36 × 10-4).

CONCLUSIONS

While still limited in sample size, this study highlights the utility of including analysis of SVs for elucidating mechanisms underlying GWAS loci and provides a valuable resource for the characterization of the effects of SVs in neurodegenerative disease pathogenesis.

Impacts of reproductive systems on grapevine genome and breeding

Diversified reproductive systems can be observed in the plant kingdom and applied in crop breeding; however, their impacts on crop genomic variation and breeding remain unclear. Grapevine (Vitis vinifera L.), a widely planted fruit tree, underwent a shift from dioecism to monoecism during domestication and involves crossing, self-pollination, and clonal propagation for its cultivation. In this study, we discover that the reproductive types, namely, crossing, selfing, and cloning, dramatically impact genomic landscapes and grapevine breeding based on comparative genomic and population genetics of wild grapevine and a complex pedigree of Pinot Noir. The impacts are widely divergent, which show interesting patterns of genomic purging and the Hill-Robertson interference. Selfing reduces genomic heterozygosity, while cloning increases it, resulting in a "double U-shaped" site frequency spectrum (SFS). Crossing and cloning conceal while selfing purges most deleterious and structural burdens. Moreover, the close leakage of large-effect deleterious and structural variations in repulsion phases maintains heterozygous genomic regions in 4.3% of the grapevine genome after successive selfing for nine generations. Our study provides new insights into the genetic basis of clonal propagation and genomic breeding of clonal crops by purging deleterious variants while integrating beneficial variants through various reproductive systems.

Application of Digital Polymerase Chain Reaction (dPCR) in Non-Invasive Prenatal Testing (NIPT)

This article reviews the current applications of the digital polymerase chain reaction (dPCR) in non-invasive prenatal testing (NIPT) and explores its potential to complement or surpass the capabilities of Next-Generation Sequencing (NGS) in prenatal testing. The growing incidence of genetic disorders in maternal-fetal medicine has intensified the demand for precise and accessible NIPT options, which aim to minimize the need for invasive prenatal diagnostic procedures. Cell-free fetal DNA (cffDNA), the core analyte in NIPT, is influenced by numerous factors such as maternal DNA contamination, placental health, and fragment degradation. dPCR, with its inherent precision and ability to detect low-abundance targets, demonstrates robustness against these interferences. Although NGS remains the gold standard due to its comprehensive diagnostic capabilities, its high costs limit widespread use, particularly in resource-limited settings. In contrast, dPCR provides comparable accuracy with lower complexity and expense, making it a promising alternative for prenatal testing.