Single Gene Research Articles

CGMega: explainable graph neural network framework with attention mechanisms for cancer gene module dissection

Cancer is rarely the straightforward consequence of an abnormality in a single gene, but rather reflects a complex interplay of many genes, represented as gene modules. Here, we leverage the recent advances of model-agnostic interpretation approach and develop CGMega, an explainable and graph attention-based deep learning framework to perform cancer gene module dissection. CGMega outperforms current approaches in cancer gene prediction, and it provides a promising approach to integrate multi-omics information. We apply CGMega to breast cancer cell line and acute myeloid leukemia (AML) patients, and we uncover the high-order gene module formed by ErbB family and tumor factors NRG1, PPM1A and DLG2. We identify 396 candidate AML genes, and observe the enrichment of either known AML genes or candidate AML genes in a single gene module. We also identify patient-specific AML genes and associated gene modules. Together, these results indicate that CGMega can be used to dissect cancer gene modules, and provide high-order mechanistic insights into cancer development and heterogeneity.

Localized in vivo gene editing of murine cancer-associated fibroblasts.

Discovering the role of fibroblasts residing in the tumor microenvironment (TME) requires controlled, localized perturbations because fibroblasts play critical roles in regulating immunity and tumor biology at multiple sites. Systemic perturbations can lead to unintended, confounding secondary effects, and methods to locally genetically engineer fibroblasts are lacking. To specifically investigate murine stromal cell perturbations restricted to the TME, we developed an adeno-associated virus (AAV)-based method to target any gene-of-interest in fibroblasts at high efficiency (>80%). As proof of concept, we generated single (sKO) and double gene KOs (dKO) of Osmr , Tgfbr2 , and Il1r1 in cancer-associated fibroblasts (CAFs) and investigated how their cell states and those of other cells of the TME subsequently change in mouse models of melanoma and pancreatic ductal adenocarcinoma (PDAC). Furthermore, we developed an in vivo knockin-knockout (KIKO) strategy to achieve long-term tracking of CAFs with target gene KO via knocked-in reporter gene expression. This validated in vivo gene editing toolbox is fast, affordable, and modular, and thus holds great potential for further exploration of gene function in stromal cells residing in tumors and beyond.

Confined placental mosaicism is a diagnostic pitfall in dystrophinopathies: a clinical report.

Single-gene copy number variants (CNVs) limited to placenta although rarely identified may have clinical implications. We describe a pregnant woman referred for chorionic villus sampling due to increased fetal nuchal translucency. Incident intragenic deletion of Duchenne muscular dystrophy (DMD) gene, affecting exons 56 and 57, was identified in a male fetus in ~23-30% of placental cells by chromosomal microarray and confirmed using multiplex ligation-dependent probe amplification (MLPA). Rapid aneuploidy testing showed normal results and the deletion was not detected in the mother. Subsequent analyses on amniotic cells yielded a normal DMD gene result, corroborating the confined placental nature of the mosaicism. Hence, this report emphasizes the importance of conducting amniocentesis following detection of mosaicism for single gene CNVs on chorionic villi, in order to preclude confined placental mosaicism (CPM). As far as we know, this report marks only the second documented situation of CPM involving an intragenic DMD deletion.

Mapping of Aegilops speltoides derived leaf rust and stripe rust resistance genes using 35K SNP array

