Miracle-dz/UniFAR-FTUs huggingface dataset:
Global gene-expression analysis reveals the molecular processes underlying ClC-5 loss-of-function in novel Dent Disease 1 cellular models

Dent disease 1 (DD1) is a rare X-linked renal proximal tubulopathy characterized by low molecular weight proteinuria (LMWP) and variable degree of hypercalciuria, nephrocalcinosis and/or nephrolithiasis with progression to chronic kidney disease (CKD). Although loss-of-function mutations in the gene CLCN5 encoding the electrogenic Cl−/H+ antiporter ClC-5, which impair endocytic uptake in proximal tubule cells, cause the disease, there is poor genotype-phenotype correlation and their contribution to proximal tubule dysfunction remains unclear. Here, in order to discover the mechanisms leading to proximal tubule dysfunction due to ClC-5 loss-of-function, we have generated and characterized new human cellular models of DD1 by silencing CLCN5 and introducing the ClC-5 pathogenic mutants V523del, E527D and I524K into the human proximal tubule-derived cell line RPTEC/TERT1. Depletion of CLCN5 or expression of mutant ClC-5 impairs albumin endocytosis, increases substrate adhesion and decreases collective migration, which correlates with a less differentiated epithelial phenotype. Interestingly, although all conditions compromised the endocytic capacity in a similar way, their impact on gene expression profiles was different. Our DNA microarray studies show that ClC-5 silencing or mutant re-introduction alter pathways related to nephron development, anion homeostasis, organic acid transport, extracellular matrix organization and cell migration, compared to control cells. Cells carrying the V523del ClC-5 mutation show the largest differences in gene expression vs WT cells, which is in agreement with the more aggressive clinical phenotype observed in some DD1 patients. Overall, this work emphasizes the use of human proximal tubule derived cell models to identify the molecular processes underlying ClC-5 deficiency.
Methylation of Dual Specificity Phosphatase 4 Controls Cell Differentiation

A collection of signaling and epigenetic events needs to be orchestrated for normal development of hematopoietic lineages. While mitogen-activated protein (MAP) kinases (MAPKs) and multiple epigenetic modulators have been implicated in the megakaryocytic (Mk) cell differentiation, the underlying molecular mechanisms of signaling-epigenetic crosstalk remain unclear. MAPKs are in general inactivated by dual specificity phosphatases (DUSPs), whose activities are tightly regulated by various posttranslational modifications. Using knockdown screening and single-cell transcriptional analysis, we determined that DUSP4 is the phosphatase that inactivates p38 MAPK in hematopoietic cells and serves as a key regulator to promote Mk differentiation. With the nextgeneration Bioorthogonal Profiling of Protein Methylation technology for live cells, we identified DUSP4 as a PRMT1 substrate. Mechanistically, PRMT1-mediated Arg351 methylation of DUSP4 triggers its ubiquitinylation by HUWE1 (an E3 ligase) and then degradation, which results in p38 MAPK activation and inhibition of Mk differentiation in vitro and in vivo. Interestingly, the mechanistic axis of the DUSP4 degradation and p38 activation is also associated with a transcriptional signature of immune activation and thus argues immunological roles of Mk cells. Collectively, these results demonstrate a critical role of PRMT1-mediated posttranslational modification of DUSP4 in regulation of Mk differentiation and maturation. In the context of thrombocytopenia observed in myelodysplastic syndromes (MDS), we demonstrated that high levels of p38 MAPK and PRMT1 are associated with low platelet counts and adverse prognosis, while pharmacological inhibition of p38 MAPK or PRMT1 stimulates megakaryopoiesis in MDS samples. These findings provide novel mechanistic insights into the role of the PRMT1-DUSP4-p38 axis on Mk differentiation and present a targeting strategy for treatment of thrombocytopenia associated with myeloid malignancies such as MDS.
Oligomerization of the Human Adenosine A2A Receptor is Driven by the Intrinsically Disordered C-terminus

