Miracle-dz/UniFAR-FTUs huggingface dataset:
Differences between apoplastic and cytosolic reactive oxygen species production in Arabidopsis during pattern-triggered immunity

Despite an ever-increasing interest in reactive oxygen species (ROS) signalling during plant-microbe interactions, very little information exists, mainly for technical reasons, on the molecular mechanisms regulating intracellular hydrogen peroxide (H2O2) signalling during PAMP-triggered immunity. Here, we used a sensitive fluorimetry method and the H2O2 sensor roGFP2-Orp1, which revealed unsuspected features on the regulation of cytoplasmic H2O2 and thiol redox dynamics upon pathogen-associated molecular patterns (PAMPs) perception by Arabidopsis thaliana. Extended PAMP-induced cytosolic roGFP2-Orp1 oxidation was distinct from the transient oxidative burst in the apoplast measured by luminol oxidation. Pharmacological and genetic analyses indicate that the prolonged PAMP-induced H2O2 increase in the cytoplasm was largely independent on NADPH oxidases and apoplastic peroxidases. By contrast, the NADPH oxidase mutant rbohF was hyper-sensitive to roGFP2-Orp1 oxidation by H2O2 and PAMP indicating a lower antioxidant capacity. Unlike previous reports, the rbohF mutant, but not rbohD, was impaired in PAMP-triggered stomatal closure and ROS production measured by a fluorescein-based probe in guard cells resulting in defects in stomatal defences against bacteria. However, stomatal closure was not correlated with an increase in roGFP2-Orp1 oxidation in guard cells. Interestingly, RBOHF also participated in PAMP-induced apoplastic alkalinisation. Altogether, our results provide novel insights on the interplay between apoplastic and cytosolic ROS dynamics and highlight the importance of RBOHF in plant immunity. Significance statement Plants mount defence responses to pathogens by detecting pathogen-associated molecular patterns (PAMPs). One response is a rapid and transient burst of reactive oxygen species (ROS, e.g. superoxide and hydrogen peroxide) in the cell wall (apoplast) produced by NADPH oxidases and cell wall peroxidases. Using a genetically-encoded hydrogen peroxide sensor roGFP2-Orp1, we found that, in contrast to the transient apoplastic ROS burst, there is also prolonged hydrogen peroxide production in the cytosol upon PAMP perception which is independent of NADPH oxidase and cell wall peroxidases. Our results suggest that apoplastic ROS rather than intracellular hydrogen peroxide is a signal triggering stomatal closure during PAMP-triggered immunity. Additionally, we re-address the relative contribution of the NADPH oxidases D and F in stomatal immunity.
Whole Brain Correlates of Individual Differences in Skin Conductance Responses during Human Fear Conditioning

Understanding the neural basis for individual differences in the conditioned skin conductance response (SCR) may inform on autonomic regulation in fear-related psychopathology. Previous region-of-interest (ROI) analyses have implicated the amygdala in regulating conditioned SCR, but whole brain analyses are lacking. This study examined correlations between individual differences in conditioned SCR and neural activity throughout the whole brain by using data from a large functional magnetic resonance imaging study (N = 285). Results show that conditioned SCR correlates with activity in the dorsal anterior cingulate cortex/anterior midcingulate cortex, anterior insula, bilateral temporoparietal junction, right frontal operculum, bilateral dorsal premotor cortex, right superior parietal lobe and midbrain. A ROI analysis replicated a correlation between amygdala activity and conditioned SCR, but amygdala contribution to SCR was modest compared with other regions. We suggest that implicated neural regions belong to a large-scale midcingulo-insular network related to salience detection and autonomic-interoceptive processing. Altered activity within this network may underlie individual differences in conditioned SCR and autonomic aspects of psychopathology.
Locomotion induced by medial septal glutamatergic neurons is linked to intrinsically generated persistent firing

Medial septal glutamatergic neurons are active during theta oscillations and locomotor activity. Prolonged optogenetic activation of medial septal glutamatergic neurons drives theta oscillations and locomotion for extended periods of time outlasting the stimulus duration. However, the cellular and circuit mechanisms supporting the maintenance of both theta oscillations and locomotion remain elusive. Specifically, it remains unclear whether the presence of theta oscillations is a necessary prerequisite for locomotion, and whether neuronal activity within the medial septum underlies its persistence. Here we show that a persistent theta oscillation can be induced by a brief transient activation of glutamatergic neurons. Moreover, persistent locomotion is initiated even if the theta oscillation is abolished by blocking synaptic transmission in the medial septum. We observe persistent spiking of medial septal neurons that outlasts the stimulus for several seconds, both in vivo and in vitro. This persistent activity is driven by intrinsic excitability of glutamatergic neurons.
Comparative anatomy and phylogenetic contribution of intracranial osseous canals and cavities in armadillos and glyptodonts (Xenarthra, Cingulata)