Wheat is an essential food commodity cultivated throughout the world. However, this crop faces continuous threats from fungal pathogens, leaf rust (LR) and stripe rust (YR). To continue feeding the growing population, these major destructors of wheat must be effectively countered by enhancing the genetic diversity of cultivated germplasm. In this study, an introgression line with hexaploid background (ILsp3603) carrying resistance against Pt pathotypes 77−5 (121R63-1), 77−9 (121R60-1) and Pst pathotypes 46S119 (46E159), 110S119 (110E159), 238S119 (238E159) was developed from donor wheat wild progenitor, Aegilops speltoides acc pau 3603. To understand the genetic basis of resistance and map these genes (named Lrsp3603 and Yrsp3603), inheritance studies were carried out in F6 and F7 mapping population, developed by crossing ILsp3603 with LR and YR susceptible cultivar WL711, which revealed a monogenic (single gene) inheritance pattern for each of these traits. Bulk segregant analysis combined with 35 K Axiom SNP array genotyping mapped both genes as separate entities on the short arm of chromosome 6B. A genetic linkage map, comprising five markers, 1 SNP, 1 PLUG and three gene based SSRs, covered a genetic distance of 12.65 cM. Lrsp3603 was flanked by markers Tag-SSR14 (located proximally at 2.42 cM) and SNP AX-94542331 (at 3.28 cM) while Yrsp3603 was mapped at one end closest to AX-94542331 at 6.62 cM distance. Functional annotation of Lrsp3603 target region (∼ 1 Mbp) revealed 10 gene IDs associated with disease resistance mechanisms including three encoding typical R gene domains.

Revealing the novel metabolism-related genes in the ossification of the ligamentum flavum based on whole transcriptomic data.

The ossification of the ligamentum flavum (OLF) is one of the major causes of thoracic myelopathy. Previous studies indicated there might be a potential link between metabolic disorder and pathogenesis of OLF. The aim of this study was to determine the potential role of metabolic disorder in the pathogenesis of OLF using the strict bioinformatic workflow for metabolism-related genes and experimental validation. A series of bioinformatic approaches based on metabolism-related genes were conducted to compare the metabolism score between OLF tissues and normal ligamentum flavum (LF) tissues using the single sample gene set enrichment analysis. The OLF-related and metabolism-related differentially expressed genes (OMDEGs) were screened out, and the biological functions of OMDEGs were explored, including the Gene Ontology enrichment analysis, Kyoto Encyclopedia of Genes and Genomes enrichment analysis, and protein-protein interaction. The competing endogenous RNA (ceRNA) network based on pairs of miRNA-hub OMDEGs was constructed. The correlation analysis was conducted to explore the potential relationship between metabolic disorder and immunity abnormality in OLF. In the end, the cell experiments were performed to validate the roles of GBE1 and TNF-α in the osteogenic differentiation of LF cells. There was a significant difference of metabolism score between OLF tissues and normal LF tissues. Forty-nine OMDEGs were screened out and their biological functions were determined. The ceRNA network containing three hub OMDEGs and five differentially expressed miRNAs (DEmiRNAs) was built. The correlation analysis between hub OMDEGs and OLF-related infiltrating immune cells indicated that metabolic disorder might contribute to the OLF via altering the local immune status of LF tissues. The cell experiments determined the important roles of GBE1 expression and TNF-α in the osteogenic differentiation of LF cells. This research, for the first time, preliminarily illustrated the vital role of metabolic disorder in the pathogenesis of OLF using strict bioinformatic algorithms and experimental validation for metabolism-related genes, which could provide new insights for investigating disease mechanism and screening effective therapeutic targets of OLF in the future.

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment—social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Genetic parameters, genome-wide association study, and selection perspective on gestation length in 16 French cattle breeds

