Variant Review with the Integrative Genomics Viewer

Insights into the role of JAK2-I724T variant in myeloproliferative neoplasms from a unique cohort of New Zealand patients

OBJECTIVES

This study aimed to compile bioinformatic and experimental information for JAK2 missense variants previously reported in myeloproliferative neoplasms (MPN) and determine if germline JAK2-I724T, recently found to be common in New Zealand Polynesians, associates with MPN.

METHODS

For all JAK2 variants found in the literature, gnomAD_exome allele frequencies were extracted and REVEL scores were calculated using the dbNSFP database. We investigated the prevalence of JAK2-I724T in a cohort of 111 New Zealand MPN patients using a TaqMan assay, examined its allelic co-occurrence with JAK2-V617F using Oxford Nanopore sequencing, and modelled the impact of I724T on JAK2 using I-Mutant and ChimeraX software.

RESULTS

Several non-V617F JAK2 variants previously reported in MPN had REVEL scores greater than 0.5, suggesting pathogenicity. JAK2-I724T (REVEL score 0.753) was more common in New Zealand Polynesian MPN patients (n = 2/27; 7.4%) than in other New Zealand patients (n = 0/84; 0%) but less common than expected for healthy Polynesians (n = 56/377; 14.9%). Patients carrying I724T (n = 2), one with polycythaemia vera and one with essential thrombocythaemia, had high-risk MPN. Both patients with JAK2-I724T were also positive for JAK2-V617F, found on the same allele as I724T, as well as separately. In silico modelling did not identify noticeable structural changes that would give JAK2-I724T a gain-of-function.

CONCLUSION

Several non-canonical JAK2 variants with high REVEL scores have been reported in MPN, highlighting the need to further understand their relationship with disease. The JAK2-I724T variant does not drive MPN, but additional investigations are required to exclude any potential modulatory effect on the MPN phenotype.

Exploring the tumor genomic landscape of aggressive prostate cancer by whole-genome sequencing of tissue or liquid biopsies

Treatment resistance remains a major issue in aggressive prostate cancer (PC), and novel genomic biomarkers may guide better treatment selection. Circulating tumor DNA (ctDNA) can provide minimally invasive information about tumor genomes, but the genomic landscape of aggressive PC based on whole-genome sequencing (WGS) of ctDNA remains incompletely characterized. Thus, we here performed WGS of tumor tissue (n = 31) or plasma ctDNA (n = 10) from a total of 41 aggressive PC patients, including 11 hormone-naïve, 15 hormone-sensitive, and 15 castration-resistant patients. Across all variant types, we found progressively more altered tumor genomic profiles in later stages of aggressive PC. The potential driver genes most frequently affected by single-nucleotide variants or insertions/deletions included the known PC-related genes TP53, CDK12, and PTEN and the novel genes COL13A1, KCNH3, and SENP3. Etiologically, aggressive PC was associated with age-related and DNA repair-related mutational signatures. Copy number variants most frequently affected 14q11.2 and 8p21.2, where no well-recognized PC-related genes are located, and also frequently affected regions near the known PC-related genes MYC, AR, TP53, PTEN, and BRCA1. Structural variants most frequently involved not only the known PC-related genes TMPRSS2 and ERG but also the less extensively studied gene in this context, PTPRD. Finally, clinically actionable variants were detected throughout all stages of aggressive PC and in both plasma and tissue samples, emphasizing the potential clinical applicability of WGS of minimally invasive plasma samples. Overall, our study highlights the feasibility of using liquid biopsies for comprehensive genomic characterization as an alternative to tissue biopsies in advanced/aggressive PC.

Frequent protein kinase A regulatory subunit A1 mutations but no GNAS mutations as potential driver in sporadic cardiac myxomas

PURPOSE

Cardiac myxomas (CMs) are the second most common benign primary cardiac tumors, mainly originating within the left atrium. Approximately 5% of CM cases are associated with Carney Complex (CNC), an autosomal dominant multiple neoplasia syndrome often caused by germline mutations in the protein kinase A regulatory subunit 1A (PRKAR1A). Data concerning PRKAR1A alterations in sporadic myxomas are variable and sparse, with PRKAR1A mutations reported to range from 0% to 87%. Therefore, we investigated the frequency of PRKAR1A mutations in sporadic CM using next-generation sequencing (NGS). Additionally, we explored mutations in the catalytic domain of the Protein Kinase A complex (PRKACA) and examined the presence of GNAS mutations as another potential driver.

METHODS AND RESULTS

This study retrospectively collected histological and clinical data from 27 patients with CM. First, we ruled out the possibility of underlying CNC through clinical evaluations and standardized interviews for each patient. Second, we performed PRKAR1A immunohistochemistry (IHC) analysis and graded the reactivity of myxoma cells semi-quantitatively. NGS was then applied to analyze the coding regions of PRKAR1A, PRKACA, and GNAS in all 27 cases. Of the 27 sporadic CM cases, 13 (48%) harbored mutations in PRKAR1A. Among these 13 mutant cases, six displayed more than one mutation in PRKAR1A. Most of the identified mutations resulted in premature stop codons or affected splicing. In PRKAR1A mutant CM cases, the loss of PRKAR1A protein expression was significantly more common. In two cases with missense mutations, protein expression remained preserved. Furthermore, a single mutation was detected in the catalytic domain of the protein kinase A complex, while no GNAS mutations were found.

CONCLUSION

We identified a relatively high frequency of PRKAR1A mutations in sporadic CM. These PRKAR1A mutations may also represent an important oncogenic mechanism in sporadic myxomas, as already known in CM cases associated with CNC.

Genetic basis of clarithromycin resistance in Bacillus anthracis

The high-consequence pathogen Bacillus anthracis causes human anthrax and often results in lethal infections without the rapid administration of effective antimicrobial treatment. Antimicrobial resistance profiling is therefore critical to inform post-exposure prophylaxis and treatment decisions, especially during emergencies such as outbreaks or where intentional release is suspected. Whole-genome sequencing using a rapid long-read sequencer can uncover antimicrobial resistance patterns if genetic markers of resistance are known. To identify genomic markers associated with antimicrobial resistance, we isolated B. anthracis derived from the avirulent Sterne strain with elevated minimal inhibitory concentrations to clarithromycin. Mutants were characterized both phenotypically through broth microdilution susceptibility testing and observations during culturing, as well as genotypically with whole-genome sequencing. We identified two different in-frame insertions in the L22 ribosomal protein-encoding gene rplV, which were subsequently confirmed to be involved in clarithromycin resistance through the reversion of the mutant gene to the parent (drug-susceptible) sequence. Detection of the rplV insertions was possible with rapid long-read sequencing, with a time-to-answer within 3 h. The mutations associated with clarithromycin resistance described here will be used in conjunction with known genetic markers of resistance for other antimicrobials to strengthen the prediction of antimicrobial resistance in B. anthracis.IMPORTANCEThe disease anthrax, caused by the pathogen Bacillus anthracis, is extremely deadly if not treated quickly and appropriately. Clarithromycin is an antibiotic recommended for the treatment and post-exposure prophylaxis of anthrax by the Centers for Disease Control and Prevention; however, little is known about the ability of B. anthracis to develop resistance to clarithromycin or the mechanism of that resistance. The characterization of clarithromycin-resistant isolates presented here provides valuable information for researchers and clinicians in the event of a release of the resistant strain. Additionally, knowledge of the genetic basis of resistance provides a foundation for susceptibility prediction through rapid genome sequencing to inform timely treatment decisions.

