SAMBLASTER: fast duplicate marking and structural variant read extraction

A Novel Gene ARHGAP44 for Longitudinal Changes in Glycated Hemoglobin (HbA1c) in Subjects without Type 2 Diabetes: Evidence from the Long Life Family Study (LLFS) and the Framingham Offspring Study (FOS)

Glycated hemoglobin (HbA1c) indicates average glucose levels over three months and is associated with insulin resistance and type 2 diabetes (T2D). Longitudinal changes in HbA1c (ΔHbA1c) are also associated with aging processes, cognitive performance, and mortality. We analyzed ΔHbA1c in 1,886 non-diabetic Europeans from the Long Life Family Study to uncover gene variants influencing ΔHbA1c. Using growth curve modeling adjusted for multiple covariates, we derived ΔHbA1c and conducted linkage-guided sequence analysis. Our genome-wide linkage scan identified a significant locus on 17p12. In-depth analysis of this locus revealed a variant rs56340929 (explaining 27% of the linkage peak) in the ARHGAP44 gene that was significantly associated with ΔHbA1c. RNA transcription of ARHGAP44 was associated with ΔHbA1c. The Framingham Offspring Study data further supported these findings on the gene level. Together, we found a novel gene ARHGAP44 for ΔHbA1c in family members without T2D. Follow-up studies using longitudinal omics data in large independent cohorts are warranted.

Analytical and Clinical Validation of a Target-Enhanced Whole Genome Sequencing-Based Comprehensive Genomic Profiling Test

Evaluation of the test performance of the Target enhanced whole-genome sequencing (TE-WGS) assay for comprehensive oncology genomic profiling. The analytical validation of the assay included sensitivity and specificity for single nucleotide variants (SNVs), insertions/deletions (indels), and structural variants (SVs), revealing a revealed a sensitivity of 99.8% for SNVs and 99.2% for indels. The positive predictive value (PPV) was 99.3% SNVs and 98.7% indels. Clinical validation was benchmarked against established orthogonal methods and demonstrated high concordance with reference methods. TE-WGS provides insights beyond targeted panels by comprehensive analysis of key biomarkers and the entire genome encompassing both germline and somatic findings.

Circulating tumour DNA detects somatic variants contributing to spatial and temporal intra-tumoural heterogeneity in head and neck squamous cell carcinoma

Background

As circulating tumour DNA (ctDNA) liquid biopsy analysis is increasingly incorporated into modern oncological practice, establishing the impact of genomic intra-tumoural heterogeneity (ITH) upon data output is paramount. Despite advances in other cancer types the evidence base in head and neck squamous cell carcinoma (HNSCC) remains poor. We sought to investigate the utility of ctDNA to detect ITH in HNSCC.

Methods

In a pilot cohort of 9 treatment-naïve HNSCC patients, DNA from two intra-tumoural sites (core and margin) was whole-exome sequenced. A 9-gene panel was designed to perform targeted sequencing on pre-treatment plasma cell-free DNA and selected post-treatment samples.

Results

Rates of genomic ITH among the 9 patients was high. COSMIC variants from 19 TCGA HNSCC genes demonstrated an 86.9% heterogeneity rate (present in one tumour sub-site only). Across all patients, cell-free DNA (ctDNA) identified 12.9% (range 7.5-19.8%) of tumour-specific variants, of which 55.6% were specific to a single tumour sub-site only. CtDNA identified 79.0% (range: 55.6-90.9%) of high-frequency variants (tumour VAF>5%). Analysis of ctDNA in serial post-treatment blood samples in patients who suffered recurrence demonstrated dynamic changes in both tumour-specific and acquired variants that predicted recurrence ahead of clinical detection.

Conclusion

We demonstrate that a ctDNA liquid biopsy identified spatial genomic ITH in HNSCC and reliably detected high-frequency driver mutations. Serial sampling allowed post-treatment surveillance and early identification of treatment failure.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Genetic Diagnosis of Retinoblastoma Using Aqueous Humour-Findings from an Extended Cohort

The identification of somatic RB1 variation is crucial to confirm the heritability of retinoblastoma. We and others have previously shown that, when tumour DNA is unavailable, cell-free DNA (cfDNA) derived from aqueous humour (AH) can be used to identify somatic RB1 pathogenic variation. Here we report RB1 pathogenic variant detection, as well as cfDNA concentration in an extended cohort of 75 AH samples from 68 patients. We show cfDNA concentration is highly variable and significantly correlated with the collection point of the AH. Cell-free DNA concentrations above 5 pg/µL enabled the detection of 93% of known or expected RB1 pathogenic variants. In AH samples collected during intravitreal chemotherapy treatment (Tx), the yield of cfDNA above 5 pg/µL and subsequent variant detection was low (≤46%). However, AH collected by an anterior chamber tap after one to three cycles of primary chemotherapy (Dx1+) enabled the detection of 75% of expected pathogenic variants. Further limiting our analysis to Dx1+ samples taken after ≤2 cycles (Dx ≤ 2) provided measurable levels of cfDNA in all cases, and a subsequent variant detection rate of 95%. Early AH sampling is therefore likely to be important in maximising cfDNA concentration and the subsequent detection of somatic RB1 pathogenic variants in retinoblastoma patients undergoing conservative treatment.

BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies

Data from both bulk and single-cell whole-genome DNA methylation experiments are under-utilized in many ways. This is attributable to inefficient mapping of methylation sequencing reads, routinely discarded genetic information, and neglected read-level epigenetic and genetic linkage information. We introduce the BISulfite-seq Command line User Interface Toolkit (BISCUIT) and its companion R/Bioconductor package, biscuiteer, for simultaneous extraction of genetic and epigenetic information from bulk and single-cell DNA methylation sequencing. BISCUIT's performance, flexibility and standards-compliant output allow large, complex experimental designs to be characterized on clinical timescales. BISCUIT is particularly suited for processing data from single-cell DNA methylation assays, with its excellent scalability, efficiency, and ability to greatly enhance mappability, a key challenge for single-cell studies. We also introduce the epiBED format for single-molecule analysis of coupled epigenetic and genetic information, facilitating the study of cellular and tissue heterogeneity from DNA methylation sequencing.

Leaf gene expression trajectories during the growing season are consistent between sites and years in American beech