Lifetime productivity is a trait of great importance to dairy cattle populations as it combines information from production and longevity variables. Therefore, we investigated the genetic background of lifetime productivity in high-producing dairy cattle by integrating genomics and transcriptomics data sets. A total of 3,365,612 test-day milk yield records from 134,029 Chinese Holstein cows were used to define 6 lifetime productivity traits, including lifetime milk yield covering full lifespan and 5 cumulative milk yield traits covering partial lifespan. Genetic parameters were estimated based on univariate and bivariate linear animal models and the Restricted Maximum Likelihood (REML) method. Genome-wide association studies (GWAS) and weighted gene co-expression network analyses (WGCNA) were performed to identify candidate genes associated with lifetime productivity based on genomic data from 3,424 cows and peripheral blood RNA-seq data from 23 cows, respectively. Lifetime milk yield averaged 24,800.8 ± 14,396.6 kg (mean ± SD) across an average of 2.4 parities in Chinese Holstein population. The heritability estimates for lifetime productivity traits ranged from 0.05 (±0.01 for SE) to 0.10 (±0.02 for SE). The estimate of genetic correlation between lifetime milk yield and productive life is 0.88 (±0.3 for SE) while the genetic correlation with 305d milk yield in the first lactation was 0.49 (±0.08 for SE). Absolute values for most genetic correlation estimates between lifetime productivity and type traits were lower than 0.30. Moderate genetic correlations were found between udder related traits and lifetime productivity, such as with udder depth (0.33), rear udder attachment height (0.33), and udder system (0.34). Some single nucleotide polymorphisms and gene co-expression modules significantly associated with lifetime milk yield were identified based on GWAS and WGCNA analyses, respectively. Functional enrichment analyses of the candidate genes identified revealed important pathways related to immune system, longevity, energy utilization and metabolism, and FoxO signaling. The genes NTMT1, FNBP1, and S1PR1 were considered to be the most important candidate genes influencing lifetime productivity in Holstein cows. Overall, our findings indicate that lifetime productivity is heritable in Chinese Holstein cattle and important candidate genes were identified by integrating genomic and transcriptomic data sets.

Early-onset lupus nephritis.

Early-onset systemic lupus erythematous (SLE) is a distinct clinical entity characterized by the onset of disease manifestations during childhood. Despite some similarities to patients who are diagnosed during adulthood, early-onset SLE typically displays a greater disease severity, with aggressive multiorgan involvement, lower responsiveness to classical therapies, and more frequent flares. Lupus nephritis is one of the most severe complications of SLE and represents a major risk factor for long-term morbidity and mortality, especially in children. This review focuses on the clinical and histological aspects of early-onset lupus nephritis, aiming at highlighting relevant differences with adult patients, emphasizing long-term outcomes and discussing the management of long-term complications. We also discuss monogenic lupus, a spectrum of conditions caused by single gene variants affecting the complement cascade, extracellular and intracellular nucleic acid sensing and processing, and occasionally other metabolic pathways. These monogenic forms typically develop early in life and often have clinical manifestations that resemble sporadic SLE, whereas their response to standard treatments is poor.

Two C-terminal isoforms of Aplysia tachykinin–related peptide receptors exhibit phosphorylation-dependent and phosphorylation-independent desensitization mechanisms

Diversity, a hallmark of G protein-coupled receptor (GPCR) signaling, partly stems from alternative splicing of a single gene generating more than one isoform for a receptor. Additionally, receptor responses to ligands can be attenuated by desensitization upon prolonged or repeated ligand exposure. Both phenomena have been demonstrated and exemplified by the deuterostome tachykinin (TK) signaling system, although the role of phosphorylation in desensitization remains a subject of debate. Here, we describe the signaling system for tachykinin-related peptides (TKRPs) in a protostome, mollusk Aplysia. We cloned the Aplysia TKRP precursor, which encodes three TKRPs (apTKRP-1, apTKRP-2a, and apTKRP-2b) containing the FXGXR-amide motif. In situ hybridization and immunohistochemistry showed predominant expression of TKRP mRNA and peptide in the cerebral ganglia. TKRPs and their post-translational modifications were observed in extracts of CNS ganglia using mass spectrometry. We identified two Aplysia TKRP receptors (TKRPRs), named apTKRPR-A and apTKRPR-B. These receptors are two isoforms generated through alternative splicing of the same gene and differ only in their intracellular C-termini. Structure-activity relationship analysis of apTKRP-2b revealed that both C-terminal amidation and conserved residues of the ligand are critical for receptor activation. C-terminal truncates and mutants of apTKRPRs suggested that there is a C-terminal phosphorylation-independent desensitization for both receptors. Moreover, apTKRPR-B also exhibits phosphorylation-dependent desensitization through the phosphorylation of C-terminal Ser/Thr residues. This comprehensive characterization of the Aplysia TKRP signaling system underscores the evolutionary conservation of the TKRP and TK signaling systems, while highlighting the intricacies of receptor regulation through alternative splicing and differential desensitization mechanisms.

