Miracle-dz/UniFAR-FTUs huggingface dataset:
In vivo stabilization of endogenous chloroplast RNAs by customized artificial pentatricopeptide repeat proteins

Pentatricopeptide repeat (PPR) proteins are helical repeat-proteins that bind RNA in a modular fashion with a sequence-specificity that can be manipulated by the use of an amino acid code. As such, PPR repeats are promising scaffolds for the design of RNA binding proteins for synthetic biology applications. However, the in vivo functional capabilities of artificial PPR proteins built from consensus PPR motifs are just starting to be explored. Here, we report in vivo functions of an artificial PPR protein, dPPRrbcL, made of consensus PPR motifs that were designed to bind a sequence near the 5’ end of rbcL transcripts in Arabidopsis chloroplasts. We used a functional complementation assay to demonstrate that this protein bound its intended RNA target with specificity in vivo and that it substituted for a natural PPR protein by stabilizing processed rbcL mRNA. We targeted a second protein of analogous design to the petL 5’ UTR, where it substituted for the native stabilizing PPR protein PGR3, albeit inefficiently. These results showed that artificial PPRs can be engineered to functionally mimic the class of native PPR proteins that serve as physical barriers against exoribonucleases.
Inferring a spatial code of cell-cell interactions across a whole animal body

Cell-cell interactions shape cellular function and ultimately organismal phenotype. Interacting cells can sense their mutual distance using combinations of ligand-receptor pairs, suggesting the existence of a spatial code, i.e., signals encoding spatial properties of cellular organization. However, this code driving and sustaining the spatial organization of cells remains to be elucidated. Here we present a computational framework to infer the spatial code underlying cell-cell interactions from the transcriptomes of the cell types across the whole body of a multicellular organism. As core of this framework, we introduce our tool cell2cell, which uses the coexpression of ligand-receptor pairs to compute the potential for intercellular interactions, and we test it across the Caenorhabditis elegans’ body. Leveraging a 3D atlas of C. elegans’ cells, we also implement a genetic algorithm to identify the ligand-receptor pairs most informative of the spatial organization of cells across the whole body. Validating the spatial code extracted with this strategy, the resulting intercellular distances are negatively correlated with the inferred cell-cell interactions. Furthermore, for selected cell-cell and ligand-receptor pairs, we experimentally confirm the communicatory behavior inferred with cell2cell and the genetic algorithm. Thus, our framework helps identify a code that predicts the spatial organization of cells across a whole-animal body.
Temporal dynamics of animacy categorization in the brain of patients with mild cognitive impairment

Electroencephalography (EEG) has been commonly used to measure brain alterations in Alzheimer’s Disease (AD). However, reported changes are limited to those obtained from using univariate measures, including activation level and frequency bands. To look beyond the activation level, we used multivariate pattern analysis (MVPA) to elicit the pattern of information processing from EEG responses to an animacy categorization task. Comparing healthy controls (HC) with patients with mild cognitive impairment (MCI), we found that the neural speed of animacy information processing is decreased in MCI patients. Moreover, we found critical time-points during which the representational pattern of animacy for MCI patients was significantly discriminable from that of HC, while the activation level remained unchanged. Together, these results suggest that the speed and pattern of animacy information processing provide clinically useful information as a potential biomarker for detecting early changes in MCI and AD patients.
Basal neural stem cells in the subventricular zone drive postnatal neurogenesis with apical stem cells acting as proliferation gate-keepers

Neural stem cells (NSCs) in the ventricular-subventricular zone (V-SVZ) contribute to olfaction by being the origin of most adult-born olfactory bulb (OB) interneurons. The current consensus maintains that adult NSCs are radial glialike progenitors apically contacting the lateral ventricle and generating intermediate progenitors migrating at the basal V-SVZ. Whether basal NSCs are present in the V-SVZ is unknown. We here used genetic tagging of NSCs in vivo and additional labelling approaches to reveal that basal NSCs lacking apical attachment represent the largest NSC type in the postnatal V-SVZ from birth onwards. Despite dividing faster than their apical counterpart, basal NSCs still undergo long-term self-renewal and quiescence. Unlike apical NSCs, they are largely devoid of primary cilia and Prominin-1, Nestin and glial fibrillary acidic protein (GFAP) immunoreactivity. Six weeks after viral tagging of apical cells, few descendant cells were detected in the basal V-SVZ, including Sox9+ progenitors and GFAP+ astrocytes, and very rare new neurons in the OB, indicating that adult-born OB neurons originate from basal and not apical NSCs. Consistent with this, we found that pregnancy, a physiological modulator of adult OB neurogenesis, selectively increases the number of basal but not apical NSCs. Lastly, we find that apical NSCs display the highest levels of Notch activation in the neural lineage, and that selective apical downregulation of Notch-signaling effector Hes1 decreases Notch activation while increasing proliferation across the V-SVZ. Thus, apical NSCs act essentially as neurogenesis gatekeepers by modulating Notch-mediated lateral inhibition of proliferation in the adult V-SVZ. Graphical Abstract Highlights Basal NSCs are the most abundant stem cell type in the adult V-SVZ from birth onwards. Apical and basal NSCs display distinct characteristics and cell cycle progression dynamics. Apical NSCs are not the main source of newly generated adult OB interneurons. Apical NSCs regulate intermediate progenitor proliferation by orchestrating Notch-mediated lateral inhibition.
Resources allocation explains the differential roles of RBS and promoter strengths in cell mass distribution and optimal protein expression productivity

