title: Megaherbivores modify forest structure and increase carbon stocks through multiple pathways

abstract: Megaherbivores have pervasive ecological effects. In African rainforests, elephants can increase aboveground carbon, though the mechanisms are unclear. Here we combine a large unpublished dataset of forest elephant feeding with published browsing preferences totaling > 120,000 records covering 700 plant species, including nutritional data for 102 species. Elephants increase carbon stocks by: 1) promoting high wood density tree species via preferential browsing on leaves from low wood density species, which are more digestible; 2) dispersing seeds of trees that are relatively large and have the highest average wood density among tree guilds based on dispersal mode. Loss of forest elephants could cause a 5-12% decline in carbon stocks due to regeneration failure of elephant-dispersed trees and an increase in abundance of low wood density trees. These results show the major importance of megaherbivores in maintaining diverse, high-carbon tropical forests. Successful elephant conservation will contribute to climate mitigation at a scale of global relevance.


title: Treated like dirt: Robust forensic and ecological inferences from soil eDNA after challenging sample storage

abstract:
Biodiversity of soil microbiota is routinely assessed with environmental DNA-based methods, among which amplification and massive parallel sequencing of marker genes (eDNA metabarcoding) is the most common. Soil microbiota may for example be investigated in relation to biodiversity research or as a tool in forensic investigations.

After sampling, the taxonomic composition of soil biotic communities may change. In order to minimize community changes after sampling, it is desirable to reduce biological activity, e.g. by freezing immediately after sampling. However, this may be impossible due to remoteness of study sites or, in forensic cases, where soil has been attached to a questioned item for protracted periods of time.

Here we investigated the effect of storage duration and conditions on the assessment of the soil biota with eDNA metabarcoding. We extracted eDNA from freshly collected soil samples and again from the same samples after storage under contrasting temperature conditions.

We used five different primer sets targeting bacteria, fungi, protists (cercozoans), general eukaryotes, and plants. For these groups, we quantified differences in richness, evenness and community composition. Subsequently, we tested whether we could correctly infer habitat type and original sample identity after storage using a large reference dataset.

We found increased community composition differences with extended storage time and with higher storage temperature. However, for samples stored less than 28 days at a maximum of 20{degrees}C, changes were generally insignificant. Classification models could successfully assign most stored samples to their exact location of origin and correct habitat type even after weeks of storage. Even samples showing larger compositional changes generally retained the original sample as the best match (relative similarity).

Our results show that for most biodiversity and forensic applications, storage of samples for days and even several weeks may not be a problem, if storage temperature does not exceed 20{degrees}C. Even after suboptimal storage conditions, significant patterns can be reproduced.

category: ecology


title: Recombination suppression and selection affect local ancestries in genomes of a migratory songbird

abstract: The patterns of genetic relatedness among individuals vary along the genome, representing fluctuation of local ancestry. The factors responsible for this variation have not been well studied in wild animals with ecological and behavioural relevance. Here, we characterise the genomic architecture of genetic relatedness in the Eurasian blackcap, an iconic songbird species in ecology and quantitative genetics of migratory behaviour. We identify 23 genomic regions with deviated local relatedness patterns, using a chromosome-level de novo assembly of the blackcap genome and whole-genome resequencing data of 179 individuals from nine populations with diverse migratory phenotypes. Five genomic regions show local relatedness patterns of polymorphic inversions, three of which are syntenic to polymorphic inversions known in the zebra finch. Phylogenetic analysis reveals these three polymorphic inversions evolved independently in the blackcap and zebra finch indicating convergence of polymorphic inversions. Population genetic analyses in these three inversions in the blackcap suggest balancing selection between two haplotypes in one locus and background selection in the other two loci. One genomic region with deviated local relatedness is under selection against gene flow by population-specific reduction in recombination rate. Other genomic islands including 11 pericentromeric regions consist of evolutionarily conserved and non-conserved recombination cold-spots under background selection. Two of these regions with non-conserved recombination suppression are known to be associated with population-specific migratory phenotypes, where local relatedness patterns support additional effects of population-specific selection. These results highlight how different forms of recombination suppression and selection jointly affect heterogeneous genomic landscape of local ancestries.

category: evolutionary biology


title: A novel zebrafish model of SPEG-related centronuclear myopathy (CNM): characterization and comparison with other CNM model zebrafish

abstract:
Centronuclear myopathy (CNM) is a congenital neuromuscular disorder caused by pathogenic variation in genes associated with membrane trafficking and excitation-contraction coupling (ECC). Bi-allelic autosomal recessive mutations in striated muscle enriched protein kinase (SPEG) account for a subset of CNM patients. Previous research has been limited by the perinatal lethality of Speg knockout mice. Thus, the precise biological role of SPEG in skeletal muscle remains unknown. To address this issue, we generated zebrafish spega, spegb, and spega/spegb (speg-DKO) mutant lines. We demonstrate that speg-DKO zebrafish faithfully recapitulate multiple phenotypes associated with human CNM, including disruption of the ECC protein machinery, dysregulation of calcium homeostasis during ECC, and impairment of muscle performance. Taking advantage of the availability of zebrafish models of multiple CNM genetic subtypes, we compared novel and known disease markers in speg-DKO with mtm1-KO and DNM2-S619L transgenic zebrafish. We observed desmin (DES) accumulation common to all CNM subtypes, and DNM2 upregulation in muscle of both speg-DKO and mtm1-KO zebrafish. In all, we establish a new model of SPEG-related CNM, and identify abnormalities in this model suitable for defining disease pathomechanisms and evaluating potential therapies.

Summary StatementWe created a novel zebrafish speg mutant model of centronuclear myopathy that recapitulates key features of the human disorder and provides insight into pathomechanisms of the disease.

category: genetics


title: The antique genetic plight of the Mediterranean monk seal (Monachus monachus)

abstract: Disentangling the impact of Late Quaternary climate change from human activities can have crucial implications on the conservation of endangered species. We investigated the population genetics and demography of the Mediterranean monk seal (Monachus monachus), one of the worlds most endangered marine mammals, through an unprecedented dataset encompassing historical (extinct) and extant populations from the eastern North Atlantic to the entire Mediterranean Basin. We show that Cabo Blanco (Western Sahara/Mauritania), Madeira, Western Mediterranean (historical range), and Eastern Mediterranean regions segregate into four populations. This structure is likely the consequence of recent drift, combined with long-term isolation by distance (R2 = 0.7), resulting from prevailing short-distance (< 500 km) and infrequent long-distance dispersal (< 1,500 km). All populations (Madeira especially), show high levels of inbreeding and low levels of genetic diversity, seemingly declining since historical time, but surprisingly not being impacted by the 1997 massive die-off in Cabo Blanco. Approximate Bayesian Computation analyses support scenarios combining local extinctions and a major effective population size decline in all populations during Antiquity. Our results suggest that the early densification of human populations around the Mediterranean Basin coupled with the development of seafaring techniques were the main drivers of the decline of Mediterranean monk seals.


title: Heterogeneity and genomic loci of ubiquitous Cre reporter transgenes in zebrafish

abstract: The most-common strategy for zebrafish Cre/lox-mediated lineage labeling experiments combines ubiquitously expressed, lox-based Switch reporter transgenes with tissue-specific Cre or 4-OH-Tamoxifen-inducible CreERT2 driver lines. Although numerous Cre driver lines have been produced, only a few broadly expressed Switch reporters exist in zebrafish and their generation by random transgene integration has been challenging due to position-effect sensitivity of the lox-flanked recombination cassettes. Here, we compare commonly used Switch reporter lines for their recombination efficiency and reporter expression pattern during zebrafish development. Using different experimental setups, we show that ubi:Switch and hsp70l:Switch outperform current generations of two additional Switch reporters due to favorable transgene integration sites. Our comparisons also document preferential Cre-dependent recombination of ubi:Switch and hsp70l:Switch in distinct zebrafish tissues at early developmental stages. To investigate what genomic features may influence Cre accessibility and lox recombination efficiency in highly functional Switch lines, we mapped these transgenes and charted chromatin dynamics at their integration sites. Our data documents the heterogeneity among lox-based Switch transgenes towards informing suitable transgene selection for lineage labeling experiments. Our work further proposes that ubi:Switch and hsp70l:Switch define genomic integration sites suitable for universal transgene or switch reporter knock-in in zebrafish.


title: The first sheep graph pan-genome reveals the spectrum of structural variations and their effects on different tail phenotypes

abstract: Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome.

category: genomics


title: The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes


title: Dead infant carrying by chimpanzee mothers in the Budongo Forest

abstract: It has been suggested that non-human primates can respond to deceased conspecifics in ways that suggest they experience psychological states not unlike humans, some of which could indicate they exhibit a notion of death. Here, we report long-term demographic data from two East African chimpanzee groups. During a combined 40-year observation period we recorded 191 births of which 68 died in infancy, mostly within the first year. We documented the post-mortem behaviour of the mothers and found that Budongo chimpanzee mothers carried infants after their death on nine occasions, usually until the body started to decompose after a few days. However, we also observed three cases of extended carrying lasting for more than two weeks, one of which was followed by the unusual extended carrying of an object and another which lasted three months. In each case, the corpses mummified. In addition, we report four instances of recurring dead infant carrying by mothers, three of whom carried the corpse for longer during the second instance. We discuss these observations in view of functional hypotheses of dead infant carrying in primates and the potential proximate mechanisms involved in this behaviour.

category: animal behavior and cognition


title: Failure to account for behavioral variability overestimates bonobo populations and compromises accuracy in population monitoring

