Variable Expression Research Articles

MiRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs.

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE.

Levodopa-responsive dystonia, parkinsonism, and treatment-resistant schizoaffective disorder in Williams syndrome.

Williams syndrome (WS; chromosome 7q11.23 deletion) is a rare, multisystemic, neurodevelopmental disorder with variable penetrance and expressivity. Although movement and psychiatric disorders are known to occur in individuals with WS, parkinsonism, dystonia, and treatment-resistant schizoaffective disorder have not been formally described. We present two unrelated cases of adults with molecularly confirmed WS and typical histories of developmental delays, intellectual/learning disabilities, and treatment-responsive anxiety/mood disorder who developed similar noteworthy neuropsychiatric expressions. We reviewed detailed neuropsychiatric histories, laboratory investigations, neuroimaging, and treatment responses and compared data for the two cases. Both individuals developed treatment-resistant schizoaffective disorder in adulthood requiring multiple trials of antipsychotic treatments. While on clozapine, both patients developed parkinsonism and generalized dystonia with truncal involvement that responded to trials of low-dose levodopa without exacerbating underlying psychotic or affective symptoms. This report illustrates the novel occurrence of levodopa-responsive movement disorders and treatment-resistant schizoaffective disorder in individuals with WS, adding to the expanding neuropsychiatric phenotypes, and highlighting potential shared underlying mechanisms. The observed treatment response suggests that levodopa, in relatively low doses, may be safe and useful in ameliorating presumed antipsychotic-associated parkinsonism and tardive dystonia in WS.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BackgroundInteractions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.ResultsWe comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.ConclusionsOur study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.

Dynamic behavior of the orthotropic hollow cylinder under a magnetic field: a thermoelastic analysis

The present article focuses on the complete analytical solution of an orthotropic hollow cylinder using a proposed magneto-thermoelastic model and finite Hankel transformation. The study considers the Lord–Shulman and Green–Lindsay models of hyperbolic generalized thermoelasticity within the context of thermoelasticity theory. The analytical expressions of the field variables of the cylinder in a magnetic field were obtained under thermo-mechanical loadings. The study also offered closed-form field functions and measured the thermoelastic behavior concerning the magnetic field intensity of the cylinder based on thermoelastic theories involving a second sound. The inclusion of hyperbolic-type heat conduction models enhanced the visualization of the reflection of thermal waves into the cylinder. The study also investigated the impacts of orthotropic material properties on the considered medium and compared the results with the existing literature.

Multiorgan manifestations of COL4A1 and COL4A2 variants and proposal for a clinical management protocol.

COL4A1/2 variants are associated with highly variable multiorgan manifestations. Depicting the whole clinical spectrum of COL4A1/2-related manifestations is challenging, and there is no consensus on management and preventative strategies. Based on a systematic review of current evidence on COL4A1/2-related disease, we developed a clinical questionnaire that we administered to 43 individuals from 23 distinct families carrying pathogenic variants. In this cohort, we extended ophthalmological and cardiological examinations to asymptomatic individuals and those with only limited or mild, often nonspecific, clinical signs commonly occurring in the general population (i.e., oligosymptomatic). The most frequent clinical findings emerging from both the literature review and the questionnaire included stroke (203/685, 29.6%), seizures or epilepsy (199/685, 29.0%), intellectual disability or developmental delay (168/685, 24.5%), porencephaly/schizencephaly (168/685, 24.5%), motor impairment (162/685, 23.6%), cataract (124/685, 18.1%), hematuria (63/685, 9.2%), and retinal arterial tortuosity (58/685, 8.5%). In oligosymptomatic and asymptomatic carriers, ophthalmological investigations detected retinal vascular tortuosity (5/13, 38.5%), dysgenesis of the anterior segment (4/13, 30.8%), and cataract (2/13, 15.4%), while cardiological investigations were unremarkable except for mild ascending aortic ectasia in 1/8 (12.5%). Our multimodal approach confirms highly variable penetrance and expressivity in COL4A1/2-related conditions, even at the intrafamilial level with neurological involvement being the most frequent and severe finding in both children and adults. We propose a protocol for prevention and management based on individualized risk estimation and periodic multiorgan evaluations.

Majority of the highly variable NLRs in maize share genomic location and contain additional target-binding domains.