Transcriptomics provides a versatile tool for ecological monitoring. Here, through genome-guided profiling of transcripts mapping to 33 042 gene models, expression differences can be discerned among multi-year and seasonal leaf samples collected from American beech trees at two latitudinally separated sites. Despite a bottleneck due to post-Columbian deforestation, the single nucleotide polymorphism-based population genetic background analysis has yielded sufficient variation to account for differences between populations and among individuals. Our expression analyses during spring-summer and summer-autumn transitions for two consecutive years involved 4197 differentially expressed protein coding genes. Using Populus orthologues we reconstructed a protein-protein interactome representing leaf physiological states of trees during the seasonal transitions. Gene set enrichment analysis revealed gene ontology terms that highlight molecular functions and biological processes possibly influenced by abiotic forcings such as recovery from drought and response to excess precipitation. Further, based on 324 co-regulated transcripts, we focused on a subset of GO terms that could be putatively attributed to late spring phenological shifts. Our conservative results indicate that extended transcriptome-based monitoring of forests can capture diverse ranges of responses including air quality, chronic disease, as well as herbivore outbreaks that require activation and/or downregulation of genes collectively tuning reaction norms maintaining the survival of long living trees such as the American beech.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed

BACKGROUND

The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown.

RESULTS

Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness.

CONCLUSIONS

Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Placentae of small appropriately-grown-for-gestational-age neonates exhibit sexually dimorphic transcriptomic changes representative of placental insufficiency

INTRODUCTION

Previous studies have reported that neonates less than the 25th BWC especially if they were male, were more likely to be associated with birth complications suggesting small neonates often identified as appropriately grown are at risk of adverse outcomes. We have questioned whether smaller neonates not typically categorized as "small for gestational age" may not reach their genetically determined growth due to placental insufficiency.

METHODS

RNA-Seq was performed on the Illumina NovaSeq 600 using term placentae from neonates that were less than the 10th birthweight centile (BWC) (n = 39), between the 10th and the 30th BWC (n = 15) or greater than the 30th BWC (n = 23). Bioinformatic analyses were conducted and statistical significance was assessed at a level of P < 0.05 for single comparisons or FDR <0.05 unless otherwise noted.

RESULTS

Gene set enrichment analysis revealed differences between BWC groups and in relation to the sex of the placenta. Genes associated with hypoxia, inflammatory responses, estrogen responsive genes, and androgen responsive genes were enriched (FDR <0.1) for in placentae of neonates <10th BWC regardless of sex and also in male placentae of neonates between the 10th-30th BWC. Female placenta of neonates between the 10th-30th BWC were comparable to placentae of neonates >30th BWC.

DISCUSSION

These findings provide evidence that small male neonates may be at a greater risk of an adverse outcome than females due to changes in gene expression that are associated with placental dysfunction. The current data raises questions of whether placental pathology for smaller appropriately grown neonates should be scientifically and clinically examined in more depth.

Select EZH2 inhibitors enhance viral mimicry effects of DNMT inhibition through a mechanism involving NFAT:AP-1 signaling

DNA methyltransferase inhibitor (DNMTi) efficacy in solid tumors is limited. Colon cancer cells exposed to DNMTi accumulate lysine-27 trimethylation on histone H3 (H3K27me3). We propose this Enhancer of Zeste Homolog 2 (EZH2)-dependent repressive modification limits DNMTi efficacy. Here, we show that low-dose DNMTi treatment sensitizes colon cancer cells to select EZH2 inhibitors (EZH2is). Integrative epigenomic analysis reveals that DNMTi-induced H3K27me3 accumulates at genomic regions poised with EZH2. Notably, combined EZH2i and DNMTi alters the epigenomic landscape to transcriptionally up-regulate the calcium-induced nuclear factor of activated T cells (NFAT):activating protein 1 (AP-1) signaling pathway. Blocking this pathway limits transcriptional activating effects of these drugs, including transposable element and innate immune response gene expression involved in viral defense. Analysis of primary human colon cancer specimens reveals positive correlations between DNMTi-, innate immune response-, and calcium signaling-associated transcription profiles. Collectively, we show that compensatory EZH2 activity limits DNMTi efficacy in colon cancer and link NFAT:AP-1 signaling to epigenetic therapy-induced viral mimicry.

A phosphatase gene is linked to nectar dihydroxyacetone accumulation in mānuka (Leptospermum scoparium)

Floral nectar composition beyond common sugars shows great diversity but contributing genetic factors are generally unknown. Mānuka (Leptospermum scoparium) is renowned for the antimicrobial compound methylglyoxal in its derived honey, which originates from the precursor, dihydroxyacetone (DHA), accumulating in the nectar. Although this nectar trait is highly variable, genetic contribution to the trait is unclear. Therefore, we investigated key gene(s) and genomic regions underpinning this trait. We used RNAseq analysis to identify nectary-associated genes differentially expressed between high and low nectar DHA genotypes. We also used a mānuka high-density linkage map and quantitative trait loci (QTL) mapping population, supported by an improved genome assembly, to reveal genetic regions associated with nectar DHA content. Expression and QTL analyses both pointed to the involvement of a phosphatase gene, LsSgpp2. The expression pattern of LsSgpp2 correlated with nectar DHA accumulation, and it co-located with a QTL on chromosome 4. The identification of three QTLs, some of the first reported for a plant nectar trait, indicates polygenic control of DHA content. We have established plant genetics as a key influence on DHA accumulation. The data suggest the hypothesis of LsSGPP2 releasing DHA from DHA-phosphate and variability in LsSgpp2 gene expression contributing to the trait variability.

Aberrant non-canonical NF-κB signalling reprograms the epigenome landscape to drive oncogenic transcriptomes in multiple myeloma

In multiple myeloma, abnormal plasma cells establish oncogenic niches within the bone marrow by engaging the NF-κB pathway to nurture their survival while they accumulate pro-proliferative mutations. Under these conditions, many cases eventually develop genetic abnormalities endowing them with constitutive NF-κB activation. Here, we find that sustained NF-κB/p52 levels resulting from such mutations favours the recruitment of enhancers beyond the normal B-cell repertoire. Furthermore, through targeted disruption of p52, we characterise how such enhancers are complicit in the formation of super-enhancers and the establishment of cis-regulatory interactions with myeloma dependencies during constitutive activation of p52. Finally, we functionally validate the pathological impact of these cis-regulatory modules on cell and tumour phenotypes using in vitro and in vivo models, confirming RGS1 as a p52-dependent myeloma driver. We conclude that the divergent epigenomic reprogramming enforced by aberrant non-canonical NF-κB signalling potentiates transcriptional programs beneficial for multiple myeloma progression.