Evidence of gene-nutrient interaction association with waist circumference, cross-sectional analysis

Background and aimsWaist circumference (WC) is a significant indicator of body adiposity and is associated with increased mortality and morbidity of cardiovascular diseases. Although, single nutrient intake and candidate genes were previously associated with WC. Little is known about WC association with overall diet quality, genetic risk score and gene-nutrient interaction. This study aims to investigate the influence of overall diet quality and multiple WC-associated single nucleotide polymorphisms on WC. In addition to investigating gene-nutrient interaction association with WC.MethodsThis study explored cross-sectional data from two large sample-size studies, to provide reproducible results. As a representation of the UK population, the Airwave Health Monitoring Study (n = 6,502) and the UK-Biobank Cohort Study (n = 171,129) were explored for factors associated with WC. Diet quality was evaluated based on the Mellen Index for Dietary Approaches to Stop Hypertension (Mellen-DASH). The genetic risk score for WC (GRS-Waist) was calculated by screening the population genotype for WC-associated single nucleotide polymorphisms. Multivariate linear regression models were built to explore WC association with diet quality and genetic risk score. Gene-nutrient interaction was explored by introducing the interaction term (GRS-Waist X Mellen-DASH score) to multivariate linear regression analysis.ResultsThe prevalence of high WC (Female > 80 cm, Male > 94 cm) was 46.5% and 51.7% in both populations. Diet quality and genetic risk score of WC were significantly associated with WC. There was no evidence of interaction between GRS-Waist, DASH diet scores and nutrient intake on WC.ConclusionThis study’s findings provided reproducible results on waist circumference association with diet and genetics and tested the possibility of gene-nutrient interaction. These reproducible results are successful in building the foundation for using diet and genetics for early identification of those at risk of having high WC and WC-associated diseases. In addition, evidence on gene-diet interactions on WC is limited and lacks replication, therefore our findings may guide future research in investigating this interaction and investigating its application in precision nutrition.

Pathogenomics analysis of high-risk clone ST147 multidrug-resistant Klebsiella pneumoniae isolated from a patient in Egypt

BackgroundThe emergence of multi-drug-resistant Klebsiella pneumoniae (MDR-KP) represents a serious clinical health concern. Antibiotic resistance and virulence interactions play a significant role in the pathogenesis of K. pneumoniae infections. Therefore, tracking the clinical resistome and virulome through monitoring antibiotic resistance genes (ARG) and virulence factors in the bacterial genome using computational analysis tools is critical for predicting the next epidemic.MethodsIn the current study, one hundred extended spectrum β-lactamase (ESBL)-producing clinical isolates were collected from Mansoura University Hospital, Egypt, in a six-month period from January to June 2022. One isolate was selected due to the high resistance phenotype, and the genetic features of MDR-KP recovered from hospitalized patient were investigated. Otherwise, the susceptibility to 25 antimicrobials was determined using the DL Antimicrobial Susceptibility Testing (AST) system. Whole genome sequencing (WGS) using Illumina NovaSeq 6000 was employed to provide genomic insights into K. pneumoniae WSF99 clinical isolate.ResultsThe isolate K. pneumoniae WSF99 was phenotypically resistant to the antibiotics under investigation via antibiotic susceptibility testing. WGS analysis revealed that WSF99 total genome length was 5.7 Mb with an estimated 5,718 protein-coding genes and a G + C content of 56.98 mol%. Additionally, the allelic profile of the WSF99 isolate was allocated to the high-risk clone ST147. Furthermore, diverse antibiotic resistance genes were determined in the genome that explain the high-level resistance phenotypes. Several β-lactamase genes, including blaCTX−M−15, blaTEM−1, blaTEM−12, blaSHV−11, blaSHV−67, and blaOXA−9, were detected in the WSF99 isolate. Moreover, a single carbapenemase gene, blaNDM−5, was predicted in the genome, positioned within a mobile cassette. In addition, other resistance genes were predicted in the genome including, aac(6’)-Ib, aph(3’)-VI, sul1, sul2, fosA, aadA, arr-2, qnrS1, tetA and tetC. Four plasmid replicons CoIRNAI, IncFIB(K), IncFIB(pQil), and IncR were predicted in the genome. The draft genome analysis revealed the occurrence of genetic mobile elements positioned around the ARGs, suggesting the ease of dissemination via horizontal gene transfer.ConclusionsThis study reports a comprehensive pathogenomic analysis of MDR-KP isolated from a hospitalized patient. These findings could be relevant for future studies investigating the diversity of antimicrobial resistance and virulence in Egypt.