Nucleotide-binding, Leucine Rich Repeat proteins (NLRs) are a major class of immune receptors in plants. NLRs include both conserved and rapidly evolving members, however their evolutionary trajectory in crops remains understudied. Availability of crop pan-genomes enables analysis of the recent events in the evolution of this highly complex gene family within domesticated species. Here, we investigated the NLR complement of 26 nested association mapping (NAM) founder lines of maize. We found that maize has just four main subfamilies containing rapidly evolving highly variable NLR (hvNLR) receptors. Curiously, three of these phylogenetically distinct hvNLR lineages are located in adjacent clusters on chromosome 10. Members of the same hvNLR clade show variable expression and methylation across lines and tissues, consistent with their rapid evolution. By combining sequence diversity analysis and AlphaFold2 computational structure prediction we predicted ligand binding sites in the hvNLRs. We also observed novel insertion domains in the LRR regions of two hvNLR subfamilies that likely contribute to target recogniton. To make this analysis accessible, we created NLRCladeFinder, a Google Colaboratory notebook, that accepts any newly identified NLR sequence, places it in the evolutionary context of the maize pan-NLRome, and provides an updated clade alignment, phylogenetic tree, and sequence diversity information for the gene of interest.

Undifferentiated Round Cell Sarcoma With CRTC1::SS18 Fusion: Expanding Clinicopathologic Features of a Rare Translocation Sarcoma With Prominent Desmoplastic Stroma

Undifferentiated round cell sarcomas (URCS) represent a diverse group of tumors, including conventional Ewing sarcoma, round cell sarcoma with EWSR1/FUS–non-ETS fusions, CIC-rearranged sarcoma, and sarcoma with BCOR alterations. Since 2018, 3 cases of URCS with a novel CRTC1::SS18 gene fusion have been reported in the literature. Herein, we report 3 additional cases of CRTC1::SS18 sarcoma, thereby doubling the number of described cases and expanding the clinicopathologic features of this rare translocation sarcoma. Together with the previously reported cases, we show that the male-to-female ratio is 1:2 with a median age of 34 years (range, 12-42 years). Tumors occurred primarily in intramuscular locations involving the lower extremity. Histologically, all tumors contained uniform round-to-epithelioid cells with a moderate amount of eosinophilic cytoplasm growing in sheets and nests with prominent desmoplastic stroma reminiscent of desmoplastic small round cell tumor. Immunohistochemical results were nonspecific, demonstrating variable expression of CD99 (patchy), ALK, GATA3, and cyclin D1. RNA sequencing revealed CRTC1::SS18 gene fusions in all cases, involving exons 1 to 2 of CRTC1 (the 5′ partner gene) on chromosome 19 and either exon 2 or exon 4 of SS18 (the 3′ partner gene) on chromosome 18. The clinical course was variable. Although 1 previously reported case demonstrated aggressive behavior with a fatal outcome, 2 others had a relatively indolent course with gradual growth for 6 to 7 years prior to resection. Two cases developed metastatic disease, including 1 case with bilateral lung metastasis and 1 with locoregional spread to a lymph node. By analyzing the clinicopathologic features, we aimed to improve recognition of this rare translocation sarcoma to better understand its biologic potential, optimize patient management, and expand the current classification of URCS.

Integrated Transcriptomic and Metabolomic Analysis Revealed Abscisic Acid-Induced Regulation of Monoterpene Biosynthesis in Grape Berries.