Intrinsic endothelial hyperresponsiveness to inflammatory mediators drives acute episodes in models of Clarkson disease

Clarkson disease, or monoclonal gammopathy-associated idiopathic systemic capillary leak syndrome (ISCLS), is a rare, relapsing-remitting disorder featuring the abrupt extravasation of fluids and proteins into peripheral tissues, which in turn leads to hypotensive shock, severe hemoconcentration, and hypoalbuminemia. The specific leakage factor(s) and pathways in ISCLS are unknown, and there is no effective treatment for acute flares. Here, we characterize an autonomous vascular endothelial defect in ISCLS that was recapitulated in patient-derived endothelial cells (ECs) in culture and in a mouse model of disease. ISCLS-derived ECs were functionally hyperresponsive to permeability-inducing factors like VEGF and histamine, in part due to increased endothelial nitric oxide synthase (eNOS) activity. eNOS blockade by administration of N(γ)-nitro-l-arginine methyl ester (l-NAME) ameliorated vascular leakage in an SJL/J mouse model of ISCLS induced by histamine or VEGF challenge. eNOS mislocalization and decreased protein phosphatase 2A (PP2A) expression may contribute to eNOS hyperactivation in ISCLS-derived ECs. Our findings provide mechanistic insights into microvascular barrier dysfunction in ISCLS and highlight a potential therapeutic approach.

Long-read single-cell RNA sequencing enables the study of cancer subclone-specific genotype and phenotype in chronic lymphocytic leukemia

Bruton's tyrosine kinase (BTK) inhibitors are effective for the treatment of chronic lymphocytic leukemia (CLL) due to BTK's role in B cell survival and proliferation. Treatment resistance is most commonly caused by the emergence of the hallmark BTKC481S mutation that inhibits drug binding. In this study, we aimed to investigate whether the presence of additional CLL driver mutations in cancer subclones harboring a BTKC481S mutation accelerates subclone expansion. In addition, we sought to determine whether BTK-mutated subclones exhibit distinct transcriptomic behavior when compared to other cancer subclones. To achieve these goals, we employ our recently published method (Qiao et al. 2024) that combines bulk DNA sequencing and single-cell RNA sequencing (scRNA-seq) data to genotype individual cells for the presence or absence of subclone-defining mutations. While the most common approach for scRNA-seq includes short-read sequencing, transcript coverage is limited due to the vast majority of the reads being concentrated at the priming end of the transcript. Here, we utilized MAS-seq, a long-read scRNAseq technology, to substantially increase transcript coverage across the entire length of the transcripts and expand the set of informative mutations to link cells to cancer subclones in six CLL patients who acquired BTKC481S mutations during BTK inhibitor treatment. We found that BTK-mutated subclones often acquire additional mutations in CLL driver genes, leading to faster subclone proliferation. When examining subclone-specific gene expression, we found that in one patient, BTK-mutated subclones are transcriptionally distinct from the rest of the malignant B cell population with an overexpression of CLL-relevant genes.

Increased genomic instability and reshaping of tissue microenvironment underlie oncogenic properties of Arid1a mutations

Oncogenic mutations accumulating in many chromatin-associated proteins have been identified in different tumor types. With a mutation rate from 10 to 57%, ARID1A has been widely considered a tumor suppressor gene. However, whether this role is mainly due to its transcriptional-related activities or its ability to preserve genome integrity is still a matter of intense debate. Here, we show that ARID1A is largely dispensable for preserving enhancer-dependent transcriptional regulation, being ARID1B sufficient and required to compensate for ARID1A loss. We provide in vivo evidence that ARID1A is mainly required to preserve genomic integrity in adult tissues. ARID1A loss primarily results in DNA damage accumulation, interferon type I response activation, and chronic inflammation leading to tumor formation. Our data suggest that in healthy tissues, the increased genomic instability that follows ARID1A mutations and the selective pressure imposed by the microenvironment might result in the emergence of aggressive, possibly immune-resistant, tumors.

Human Coronavirus Infection Reorganizes Spatial Genomic Architecture in Permissive Lung Cells

Chromatin conformation capture followed by next-generation sequencing in combination with large-scale polymer simulations (4DHiC) produces detailed information on genomic loci interactions, allowing for the interrogation of 3D spatial genomic structures. Here, Hi-C data was acquired from the infection of fetal lung fibroblast (MRC5) cells with α-coronavirus 229E (CoV229E). Experimental Hi-C contact maps were used to determine viral-induced changes in genomic architecture over a 48-hour time period following viral infection, revealing substantial alterations in contacts within chromosomes and in contacts between different chromosomes. To gain further structural insight and quantify the underlying changes, we applied the 4DHiC polymer simulation method to reconstruct the 3D genomic structures and dynamics corresponding to the Hi-C maps. The models successfully reproduced experimental Hi-C data, including the changes in contacts induced by viral infection. Our 3D spatial simulations uncovered widespread chromatin restructuring, including increased chromosome compactness and A-B compartment mixing arising from infection. Our model also suggests increased spatial accessibility to regions containing interferon-stimulated genes upon infection with CoV229E, followed by chromatin restructuring at later time points, potentially inducing the migration of chromatin into more compact regions. This is consistent with previously observed suppression of gene expression. Our spatial genomics study provides a mechanistic structural basis for changes in chromosome architecture induced by coronavirus infection in lung cells.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Lipocalin 2 (LCN2) confers acquired resistance to almonertinib in NSCLC through LCN2-MMP-9 signaling pathway