abstract: Indirect wildlife population surveying largely depends upon counts of artefacts of behavior (e.g., nests or dung). Likelihood to encounter these artefacts is derived from both artefact production and decay, and variability in artefact production behavior is considered to contribute minimally to inaccuracy in wildlife estimation. Here we demonstrate how ignoring behavioral variability contributes to significant population misestimation, using an example of an endangered ape, the bonobo (Pan paniscus). Until now, a single estimate of (i.e., nest construction) has been used to extrapolate bonobo densities, assumed to be broadly representative of bonobo sign production behavior. We estimated nest construction rate across seasons and social groups at the Kokolopori Bonobo Reserve, DRC and find that nest construction rates in bonobos to be highly variable across populations as well as seasonal. Failure to account for this variability led to potentially severe degradation in accuracy of bonobo population estimates of abundance, accounting for a likely overestimation of bonobo numbers by 34%, and at worst as high as 80%. With this example, we demonstrate that failure to account for inter- and intra-population behavioral variation compromises the ability to estimate both relative and absolute wildlife abundances. We argue that variation in sign production is but one of several potential ways that behavioral variability can affect conservation monitoring, should be measured across contexts whenever possible, and must be considered in population estimation confidence intervals. With increasing attention to behavioral variability as a potential tool for conservation, conservationists must also account for the impact that behavioral variability can play upon wildlife population estimation. Our results underline the importance of observational research to wildlife monitoring schemes as a critical but overlooked component of conservation management. We extensively discuss the various avenues through which behavioral variability is likely to impact wildlife monitoring accuracy and precision and propose potential avenues for accounting for behavioral variability in wildlife monitoring.


title: Failure to account for behavioral variability significantly compromises accuracy in indirect population monitoring


title: Cultural diffusion dynamics depend on behavioural production rules

abstract: Culture is an outcome of the acquisition of knowledge about behaviour through social transmission, and its subsequent production. Transmission and production are often discussed interchangeably or modeled separately, yet to date, no study has accounted for both processes and explored their interaction. We present a generative model that integrates the two in order to explore how variation in either might shape cultural diffusion dynamics. Agents make behavioural choices that change as they learn from their behavioural productions. Their repertoires also change over time, and the social transmission of behaviours depends on their frequency. We diffuse a novel behaviour through social networks across a large parameter space to demonstrate how accounting for both transmission and production reveals dependencies between individual-level behavioural production rules and population-level diffusion dynamics. We then investigate how such dependencies might affect the performance of two commonly used inferential models for social learning; Network-based Diffusion Analysis (NBDA), and Experienced Weighted Attraction models (EWA). By clarifying the distinction between acquisition and usage, we illuminate often-overlooked theoretical differences between social learning and social influence. These distinctions yield consequences and new considerations for how inferential methods are applied to empirical studies of culture.


title: Epilipidomics platform for holistic profiling of oxidized complex lipids in blood plasma of obese individuals

abstract: Lipids are a structurally diverse class of biomolecules which can undergo a variety of chemical modifications. Among them, lipid (per)oxidation attracts most of the attention due to its significance in regulation of inflammation, cell proliferation and death programs. Despite their apparent regulatory significance, the molecular repertoire of oxidized lipids remains largely elusive as accurate annotation of lipid modifications is challenged by their low abundance and largely unknown, biological context-dependent structural diversity. Here we provide a holistic workflow based on the combination of bioinformatics and LC-MS/MS technologies to support identification and relative quantification of oxidized complex lipids in a modification type- and position-specific manner. The developed methodology was used to identify epilipidomics signatures of lean and obese individuals with and without type II diabetes. Characteristic signature of lipid modifications in lean individuals, dominated by the presence of modified octadecanoid acyl chains in phospho- and neutral lipids, was drastically shifted towards lipid peroxidation-driven accumulation of oxidized eicosanoids, suggesting significant alteration of endocrine signalling by oxidized lipids in metabolic disorders.

category: biochemistry


title: Functional design of bacterial superoxide:quinone oxidoreductase

abstract: The superoxide anion - molecular oxygen reduced by a single electron - is produced in large amounts by enzymatic and adventitious reactions and can perform a range of cellular functions, including bacterial warfare and iron uptake, signalling and host immune response in eukaryotes. However, it also serves as precursor for more deleterious species such as the hydroxyl anion or peroxynitrite and therefore, cellular defense mechanisms for superoxide neutralization have evolved. In addition to the soluble proteins superoxide dismutase and superoxide reductase, recently the membrane embedded diheme cytochrome b561(CybB) from E. coli has been proposed to act as a superoxide:quinone oxidoreductase. Here, we confirm superoxide and cellular ubiquinones or menaquinones as natural substrates and show that quinone binding to the enzyme accelerates the reaction with superoxide. The reactivity of the substrates is in accordance with the here determined midpoint potential of the two b hemes (+48 and -23 mV / NHE). Our data suggest that the enzyme can work near the diffusion limit in the forward direction and can also catalyse the reverse reaction efficiently under physiological conditions. The data is discussed in context of described cytochrome b561 proteins and potential physiological roles of CybB.


title: Structure and mechanism of a cytosine transport protein

abstract: CodB is a cytosine transporter from the Nucleobase-Cation-Symport-1 (NCS1) transporter family, a member of the widespread LeuT superfamily. Previous experiments with the nosocomial pathogen Pseudomonas aeruginosa have shown CodB also to be important in the uptake of 5-fluorocytosine, which has been suggested as a novel drug to combat antimicrobial resistance by suppressing virulence in the organism. Here we solve the crystal structure of CodB from Proteus vulgaris, at 2.4[A] resolution in complex with cytosine. We show that the protein carries out the sodium-dependent uptake of cytosine and can bind 5-fluorocytosine. Comparison of the substrate-bound structures of CodB and the hydantoin transporter Mhp1, the only other NCS1 family member for which the structure is known, highlight the importance of the hydrogen bonds that the substrates make with the main chain at the breakpoint in the discontinuous helix, TM6. In contrast to other LeuT superfamily members, neither CodB nor Mhp1 make specific interactions with residues on TM1. Comparison of the structures provides insight into the intricate mechanisms of how these proteins transport substrates across the plasma membrane.


title: Stain-free nucleus identification in holographic learning flow cyto-tomography

abstract: Quantitative Phase Imaging (QPI) has gained popularity because it can avoid the staining step, which in some cases is difficult or impossible. However, QPI does not provide the well-known specificity to various parts of the cell (e.g., organelles, membrane). Here we show a novel computational segmentation method based on statistical inference that bridges the gap between the specificity of Fluorescence Microscopy (FM) and the label-free property of QPI techniques to identify the cell nucleus. We demonstrate application to stain-free cells reconstructed through the holographic learning and in flow cyto-tomography modality. In particular, by means of numerical simulations and two cancer cell lines, we demonstrate that the nucleus-like regions can be accurately distinguished within the stain-free tomograms. We show that our experimental results are consistent with confocal FM data and microfluidic cytofluorimeter outputs. This is a significant step towards extracting the three-dimensional (3D) intracellular specificity directly from the phase-contrast data in a typical flow cytometry configuration.

category: bioengineering


title: bulkAnalyseR: An accessible, interactive pipeline for analysing and sharing bulk sequencing results

abstract:
Bulk sequencing experiments (single- and multi-omics) are essential for exploring wide-ranging biological questions. To facilitate interactive, exploratory tasks, coupled with the sharing of easily accessible information, we present bulkAnalyseR, a package integrating state-of-the-art approaches using an expression matrix as the starting point (pre-processing functions are available as part of the package). Static summary images are replaced with interactive panels illustrating quality-checking, differential expression analysis (with noise detection) and biological interpretation (enrichment analyses, identification expression patterns, followed by inference and comparison of regulatory interactions). bulkAnalyseR can handle different modalities, facilitating robust integration and comparison of cis-, trans- and customised regulatory networks.