The complete sequence and comparative analysis of ape sex chromosomes

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Genome-wide analyses reveals an association between invasive urothelial carcinoma in the Shetland sheepdog and NIPAL1

Naturally occurring canine invasive urinary carcinoma (iUC) closely resembles human muscle invasive bladder cancer in terms of histopathology, metastases, response to therapy, and low survival rate. The heterogeneous nature of the disease has led to the association of large numbers of risk loci in humans, however most are of small effect. There exists a need for new and accurate animal models of invasive bladder cancer. In dogs, distinct breeds show markedly different rates of iUC, thus presenting an opportunity to identify additional risk factors and overcome the locus heterogeneity encountered in human mapping studies. In the association study presented here, inclusive of 100 Shetland sheepdogs and 58 dogs of other breeds, we identify a homozygous protein altering point mutation within the NIPAL1 gene which increases risk by eight-fold (OR = 8.42, CI = 3.12-22.71), accounting for nearly 30% of iUC risk in the Shetland sheepdog. Inclusion of six additional loci accounts for most of the disease risk in the breed and explains nearly 75% of the phenotypes in this study. When combined with sequence data from tumors, we show that variation in the MAPK signaling pathway is an overarching cause of iUC susceptibility in dogs.

Ancient Mycobacterium leprae genome reveals medieval English red squirrels as animal leprosy host

Leprosy, one of the oldest recorded diseases in human history, remains prevalent in Asia, Africa, and South America, with over 200,000 cases every year.1,2 Although ancient DNA (aDNA) approaches on the major causative agent, Mycobacterium leprae, have elucidated the disease's evolutionary history,3,4,5 the role of animal hosts and interspecies transmission in the past remains unexplored. Research has uncovered relationships between medieval strains isolated from archaeological human remains and modern animal hosts such as the red squirrel in England.6,7 However, the time frame, distribution, and direction of transmissions remains unknown. Here, we studied 25 human and 12 squirrel samples from two archaeological sites in Winchester, a medieval English city well known for its leprosarium and connections to the fur trade. We reconstructed four medieval M. leprae genomes, including one from a red squirrel, at a 2.2-fold average coverage. Our analysis revealed a phylogenetic placement of all strains on branch 3 as well as a close relationship between the squirrel strain and one newly reconstructed medieval human strain. In particular, the medieval squirrel strain is more closely related to some medieval human strains from Winchester than to modern red squirrel strains from England, indicating a yet-undetected circulation of M. leprae in non-human hosts in the Middle Ages. Our study represents the first One Health approach for M. leprae in archaeology, which is centered around a medieval animal host strain, and highlights the future capability of such approaches to understand the disease's zoonotic past and current potential.

Exon-junction complex association with stalled ribosomes and slow translation-independent disassembly

Exon junction complexes are deposited at exon-exon junctions during splicing. They are primarily known to activate non-sense mediated degradation of transcripts harbouring premature stop codons before the last intron. According to a popular model, exon-junction complexes accompany mRNAs to the cytoplasm where the first translating ribosome pushes them out. However, they are also removed by uncharacterized, translation-independent mechanisms. Little is known about kinetic and transcript specificity of these processes. Here we tag core subunits of exon-junction complexes with complementary split nanoluciferase fragments to obtain sensitive and quantitative assays for complex formation. Unexpectedly, exon-junction complexes form large stable mRNPs containing stalled ribosomes. Complex assembly and disassembly rates are determined after an arrest in transcription and/or translation. 85% of newly deposited exon-junction complexes are disassembled by a translation-dependent mechanism. However as this process is much faster than the translation-independent one, only 30% of the exon-junction complexes present in cells at steady state require translation for disassembly. Deep RNA sequencing shows a bias of exon-junction complex bound transcripts towards microtubule and centrosome coding ones and demonstrate that the lifetimes of exon-junction complexes are transcript-specific. This study provides a dynamic vision of exon-junction complexes and uncovers their unexpected stable association with ribosomes.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Somatic variants as a cause of drug-resistant epilepsy including mesial temporal lobe epilepsy with hippocampal sclerosis

OBJECTIVE

The contribution of somatic variants to epilepsy has recently been demonstrated, particularly in the etiology of malformations of cortical development. The aim of this study was to determine the diagnostic yield of somatic variants in genes that have been previously associated with a somatic or germline epilepsy model, ascertained from resected brain tissue from patients with multidrug-resistant focal epilepsy.

METHODS

Forty-two patients were recruited across three categories: (1) malformations of cortical development, (2) mesial temporal lobe epilepsy with hippocampal sclerosis, and (3) nonlesional focal epilepsy. Participants were subdivided based on histopathology of the resected brain. Paired blood- and brain-derived DNA samples were sequenced using high-coverage targeted next generation sequencing to high depth (585× and 1360×, respectively). Variants were identified using Genome Analysis ToolKit (GATK4) MuTect-2 and confirmed using high-coverage Amplicon-EZ sequencing.

RESULTS

Sequence data on 41 patients passed quality control. Four somatic variants were validated following amplicon sequencing: within CBL, ALG13, MTOR, and FLNA. The diagnostic yield across 41 patients was 10%, 9% in mesial temporal lobe epilepsy with hippocampal sclerosis and 20% in malformations of cortical development.

SIGNIFICANCE

This study provides novel insights into the etiology of mesial temporal lobe epilepsy with hippocampal sclerosis, highlighting a potential pathogenic role of somatic variants in CBL and ALG13. We also report candidate diagnostic somatic variants in FLNA in focal cortical dysplasia, while providing further insight into the importance of MTOR and related genes in focal cortical dysplasia. This work demonstrates the potential molecular diagnostic value of variants in both germline and somatic epilepsy genes.

Panel Comparative Analysis Tool: An Open-Source Comparative Analysis Tool for Next-Generation Sequencing Panel Target Regions

Multigene next-generation sequencing (NGS) panels have become a routine diagnostic method in the contemporary practice of personalized medicine. To avoid inadequate test choice or interpretation, a detailed understanding of the precise panel target regions is required. However, the necessary bioinformatic expertise is not always available, and publicly accessible and easily interpretable analyses of target regions are scarce. To address this critical knowledge gap, we present the Panel Comparative Analysis Tool (PanelCAT), an open-source application to analyze, visualize, and compare NGS panel DNA target regions. PanelCAT uses Reference Sequence, ClinVar, and Catalogue of Somatic Mutations in Cancer mutation census databases to quantify the exon and mutation coverage of target regions and provides interactive graphical representations and search functions to inspect the results. We demonstrate the utility of PanelCAT by analyzing two large NGS panels (TruSight Oncology 500 and Human Pan Cancer Panel) to validate the advertised target genes, quantify targeted exons and mutations, and identify differences between panels. PanelCAT will enable institutions and researchers to catalog and visualize NGS panel target regions independent of the manufacturer, promote transparency of panel limitations, and share this information with employees and requisitioners.

Next-generation DNA sequencing of Panax samples revealed new genotypes: Burrows-Wheeler Aligner, Python-based abundance and clustering analysis

Background

There are two major species of the Panax genus, namely Panax ginseng and Panax quinquefolius. Other than the nucleic acid test and nucleic acid amplification test, DNA sequencing can be used to authenticate the species of ginseng samples, especially when their physical forms cannot be used for differentiation.

Method

In this work, next generation sequencing was used to obtain millions of reads from fourteen ginseng samples (root, powder, and granule). Then Gaussian Mixture clustering analysis was applied to analyze the reads from each sample.