Monoterpenes are a class of volatile organic compounds that play crucial roles in imparting floral and fruity aromas to Muscat-type grapes. However, our understanding of the regulatory mechanisms underpinning monoterpene biosynthesis in grapes, particularly following abscisic acid (ABA) treatment, remains elusive. This study aimed to explore the impact of exogenous ABA on monoterpene biosynthesis in Ruiduhongyu grape berries by employing Headspace Solid-Phase Micro-Extraction Gas Chromatography-Mass Spectrometry (HS-SPME/GC-MS) analysis and transcriptome sequencing. The results suggested significant differences in total soluble solids (TSS), pH, and total acid content. ABA treatment resulted in a remarkable increase in endogenous ABA levels, with concentrations declining from veraison to ripening stages. ABA treatment notably enhanced monoterpene concentrations, particularly at the E_L37 and E_L38 stages, elevating the overall floral aroma of grape berries. According to the variable gene expression patterns across four developmental stages in response to ABA treatment, the E_L37 stage had the largest number of differential expressed genes (DEGs), which was correlated with a considerable change in free monoterpenes. Furthermore, functional annotation indicated that the DEGs were significantly enriched in primary and secondary metabolic pathways, underlining the relationship between ABA, sugar accumulation, and monoterpene biosynthesis. ABA treatment upregulated key genes involved in the methylerythritol phosphate (MEP) pathway, enhancing carbon allocation and subsequently impacting terpene synthesis. This study also identified transcription factors, including MYB and AP2/ERF families, potentially modulating monoterpene and aroma-related genes. Weighted gene co-expression network analysis (WGCNA) linked ABA-induced gene expression to monoterpene accumulation, highlighting specific modules enriched with genes associated with monoterpene biosynthesis; one of these modules (darkgreen) contained genes highly correlated with most monoterpenes, emphasizing the role of ABA in enhancing grape quality during berry maturation. Together, these findings provide valuable insights into the multifaceted effects of exogenous ABA on monoterpene compounds and grape berry flavor development, offering potential applications in viticulture and enology.

Reduction in ACE2 expression in peripheral blood mononuclear cells during COVID-19 – implications for post COVID-19 conditions

BackgroundSevere COVID-19 is uncommon, restricted to 19% of the total population. In response to the first virus wave (alpha variant of SARS-CoV-2), we investigated whether a biomarker indicated severity of disease and, in particular, if variable expression of angiotensin converting enzyme 2 (ACE2) in blood might clarify this difference in risk and of post COVID -19 conditions (PCC).MethodsThe IRB-approved study compared patients hospitalized with severe COVID-19 to healthy controls. Severe infection was defined requiring oxygen or increased oxygen need from baseline at admission with positive COVID-19 PCR. A single blood sample was obtained from patients within a day of admission. ACE2 RNA expression in blood cells was measured by an RT-PCR assay. Plasma ACE1 and ACE2 enzyme activities were quantified by fluorescent peptides. Plasma TIMP-1, PIIINP and MMP-9 antigens were quantified by ELISA. Data were entered into REDCap and analyzed using STATA v 14 and GraphPad Prism v 10.ResultsForty-eight patients and 72 healthy controls were recruited during the pandemic. ACE2 RNA expression in peripheral blood mononuclear cells (PBMC) was rarely detected acutely during severe COVID-19 but common in controls (OR for undetected ACE2: 12.4 [95% CI: 2.62-76.1]). ACE2 RNA expression in PBMC did not determine plasma ACE1 and ACE2 activity, suggesting alternative cell-signaling pathways. Markers of fibrosis (TIMP-1 and PIIINP) and vasculopathy (MMP-9) were additionally elevated. ACE2 RNA expression during severe COVID-19 often responded within hours to convalescent plasma. Analogous to oncogenesis, we speculate that potent, persistent, cryptic processes following COVID-19 (the renin-angiotensin system (RAS), fibrosis and vasculopathy) initiate or promote post-COVID-19 conditions (PCC) in susceptible individuals.ConclusionsThis work elucidates biological and temporal plausibility for ACE2, TIMP1, PIIINP and MMP-9 in the pathogenesis of PCC. Intersection of these independent systems is uncommon and may in part explain the rarity of PCC.

Rare Variant in MRC2 Associated With Familial Supraventricular Tachycardia and Wolff-Parkinson-White Syndrome.

Accessory pathways are a common cause of supraventricular tachycardia (SVT) and can lead to sudden cardiac death in otherwise healthy children and adults when associated with Wolff-Parkinson-White syndrome. The goal of this study was to identify genetic variants within a large family with structurally normal hearts affected by SVT and Wolff-Parkinson-White syndrome and determine causality of the gene deficit in a corresponding mouse model. Whole exome sequencing performed on 2 distant members of a 3-generation family in which multiple members were affected by SVT or Wolff-Parkinson-White pattern (preexcitation) on ECG identified MRC2 as a candidate gene. Serial electrocardiograms, intracardiac electrophysiology studies, echocardiography, optical mapping studies, and histology were performed on both Mrc2 mutant and WT (wild-type) mice. A rare HET (heterozygous) missense variant c.2969A>G;p.Glu990Gly (E990G) in MRC2 was identified as the leading candidate gene variant segregating with the cardiac phenotype following an autosomal-dominant Mendelian trait segregation pattern with variable expressivity. In vivo electrophysiology studies revealed reentrant SVT in E990G mice. Optical mapping studies in E990G mice demonstrated abnormal retrograde conduction, suggesting the presence of an accessory pathway. Histological analysis of E990G mouse hearts showed a disordered ECM (extracellular matrix) in the annulus fibrosus. Finally, Mrc2 knockdown in human cardiac fibroblasts enhanced accelerated cell migration. This study identified a rare nonsynonymous variant in the MRC2 gene in individuals with familial reentrant SVT, Wolff-Parkinson-White ECG pattern, and structurally normal hearts. Furthermore, Mrc2 knock-in mice revealed an increased incidence of reentrant SVT and bypass tract formation in the setting of preserved cardiac structure and function.