bulkAnalyseR is available on CRAN and GitHub, with extensive documentation and usage examples (https://github.com/Core-Bioinformatics/bulkAnalyseR, https://cran.r-project.org/web/packages/bulkAnalyseR/)

ContactIrina Mohorianu iim22@cam.ac.uk

category: bioinformatics


title: bulkAnalyseR: An accessible, interactive pipeline for analysing and sharing bulk multi-modal sequencing data


title: ePat: extended PROVEAN annotation tool

abstract: The  ePat (extended PROVEAN annotation tool) is a software tool that extends the functionality of PROVEAN: a software tool for predicting whether amino acid substitutions and indels will affect the biological function of proteins. The  ePat extends the conventional PROVEAN to enable the following two things, which the conventional PROVEAN could not calculate the pathogenicity of these variants. First is to calculate and score the pathogenicity of indel mutations with frameshift and variants near splice junctions. Second is to use batch processing to calculate the pathogenicity of multiple variants into a variants list (VCF file) in a single step. ePat can help extract variants that affect biological functions by utilizing not only point mutations, and indel mutations that does not cause frameshift, but also frameshift, stop gain, and splice variants. These extended features will increase detection rate and improve diagnostic of inherited diseases or associate specific variant to phenotype.


title: A Deep Learning Approach for NeoAG-Specific Prediction Considering Both HLA-Peptide Binding and Immunogenicity: Finding Neoantigens to Making T-Cell Products More Personal

abstract:
BackgroundAn emerging type of cancer treatment, known as cell immunotherapy, is gaining popularity over chemotherapy or other radiation therapy that causes mass destruction to our body. One favourable approach in cell immunotherapy is the use of neoantigens as targets that help our body immune system identify the cancer cells from healthy cells. Neoantigens, which are non-autologous proteins with individual specificity, are generated by non-synonymous mutations in the tumor cell genome. Owing to its strong immunogenicity and lack of expression in normal tissues, it is now an important target for tumor immunotherapy. Neoantigens are some form of special protein fragments excreted as a by-product on the surface of cancer cells during the DNA mutation at the tumour. In cancer immunotherapies, certain neoantigens which exist only on cancer cells elicit our white blood cells (bodys defender, anti-cancer T-cell) responses that fight the cancer cells while leaving healthy cells alone. Personalized cancer vaccines therefore can be designed de novo for each individual patient, when the specific neoantigens are found to be relevant to his/her tumour. The vaccine which is usually coded in synthetic long peptides, RNA or DNA representing the neoantigens trigger an immune response in the body to destroy the cancer cells (tumour). The specific neoantigens can be found by a complex process of biopsy and genome sequencing. Alternatively, modern technologies nowadays tap on AI to predict the right neoantigen candidates using algorithms. However, determining the binding and non-binding of neoantigens on T-cell receptors (TCR) is a challenging computational task due to its very large search space.

ObjectiveTo enhance the efficiency and accuracy of traditional deep learning tools, for serving the same purpose of finding potential responsiveness to immunotherapy through correctly predicted neoantigens. It is known that deep learning is possible to explore which novel neoantigens bind to T-cell receptors and which ones dont. The exploration may be technically expensive and time-consuming since deep learning is an inherently computational method. one can use putative neoantigen peptide sequences to guide personalized cancer vaccines design.

MethodsThese models all proceed through complex feature engineering, including feature extraction, dimension reduction and so on. In this study, we derived 4 features to facilitate prediction and classification of 4 HLA-peptide binding namely AAC and DC from the global sequence, and the LAAC and LDC from the local sequence information. Based on the patterns of sequence formation, a nested structure of bidirectional long-short term memory neural network called local information module is used to extract context-based features around every residue. Another bilstm network layer called global information module is introduced above local information module layer to integrate context-based features of all residues in the same HLA-peptide binding chain, thereby involving inter-residue relationships in the training process. introduced

ResultsFinally, a more effective model is obtained by fusing the above two modules and 4 features matric, the method performs significantly better than previous prediction schemes, whose overall r-square increased to 0.0125 and 0.1064 on train and increased to 0.0782 and 0.2926 on test datasets. The RMSE for our proposed models trained decreased to approximately 0.0745 and 1.1034, respectively, and decreased to 0.6712 and 1.6506 on test dataset.

ConclusionOur work has been actively refining a machine-learning model to improve neoantigen identification and predictions with the determinants for Neoantigen identification. The final experimental results show that our method is more effective than existing methods for predicting peptide types, which can help laboratory researchers to identify the type of novel HLA-peptide binding.


title: In situ distance measurements in a membrane transporter using maleimide functionalized orthogonal spin labels and 5-pulse electron double resonance spectroscopy

abstract: Spectroscopic investigation of membrane proteins in their native environment is a challenging task. Earlier we demonstrated the feasibility to measure precise distances within outer membrane proteins in E. coli and native membranes using methanethiosulfonate (MTS) functionalized labels combined with pulsed electron double resonance spectroscopy. Here we show the application of maleimide functionalized Gd(III), nitroxide, and trityl labels for in situ distance measurement using the cobalamin transporter BtuB. These labels enabled distance measurements for BtuB in E. coli and native outer membranes and in the membranes maleimide-Gd-DOTA also is effective. Further, we show that the observable dipolar evolution time can be significantly prolonged in the native environments using the Carr-Purcell 5-pulse electron double resonance sequence. For a nitroxide-nitroxide pair, application of sech/tanh inversion pulses substantially suppressed the 4-pulse artifact at the Q-band frequency. In the case of a nitroxide-trityl pair, Gaussian pump pulses of varying amplitude are sufficient to suppress the artifact to the typical noise level. The feasibility of a range of bioresistant spin labels and the 5-pulse electron double resonance offers promising tools for investigating heterooligomeric membrane protein complexes in their native environment.

category: biophysics


title: Multi-omics Topic Modeling for Breast Cancer Classification

abstract:
The integration of transcriptional data with other layers of information, such as the post-transcriptional regulation mediated by microRNAs, can be crucial to identify the driver genes and the subtypes of complex and heterogeneous diseases such as cancer. This paper presents an approach based on topic modeling to accomplish this integration task. More specifically, we show how an algorithm based on a hierarchical version of stochastic block modeling can be naturally extended to integrate any combination of  omics data. We test this approach on breast cancer samples from the TCGA database, integrating data on messenger RNA, microRNAs and copy number variations. We show that the inclusion of the microRNA layer significantly improves the accuracy of subtype classification. Moreover, some of the hidden structures or "topics" that the algorithm extracts actually correspond to genes and microRNAs involved in breast cancer development, and are associated to the survival probability.

Simple SummaryTopic models are algorithms introduced for discovering hidden topics or latent variables in large unstructured text corpora. Leveraging on analogies between texts and gene expression profiles, these algorithms can be used to find structures in expression data. This work presents an application of topic modeling techniques to the identification of breast cancer subtypes. In particular, we extended a specific class of topic models to allow a multi-omics approach. As an illustrative example, considering both messenger RNA and microRNA expression levels, we were able to clearly distinguish healthy from tumor samples, as well as the different breast cancer subtypes. The integration of different layers of information is crucial for the observed classification accuracy. Our approach naturally provides the genes and the microRNAs associated to the specific topics that are used for the sample organization. We show that indeed these topics often contain genes involved in breast cancer development, and are associated to different survival probabilities.

category: cancer biology


title: Comparative assessment of electronic nicotine delivery systems aerosol and cigarette smoke on endothelial cell migration: the Replica Project

abstract:
Cigarette smoking is associated with impairment of repair mechanisms necessary for vascular endothelium homeostasis. Reducing the exposure to smoke toxicants may result in the mitigation of the harmful effect on the endothelium and cardiovascular disease development. Previous investigations performed by the tobacco industries evaluated in vitro the effect of electronic cigarette (e-cig) compared to cigarette smoke demonstrating a significant reduction in endothelial cell migration inhibition following e-cig aerosol exposure. In the present study, we replicated one of these studies, evaluating the effects of cigarette smoke on endothelial cell migration compared to e-cig and heated tobacco products. We used a multi-center approach (ring-study) to verify the robustness and reliability of the results obtained in the replicated study. Consistently with the original study, we observed a substantial reduction of the effects of e-cig and tobacco heated products on endothelial cell migration compared to cigarette smoke. In conclusion, our study further confirms the importance of e-cig and tobacco heated products as a possible harm reduction strategy for cardiovascular diseases development in smokers.

HighlightsO_LICigarette smoking is strictly related to impairment of vascular repair mechanisms
C_LIO_LIReducing the exposure to toxicants in smoke, could reduce the harm to endothelium
C_LIO_LIENDS showed a reduced effect on endothelial cell migration compared to cigarette
C_LIO_LIThese data demonstrated the reduced toxicity of ENDS compared to cigarettes.
C_LI

category: pharmacology and toxicology


title: Lack of pulmonary fibrogenicity and carcinogenicity of titanium dioxide nanoparticles in 26-week inhalation study in rasH2 mouse model

abstract:
BackgroundWith the rapid development of alternative methods based on the spirit of animal welfare, the publications of animal studies evaluating endpoints such as cancer have been extremely reduced. There have been no systemic inhalation exposure studies of titanium dioxide nanoparticles (TiO2 NPs) using CByB6F1-Tg(HRAS)2Jic (rasH2) 26-week study mice model for detecting carcinogenicity.

MethodsMale and female rasH2 mice were exposed to 2, 8 or 32 mg/m3 of TiO2 NPs for 6 hours/day, 5 days/week for 26 weeks using a whole-body inhalation exposure system, with reference to the Organization for Economic Co-operation and Development principles of Good Laboratory Practice. All tissues including lungs, and blood were collected and subjected to biological and histopathological analyses. Additionally, Ki67 positive index were evaluated in mice lung alveolar epithelial type 2 cell (AEC2).

ResultsThis study established a stable method for generating and exposing TiO2 NPs aerosol, and clarified the dose-response relationship by TiO2 NPs inhalation to rasH2 mice. TiO2 NPs exposure induced deposition of particles in lungs and mediastinal lymph nodes in a dose-dependent manner in each exposure group. Additionally, alveolar inflammation was only observed in 32 mg/m3 exposure group in both the sexes. Exposure to TiO2 NPs, as well as other organs, did not increase the incidence of lung tumors in any group, and pulmonary fibrosis and pre-neoplastic lesions were not observed in all groups. Finally, the cell proliferative activity of AEC2 was examined, and it was not increased by exposure to TiO2 NPs.

ConclusionsThis is the first report showing the lack of pulmonary fibrogenicity and carcinogenicity (no evidence of carcinogenic activity) of TiO2 NPs in 26-week inhalation study in rasH2 mice exposed up to 32 mg/m3, which is considered to be a high concentration. Macrophages undergoing phagocytosis due to TiO2 NPs exposure formed inflammatory foci in the alveolar regions of exposed mice but did not develop fibrosis or hyperplasia or tumors. Moreover, the cell proliferative ability of AEC2 in lesions was not increased. In addition, no carcinogenicity was observed for any organs other than the lungs in this study.


title: ROOT PENETRATION INDEX 3, a major quantitative trait locus (QTL) associated with root system penetrability in Arabidopsis.

abstract:
Soil mechanical impedance precludes root penetration, confining root system development to shallow soil horizons where mobile nutrients are scarce. Using a two-phase-agar system, we characterized Arabidopsis thaliana responses to low and high mechanical impedance at three root penetration stages. We found that seedlings whose roots fail to penetrate agar barriers show drastic changes in shoot and root morphology, while those capable of penetrating have only minor morphological effects. The assessment of 21 Arabidopsis accessions revealed that primary root penetrability (PRP) varies widely among accessions. To search for quantitative trait loci (QTLs) associated to root system penetrability, we evaluated a recombinant inbred population (RIL) derived from Landsberg erecta (Ler-0, with a high PRP) and Shahdara (Sha, with a low PRP) accessions. QTL analysis revealed a major-effect QTL localized in chromosome 3 (q-RPI3), which accounted for 29.98% (LOD = 8.82) of the total phenotypic variation. Employing an introgression line (IL-321), with a homozygous q-RPI3 region from Sha in the Ler-0 genetic background, we demonstrated that q-RPI3 plays a crucial role in root penetrability. This multiscale study revels new insights into root plasticity during the penetration process in hard agar layers, natural variation and genetic architecture behind primary root penetrability in Arabidopsis.

HighlightWe found a wide natural variation in the capacity of Arabidopsis accessions to penetrate hard agar layers. Using a Ler-0 x Sha recombinant inbred population, a major-effect QTL (q-RPI3) strongly associated with root penetrability of compact agar layers was identified.

category: plant biology


title: A synthetic membrane shaper for controlled liposome deformation

abstract: Shape defines the structure and function of cellular membranes. In cell division, the cell membrane deforms into a  dumbbell shape, while organelles such as the autophagosome exhibit  stomatocyte shapes. Bottom-up in vitro reconstitution of protein machineries that stabilize or resolve the membrane necks in such deformed liposome structures is of considerable interest to characterize their function. Here we develop a DNA-nanotechnology-based approach that we call Synthetic Membrane Shaper (SMS), where cholesterol-linked DNA structures attach to the liposome membrane to reproducibly generate high yields of stomatocytes and dumbbells. In silico simulations confirm the shape-stabilizing role of the SMS. We show that the SMS is fully compatible with protein reconstitution by assembling bacterial divisome proteins (DynaminA, FtsZ:ZipA) at the catenoidal neck of these membrane structures. The SMS approach provides a general tool for studying protein binding to complex membrane geometries that will greatly benefit synthetic cell research.

category: synthetic biology


title: Sustained TNF signaling is required for the synaptic and behavioral response to acute stress

abstract: Acute stress triggers plasticity of forebrain synapses as well as behavioral changes. Here we reveal that Tumor Necrosis Factor  (TNF) is a required downstream mediator of the stress response in mice, necessary for stress-induced synaptic potentiation in the ventral hippocampus and for an increase in anxiety-like behaviour. Acute stress is sufficient to activate microglia, triggering the long-term release TNF. Critically, on-going TNF signaling in the ventral hippocampus is necessary to sustain both the stress-induced synaptic and behavioral changes, as these could be reversed hours after induction by antagonizing TNF signaling. This demonstrates that TNF maintains the synaptic and behavioral stress response in vivo, making TNF a potential novel therapeutic target for stress disorders.

category: neuroscience


title: Correlation of receptor density and mRNA expression patterns in the human cerebral cortex

abstract:
Changes in distribution of associated molecular targets have been reported across several neuropsychiatric disorders. However, the high-resolution topology of most proteins is unknown and simultaneous in vivo measurement in multi-receptor systems is complicated. To account for the missing proteomic information, mRNA transcripts are typically used as a surrogate. Nonetheless, post-transcriptional and post-translational processes might cause the discrepancy between the final distribution of proteins and gene expression patterns. Therefore, this study aims to investigate ex vivo links between mRNA expression and corresponding receptor density in the human cerebral cortex.

To this end, autoradiography data on the density of 15 different receptors in 38 brain regions were correlated with the expression patterns of 50 associated genes derived from microarray data (mA), RNA sequencing data (RNA-Seq) provided by the Allen Human Brain Atlas and predicted mRNA expression patterns (pred-mRNA). Spearmans rank correlation was used to evaluate the possible links between proteomic data and mRNA expression patterns.

Correlations between mRNA and protein density varied greatly between targets: Positive associations were found for e.g. the serotonin 1A (pred-mRNA: rs = 0.708; mA: rs = 0.601) or kainate receptor (pred-mRNA: rs = 0.655; mA: rs = 0.601; RNA-Seq: rs = 0.575), while most of the investigated target receptors showed low or negative correlations.

The high variability in the correspondence of mRNA expression and receptor warrants caution when inferring the topology of molecular targets in the brain from transcriptome data. This highlights the longstanding value of molecular imaging data and the need for comprehensive proteomic data.


title: Remote contextual fear retrieval engages activity from salience network regions in rats

abstract: The ability to retrieve contextual fear memories depends on the coordinated activation of a brain-wide circuitry. Transition from recent to remote memories seems to involve the reorganization of this circuitry, a process called systems consolidation that has been associated with time-dependent fear generalization. However, it is unknown whether emotional memories acquired under different stress levels can undergo different systems consolidation processes. Here, we explored the activation pattern and functional connectivity of key brain regions associated with contextual fear conditioning (CFC) retrieval after recent (2 days) or remote (28 days) memory tests performed in rats submitted to strong (1.0mA footshock) or mild (0.3mA footshock) training. We used brain tissue from Wistar rats from a previous study, where we observed that increasing training intensity promotes fear memory generalization over time, possibly due to an increase in corticosterone (CORT) levels during memory consolidation. Analysis of Fos expression across 8 regions of interest (ROIs) allowed us to identify coactivation between them at both timepoints following memory recall. Our results showed that strong CFC elicits higher Fos activation in the anterior insular and prelimbic cortices during remote retrieval, which was positively correlated with freezing along with the basolateral amygdala. Rats trained either with mild or strong CFC showed broad functional connectivity at the recent timepoint whereas only animals submitted to the strong CFC showed a widespread loss of coactivation during remote retrieval. Post-training plasma CORT levels are positively correlated with FOS expression during recent retrieval in strong CFC, but negatively correlated with FOS expression during remote retrieval in mild CFC. Our findings suggest that increasing training intensity results in differential processes of systems consolidation, possibly associated with increased post-training CORT release, and that strong CFC engages activity from the aIC, BLA and PrL - areas associated with the Salience Network in rats - during remote retrieval.


title: Multimodal lesion network mapping to predict sensorimotor behavior in stroke patients

abstract: Lesion network mapping (LNM) has proved to be a successful technique to map symptoms to brain networks after acquired brain injury. Beyond the characteristics of a lesion, such as its etiology, size or location, LNM has shown that common symptoms in patients after injury may reflect the effects of their lesions on the same circuits, thereby linking symptoms to specific brain networks. Here, we extend LNM to its multimodal form, using a combination of functional and structural connectivity maps drawn from data from 1000 healthy participants in the Human Connectome Project. We applied the multimodal LNM to a cohort of 54 stroke patients with the aim of predicting sensorimotor behavior, as assessed through a combination of motor and sensory tests. Test scores were predicted using a Canonical Correlation Analysis with multimodal brain maps as independent variables, and cross-validation strategies were employed to overcome overfitting. The results obtained led us to draw three conclusions. First, the multimodal analysis reveals how functional connectivity maps contribute more than structural connectivity maps in the optimal prediction of sensorimotor behavior. Second, the maximal association solution between the behavioral outcome and multimodal lesion connectivity maps suggests an equal contribution of sensory and motor coefficients, in contrast to the unimodal analyses where the sensory contribution dominates in both structural and functional maps. Finally, when looking at each modality individually, the performance of the structural connectivity maps strongly depends on whether sensorimotor performance was corrected for lesion size, thereby eliminating the effect of larger lesions that produce more severe sensorimotor dysfunction. By contrast, the maps of functional connectivity performed similarly irrespective of any correction for lesion size. Overall, these results support the extension of LNM to its multimodal form, highlighting the synergistic and additive nature of different types of imaging modalities, and the influence of their corresponding brain networks on behavioral performance after acquired brain injury.


title: Post-error slowing reflects the joint impact of adaptive and maladaptive processes during decision making

abstract: Errors and their consequences are typically studied by investigating changes in decision speed and accuracy in trials that follow an error, commonly referred to as "post-error adjustments". Many studies have reported that subjects slow down following an error, a phenomenon called "post-error slowing" (PES). However, the functional significance of PES is still a matter of debate as it is not always adaptive. That is, it is not always associated with a gain in performance and can even occur with a decline in accuracy. Here, we hypothesized that the nature of PES is influenced by ones speed-accuracy tradeoff policy, which determines the overall level of choice accuracy in the task at hand. To test this hypothesis, we investigated post-error adjustments in subjects performing the same task while they were required to either emphasize speed (low accuracy) or cautiousness (high accuracy) in two distinct contexts (hasty and cautious contexts, respectively) experienced on separate days. Accordingly, our data indicate that post-error adjustments varied according to the context in which subjects performed the task, with PES being solely significant in the hasty context. In addition, we only observed a gain in performance after errors in a specific trial type, suggesting that post-error adjustments depend on a complex combination of processes that affect the speed of ensuing actions as well as the degree to which such PES comes with a gain in performance.


title: Dual Passive Reactive Brain Computer Interface: a Novel Approach to Human-Machine Symbiosis

abstract: The present study proposes a novel concept of neuroadaptive technology, namely a dual passive-reactive Brain-Computer Interface (BCI), that enables bi-directional interaction between humans and machines. We have implemented such a system in a realistic flight simulator using the NextMind classification algorithms and framework to decode pilots intention (reactive BCI) and to infer their level of attention (passive BCI). Twelve pilots used the reactive BCI to perform checklists along with an anti-collision radar monitoring task that was supervised by the passive BCI. The latter simulated an automatic avoidance maneuver when it detected that pilots missed an incoming collision. The reactive BCI reached 100% classification accuracy with a mean reaction time of 1.6s when exclusively performing the checklist task. Accuracy was up to 98.5% with a mean reaction time of 2.5s when pilots also had to fly the aircraft and monitor the anti-collision radar. The passive BCI achieved a F1 - score of 0.94. This first demonstration shows the potential of a dual BCI to improve human-machine teaming which could be applied to a variety of applications.


title: Theta oscillations in the prefrontal-hippocampal circuit do not couple to respiration-related oscillations

abstract: Oscillatory activity is thought to coordinate neural computations across brain regions, and theta oscillations are critical for learning and memory. Because the frequency of respiratory-related oscillations (RROs) in rodents can overlap with the frequency of theta in the prefrontal cortex (PFC) and the hippocampus, we asked whether odor-cued working memory may be supported by coupling between these two oscillations. We first confirmed that RROs are propagated to the hippocampus and PFC and that RRO frequency overlaps with canonical theta frequency. However, we found low coherence between RROs and local theta oscillations in the hippocampus-PFC network when the two types of oscillations overlapped in frequency. This effect was observed during all behavioral phases including during movement and while odors were actively sampled when stationary. Despite the similarity in frequency, RROs and theta oscillations therefore appear to be limited to supporting computation in distinct networks, which suggests that sustained long-range coordination between oscillation patterns that depend on separate pacemakers is not necessary to support at least one type of working memory.


title: Spherical harmonic based noise rejection and neuronal sampling with multi-axis OPMs

abstract: In this study we explore the interference rejection and spatial sampling properties of multi-axis Optically Pumped Magnetometer (OPM) data. We use both vector spherical harmonics and eigenspectra to quantify how well an array can separate neuronal signal from environmental interference while adequately sampling the entire cortex. We found that triaxial OPMs have superb noise rejection properties allowing for very high orders of interference (L=6) to be accounted for while minimally affecting the neural space (2dB attenuation for a 60-sensor triaxial system). To adequately model the signals arising from the cortex, we show that at least 11th order (143 spatial degrees of freedom) irregular solid harmonics or 95 eigenvectors of the lead field are needed to model the neural space for OPM data (regardless of number of axes measured). This can be adequately sampled with 75-100 equidistant triaxial sensors (225-300 channels) or 200 equidistant radial channels. In other words, ordering the same number of channels in triaxial (rather than purely radial) configuration gives significant advantages not only in terms of external noise rejection but also minimizes cost, weight and cross-talk.


title: Structural basis for antibody resistance to SARS-CoV-2 omicron variant

abstract:
The recently reported B.1.1.529 Omicron variant of SARS-CoV-2 includes 34 mutations in the spike protein relative to the Wuhan strain that initiated the COVID-19 pandemic, including 15 mutations in the receptor binding domain (RBD). Functional studies have shown omicron to substantially escape the activity of many SARS-CoV-2-neutralizing antibodies. Here we report a 3.1 [A] resolution cryo-electron microscopy (cryo-EM) structure of the Omicron spike protein ectodomain. The structure depicts a spike that is exclusively in the 1-RBD-up conformation with increased mobility and inter-protomer asymmetry. Many mutations cause steric clashes and/or altered interactions at antibody binding surfaces, whereas others mediate changes of the spike structure in local regions to interfere with antibody recognition. Overall, the structure of the omicron spike reveals how mutations alter its conformation and explains its extraordinary ability to evade neutralizing antibodies.

HighlightsO_LISARS-CoV-2 omicron spike exclusively adopts 1-RBD-up conformation
C_LIO_LIOmicron substitutions alter conformation and mobility of RBD
C_LIO_LIA subset of omicron mutations change the local conformation of spike
C_LIO_LIThe structure reveals the basis of antibody neutralization escape
C_LI

category: immunology


title: Anti-SARS-CoV-2 human antibodies retaining neutralizing activity against SARS-CoV-2 B.1.1.529 (omicron)

abstract: SARS-CoV-2 B.1.1.529, designated omicron, was recently identified as a new variant of concern by WHO and is rapidly replacing SARS-CoV-2 delta as the most dominant variant in many countries. Unfortunately, because of the high number of mutations present in the spike of SARS-CoV-2 omicron, most monoclonal antibodies (mAbs) currently approved for treatment of COVID-19 lose their in vitro neutralizing activity against this variant. We recently described a panel of human anti-SARS-CoV-2 mAbs that potently neutralize SARS-CoV-2 Wuhan, D614G and variants alpha, beta, gamma and delta. In this work, we evaluated our mAb panel for potential in vitro activity against SARS-CoV-2 delta and omicron. Three mAbs from our panel retain neutralizing activity against both delta and omicron, with mAb 3B8 still resulting in complete neutralization at a concentration as low as 0.02 g/ml for both variants. Overall, our data indicate that mAb 3B8 may have the potential to become a game-changer in the fight against the continuously evolving SARS-CoV-2.


title: High throughput hemogram of T cells using digital holographic microscopy and deep learning

abstract: T cells of the adaptive immune system provide effective protection to the human body against numerous pathogenic challenges. Current labelling methods of detecting these cells, such as flow cytometry or magnetic bead labelling, are time consuming and expensive. To overcome these limitations, the label-free method of digital holographic microscopy (DHM) combined with deep learning has recently been introduced which is both time and cost effective. In this study, we demonstrate the application of digital holographic microscopy with deep learning to classify the key CD4+ and CD8+ T cell subsets. We show that combining DHM of varying fields of view, with deep learning, can potentially achieve a classification throughput rate of 78,000 cells per second with an accuracy of 76.2% for these morphologically similar cells. This throughput rate is 100 times faster than the previous studies and proves to be an effective replacement for labelling methods.


title: MAVS mediates a protective immune response in the brain to Rift Valley fever virus

abstract:
Rift Valley fever virus (RVFV) is a highly pathogenic mosquito-borne virus capable of causing hepatitis, encephalitis, blindness, hemorrhagic syndrome, and death in humans and livestock. Upon aerosol infection with RVFV, the brain is a major site of viral replication and tissue damage, yet pathogenesis in this organ has been understudied. Here, we investigated the immune response in the brain of RVFV infected mice. In response to infection, microglia initiate robust transcriptional upregulation of antiviral immune genes, as well as increased levels of activation markers and cytokine secretion that is dependent on mitochondrial antiviral-signaling protein (MAVS) and independent of toll-like receptors 3 and 7. In vivo, Mavs-/- mice displayed enhanced susceptibility to RVFV as determined by increased brain viral burden and higher mortality. Single-cell RNA sequence analysis identified microglia-specific defects in type I interferon and interferon responsive gene expression in Mavs-/- mice, as well as dysregulated lymphocyte infiltration. The results of this study provide a crucial step towards understanding the precise molecular mechanisms by which RVFV infection is controlled in the brain and will help inform the development of vaccines and antiviral therapies that are effective in preventing encephalitis.

Author SummaryRift Valley fever virus causes severe disease in humans and livestock and in some cases can be fatal. There is concern about the use of Rift Valley fever virus as a bioweapon since it can be transmitted through the air, and there are no vaccines or antiviral treatments. Airborne transmission of the virus causes severe inflammation of the brain, yet little is known about the immune response against the virus in this organ. Here, we investigated the immune response in the brain to Rift Valley fever virus following intranasal infection. We determined that microglia, the resident immune cells of the brain, initiate a robust response to Rift Valley fever virus infection and identified a key immune pathway that is critical for the ability of microglia to respond to infection. When this immune pathway is rendered non-functional, mice have a dysregulated response to infection in the brain. This study provides insight into how the immune response can control Rift Valley fever virus infection of the brain.


title: Construction of a potent pan-vaccine based on the evolutionary tendency of SARS-CoV-2 spike protein.

abstract: SARS-CoV-2 continued to spread globally along with different variants. Here, we systemically analyzed viral infectivity and immune-resistance of SARS-CoV-2 variants to explore the underlying rationale of viral mutagenesis. We found that the Beta variant harbors both high infectivity and strong immune resistance, while the Delta variant is the most infectious with only a mild immune-escape ability. Remarkably, the Omicron variant is even more immune-resistant than the Beta variant, but its infectivity increases only in Vero E6 cells implying a probable preference for the endocytic pathway. A comprehensive analysis revealed that SARS-CoV-2 spike protein evolved into distinct evolutionary paths of either high infectivity plus low immune resistance or low infectivity plus high immune resistance, resulting in a narrow spectrum of the current single-strain vaccine. In light of these findings and the phylogenetic analysis of 2674 SARS-CoV-2 S-protein sequences, we generated a consensus antigen (S6) taking the most frequent mutations as a pan-vaccine against heterogeneous variants. As compared to the ancestry SWT vaccine with significantly declined neutralizations to emerging variants, the S6 vaccine elicits broadly neutralizing antibodies and full protections to a wide range of variants. Our work highlights the importance and feasibility of a universal vaccine strategy to fight against antigen drift of SARS-CoV-2.

category: microbiology


title: Nucleocapsid 203 mutations enhance SARS-CoV-2 immune evasion

abstract: Previous work indicated that the nucleocapsid 203 mutation increase the virulence and transmission of the SARS-CoV-2 Alpha variant. However, Delta later outcompeted Alpha and other lineages, promoting a new wave of infections. Delta also possesses a nucleocapsid 203 mutation, R203M. Large-scale epidemiological analyses suggest a synergistic effect of the 203 mutation and the spike L452R mutation, associated with Delta expansion. Viral competition experiments demonstrate the synergistic effect in fitness and infectivity. More importantly, we found that the combination of R203M and L452R brings in a 3.2-fold decrease in neutralizing titers to the neutralizing serum relative to L452R-only virus. R203M/L452R show an increased fitness after the initiation of global vaccination programmes, possibly associated with the enhanced immune evasion. Another rapidly emerging variant Omicron also bears the 203 mutation. Thus, we proposed that nucleocapsid mutations play an essential role for the rise and predominance of variants in concern.


title: Direct RNA sequencing of Respiratory Syncytial Virus infected human cells generates a detailed overview of RSV polycistronic mRNA and transcript abundance

abstract: To characterize species of viral mRNA transcripts generated during respiratory syncytial virus (RSV) infection, human fibroblast-like MRC5 lung cells were infected with subgroup A RSV for 6, 16 and 24 hours. Total RNA was harvested and polyadenylated mRNA was enriched and sequenced by direct RNA sequencing on an Oxford nanopore device. This yielded over 150,000 direct mRNA transcript reads which were mapped to the viral genome and analysed to determine relative mRNA levels of viral genes using our in-house ORF-centric pipeline. We were also able to examine frequencies with which polycistronic readthrough mRNAs were generated and to assess the length of the polyadenylated tails for each group of transcripts. We show that there is a general but non-linear decline in gene transcript abundance across the viral genome, as predicted by the model of RSV gene transcription. However, the decline in transcript abundance is not consistent. We show that the polyadenylate tails generated by the viral polymerase are similar in length to those generated by the host cells polyadenylation machinery and broadly declined in length for most transcripts as infection progressed. Finally, we observed that the steady state abundance of transcripts with very short polyadenylate tails is much less for N, SH and G transcripts compared to NS1, NS2, P, M, F and M2 which may reflect differences in mRNA stability and/or translation rates.


title: Effects of α-crystallin gene knockout on zebrafish lens development

abstract: The -crystallin small heat shock proteins contribute to the transparency and refractive properties of the vertebrate eye lens and prevent the protein aggregation that would otherwise produce lens cataract, the leading cause of human blindness. There are conflicting data in the literature as to what role the -crystallins may play in early lens development. In this study we used CRISPR gene editing to produce zebrafish lines with null mutations for each of the three -crystallin genes (cryaa, cryaba and cryabb). Absence of protein was confirmed by mass spectrometry and lens phenotypes were assessed with differential interference contrast microscopy and histology. Loss of A-crystallin produced a variety of lens defects with varying severity in larval lenses at 3 and 4 dpf, but little significant change in normal fiber cell denucleation. Loss of either Ba- or Bb-crystallin produced no significant lens defects. Mutation of each -crystallin gene did not alter the expression levels of the remaining two, suggesting a lack of genetic compensation. These data confirm a developmental role for A-crystallin in lens development, but the range of phenotype severity suggests its loss simply increases the chance for defect, and that the protein is not essential. Our finding that cryaba and cryabb null mutants lack noticeable lens defects is congruent with insignificant transcript levels in lens epithelial and fiber cells. Future experiments can explore the molecular consequences of cryaa mutation and causes of lens defects in this null mutant, as well as the roles of other genes in lens development and function.

category: developmental biology


title: A single-cell atlas of human fetal lung development between 14 and 19 weeks of gestation

abstract:
RationaleHuman lung development has been mainly described in morphologic studies and the potential underlying molecular mechanisms were extrapolated from animal models. Therefore, there is a need to gather knowledge from native human lung tissue. In this study we describe changes at a single-cell level in human fetal lungs during the pseudoglandular stage.

MethodsWe report the cellular composition, cell trajectories and cell-to-cell communication in developing human lungs with single-nuclei RNA sequencing (snRNA-seq) on 23,251 nuclei isolated from nine human fetuses with gestational ages between 14 to 19 weeks of gestation.

ResultsWe identified nine different cell types, including a rare pulmonary neuroendocrine cells population. For each cell type, marker genes are reported, and selected marker genes are used for spatial validation with fluorescent RNA in situ hybridization. Enrichment and developmental trajectory analysis provide insight into molecular mechanisms and signaling pathways within individual cell clusters according to gestational age. Lastly, ligand-receptor analysis highlights determinants of cell-to-cell communication among the different cell types through the pseudoglandular stage, including general developmental pathways (NOTCH and TGFB), as well as more specific pathways involved in vasculogenesis, neurogenesis, and immune system regulation.

ConclusionThese findings provide a clinically relevant background for research hypotheses generation in projects studying normal or impaired lung development and help to develop and validate surrogate models to study human lung development, such as human lung organoids.

TAKE HOME MESSAGEUsing a single-cell transcriptomic approach (single-nuclei RNA sequencing), we describe here, for the first time, the cellular landscape, cell developmental trajectories, and cell-to-cell communication in the developing human lung during the pseudoglandular stage.


title: Restricted proliferation during neurogenesis contributes to regionalization of the amphioxus nervous system.

abstract: The central nervous system of the cephalochordate amphioxus consists of a dorsal neural tube with an anterior brain. Two decades of gene expression analyses in developing amphioxus embryos have shown that despite the lack of overt segmentation the amphioxus neural tube is highly regionalized at the molecular level. However, little is known about the mechanisms that generate such precise regionalization. Proliferation is a key driver of pattern formation and cell type diversification, but in amphioxus it has never been studied in detail nor in the specific context of neurogenesis. Here, we describe the dynamics of cell division during the formation of the central nervous system in amphioxus embryos and its contributions to the regionalization of the neural axis. We show that after gastrulation, proliferation pauses to become spatially restricted to the anterior and posterior ends of the neural tube at neurula stages. Only at the onset of larval life, proliferation resumes in the central part of the nervous system. By marking specific populations and inhibiting cell division during neurulation, we demonstrate that proliferation in the anterior cerebral vesicle is required to establish the full cell type repertoire of the frontal eye complex and the putative hypothalamic region of the amphioxus brain, while posterior proliferating progenitors, which were found here to derive from the dorsal lip of the blastopore, contribute to elongate the caudal floor plate. Between these proliferative domains, we find trunk nervous system differentiation is independent from cell division, which decreases during neurulation and resumes at the early larval stage. Taken together, our results highlight multiple roles for proliferation in shaping the amphioxus nervous system.


title: Applying a Modified Metabarcoding Approach for the Sequencing of Macrofungal Specimens from Fungarium Collections

abstract:
O_LIPremise: Fungaria are a largely untapped source for understanding fungal biodiversity. The effort and cost in producing DNA barcode sequence data for large numbers of fungal specimens can be prohibitive. This study applies a modified metabarcoding approach that provides a labor and cost-effective solution for sequencing the fungal DNA barcode from hundreds of specimens at once.
C_LIO_LIMethods: A two-step PCR approach uses nested barcoded primers to nrITS2 sequence data. We applied this to 766 macrofungal specimens that represent a broad taxonomic sampling of the Dikarya, of which 382 Lactarius specimens are used to identify molecular operational taxonomic units (MOTUs) through a phylogenetic approach. Scripts in Python and R were used to organize sequence data and execute packages CutAdapt and DADA2 were used for primer removal and assessing sequence quality. Sequences were compared to NCBI and UNITE databases and Sanger-produced sequences.
C_LIO_LIResults: Specimen taxonomic identities from nrITS2 sequence data are >90% accurate across all specimens sampled. Phylogenetic analysis of Lactarius sequences identified 20 MOTUs.
C_LIO_LIDiscussion: The results demonstrate the capacity of these methods to produce nrITS2 sequences from large numbers of fungarium specimens. This provides an opportunity to more effectively use fungarium collections in advancing fungal diversity identification and documentation.
C_LI


title: Weakly deleterious natural genetic variation amplifies probability of resistance in multiplexed gene drive systems

abstract: Evolution of resistance is a major barrier to successful deployment of gene drive systems to suppress natural populations, which could greatly reduce the burden of many vector borne diseases. Multiplexed guide RNAs that require resistance mutations in all target cut sites is a promising anti-resistance strategy, since in principle resistance would only arise in unrealistically large populations. Using novel stochastic simulations that accurately model evolution at very large population sizes, we explore the probability of resistance due to three important mechanisms: 1) non-homologous end-joining mutations, 2) single nucleotide mutants arising de novo or, 3) single nucleotide polymorphisms pre-existing as standing variation. Our results explore the relative importance of these mechanisms and highlight a complexity of the mutation-selection-drift balance between haplotypes with complete resistance and those with an incomplete number of resistant alleles. We find this leads to a qualitatively new phenomenon where weakly deleterious naturally occurring variants greatly amplify the probability of multi-site resistance. This challenges the intuition that many target sites would guarantee prevention of resistance, where in the face of standing genetic variation, it can be probable even in not very large populations. This result has broad application to resistance arising in many multi-site evolutionary scenarios including multi-drug resistance to antibiotics, antivirals and cancer treatments, as well as the evolution of vaccine escape mutations in large populations.