Results and Discussion

A new genotype has been revealed in this study. Two samples have been authenticated with certainty, while the others may be hybrid in nature as revealed by the clustering results.

Transcriptomic Signature and Pro-Osteoclastic Secreted Factors of Abnormal Bone-Marrow Stromal Cells in Fibrous Dysplasia

Fibrous dysplasia (FD) is a mosaic skeletal disorder caused by somatic activating variants of GNAS encoding for Gαs and leading to excessive cyclic adenosine monophosphate signaling in bone-marrow stromal cells (BMSCs). The effect of Gαs activation in the BMSC transcriptome and how it influences FD lesion microenvironment are unclear. We analyzed changes induced by Gαs activation in the BMSC transcriptome and secretome. RNAseq analysis of differential gene expression of cultured BMSCs from patients with FD and healthy volunteers, and from an inducible mouse model of FD, was performed, and the transcriptomic profiles of both models were combined to build a robust FD BMSC genetic signature. Pathways related to Gαs activation, cytokine signaling, and extracellular matrix deposition were identified. To assess the modulation of several key secreted factors in FD pathogenesis, cytokines and other factors were measured in culture media. Cytokines were also screened in a collection of plasma samples from patients with FD, and positive correlations of several cytokines to their disease burden score, as well as to one another and bone turnover markers, were found. These data support the pro-inflammatory, pro-osteoclastic behavior of FD BMSCs and point to several cytokines and other secreted factors as possible therapeutic targets and/or circulating biomarkers for FD.

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Self-incompatibility phenotypes of SRK mutants can be predicted with high accuracy

Only very limited information is available on why some non-synonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on S-locus receptor kinase, SRK, which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 299 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed and analyzed to determine the genotype and self-incompatibility phenotype of each transformant. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. These mutations were likely to underlie the self-incompatibility defect as they caused significant changes to amino acid properties. Such findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of function. In addition, this study showed the RandomForest and Extreme Gradient Boosting methods could predict self-incompatibility phenotypes of SRK mutants with high accuracy.

A chromosome-scale fishing cat reference genome for the evaluation of potential germline risk variants

The fishing cat, Prionailurus viverrinus, faces a population decline, increasing the importance of maintaining healthy zoo populations. Unfortunately, zoo-managed individuals currently face a high prevalence of transitional cell carcinoma (TCC), a form of bladder cancer. To investigate the genetics of inherited diseases among captive fishing cats, we present a chromosome-scale assembly, generate the pedigree of the zoo-managed population, reaffirm the close genetic relationship with the Asian leopard cat (Prionailurus bengalensis), and identify 7.4 million single nucleotide variants (SNVs) and 23,432 structural variants (SVs) from whole genome sequencing (WGS) data of healthy and TCC cats. Only BRCA2 was found to have a high recurrent number of missense mutations in fishing cats diagnosed with TCC when compared to inherited human cancer risk variants. These new fishing cat genomic resources will aid conservation efforts to improve their genetic fitness and enhance the comparative study of feline genomes.

Multi-omics analysis reveals that linoleic acid metabolism is associated with variations of trained immunity induced by distinct BCG strains

Trained immunity is one of the mechanisms by which BCG vaccination confers persistent nonspecific protection against diverse diseases. Genomic differences between the different BCG vaccine strains that are in global use could result in variable protection against tuberculosis and therapeutic effects on bladder cancer. In this study, we found that four representative BCG strains (BCG-Russia, BCG-Sweden, BCG-China, and BCG-Pasteur) covering all four genetic clusters differed in their ability to induce trained immunity and nonspecific protection. The trained immunity induced by BCG was associated with the Akt-mTOR-HIF1α axis, glycolysis, and NOD-like receptor signaling pathway. Multi-omics analysis (epigenomics, transcriptomics, and metabolomics) showed that linoleic acid metabolism was correlated with the trained immunity-inducing capacity of different BCG strains. Linoleic acid participated in the induction of trained immunity and could act as adjuvants to enhance BCG-induced trained immunity, revealing a trained immunity-inducing signaling pathway that could be used in the adjuvant development.

Genomic Characterization of IS6110 Insertions in Mycobacterium orygis

Mycobacterium orygis, a subspecies of the Mycobacterium tuberculosis complex (MTBC), has emerged as a significant concern in the context of One Health, with implications for zoonosis or zooanthroponosis or both. MTBC strains are characterized by the unique insertion element IS6110, which is widely used as a diagnostic marker. IS6110 transposition drives genetic modifications in MTBC, imparting genome plasticity and profound biological consequences. While IS6110 insertions are customarily found in the MTBC genomes, the evolutionary trajectory of strains seems to correlate with the number of IS6110 copies, indicating enhanced adaptability with increasing copy numbers. Here, we present a comprehensive analysis of IS6110 insertions in the M. orygis genome, utilizing ISMapper, and elucidate their genetic consequences in promoting successful host adaptation. Our study encompasses a panel of 67 paired-end reads, comprising 11 isolates from our laboratory and 56 sequences downloaded from public databases. Among these sequences, 91% exhibited high-copy, 4.5% low-copy, and 4.5% lacked IS6110 insertions. We identified 255 insertion loci, including 141 intragenic and 114 intergenic insertions. Most of these loci were either unique or shared among a limited number of isolates, potentially influencing strain behavior. Furthermore, we conducted gene ontology and pathway analysis, using eggNOG-mapper 5.0, on the protein sequences disrupted by IS6110 insertions, revealing 63 genes involved in diverse functions of Gene Ontology and 45 genes participating in various KEGG pathways. Our findings offer novel insights into IS6110 insertions, their preferential insertion regions, and their impact on metabolic processes and pathways, providing valuable knowledge on the genetic changes underpinning IS6110 transposition in M. orygis.

Multiple independent de novo mutations are associated with the development of schistosoma reflexum, a lethal syndrome in cattle

Schistosoma reflexum (SR) is a lethal congenital syndrome characterized by U-shaped dorsal retroflexion of the spine and exposure of abdominal viscera. SR is usually associated with severe dystocia. The syndrome is thought to be inherited as a Mendelian trait. We collected a series of 23 SR-affected calves from four breeds (20 Holstein, one Red Danish, one Limousin, one Romagnola) and performed whole-genome sequencing (WGS). WGS was performed on 51 cattle, including 14 cases with parents (trio-based; Group 1) and nine single cases (solo-based; Group 2). Sequencing-based genome-wide association studies with 20 Holstein cases and 154 controls showed no association (above Bonferroni threshold; P-value<3 ×10-09). Assuming a monogenic recessive inheritance, no region of shared homozygosity was observed, suggesting heterogeneity. Alternatively, the presence of possible dominant acting de novo mutations were assessed. In Group 1, heterozygous private variants, absent in both parents, were found in seven cases. These involved the ACTL6A, FLNA, GLG1, IQSEC2, MAST3, MBTPS2, and MLLT1 genes. In addition, heterozygous private variants affecting the genes DYNC1LI1, PPP2R2B, SCAF8, SUGP1, and UBP1 were identified in five cases from Group 2. The detected frameshift and missense variants are predicted to cause haploinsufficiency. Each of these 12 affected genes belong to the class of haploinsufficient loss-of-function genes or are involved in embryonic and pre-weaning lethality or are known to be associated with severe malformation syndromes in humans and/or mice. This study presents for the first time a detailed genomic evaluation of bovine SR, suggesting that independent de novo mutations may explain the sporadic occurrence of SR in cattle.