French National Protocol for diagnosis and care of facioscapulohumeral muscular dystrophy (FSHD).

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most common genetically inherited myopathies in adults. It is characterized by incomplete penetrance and variable expressivity. Typically, FSHD patients display asymmetric weakness of facial, scapular, and humeral muscles that may progress to other muscle groups, particularly the abdominal and lower limb muscles. Early-onset patients display more severe muscle weakness and atrophy, resulting in a higher frequency of associated skeletal abnormalities. In these patients, multisystem involvement, including respiratory, ocular, and auditory, is more frequent and severe and may include the central nervous system. Adult-onset FSHD patients may also display some degree of multisystem involvement which mainly remains subclinical. In 95% of cases, FSHD patients carry a pathogenic contraction of the D4Z4 repeat units (RUs) in the subtelomeric region of chromosome 4 (4q35), which leads to the expression of DUX4 retrogene, toxic for muscles (FSHD1). Five percent of patients display the same clinical phenotype in association with a mutation in the SMCHD1 gene located in chromosome 18, inducing epigenetic modifications of the 4q D4Z4 repeated region and expression of DUX4 retrogene. This review highlights the complexities and challenges of diagnosing and managing FSHD, underscoring the importance of standardized approaches for optimal patient outcomes. It emphasizes the critical role of multidisciplinary care in addressing the diverse manifestations of FSHD across different age groups, from skeletal abnormalities in early-onset cases to the often-subclinical multisystem involvement in adults. With no current cure, the focus on alleviating symptoms and slowing disease progression through coordinated care is paramount.

Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

Exome sequencing identifies novel genes underlying primary congenital glaucoma in the National Birth Defects Prevention Study.

Primary congenital glaucoma (PCG) affects approximately 1 in 10,000 live born infants in the United States (U.S.). PCG has a autosomal recessive inheritance pattern, and variable expressivity and reduced penetrance have been reported. Likely causal variants in the most commonly mutated gene, CYP1B1, are less prevalent in the U.S., suggesting that alternative genes may contribute to the condition. This study utilized exome sequencing to investigate the genetic architecture of PCG in the U.S. and to identify novel genes and variants. We studied 37 family trios where infants had PCG and were part of the National Birth Defects Prevention Study (births 1997-2011), a U.S. multicenter study of birth defects. Samples underwent exome sequencing and sequence reads were aligned to the human reference sample (NCBI build 37/hg19). Variant filtration was conducted under de novo and Mendelian inheritance models using GEMINI. Among candidate variants, CYP1B1 was most represented (five trios, 13.5%). Twelve probands (32%) had potentially pathogenic variants in other genes not previously linked to PCG but important in eye development and/or to underlie Mendelian conditions with potential phenotypic overlap (e.g., CRYBB2, RXRA, GLI2). Variation in the genes identified in this population-based study may help to further explain the genetics of PCG.

Expanding the Spectrum of NUTM1 -Rearranged Sarcoma : A Clinicopathologic and Molecular Genetic Study of 8 Cases.