The evolutionary history of the Cingulata, as for many groups, remains a highly debated topic to this day, particularly for one of their most emblematic representatives: the glyptodonts. There is no consensus among morphological and molecular phylogenies relative to their position within Cingulata. As demonstrated by recent works, the study of the internal anatomy constitutes a promising path for enriching morphological matrices for the phylogenetic study of armadillos. However, internal cranial anatomy remains under-studied in the Cingulata. Here we explored and compared the anatomy of intracranial osseous canals and cavities in a diverse sample of extant and extinct cingulates, including the earliest well-preserved glyptodont crania. The virtual 3D reconstruction (using X-ray microtomography) of selected canals, i.e., the nasolacrimal canal, the palatine canal, the sphenopalatine canal, the canal for the frontal diploic vein, the transverse canal, the orbitotemporal canal, the canal for the capsuloparietal emissary vein and the posttemporal canal, and alveolar cavities related to cranial vascularization, innervation or tooth insertion allowed us to compare the locations, trajectories and shape of these structures and to discuss their potential interest for cingulate systematics. We tentatively reconstructed evolutionary scenarios for eight selected traits on these structures, in which glyptodonts often showed a greater resemblance to pampatheres, to the genus Proeutatus and/or to chlamyphorines. This latter pattern was partly congruent with recent molecular hypotheses, but more research is needed on these resemblances and on the potential effects of development and allometry on the observed variations. Overall, these comparisons have enabled us to highlight new anatomical variation that may be of great interest to further explore the evolutionary history of cingulates and the origins of glyptodonts on a morphological basis.
The effective population size modulates the strength of GC biased gene conversion in two passerines

Understanding the determinants of genomic base composition is fundamental to understanding genome evolution. GC biased gene conversion (gBGC) is a key driving force behind genomic GC content, through the preferential incorporation of GC alleles over AT alleles during recombination, driving them towards fixation. The majority of work on gBGC has focussed on its role in coding regions, largely to address how it confounds estimates of selection. Non-coding regions have received less attention, particularly in regard to the interaction of gBGC and the effective population size (Ne) within and between species. To address this, we investigate how the strength of gBGC (B = 4Neb, where b is the conversion bias) varies within the non-coding genome of two wild passerines. We use a dataset of published high coverage genomes (10 great tits and 10 zebra finches) to estimate B, nucleotide diversity, changes in Ne, and crossover rates from linkage maps, in 1Mb homologous windows in each species. We demonstrate remarkable conservation of both B and crossover rate between species. We show that the mean strength of gBGC in the zebra finch is more than double that in the great tit, consistent with its twofold greater effective population size. B also correlates with both crossover rate and nucleotide diversity in each species. Finally, we estimate equilibrium GC content from both divergence and polymorphism data, which indicates that B has been increasing in both species, and provide support for population expansion explaining a large proportion of this increase in the zebra finch. Significance statement Understanding the forces that change the nucleotide base composition of genomes is central to understanding their evolution. One such force is GC biased gene conversion, a process that during recombination converts some heterozygous base positions to homozygous. This process is more likely to convert adenine and thymine bases to guanine and cytosine bases than the other way around, hence is GC biased. This increases the frequency of GC alleles in a way similar to positive selection. This process has largely been studied within protein coding regions, and not often compared between species. We measure its strength in the non-coding areas of the genomes of two bird species, showing it to be stronger in the species with the larger population size.
Deep learning identifies pathological abnormalities predictive of graft loss in kidney transplant biopsies