Almonertinib, a third-generation epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor, is highly selective for EGFR-activating mutations as well as the EGFR T790M mutation in patients with advanced non-small cell lung cancer (NSCLC). However, the development of resistance inevitably occurs and poses a major obstacle to the clinical efficacy of almonertinib. Therefore, a clear understanding of the mechanism is of great significance to overcome drug resistance to almonertinib in the future. In this study, NCI-H1975 cell lines resistant to almonertinib (NCI-H1975 AR) were developed by concentration-increasing induction and were employed for clarification of underlying mechanisms of acquired resistance. Through RNA-seq analysis, the HIF-1 and TGF-β signaling pathways were significantly enriched by gene set enrichment analysis. Lipocalin-2 (LCN2), as the core node in these two signaling pathways, were found to be positively correlated to almonertinib-resistance in NSCLC cells. The function of LCN2 in the drug resistance of almonertinib was investigated through knockdown and overexpression assays in vitro and in vivo. Moreover, matrix metalloproteinases-9 (MMP-9) was further identified as a critical downstream effector of LCN2 signaling, which is regulated via the LCN2-MMP-9 axis. Pharmacological inhibition of MMP-9 could overcome resistance to almonertinib, as evidenced in both in vitro and in vivo models. Our findings suggest that LCN2 was a crucial regulator for conferring almonertinib-resistance in NSCLC and demonstrate the potential utility of targeting the LCN2-MMP-9 axis for clinical treatment of almonertinib-resistant lung adenocarcinoma.

Resistance of HNSCC cell models to pan-FGFR inhibition depends on the EMT phenotype associating with clinical outcome

BACKGROUND

Focal adhesion signaling involving receptor tyrosine kinases (RTK) and integrins co-controls cancer cell survival and therapy resistance. However, co-dependencies between these receptors and therapeutically exploitable vulnerabilities remain largely elusive in HPV-negative head and neck squamous cell carcinoma (HNSCC).

METHODS

The cytotoxic and radiochemosensitizing potential of targeting 10 RTK and β1 integrin was determined in up to 20 3D matrix-grown HNSCC cell models followed by drug screening and patient-derived organoid validation. RNA sequencing and protein-based biochemical assays were performed for molecular characterization. Bioinformatically identified transcriptomic signatures were applied to patient cohorts.

RESULTS

Fibroblast growth factor receptor (FGFR 1-4) targeting exhibited the strongest cytotoxic and radiosensitizing effects as monotherapy and combined with β1 integrin inhibition, exceeding the efficacy of the other RTK studied. Pharmacological pan-FGFR inhibition elicited responses ranging from cytotoxicity/radiochemosensitization to resistance/radiation protection. RNA sequence analysis revealed a mesenchymal-to-epithelial transition (MET) in sensitive cell models, whereas resistant cell models exhibited a partial epithelial-to-mesenchymal transition (EMT). Accordingly, inhibition of EMT-associated kinases such as EGFR caused reduced adaptive resistance and enhanced (radio)sensitization to FGFR inhibition cell model- and organoid-dependently. Transferring the EMT-associated transcriptomic profiles to HNSCC patient cohorts not only demonstrated their prognostic value but also provided a conclusive validation of the presence of EGFR-related vulnerabilities that can be strategically exploited for therapeutic interventions.

CONCLUSIONS

This study demonstrates that pan-FGFR inhibition elicits a beneficial radiochemosensitizing and a detrimental radioprotective potential in HNSCC cell models. Adaptive EMT-associated resistance appears to be of clinical importance, and we provide effective molecular approaches to exploit this therapeutically.

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.

Genomic alterations in oral multiple primary cancers

Oral squamous cell carcinoma (OSCC) is the predominant type of oral cancer, while some patients may develop oral multiple primary cancers (MPCs) with unclear etiology. This study aimed to investigate the clinicopathological characteristics and genomic alterations of oral MPCs. Clinicopathological data from patients with oral single primary carcinoma (SPC, n = 202) and oral MPCs (n = 34) were collected and compared. Copy number alteration (CNA) analysis was conducted to identify chromosomal-instability differences among oral MPCs, recurrent OSCC cases, and OSCC patients with lymph node metastasis. Whole-exome sequencing was employed to identify potential unique gene mutations in oral MPCs patients. Additionally, CNA and phylogenetic tree analyses were used to gain preliminary insights into the molecular characteristics of different primary tumors within individual patients. Our findings revealed that, in contrast to oral SPC, females predominated the oral MPCs (70.59%), while smoking and alcohol use were not frequent in MPCs. Moreover, long-term survival outcomes were poorer in oral MPCs. From a CNA perspective, no significant differences were observed between oral MPCs patients and those with recurrence and lymph node metastasis. In addition to commonly mutated genes such as CASP8, TP53 and MUC16, in oral MPCs we also detected relatively rare mutations, such as HS3ST6 and RFPL4A. Furthermore, this study also demonstrated that most MPCs patients exhibited similarities in certain genomic regions within individuals, and distinct differences of the similarity degree were observed between synchronous and metachronous oral MPCs.

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. Nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). Here, we profile contributions of UHRF1 and DNMT1 to genome-wide DNA methylation inheritance and dissect specific roles for ubiquitin signaling in this process. We reveal DNA methylation maintenance at low-density CpGs is vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMD), a methylation signature observed across human cancers. Furthermore, disrupting DNMT1 UIM2 function abolishes DNA methylation maintenance. Collectively, we show DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to the development of PMDs in human cancers.

Parkinson's disease risk enhancers in microglia

Genome-wide association studies have identified thousands of single nucleotide polymorphisms that associate with increased risk for Parkinson's disease (PD), but the functions of most of them are unknown. Using assay for transposase-accessible chromatin (ATAC) and H3K27ac chromatin immunoprecipitation (ChIP) sequencing data, we identified 73 regulatory elements in microglia that overlap PD risk SNPs. To determine the target genes of a "risk enhancer" within intron two of SNCA, we used CRISPR-Cas9 to delete the open chromatin region where two PD risk SNPs reside. The loss of the enhancer led to reduced expression of multiple genes including SNCA and the adjacent gene MMRN1. It also led to expression changes of genes involved in glucose metabolism, a process that is known to be altered in PD patients. Our work expands the role of SNCA in PD and provides a connection between PD-associated genetic variants and underlying biology that points to a risk mechanism in microglia.

PCGF6 controls murine Tuft cell differentiation via H3K9me2 modification independently of Polycomb repression