Design of synthetic genetic circuits without considering the impact of host–circuit interactions results in an inefficient design process and lengthy trial-and-error iterations to appropriately tune the expression levels. Microorganisms have evolved to reach an optimal use of cellular resources. This balance is perturbed by circuit-host interactions resulting from the interaction among the cell environment from which the cell takes substrates, its metabolism, and the needs of exogenous synthetic genetic circuit introduced in the cell host. The resulting competition for common shared cell resources introduces spurious dynamics leading to problems of malfunctioning of the synthetic circuit due to lack of enough cellular resources. Therefore, there is an increasing interest in development of methods for model-based design of synthetic gene circuits considering host-circuit interactions. Here we present a small-size model of gene expression dynamics in bacterial cells accounting for host-circuit interactions. For each gene, the model defines the cellular resources recruitment strength as the key functional coefficient that allows to explain the distribution of resources among the host and the genes of interest and the relationship between the usage of resources and cell growth. This functional coefficient explicitly takes into account the availability of resources and lab-accessible gene characteristics, such as promoter and ribosome binding site (RBS) strengths and capture their interplay with the availability of free cell resources. In spite of its simplicity, the model is able to explain the differential role of promoter and RBS strengths in the distribution of protein mass and the optimal protein expression productivity with remarkable fit to the experimental data from the literature for E. coli. This makes the model amenable for model-based circuit design purposes. Moreover, the model also allows to understand why endogenous ribosomal and non-ribosomal genes have evolved different strategies in the expression space.
TimeCycle: topology inspired method for the detection of cycling transcripts in circadian time-series data

The circadian rhythm drives the oscillatory expression of thousands of genes across all tissues. The recent revolution in high-throughput transcriptomics, coupled with the significant implications of the circadian clock for human health, has sparked an interest in circadian profiling studies to discover genes under circadian control. Here we present TimeCycle: a topology-based rhythm detection method designed to identify cycling transcripts. For a given time-series, the method reconstructs the state space using time-delay embedding, a data transformation technique from dynamical systems. In the embedded space, Takens’ theorem proves that the dynamics of a rhythmic signal will exhibit circular patterns. The degree of circularity of the embedding is calculated as a persistence score using persistent homology, an algebraic method for discerning the topological features of data. By comparing the persistence scores to a bootstrapped null distribution, cycling genes are identified. Results in both synthetic and biological data highlight TimeCycle’s ability to identify cycling genes across a range of sampling schemes, number of replicates, and missing data. Comparison to competing methods highlights their relative strengths, providing guidance as to the optimal choice of cycling detection method.
Towards Personalized Auditory Models: Predicting Individual Sensorineural Hearing-Loss Profiles From Recorded Human Auditory Physiology

Over the past decades, different types of auditory models have been developed to study the functioning of normal and impaired auditory processing. Several models can simulate frequency-dependent sensorineural hearing loss (SNHL) and can in this way be used to develop personalized audio-signal processing for hearing aids. However, to determine individualized SNHL profiles, we rely on indirect and noninvasive markers of cochlear and auditory-nerve (AN) damage. Our progressive knowledge of the functional aspects of different SNHL subtypes stresses the importance of incorporating them into the simulated SNHL profile, but has at the same time complicated the task of accomplishing this on the basis of noninvasive markers. In particular, different auditory-evoked potential (AEP) types can show a different sensitivity to outer-hair-cell (OHC), inner-hair-cell (IHC), or AN damage, but it is not clear which AEP-derived metric is best suited to develop personalized auditory models. This study investigates how simulated and recorded AEPs can be used to derive individual AN- or OHC-damage patterns and personalize auditory processing models. First, we individualized the cochlear model parameters using common methods of frequency-specific OHC-damage quantification, after which we simulated AEPs for different degrees of AN damage. Using a classification technique, we determined the recorded AEP metric that best predicted the simulated individualized cochlear synaptopathy profiles. We cross-validated our method using the data set at hand, but also applied the trained classifier to recorded AEPs from a new cohort to illustrate the generalizability of the method.
TEMPORAL DYNAMICS OF IMPLICIT MEMORY 1 Temporal dynamics of implicit memory underlying serial dependence