Multi-omics panoramic analysis of HBV integration, transcriptional regulation, and epigenetic modifications in PLC/PRF/5 cell line

The clearance or transcriptional silencing of integrated HBV DNA is crucial for achieving a functional cure in patients with chronic hepatitis B and reducing the risk of hepatocellular carcinoma development. The PLC/PRF/5 cell line is commonly used as an in vitro model for studying HBV integration. In this study, we employed a range of multi-omics techniques to gain a panoramic understanding of the characteristics of HBV integration in PLC/PRF/5 cells and to reveal the transcriptional regulatory mechanisms of integrated HBV DNA. Transcriptome long-read sequencing (ONT) was conducted to analyze and characterize the transcriptional activity of different HBV DNA integration sites in PLC/PRF/5 cells. Additionally, we collected data related to epigenetic regulation, including whole-genome bisulfite sequencing (WGBS), histone chromatin immunoprecipitation sequencing (ChIP-seq), and assays for transposase-accessible chromatin using sequencing (ATAC-seq), to explore the potential mechanisms involved in the transcriptional regulation of integrated HBV DNA. Long-read RNA sequencing analysis revealed significant transcriptional differences at various integration sites in the PLC/PRF/5 cell line, with higher HBV DNA transcription levels at integration sites on chr11, chr13, and the chr13/chr5 fusion chromosome t (13:5). Combining long-read DNA and RNA sequencing results, we found that transcription of integrated HBV DNA generally starts downstream of the SP1, SP2, or XP promoters. ATAC-seq data confirmed that chromatin accessibility has limited influence on the transcription of integrated HBV DNA in the PLC/PRF/5 cell line. Analysis of WGBS data showed that the methylation intensity of integrated HBV DNA was highly negatively correlated with its transcription level (r = -0.8929, p = 0.0123). After AzaD treatment, the transcription level of integrated HBV DNA significantly increased, especially for the integration chr17, which had the highest level of methylation. Through ChIP-seq data, we observed the association between histone modification of H3K4me3 and H3K9me3 with the transcription of integrated HBV DNA. Our findings suggest that the SP1, SP2 and XP in integrated HBV DNA, methylation level of surrounding host chromosome, and histone modifications affect the transcription of integrated HBV DNA in PLC/PRF/5 cells. This provides important clues for future studies on the expression and regulatory mechanisms of integrated HBV.

Whole genome resequencing reveals candidate genes for postaxial polydactyly in Large White pigs

Polydactyly is a genetic abnormality that affects both pig welfare and industry profits. Despite efforts to explore the genetic basis of pig polydactyly, progress remains limited. In this study, we analyzed a group of Large White pigs with postaxial polydactyly, including 29 cases and 79 controls from 24 families. High-depth sequencing was performed on 20 pigs, while low-depth sequencing was improved through imputation for the remaining pigs. A genome-wide association study (GWAS) and genetic differentiation were conducted using the resequencing dataset, resulting in the identification of 48 significantly associated SNPs and 27 candidate regions. The genetic differentiation regions on chromosomes 5 and 18, which harbored GWAS-identified SNPs, were delineated as confidence regions. The confidence region at Chr18: 1.850-1.925 Mb covers the fifth intron of LMBR1, a gene that contains an important regulatory element for SHH, known as ZRS. Mutations in this ZRS have been found to cause polydactyly in animals and humans. Therefore, we propose LMBR1 as a prospective candidate gene for postaxial polydactyly. These findings emphasize the importance of exploring the role of ZRS within LMBR1 in the pathogenesis of polydactyly in pigs.

Histologic and Genomic Analysis of Conjunctival SCC in African and American Cohorts Reveal UV Light and HPV Signatures and High Tumor Mutation Burden

Purpose

Conjunctival squamous cell carcinoma (conjSCC) is more prevalent and aggressive in sub-Saharan African countries compared with the rest of the world. This study aims to compare the genomic, immunophenotypic, and histologic features between patients from the United States and Ethiopia, to identify etiopathogenic mechanisms and unveil potential treatment strategies.

Methods

We compared histologic features and mutational profiles using whole exome sequencing, high-risk human papillomavirus (HPV) status, PD-L1 expression, and tumor-infiltrating lymphocytes in conjSCC tumors of patients from Ethiopia (ETH; n = 25) and the United States (from MD Anderson [the MDA cohort]; n = 29). Genomic alterations were compared with SCCs from other anatomic sites using data from The Cancer Genome Atlas.

Results

Solar elastosis was seen in 78% of ETH and 10% of MDA samples. Thicker tumors had higher density of CD8+ and CD3+ cells. HPV status was similar between the cohorts (ETH = 21% and MDA = 28%). The mean tumor mutation burden (TMB) was significantly higher in conjSCC (3.01/Mb, log10) and cutaneous SCC compared other SCC subtypes. ETH samples had higher TMB compared to the MDA cohort (3.34 vs. 2.73). Mutations in genes associated with ultraviolet light (UV) signature were most frequently encountered (SBS7b = 74% and SBS7a = 72%), with higher prevalence in the ETH cohort, whereas SBS2 and SBS13 signatures were more common among MDA HPV+ conjSCCs.

Conclusions

Our findings suggest that UV exposure may play a major role in conjSCC, with a higher prevalence in the ETH cohort compared with the MDA cohort, where HPV also contributes.

Three Modes of Viral Adaption by the Heart

Viruses elicit long-term adaptive responses in the tissues they infect. Understanding viral adaptions in humans is difficult in organs such as the heart, where primary infected material is not routinely collected. In search of asymptomatic infections with accompanying host adaptions, we mined for cardio-pathogenic viruses in the unaligned reads of nearly one thousand human hearts profiled by RNA sequencing. Among virus-positive cases (~20%), we identified three robust adaptions in the host transcriptome related to inflammatory NFκB signaling and post-transcriptional regulation by the p38-MK2 pathway. The adaptions are not determined by the infecting virus, and they recur in infections of human or animal hearts and cultured cardiomyocytes. Adaptions switch states when NFκB or p38-MK2 are perturbed in cells engineered for chronic infection by the cardio-pathogenic virus, coxsackievirus B3. Stratifying viral responses into reversible adaptions adds a targetable systems-level simplification for infections of the heart and perhaps other organs.

Development, testing and validation of a targeted NGS-panel for the detection of actionable mutations in lung cancer (NSCLC) using anchored multiplex PCR technology in a multicentric setting

Lung cancer is a paradigm for a genetically driven tumor. A variety of drugs were developed targeting specific biomarkers requiring testing for tumor genetic alterations in relevant biomarkers. Different next-generation sequencing technologies are available for library generation: 1) anchored multiplex-, 2) amplicon based- and 3) hybrid capture-based-PCR. Anchored multiplex PCR-based sequencing was investigated for routine molecular testing within the national Network Genomic Medicine Lung Cancer (nNGM). Four centers applied the anchored multiplex ArcherDX-Variantplex nNGMv2 panel to re-analyze samples pre-tested during routine diagnostics. Data analyses were performed by each center and compiled centrally according to study design. Pre-defined standards were utilized, and panel sensitivity was determined by dilution experiments. nNGMv2 panel sequencing was successful in 98.9% of the samples (N = 90). With default filter settings, all but two potential MET exon 14 skipping variants were identified at similar allele frequencies. Both MET variants were found with an adapted calling filter. Three additional variants (KEAP1, STK11, TP53) were called that were not identified in pre-testing analyses. Only total DNA amount but not a qPCR-based DNA quality score correlated with average coverage. Analysis was successful with a DNA input as low as 6.25 ng. Anchored multiplex PCR-based sequencing (nNGMv2) and a sophisticated user-friendly Archer-Analysis pipeline is a robust and specific technology to detect tumor genetic mutations for precision medicine of lung cancer patients.