Cell fate is determined by specific transcription programs that are essential for tissue homeostasis and regeneration. The E3-ligases RING1A and B represent the core activity of the Polycomb repressive complex 1 (PRC1) that deposits repressive histone H2AK119 mono-ubiquitination (H2AK119ub1), which is essential for mouse intestinal homeostasis by preserving stem cell functions. However, the specific role of different PRC1 forms, which are defined by the six distinct PCGF1-6 paralogs, remains largely unexplored in vivo. We report that PCGF6 regulates mouse intestinal Tuft cell differentiation independently of H2AK119ub1 deposition. We show that PCGF6 chromatin occupancy expands outside Polycomb repressive domains, associating with unique promoter and distal regulatory elements. This occurs in the absence of RING1A/B and involves MGA-mediated E-BOX recognition and specific H3K9me2 promoter deposition. PCGF6 inactivation induces an epithelial autonomous accumulation of Tuft cells that was not phenocopied by RING1A/B loss. This involves direct PCGF6 association with a Tuft cell differentiation program that identified Polycomb-independent properties of PCGF6 in adult tissues homeostasis.

Genomic vulnerability of a freshwater salmonid under climate change

Understanding the adaptive potential of populations and species is pivotal for minimizing the loss of biodiversity in this era of rapid climate change. Adaptive potential has been estimated in various ways, including based on levels of standing genetic variation, presence of potentially beneficial alleles, and/or the severity of environmental change. Kokanee salmon, the non-migratory ecotype of sockeye salmon (Oncorhynchus nerka), is culturally and economically important and has already been impacted by the effects of climate change. To assess its climate vulnerability moving forward, we integrated analyses of standing genetic variation, genotype-environment associations, and climate modeling based on sequence and structural genomic variation from 224 whole genomes sampled from 22 lakes in British Columbia and Yukon (Canada). We found that variables for extreme temperatures, particularly warmer temperatures, had the most pervasive signature of selection in the genome and were the strongest predictors of levels of standing variation and of putatively adaptive genomic variation, both sequence and structural. Genomic offset estimates, a measure of climate vulnerability, were significantly correlated with higher increases in extreme warm temperatures, further highlighting the risk of summer heat waves that are predicted to increase in frequency in the future. Levels of standing genetic variation, an important metric for population viability and resilience, were not correlated with genomic offset. Nonetheless, our combined approach highlights the importance of integrating different sources of information and genomic data to formulate more comprehensive and accurate predictions on the vulnerability of populations and species to future climate change.

Genetic Risk Factors for Early-Onset Merkel Cell Carcinoma

Importance

Merkel cell carcinoma (MCC) is a rare, aggressive neuroendocrine skin cancer. Of the patients who develop MCC annually, only 4% are younger than 50 years.

Objective

To identify genetic risk factors for early-onset MCC via genomic sequencing.

Design, Setting, and Participants

The study represents a multicenter collaboration between the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS), the National Institute of Allergy and Infectious Diseases (NIAID), and the University of Washington. Participants with early-onset and later-onset MCC were prospectively enrolled in an institutional review board-approved study at the University of Washington between January 2003 and May 2019. Unrelated controls were enrolled in the NIAID Centralized Sequencing Program (CSP) between September 2017 and September 2021. Analysis was performed from September 2021 and March 2023. Early-onset MCC was defined as disease occurrence in individuals younger than 50 years. Later-onset MCC was defined as disease occurrence at age 50 years or older. Unrelated controls were evaluated by the NIAID CSP for reasons other than familial cancer syndromes, including immunological, neurological, and psychiatric disorders.

Results

This case-control analysis included 1012 participants: 37 with early-onset MCC, 45 with later-onset MCC, and 930 unrelated controls. Among 37 patients with early-onset MCC, 7 (19%) had well-described variants in genes associated with cancer predisposition. Six patients had variants associated with hereditary cancer syndromes (ATM = 2, BRCA1 = 2, BRCA2 = 1, and TP53 = 1) and 1 patient had a variant associated with immunodeficiency and lymphoma (MAGT1). Compared with 930 unrelated controls, the early-onset MCC cohort was significantly enriched for cancer-predisposing pathogenic or likely pathogenic variants in these 5 genes (odds ratio, 30.35; 95% CI, 8.89-106.30; P < .001). No germline disease variants in these genes were identified in 45 patients with later-onset MCC. Additional variants in DNA repair genes were also identified among patients with MCC.

Conclusions and Relevance

Because variants in certain DNA repair and cancer predisposition genes are associated with early-onset MCC, genetic counseling and testing should be considered for patients presenting at younger than 50 years.

Chromosome-level genome assemblies of Nicotiana tabacum, Nicotiana sylvestris, and Nicotiana tomentosiformis

The Solanaceae species Nicotiana tabacum, an economically important crop plant cultivated worldwide, is an allotetraploid species that appeared about 200,000 years ago as the result of the hybridization of diploid ancestors of Nicotiana sylvestris and Nicotiana tomentosiformis. The previously published genome assemblies for these three species relied primarily on short-reads, and the obtained pseudochromosomes only partially covered the genomes. In this study, we generated annotated de novo chromosome-level genomes of N. tabacum, N. sylvestris, and N. tomentosiformis, which contain 3.99 Gb, 2.32 Gb, and 1.74 Gb, respectively of sequence data, with 97.6%, 99.5%, and 95.9% aligned in chromosomes, and represent 99.2%, 98.3%, and 98.5% of the near-universal single-copy orthologs Solanaceae genes. The completion levels of these chromosome-level genomes for N. tabacum, N. sylvestris, and N. tomentosiformis are comparable to other reference Solanaceae genomes, enabling more efficient synteny-based cross-species research.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases

PURPOSE

Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION).

METHODS

PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation.

RESULTS

A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants.

CONCLUSION

GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

Multi-omic dataset of patient-derived tumor organoids of neuroendocrine neoplasms

BACKGROUND

Organoids are 3-dimensional experimental models that summarize the anatomical and functional structure of an organ. Although a promising experimental model for precision medicine, patient-derived tumor organoids (PDTOs) have currently been developed only for a fraction of tumor types.

RESULTS

We have generated the first multi-omic dataset (whole-genome sequencing [WGS] and RNA-sequencing [RNA-seq]) of PDTOs from the rare and understudied pulmonary neuroendocrine tumors (n = 12; 6 grade 1, 6 grade 2) and provide data from other rare neuroendocrine neoplasms: small intestine (ileal) neuroendocrine tumors (n = 6; 2 grade 1 and 4 grade 2) and large-cell neuroendocrine carcinoma (n = 5; 1 pancreatic and 4 pulmonary). This dataset includes a matched sample from the parental sample (primary tumor or metastasis) for a majority of samples (21/23) and longitudinal sampling of the PDTOs (1 to 2 time points), for a total of n = 47 RNA-seq and n = 33 WGS. We here provide quality control for each technique and the raw and processed data as well as all scripts for genomic analyses to ensure an optimal reuse of the data. In addition, we report gene expression data and somatic small variant calls and describe how they were generated, in particular how we used WGS somatic calls to train a random forest classifier to detect variants in tumor-only RNA-seq. We also report all histopathological images used for medical diagnosis: hematoxylin and eosin-stained slides, brightfield images, and immunohistochemistry images of protein markers of clinical relevance.