title: Weakly deleterious natural genetic variation greatly amplifies probability of resistance in multiplexed gene drive systems


title: Adjustment for SNPs instead of demographic and clinical covariates in GWAS

abstract: With capabilities of sequencing ancient DNA to high coverage often limited by sample quality or cost, imputation of missing genotypes presents a possibility to increase power of inference as well as cost-effectiveness for the analysis of ancient data. However, the high degree of uncertainty often associated with ancient DNA poses several methodological challenges, and performance of imputation methods in this context has not been fully explored. To gain further insights, we performed a systematic evaluation of imputation of ancient data using Beagle 4.0 and reference data from phase 3 of the 1000 Genomes project, investigating the effects of coverage, phased reference and study sample size. Making use of five ancient samples with high-coverage data available, we evaluated imputed data with respect to accuracy, reference bias and genetic affinities as captured by PCA. We obtained genotype concordance levels of over 99% for data with 1x coverage, and similar levels of accuracy and reference bias at levels as low as 0.75x. Our findings suggest that using imputed data can be a realistic option for various population genetic analyses even for data in coverage ranges below 1x. We also show that a large and varied phased reference set as well as the inclusion of low-to moderate-coverage ancient samples can increase imputation performance, particularly for rare alleles. In-depth analysis of imputed data with respect to genetic variants and allele frequencies gave further insight into the nature of errors arising during imputation, and can provide practical guidelines for post-processing and validation prior to downstream analysis.