An all-out assault on a dominant resistance gene: Local emergence, establishment, and spread of strains of tomato spotted wilt orthotospovirus (TSWV) that overcome Sw-5b-mediated resistance in fresh market and processing tomatoes in California.

Tomato spotted wilt orthotospovirus (TSWV) causes substantial economic loss to tomato production, and the Sw-5b resistance gene is widely deployed for management. Here, we show (i) the emergence of resistance-breaking (RB) TSWV strains in processing and fresh market tomato production in California over the past ten years, and (ii) evolutionary relationships with RB strains from other areas. A specific RT-PCR test was used to show the C118Y RB strain that emerged in Fresno County in 2016 quickly became predominant in the central production area and remained so through this study. In 2021, the C118Y strain was detected in the Northern production area, and was predominant in 2022. However, in 2023, the C118Y strain was unexpectedly detected in fewer spotted wilt samples from resistant varieties. This was due to emergence of the T120N RB strain, previously known to occur in Spain. A specific RT-PCR test was developed and used to show that the T120N RB strain was predominant in Colusa and Sutter counties (detected in 75-80% of samples), and detected in ~50% of samples from Yolo County. Pathogenicity tests confirmed California isolates of the T120N strain infected Sw-5b tomato varieties and induced severe symptoms. Another RB strain, C118F, was associated with spotted wilt samples of Sw-5 varieties from fresh market tomato production in southern California. Phylogenetic analyses with complete NSm sequences revealed that the C118Y and T120N RB strains infecting resistant processing tomato in California emerged locally, whereas those from fresh market production were more closely related to isolates from Mexico. Thus, widespread deployment of this single dominant resistance gene in California has driven the local emergence of multiple RB strains in different tomato production areas and types. These results further emphasize the need for ongoing monitoring for RB strains, and identification of sources of resistance to these strains.

Dominant suppressor genes of p53-induced apoptosis in Drosophila melanogaster.

One of the major functions of programmed cell death (apoptosis) is the removal of cells that suffered oncogenic mutations, thereby preventing cancerous transformation. By making use of a Double-Headed-EP (DEP) transposon, a P element derivative made in our laboratory, we made an insertional mutagenesis screen in Drosophila melanogaster to identify genes that, when overexpressed, suppress the p53-activated apoptosis. The DEP element has Gal4-activatable, outward-directed UAS promoters at both ends, which can be deleted separately in vivo. In the DEP insertion mutants, we used the GMR-Gal4 driver to induce transcription from both UAS promoters and tested the suppression effect on the apoptotic rough eye phenotype generated by an activated UAS-p53 transgene. By DEP insertions, 7 genes were identified, which suppressed the p53-induced apoptosis. In 4 mutants, the suppression effect resulted from single genes activated by 1 UAS promoter (Pka-R2, Rga, crol, and Spt5). In the other 3 (Orct2, Polr2M, and stg), deleting either UAS promoter eliminated the suppression effect. In qPCR experiments, we found that the genes in the vicinity of the DEP insertion also showed an elevated expression level. This suggested an additive effect of the nearby genes on suppressing apoptosis. In the eukaryotic genomes, there are coexpressed gene clusters. Three of the DEP insertion mutants are included, and 2 are in close vicinity of separate coexpressed gene clusters. This raises the possibility that the activity of some of the genes in these clusters may help the suppression of the apoptotic cell death.