G protein-coupled receptors (GPCRs) have long been shown to exist as oligomers with functional properties distinct from those of the monomeric counterparts, but the driving factors of GPCR oligomerization remain relatively unexplored. In this study, we focus on the human adenosine A2A receptor (A2AR), a model GPCR that forms oligomers both in vitro and in vivo. Combining experimental and computational approaches, we discover that the intrinsically disordered C-terminus of A2AR drives the homo-oligomerization of the receptor. The formation of A2AR oligomers declines progressively and systematically with the shortening of the C-terminus. Multiple interaction sites and types are responsible for A2AR oligomerization, including disulfide linkages, hydrogen bonds, electrostatic interactions, and hydrophobic interactions. These interactions are enhanced by depletion interactions along the C-terminus, forming a tunable network of bonds that allow A2AR oligomers to adopt multiple interfaces. This study uncovers the disordered C-terminus as a prominent driving factor for the oligomerization of a GPCR, offering important guidance for structure-function studies of A2AR and other GPCRs.
Grammatical category and the neural processing of phrases

The interlocking roles of lexical, syntactic and semantic processing in language comprehension has been the subject of longstanding debate. Recently, the cortical response to a frequency-tagged linguistic stimulus has been shown to track the rate of phrase and sentence, as well as syllable, presentation. This could be interpreted as evidence for the hierarchical processing of speech, or as a response to the repetition of grammatical category. To examine the extent to which hierarchical structure plays a role in language processing we recorded EEG from human participants as they listen to isochronous streams of monosyllabic words. Comparing responses to sequences in which grammatical category is strictly alternating and chosen such that two-word phrases can be grammatically constructed—cold food loud room—or is absent—rough give ill tell—showed cortical entrainment at the two-word phrase rate was only present in the grammatical condition. Thus, grammatical category repetition alone does not yield entertainment at higher level than a word. On the other hand, cortical entrainment was reduced for the mixed-phrase condition that contained two-word phrases but no grammatical category repetition—that word send less—which is not what would be expected if the measured entrainment reflected purely abstract hierarchical syntactic units. Our results support a model in which word-level grammatical category information is required to build larger units.
The jasmine (Jasminum sambac) genome and flower fragrances

Jasminum sambac, a world-renowned plant appreciated for its exceptional flower fragrance, is of cultural and economic importance. However, the genetic basis of its fragrance is largely unknown. Here, we present the first de novo genome of J. sambac with 550.12 Mb (scaffold N50 = 40.1 Mb) assembled into 13 pseudochromosomes. Terpene synthase genes associated with flower fragrance are significantly amplified in the form of gene clusters through tandem duplications in the genome. Eleven homolog genes within the SABATH super-family were identified as related to phenylpropanoid/benzenoid compounds. Several key genes regulating jasmonate biosynthesis were duplicated causing increased copy numbers. Furthermore, multi-omics analyses identified various aromatic compounds and the key genes involved in fragrance biosynthesis pathways. Our genome of J. sambac offers a basic genetic resource for studying floral scent biosynthesis and provides an essential foundation for functional genomic research and variety improvements in Jasminum.
Higher-order structures of the foot-and-mouth disease virus RNA-dependent RNA polymerase required for dynamic inter-molecular interactions involved in viral genome replication

Replication of many positive-sense RNA viruses occurs within intracellular membrane-associated compartments. These are believed to provide a favourable environment for replication to occur, concentrating essential viral structural and non-structural components, as well as protecting these components from host-cell pathogen recognition and innate immune responses. However, the details of the molecular interactions and dynamics within these structures is very limited. One of the key components of the replication machinery is the RNA-dependent RNA polymerase, RdRp. This enzyme has been shown to form higher-order fibrils in vitro. Here, using the RdRp from foot-and-mouth disease virus (termed 3Dpol), we report fibril structures, solved at ~7-9 Å resolution by cryo-EM, revealing multiple conformations of a flexible assembly. Fitting high-resolution coordinates led to the definition of potential intermolecular interactions. We employed mutagenesis using a sub-genomic replicon system to probe the importance of these interactions for replication. We use these data to propose models for the role of higher order 3Dpol complexes as a dynamic scaffold within which RNA replication can occur.
Genotype-phenotype map of an RNA-ligand complex