Apart from the lethal midline carcinoma (NUT carcinoma), NUTM1 translocation has also been reported in mesenchymal tumors, but is exceedingly rare. Here, we describe a series of 8 NUTM1 -rearranged sarcomas to further characterize the clinicopathologic features of this emerging entity. This cohort included 2 males and 6 females with age ranging from 24 to 64 years (mean: 51y; median: 56y). Tumors occurred in the colon (2), abdomen (2), jejunum (1), esophagus (1), lung (1) and infraorbital region (1). At diagnosis, 6 patients presented with metastatic disease. Tumor size ranged from 1 to 10.5cm (mean: 6cm; median: 5.5 cm). Histologically, 4 tumors were composed of primitive small round cells to epithelioid cells intermixed with variable spindle cells, while 3 tumors consisted exclusively of small round cells to epithelioid cells and 1 tumor consisted predominantly of high-grade spindle cells. The neoplastic cells were arranged in solid sheets, nests, or intersecting fascicles. Mitotic activity ranged from 1 to 15/10 HPF (median: 5/10 HPF). Other features included rhabdoid phenotype (4/8), pronounced nuclear convolutions (2/8), prominent stromal hyalinization (2/8), focally myxoid stroma (1/8), foci of osteoclasts (1/8), and necrosis (1/8). By immunohistochemistry, all tumors showed diffuse and strong nuclear staining of NUT protein, with variable expression of pancytokeratin (AE1/AE3) (2/8), CK18 (1/8), CD99 (3/8), NKX2.2 (2/8), cyclin D1 (2/8), desmin (2/8), BCOR (2/8), S100 (1/8), TLE1 (1/8), and synaptophysin (1/8). Seven of 8 tumors demonstrated NUTM1 rearrangement by fluorescence in situ hybridization analysis. RNA-sequencing analysis identified MXD4::NUTM1 (3/7), MXI1::NUTM1 (3/7), and MGA::NUTM1 (1/7) fusions, respectively. DNA-based methylation profiling performed in 2 cases revealed distinct methylation cluster differing from those of NUT carcinoma and undifferentiated small round cell and spindle cell sarcomas. At follow-up (range: 4 to 24mo), 1 patient experienced recurrence at 8.5 months, 4 patients were alive with metastatic disease (5, 10, 11, and 24mo after diagnosis), 3 patients remained well with no signs of recurrence or metastasis (4, 6, and 12mo after diagnosis). Our study further demonstrated that NUTM1 -rearranged sarcoma had a broad range of clinicopathologic spectrum. NUT immunohistochemistry should be included in the diagnostic approach of monotonous undifferentiated small round, epithelioid to high-grade spindle cell malignancies that difficult to classify by conventional means. DNA-based methylation profiling might provide a promising tool in the epigenetic classification of undifferentiated sarcomas.

A parametric approach to flow loss with incidence and tip clearance variables for a compressor linear cascade

A study was conducted using theoretical and numerical calculations to analyze the flow losses of a compressor linear cascade called SMU37-RL02. Its main goal is to determine the parametric expressions for vortex structure weight loss at the 103% axial chord position on the quasi-S3 section of the cascade. A flow loss distribution band delineation extraction method is proposed to improve data acquisition efficiency. The method expands the database by single-parameter variable expressions, and the simulated annealing algorithm solves the full-condition optimized tip clearance range for the cascade. The results show that the dominant flow loss weight of the vortex structures is transferred from the passage vortex to the tip leakage vortex following the increased tip clearance. Making the tip clearance Ct = 0.19%H as the optimization objective has the most stable and best performance at the full incoming flow conditions compared to the cascade with an equidistant tip clearance. The maximum loss reduction is 13.16% relative to the original cascade loss at the i = +1.9° incoming flow condition. The optimal cascade achieves relatively stable optimization over the full range of incoming flow conditions.

Treacher Collins Syndrome (mandibulofacial dysostosis) – A case report

Background. Treacher-Collins Syndrome (TCS) (or mandibulofacial dysostosis) is a rare autosomal dominant genetic disorder involving 1st and 2nd branchial arches present with craniofacial deformities with variable expressivity. TCS commonly manifests as facial dysmorphism with varied presentations. In this case report, we demonstrate the typical CT (Computed Tomography) findings of the characteristic features. Case report. Here, we report an 8-year-old boy who presented for evaluation of uneven facial features. Upon a multidisciplinary evaluation, he was diagnosed with TCS, and the characteristic radiological features are described in this report. Conclusions. TCS has a variable degree of phenotypic expression with no related developmental delay or neurologic disorder. Proper counselling, planning, and surgical procedures are necessary for more promising results.

A prospective study of HER3 expression pre and post neoadjuvant therapy of different breast cancer subtypes: implications for HER3 imaging therapy guidance