CONCLUSIONS

This dataset will be critical to future studies relying on this PDTO biobank, such as drug screens for novel therapies and experiments investigating the mechanisms of carcinogenesis in these understudied diseases.

Fragmentstein-facilitating data reuse for cell-free DNA fragment analysis

SUMMARY

Method development for the analysis of cell-free DNA (cfDNA) sequencing data is impeded by limited data sharing due to the strict control of sensitive genomic data. An existing solution for facilitating data sharing removes nucleotide-level information from raw cfDNA sequencing data, keeping alignment coordinates only. This simplified format can be publicly shared and would, theoretically, suffice for common functional analyses of cfDNA data. However, current bioinformatics software requires nucleotide-level information and cannot process the simplified format. We present Fragmentstein, a command-line tool for converting non-sensitive cfDNA-fragmentation data into alignment mapping (BAM) files. Fragmentstein complements fragment coordinates with sequence information from a reference genome to reconstruct BAM files. We demonstrate the utility of Fragmentstein by showing the feasibility of copy number variant (CNV), nucleosome occupancy, and fragment length analyses from non-sensitive fragmentation data.

AVAILABILITY AND IMPLEMENTATION

Implemented in bash, Fragmentstein is available at https://github.com/uzh-dqbm-cmi/fragmentstein, licensed under GNU GPLv3.

Whole‑exome sequencing reveals Lewis lung carcinoma is a hypermutated Kras/Nras-mutant cancer with extensive regional mutation clusters in its genome

Lewis lung carcinoma (LLC), as a widely used preclinical cancer model, has still not been genetically and genomically characterized. Here, we performed a whole-exome sequencing analysis on the LLC cell line to elucidate its molecular characteristics and etiologies. Our data showed that LLC originated from a male mouse belonging to C57BL/6L (a transitional strain between C57BL/6J and C57BL/6N) and contains substantial somatic SNV and InDel mutations (> 20,000). Extensive regional mutation clusters are present in its genome, which were caused mainly by the mutational processes underlying the SBS1, SBS5, SBS15, SBS17a, and SBS21 signatures during frequent structural rearrangements. Thirty three deleterious mutations are present in 30 cancer genes including Kras, Nras, Trp53, Dcc, and Cacna1d. Cdkn2a and Cdkn2b are biallelically deleted from the genome. Five pathways (RTK/RAS, p53, cell cycle, TGFB, and Hippo) are oncogenically deregulated or affected. The major mutational processes in LLC include chromosomal instability, exposure to metabolic mutagens, spontaneous 5-methylcytosine deamination, defective DNA mismatch repair, and reactive oxygen species. Our data also suggest that LLC is a lung cancer similar to human lung adenocarcinoma. This study lays a molecular basis for the more targeted application of LLC in preclinical research.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Familial gigantiform cementoma with recurrent ANO5 p.Cys356Tyr mutations: Clinicopathological and genetic study with literature review

BACKGROUND

Familial gigantiform cementoma (FGC) is a rare tumor characterized by the early onset of multi-quadrant fibro-osseous lesions in the jaws, causing severe maxillofacial deformities. Its clinicopathological features overlap with those of other benign fibro-osseous lesions. FGC eventually exhibits progressively rapid growth, but no suspected causative gene has been identified.

METHODS

In this study, three patients with FGC were recruited, and genomic DNA from the tumor tissue and peripheral blood was extracted for whole-exome sequencing.

RESULTS

Results showed that all three patients harbored the heterozygous mutation c.1067G > A (p.Cys356Tyr) in the ANO5 gene. Furthermore, autosomal dominant mutations in ANO5 at this locus have been identified in patients with gnathodiaphyseal dysplasia (GDD) and are considered a potential causative agent, suggesting a genetic association between FGC and GDD. In addition, multifocal fibrous bone lesions with similar clinical presentations were detected, including five cases of florid cemento-osseous dysplasia, five cases of polyostotic fibrous dysplasia, and eight cases of juvenile ossifying fibromas; however, none of them harbored mutations in the ANO5 gene.

CONCLUSION

Our findings indicate that FGC may be an atypical variant of GDD, providing evidence for the feasibility of ANO5 gene testing as an auxiliary diagnostic method for complex cases with multiple quadrants.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Expanding Genotype-Phenotype Correlation of CLCNKA and CLCNKB Variants Linked to Hearing Loss

The ClC-K channels CLCNKA and CLCNKB are crucial for the transepithelial transport processes required for sufficient urinary concentrations and sensory mechanoelectrical transduction in the cochlea. Loss-of-function alleles in these channels are associated with various clinical phenotypes, ranging from hypokalemic alkalosis to sensorineural hearing loss (SNHL) accompanied by severe renal conditions, i.e., Bartter's syndrome. Using a stepwise genetic approach encompassing whole-genome sequencing (WGS), we identified one family with compound heterozygous variants in the ClC-K channels, specifically a truncating variant in CLCNKA in trans with a contiguous deletion of CLCNKA and CLCNKB. Breakpoint PCR and Sanger sequencing elucidated the breakpoint junctions derived from WGS, and allele-specific droplet digital PCR confirmed one copy loss of the CLCNKA_CLCNKB contiguous deletion. The proband that harbors the CLCNKA_CLCNKB variants is characterized by SNHL without hypokalemic alkalosis and renal anomalies, suggesting a distinct phenotype in the ClC-K channels in whom SNHL predominantly occurs. These results expanded genotypes and phenotypes associated with ClC-K channels, including the disease entities associated with non-syndromic hearing loss. Repeated identification of deletions across various extents of CLCNKA_CLCNKB suggests a mutational hotspot allele, highlighting the need for an in-depth analysis of the CLCNKA_CLCNKB intergenic region, especially in undiagnosed SNHL patients with a single hit in CLCNKA.