Background Interstitial fibrosis, tubular atrophy, and inflammation are major contributors to renal allograft failure. Here we seek an objective, quantitative pathological assessment of these lesions to improve predictive utility. Methods We constructed a deep-learning-based pipeline recognizing normal vs. abnormal kidney tissue compartments and mononuclear leukocyte (MNL) infiltrates from Periodic acid-Schiff (PAS) stained slides of transplant biopsies (training: n=60, testing: n=33) that quantified pathological lesions specific for interstitium, tubules and MNL infiltration. The pipeline was applied to 789 whole slide images (WSI) from baseline (n=478, pre-implantation) and 12-month post-transplant (n=311) protocol biopsies in two independent cohorts (GoCAR: 404 patients, AUSCAD: 212 patients) of transplant recipients to correlate composite lesion features with graft loss. Results Our model accurately recognized kidney tissue compartments and MNLs. The digital features significantly correlated with Banff scores, but were more sensitive to subtle pathological changes below the thresholds in Banff scores. The Interstitial and Tubular Abnormality Score (ITAS) in baseline samples was highly predictive of 1-year graft loss (p=2.8e-05), while a Composite Damage Score (CDS) in 12-month post-transplant protocol biopsies predicted later graft loss (p=7.3e-05). ITAS and CDS outperformed Banff scores or clinical predictors with superior graft loss prediction accuracy. High/intermediate risk groups stratified by ITAS or CDS also demonstrated significantly higher incidence of eGFR decline and subsequent graft damage. Conclusions This deep-learning approach accurately detected and quantified pathological lesions from baseline or post-transplant biopsies, and demonstrated superior ability for prediction of posttransplant graft loss with potential application as a prevention, risk stratification or monitoring tool.
GC-biased gene conversion in X-chromosome palindromes conserved in human, chimpanzee, and rhesus macaque

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich than the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that additional factors are required to explain elevated GC content in palindrome arms. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion, and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.
Dynamics of a hepatocyte-cholangiocyte decision-making gene regulatory network during liver development and regeneration

Liver is one of the few organs with immense regenerative potential even at adulthood in mammals. It is composed of primarily two cell types: hepatocytes and cholangiocytes, that can trans-differentiate to one another either directly or through intermediate progenitor states, contributing to remarkable regenerative potential of the liver. However, the dynamical features of decision-making between these cell-fates during liver development and regeneration remains elusive. Here, we identify a core gene regulatory network comprising c/EBPα, TGFBR2 and SOX9 that underlies liver development and injury-induced reprogramming. Dynamic simulations for this network reveal its multistable nature, enabling three distinct cell states – hepatocytes, cholangiocytes and liver progenitor cells (hepatoblasts/oval cells) – and stochastic switching among them. Predicted expression signature for these three states are validated through multiple bulk and single-cell transcriptomic datasets collected across developmental stages and injury-induced liver repair. This network can also explain the experimentally observed spatial organisation of phenotypes in liver parenchyma and predict strategies for efficient cellular reprogramming among these cell-fates. Our analysis elucidates how the emergent multistable dynamics of underlying gene regulatory networks drive diverse cell-state decisions in liver development and regeneration.
Exploring the genomic diversity and antimicrobial susceptibility of Bifidobacterium pseudocatenulatum in the Vietnamese population to aid probiotic design

Bifidobacterium pseudocatenulatum is a member of the human gut microbiota, and has previously been used as a probiotic to improve gut integrity and reduce inflammatory responses. We showed previously that B. pseudocatenulatum was significantly depleted during dysenteric diarrhea, suggesting the organism may aid in recovery from diarrhea. Here, in order to investigate its probiotic potential, we aimed to assess the genomic diversity and predicted metabolic profiles of B. pseudocatenulatum found colonizing the gut of healthy Vietnamese adults and children. We found that the population of B. pseudocatenulatum from each individual was distinct and highly diverse, with intra-clonal variation attributed to gain or loss of carbohydrate utilizing enzymes. The B. pseudocatenulatum genomes were enriched with glycosyl hydrolases that target plant-based non-digestible carbohydrates (GH13, GH43), but not host-derived glycans. Notably, the exopolysaccharide biosynthesis region from organisms isolated from healthy children showed greater genetic diversity, and was subject to a high degree of genetic modification. Antimicrobial susceptibility testing revealed that the Vietnamese B. pseudocatenulatum were uniformly susceptible to beta-lactams, but exhibited variable resistance to azithromycin, tetracycline, ciprofloxacin and metronidazole. The genomic presence of ermX and tet variants conferred resistance against azithromycin and tetracycline, respectively; ciprofloxacin resistance was associated with mutation(s) in the quinolone resistance determining region (GyrA, S115 and/or D119). Our work provides the first detailed genomic and antimicrobial resistance characterization of B. pseudocatenulatum found in the Vietnamese population, which could inform the next phase of rational probiotic design. Importance Bifidobacterium pseudocatenulatum is a probiotic candidate with potential applications in several health conditions, but its efficacy is largely strain-dependent and associated with distinct genomic and biochemical features. However, most commercial probiotics have been developed by Western institutions, which may not have ideal efficacy when administered in developing countries. This study taps into the underexplored diversity of the organism in Vietnam, and provides more understanding to its lifestyles and antimicrobial susceptibility. These data are key for selecting an optimal probiotic candidate, from our established collection, for downstream investigations and validation. Thus, our work represents a model in identifying and characterizing bespoke probiotics from an indigenous population in a developing setting.
GPT-2’s activations predict the degree of semantic comprehension in the human brain