Serial dependence is the effect in which the immediately preceding trial influences participants' responses to the current stimulus. But for how long does this bias last in the absence of interference from other stimuli? Here, we had 20 healthy young adult participants (12 women) perform a coincident timing task using different inter-trial intervals to characterize the serial dependence effect as the time between trials increases. Our results show that serial dependence abruptly decreases from 0.1 s to 1 s inter-trial interval, but it remains pronounced after that for up to 8 s. In addition, participants' response variability slightly decreases over longer intervals. We discuss these results in light of recent models suggesting that serial dependence might rely on a short-term memory trace kept through changes in synaptic weights, which might explain its long duration and apparent stability over time.
The early endosomal protein Rab21 is critical for enterocyte functions and intestinal homeostasis

Membrane trafficking is defined as the vesicular transport of molecules into, out of, and throughout the cell. In intestinal enterocytes, defects in endocytic/recycling pathways result in impaired function and are linked to genetic diseases. However, how these trafficking pathways regulate intestinal tissue homeostasis is poorly understood. Using the Drosophila intestine as an in vivo model system, we investigated enterocyte-specific functions for the early endosomal trafficking machinery in gut homeostasis. We focused on the small GTPase Rab21, which regulates specific steps in early endosomal trafficking. Rab21-depleted guts showed severe abnormalities in intestinal morphology, with deregulated homeostasis associated with a gain in mitotic cells and increased cell death. Increases in both apoptosis and yorkie signaling were responsible for compensatory proliferation and tissue inflammation. Using a RNA interference screen, we identified specific regulators of autophagy and membrane trafficking that phenocopied Rab21 loss. We further showed that Rab21-induced hyperplasia was rescued by inhibition of epidermal growth factor receptor signaling, and identified improperly trafficked cargoes in Rab21-depleted enterocytes. Our data shed light on an important role for early endosomal trafficking, and particularly Rab21, in enterocyte-mediated intestinal homeostasis.
Interactions of anti-COVID-19 drug candidates with multispecific ABC and OATP drug transporters

In the COVID-19 epidemic, several repurposed drugs have been proposed to alleviate the major health effects of the disease. These drugs are often applied together with analgesics or non-steroid anti-inflammatory compounds, and co-morbid patients may also be treated with anticancer, cholesterol-lowering or antidiabetic agents. Since drug ADME-tox properties may be significantly affected by multispecific transporters, here we examined the interactions of the repurposed drugs with the key human multidrug transporters, present in the major tissue barriers and strongly affecting pharmacokinetics. Our in vitro studies, using a variety of model systems, explored the interactions of the antimalarial agents chloroquine and hydroxychloroquine, the antihelmintic ivermectin, and the proposed antiviral compounds, ritonavir, lopinavir, favipiravir and remdesivir with the ABCB1/Pgp, ABCG2/BCRP and ABCC1/MRP1 exporters, as well as the OATP2B1 and OATP1A2 uptake transporters. The results presented here show numerous pharmacologically relevant transporter interactions and may provide a warning for the potential toxicities of these repurposed drugs, especially in drug combinations at the clinic.
Accurate genomic variant detection in single cells with primary template-directed amplification

Significance Primary template-directed amplification (PTA) is a major improvement in whole genome amplification, which is required to study intratissue cellular evolution. As presented in the manuscript, PTA produces significantly improved and reproducible genome sequencing coverage and variant detection from a single genome of a single cell. Applications of measuring genetic diversity from single cells with PTA presented in the manuscript include examining the acquisition of genetic changes during normal development and aging, measuring the consequences of specific perturbations such as genome editing, and characterizing the evolution of clonal populations during cancer formation. Improvements in whole genome amplification (WGA) would enable new types of basic and applied biomedical research, including studies of intratissue genetic diversity that require more accurate single-cell genotyping. Here, we present primary template-directed amplification (PTA), an isothermal WGA method that reproducibly captures >95% of the genomes of single cells in a more uniform and accurate manner than existing approaches, resulting in significantly improved variant calling sensitivity and precision. To illustrate the types of studies that are enabled by PTA, we developed direct measurement of environmental mutagenicity (DMEM), a tool for mapping genome-wide interactions of mutagens with single living human cells at base-pair resolution. In addition, we utilized PTA for genome-wide off-target indel and structural variant detection in cells that had undergone CRISPR-mediated genome editing, establishing the feasibility for performing single-cell evaluations of biopsies from edited tissues. The improved precision and accuracy of variant detection with PTA overcomes the current limitations of accurate WGA, which is the major obstacle to studying genetic diversity and evolution at cellular resolution.
Common allotypes of ER aminopeptidase 1 have substrate-dependent and highly variable enzymatic properties