Dupsifter: a lightweight duplicate marking tool for whole genome bisulfite sequencing

SUMMARY

In whole genome sequencing data, polymerase chain reaction amplification results in duplicate DNA fragments coming from the same location in the genome. The process of preparing a whole genome bisulfite sequencing (WGBS) library, on the other hand, can create two DNA fragments from the same location that should not be considered duplicates. Currently, only one WGBS-aware duplicate marking tool exists. However, it only works with the output from a single tool, does not accept streaming input or output, and requires a substantial amount of memory relative to the input size. Dupsifter provides an aligner-agnostic duplicate marking tool that is lightweight, has streaming capabilities, and is memory efficient.

AVAILABILITY AND IMPLEMENTATION

Source code and binaries are freely available at https://github.com/huishenlab/dupsifter under the MIT license. Dupsifter is implemented in C and is supported on macOS and Linux.

β1 integrin signaling governs necroptosis via the chromatin-remodeling factor CHD4

Fibrosis, characterized by sustained activation of myofibroblasts and excessive extracellular matrix (ECM) deposition, is known to be associated with chronic inflammation. Receptor-interacting protein kinase 3 (RIPK3), the central kinase of necroptosis signaling, is upregulated in fibrosis and contributes to tumor necrosis factor (TNF)-mediated inflammation. In bile-duct-ligation-induced liver fibrosis, we found that myofibroblasts are the major cell type expressing RIPK3. Genetic ablation of β1 integrin, the major profibrotic ECM receptor in fibroblasts, not only abolished ECM fibrillogenesis but also blunted RIPK3 expression via a mechanism mediated by the chromatin-remodeling factor chromodomain helicase DNA-binding protein 4 (CHD4). While the function of CHD4 has been conventionally linked to the nucleosome-remodeling deacetylase (NuRD) and CHD4-ADNP-HP1(ChAHP) complexes, we found that CHD4 potently repressed a set of genes, including Ripk3, with high locus specificity but independent of either the NuRD or the ChAHP complex. Thus, our data uncover that β1 integrin intrinsically links fibrotic signaling to RIPK3-driven inflammation via a novel mode of action of CHD4.

Hi-C metagenome sequencing reveals soil phage-host interactions

Bacteriophages are abundant in soils. However, the majority are uncharacterized, and their hosts are unknown. Here, we apply high-throughput chromosome conformation capture (Hi-C) to directly capture phage-host relationships. Some hosts have high centralities in bacterial community co-occurrence networks, suggesting phage infections have an important impact on the soil bacterial community interactions. We observe increased average viral copies per host (VPH) and decreased viral transcriptional activity following a two-week soil-drying incubation, indicating an increase in lysogenic infections. Soil drying also alters the observed phage host range. A significant negative correlation between VPH and host abundance prior to drying indicates more lytic infections result in more host death and inversely influence host abundance. This study provides empirical evidence of phage-mediated bacterial population dynamics in soil by directly capturing specific phage-host interactions.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Allelic hierarchy for USH2A influences auditory and visual phenotypes in South Korean patients

When medical genetic syndromes are influenced by allelic hierarchies, mutant alleles have distinct effects on clinical phenotypes. Genotype-phenotype correlations for Usher syndrome type 2 (USH2) suggest that the USH2A gene exhibits an allelic hierarchy. Here, we analyzed the phenotypes and genotypes of 16 South Korean patients with USH2A biallelic variants to investigate an allelic hierarchy from audiological and ophthalmological perspectives. Using whole exome and genome sequencing, 18 mutant alleles, including 4 novel alleles, were identified and implicated in USH2A-related disorders. Truncated alleles were linked to earlier onset of subjective hearing loss and more severe thresholds; biallelic truncated alleles had more severe effects. Truncated alleles were also associated with retinal structure degeneration and severe functional deterioration. However, younger patients (aged < 16 years) did not exhibit overt retinitis pigmentosa even when they had biallelic truncated alleles, suggesting that USH2A-related USH2 can mimic nonsyndromic hearing loss. For truncated alleles, there was a clear correlation between mean hearing threshold and 30-Hz flicker electroretinography implicit time. This study provides the first evidence of an USH2A-related allelic hierarchy among South Korean patients; our data yield valuable insights concerning the natural courses of clinical phenotypes and how genotype-based therapies may be used.

Vaginal microbes alter epithelial transcriptomic and epigenomic modifications providing insight into the molecular mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health with microbial communities dominated by Lactobacillus spp. being considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes including Gardnerella vaginalis , have been associated with adverse reproductive outcomes. However, the molecular pathways modulated by microbe interactions with the cervicovaginal epithelia remain unclear. Using RNA-sequencing, we characterize the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulated genes were associated with an activated innate immune response including anti-microbial peptides and inflammasome pathways, represented by NLRP3-mediated increases in caspase-1, IL-1β and cell death. Cervicovaginal epithelial cells exposed to L. crispatus showed limited transcriptomic changes, while exposure to L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells. ATAC-sequencing confirmed epigenetic changes with reduced chromatin accessibility. This study reveals new insight into host-microbe interactions in the lower reproductive tract and suggest potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Variants of IFNL4 Gene in Amazonian and Northern Brazilian Populations

Since the discovery of the polymorphic nature of the IFNL4 gene, its variants have been investigated and associated with several viral diseases, with an emphasis on hepatitis C. However, the impacts of these variants on mixed-race and native populations in the northern region of Brazil are scarce. We investigated three variants of the IFNL4 gene in populations from this location, which were among the 14 most frequent variants in worldwide populations, and compared the frequencies obtained to populational data from the 1000 Genomes Project, gnomAD and ABraOM databases. Our results demonstrate that mixed-race and native populations from the northern region of Brazil present frequencies like those of European and Asian groups for the rs74597329 and rs11322783 variants, and like all populations presented for the rs4803221 variant. These data reinforce the role of world populations in shaping the genetic profile of Brazilian populations, indicate patterns of illness according to the expressed genotype, and infer an individual predisposition to certain diseases.

A reference assembly for the legume cover crop hairy vetch (Vicia villosa)