Combination of error-prone PCR (epPCR) and Circular Polymerase Extension Cloning (CPEC) for improving the coverage of random mutagenesis libraries

Random mutagenesis, such as error-prone PCR (epPCR), is a technique capable of generating a wide variety of a single gene. However, epPCR can produce a large number of mutated gene variants, posing a challenge in ligating these mutated PCR products into plasmid vectors. Typically, the primers for mutagenic PCRs incorporate artificial restriction enzyme sites compatible with chosen plasmids. Products are cleaved and ligated to linearized plasmids, then recircularized by DNA ligase. However, this cut-and-paste method known as ligation-dependent process cloning (LDCP), has limited efficiency, as the loss of potential mutants is inevitable leading to a significant reduction in the library’s breadth. An alternative to LDCP is the circular polymerase extension cloning (CPEC) method. This technique involves a reaction where a high-fidelity DNA polymerase extends the overlapping regions between the insert and vector, forming a circular molecule. In this study, our objective was to compare the traditional cut-and-paste enzymatic method with CPEC in producing a variant library from the gene encoding the red fluorescent protein (DsRed2) obtained by epPCR. Our findings suggest that CPEC can accelerate the cloning process in gene library generation, enabling the acquisition of a greater number of gene variants compared to methods reliant on restriction enzymes.

Unmutated RRAS2 emerges as a key oncogene in post-partum-associated triple negative breast cancer

BackgroundBreast cancer (BC) is the most common cancer in women, with triple negative BC (TNBC) accounting for 20% of cases. While early detection and targeted therapies have improved overall life expectancy, TNBC remains resistant to current treatments. Although parity reduces the lifetime risk of developing BC, pregnancy increases the risk of developing TNBC for years after childbirth. Although numerous gene mutations have been associated with BC, no single gene alteration has been identified as a universal driver. RRAS2 is a RAS-related GTPase rarely found mutated in cancer.MethodsConditional knock-in mice were generated to overexpress wild type human RRAS2 in mammary epithelial cells. A human sample cohort was analyzed by RT-qPCR to measure RRAS2 transcriptional expression and to determine the frequency of both a single-nucleotide polymorphism (SNP rs8570) in the 3’UTR region of RRAS2 and of genomic DNA amplification in tumoral and non-tumoral human BC samples.ResultsHere we show that overexpression of wild-type RRAS2 in mice is sufficient to develop TNBC in 100% of females in a pregnancy-dependent manner. In human BC, wild-type RRAS2 is overexpressed in 68% of tumors across grade, location, and molecular type, surpassing the prevalence of any previously implicated alteration. Still, RRAS2 overexpression is notably higher and more frequent in TNBC and young parous patients. The increased prevalence of the alternate C allele at the SNP position in tumor samples, along with frequent RRAS2 gene amplification in both tumors and blood of BC patients, suggests a cause-and-effect relationship between RRAS2 overexpression and breast cancer.ConclusionsHigher than normal expression of RRAS2 not bearing activating mutations is a key driver in the majority of breast cancers, especially those of the triple-negative type and those linked to pregnancy.

Water Versus Electrolyte Rehydration on Vocal Fold Osmotic and Oxidative Stress Gene Expression.