title: Electron relay via diffusible and protein-bound flavin cofactor in human six-transmembrane epithelial antigen of the prostate (STEAP)

abstract: Six transmembrane epithelial antigen of the prostate (STEAP) is a family of membrane-embedded metal ion reductases that transfer electrons across the cell membranes. STEAPs are unique to mammals and implicated in metabolic and inflammatory responses and are significantly upregulated in many types of cancer cells. There are four members in the family, STEAP1 - 4, and all STEAPs have a transmembrane domain (TMD) that has a conserved heme binding site, and STEAP2 - 4, but not STEAP1, have an intracellular reductase domain (RED) that binds to NADPH and FAD. NADPH, FAD, and heme form an electron transfer chain that allows electron flow across the cell membranes, however, the mechanism of the stepwise crossmembrane electron transfer remains unclear. It is also unclear how STEAP1, which does not have a RED, acquires and transfers electrons. We expressed and purified human STEAP2 (hSTEAP2), and constructed the electron transfer chain in vitro. Purified hSTEAP2 mediates electron transfer from NADPH to FAD and to heme, with a NADPH oxidation rate of 2.6 x 10-3 s-1. The time course for reduction of heme is more complex with an initial rate of ~ 1.6 x 10-4 s-1. We also found that the heme in hSTEAP2 has a low-spin electron structure and thus a rigid coordination, which is consistent with its high occupancy in the purified protein and its role as part of the electron transfer chain. We then determined the structure of hSTEAP2 in complex with NADP+, FAD, and heme by cryo-electron microscopy to 3.2 [A]. Human STEAP2 forms a homotrimer and its structure is similar to that of hSTEAP4. NADP+, FAD, and heme are well-resolved in the structure, and while the current conformation would allow electron transfer from FAD to heme, the FAD isoalloxazine ring is ~ 19 [A] away from NADPH and does not support hydride transfer. Significant structural changes are required to accommodate dissociation of the FAD isoalloxazine ring from the TMD such that the FAD may become diffusible after its reduction. To test this hypothesis and also to find out how STEAP1 may transfer electrons, we reconstructed an electron transfer chain for STEAP1 and found that the heme in STEAP1 can be reduced by FAD produced either by the full-length STEAP2 or by the soluble RED domain from STEAP4. These results support a "diffusible FAD" mechanism and demonstrate that STEAP1 is capable of mediating electron transfer across the cell membranes. In summary, our study established a structural and functional framework for further analyses for resolving the mechanism of electron transfer in STEAPs.