Mutation analysis and clinical profile of South African patients with Neurofibromatosis type 1 (NF1) phenotype

Neurofibromatosis type 1 (NF1) is an autosomal dominant genetic condition with complete age-dependent penetrance, variable expressivity and a global prevalence of ∼1/3,000. It is characteriszed by numerous café-au-lait macules, skin freckling in the inguinal or axillary regions, Lisch nodules of the iris, optic gliomas, neurofibromas, and tumour predisposition. The diagnostic testing strategy for NF1 includes testing for DNA single nucleotide variants (SNVs), copy number variants (CNVs) as well as RNA analysis for deep intronic and splice variants, which can cumulatively identify the causative variant in 95% of patients. In the present study, NF1 patients were screened using a next-generation sequencing (NGS) assay targeting NF1 exons and intron/exon boundaries for SNV and NF1 multiple ligation-dependent probe amplification (MLPA) analysis for CNV detection. Twenty-six unrelated Southern African patients clinically suspected of having NF1, based on the clinical diagnostic criteria developed by the National Institute of Health (NIH), were included in the current study. A detection rate of 58% (15/26) was obtained, with SNVs identified in 80% (12/15) using a targeted gene panel and NF1 gene deletion in 20% (3/15) identified using MLPA. Ten patients (38%) had no variants identified, although they met NF1 diagnostic criteria. One VUS was identified in this study in a patient that met NF1 diagnostic criteria, however there was no sufficient information to classify variant as pathogenic. The clinical features of Southern African patients with NF1 are similar to that of the known NF1 phenotype, with the exception of a lower frequency of plexiform neurofibromas and a higher frequency of developmental/intellectual disability compared to other cohorts. This is the first clinical and molecular characterisation of a Southern African ancestry NF1 cohort using both next-generation sequencing and MLPA analysis. A significant number of patients remained without a diagnosis following DNA-level testing. The current study offers a potential molecular testing strategy for our low resource environment that could benefit a significant proportion of patients who previously only received a clinical diagnosis without molecular confirmation.

PRMT1 promotes pancreatic cancer development and resistance to chemotherapy

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal types of cancer, and novel treatment regimens are direly needed. Epigenetic regulation contributes to the development of various cancer types, but its role in the development of and potential as a therapeutic target for PDAC remains underexplored. Here, we show that PRMT1 is highly expressed in murine and human pancreatic cancer and is essential for cancer cell proliferation and tumorigenesis. Deletion of PRMT1 delays pancreatic cancer development in a KRAS-dependent mouse model, and multi-omics analyses reveal that PRMT1 depletion leads to global changes in chromatin accessibility and transcription, resulting in reduced glycolysis and a decrease in tumorigenic capacity. Pharmacological inhibition of PRMT1 in combination with gemcitabine has a synergistic effect on pancreatic tumor growth in vitro and in vivo. Collectively, our findings implicate PRMT1 as a key regulator of pancreatic cancer development and a promising target for combination therapy.

Molecular and clinicopathologic characteristics of CNS embryonal tumors with BRD4::LEUTX fusion

Central nervous system (CNS) embryonal tumors are a heterogeneous group of high-grade malignancies, and the increasing clinical use of methylation profiling and next-generation sequencing has led to the identification of molecularly distinct subtypes. One proposed tumor type, CNS tumor with BRD4::LEUTX fusion, has been described. As only a few CNS tumors with BRD4::LEUTX fusions have been described, we herein characterize a cohort of 9 such cases (4 new, 5 previously published) to further describe their clinicopathologic and molecular features. We demonstrate that CNS embryonal tumor with BRD4::LEUTX fusion comprises a well-defined methylation class/cluster. We find that patients are young (4 years or younger), with large tumors at variable locations, and frequently with evidence of leptomeningeal/cerebrospinal fluid (CSF) dissemination. Histologically, tumors were highly cellular with high-grade embryonal features. Immunohistochemically, 5/5 cases showed synaptophysin and 4/5 showed OLIG2 expression, thus overlapping with CNS neuroblastoma, FOXR2-activated. DNA copy number profiles were generally flat; however, two tumors had chromosome 1q gains. No recurring genomic changes, besides the presence of the fusion, were found. The LEUTX portion of the fusion transcript was constant in all cases assessed, while the BRD4 portion varied but included a domain with proto-oncogenic activity in all cases. Two patients with clinical follow up available had tumors with excellent response to chemotherapy. Two of our patients were alive without evidence of recurrence or progression after gross total resection and chemotherapy at 16 and 33 months. One patient relapsed, and the last of our four patients died of disease one month after diagnosis. Overall, this case series provides additional evidence for this as a distinct tumor type defined by the presence of a specific fusion as well as a distinct DNA methylation signature. Studies on larger series are required to further characterize these tumors.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

Repeat expansions in AR, ATXN1, ATXN2 and HTT in Norwegian patients diagnosed with amyotrophic lateral sclerosis

Genetic repeat expansions cause neuronal degeneration in amyotrophic lateral sclerosis as well as other neurodegenerative disorders such as spinocerebellar ataxia, Huntington's disease and Kennedy's disease. Repeat expansions in the same gene can cause multiple clinical phenotypes. We aimed to characterize repeat expansions in a Norwegian amyotrophic lateral sclerosis cohort. Norwegian amyotrophic lateral sclerosis patients (n = 414) and neurologically healthy controls adjusted for age and gender (n = 713) were investigated for repeat expansions in AR, ATXN1, ATXN2 and HTT using short read exome sequencing and the ExpansionHunter software. Five amyotrophic lateral sclerosis patients (1.2%) and two controls (0.3%) carried ≥36 repeats in HTT (P = 0.032), and seven amyotrophic lateral sclerosis patients (1.7%) and three controls (0.4%) carried ≥29 repeats in ATXN2 (P = 0.038). One male diagnosed with amyotrophic lateral sclerosis carried a pathogenic repeat expansion in AR, and his diagnosis was revised to Kennedy's disease. In ATXN1, 50 amyotrophic lateral sclerosis patients (12.1%) and 96 controls (13.5%) carried ≥33 repeats (P = 0.753). None of the patients with repeat expansions in ATXN2 or HTT had signs of Huntington's disease or spinocerebellar ataxia type 2, based on a re-evaluation of medical records. The diagnosis of amyotrophic lateral sclerosis was confirmed in all patients, with the exception of one patient who had primary lateral sclerosis. Our findings indicate that repeat expansions in HTT and ATXN2 are associated with increased likelihood of developing amyotrophic lateral sclerosis. Further studies are required to investigate the potential relationship between HTT repeat expansions and amyotrophic lateral sclerosis.

The transcriptomic insight into the differential susceptibility of African Swine Fever in inbred pigs