RNA-ligand interactions play important roles in biology and biotechnology, but they often involve complex three-dimensional folding of RNA and are difficult to predict. To systematically explore the phenotypic landscape of an RNA-ligand complex, we used microarrays to investigate all possible single and double mutants of the 49-nt RNA aptamer Broccoli bound to the fluorophore DFHBI-1T. We collected more than seven million fluorescence measurements in varying conditions, and inferred dissociation rate constants, spectral shifts, and intragenic epistasis. Our results reveal an unexpectedly complex phenotypic landscape, in which mutations near the fluorophore binding pocket modulated magnesium-, potassium- and fluorophore-binding and fluorescence spectra, while distal mutations influenced structural stability and fluorescence intensity. We trained a machine learning model that accurately predicted RNA secondary structure from local epistatic interactions, despite the presence of G-quadruplexes and other noncanonical structures. Our experimental platform will facilitate the discovery and analysis of new RNA-ligand interactions.
Pan-human consensus genome significantly improves the accuracy of RNA-seq analyses

The Human Reference Genome serves as the foundation for modern genomic analyses. However, in its present form, it does not adequately represent the vast genetic diversity of the human population. In this study, we explored the consensus genome as a potential successor of the current Reference genome, and assessed its effect on the accuracy of RNA-seq read alignment. In order to find the best haploid genome representation, we constructed consensus genomes at the Pan-human, Super-population and Population levels, utilizing variant information from the 1000 Genomes project. Using personal haploid genomes as the ground truth, we compared mapping errors for real RNA-seq reads aligned to the consensus genomes versus the Reference genome. For reads overlapping homozygous variants, we found that the mapping error decreased by a factor of ~2-3 when the Reference was replaced with the Pan-human consensus genome. Interestingly, we also found that using more population-specific consensuses resulted in little to no increase over using the Pan-human consensus, suggesting a limit in the utility of incorporating more specific genomic variation. To assess the functional impact, we compared splice junction expression in the different genomes and found that the Pan-human consensus increases accuracy of splice junction quantification for hundreds of splice junctions.
Collapse transitions in protein‐like lattice polymers: The effect of sequence patterns

The collapse transition of lattice protein‐like heteropolymers has been studied by means of the Monte Carlo method. The protein model has been reduced to the α‐carbon trace restricted to a high coordination lattice. The sequences of model heteropolymers contain two types of mers: hydrophobic/nonpolar (H) and hydrophilic/polar (P). Interactions of HH and PP pairs were assumed to be negative (weaker attractions of PP pairs) while the contact energy for HP pairs was equal to zero. All sequence‐specific short‐range interactions have been neglected in the present studies. It has been found that homopolymeric chains undergo a smooth collapse transition to a dense globular state. The globule lacks any signatures of local ordering that could be interpreted as a model of protein secondary structure. Hetero‐polymers with the sequences of hydrophilic and hydrophobic residues characteristic for α‐ and β‐type proteins undergo a somewhat sharper (though continuous) collapse transition to a dense globular state with elements of local ordering controlled by the sequence. The helical pattern induces more secondary structure than the β‐type pattern. For all examined sequences the level of local ordering was lower than the average secondary structure content of globular proteins. The results are compared with other theoretical work and with known experimental facts. The implications for the reduced modeling of protein systems are briefly discussed. © 1997 John Wiley & Sons, Inc. Biopoly 42: 537–548, 1997
The Urinary Microbiome in Postmenopausal Women with Recurrent Urinary Tract Infections