PurposeHER3, a member of the EGFR receptor family, plays a central role in driving oncogenic cell proliferation in breast cancer. Novel HER3 therapeutics are showing promising results while recently developed HER3 PET imaging modalities aid in predicting and assessing early treatment response. However, baseline HER3 expression, as well as changes in expression while on neoadjuvant therapy, have not been well-characterized. We conducted a prospective clinical study, pre- and post-neoadjuvant/systemic therapy, in patients with newly diagnosed breast cancer to determine HER3 expression, and to identify possible resistance mechanisms maintained through the HER3 receptor.Experimental designThe study was conducted between May 25, 2018 and October 12, 2019. Thirty-four patients with newly diagnosed breast cancer of any subtype (ER ± , PR ± , HER2 ±) were enrolled in the study. Two core biopsy specimens were obtained from each patient at the time of diagnosis. Four patients underwent a second research biopsy following initiation of neoadjuvant/systemic therapy or systemic therapy which we define as neoadjuvant therapy. Molecular characterization of HER3 and downstream signaling nodes of the PI3K/AKT and MAPK pathways pre- and post-initiation of therapy was performed. Transcriptional validation of finings was performed in an external dataset (GSE122630).ResultsVariable baseline HER3 expression was found in newly diagnosed breast cancer and correlated positively with pAKT across subtypes (r = 0.45). In patients receiving neoadjuvant/systemic therapy, changes in HER3 expression were variable. In a hormone receptor-positive (ER +/PR +/HER2-) patient, there was a statistically significant increase in HER3 expression post neoadjuvant therapy, while there was no significant change in HER3 expression in a ER +/PR +/HER2+ patient. However, both of these patients showed increased downstream signaling in the PI3K/AKT pathway. One subject with ER +/PR −/HER2− breast cancer and another subject with ER +/PR +/HER2 + breast cancer showed decreased HER3 expression. Transcriptomic findings, revealed an immune suppressive environment in patients with decreased HER3 expression post therapy.ConclusionThis study demonstrates variable HER3 expression across breast cancer subtypes. HER3 expression can be assessed early, post-neoadjuvant therapy, providing valuable insight into cancer biology and potentially serving as a prognostic biomarker. Clinical translation of neoadjuvant therapy assessment can be achieved using HER3 PET imaging, offering real-time information on tumor biology and guiding personalized treatment for breast cancer patients.

A Comprehensive Clinicopathologic and Molecular Reappraisal of GLI1 -altered Mesenchymal Tumors with Pooled Outcome Analysis Showing Poor Survival in GLI1 - amplified Versus GLI1- rearranged Tumors.

GLI1 -altered mesenchymal tumor is a recently described distinct pathologic entity with an established risk of malignancy, being defined molecularly by either GLI1 gene fusions or amplifications. The clinicopathologic overlap of tumors driven by the 2 seemingly distinct mechanisms of GLI1 activation is still emerging. Herein, we report the largest series of molecularly confirmed GLI1 -altered mesenchymal neoplasms to date, including 23 GLI1- amplified and 15 GLI1 -rearranged new cases, and perform a comparative clinicopathologic, genomic, and survival investigation. GLI1- rearranged tumors occurred in younger patients (42 vs. 52y) and were larger compared with GLI1 -amplified tumors (5.6cm vs. 1.5cm, respectively). Histologic features were overall similar between the 2 groups, showing a multinodular pattern and a nested architecture of epithelioid, and less commonly spindle cells, surrounded by a rich capillary network. A distinct whorling pattern was noted among 3 GLI1 -amplified tumors. Scattered pleomorphic giant cells were rarely seen in both groups. The immunoprofile showed consistent expression of CD56, with variable S100, CD10 and SMA expression. Genomically, both groups had overall low mutation burdens, with rare TP53 mutations seen only in GLI1- amplified tumors. GLI1 -amplified mesenchymal tumors exhibit mostly a single amplicon at the 12q13-15 locus, compared with dedifferentiated liposarcoma, which showed a 2-peak amplification centered around CDK4 (12q14.1) and MDM2 (12q15). GLI1 -amplified tumors had a significantly higher GLI1 mRNA expression compared with GLI1 -rearranged tumors. Survival pooled analysis of current and published cases (n=83) showed a worse overall survival in GLI1 -amplified patients, with 16% succumbing to disease compared with 1.7% in the GLI1- rearranged group. Despite comparable progression rates, GLI1 -amplified tumors had a shorter median progression-free survival compared with GLI1 -rearranged tumors (25mo vs. 77mo). Univariate analysis showed that traditional histologic predictors of malignancy (mitotic count ≥4/10 high-power fields, presence of necrosis, and tumor size ≥5cm) are associated with worse prognosis among GLI1 -altered mesenchymal tumors.