African swine fever (ASF) is a global threat to animal health and food security. ASF is typically controlled by strict biosecurity, rapid diagnosis, and culling of affected herds. Much progress has been made in developing modified live virus vaccines against ASF. There is host variation in response to ASF infection in the field and under controlled conditions. To better understand the dynamics underlying this host differential morbidity, whole transcriptome profiling was carried out in twelve immunized and five sham immunized pigs. Seventeen MHC homozygous inbred Large white Babraham pigs were sampled at three time points before and after the challenge. The changes in the transcriptome profiles of infected animals were surveyed over time. In addition, the immunization effect on the host response was studied as well among the contrasts of all protection subgroups. The results showed two promising candidate genes to distinguish between recovered and non-recovered pigs after infection with a virulent African swine fever virus (ASFV) pre-infection: HTRA3 and GFPT2 (padj < 0.05). Variant calling on the transcriptome assemblies showed a two-base pair insertion into the ACOX3 gene closely located to HTRA3 that may regulate its expression as a putative genomic variant for ASF. Several significant DGEs, enriched gene ontology (GO) terms, and KEGG pathways at 1 day and 7 days post-infection, compared to the pre-infection, indicate a significant inflammation response immediately after ASF infection. The presence of the virus was confirmed by the mapping of RNA-Seq reads on two whole viral genome sequences. This was concordant with a higher virus load in the non-recovered animals 7 days post-infection. There was no transcriptome signature on the immunization at pre-infection and 1 day post-infection. More samples and data from additional clinical trials may support these findings.

Erg251 has complex and pleiotropic effects on azole susceptibility, filamentation, and stress response phenotypes

Ergosterol is essential for fungal cell membrane integrity and growth, and numerous antifungal drugs target ergosterol. Inactivation or modification of ergosterol biosynthetic genes can lead to changes in antifungal drug susceptibility, filamentation and stress response. Here, we found that the ergosterol biosynthesis gene ERG251 is a hotspot for point mutations during adaptation to antifungal drug stress within two distinct genetic backgrounds of Candida albicans. Heterozygous point mutations led to single allele dysfunction of ERG251 and resulted in azole tolerance in both genetic backgrounds. This is the first known example of point mutations causing azole tolerance in C. albicans. Importantly, single allele dysfunction of ERG251 in combination with recurrent chromosome aneuploidies resulted in bona fide azole resistance. Homozygous deletions of ERG251 caused increased fitness in low concentrations of fluconazole and decreased fitness in rich medium, especially at low initial cell density. Dysfunction of ERG251 resulted in transcriptional upregulation of the alternate sterol biosynthesis pathway and ZRT2, a Zinc transporter. Notably, we determined that overexpression of ZRT2 is sufficient to increase azole tolerance in C. albicans. Our combined transcriptional and phenotypic analyses revealed the pleiotropic effects of ERG251 on stress responses including cell wall, osmotic and oxidative stress. Interestingly, while loss of either allele of ERG251 resulted in similar antifungal drug responses, we observed functional divergence in filamentation regulation between the two alleles of ERG251 (ERG251-A and ERG251-B) with ERG251-A exhibiting a dominant role in the SC5314 genetic background. Finally, in a murine model of systemic infection, homozygous deletion of ERG251 resulted in decreased virulence while the heterozygous deletion mutants maintain their pathogenicity. Overall, this study provides extensive genetic, transcriptional and phenotypic analysis for the effects of ERG251 on drug susceptibility, fitness, filamentation and stress responses.

Multi-platform whole genome sequencing for tuberculosis clinical and surveillance applications

Whole genome sequencing (WGS) of Mycobacterium tuberculosis offers valuable insights for tuberculosis (TB) control. High throughput platforms like Illumina and Oxford Nanopore Technology (ONT) are increasingly used globally, although ONT is known for higher error rates and is less established for genomic studies. Here we present a study comparing the sequencing outputs of both Illumina and ONT platforms, analysing DNA from 59 clinical isolates in highly endemic TB regions of Thailand. The resulting sequence data were used to profile the M. tuberculosis pairs for their lineage, drug resistance and presence in transmission chains, and were compared to publicly available WGS data from Thailand (n = 1456). Our results revealed isolates that are predominantly from lineages 1 and 2, with consistent drug resistance profiles, including six multidrug-resistant strains; however, analysis of ONT data showed longer phylogenetic branches, emphasising the technologies higher error rate. An analysis incorporating the larger dataset identified fifteen of our samples within six potential transmission clusters, including a significant clade of 41 multi-drug resistant isolates. ONT's extended sequences also revealed strain-specific structural variants in pe/ppe genes (e.g. ppe50), which are candidate loci for vaccine development. Despite some limitations, our results show that ONT sequencing is a promising approach for TB genomic research, supporting precision medicine and decision-making in areas with less developed infrastructure, which is crucial for tackling the disease's significant regional burden.

An atlas of epithelial cell states and plasticity in lung adenocarcinoma

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

Correlation between variant call accuracy and quality parameters in comprehensive cancer genomic profiling tests

Background

Comprehensive genomic profiling (CGP) tests have been widely utilized in clinical practice. In this test, the variant list automatically output from the data analysis pipeline often contains false-positive variants, although the correlation between the quality parameters and prevalence of false-positive variants remains unclear.

Methods

We analyzed 125 CGP tests performed in our laboratory. False-positive variants were manually detected via visual inspection. The quality parameters of both wet and dry processes were also analyzed.

Results

Among the 125 tests, 52 (41.6%) required more than one correction of the called variants, and 21 (16.8%) required multiple corrections. A significant correlation was detected between somatic false-positive variants and quality parameters in the wet (ΔΔCq, pre-capture library peak size, pre-capture library DNA amount, capture library peak size, and capture library concentration) and dry processes (total reads, mapping rates, duplication rates, mean depth, and depth coverage). Capture library concentration and mean depth were strong independent predictors of somatic false-positive variants.

Conclusions

We demonstrated a correlation between somatic false-positive variants and quality parameters in the CGP test. This study facilitates gaining a better understanding of CGP test quality management.

Characterization of CYP2B6 and CYP2A6 Pharmacogenetic Variation in Sub-Saharan African Populations

Genetic variation in CYP2B6 and CYP2A6 is known to impact interindividual response to antiretrovirals, nicotine, and bupropion, among other drugs. However, the full catalogue of clinically relevant pharmacogenetic variants in these genes is yet to be established, especially across African populations. This study therefore aimed to characterize the star allele (haplotype) distribution in CYP2B6 and CYP2A6 across diverse and understudied sub-Saharan African (SSA) populations. We called star alleles from 961 high-depth full genomes using StellarPGx, Aldy, and PyPGx. In addition, we performed CYP2B6 and CYP2A6 star allele frequency comparisons between SSA and other global biogeographical groups represented in the new 1000 Genomes Project high-coverage dataset (n = 2,000). This study presents frequency information for star alleles in CYP2B6 (e.g., *6 and *18; frequency of 21-47% and 2-19%, respectively) and CYP2A6 (e.g., *4, *9, and *17; frequency of 0-6%, 3-10%, and 6-20%, respectively), and predicted phenotypes (for CYP2B6), across various African populations. In addition, 50 potentially novel African-ancestry star alleles were computationally predicted by StellarPGx in CYP2B6 and CYP2A6 combined. For each of these genes, over 4% of the study participants had predicted novel star alleles. Three novel star alleles in CYP2A6 (*54, *55, and *56) and CYP2B6 apiece, and several suballeles were further validated via targeted Single-Molecule Real-Time resequencing. Our findings are important for informing the design of comprehensive pharmacogenetic testing platforms, and are highly relevant for personalized medicine strategies, especially relating to antiretroviral medication and smoking cessation treatment in Africa and the African diaspora. More broadly, this study highlights the importance of sampling diverse African ethnolinguistic groups for accurate characterization of the pharmacogene variation landscape across the continent.