Systemic dehydration may induce osmotic and oxidative stress in the vocal folds, but our knowledge of the biology and mitigation with rehydration is limited. The purpose of this experiment was to evaluate whether systemic dehydration induces vocal fold oxidative and osmotic stress and to compare the impact of rehydration by water intake versus electrolyte intake on osmotic and oxidative stress-related gene expression. Four-month-old male Sprague-Dawley rats (N = 32) underwent water restriction. Rehydration was achieved with ad libitum access to water or electrolytes for 24 hours. Rats were divided into four groups: euhydration control, dehydration-only, dehydration followed by either water or electrolyte rehydration (n = 8/group). Gene expression was assessed via RT2 Gene Expression Profiler arrays. With respect to oxidative stress, 10 genes were upregulated and 2 were downregulated after vocal fold dehydration compared with the euhydrated control. Concerning osmotic stress, six genes were upregulated with dehydration only, six genes were upregulated following rehydration with water, whereas a single gene was upregulated with electrolyte rehydration. All genes with significantly different expression between the rehydration groups showed lower expression with electrolytes compared with water. The results support a potential role of oxidative and osmotic stresses in vocal folds related to systemic dehydration. The differences in stress-related gene expression in vocal fold tissue between rehydration with electrolytes or water, albeit modest, suggest that both rehydration options offer clinical utility to subjects experiencing vocal fold dehydration with preliminary evidence that electrolytes may be more effective than water in resolving osmotic stress. NA (prospective animal study) Laryngoscope, 134:4636-4641, 2024.

Molecular characterization of the salivary adenoid cystic carcinoma immune landscape by anatomic subsites

Adenoid cystic carcinoma (AdCC) is a slow-growing salivary gland malignancy that relapses frequently. AdCCs of the submandibular gland exhibit unique differences in prognosis and treatment response to adjuvant radiotherapy compared to other sites, yet the role of tumor anatomic subsite on gene expression and tumor immune microenvironment (TIME) composition remains unclear. We used 87 samples, including 48 samples (27 AdCC and 21 normal salivary gland tissue samples) from 4 publicly available AdCC RNA sequencing datasets, a validation set of 33 minor gland AdCCs, and 39 samples from an in-house cohort (30 AdCC and 9 normal salivary gland samples). RNA sequencing data were used for single sample gene set enrichment analysis and TIME deconvolution. Quantitative PCR and multiplex immunofluorescence were performed on the in-house cohort. Wilcoxon rank-sum, nonparametric equality-of-medians tests and linear regression models were used to evaluate tumor subsite differences. AdCCs of different anatomic subsites including parotid, submandibular, sublingual, and minor salivary glands differed with respect to expression of several key tumorigenic pathways. Among the three major salivary glands, the reactive oxygen species (ROS)/nuclear factor erythroid 2-related factor 2 (NRF2) pathway signature was significantly underexpressed in AdCC of submandibular compared to parotid and sublingual glands while this association was not observed among normal glands. Additionally, the NRF2 pathway, whose expression was associated with favorable overall survival, was overexpressed in AdCCs of parotid gland compared to minor and submandibular glands. The TIME deconvolution identified differences in CD4+ T cell populations between AdCC of major and minor glands and natural killer (NK) cells among AdCC of minor, submandibular, and parotid glands while plasma cells were enriched in normal submandibular glands compared to other normal gland controls. Our data reveal key molecular differences in AdCC of different anatomic subsites. The ROS and NRF2 pathways are underexpressed in submandibular and minor AdCCs compared to parotid gland AdCCs, and NRF2 pathway expression is associated with favorable overall survival. The CD4+ T, NK, and plasma cell populations also vary by tumor subsites, suggesting that the observed submandibular AdCC tumor-intrinsic pathway differences may be responsible for influencing the TIME composition and survival differences.

Accurate estimation of pathway activity in single cells for clustering and differential analysis.