title: Growth kinetics and power laws indicate distinct mechanisms of cell-cell interactions in the aggregation process

abstract: Cellular aggregation is a complex process orchestrated by various kinds of interactions depending on its environments. Different interactions give rise to different pathways of cellular rearrangement and the development of specialized tissues. To distinguish the underlying mechanisms, in this theoretical work, we investigate the spontaneous emergence of tissue patterns from an ensemble of single cells on a substrate following three leading pathways of cell-cell interactions, namely, direct cell adhesion contacts, matrix mediated mechanical interaction, and chemical signalling. Our analysis shows that the growth kinetics of the aggregation process is distinctly different for each pathway and bears the signature of the specific cell-cell interactions. Interestingly, we find that the average domain size and the mass of the clusters exhibit a power law growth in time under certain interaction mechanisms hitherto unexplored. Further, as observed in experiments, the cluster size distribution can be characterized by stretched exponential functions showing distinct cellular organization processes.


title: Transport receptor occupancy in Nuclear Pore Complex mimics

abstract: Nuclear Pore Complexes (NPCs) regulate all molecular transport between the nucleus and the cytoplasm in eukaryotic cells. Intrinsically disordered Phe-Gly nucleoporins (FG Nups) line the central conduit of NPCs to impart a selective barrier where large proteins are excluded unless bound to a transport receptor (karyopherin; Kap). Here, we assess  Kap-centric NPC models, which postulate that Kaps participate in establishing the selective barrier. We combine biomimetic nanopores, formed by tethering Nsp1 to the inner wall of a solid-state nanopore, with coarse-grained modeling to show that yeast Kap95 exhibits two populations in Nsp1-coated pores: one population that is transported across the pore in milliseconds, and a second population that is stably assembled within the FG mesh of the pore. Ionic current measurements show a conductance decrease for increasing Kap concentrations and noise data indicate an increase in rigidity of the FG-mesh. Modeling reveals an accumulation of Kap95 near the pore wall, yielding a conductance decrease. We find that Kaps only mildly affect the conformation of the Nsp1 mesh and that, even at high concentrations, Kaps only bind at most 8% of the FG-motifs in the nanopore, indicating that Kap95 occupancy is limited by steric constraints rather than by depletion of available FG-motifs. Our data provide an alternative explanation of the origin of bimodal NPC binding of Kaps, where a stable population of Kaps binds avidly to the NPC periphery, while fast transport proceeds via a central FG-rich channel through lower affinity interactions between Kaps and the cohesive domains of Nsp1.