Dissecting genetics of spectrum of epilepsies with eyelid myoclonia by exome sequencing

OBJECTIVE

Epilepsy with eyelid myoclonia (EEM) spectrum is a generalized form of epilepsy characterized by eyelid myoclonia with or without absences, eye closure-induced seizures with electroencephalographic paroxysms, and photosensitivity. Based on the specific clinical features, age at onset, and familial occurrence, a genetic cause has been postulated. Pathogenic variants in CHD2, SYNGAP1, NEXMIF, RORB, and GABRA1 have been reported in individuals with photosensitivity and eyelid myoclonia, but whether other genes are also involved, or a single gene is uniquely linked with EEM, or its subtypes, is not yet known. We aimed to dissect the genetic etiology of EEM.

METHODS

We studied a cohort of 105 individuals by using whole exome sequencing. Individuals were divided into two groups: EEM- (isolated EEM) and EEM+ (EEM accompanied by intellectual disability [ID] or any other neurodevelopmental/psychiatric disorder).

RESULTS

We identified nine variants classified as pathogenic/likely pathogenic in the entire cohort (8.57%); among these, eight (five in CHD2, one in NEXMIF, one in SYNGAP1, and one in TRIM8) were found in the EEM+ subcohort (28.57%). Only one variant (IFIH1) was found in the EEM- subcohort (1.29%); however, because the phenotype of the proband did not fit with published data, additional evidence is needed before considering IFIH1 variants and EEM- an established association. Burden analysis did not identify any single burdened gene or gene set.

SIGNIFICANCE

Our results suggest that for EEM, as for many other epilepsies, the identification of a genetic cause is more likely with comorbid ID and/or other neurodevelopmental disorders. Pathogenic variants were mostly found in CHD2, and the association of CHD2 with EEM+ can now be considered a reasonable gene-disease association. We provide further evidence to strengthen the association of EEM+ with NEXMIF and SYNGAP1. Possible new associations between EEM+ and TRIM8, and EEM- and IFIH1, are also reported. Although we provide robust evidence for gene variants associated with EEM+, the core genetic etiology of EEM- remains to be elucidated.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia

BACKGROUND

Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus.

METHODS

Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region.

RESULTS

Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements.

CONCLUSION

Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

Sirt6 regulates the proliferation of neural precursor cells and cortical neurogenesis in mice

Sirt6, a member of the class III histone deacetylases (HDACs), functions in the regulation of genomic stability, DNA repair, cancer, metabolism and aging. Sirt6 deficiency is lethal, and newborn SIRT6-null cynomolgus monkeys show unfinished brain development. After the generation of a cortex-specific Sirt6 conditional knockout mouse model, we investigated the specific deletion of Sirt6 in NPCs at E10.5. This study found that Sirt6 deficiency causes excessive proliferation of neural precursor cells (NPCs) and retards differentiation. The results suggest that endogenous Sirt6 in NPCs regulates histone acetylation and limits stemness-related genes, including Notch1, in order to participate in NPC fate determination. These findings help elucidate Sirt6's role in brain development and in NPC fate determination while providing data on species generality and differentiation.

A comprehensive review of deep learning-based variant calling methods

Genome sequencing data have become increasingly important in the field of personalized medicine and diagnosis. However, accurately detecting genomic variations remains a challenging task. Traditional variation detection methods rely on manual inspection or predefined rules, which can be time-consuming and prone to errors. Consequently, deep learning-based approaches for variation detection have gained attention due to their ability to automatically learn genomic features that distinguish between variants. In our review, we discuss the recent advancements in deep learning-based algorithms for detecting small variations and structural variations in genomic data, as well as their advantages and limitations.

Homozygous CNP Mutation and Neurodegeneration in Weimaraners: Myelin Abnormalities and Accumulation of Lipofuscin-like Inclusions

A progressive neurological disorder was observed in a male neutered Weimaraner. Clinical signs included fecal incontinence, lethargy, moderate paraparesis, proprioceptive pelvic limb ataxia, falling, cognitive decline, incoordination, decreased interest in food, changes in posture, and episodes of trance-like behavior. Neurologic signs were first observed at approximately 4 years, 10 months of age and progressed slowly. Magnetic resonance imaging showed generalized brain atrophy with areas of white matter pathology. Humane euthanasia was elected at 6 years, 7 months of age due to increasing severity of the neurological signs. Autofluorescent intracellular granules were observed in the cerebral and cerebellar cortexes, optic nerve, and cardiac muscle of the affected dog. These abnormal inclusions in the cerebral cortex and cardiac muscle immunolabeled with antibodies to mitochondrial ATP synthase subunit c protein, like that observed in the neuronal ceroid lipofuscinosis group of lysosomal storage diseases. Immunolabeling also demonstrated pronounced neuroinflammation in brain tissues. The ultrastructural appearances of the disease-related inclusion bodies in the brain and optic nerve were quite variable. The ultrastructure and locations of many of the inclusions in the nervous tissues suggested that they were derived, at least in part, from the myelin surrounding axons. The storage bodies in the cardiac muscle were located in mitochondria-rich regions and consisted of parallel arrays of membrane-like components interspersed with electron-dense flocculent material. The disease was characterized by pronounced abnormalities in the myelin of the brain and optic nerve consisting of distinctive areas of ballooning between the layers of myelin. The whole genome sequence generated from the affected dog contained a homozygous G-to-A missense mutation in CNP, which encodes proteins with CNPase enzyme activity and a structural role in myelin. The mutation predicts a Thr42Met amino acid sequence substitution. Genotyping of archived Weimaraner DNA samples identified an additional G > A variant homozygote with a clinical history and brain lesions similar to those of the proband. Of 304 Weimaraners and over 4000 other dogs of various breeds, the proband and the other Weimaraner that exhibited similar signs were the only two that were homozygous for the CNP missense variant. CNPase immunolabeling was widespread in brain tissues from normal dogs but was undetectable in the same tissues from the proband. Based on the clinical history, fluorescence and electron-microscopy, immunohistochemistry, and molecular genetic findings, the late-onset Weimaraner disorder likely results from the missense mutation that results in CNPase deficiency, leading to myelin abnormalities, accumulation of lysosomal storage bodies, and brain atrophy. Similar disorders have been associated with different CNP variants in Dalmatians and in human subjects.

Unselfish meiotic drive maintains heterozygosity in a parthenogenetic ant

According to Mendel's second law, chromosomes segregate randomly in meiosis. Nonrandom segregation is primarily known for cases of selfish meiotic drive in females, in which particular alleles bias their own transmission into the oocyte1,2. Here, we report a rare example of unselfish meiotic drive for crossover inheritance in the clonal raider ant, Ooceraea biroi. This species produces diploid offspring parthenogenetically via fusion of two haploid nuclei from the same meiosis3. This process should cause rapid genotypic degeneration due to loss of heterozygosity, which results if crossover recombination is followed by random (Mendelian) segregation of chromosomes4,5. However, by comparing whole genomes of mothers and daughters, we show that loss of heterozygosity is exceedingly rare, raising the possibility that crossovers are infrequent or absent in O. biroi meiosis. Using a combination of cytology and whole genome sequencing, we show that crossover recombination is, in fact, common, but that loss of heterozygosity is avoided because crossover products are faithfully co-inherited. This results from a programmed violation of Mendel's law of segregation, such that crossover products segregate together rather than randomly. This discovery highlights an extreme example of cellular "memory" of crossovers, which could be a common yet cryptic feature of chromosomal segregation.