Recurrent urinary tract infections (UTI) are highly prevalent in postmenopausal women, where vaginal estrogen and prophylactic antibiotics are used for treatment. The etiology of recurrent UTIs is not completely known, but the urinary microbiome is thought to be implicated. Thus, we aimed to compare the “steady state” urinary microbiome in three groups of menopausal women who were all using topically-applied vaginal estrogen: 1) women with recurrent UTIs on daily antibiotic prophylaxis; 2) women with recurrent UTIs not on antibiotic prophylaxis; and 3) age-matched controls without recurrent UTIs. Here we present a cross-sectional analysis of baseline data from 64 women enrolled in a longitudinal cohort study. Catheterized urine samples were collected > 4 weeks after last treatment for UTI. Samples were evaluated using expanded quantitative urine culture (EQUC) and 16S rRNA gene sequencing. With EQUC techniques, there were no significant differences in the median numbers of microbial species isolated among groups (p=0.96), even when considering Lactobacilli (p=0.72). However, there were trends towards different Lactobacillus species between groups. With sequencing the overwhelming majority of urinary samples contained Lactobacilli, with non-significant trends in relative abundance of Lactobacilli among groups. Using a Bayesian regression analysis for compositional data, we identified significant differences in anaerobic taxa that were associated with phenotypic groups. Most of these differences centered on Bacteroidales and the family Prevotellaceae, though differences were also noted in Actinobacteria and certain genera of Clostridiales. Associations between anaerobes within the urinary microbiome and recurrent UTI warrants further investigation. IMPORTANCE In menopausal women with recurrent urinary tract infections (UTIs) compared to those without, the abundance of Lactobacillus within the urinary microbiome is not significantly different when vaginal estrogen is regularly used. In this population, Lactobacillaceae were identified in 97% of urine samples using culture-independent techniques. However, with expanded urine cultures, women with recurrent UTIs taking daily antibiotics had a disproportionately low amount of L. gasseri/L. acidophilus compared to the other phenotypic groups. These findings support the theory that certain Lactobacillus species may be more important than others in the pathophysiology of postmenopausal recurrent UTIs. Furthermore, when using culture-independent techniques to explore urinary microbiota across phenotypic groups, we identified differences in multiple anaerobic taxa. Taken together, these results suggest that altered ratios of anaerobes and certain Lactobacillus species within the urinary microbiome may be implicated in postmenopausal recurrent UTI.
Comprehensive Analysis of hsa-miR-654-5p’s Tumor-Suppressing Functions

The pivotal roles of miRNAs in carcinogenesis, metastasis, and prognosis have been demonstrated recently in various cancers. This study intended to investigate the specific roles of hsa-miR-654-5p in lung cancer, which was poorly discussed. A series of closed-loop bioinformatic functional analyses were integrated with in vitro experimental validation to explore the overall biological functions and pan-cancer regulation pattern of miR-654-5p. We found that miR-654-5p abundance was significantly elevated in LUAD tissues and correlated with patients’ survival. 275 potential targets of miR-654-5p were then identified and the miR-654-5p-RNF8 regulation axis was validated in vitro as a proof of concept. Furthermore, we revealed the tumor-suppressing roles of miR-654-5p and demonstrated that miR-654-5p inhibited lung cancer cell epithelial-mesenchymal transition (EMT) process, cell proliferation, and migration using target-based, abundance-based and ssGSEA-based bioinformatic methods and in vitro validation. Following the construction of a protein-protein interaction network, 11 highly interconnected hub genes were identified and a five genes risk scoring model was developed to assess their potential prognostic ability. Our study will not only provide a basic miRNA-mRNA-phenotypes reference map for understanding the function of miR-654-5p in different cancers but also reveal the tumor-suppressing roles and prognostic values of miR-654-5p in lung cancer. Graphical abstract
Transcriptional analysis of immune modulatory genes in melanoma treated with PD-1 blockade