Inferring which and how biological pathways and gene sets change is a key question in many studies that utilize single-cell RNA sequencing. Typically, these questions are addressed by quantifying the enrichment of known gene sets in lists of genes derived from global analysis. Here we offer SiPSiC, a new method to infer pathway activity in every single cell. This allows more sensitive differential analysis and utilization of pathway scores to cluster cells and compute UMAP or other similar projections. We apply our method to COVID-19, lung adenocarcinoma and glioma data sets, and demonstrate its utility. SiPSiC analysis results are consistent with findings reported in previous studies in many cases, but SiPSiC also reveals the differential activity of novel pathways, enabling us to suggest new mechanisms underlying the pathophysiology of these diseases and demonstrating SiPSiC's high accuracy and sensitivity in detecting biological function and traits. In addition, we demonstrate how it can be used to better classify cells based on activity of biological pathways instead of single genes and its ability to overcome patient-specific artifacts.

A single gene mutation underpins metabolic adaptation and acquisition of filamentous competence in the emerging fungal pathogen Candida auris.

Filamentous cell growth is a vital property of fungal pathogens. The mechanisms of filamentation in the emerging multidrug-resistant fungal pathogen Candida auris are poorly understood. Here, we show that exposure of C. auris to glycerol triggers a rod-like filamentation-competent (RL-FC) phenotype, which forms elongated filamentous cells after a prolonged culture period. Whole-genome sequencing analysis reveals that all RL-FC isolates harbor a mutation in the C2H2 zinc finger transcription factor-encoding gene GFC1 (Gfc1 variants). Deletion of GFC1 leads to an RL-FC phenotype similar to that observed in Gfc1 variants. We further demonstrate that GFC1 mutation causes enhanced fatty acid β-oxidation metabolism and thereby promotes RL-FC/filamentous growth. This regulation is achieved through a Multiple Carbon source Utilizer (Mcu1)-dependent mechanism. Interestingly, both the evolved RL-FC isolates and the gfc1Δ mutant exhibit an enhanced ability to colonize the skin. Our results reveal that glycerol-mediated GFC1 mutations are beneficial during C. auris skin colonization and infection.

LASSO-derived nomogram for early identification of pediatric monogenic lupus.

Monogenic lupus is defined as systemic lupus erythematosus (SLE)/SLE-like patients with either dominantly or recessively inherited pathogenic variants in a single gene with high penetrance. However, because the clinical phenotype of monogenic SLE is extensive and overlaps with that of classical SLE, it causes a delay in diagnosis and treatment. Currently, there is a lack of early identification models for clinical practitioners to provide early clues for recognition. Our goal was to create a clinical model for the early identification of pediatric monogenic lupus, thereby facilitating early and precise diagnosis and treatment for patients. This retrospective cohort study consisted of 41 cases of monogenic lupus treated at the Department of Pediatrics at Peking Union Medical College Hospital from June 2012 to December 2022. The control group consisted of classical SLE patients recruited at a 1:2 ratio. Patients were randomly divided into a training group and a validation group at a 7:3 ratio. A logistic regression model was established based on the least absolute shrinkage and selection operator to generate the coefficient plot. The predictive ability of the model was evaluated using receiver operator characteristic curves and the area under the curve (AUC) index. A total of 41 cases of monogenic lupus patients and 82 cases of classical SLE patients were included. Among the monogenic lupus cases (with a male-to-female ratio of 1:1.05 and ages of onset ranging from birth to 15years), a total of 18 gene mutations were identified. The variables included in the coefficient plot were age of onset, recurrent infections, intracranial calcifications, growth and developmental delay, abnormal muscle tone, lymphadenopathy/hepatosplenomegaly, and chilblain-like skin rash. Our model demonstrated satisfactory diagnostic performance through internal validation, with an AUC value of 0.97 (95% confidence interval = 0.92-0.97). We summarized and analyzed the clinical characteristics of pediatric monogenic lupus and developed a predictive model for early identification by clinicians. Clinicians should exercise high vigilance for monogenic lupus when the score exceeds - 9.032299.