Vicia villosa is an incompletely domesticated annual legume of the Fabaceae family native to Europe and Western Asia. V. villosa is widely used as a cover crop and forage due to its ability to withstand harsh winters. Here, we generated a reference-quality genome assembly (Vvill1.0) from low error-rate long-sequence reads to improve the genetic-based trait selection of this species. Our Vvill1.0 assembly includes seven scaffolds corresponding to the seven estimated linkage groups and comprising approximately 68% of the total genome size of 2.03 Gbp. This assembly is expected to be a useful resource for genetically improving this emerging cover crop species and provide useful insights into legume genomics and plant genome evolution.

ATF4 renders human T-cell acute lymphoblastic leukemia cell resistance to FGFR1 inhibitors through amino acid metabolic reprogramming

Abnormalities of FGFR1 have been reported in multiple malignancies, suggesting FGFR1 as a potential target for precision treatment, but drug resistance remains a formidable obstacle. In this study, we explored whether FGFR1 acted a therapeutic target in human T-cell acute lymphoblastic leukemia (T-ALL) and the molecular mechanisms underlying T-ALL cell resistance to FGFR1 inhibitors. We showed that FGFR1 was significantly upregulated in human T-ALL and inversely correlated with the prognosis of patients. Knockdown of FGFR1 suppressed T-ALL growth and progression both in vitro and in vivo. However, the T-ALL cells were resistant to FGFR1 inhibitors AZD4547 and PD-166866 even though FGFR1 signaling was specifically inhibited in the early stage. Mechanistically, we found that FGFR1 inhibitors markedly increased the expression of ATF4, which was a major initiator for T-ALL resistance to FGFR1 inhibitors. We further revealed that FGFR1 inhibitors induced expression of ATF4 through enhancing chromatin accessibility combined with translational activation via the GCN2-eIF2α pathway. Subsequently, ATF4 remodeled the amino acid metabolism by stimulating the expression of multiple metabolic genes ASNS, ASS1, PHGDH and SLC1A5, maintaining the activation of mTORC1, which contributed to the drug resistance in T-ALL cells. Targeting FGFR1 and mTOR exhibited synergistically anti-leukemic efficacy. These results reveal that FGFR1 is a potential therapeutic target in human T-ALL, and ATF4-mediated amino acid metabolic reprogramming contributes to the FGFR1 inhibitor resistance. Synergistically inhibiting FGFR1 and mTOR can overcome this obstacle in T-ALL therapy.

Comparative analysis of genome-wide copy number variations between Tibetan sheep and White Suffolk sheep

The DNA copy number variations (CNVs) are widely involved in affecting various kinds of biological functions, such as environmental adaptation. Tibetan sheep and White Suffolk sheep are two representative indigenous and exotic breeds raised in Sichuan, China, and both of them have many contrasting biological characteristics. In this study, we employed high-throughput sequencing approach to investigate genome-wide CNVs between the two sheep breeds. A total of 11,135 CNV regions (CNVRs) consisting of 6,488 deletions and 4,647 duplications were detected, whose length ranged from 1,599 bp to 0.56 Mb with the mean of 4,658 bp. There were 281 CNVRs segregated between Tibetan sheep and White Suffolk sheep, and 18 of them have been fixed within both breeds. Functional analyses of candidate genes within the segregating CNVRs revealed the thyroid hormone signaling pathway and CTNNB1 gene that would be responsible for differential biological characteristics of breeds, such as energy metabolism, seasonal reproduction, and litter size. Furthermore, the segregating CNVRs identified in this study were overlapped with many known quantitative trait loci that are associated with growth, testis weight, and reproductive seasonality. In conclusion, these results help us better understanding differential biological characteristics between Tibetan sheep and White Suffolk sheep.

Leishmania genetic exchange is mediated by IgM natural antibodies

Host factors that mediate Leishmania genetic exchange are not well defined. Here we demonstrate that natural IgM (IgMn)1-4 antibodies mediate parasite genetic exchange by inducing the transient formation of a spherical parasite clump that promotes parasite fusion and hybrid formation. We establish that IgMn from Leishmania-free animals binds to the surface of Leishmania parasites to induce significant changes in the expression of parasite transcripts and proteins. Leishmania binding to IgMn is partially lost after glycosidase treatment, although parasite surface phosphoglycans, including lipophosphoglycan, are not required for IgMn-induced parasite clumping. Notably, the transient formation of parasite clumps is essential for Leishmania hybridization in vitro. In vivo, we observed a 12-fold increase in hybrid formation in sand flies provided a second blood meal containing IgMn compared with controls. Furthermore, the generation of recombinant progeny from mating hybrids and parental lines were only observed in sand flies provided with IgMn. Both in vitro and in vivo IgM-induced Leishmania crosses resulted in full genome hybrids that show equal patterns of biparental contribution. Leishmania co-option of a host natural antibody to facilitate mating in the insect vector establishes a new paradigm of parasite-host-vector interdependence that contributes to parasite diversity and fitness by promoting genetic exchange.

Structural variants and short tandem repeats impact gene expression and splicing in bovine testis tissue

Structural variants (SVs) and short tandem repeats (STRs) are significant sources of genetic variation. However, the impacts of these variants on gene regulation have not been investigated in cattle. Here, we genotyped and characterized 19,408 SVs and 374,821 STRs in 183 bovine genomes and investigated their impact on molecular phenotypes derived from testis transcriptomes. We found that 71% STRs were multiallelic. The vast majority (95%) of STRs and SVs were in intergenic and intronic regions. Only 37% SVs and 40% STRs were in high linkage disequilibrium (LD) (R2 > 0.8) with surrounding SNPs/insertions and deletions (Indels), indicating that SNP-based association testing and genomic prediction are blind to a nonnegligible portion of genetic variation. We showed that both SVs and STRs were more than 2-fold enriched among expression and splicing QTL (e/sQTL) relative to SNPs/Indels and were often associated with differential expression and splicing of multiple genes. Deletions and duplications had larger impacts on splicing and expression than any other type of SV. Exonic duplications predominantly increased gene expression either through alternative splicing or other mechanisms, whereas expression- and splicing-associated STRs primarily resided in intronic regions and exhibited bimodal effects on the molecular phenotypes investigated. Most e/sQTL resided within 100 kb of the affected genes or splicing junctions. We pinpoint candidate causal STRs and SVs associated with the expression of SLC13A4 and TTC7B and alternative splicing of a lncRNA and CAPP1. We provide a catalog of STRs and SVs for taurine cattle and show that these variants contribute substantially to gene expression and splicing variation.