Local Genomic Instability of the SpTransformer Gene Family in the Purple Sea Urchin Inferred from BAC Insert Deletions

The SpTransformer (SpTrf) gene family in the purple sea urchin, Strongylocentrotus purpuratus, encodes immune response proteins. The genes are clustered, surrounded by short tandem repeats, and some are present in genomic segmental duplications. The genes share regions of sequence and include repeats in the coding exon. This complex structure is consistent with putative local genomic instability. Instability of the SpTrf gene cluster was tested by 10 days of growth of Escherichia coli harboring bacterial artificial chromosome (BAC) clones of sea urchin genomic DNA with inserts containing SpTrf genes. After the growth period, the BAC DNA inserts were analyzed for size and SpTrf gene content. Clones with multiple SpTrf genes showed a variety of deletions, including loss of one, most, or all genes from the cluster. Alternatively, a BAC insert with a single SpTrf gene was stable. BAC insert instability is consistent with variations in the gene family composition among sea urchins, the types of SpTrf genes in the family, and a reduction in the gene copy number in single coelomocytes. Based on the sequence variability among SpTrf genes within and among sea urchins, local genomic instability of the family may be important for driving sequence diversity in this gene family that would be of benefit to sea urchins in their arms race with marine microbes.

Bioanthropological analysis of human remains from the archaic and classic period discovered in Puyil cave, Mexico

OBJECTIVES

Determine the geographic place of origin and maternal lineage of prehistoric human skeletal remains discovered in Puyil Cave, Tabasco State, Mexico, located in a region currently populated by Olmec, Zoque and Maya populations.

MATERIALS AND METHODS

All specimens were radiocarbon (14 C) dated (beta analytic), had dental modifications classified, and had an analysis of 13 homologous reference points conducted to evaluate artificial cranial deformation (ACD). Following DNA purification, hypervariable region I (HVR-1) of the mitogenome was amplified and Sanger sequenced. Finally, Next Generation Sequencing (NGS) was performed for total DNA. Mitochondrial DNA (mtDNA) variants and haplogroups were determined using BioEdit 7.2 and IGV software and confirmed with MITOMASTER and WebHome softwares.

RESULTS

Radiocarbon dating (14 C) demonstrated that the inhabitants of Puyil Cave lived during the Archaic and Classic Periods and displayed tabular oblique and tabular mimetic ACD. These pre-Hispanic remains exhibited five mtDNA lineages: A, A2, C1, C1c and D4. Network analysis revealed a close genetic affinity between pre-Hispanic Puyil Cave inhabitants and contemporary Maya subpopulations from Mexico and Guatemala, as well as individuals from Bolivia, Brazil, the Dominican Republic, and China.

CONCLUSIONS

Our results elucidate the dispersal of pre-Hispanic Olmec and Maya ancestors and suggest that ACD practices are closely related to Olmec and Maya practices. Additionally, we conclude that ACD has likely been practiced in the region since the Middle-Archaic Period.

ZNF746 plays cardinal roles on colorectal cancer (CRC) cell invasion and migration and regulates mitochondrial dynamics and morphological changes of CRC cells-Role of combined melatonin and 5-FU regimen

The underlying mechanism of colorectal cells developing into cancer cells has been extensively investigated, yet is still not fully delineated, resulting in the treatment of advanced colorectal cancer (CRC) remains regrettably an unmet need. Zinc Finger Protein 746/Parkin-interacting substrate (ZNF746/PARIS) has previously been identified to play a fundamental role on bladder cancer cell proliferation and metastasis that were effectively inhibited by melatonin (Mel). In this study, we utilized ex vivo/in vivo studies to verify whether the ZNF746 signaling was also crucial in CRC growth/invasion/migration. Tissue-bank specimens showed that the protein expression of ZNF746 was significantly increased in CRC than that of healthy colorectal tissues (p < 0.001). Additionally, in vitro study demonstrated that excessive expression of ZNF746 significantly inhibited mitochondrial activity via (1) interfering with the dynamic balance of mitochondrial fusion/fission and (2) inhibiting the protein expression of MFN1/MFN2/PGC1a (all p < 0.001). Furthermore, we identified that inhibition of ZNF746 protein expression significantly reduced the resistance of CRC cell lines to the anticancer drug of 5-FU (p < 0.001), whereas overexpression of ZNF746 significantly augmented resistance of CRC cells to 5-FU (all p < 0.001). Finally, using the cell culture method, we found that combined Mel and 5-FU was superior to merely one on promoting the CRC cell apoptosis (p < 0.001). Our results confirmed that ZNF746 signaling played a cardinal role of CRC cell proliferation/survival and combined Mel and 5-FU treatment attenuated the resistance of CRC cells to the drug mainly through suppressing this signaling.

A feasible molecular diagnostic strategy for rare genetic disorders within resource-constrained environments

Timely and accurate diagnosis of rare genetic disorders is critical, as it enables improved patient management and prognosis. In a resource-constrained environment such as the South African State healthcare system, the challenge is to design appropriate and cost-effective assays that will enable accurate genetic diagnostic services in patients of African ancestry across a broad disease spectrum. Next-generation sequencing (NGS) has transformed testing approaches for many Mendelian disorders, but this technology is still relatively new in our setting and requires cost-effective ways to implement. As a proof of concept, we describe a feasible diagnostic strategy for genetic disorders frequently seen in our genetics clinics (RASopathies, Cornelia de Lange syndrome, Treacher Collins syndrome, and CHARGE syndrome). The custom-designed targeted NGS gene panel enabled concurrent variant screening for these disorders. Samples were batched during sequencing and analyzed selectively based on the clinical phenotype. The strategy employed in the current study was cost-effective, with sequencing and analysis done at USD849.68 per sample and achieving an overall detection rate of 54.5%. The strategy employed is cost-effective as it allows batching of samples from patients with different diseases in a single run, an approach that can be utilized with rare and less frequently ordered molecular diagnostic tests. The subsequent selective analysis pipeline allowed for timeous reporting back of patients results. This is feasible with a reasonable yield and can be employed for the molecular diagnosis of a wide range of rare monogenic disorders in a resource-constrained environment.

IDR-targeting compounds suppress HPV genome replication via disruption of phospho-BRD4 association with DNA damage response factors

Compounds binding to the bromodomains of bromodomain and extra-terminal (BET) family proteins, particularly BRD4, are promising anticancer agents. Nevertheless, side effects and drug resistance pose significant obstacles in BET-based therapeutics development. Using high-throughput screening of a 200,000-compound library, we identified small molecules targeting a phosphorylated intrinsically disordered region (IDR) of BRD4 that inhibit phospho-BRD4 (pBRD4)-dependent human papillomavirus (HPV) genome replication in HPV-containing keratinocytes. Proteomic profiling identified two DNA damage response factors-53BP1 and BARD1-crucial for differentiation-associated HPV genome amplification. pBRD4-mediated recruitment of 53BP1 and BARD1 to the HPV origin of replication occurs in a spatiotemporal and BRD4 long (BRD4-L) and short (BRD4-S) isoform-specific manner. This recruitment is disrupted by phospho-IDR-targeting compounds with little perturbation of the global transcriptome and BRD4 chromatin landscape. The discovery of these protein-protein interaction inhibitors (PPIi) not only demonstrates the feasibility of developing PPIi against phospho-IDRs but also uncovers antiviral agents targeting an epigenetic regulator essential for virus-host interaction and cancer development.

Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples

BACKGROUND

Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold.

RESULTS

To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value.

CONCLUSIONS

Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.