Prediction by genetic MATS of 4CMenB vaccine strain coverage of invasive meningococcal serogroup B isolates circulating in Taiwan between 2003 and 2020.

Neisseria meningitidis serogroup B (NmB) strains have diverse antigens, necessitating methods for predicting meningococcal serogroup B (MenB) vaccine strain coverage. The genetic Meningococcal Antigen Typing System (gMATS), a correlate of MATS estimates, predicts strain coverage by the 4-component MenB (4CMenB) vaccine in cultivable and non-cultivable NmB isolates. In Taiwan, 134 invasive, disease-causing NmB isolates were collected in 2003-2020 (23.1%, 4.5%, 5.2%, 29.8%, and 37.3% from individuals aged ≤11 months, 12-23 months, 2-4 years, 5-29 years, and ≥30 years, respectively). NmB isolates were characterized by whole-genome sequencing and vaccine antigen genotyping, and 4CMenB strain coverage was predicted using gMATS. Analysis of phylogenetic relationships with 502 global NmB genomes showed that most isolates belonged to three global hyperinvasive clonal complexes: ST-4821 (27.6%), ST-32 (23.9%), and ST-41/44 (14.9%). Predicted strain coverage by gMATS was 62.7%, with 27.6% isolates covered, 2.2% not covered, and 66.4% unpredictable by gMATS. Age group coverage point estimates ranged from 42.9% (2-4 years) to 66.1% (≤11 months). Antigen coverage estimates and percentages predicted as covered/not covered were highly variable, with higher estimates for isolates with one or more gMATS-positive antigens than for isolates positive for one 4CMenB antigen. In conclusion, this first study on NmB strain coverage by 4CMenB in Taiwan shows 62.7% coverage by gMATS, with predictable coverage for 29.8% of isolates. These could be underestimated since the gMATS calculation does not consider synergistic mechanisms associated with simultaneous antibody binding to multiple targets elicited by multicomponent vaccines or the contributions of minor outer membrane vesicle vaccine components.IMPORTANCEMeningococcal diseases, caused by the bacterium Neisseria meningitidis (meningococcus), include meningitis and septicemia. Although rare, invasive meningococcal disease is often severe and can be fatal. Nearly all cases are caused by six meningococcal serogroups (types), including meningococcal serogroup B. Vaccines are available against meningococcal serogroup B, but the antigens targeted by these vaccines have highly variable genetic features and expression levels, so the effectiveness of vaccination may vary depending on the strains circulating in particular countries. It is therefore important to test meningococcal serogroup B strains isolated from specific populations to estimate the percentage of bacterial strains that a vaccine can protect against (vaccine strain coverage). Meningococcal isolates were collected in Taiwan between 2003 and 2020, of which 134 were identified as serogroup B. We did further investigations on these isolates, including using a method (called gMATS) to predict vaccine strain coverage by the 4-component meningococcal serogroup B vaccine (4CMenB).

SEraster: a rasterization preprocessing framework for scalable spatial omics data analysis.

Spatial omics data demand computational analysis but many analysis tools have computational resource requirements that increase with the number of cells analyzed. This presents scalability challenges as researchers use spatial omics technologies to profile millions of cells. To enhance the scalability of spatial omics data analysis, we developed a rasterization preprocessing framework called SEraster that aggregates cellular information into spatial pixels. We apply SEraster to both real and simulated spatial omics data prior to spatial variable gene expression analysis to demonstrate that such preprocessing can reduce computational resource requirements while maintaining high performance, including as compared to other down-sampling approaches. We further integrate SEraster with existing analysis tools to characterize cell-type spatial co-enrichment across length scales. Finally, we apply SEraster to enable analysis of a mouse pup spatial omics dataset with over a million cells to identify tissue-level and cell-type-specific spatially variable genes as well as spatially co-enriched cell types that recapitulate expected organ structures. SEraster is implemented as an R package on GitHub (https://github.com/JEFworks-Lab/SEraster) with additional tutorials at https://JEF.works/SEraster.