title: Optimizing a High Performing Multiplex-CRISPRi P. putida strain with Integrated Metabolomics and 13C-Metabolic Flux Analyses

abstract: Microbial cell factory development often faces bottlenecks after initial rounds of design-build-test-learn (DBTL) cycles as engineered producers respond unpredictably to further genetic modifications. Thus, deciphering metabolic flux and correcting bottlenecks are key components of DBTL cycles. Here, a 14-gene edited Pseudomonas putida KT2440 strain for heterologous indigoidine production was examined using both 13C-metabolic flux analysis (13C-MFA) and metabolite measurements. The results indicated the conservation of the cyclic Entner-Doudoroff (ED)-EMP pathway flux, downregulation of the TCA cycle and pyruvate shunt, and glyoxylate shunt activation. At the metabolite level, the CRISPR/dCpf1-interference mediated multiplex repression decreased gluconate/2-ketogluconate secretion and altered several intracellular TCA metabolite concentrations, leading to succinate overflow. Further strain engineering based on the metabolic knowledge first employed an optimal ribosome binding site (RBS) to achieve stronger product-substrate growth coupling (1.6-fold increase). Then, deletion strains were constructed using ssDNA recombineering. Of the five strains tested, deletion of the PHA operon ({Delta}phaAZC-IID) resulted in a 2.2-fold increase in growth phase production compared to the optimal RBS construct. After 72 h of batch cultivation, the {Delta}phaAZC-IID strain had 1.5-fold and 1.8-fold increases of indigoidine titer compared to the improved RBS construct and the original strain, respectively. Overall, the findings provided actionable DBTL targets as well as insights into physiological responses and flux buffering when new recombineering tools were used for engineering P. putida KT2440.


title: Identifying control ensembles for decision-making within the cortico-basal ganglia-thalamic circuit

abstract:
In situations featuring uncertainty about action-reward contingencies, mammals exhibit a high degree of flexibility in adopting strategies for decision-making that are tuned in response to the conditions that they encounter. Although the cortico-basal ganglia thalamic (CBGT) network is implicated in information processing during decisions, it features a complex synaptic architecture, comprised of multiple feed-forward, reciprocal, and feedback pathways, that complicate efforts to elucidate the roles of specific CBGT populations in the process of evidence accumulation. In this paper we apply a strategic sampling approach, based on Latin hypercube sampling, to explore how variations in CBGT network properties, including subpopulation firing rates and synaptic weights, map to variability of parameters in a normative drift diffusion model (DDM), representing algorithmic aspects of information processing during decision-making. Through the application of canonical correlation analysis, we find that this relationship can be characterized in terms of three low-dimensional control ensembles within the CBGT network that impact specific qualities of the emergent decision policy: responsiveness (associated with overall activity in corticothalamic and direct pathways), pliancy (associated largely with overall activity in components of the indirect pathway of the basal ganglia), and choice (associated with differences in direct and indirect pathways across action channels). These analyses provide key mechanistic predictions about the roles of specific CBGT network elements in tuning information processing dynamics during decisions.

Author summaryMammals are continuously subjected to uncertain situations in which they have to choose among behavioral options. The cortico-basal ganglia-thalamic (CBGT) circuit is a complicated collection of interconnected nuclei believed to strongly influence the ability to adapt to environmental changes. The roles of specific CBGT components in controlling information during decisions remains unclear. At a more phenomenological, algorithmic level, drift-diffusion models have been shown to be able to reproduce behavioral data (action selection probabilities and the time needed to make a decision) obtained experimentally from mammals and to provide an abstract representation of a decision policy. In this work, we use simulated decision-making to establish a mapping from neural activity in the CBGT circuit to behavioral outcomes. This mapping illuminates the importance of three core sets of CBGT subnetworks in the action selection process and how they are involved in adapting decision policies across exploitative and exploratory situations.


title: Identifying control ensembles for information processing within the cortico-basal ganglia-thalamic circuit


title: Crosstalk between biochemical signaling and trafficking determines AMPAR dynamics in synaptic plasticity

abstract: Synaptic plasticity involves the modification of both biochemical and structural components of neurons. Many studies have revealed that the change in the number density of the glutamatergic receptor AMPAR at the synapse is proportional to synaptic weight update; increase in AMPAR corresponds to strengthening of synapses while decrease in AMPAR density weakens synaptic connections. The dynamics of AMPAR are thought to be regulated by upstream signaling, primarily the calcium-CaMKII pathway, trafficking to and from the synapse, and influx from extrasynaptic sources. Here, we have developed a set of models using compartmental ordinary differential equations to systematically investigate contributions of signaling and trafficking variations on AMPAR dynamics at the synaptic site. We find that the model properties including network architecture and parameters significantly affect the integration of fast upstream species by slower downstream species. Furthermore, we predict that the model outcome, as determined by bound AMPAR at the synaptic site, depends on (a) the choice of signaling model (bistable CaMKII or monostable CaMKII dynamics), (b) trafficking versus influx contributions, and (c) frequency of stimulus. Therefore, AMPAR dynamics can have unexpected dependencies when upstream signaling dynamics (such as CaMKII and PP1) are coupled with trafficking modalities.


title: Subject-specific features of excitation/inhibition profiles in neurodegenerative diseases