Language transformers, like GPT-2, have demonstrated remarkable abilities to process text, and now constitute the backbone of deep translation, summarization and dialogue algorithms. However, whether these models encode information that relates to human comprehension remains controversial. Here, we show that the representations of GPT-2 not only map onto the brain responses to spoken stories, but also predict the extent to which subjects understand narratives. To this end, we analyze 101 subjects recorded with functional Magnetic Resonance Imaging while listening to 70 min of short stories. We then fit a linear model to predict brain activity from GPT-2’s activations, and correlate this mapping with subjects’ comprehension scores as assessed for each story. The results show that GPT-2’s brain predictions significantly correlate with semantic comprehension. These effects are bilaterally distributed in the language network and peak with a correlation of R=0.50 in the angular gyrus. Overall, this study paves the way to model narrative comprehension in the brain through the lens of modern language algorithms.
Unlocking inaccessible historical genomes preserved in formalin

Background Museum specimens represent an unparalleled record of historical genomic data. However, the wide-spread practice of formalin preservation has thus far impeded genomic analysis of a large proportion of specimens. Limited DNA sequencing from formalin-preserved specimens has yielded low genomic coverage with unpredictable success. We set out to refine sample processing methods and to identify specimen characteristics predictive of sequencing success. With a set of taxonomically diverse specimens collected between 1936 and 2015 and ranging in preservation quality, we compared the efficacy of several end-to-end whole genome sequencing workflows alongside a k-mer-based trimming-free read alignment approach to maximize mapping of endogenous sequence. Results We recovered complete mitochondrial genomes and up to 3X nuclear genome coverage from formalin-fixed tissues. Hot alkaline lysis coupled with phenol-chloroform extraction out- performed proteinase K digestion in recovering DNA, while library preparation method had little impact on sequencing success. The strongest predictor of DNA yield was overall specimen condition, which additively interacts with preservation conditions to accelerate DNA degradation. Conclusions We demonstrate a significant advance in capability beyond limited recovery of a small number of loci via PCR or target-capture sequencing. To facilitate strategic selection of suitable specimens for genomic sequencing, we present a decision-making framework that utilizes independent and non-destructive assessment criteria. Sequencing of formalin-fixed specimens will contribute to a greater understanding of temporal trends in genetic adaptation, including those associated with a changing climate. Our work enhances the value of museum collections worldwide by unlocking genomes of specimens that have been disregarded as a valid molecular resource.
Evolutionary rates of testes-expressed genes differ in monogamous and promiscuous Peromyscus species

Reproductive proteins, including those expressed in the testes, are among the fastest evolving proteins across the tree of life. Sexual selection on traits involved in sperm competition is thought to be a primary driver of testes gene evolution and is expected to differ between promiscuous and monogamous species due to intense competition between males to fertilize females in promiscuous lineages and lack thereof in monogamous ones. Here, we employ the rodent genus Peromyscus as a model to explore differences in evolutionary rates between testes-expressed genes in monogamous and promiscuous species. We find candidate genes that may be associated with increased sperm production in promiscuous species and gene ontology categories that show patterns of molecular convergence associated with phenotypic convergence in independently evolved monogamous species. Overall, our results highlight the possible molecular consequences of differences in mating system, likely due to differences in selective pressures.
Connectomic assessment of injury burden and longitudinal structural network alterations in moderate‐to‐severe traumatic brain injury

Traumatic brain injury (TBI) is a major public health problem. Caused by external mechanical forces, a major characteristic of TBI is the shearing of axons across the white matter, which causes structural connectivity disruptions between brain regions. This diffuse injury leads to cognitive deficits, frequently requiring rehabilitation. Heterogeneity is another characteristic of TBI as severity and cognitive sequelae of the disease have a wide variation across patients, posing a big challenge for treatment. Thus, measures assessing network-wide structural connectivity disruptions in TBI are necessary to quantify injury burden of individuals, which would help in achieving personalized treatment, patient monitoring, and rehabilitation planning. Despite TBI being a disconnectivity syndrome, connectomic assessment of structural disconnectivity has been very scarce. In this study, we propose a novel connectomic measure that we call network anomaly score (NAS) to capture the integrity of structural connectivity in TBI patients by leveraging two major characteristics of the disease: diffuseness of axonal injury and heterogeneity of the disease. Over a longitudinal cohort of moderate-to-severe TBI patients, we demonstrate that structural network topology of patients are more heterogeneous and are significantly different than that of healthy controls at 3 months post-injury, where dissimilarity further increases up to 12 months. We also show that NAS captures injury burden as quantified by post-traumatic amnesia and that alterations in the structural brain network is not related to cognitive recovery. Finally we compare NAS to major graph theory measures used in TBI literature and demonstrate the superiority of NAS in characterizing the disease.