Objective Polymorphic variation of immune system proteins can drive variability of individual immune responses. ER aminopeptidase 1 (ERAP1) generates antigenic peptides for presentation by MHC class I molecules. Coding single nucleotide polymorphisms (SNPs) in ERAP1 have been associated with predisposition to inflammatory rheumatic disease and shown to affect functional properties of the enzyme, but the interplay between combinations of these SNPs as they exist in allotypes, has not been thoroughly explored. Methods We used phased genotype data to estimate ERAP1 allotype frequency in 2,504 individuals across five major human populations, generated highly pure recombinant enzymes corresponding to the 10 most common ERAP1 allotypes and systematically characterized their in vitro enzymatic properties. Results We find that ERAP1 allotypes possess a wide range of enzymatic activities, up to 60-fold, whose ranking is substrate-dependent. Strikingly, allotype 10, previously associated with Behçet’s disease, is consistently a low-activity outlier, suggesting that a significant percentage of individuals carry a sub-active ERAP1 gene. Enzymatic analysis revealed that ERAP1 allotypes can differ in both catalytic efficiency and substrate affinity, differences that can change intermediate accumulation in multi-step trimming reactions. Alterations in efficacy of an allosteric inhibitor that targets the regulatory site suggest that allotypic variation influences the communication between the regulatory and the active site. Conclusion Our work defines the wide landscape of ERAP1 activity in human populations and demonstrates how common allotypes can induce substrate-dependent variability in antigen processing, thus contributing, in synergy with MHC haplotypes, to immune response variability and to predisposition to chronic inflammatory conditions.
Admixture can readily lead to the formation of supergenes

Supergenes are genetic architectures allowing the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of sets of alleles at linked loci is determined by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, color or behavioral polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the genetic origin of supergenes remains poorly understood. Here, using a population genetic two-island model, we explore how the evolution of overdominant chromosomal inversions may lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations connected by gene flow leads to an increase in frequency of poorly adapted, immigrant haplotypes. When inversions are associated with recessive fitness cost hampering their fixation (such as a mutational load), this results in the formation of supergenes. These results provide a realistic scenario for the evolution of supergenes and inversion polymorphisms in general, and bring new light into the importance of admixture in the formation of new genetic architectures.
Shifting gradients of macroscale cortical organization mark the transition from childhood to adolescence

Significance Here, we describe age-dependent shifts in the macroscale organization of cortex in childhood and adolescence. The characterization of functional connectivity patterns in children revealed an overarching organizational framework anchored within the unimodal cortex, between somatosensory/motor and visual regions. Conversely, in adolescents, we observed a transition into an adult-like gradient, situating the default network at the opposite end of a spectrum from primary somatosensory/motor regions. This spatial framework emerged gradually with age, reaching a sharp inflection point at the transition from childhood to adolescence. These data reveal a developmental change from a functional motif first dominated by the distinction between sensory and motor systems and then balanced through interactions with later-maturing aspects of association cortex that support more abstract cognitive functions. The transition from childhood to adolescence is marked by pronounced shifts in brain structure and function that coincide with the development of physical, cognitive, and social abilities. Prior work in adult populations has characterized the topographical organization of the cortex, revealing macroscale functional gradients that extend from unimodal (somatosensory/motor and visual) regions through the cortical association areas that underpin complex cognition in humans. However, the presence of these core functional gradients across development as well as their maturational course have yet to be established. Here, leveraging 378 resting-state functional MRI scans from 190 healthy individuals aged 6 to 17 y old, we demonstrate that the transition from childhood to adolescence is reflected in the gradual maturation of gradient patterns across the cortical sheet. In children, the overarching organizational gradient is anchored within the unimodal cortex, between somatosensory/motor and visual territories. Conversely, in adolescence, the principal gradient of connectivity transitions into an adult-like spatial framework, with the default network at the opposite end of a spectrum from primary sensory and motor regions. The observed gradient transitions are gradually refined with age, reaching a sharp inflection point in 13 and 14 y olds. Functional maturation was nonuniformly distributed across cortical networks. Unimodal networks reached their mature positions early in development, while association regions, in particular the medial prefrontal cortex, reached a later peak during adolescence. These data reveal age-dependent changes in the macroscale organization of the cortex and suggest the scheduled maturation of functional gradient patterns may be critically important for understanding how cognitive and behavioral capabilities are refined across development.