abstract: Brain pathologies are based on microscopic changes in neurons and synapses that reverberate into large scale networks altering brain dynamics and functional states. An important yet unresolved issue concerns the impact of patients excitation/inhibition profiles on neurodegenerative diseases including Alzheimers disease, Frontotemporal Dementia and Amyotrophic Lateral Sclerosis. In this work we used a simulation platform, The Virtual Brain, to simulate brain dynamics in healthy controls and in Alzheimers disease, Frontotemporal Dementia and Amyotrophic Lateral Sclerosis patients. The brain connectome and functional connectivity were extracted from 3T-MRI scans and The Virtual Brain nodes were represented by a Wong-Wang neural mass model endowing an explicit representation of the excitatory/inhibitory balance. The integration of cerebro-cerebellar loops improved the correlation between experimental and simulated functional connectivity, and hence The Virtual Brain predictive power, in all pathological conditions. The Virtual Brain biophysical parameters differed between clinical phenotypes, predicting higher global coupling and inhibition in Alzheimers disease and stronger NMDA (N-methyl-D-aspartate) receptor-dependent excitation in Amyotrophic Lateral Sclerosis. These physio-pathological parameters allowed an advanced analysis of patients state. In backward regressions, The Virtual Brain parameters significantly contributed to explain the variation of neuropsychological scores and, in discriminant analysis, the combination of The Virtual Brain parameters and neuropsychological scores significantly improved discriminative power between clinical conditions. Eventually, cluster analysis provided a unique description of the excitatory/inhibitory balance in individual patients. In aggregate, The Virtual Brain simulations reveal differences in the excitatory/inhibitory balance of individual patients that, combined with cognitive assessment, can promote the personalized diagnosis and therapy of neurodegenerative diseases.


title: Distinct hippocampal network states support theta phase precession and theta sequences in CA1

abstract: The hippocampus likely uses temporal coding to represent complex memories via mechanisms such as theta phase precession and theta sequences. Theta sequences are rapid sweeps of spikes from multiple place cells, encoding past or planned trajectories or non-spatial information. Phase precession, the correlation between a place cells theta firing phase and animal position has been suggested to facilitate sequence emergence. We find that CA1 phase precession varies strongly across cells and environmental contingencies. Phase precession depends on the CA1 network state, and is only present when the medium gamma oscillation (60-90 Hz, linked to Entorhinal inputs) dominates. Conversely, theta sequences are most evident for non-precessing cells or with leading slow gamma (20-45 Hz, linked to CA3 inputs). These results challenge the view that phase precession is the mechanism underlying the emergence of theta sequences and point at a "dual network states" model for hippocampal temporal code, potentially supporting merging of memory and exogenous information in CA1.


title: Toll-like receptor 5 as a novel receptor for fungal zymosan

abstract: Microbial pathogens carry specific structural patterns which were termed pathogen-associated molecular patterns (PAMPs). Toll-like receptors (TLRs) as key elements for the recognition of microbial pathogens are necessary for the activation of innate immune pathways. TLRs are activated by binding PAMPs of bacteria, viruses and fungi and initiate a signaling pathway resulting in the activation of transcription factors which modulate the production of various proinflammatory cytokines. It is not fully clear in detail which microbial pattern is recognized by which TLR. Here we show for the first time that TLR5 is a strong receptor for the yeast particle zymosan. We have generated stable human cell lines with combinations of TLR2 and TLR5 knock in/knock out together with stable nuclear factor kappaB (NF-{kappa}B) luciferase reporters. We found that both receptors TLR5 and TLR2 lead to an independent activation of the NF-{kappa}B pathway when simulated with zymosan. Our results demonstrate that TLR5 is a receptor for the fungal particle zymosan in addition to bacterial fragments like flagellin. Distinct cytokine patterns might suggest that TLR5 is potentially important for the differentiation in the recognition of the specific type of the foreign microorganisms and in the specific host defense response.


title: sRNAs enriched in outer membrane vesicles of pathogenic Flavobacterium psychrophilum interact with immune genes of rainbow trout

abstract:
Outer membrane vesicles (OMVs) released by gram-negative bacteria during host-pathogen interactions harbor cargos, such as DNA, RNA, toxins, and virulence factors. We hypothesized that sRNAs carried within OMVs of Flavobacterium psychrophilum interact with host immune genes and affect their expression. OMVs were isolated from F. psychrophilum and visualized using transmission electron microscopy (TEM). RNA-Seq datasets generated from whole-cell F. psychrophilum and their OMVs indicated enrichment of specific sRNAs in the OMVs compared to the parent cell. Fluorescent in situ hybridization (FISH) and confocal microscopy confirmed the expression of a randomly chosen sRNA.

Integrated RNA-Seq analyses of host transcriptome and bacterial sRNAs on day 5 post-infection of F. psychrophilum-resistant and -susceptible rainbow trout genetic lines revealed 516 protein-coding, 595 lncRNA, and 116 bacterial sRNA differentially expressed (DE) transcripts. Integrated and network analyses of these DE transcripts revealed immune genes targeted by bacterial sRNAs. On the top of these genes, an isoform encoding anaphase-promoting complex subunit 13 (ANAPC13_1) was highly upregulated and exhibited interaction and reciprocal expression with 21 DE sRNAs enriched in OMVs and/or located in pathogenicity islands (PAIs). In vitro treatment of the rainbow trout epithelial cell line RTgill-W1 with OMVs showed signs of cell autolysis accompanied by dynamic changes in expression of host genes when profiled 24h following treatment. The OMV-enriched sRNAs, soFE013584 and soFE002123, showed high interactions with the protection of telomeres 1 gene (POT1); essential for chromosome stability and cellular viability. Modulation of the host gene expression following OMV-treatment, which favors elements from the phagocytic, endocytic, and antigen presentation pathways in addition to HSP70, HSP90, and cochaperone proteins, provided evidence for a potential role of OMVs in boosting the host immune response. In conclusion, our work identified novel microbial targets and inherent characteristics of OMVs that could open up new avenues of treatment and prevention of fish infections.


title: Efficacy of antiviral drugs against the omicron variant of SARS-CoV-2.

abstract: The Omicron variant of the SARS-CoV-2 virus was first detected in South Africa in November 2021. The analysis of the sequence data in the context of earlier variants suggested that it may show very different characteristics, including immune evasion and increased transmission. These assumptions were partially confirmed, and the reduction in protection in convalescent patients and vaccinated individuals have been confirmed. Here, we have evaluated the efficacy of antivirals against SARS-CoV-2 variants, Omicron, Delta, and the early 2020 isolate.


title: New cell biological explanations for kinesin-linked axon degeneration

abstract: Axons are the slender, up to meter-long projections of neurons that form the biological cables wiring our bodies. Most of these delicate structures must survive for an organisms lifetime, meaning up to a century in humans. Axon maintenance requires life-sustaining motor protein-driven transport distributing materials and organelles from the distant cell body. It seems logic that impairing this transport causes systemic deprivation linking to axon degeneration. But the key steps underlying these pathological processes are little understood. To investigate mechanisms triggered by motor protein aberrations, we studied more than 40 loss- and gain-of-function conditions of motor proteins, cargo linkers or further genes involved in related processes of cellular physiology. We used one standardised Drosophila primary neuron system and focussed on the organisation of axonal microtubule bundles as an easy to assess readout reflecting axon integrity. We found that bundle disintegration into curled microtubules is caused by the losses of Dynein heavy chain and the Kif1 and Kif5 homologues Unc-104 and Kinesin heavy chain (Khc). Using point mutations of Khc and functional loss of its linker proteins, we studied which of Khcs sub-functions might link to microtubule curling. One cause was emergence of harmful reactive oxygen species through loss of Milton/Miro-mediated mitochondrial transport. In contrast, loss of the Kinesin light chain linker caused microtubule curling through an entirely different mechanism appearing to involve increased mechanical challenge to microtubule bundles through de-inhibition of Khc. The wider implications of our findings for the understanding of axon maintenance and pathology are discussed.

category: cell biology


title: Single-cell RNA sequencing reveals a role for reactive oxygen species and peroxiredoxins in fatty-acid-induced rat β-cell proliferation

abstract: The functional mass of insulin-secreting pancreatic {beta} cells expands to maintain glucose homeostasis in the face of nutrient excess, in part via replication of existing {beta} cells. To decipher the underlying molecular mechanisms, we assessed {beta}-cell proliferation in isolated rat islets exposed to glucose and oleate or palmitate for 48 h and analyzed the transcriptional response by single-cell RNA sequencing. Unsupervised clustering of pooled {beta} cells identified subpopulations, including proliferating {beta} cells. {beta}-cell proliferation increased in response to oleate but not palmitate. Both fatty acids enhanced the expression of genes involved energy metabolism and mitochondrial activity. Comparison of proliferating vs. non-proliferating {beta} cells and pseudotime ordering suggested the involvement of reactive oxygen species (ROS) and peroxiredoxin signaling. Accordingly, the antioxidant N-acetyl cysteine and the peroxiredoxin inhibitor Conoidin A both blocked oleate-induced {beta}-cell proliferation. Our data reveal a key role for ROS signaling in {beta}-cell proliferation in response to nutrients.


title: Baltica: integrated splice junction usage analysis

abstract: Alternative splicing is a tightly regulated co- and post-transcriptional process contributing to the transcriptome diversity observed in eukaryotes. Several methods for detecting differential junction usage (DJU) from RNA sequencing (RNA-seq) datasets exist. Yet, efforts to integrate the results from DJU methods are lacking. Here, we present Baltica, a framework that provides workflows for quality control, de novo transcriptome assembly with StringTie2, and currently 4 DJU methods: rMATS, JunctionSeq, Majiq, and LeafCutter. Baltica puts the results from different DJU methods into context by integrating the results at the junction level. We present Baltica using 2 datasets, one containing known artificial transcripts (SIRVs) and the second dataset of paired Illumina and Oxford Nanopore Technologies RNA-seq. The data integration allows the user to compare the performance of the tools and reveals that JunctionSeq outperforms the other methods, in terms of F1 score, for both datasets. Finally, we demonstrate for the first time that meta-classifiers trained on scores of multiple methods outperform classifiers trained on scores of a single method, emphasizing the application of our data integration approach for differential splicing identification. Baltica is available at https://github.com/dieterich-lab/Baltica under MIT license.