We aimed to characterize immunological features of melanoma patients treated with PD-1 blockade using tumor transcriptomic datasets. Response-dependent and response-independent predictors based on biological knowledge were investigated. Domain knowledge-driven regression-based analysis identified CEACAM1, CD40, B7-H3, and CD112 as key genes that determine the melanoma immune status. We devised the transcriptional deviance score (TDS) representing the individual sample-wise contribution to the immune network. The TDS not only showed good predictive power for immune checkpoint inhibitor (ICI) responses but also suggested specific gene interactions that determine ICI responses. Dynamic TDS changes following ICI treatment were related to long survival, indicating immune network modulation by ICIs occurred in responders. A predictive model incorporating B7-H3 and CEACAM1 expression, mutational status, clinical features, and the TDS showed excellent performance for ICI response. Thus, our approaches suggest a novel measure for the tumor immune temperature and provide insight into melanoma immunobiology. Highlights We applied outcome-independent and outcome-dependent methods to investigate melanoma immunobiology. CEACAM1, CD40, B7-H3, and CD112 expression levels are key determinants of immune status. We devised a TDS that could measure tumor immune network status at the individual level. Incorporating regression and correlation approaches greatly improves predictive power.
Multiple overlapping dynamic patterns of the visual sensory network in schizophrenia

Although visual processing impairments have been explored in schizophrenia (SZ), their underlying neurobiology of the visual processing impairments has not been widely studied. Also, while some research has hinted at differences in information transfer and flow in SZ, there are few investigations of the dynamics of functional connectivity within visual networks. In this study, we analyzed resting-state fMRI data of the visual sensory network (VSN) in 160 healthy control (HC) subjects and 151 SZ subjects. We estimated 9 independent components within the VSN. Then, we calculated the dynamic functional network connectivity (dFNC) using the Pearson correlation. Next, using k-means clustering, we partitioned the dFNCs into five distinct states, and then we calculated the portion of time each subject spent in each state, that we termed the occupancy rate (OCR). Using OCR, we compared HC with SZ subjects and investigated the link between OCR and visual learning in SZ subjects. Besides, we compared the VSN functional connectivity of SZ and HC subjects in each state. We found that this network is indeed highly dynamic. Each state represents a unique connectivity pattern of fluctuations in VSN FNC, and all states showed significant disruption in SZ. Overall, HC showed stronger connectivity within the VSN in states. SZ subjects spent more time in a state in which the connectivity between the middle temporal gyrus and other regions of VNS is highly negative. Besides, OCR in a state with strong positive connectivity between middle temporal gyrus and other regions correlated significantly with visual learning scores in SZ.
Development of functional corticomuscular connectivity from childhood to adulthood during precision grip with the dominant and non-dominant hand

How does the neural control of manual movements mature from childhood to adulthood? Here, we investigated developmental differences in functional corticomuscular connectivity using coherence techniques in 91 individuals recruited from four different age groups covering the age range 8-30y. EEG and EMG were recorded while participants performed a unimanual visual force-tracing task requiring fine control of the force produced in a precision grip with both the dominant and non-dominant hand. Using beamforming methods, we reconstructed source activity from EEG data displaying peak coherence with the active FDI muscle during the task in order to assess functional corticomuscular connectivity. Our results revealed that coherence was greater in adolescents and adults than in children and that the difference in coherence between children and adults was driven by a greater magnitude of descending (cortex-to-muscle) coherence. This was paralleled by the observation of a posterior-to-anterior shift in the cortical sources displaying corticomuscular coherence within the contralateral hemisphere from late adolescence. Finally, we observed that corticomuscular coherence was higher on the non-dominant compared to the dominant hand across age groups. These findings provide a detailed characterization of differences in task-related corticomuscular connectivity for individuals at different stages of typical ontogenetic development that may be related to the maturational refinement of dexterous motor control. Key points ‐ Fine motor control is gradually refined during human motor development, but little is known about the underlying neurophysiological mechanisms. ‐ Here, we used EEG and EMG to investigate functional corticomuscular connectivity during a precision grip with the dominant and non-dominant hand in 91 typically developed children, adolescents and adults (age range 8-30y). ‐ We show that older adolescents and adults are characterized by greater levels of corticomuscular coherence compared to children and that this is mainly driven by greater magnitudes of descending coherence (cortex-to-muscle). This is paralleled by a more anterior cortical site of coherence in older adolescents and adults compared to younger individuals. ‐ These results help us better understand the ontogenetic development of task-related functional connectivity in sensorimotor networks.