title: A cytoplasmic protein kinase in Chlamydomonas couples engagement of ciliary receptors to rapid cellular responses

abstract: The principal function of the primary cilium is to convert cues from the extracellular milieu into changes in cyclic nucleotide concentration and cytoplasmic responses, but fundamental questions remain about the mechanisms of transmission of cilium-to-cytoplasm signals. During fertilization in Chlamydomonas reinhardtii, ciliary adhesion between plus and minus gametes triggers an immediate [~]10-fold increase in cellular cAMP and activation for cell fusion. Here, we identify Gamete-Specific Protein Kinase (GSPK) as an essential link between cilary receptor engagement and gamete activation. The ciiary adhesion-induced increase in cAMP and cell fusion are severely impaired in gspk mutants but fusion is rescued by a cell-permeable form of cAMP, indicating that GSPK functions upstream of the cAMP increase. GSPK is cytoplasmic, and, remarkably, the entire cellular complement is phosphorylated in less than 60 seconds after ciliary contact. Thus, a cytoplasmic protein kinase rapidly converts a ciliary membrane cue into a global cellular response.

category: cell biology


title: Glucagon Receptor-mediated Regulation of Gluconeogenic Gene Transcription is Endocytosis-dependent in Primary Hepatocytes

abstract: A number of G protein-coupled receptors (GPCRs) are now thought to use endocytosis to promote cellular cAMP signaling that drives downstream transcription of cAMP-dependent genes. We tested if this is true for the Glucagon Receptor (GCGR), which mediates physiological regulation of hepatic glucose metabolism via cAMP signaling. We show that epitope-tagged GCGRs undergo clathrin and dynamin-dependent endocytosis in HEK293 cells after activation by glucagon, and transit via EEA1-marked endosomes shown previously to be sites of GPCR/Gs-stimulated production of cAMP. We further show that endocytosis potentiates cytoplasmic cAMP elevation produced by GCGR activation and promotes transcription of PCK1, the gene which encodes the enzyme catalyzing the rate-limiting step in gluconeogenesis. We verify endocytosis-dependent induction of PCK1 expression by endogenous GCGRs in primary hepatocytes, and show similar control of two other gluconeogenic genes (PGC1 and G6PC). Together, these results implicate the endosomal signaling paradigm in metabolic regulation by glucagon.


title: MAR1 links membrane adhesion to membrane merger during cell-cell fusion in Chlamydomonas

abstract: Union of two gametes to form a zygote is a defining event in the life of sexual eukaryotes, yet the mechanisms that underlie cell-cell fusion during fertilization remain poorly characterized. Here, in studies of fertilization in the green alga, Chlamydomonas, we report identification of a membrane protein on minus gametes, Minus Adhesion Receptor 1 (MAR1), that is essential for the membrane attachment with plus gametes that immediately precedes lipid bilayer merger. We show that MAR1 forms a receptor pair with previously identified receptor FUS1 on plus gametes, whose ectodomain architecture we find is identical to a sperm adhesion protein conserved throughout plant lineages. Strikingly, before fusion, MAR1 is biochemically and functionally associated with the ancient, evolutionarily conserved eukaryotic class II fusion protein HAP2 on minus gametes. Thus, the integral membrane protein MAR1 provides a molecular link between membrane adhesion and bilayer merger during fertilization in Chlamydomonas.


title: A cell atlas of the fly kidney

abstract: Like humans, insects rely on precise regulation of their internal environments to survive. The insect renal system consists of Malpighian tubules and nephrocytes that share similarities to the mammalian kidney. Studies of the Drosophila Malpighian tubules and nephrocytes have provided many insights into our understanding of the excretion of waste products, stem cell regeneration, protein reabsorption, and as human kidney disease models. Here, we analyzed single-nucleus RNA sequencing (snRNA-seq) data sets to characterize the cell types of the adult fly kidney. We identified 11 distinct clusters representing renal stem cells (RSCs), stellate cells (SCs), regionally specific principal cells (PCs), garland nephrocyte cells (GCs) and pericardial nephrocytes (PNs). Analyses of these clusters revealed many new interesting features. For example, we found a new, previously unrecognized cell cluster: lower segment PCs that express Esyt2. In addition, we find that the SC marker genes RhoGEF64c, Frq2, Prip and CG10939 regulate their unusual cell shape. Further, we identified transcription factors specific to each cluster and built a network of signaling pathways that are potentially involved in mediating cell-cell communication between Malpighian tubule cell types. Finally, cross-species analysis allowed us to match the fly kidney cell types to mouse kidney cell types and planarian protonephridia - knowledge that will help the generation of kidney disease models. To visualize this dataset, we provide a web-based resource for gene expression in single cells (https://www.flyrnai.org/scRNA/kidney/). Altogether, our study provides a comprehensive resource for addressing gene function in the fly kidney and future disease studies.

category: developmental biology


title: Wheat EARLY FLOWERING3 is a dawn-expressed circadian oscillator component that regulates heading date

abstract: Optimising the seasonal control of flowering in the major crops is an important component of breeding to match crop adaptation to the target environment. Using an eight parent Multiparent Advanced Generation Inter-Cross (MAGIC) population we investigated the contribution of variation at circadian clock-associated genes to the regulation of heading date in UK and European winter wheat varieties. We identified homoeologues of EARLY FLOWERING 3 (ELF3) as candidate genes for the Earliness per se (Eps) D1 and B1 loci in field conditions. We confirmed that a SNP within the coding region of TaELF3-B1 is the likely causal polymorphism underlying the Eps-B1 locus. We also identified that a reported deletion at the Eps-D1 locus encompassing TaELF3-D1, is in fact a novel allele that lies within an introgression region that contains an inversion relative to the Chinese Spring D genome. Our findings that ELF3 might be associated with the regulation of heading date prompted us to investigate whether ELF3 is a circadian oscillator gene in wheat, as it is in Arabidopsis. Using T. turgidum cv. Kronos carrying loss of function alleles for both copies of TtELF3, we found that circadian rhythms were severely disrupted. Furthermore, in T. aestivum, we found that loss of functional LUX ARRHYTHMO (LUX), an orthologue of the protein partner of ELF3 in Arabidopsis, also severely disrupted circadian rhythms. Whilst these data suggest a function for both ELF3 and LUX in the wheat circadian oscillator, that oscillator might be structured differently to that of Arabidopsis because wheat ELF3 and LUX transcripts are maximal at the end of the night and day respectively, rather than co-expressed at dusk as they are in Arabidopsis. We conclude that there is sufficient allelic diversity within the three wheat ELF3 homoeologues for selection to delay or advance heading, and that this can be achieved without pleiotropic deleterious alterations to circadian rhythms.

category: plant biology


title: Assessing the effectiveness of honey bee pollinators for cultivated blueberries in South Africa

abstract: Globally, agricultural crops are often dependent on insect pollination. Blueberries are an example of such a crop and owing to their proposed health benefits they are grown around the world, including locations where their native bumble bee pollinators do not occur. In the absence of bumble bees, blueberry pollination in South Africa and many other parts of the non-native, commercial range is performed primarily by honey bees. Despite this, the effectiveness of honey bee pollination on blueberries remains understudied. This study determined the effect of honey bee pollination on components of fruit yield (fruit set and mass) of five blueberry varieties that are extensively planted in South Africa. For each variety, two metrics were calculated: 1) the benefit of bees -- a comparison of fruit yields after exposure to honey bees and fruit yields after honey bee exclusion, 2) the pollination deficit -- the difference in yield between hand pollination (maximum yield potential) and yields after exposure to honey bee pollinators. Honey bee pollination consistently resulted in improved yields, although the magnitude of this improvement (i.e., the benefit of bees) was dependent on the variety considered. Similarly, the pollination deficit also varied considerably across varieties and while some varieties appeared to perform close to maximum potential (small pollination deficit), others yielded well below their maximum potential under honey bee pollination. This study demonstrates that honey bees are functional pollinators of blueberries in areas where native blueberry pollinators are absent. However, in such areas, it is important that special focus be given to selecting blueberry varieties that perform well with honey bees as their sole pollinator. Further research is necessary to determine how the pollination deficit of blueberry varieties can be decreased as well as how to increase the effectiveness of honey bee pollination.


title: Genomic evidence for divergent co-infections of SARS-CoV-2 lineages

abstract: Recently, patients co-infected by two SARS-CoV-2 lineages have been sporadically reported. Concerns are raised because previous studies have demonstrated co-infection may contribute to the recombination of RNA viruses and cause severe clinic symptoms. In this study, we have estimated the compositional lineage(s), tendentiousness, and frequency of co-infection events in population from a large-scale genomic analysis for SARS-CoV-2 patients. SARS-CoV-2 lineage(s) infected in each sample have been recognized from the assignment of within-host site variations into lineage-defined feature variations by introducing a hypergeometric distribution method. Of all the 29,993 samples, 53 (~0.18%) co-infection events have been identified. Apart from 52 co-infections with two SARS-CoV-2 lineages, one sample with co-infections of three SARS-CoV-2 lineages was firstly identified. As expected, the co-infection events mainly happened in the regions where have co-existed more than two dominant SARS-CoV-2 lineages. However, co-infection of two sub-lineages in Delta lineage were detected as well. Our results provide a useful reference framework for the high throughput detecting of SARS-CoV-2 co-infection events in the Next Generation Sequencing (NGS) data. Although low in average rate, the co-infection events showed an increasing tendency with the increased diversity of SARS-CoV-2. And considering the large base of SARS-CoV-2 infections globally, co-infected patients would be a nonnegligible population. Thus, more clinical research is urgently needed on these patients.

category: microbiology


title: Network analysis outlines strengths and weaknesses of emerging SARS-CoV-2 Spike variants

abstract: The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to dissect and understand the structure and dynamics of this complex pathogen. The Spike glycoprotein of SARS-CoV-2 has received special attention as it is the means by which the virus enters the human host cells. The N-terminal domain (NTD) is one of the targeted regions of the Spike protein for therapeutics and neutralizing antibodies against COVID-19. Though its function is not well-understood, the NTD is reported to acquire mutations and deletions that can accelerate the evolutionary adaptation of the virus driving antibody escape. Cellular processes are known to be regulated by complex interactions at the molecular level, which can be characterized by means of a graph representation facilitating the identification of key residues and critical communication pathways within the molecular complex. From extensive all-atom molecular dynamics simulations of the entire Spike for the wild-type and the dominant variant, we derive a weighted graph representation of the protein in two dominant conformations of the receptor-binding-domain; all-down and one-up. We implement graph theory techniques to characterize the relevance of specific residues at facilitating roles of communication and control, while uncovering key implications for fitness and adaptation. We find that many of the reported high-frequency mutations tend to occur away from the critical residues highlighted by our graph theory analysis, implying that these mutations tend to avoid targeting residues that are most critical for protein allosteric communication. We propose that these critical residues could be candidate targets for novel antibody therapeutics. In addition, our analysis provides quantitative insights of the critical role of the NTD and furin cleavage site and their wide-reaching influence over the protein at large. Many of our conclusions are supported by empirical evidence while others point the way towards crucial simulation-guided experiments.

category: molecular biology


title: Distinct stress-dependent signatures of cellular and extracellular tRNA-derived small RNAs (tDRs)

abstract: The cellular response to stress is an important determinant of disease pathogenesis. Uncovering the molecular fingerprints of distinct stress responses may yield novel biomarkers for different diseases, and potentially identify key signaling pathways important for disease progression. tRNAs and tRNA-derived small RNAs (tDRs) comprise one of the most abundant RNA species in cells and have been associated with cellular stress responses. The presence of RNA modifications on tDRs has been an obstacle for accurately identifying tDRs with conventional small RNA sequencing. Here, we use AlkB-facilitated methylation sequencing (ARM-seq) to uncover a comprehensive landscape of cellular and extracellular tDR expression in a variety of human and rat cells during common stress responses, including nutritional deprivation, hypoxia, and oxidative stress. We found that extracellular tDRs have a distinct fragmentation signature with a predominant length of 31-33 nts and a highly specific termination position when compared with intracellular tDRs. Importantly, we found these signatures are better discriminators of different cellular stress responses compared to extracellular miRNAs. Distinct extracellular tDR signatures for each profiled stressor are elucidated in four different types of cells. This distinct extracellular tDR fragmentation pattern is also noted in plasma extracellular RNAs from patients on cardiopulmonary bypass. The observed overlap of these patient tDR signatures with the signatures of nutritional deprivation and oxidative stress in our cellular models provides preliminary in vivo corroboration of our findings and demonstrates the potential to establish novel extracellular tDR biomarkers in human disease models.


title: A separation-of-function ZIP4 wheat mutant allows crossover between related chromosomes and is meiotically stable

abstract: Many species, including most flowering plants, are polyploid, possessing multiple genomes. During polyploidisation, fertility is preserved via the evolution of mechanisms to control the behaviour of these multiple genomes during meiosis. On the polyploidisation of wheat, the major meiotic gene ZIP4 duplicated and diverged, with the resulting new gene TaZIP4-B2 being inserted into chromosome 5B. Previous studies showed that this TaZIP4-B2 promotes pairing and synapsis between wheat homologous chromosomes, whilst suppressing crossover between related (homoeologous) chromosomes. Moreover, in wheat, the presence of TaZIP4-B2 preserves up to 50% of grain number. The present study exploits a  separation-of-function wheat Tazip4-B2 mutant named zip4-ph1d, in which the Tazip4-B2 copy still promotes correct pairing and synapsis between homologues (resulting in the same pollen profile and fertility normally found in wild type wheat), but which also allows crossover between the related chromosomes in wheat haploids of this mutant. This suggests an improved utility for the new zip4-ph1d mutant line during wheat breeding exploitation, compared to the previously described CRISPR Tazip4-B2 and ph1 mutant lines. The results also reveal that loss of suppression of homoeologous crossover between wheat chromosomes does not in itself reduce wheat fertility when promotion of homologous pairing and synapsis by TaZIP4-B2 is preserved.

category: genetics


title: Structures of the human peroxisomal fatty acid transporter ABCD1 in a lipid environment

abstract: The peroxisomal very long chain fatty acid (VLCFA) transporter ABCD1 is central to cellular fatty acid catabolism and lipid biosynthesis. Its dysfunction underlies toxic cytosolic accumulation of VLCFAs, progressive nervous system demyelination, and neurological impairments including the potentially fatal disease X-linked adrenoleukodystrophy (X-ALD). Molecular details underlying substrate recognition and transport by ABCD1 are poorly understood. Here we determined cryo-EM structures of ABCD1 in phospholipid nanodiscs in a nucleotide bound conformation open to the peroxisomal lumen and an inward facing conformation open to the cytosol at up to 3.5 [A] resolution that reveal key details of its transmembrane cavity and ATP dependent conformational transitions. We identify structural elements distinguishing ABCD1 from its closest homologs and show that coenzyme A (CoA) esters of VLCFAs modulate ABCD1 activity in a species dependent manner. Together, our data support a transport mechanism where only the CoA moieties of VLCFA-CoAs enter the hydrophilic transmembrane cavity while the acyl chains extend out into the surrounding membrane bilayer, help rationalize disease causing mutations, and provide a framework for ABCD1 targeted structure-based drug design.

category: biochemistry


title: An Extended Coatomer Binding Motif in the SARS-CoV-2 Spike Protein

abstract: {beta}-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for COPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. COPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking.


title: An Extended Motif in the SARS-CoV-2 Spike Modulates Binding and Release of Host Coatomer in Retrograde Trafficking


title: DOCKGROUND scoring benchmarks for protein docking

abstract: Protein docking protocols typically involve global docking scan, followed by re-ranking of the scan predictions by more accurate scoring functions that are either computationally too expensive or algorithmically impossible to include in the global scan. Development and validation of scoring methodologies are often performed on scoring benchmark sets (docking decoys) which offer concise and nonredundant representation of the global docking scan output for a large and diverse set of protein-protein complexes. Two such protein-protein scoring benchmarks were built for the DO_SCPLOWOCKGROUNDC_SCPLOW resource, which contains various datasets for the development and testing of protein docking methodologies. One set was generated based on the DO_SCPLOWOCKGROUNDC_SCPLOW unbound docking benchmark 4, and the other based on protein models from the DO_SCPLOWOCKGROUNDC_SCPLOW model-model benchmark 2. The docking decoys were designed to reflect the reality of the real-case docking applications (e.g., correct docking predictions defined as near-native rather than native structures), and to minimize applicability of approaches not directly related to the development of scoring functions (reducing clustering of predictions in the binding funnel and disparity in structural quality of the near-native and non-native matches). The sets were further characterized by the source organism and the function of the protein-protein complexes. The sets, freely available to the research community on the DO_SCPLOWOCKGROUNDC_SCPLOW webpage, present a unique, user-friendly resource for the developing and testing of protein-protein scoring approaches.

category: bioinformatics


title: Structural anatomy of C1 domain interactions with DAG and other agonists

abstract: Diacylglycerol (DAG) is a versatile lipid whose 1,2-sn-stereoisomer serves both as second messenger in signal transduction pathways that control vital cellular processes, and as metabolic precursor for downstream signaling lipids such as phosphatidic acid1,2. DAG effector proteins compete for available lipid using conserved homology 1 (C1) domains as DAG-sensing modules. Yet, how C1 domains recognize and capture DAG in the complex environment of a biological membrane has remained unresolved for the 40 years since the discovery of Protein Kinase C (PKC)3 as the first member of the DAG effector cohort. Herein, we report the first high-resolution crystal structures of a C1 domain (C1B from PKC{delta}) complexed to DAG and to each of four potent PKC agonists that produce different biological readouts and that command intense therapeutic interest. This structural information details the mechanisms of stereospecific recognition of DAG by the C1 domains, the functional properties of the lipid-binding site, and the identities of the key residues required for the recognition and capture of DAG and exogenous agonists. Moreover, the structures of the five C1 domain complexes provide the high-resolution guides for the design of agents that modulate the activities of DAG effector proteins.

category: biophysics


title: GLS1 is a Protective Factor Rather than a Molecular Target in ARID1A-Mutated Ovarian Clear Cell Carcinoma

abstract: Targeting glutamine metabolism has emerged as a novel therapeutic strategy for several human cancers, including ovarian cancer. The primary target of this approach is the kidney isoform of glutaminase, glutaminase 1 (GLS1), a key enzyme in glutamine metabolism that is overexpressed in several human cancers. A first-in-class inhibitor of GLS1, called CB839 (Telaglenastat), has been investigated in several clinical trials, with promising results. The first clinical trial of CB839 in platinum-resistant ovarian cancer patients is forthcoming. ARID1A-mutated ovarian clear cell carcinoma (OCCC) is a relatively indolent and chemoresistant ovarian cancer histotype. In OCCC-derived cells in vitro and mouse models, loss of ARID1A leads to upregulation of GLS1. Thus, targeting of GLS1 with CB839 has been suggested as a targeted approach for OCCC patients with tumors harboring ARID1A-mutations. Here, we investigated whether GLS1 is differentially expressed between OCCC patients whose tumors are ARID1A positive and patients whose tumors are ARID1A negative. In clinical specimens of OCCC, we found that GLS1 overexpression was not correlated with ARID1A loss. In addition, GLS1 overexpression was associated with better clinical outcomes. Our findings suggest that GLS1 expression in OCCC may be a protective factor and that caution should be taken when considering the use of CB839 to treat OCCC patients.

category: cancer biology


title: Tbet expression by Tregs is needed to protect against Th1-mediated immunopathology during Toxoplasma infection in mice

abstract: T. gondii infection has proven to be an ideal model to understand the delicate balance between protective immunity and immune-mediated pathology during infection. Lethal infection causes a collapse of Tregs mediated by loss of IL-2, and conversion of Tregs to IFN{gamma} producing cells. Importantly, these Tregs highly express the Th1 transcription factor Tbet. To determine the role of Tbet in Tregs, we infected Tbx21f/f-Foxp3YFPCre and control Foxp3YFPCre mice with the type II strain of T. gondii, ME49. The majority of Tbx21f/f-Foxp3YFPCre mice succumb to a non-lethal acute infection. Notably, parasite burden is comparable between Tbx21f/f-Foxp3YFPCre and Foxp3YFPCre control mice. We found that Tbx21f/f-Foxp3YFPCre mice have significantly higher serum levels of proinflammatory cytokines IFN{gamma} and TNF, suggestive of a heightened immune response. To test if CD4+ T cells were driving immunopathology, we treated Tbx21f/f-Foxp3YFPCre mice with anti-CD4 depleting antibody and partially rescued these mice. Broad spectrum antibiotic treatment also improved survival, demonstrating a role for commensal flora in immunopathology in Tbx21f/f-Foxp3YFPCre mice. RNA-seq analysis reinforced that Tbet regulates several key cellular pathways, including chromosome segregation, cytokine receptor activity and cell cycle progression, that help to maintain fitness in Tregs during Th1 responses. Taken together, our data shows an important role for Tbet in Tregs in preventing lethal immunopathology during Toxoplasma gondii infection, further highlighting the protective role of Treg plasticity to self and microbiota.

category: immunology


title: AbmR is a mycobacterial dual-function transcription factor and ribonucleoprotein with distinct DNA and RNA-binding determinants

abstract: The bacterium Mycobacterium tuberculosis (Mtb) must adapt to myriad host-associated stressors. A recently identified transcription factor, AbmR (ATP-binding mcr11-regulator), regulates expression of an essential stress-responsive small RNA (Mcr11) and inhibits the growth of Mtb. Previously, AbmR was found to make 39S complexes of unknown function. Here we report that AbmR 39S complexes are comprised of AbmR and co-purifying RNAs and that RNA-binding inhibits AbmRs DNA-binding function. While AbmR binds DNA and regulates gene expression in a sequence specific manner, RNA-binding is not sequence specific. Amino acid R146 is important for DNA-binding but completely dispensable for RNA-binding and 39S complex formation, establishing that the RNA- and DNA-binding functions of AbmR are distinct. RNA bound by AbmR was protected from RNase digestion, supporting an RNA modulatory function for the 39S complex. We also found that abmR is required for optimal survival during treatment with the ATP-depleting antibiotic bedaquiline, which is associated with extended RNA stability. These data establish a paradigm wherein a transcription factor assembles into large complexes to transition between mutually exclusive DNA-binding gene regulatory and RNA-binding RNA modulatory functions. Our findings indicate that AbmR is a dual-function protein that may have novel RNA regulatory roles in stress adapted Mtb.


title: When development is constantly but weakly perturbed - Canalization by microRNAs

abstract: Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of GRN (gene regulatory network) stabilization as well as developmental canalization against constant but small environmental perturbations. Since the complete removal of miRNAs is lethal, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs. We hypothesize that flies with modest miRNA reductions will gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by [~] 90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. We further show that the developmental failure is associated with GRN aberration in the larval stages even before phenotypic aberrations become observable. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be the genomes solution to weak but constant environmental perturbations.


title: Tissue scale properties of the extracellular matrix regulates nuclear shape, organisation and fate in the embryonic midline sutures

abstract: Balancing self-renewal and differentiation is a key feature of every stem cell niche and one that is tuned by mechanical interactions of cells with their neighbors and surrounding extracellular matrix. The fibrous stem cell niches that develop as sutures between skull bones must balance the complex extracellular environment that emerges to define them with self-renewal and bone production. Here, we address the role for physical stimuli in suture development by probing the relationship between nuclear shape, organization and gene expression in response to a developing collagen network in embryonic midline sutures. This work complements genetic approaches used to study sutures and provides the first quantitative analyses of physical structure in these sutures. By combining multiple imaging modalities with novel shape description, in addition to network analysis methods, we find the early emergence of a complex extracellular collagen network to have an important role in regulating morphogenesis and cell fate. We show that disrupted collagen crosslinking can alter ECM organization of midline sutures as well as stimulate expression of bone differentiation markers. Further, our findings suggest that in vivo, skeletal tissues can uncouple the response of the nuclear lamina from collagen mediated tissue stiffening seen in vitro. Our findings highlight a crucial relationship between the cellular microenvironment, tissue stiffness and geometry with gene expression in normal development and maintenance of progenitor fate in embryonic sutures.


title: Production of offspring from azoospermic mice with meiotic failure: Precise biparental meiosis within halved oocytes

abstract: While the large volume of mammalian oocytes is necessary for embryo development, it can lead to error-prone chromosomal segregation during meiosis. Conversely, we hypothesized that smaller oocytes would have a great unidentified potential to stabilize unstable meiosis and improve the development of the resultant embryos. Here, we show that reducing ooplasmic volume can rescue highly error-prone fertilization using primary spermatocytes by preventing segregation errors of chromosomes during biparental meiosis. High-resolution live-imaging analysis revealed that erroneous chromosome segregation occurred in most (90%) spermatocyte-injected oocytes of normal size, but could be ameliorated to 40% in halved oocytes. The birth rate improved remarkably from 1% to 19% (P < 0.0001). Importantly, this technique enabled the production of offspring from azoospermic mice with spermatocyte arrest caused by STX2 deficiency, an azoospermia factor also found in humans. Thus, contrary to popular opinion, oocytes inherently have a strong potential for precise meiotic divisions, which can be evoked by reduction of the ooplasmic volume. Their potential might help rescue cases of untreatable human azoospermia with spermatocyte arrest.


title: Production of Offspring from Azoospermic Mice with Meiotic Failure: Precise Biparental Meiosis within Halved Oocytes


title: Evolution and functional dynamics of dehydrins in model Brachypodium grasses

abstract: Dehydration proteins (dehydrins, DHNs) confer tolerance to water-stress deficit to plants, thus playing a fundamental role in plant response and adaptation to water-deprivation stressful environments. We have performed a comparative genomics and evolutionary study of DHN genes in four model Brachypodium grass species, and a drought-induced functional analysis in 32 ecotypes of the flagship species B. distachyon, to gain insight into the origins and dynamics of these proteins and the correlated drought-mediated phenotypic responses in ecotypes showing different hydric requirements. Genomic sequence analysis detected 10 types of dehydrin genes (Bdhn) across the Brachypodium species, totalling 47 genes. Domain and conserved motif contents of peptides encoded by Bdhn genes revealed eight protein architectures, YS{phi}K2 being the most common architecture. Bdhn genes were spread across several chromosomes and more frequent in syntenic chromosomes 3 and 4 of B. distachyon, 4 and 5 of B. stacei and 4 of B. sylvaticum. Tandem and segmental duplication events were detected for four Bdhn genes. Selection analysis indicated that all the Bdhn genes were constrained by purifying selection. Three upstream cis-regulatory motifs (BES1, MYB124, ZAT) were consistently detected in several Bdhn genes. Functional analysis in 32 natural accessions of B. distachyon demonstrated that only four Bdhn genes (Bdhn1, Bdhn2, Bdhn3, Bdhn7) were expressed in mature leaves and that all of them were significantly more highly expressed in plants under drought conditions. These genes corresponded to wheat orthologs that were also significantly more expressed under drought stress. Brachypodium dehydrin expression was significantly correlated with drought-response phenotypic traits (plant biomass, leaf carbon and proline contents and WUE increases, leaf water and nitrogen content changes) being more pronounced in drought-tolerant ecotypes. Bdhn expression, associated phenotypic trait changes and climate niche variation did not show significant phylogenetic signal when tested in the B. distachyon genealogical-species tree. By contrast, some of them showed low or marginal significant phylogenetic signal when tested in the B. distachyon Bdhn tree, suggesting that Bdhn gene evolution is partially related to adaptation to drought in this species. Our results demonstrate that dehydrin composition and regulation is a key factor determining the acquisition of water-stress tolerance in grasses.


title: Modeling CRISPR-Cas13d on-target and off-target effects using machine learning approaches

abstract: A major challenge in the application of the CRISPR-Cas13d (RfxCas13d, or CasRx) RNA editing system is to accurately predict its guide RNA (gRNA) dependent on-target and off-target effect. Here, we performed CRISPR-Cas13d proliferation screens that target protein-coding genes and long non-coding RNAs (lncRNAs), followed by a systematic modeling of Cas13d on-target efficiency and off-target viability effect. We first designed a deep learning model, named DeepCas13, to predict the on-target activity of a gRNA with high accuracy from its sequence and secondary structure. DeepCas13 outperforms existing methods and accurately predicts the efficiency of guides targeting both protein-coding and non-coding RNAs (e.g., circRNAs and lncRNAs). Next, we systematically studied guides targeting non-essential genes, and found that the off-target viability effect, defined as the unintended effect of guides on cell viability, is closely related to their on-target RNA cleavage efficiency. This finding suggests that these gRNAs should be used as negative controls in proliferation screens to reduce false positives, possibly coming from the unwanted off-target viability effect of efficient guides. Finally, we applied these models to our screens that included guides targeting 234 lncRNAs, and identified lncRNAs that affect cell viability and proliferation in multiple cell lines. DeepCas13 is freely accessible via http://deepcas13.weililab.org.


title: KIRCLE: An Analysis of Variations in KIR Genes in The Cancer Genome Atlas and UK Biobank

abstract:
BackgroundNatural killer (NK) cells represent a critical component of the innate immune systems response against cancer and viral infections, among other diseases. To distinguish healthy host cells from infected or tumor cells, killer immunoglobulin receptors (KIR) on NK cells bind and recognize Human Leukocyte Antigen (HLA) complexes on their target cells. Just like the HLAs they bind, these KIRs exhibit high allelic diversity in the human population.

ResultsIn order to better understand these immunoreceptors, we have developed KIRCLE, a novel method for genotyping individual KIR genes from whole exome sequencing data, and used it to analyze approximately 60,000 patient samples in The Cancer Genome Atlas and UK Biobank. We were able to assess population frequencies for different KIR alleles and demonstrate that, similar to HLA alleles, individuals KIR alleles correlate strongly with their ethnicities. In addition, we observed associations between different KIR alleles and HLA alleles, including HLA-B*53 with KIR3DL2*013 (Fishers Exact FDR = 7.64e-51). Finally, we showcased statistically significant associations between KIR alleles and various clinical correlates, including peptic ulcer disease (Fishers Exact FDR = 0.0429) and age of onset of atopy and various KIR alleles (Mann-Whitney-U FDR = 0.0751).

ConclusionsKIR polymorphism and NK cells play a critical role in many diseases, often through their interactions with HLA complexes. Peptic ulcer disease and atopy are just two diseases in which NK cells may play a role beyond their "classical" realm of anti-tumor and anti-viral responses.


title: MOSCATO: A Supervised Approach for Analyzing Multi-Omic Single-Cell Data

abstract: Advancements in genomic sequencing continually improve personalized medicine in complex diseases. Recent breakthroughs generate multiple types of signatures (or multi-omics) from each cell, producing different data  omic types per single-cell experiment. We introduce MOSCATO, a technique for selecting features across multi-omic single-cell datasets that relate to clinical outcomes. For example, we leverage penalization concepts often used in multi-omic network analytics to accommodate the high-dimensionality where multiple-testing is likely underpowered. We organize the data into multi-dimensional tensors where the dimensions correspond to the different  omic types. Using the outcome and the single-cell tensors, we perform regularized tensor regression to return a variable set for each  omic type that forms the clinically-associated network. Robustness is assessed over simulations based on available single-cell simulation methods. Real data comparing healthy subjects versus subjects with leukemia is also considered in order to identify genes associated with the disease. The flexibility of our approach enables future extensions on distributional assumptions and covariate adjustments. This algorithm may identify clinically-relevant genetic patterns on a cellular-level that span multiple layers of sequencing data and ultimately inform highly precise therapeutic targets in complex diseases. Code to perform MOSCATO and replicate the real data application is publicly available on GitHub at https://github.com/lorinmil/MOSCATO and https://github.com/lorinmil/MOSCATOLeukemiaExample.


title: Protein Structure Prediction Using a Maximum Likelihood Formulation of a Recurrent Geometric Network

abstract: Only [~]40% of the human proteome has structural coordinates available from experiment (i.e., X-ray crystallography, NMR spectroscopy, or cryo-EM) or homology modeling with quality templates (i.e., 30% sequence identity or greater), leaving most of the proteome structurally unsolved. Deep learning (DL) methods for predicting protein structure can help close knowledge gaps where experimental and homology models are difficult to obtain. Recent advances in these DL methods have shown promising results in expanding structural coverage to the scale of the entire human proteome, providing researchers with more complete protein structural information. Here, we improve upon an existing DL algorithm for protein structure prediction, the Recurrent Geometric Network (RGN). We first expand the training dataset to include experimental uncertainty data in the form of atomic displacement parameters, then derive a maximum likelihood loss function that incorporates this uncertainty data into model training. Compared to the original RGN, our novel maximum likelihood model improves the rate of convergence of initial model training and ultimately results in more accurate structure prediction according to the root mean square deviation (RMSD) of backbone atoms, the Global Distance Test (GDT), the Global Distance Test High Accuracy (GDT-HA), and the Template-Modeling Score (TM-Score). Our model also predicts structures with more favorable backbone torsions, which provide more accurate starting coordinates for downstream physics-based simulations. Based on these results, our maximum likelihood reformulation provides a framework for improving existing or future machine learning algorithms for protein structure prediction. The augmented dataset, data collection scripts, reformulated RGN source code, and a series of trained models are publicly available at https://github.com/SchniedersLab/likelihood-rgn.


title: Enzyme Aggregation and Fragmentation Induced by Catalysis Relevant Species

abstract: It is usually assumed that enzymes retain their native structure during catalysis. However, the aggregation and fragmentation of proteins can be difficult to detect and sometimes conclusions are drawn based on the assumption that the protein is in its native form. We have examined three model enzymes, alkaline phosphatase (AkP), hexokinase (HK) and glucose oxidase (GOx). We find that these enzymes aggregate or fragment after addition of chemical species directly related to their catalysis. We used several independent techniques to study this behavior. Specifically, we found that glucose oxidase and hexokinase fragment in the presence of D-Glucose but not L-glucose, while hexokinase aggregates in the presence of Mg2+ ion and either ATP or ADP at low pH. Alkaline phosphatase aggregates in the presence of Zn2+ ion and inorganic phosphate. The aggregation of hexokinase and alkaline phosphatase does not appear to attenuate their catalytic activity. Our study indicates that specific multimeric structures of native enzymes may not be retained during catalysis and suggests pathways for different enzymes to associate or separate over the course of substrate turnover.


title: CaV1.3 enhanced store operated calcium promotes resistance to androgen deprivation in prostate cancer

abstract:
Androgen deprivation therapy (ADT) is the main treatment for advanced prostate cancer (PCa) but resistance results in progression to terminal castrate resistant PCa (CRPC), where there is an unmet therapeutic need. Aberrant intracellular calcium (Cai2+) is known to promote neoplastic transformation and treatment resistance. There is growing evidence that expression of voltage gated calcium channels (VGCC) is increased in cancer, particularly the CACNA1D/CaV1.3 in CRPC. The aim of this study was to investigate if increased CaV1.3 drives resistance to ADT and determine its associated impact on Cai2+ and cancer biology.

Bioinformatic analysis revealed that CACNA1D gene expression is increased in ADT treated PCa patients regardless of TMPRSS2:ERG status. Corroborated in both in vivo LNCaP xenograft mouse and in vitro PCa cell line models which demonstrated a significant increase in CaV1.3 protein expression following ADT with bicalutamide. The expression was found to be a shortened 170kDA CaV1.3 isoform associated which failed to mediate calcium influx following membrane depolarisation. Instead, under ADT CaV1.3 mediated a rise in basal cytosolic calcium and an increase in store operated calcium entry (SOCE). This in turn drove both proliferation and survival of long-term ADT CRPC cells.

Overall, this study demonstrates for the first time in PCa that increased SOCE through a novel CaV1.3 mechanism which represents a novel oncogenic switch that contributes to ADT resistance and promotes CRPC biology. Highlighting aberrant intracellular calcium in CRPC as a potential area for therapeutic development to improve patient outcomes.


title: Numb is required for optimal contraction of skeletal muscle

abstract:
BackgroundThe role of Numb, a protein that is important for cell fate and development was investigated in adult skeletal muscle in mice using a conditional, inducible knockout (cKO) model.

MethodsNumb expression was evaluated by Western blot. Numb localization was determined by confocal microscopy. The effects of cKO of Numb and the closely-related gene Numb-like in skeletal muscle fibers was evaluated by in-situ physiology; transmission and focused ion beam scanning electron microscopy; 3-dimensional reconstruction of mitochondrial; lipidomics; and bulk RNA-sequencing. Additional studies using primary mouse myotubes investigated the effects the effects of Numb knockdown on cell fusion, mitochondrial function and calcium transients.

ResultsNumb protein expression was reduced by [~]70% (p < 0.01) at 24 as compared to 3 months of age. Numb was localized within muscle fibers as bands traversing fibers at regularly spaced intervals in close proximity to dihydropyridine receptors.

The cKO of Numb and Numb-like reduced specific tetanic force by 36%, p < 0.01), altered mitochondrial spatial relationships to sarcomeric structures, increased Z-line spacing by 30% (p < 0.0001), perturbed sarcoplasmic reticulum organization and reduced mitochondrial volume by over 80% (p < 0.01). Only six genes were differentially expressed in cKO mice: Itga4, Sema7a, Irgm2, Vezf1, Mib1 and Tmem132a. Several lipid mediators derived from polyunsaturated fatty acid (PUFAs) through lipoxygenases were upregulated in Numb cKO skeletal muscle; 12-HEPE was increased by [~]250% (p < 0.05) and 17,18-EpETE by [~]240% (p < 0.05). In mouse primary myotubes, Numb knock-down reduced cell fusion ([~]20%, p < 0.01) and mitochondrial function and delayed the caffeine-induced rise in cytosolic calcium concentrations by more than 100% (p < 0.01).

ConclusionsThese findings implicate Numb as a critical factor in skeletal muscle structure and function which appear to be critical for calcium release; we therefore speculate Numb plays critical roles in excitation-contraction coupling, one of the putative targets of aged skeletal muscles. These findings provide new insights into the molecular underpinnings of the loss of muscle function observed with sarcopenia.

category: physiology


title: Inexpensive and easy method for 6 fragment Golden Gate Assembly of a modular S/MARs mammalian expression vector and its variants

abstract:
BackgroundA basic requirement for synthetic biology is the availability of efficient DNA assembly methods. Numerous methods have been previously reported to accomplish this task. One such method has been reported, which allows parallel assembly of multiple DNA fragments in a one-tube reaction, called Golden Gate Assembly. This study aims to further simplify that method and make it more suitable for small labs and students.

MethodsPrior to amplification of the parental plasmids used in building the modules were domesticated using a variation of SDM (Site Directed Mutagenesis) called SPRIP. After careful design and amplification of the desired modules, using a high-fidelity polymerase, amplified PCR fragments that enter the one-step-one-pot reaction were stored in Zymo DNA/RNA Shield at -20 degrees C and thawed whenever needed to be used as fragments or modules in the assembly. The fragments were designed to posses unique overhangs using NEB Golden Gate assembly tool and Snapgene, amplification of modules was performed using a Q5 high fidelity polymerase from preexisting plasmids or gene fragments, clean-up of the PCR products (fragments) was performed in one tube per assembly using Zymo DNA Clean and Concentrator-5, assembled using BsaI and T4 ligase, DpnI digestion performed for eliminating the background plasmids that remain after the PCR reaction and the resulting assembled product was transformed into competent E.coli cells. Transformants were screened using diagnostic digest, transfected into HEK293T cells and the fluorescence was evaluated using fluorescent microscopy and flow cytometry.

ResultsHerein presented is a simple and inexpensive alternate protocol to build modular plasmids using the Golden Gate Assembly method. A total of p37 S/MARs mammalian expression vectors were designed and constructed using 6 modules previously amplified by PCR and stored in the appropriate buffer to eliminate exo- and endonuclease activity and to protect the DNA from freeze thaw cycles. The existing modules were interchangeable and new modules were easily amplified and stored for use when needed. The mammalian expression vectors constructed showed the desired restriction pattern and GFP expression in bacteria and in mammalian cells. A comparison of 7 pNoname variants was conducted using flow cytometry. Interestingly, no pNoname variant harbouring the SV40 promoter showed expression in tested HEK293T cells. It appears that using the Ef1a promoter in combination with the BGH polyA signal provides the best expression in S/MARS vectors harboring the DTS40 region, as measured by flow cytometry.

ConclusionsProvided the design steps are respected and the fragments are stored and labeled appropriately, multiple plasmid variants and combinations of the pre-designed modules can be assembled in one day, easier and using less resources than the established protocols, with good efficiency. The simplicity of the design and the affordability of the method could make modular cloning of plasmid constructs more accessible to small labs and students.

category: synthetic biology


title: In lupus nephritis, specific in situ inflammatory states are associated with refractory disease and progression to renal failure

abstract:
In human lupus nephritis (LN), tubulointerstitial inflammation (TII) on biopsy predicts progression to end stage renal disease (ESRD). However, while approximately half of patients with moderate or severe TII develop ESRD, half do not. Therefore, we hypothesized that TII is heterogeneous, with distinct inflammatory states associated with different renal outcomes. We interrogated renal biopsies from LN longitudinal and cross-sectional cohorts using both conventional and highly multiplexed confocal microscopy. To accurately segment cells across whole biopsies, and to understand their spatial relationships, we developed unique computational pipelines by training and implementing several deep learning models and other computer vision techniques. Surprisingly, across biopsies, high B cell densities were strongly associated with protection from ESRD. In contrast, elevated CD4-T cell population densities, which included CD8, {gamma}{delta} and double negative (CD4-CD8-{delta}-, DN) T cells, were associated with both acute refractory renal failure and gradual progression to ESRD. Interestingly, lymphocytes and dendritic cells were organized into discrete clusters or neighborhoods that could be characterized by the enrichment for specific cell populations. B cells were often organized into large neighborhoods with CD4+ T cells including T follicular helper-like cells. In contrast, the CD4-T cell populations formed small cellular neighborhoods whose frequency predicted subsequent progression to ESRD. These data reveal that in LN, specific in situ inflammatory states are associated with refractory disease and progression to ESRD.

One sentence summaryUsing deep machine learning to analyze confocal microscopy data, we demonstrate that in lupus nephritis, CD4-T cell populations, including CD8+ and {gamma}{delta} T cells, organize into specific spatial neighborhoods that predict progression to renal failure.


title: Aberrant localization of CDC42 C-terminal variants to the Golgi apparatus drives pyrin inflammasome-dependent autoinflammation

abstract: Mutations in the C-terminal region of the CDC42 gene cause severe neonatal-onset autoinflammation. Elevated levels of serum IL-18 in patients and effectiveness of IL-1{beta}-blocking therapy indicate that the pathology involves abnormal inflammasome activation; however, the mechanism underlying autoinflammation remains to be elucidated. Using induced-pluripotent stem cells established from patients carrying CDC42R186C, we found that patient-derived cells secreted larger amounts of IL-1{beta} in response to pyrin-activating stimuli. Aberrant palmitoylation and localization of CDC42R186C protein to the Golgi apparatus promoted pyrin inflammasome assembly downstream of pyrin dephosphorylation. Aberrant subcellular localization was the common pathological feature shared by CDC42 C-terminal variants with inflammatory phenotypes, including CDC42*192C*24 that also localizes to the Golgi apparatus. Furthermore, the level of pyrin inflammasome overactivation paralleled that of mutant protein accumulation in the Golgi apparatus, but no that of the mutant GTPase activity. These results reveal an unexpected association between CDC42 subcellular localization and pyrin inflammasome activation that could pave way for elucidating the mechanism of pyrin inflammasome formation.


title: Chromosomal-level genome assembly of the bioluminescent cardinalfishSiphamia tubifer, an emerging model for symbiosis research

abstract: The bioluminescent symbiosis between the sea urchin cardinalfish Siphamia tubifer (Kurtiformes: Apogonidae) and the luminous bacterium Photobacterium mandapamensis is an emerging vertebrate-bacteria model for the study of microbial symbiosis. However, there is little genetic data available for the host fish, limiting the scope of potential research that can be carried out with this association. In this study, we present a chromosomal-level genome assembly of S. tubifer using a combination of PacBio HiFi sequencing and Hi-C technologies. The final genome assembly was 1.2 Gb distributed on 23 chromosomes and contained 32,365 protein coding genes with a BUSCO completeness score of 99%. A comparison of the S. tubifer genome to that of another non-luminous cardinalfish revealed a high degree of synteny, whereas a similar comparison to a more distant relative in the Gobiiformes order revealed a fusion of two chromosomes in the cardinalfish genomes. An additional comparison of orthologous clusters among these three genomes revealed a set of 710 clusters that were unique to S. tubifer in which 23 GO pathways were significantly enriched, including several relating to host-microbe interactions and one involved in visceral muscle development, which could be related to the musculature involved in the gut-associated light organ of S. tubifer. We also assembled the complete mitogenome of S. tubifer and discovered both an inversion in the WANCY tRNA gene region resulting in a WACNY gene order as well as heteroplasmy in the length of the control region for this individual. A phylogenetic analysis based on the whole mitochondrial genome indicated that S. tubifer is divergent from the rest of the cardinalfish family, bringing up questions of the involvement of the bioluminescent symbiosis in the initial divergence of the ancestral Siphamia species. This draft genome assembly of S. tubifer will enable future studies investigating the evolution of bioluminescence in fishes as well as candidate genes involved in the symbiosis and will provide novel opportunities to use this system as a vertebrate-bacteria model for symbiosis research.

category: genomics


title: Plasma derived cell-free mitochondrial DNA originates mainly from circulating cell-free mitochondria

abstract: Circulating mitochondrial DNA (cir-mtDNA) could have a potential comparable to circulating nuclear DNA (cir-nDNA), with numerous applications. However, research and development in this area have fallen behind, particularly considering its origin and structural features. To tackle this, we initially combined Q-PCR and low-pass whole genome sequencing in the same analytical strategy previously and successfully used for cir-nDNA. This revealed unexplained structural patterns and led us to correlate these data with observations made during physical examinations such as filtration, and differential centrifugation in various plasma preparations. Both the integrity index and number of reads revealed a very minor proportion of low size-ranged fragments (<1000 bp) in plasma obtained with a standard preparation (0.06%). Filtration and high speed second step centrifugation revealed that 98.7 and 99.4% corresponded to extracellular mitochondria either free or in large extracellular vesicles. When avoiding platelet activation during plasma preparation, the proportion of both types of entities was still preponderant (76-80%), but the amount of detected mitochondrial DNA decreased 67-fold. In correlation with our previous study on the presence of circulating cell-free mitochondria in blood, our differential centrifugation procedure suggested that cir-mtDNA is also associated with approximately 18% small extracellular vesicles, 1.7% exosomes and 4% protein complexes.


title: A large-volume sputum dry storage and transportation device for molecular and culture-based diagnosis of tuberculosis

abstract: Technologies for preservation of specimens in the absence of cold chains are essential for optimum utilization of existing laboratory services in the developing world. We present a prototype called specimen transportation tube (SPECTRA-tube) for the collection, exposure-free drying, ambient transportation, and liquid state recovery of large-volume (>1 mL) specimens. Specimens introduced into SPECTRA-tube are dried in glass fiber membranes, which are critical for efficient liquid-state sample recovery by rehydration and centrifugation. Mycobacterium smegmatis (Msm)-spiked mock sputum dried in native Standard 17 glass fiber was stable for molecular testing after 10-day storage at 45{degrees}C, and for culture testing after 10- and 5-day storage at 37{degrees}C and 45{degrees}C, respectively. Compatibility with human sputum storage was demonstrated by dry storing Mycobacterium bovis-spiked pooled human sputum in SPECTRA-tube for 5 days at room temperature followed by successful qPCR detection. By significantly increasing the volume of samples that can be transported in the dry state and enabling recovery of the entire sample in liquid state, SPECTRA-tube presents a potential universal solution for the preservation and transportation of liquid specimens.

category: bioengineering


title: Versatile, do-it-yourself, low-cost spinning disk confocal microscope

abstract: Confocal microscopy is an invaluable tool for 3D imaging of biological specimens, however, accessibility is often limited to core facilities due to the high cost of the hardware. We describe an inexpensive do-it-yourself (DIY) spinning disk confocal microscope (SDCM) module based on a commercially fabricated chromium photomask that can be added on to a laser-illuminated epifluorescence microscope. The SDCM achieves strong performance across a wide wavelength range (~400-800 nm) as demonstrated through a series of biological imaging applications that include conventional microscopy (immunofluorescence, small-molecule stains, and fluorescence in situ hybridization) and super-resolution microscopy (single-molecule localization microscopy and expansion microscopy). This low-cost and simple DIY SDCM is well-documented and should help increase accessibility to confocal microscopy for researchers.


title: BCN057, a Modulator of GSK3b Induces KRAS G12D Mutant Pancreatic Cancer Cell Death

abstract: Effective treatment for Pancreatic Cancer remains a major challenge due to its resistance to radiation/chemotherapy and poor drug permeability. Moreover, treatment induced normal tissue toxicity, mainly to the duodenum and gastrointestinal epithelium, is common and is a dose limiting event, while toxicity to the pancreas is relatively rare1-3. Gastrointestinal toxicity, however, often results in interruption, reduction or premature withdrawal of anti-cancer therapy which is a very significant factor impacting the overall survival of patients being treated. Therefore, development of a therapeutic strategy to selectively sensitize tumor tissue without inducing normal tissue toxicity is important. In this manuscript, we show that the novel small molecule BCN057 can modulate chemo-sensitivity of oncogenic RAS pancreatic cancer cells while conversely protecting normal intestinal epithelium from off target toxicity. In particular, BCN 057 protects Lgr5 positive intestinal stem cells, thereby preserving barrier function. Further, it is demonstrated that BCN057 inhibits GSK3{beta} and thereby induces a pro-apoptotic phosphorylation pattern on c-Jun in KRAS G12D mutant pancreatic cancer cells (Panc-1) leading to the restoration of PTEN expression and consequent apoptosis. This appears to be a new mechanistic observation for the oncogenic RAS phenotype. Lastly, concurrent with its GSK3{beta} inhibition, BCN057 is a small molecule inhibitor of PD-1 expression on human T-lymphocytes co-cultured with human pancreatic cancer cells. In summary, BCN057 can promote synthetic lethality specifically to malignant cells and therefore should be considered to improve the therapeutic ratio in pancreatic and epithelial cancer treatment in conjunction with chemotherapy and radiation.


title: BCN057, a Modulator of GSK3β Induces KRAS G12D Mutant Pancreatic Cancer Cell Death


title: Native Metabolomics Identifies the Rivulariapeptolide Family of Protease Inhibitors

abstract: The identity and biological activity of most metabolites still remain unknown. A key bottleneck in the full exploration of this tremendous source of new structures and pharmaceutical activities is the compound purification needed for bioactivity assignments of individual compounds and downstream structure elucidation. To enable bioactivity-focused compound identification from complex mixtures, we developed a scalable native metabolomics approach that integrates non-targeted liquid chromatography tandem mass spectrometry, and simultaneous detection of protein binding via native mass spectrometry. While screening for new protease inhibitors from an environmental cyanobacteria community, native metabolomics revealed 30 cyclodepsipeptides as chymotrypsin binders. Mass spectrometry-guided purification then allowed for the full structure elucidation of four new specialized metabolites via tandem mass spectrometry, chemical derivatization, and nuclear magnetic resonance spectroscopy. Together with the evaluation of biological activities, our results identified the rivulariapeptolides as a family of serine protease inhibitors with nanomolar potency, highlighting native metabolomics as promising approach for drug discovery, chemical ecology, and chemical biology studies.

category: systems biology


title: Metabolic modeling of the International Space Station microbiome reveals key microbial interactions

abstract:
BackgroundRecent studies have provided insights into the persistence and succession of microbes aboard the International Space Station (ISS), notably the dominance of Klebsiella pneumoniae. However, the interactions between the various microbes aboard the ISS, and how it shapes the microbiome remain to be clearly understood. In this study, we apply a computational approach to predict possible metabolic interactions in the ISS microbiome and shed further light on its organization.

ResultsThrough a combination of a systems-based graph-theoretical approach, and a constraint-based community metabolic modelling approach, we demonstrated several key interactions in the ISS microbiome. These complementary approaches provided insights into the metabolic interactions and dependencies present amongst various microbes in a community, highlighting key interactions and keystone species. Our results showed that the presence of K. pneumoniae is beneficial to many other microorganisms it coexists with, notably those from the Pantoea genus. Species belonging to the Enterobacteriaceae family were often found to be the most beneficial for the survival of other microorganisms in the ISS microbiome. However, K. pneumoniae was found to exhibit parasitic and amensalistic interactions with Aspergillus and Penicillium species, respectively. To prove this metabolic prediction, K. pneumoniae and Aspergillus fumigatus were co-cultured under normal and simulated microgravity, where K. pneumoniae cells showed parasitic characteristics to the fungus. The electron micrography revealed that the presence of K. pneumoniae compromised the morphology of fungal conidia and its biofilm biofilm-forming structures.

ConclusionsOur study underscores the importance of K. pneumoniae in the ISS, and its potential contribution to the survival (mutualism) and eradication (parasitism) of other microbes, including potential pathogens. This integrated modelling approach, combined with experiments, demonstrates immense potential for understanding the organization of other such microbiomes, unravelling key organisms and their interdependencies.


title: Males are more sensitive to their audience than females when scent-marking in the redfronted lemur

abstract: Audience effects, i.e. changes in behaviour caused by the presence of conspecifics, have rarely been studied in the context of olfactory communication, even though they may provide important insights into the functions of olfactory signals. Functional sex differences in scent-marking behaviours are common and influenced by the social system. To date, patterns of functional sex differences in scent-marking behaviours remain unknown in species without overt dominance relationships. We investigated sex differences in intra-group audience effects on anogenital scent-marking in a wild population of redfronted lemurs (Eulemur rufifrons) by performing focal scent-marking observations. With a combination of generalised linear mixed models and exponential random graph models, we found different audience effects in both sexes. Males were overall more sensitive than females to their audience. Only males seemed to be sensitive to the presence of both members of the opposite sex and same-sex conspecifics in the audience. Females were only moderately sensitive to the presence of other females in the audience. This study offers a potential behavioural pattern associated with anogenital scent-marking that seem to differ from those described for species exhibiting female dominance, supporting the notion that the social systems co-varies with scent-marking behaviours and scent-complexity in strepsirrhines.

category: animal behavior and cognition


title: Different B cell subpopulations show distinct patterns in their IgH repertoire metrics

abstract:
BackgroundSeveral human B-cell subpopulations are recognized in the peripheral blood, which play distinct roles in the humoral immune response. These cells undergo developmental and maturational changes involving VDJ recombination, somatic hypermutation and class switch recombination, altogether shaping their immunoglobulin heavy chain (IgH) repertoire.

MethodsHere, we sequenced the IgH repertoire of naive, marginal zone, switched and plasma cells from 10 healthy adults along with matched unsorted and in silico separated CD19+ bulk B cells. We used advanced bioinformatic analysis and machine learning to thoroughly examine and compare these repertoires.

ResultsWe show that sorted B cell subpopulations are characterised by distinct repertoire characteristics on both the individual sequence and the repertoire level. Sorted subpopulations shared similar repertoire characteristics with their corresponding in silico separated subsets. Furthermore, certain IgH repertoire characteristics correlated with the position of the constant region on the IgH locus.

ConclusionOverall, this study provides unprecedented insight over mechanisms of B cell repertoire control in peripherally circulating B cell subpopulations.


title: The AmiC/NlpD pathway dominates peptidoglycan breakdown in Neisseria meningitidisand affects cell separation, NOD1 agonist production, and infection

abstract: The human-restricted pathogen Neisseria meningitidis, which is best known for causing invasive meningococcal disease, has a nonpathogenic lifestyle as an asymptomatic colonizer of the human naso- and oropharyngeal space. N. meningitidis releases small peptidoglycan (PG) fragments during growth. It was demonstrated previously that N. meningitidis releases low levels of tripeptide PG monomer, which is an inflammatory molecule recognized by the human intracellular innate immune receptor NOD1. In this present study, we demonstrated that N. meningitidis released more PG-derived peptides compared to PG monomers. Using a reporter cell line overexpressing human NOD1, we showed that N. meningitidis activates NOD1 using PG-derived peptides. Generation of such peptides required the presence of the periplasmic N-acetylmuramyl-L-alanine amidase AmiC, and the outer membrane lipoprotein, NlpD. AmiC and NlpD were found to function in cell separation, and mutation of either amiC or nlpD resulted in large clumps of unseparated N. meningitidis cells instead of the characteristic diplococci. Using stochastic optical reconstruction microscopy, we demonstrated that FLAG epitope-tagged NlpD localized to the septum, while similarly-tagged AmiC was found at the septum in some diplococci but distributed around the cell in most cases. In a human whole blood infection assay, an nlpD mutant was severely attenuated and showed particular sensitivity to complement. Thus, in N. meningitidis the cell separation proteins AmiC and NlpD are necessary for NOD1 stimulation and for survival during infection of human blood.


title: Genomic characterisation of an mcr-1 and mcr-3-producing Escherichia coli strain isolated from pigs in France

abstract: The current study is about genomic characterisation of an atypical multidrug-resistant Escherichia coli harbouring two mobilised colistin resistance (mcr) genes isolated from pigs in France. Stool samples taken from a pig farm in Avignon in the department of the Vaucluse were subjected to a molecular screening for the detection of mcr gene variants. These samples were cultured on selective LBJMR medium. Growing bacteria were identified using MALDI-TOF, followed by antibiotic susceptibility testing (AST). Whole genome sequencing (WGS) and bioinformatic genome analysis was performed. The selective culture of stools revealed the presence of an E. coli strain named Q4552 which was simultaneously positive for the mcr-1.1 and mcr-3.5 genes. This strain exhibited resistance phenotype to fourteen antibiotics, including colistin. Genome sequencing revealed a circular chromosome and eight plasmids. Genomic analysis revealed a chromosomic integration of a mobile genetic element (MGE) harbouring the mcr-1.1 gene, while the mcr-3.5 gene was plasmidic (i.e., an IncFII plasmid). Its resistome exhibited twenty-two resistance genes, explaining its multidrug resistance phenotype. The Q4552 strain is an ST-843 clone belonging to the clonal complex Cplx-568 and is the only ST type of this cplx-568 which has been isolated from animals, humans, and the environment. Here, we report the first co-occurrence of the mcr-1 and mcr-3 genes in France from a pathogenic E. coli strain isolated from a pig farm. Since this clone (ST-843) has been reported in zoonotic transmissions, programs to monitor such colistin resistant bacterium are urgently required to avoid its spread and zoonotic transmission to humans.


title: RNASeq analysis of a Pax3-expressing myoblast clone in-vitro and effect of culture surface stiffness on differentiation.

abstract: Skeletal muscle satellite cells cultured on soft surfaces (12kPa) show improved differentiation than cells cultured on stiff surfaces (approximately 100kPa). To better understand the reasons for this, we performed an RNA-Seq analysis for a single satellite cell clone (C1F) derived from the H2kb-tsA58 immortomouse, which differentiates into myotubes under tightly regulated conditions (withdrawal of {gamma}-interferon, 37{degrees}C). As expected, the largest change in overall gene expression occurred at day 1, as cells switch from proliferation to differentiation. Surprisingly, further analysis showed that proliferating C1F cells express Pax3 and not Pax7, confirmed by immunostaining, yet their subsequent differentiation into myotubes is normal, and enhanced on softer surfaces, as evidenced by significantly higher expression levels of myogenic regulatory factors, sarcomeric genes, enhanced fusion and improved myofibrillogenesis. Levels of RNA encoding extracellular matrix structural constituents and related genes were consistently upregulated on hard surfaces, suggesting that a consequence of differentiating satellite cells on hard surfaces is that they attempt to manipulate their niche prior to differentiating. This comprehensive RNASeq dataset will be a useful resource for understanding Pax3 expressing cells.

category: cell biology


title: Adipose deficiency and aberrant autophagy in a Drosophila model of MPS VII is corrected by pharmacological stimulators of mTOR

abstract: Mucopolysaccharidosis type VII (MPS VII) is a recessively inherited lysosomal storage disorder caused due to {beta}-glucuronidase ({beta}-GUS) enzyme deficiency. Prominent clinical symptoms include hydrops fetalis, musculoskeletal deformities, neurodegeneration and hepatosplenomegaly leading to premature death in most cases. Apart from these, MPS VII is also characterized as adipose storage deficiency disorder although the underlying mechanism of this lean phenotype in the patients or {beta}-GUS-deficient mice still remains a mystery. We addressed this issue using our recently developed Drosophila model of MPS VII (the CG2135-/- fly), which also exhibited a significant loss of adiposity. We report here that the lean phenotype of the CG2135-/- fly is due to fewer number of adipocytes, smaller lipid droplets and reduced adipogenesis. Our data further revealed that there is an abnormal accumulation of autophagosomes in the CG2135-/- larvae due to autophagosome-lysosome fusion defect. Decreased lysosome-mediated turnover also led to attenuated mTOR activity in the CG2135-/- flies. Interestingly, treatment of the CG2135-/- larvae with mTOR stimulators, 3BDO or glucose, led to the restoration of mTOR activity with simultaneous correction of the autophagy defect and adipose storage deficiency. Our finding thus established a hitherto unknown mechanistic link between autophagy dysfunction, mTOR downregulation and reduced adiposity in MPS VII.


title: Allosteric mechanism of signal transduction in the two-component system histidine kinase PhoQ

abstract: Transmembrane signaling proteins couple extracytosolic sensors to cytosolic effectors. Here, we examine how binding of Mg2+ to the sensor domain of an E. coli two component histidine kinase (HK), PhoQ, modulates its cytoplasmic kinase domain. We use cysteine-crosslinking and reporter-gene assays to simultaneously and independently probe the signaling state of PhoQs sensor and autokinase domains in a set of over 30 mutants. Strikingly, conservative single-site mutants distant from the sensor or catalytic site strongly influence PhoQs ligand-sensitivity as well as the magnitude and direction of the signal, endowing diverse signaling characteristics without need for epistasis. Data from 35 mutants are explained by a semi-empirical 3-domain model in which the sensor, intervening HAMP, and catalytic domains can adopt kinase-promoting or inhibiting conformations, that are in allosteric communication. The catalytic and sensor domains intrinsically favor a constitutively  kinase-on conformation, while the HAMP favors the  off state; when coupled, they create a bistable system responsive to physiological [Mg2+]. Mutants alter signaling by locally modulating these intrinsic equilibrium constants and couplings. Our model suggests signals transmit via interdomain allostery rather than propagation of a single concerted conformational change, explaining the diversity of signaling structural transitions observed in individual HK domains.


title: Structure and ion-release mechanism of PIB-4-type ATPases

abstract: Transition metals, such as zinc, are essential micronutrients in all organisms, but also highly toxic in excessive amounts. Heavy-metal transporting P-type (PIB) ATPases are crucial for homeostasis, conferring cellular detoxification and redistribution through transport of these ions across cellular membranes. No structural information is available for the PIB-4-ATPases, the subclass with the broadest cargo scope, and hence even their topology remains elusive. Here we present structures and complementary functional analyses of an archetypal PIB-4-ATPases, sCoaT from Sulfitobacter sp. NAS14-1. The data disclose the architecture, devoid of classical so-called heavy metal binding domains, and provides fundamentally new insights into the mechanism and diversity of heavy-metal transporters. We reveal several novel P-type ATPase features, including a dual role in heavy-metal release, and as an internal counter ion, of an invariant, central histidine. We also establish that the turn-over of PIB-ATPases is potassium independent, contrasting to many other P-type ATPases. Combined with new inhibitory compounds, our results open up for efforts in e.g. drug discovery, since PIB-4-ATPases function as virulence factors in many pathogens.


title: Genetic, metabolic, and molecular insights into the diverse outcomes of diet-induced obesity in mice

abstract: Overweight and obesity are increasingly common public health issues worldwide, leading to a wide range of diseases from metabolic syndrome to steatohepatitis and cardiovascular diseases. While the increase in the prevalence of obesity is partly attributable to changes in lifestyle (i.e. increased sedentarity and changes in eating behaviour), the metabolic and clinical impacts of these obesogenic conditions varies between sexes and genetic backgrounds. The conception of personalised treatments of obesity and its complications require a thorough understanding of the diversity of responses to conditions such as high-fat diet intake. By analysing nine genetically diverse mouse strains, we show that much like humans, mice respond to high-fat diet in a genetic- and sex-dependent manner. Physiological and molecular responses to high-fat diet are associated with expression of genes involved in immunity and mitochondrial function. Finally, we find that mitochondrial function may explain part of the diversity of physiological responses. By exploring the complex interactions between genetics and metabolic phenotypes via gene expression and molecular traits, we shed light on the importance of genetic background and sex in determining metabolic outcomes. In addition to providing the community with an extensive resource for optimizing future experiments, our work serves as an exemplary design for more generalizable translational studies.


title: The nuclear effector ArPEC25 from the necrotrophic fungus Ascochyta rabiei targets the chickpea transcription factor CaβLIM1a and negatively modulates lignin biosynthesis for host susceptibility

abstract:
Fungal pathogens deploy a barrage of secreted effectors to subvert host immunity, often by evading, disrupting, or altering key components of transcription, defense signaling, and metabolic pathways. However, the underlying mechanisms of effectors and their host targets are largely unexplored in necrotrophic fungal pathogens. Here, we describe the effector protein ArPEC25, which is secreted by the necrotroph Ascochyta rabiei, the causal agent of Ascochyta blight disease in chickpea (Cicer arietinum), and is indispensable for virulence. After entering host cells, ArPEC25 localizes to the nucleus and targets the host LIM transcription factor Ca{beta}LIM1a. Ca{beta}LIM1a is a transcriptional regulator of CaPAL1, which encodes phenylalanine ammonia lyase, the regulatory, gatekeeping enzyme of the phenylpropanoid pathway. ArPEC25 inhibits the transactivation of Ca{beta}LIM1a by interfering with its DNA binding ability. This results in negative regulation of the phenylpropanoid pathway and decreased levels of intermediates of lignin biosynthesis, thereby suppressing lignin production. Our findings illustrate the role of fungal effectors in enhancing virulence by targeting a key defense pathway that leads to the biosynthesis of various secondary metabolites and antifungal compounds. This study provides a template for the study of less explored necrotrophic effectors and their host target functions.

One-sentence summaryThe Ascochyta rabiei effector ArPEC25 enters the host nucleus and targets the transcription factor Ca{beta}LIM1a to manipulate phenylpropanoid pathway for negative modulation of chickpea lignin biosynthesis.


title: Robust de novo design of protein binding proteins from target structural information alone

abstract: The design of proteins that bind to a specific site on the surface of a target protein using no information other than the three-dimensional structure of the target remains an outstanding challenge. We describe a general solution to this problem which starts with a broad exploration of the very large space of possible binding modes and interactions, and then intensifies the search in the most promising regions. We demonstrate its very broad applicability by de novo design of binding proteins to 12 diverse protein targets with very different shapes and surface properties. Biophysical characterization shows that the binders, which are all smaller than 65 amino acids, are hyperstable and bind their targets with nanomolar to picomolar affinities. We succeeded in solving crystal structures of four of the binder-target complexes, and all four are very close to the corresponding computational design models. Experimental data on nearly half a million computational designs and hundreds of thousands of point mutants provide detailed feedback on the strengths and limitations of the method and of our current understanding of protein-protein interactions, and should guide improvement of both. Our approach now enables targeted design of binders to sites of interest on a wide variety of proteins for therapeutic and diagnostic applications.


title: Control of protein synthesis and memory by GluN3A-NMDA receptors through inhibition of GIT1/mTORC1 assembly

abstract: De novo protein synthesis is required for synapse modifications underlying stable memory encoding. Yet neurons are highly compartmentalized cells and how protein synthesis can be regulated at the synapse level is unknown. Here we characterize neuronal signaling complexes formed by the postsynaptic scaffold GIT1, the mTOR kinase and Raptor that couple synaptic stimuli to mTOR-dependent protein synthesis; and identify NMDA receptors containing GluN3A subunits as key negative regulators of GIT1 binding to mTOR. Disruption of GIT1/mTOR complexes by enhancing GluN3A expression or silencing GIT1 inhibits synaptic mTOR activation and restricts the mTOR-dependent translation of specific activity-regulated mRNAs. Conversely, GluN3A removal enables complex formation, potentiates mTOR-dependent protein synthesis, and facilitates the consolidation of associative and spatial memories in mice. The memory enhancement becomes evident with light or spaced training, can be achieved by selectively deleting GluN3A from excitatory neurons during adulthood, and does not compromise other aspects of cognition such as memory flexibility or extinction. Our findings provide mechanistic insight into synaptic translational control and reveal a potentially selective target for cognitive enhancement.

category: neuroscience


title: Generation of CD34+CD43+ hematopoietic progenitors to induce thymocytes from human pluripotent stem cells

abstract: Immunotherapy using primary T cells has revolutionized medical care in some pathologies in recent years but limitations associated to challenging cell genome edition, insufficient cell number production, the use of only autologous cells and lack of product standardization have limited its uses in the clinic. The alternative use of T cells generated in vitro from human pluripotent stem cells (hPSCs) offers great advantages by providing a self-renewing source of T cells that can be readily genetically modified and facilitate the use of standardized universal off-the-shelf allogeneic cell products and rapid clinic access. However, despite their potential, the feasibility and functionality of T-cells differentiated from hPSCs needs better comprehension before moving to the clinic. In this study, we generated human induced pluripotent stem cells from T-cells (T-iPSCs) allowing preservation of already recombined TCR, with the same properties as human embryonic stem cells (hESCs). Based on these cells, we differentiated with high efficiency hematopoietic progenitor stem cells (HPSCs), capable of self-renewal and differentiation into any cell blood type, and then DN3a thymic progenitors from several T-iPSC lines. To better comprehend differentiation, we analyzed the transcriptomic profiles of the different cell types and demonstrated that HPSCs differentiated from hiPSCs had a very similar profile to cord blood hematopoietic stem cells (HSCs). Furthermore, differentiated T-cell progenitors had a similar profile to thymocytes at the DN3a stage of thymic lymphopoiesis. Therefore, with this approach, we were able to regenerate precursors of therapeutic human T cells to potentially treat a wide number of diseases.

category: developmental biology


title: Intronic elements associated with insomnia and restless legs syndrome exhibit cell type-specific epigenetic features contributing to MEIS1 regulation

abstract: A highly evolutionarily conserved MEIS1 intronic region is strongly associated with restless legs syndrome (RLS) and insomnia. To understand its regulatory function, we dissected the region by analyzing chromatin accessibility, enhancer-promoter contacts, DNA methylation, and eQTLs in different human neural cell types and tissues. We observed specific activity with respect to cell type and developmental maturation, indicating a prominent role for distinct highly conserved intronic elements in forebrain inhibitory neuron differentiation. Two elements were hypomethylated in neural cells with higher MEIS1 expression, suggesting a role of enhancer demethylation in gene regulation. MEIS1 eQTLs showed a striking modular chromosomal distribution, with forebrain eQTLs clustering in intron 8/9. CRISPR interference targeting of individual elements in this region attenuated MEIS1 expression, revealing a complex regulatory interplay of distinct elements. In summary, we found that MEIS1 regulation is organized in a modular pattern. Disease-associated intronic regulatory elements control MEIS1 expression with cell type and maturation stage specificity, particularly in the inhibitory neuron lineage. The precise spatiotemporal activity of these elements likely contributes to the pathogenesis of insomnia and RLS.


title: TTYH family members form tetrameric complexes within the cell membrane

abstract: The conserved Tweety homolog (TTYH) family consists of three paralogs in vertebrates, displaying a ubiquitous expression pattern. Although considered as ion channels for almost two decades, recent structural and functional analyses refuted this role. Intriguingly, while all paralogs, studied following detergent solubilization, shared a dimeric stoichiometry, their spatial organization differed. Here, we determined the stoichiometry of intact mouse TTYH (mTTYH) complexes in cells. Using cross-linking and single-molecule fluorescence microscopy, we demonstrated that mTTYH1 and mTTYH3 form tetramers at the plasma membrane. Blue-native PAGE and fluorescence-detection size-exclusion chromatography analyses revealed that detergent solubilization results in the dissolution of tetramers into dimers, suggesting a dimer-of-dimers assembly mode. As cross-linking analysis of the soluble extracellular domains also showed tetrameric stoichiometry, we explored the effect of membrane solubilization and disulfide bridges integrity and established their contribution to tetramer stability. Future studies of the native tetrameric TTYH characterized here may illuminate their long-sought cellular function.


title: Scallop2 enables accurate assembly of multiple-end RNA-seq data

abstract: Transcript assembly (i.e., to reconstruct the full-length expressed transcripts from RNA-seq data) has been a critical but yet unsolved step in RNA-seq analysis. Modern RNA-seq protocols can produce paired-/multiple-end RNA-seq reads, where information is available that two or more reads originate from the same transcript. The long-range constraints implied in these paired-/multiple-end reads can be much beneficial in correctly phasing the complicated spliced isoforms. However, there often exist gaps among individual ends, which may even contain junctions, making the efficient use of such constraints algorithmically challenging. Here we introduce Scallop2, a new reference-based transcript assembler optimized for multiple-end (including paired-end) RNA-seq data. Scallop2 uses an algorithmic frame-work that first represents reads from the same molecule as the so-called multiple-end phasing paths in the context of a splice graph, then "bridges" each multiple-end phasing path into a long, single-end phasing path, and finally decomposes the splice graph into paths (i.e., transcripts) guided by the bridged phasing paths. An efficient bridging algorithm is designed to infer the true path connecting two consecutive ends following a novel formulation that is robust to sequencing errors and transcript noises. By observing that failing to bridge two ends is mainly due to incomplete splice graphs, we propose a new method to determine false starting/ending vertices of the splice graphs which has been showed efficient in reducing false positive rate. Evaluations on both (multiple-end) single-cell RNA-seq datasets from Smart-seq3 protocol and Illumina paired-end RNA-seq samples demonstrate that Scallop2 vastly outperforms recent assemblers including StringTie2, Scallop, and CLASS2 in assembly accuracy.


title: Intersectional, anterograde transsynaptic targeting of the neurons receiving monosynaptic inputs from two upstream regions

abstract: A brain region typically receives inputs from multiple upstream areas. However, currently, no method is available to selectively access neurons that receive monosynaptic inputs from two upstream regions. Here, we devised a method to genetically label such neurons in mice by combining the anterograde transsynaptic spread of adeno-associated virus serotype 1 (AAV1) with intersectional gene expression. Injections of AAV1s expressing either Cre or Flpo recombinases and the Cre/Flpo double-dependent AAV into two upstream regions and the downstream region, respectively, were used to label the postsynaptic neurons receiving inputs from the two upstream regions. We demonstrated this labelling in two distinct circuits: the retina/primary visual cortex to the superior colliculus and the bilateral motor cortex to the dorsal striatum. Systemic delivery of the intersectional AAV allowed for unbiased detection of the labelled neurons throughout the brain. This strategy may help analyse the interregional integration of information in the brain.


title: ISL1 is an essential determinant of structural and functional tonotopic representation of sound

abstract: A cardinal feature of the auditory pathway is frequency selectivity, represented in a tonotopic map from the cochlea to the cortex. The molecular determinants of the auditory frequency map are unknown. Here, we discovered that the transcription factor ISL1 regulates the molecular and cellular features of auditory neurons, including the formation of the spiral ganglion and peripheral and central processes that shape the tonotopic representation of the auditory map. We selectively knocked out Isl1 in auditory neurons using Neurod1Cre strategies. In the absence of Isl1, spiral ganglion neurons migrate into the central cochlea and beyond, and the cochlear wiring is profoundly reduced and disrupted. The central axons of Isl1 mutants lose their topographic projections and segregation at the cochlear nucleus. Transcriptome analysis of spiral ganglion neurons shows that Isl1 regulates neurogenesis, axonogenesis, migration, neurotransmission-related machinery, and synaptic communication patterns. We show that peripheral disorganization in the cochlea affects the physiological properties of hearing in the midbrain and auditory behavior. Surprisingly, auditory processing features are preserved despite the significant hearing impairment, revealing central auditory pathway resilience and plasticity in Isl1 mutant mice. Mutant mice have a reduced acoustic startle reflex, altered prepulse inhibition, and characteristics of compensatory neural hyperactivity centrally. Our findings show that ISL1 is one of the obligatory factors required to sculpt auditory structural and functional tonotopic maps. Still, upon Isl1 deletion, the ensuing central compensatory plasticity of the auditory pathway does not suffice to overcome developmentally induced peripheral dysfunction of the cochlea.


title: ISL1 is necessary for auditory neuron development and tonotopic organization


title: ISL1 is necessary for auditory neuron development and contributes towards tonotopic organization


title: Neural circuit dynamics of drug-context associative learning in the hippocampus

abstract: The environmental context associated with previous drug consumption serves as a potent trigger for relapse to drug use. The mechanism by which existing neural representations of context are modified to incorporate information associated with a given drug however, remains unknown. Using longitudinal calcium imaging in freely behaving mice, we reveal that drug-context associations for psychostimulants and opioids are encoded in a subset of hippocampal neurons. In these neurons, drug context pairing in a conditioned place preference task weakened their spatial coding for the nondrug-paired context, with drug-induced changes to spatial coding predictive of drug-seeking behavior. Furthermore, the dissociative drug ketamine blocked both the drug-induced changes to hippocampal coding and corresponding drug-seeking behavior. Together, this work reveals how drugs of abuse can alter the hippocampal circuit to encode drug-context associations and points to the hippocampus as a key node in the cognitive process of drug addiction and context-induced drug relapse.


title: Over-expression Screen of Interferon-Stimulated Genes Identifies RARRES3 as a Restrictor of Toxoplasma gondii Infection

abstract: Toxoplasma gondii is an important human pathogen infecting an estimated 1 in 3 people worldwide. The cytokine interferon gamma (IFN{gamma}) is induced during infection and is critical for restricting T. gondii growth in human cells. Growth restriction is presumed to be due to the induction interferon stimulated genes (ISGs) that are upregulated to protect the host from infection. Although there are hundreds of ISGs induced by IFN{gamma}, their individual roles in restricting parasite growth in human cells remain somewhat elusive. To address this deficiency, we screened a library of 414 IFN{gamma} induced ISGs to identify factors that impact T. gondii infection in human cells. In addition to IRF1, which likely acts through induction of numerous downstream genes, we identified RARRES3 as a single factor that restricts T. gondii infection by inducing premature egress of the parasite in multiple human cell lines. Overall, while we successfully identified a novel IFN{gamma} induced factor restricting T. gondii infection, the limited number of ISGs capable of restricting T. gondii infection when individually expressed suggests that IFN{gamma} mediated immunity to T. gondii infection is a complex, multifactorial process.

category: microbiology


title: Hippo pathway and Bonus control developmental cell fate decisions in the Drosophila eye

abstract: The Hippo pathway controls organ growth, however its role in cell fate determination and the underlying mechanism is not well understood. Here, we uncover the function of the Hippo pathway in developmental cell fate decisions in the Drosophila eye, exerted through the interaction of Yorkie (Yki) with a transcriptional regulator Bonus (Bon). Activation of either Bon or Yki is sufficient to promote epidermal fate at the expense of eye fate through the recruitment of multiple transcriptional and post-transcriptional regulators. Transcriptome analysis reveals that Bon and Yki jointly upregulate epidermal differentiation genes and downregulate Notch target genes that modulate the eye-to-epidermal fate switch. The Hippo pathway and Bon also control the early eye-antennal specification, with activated Yki and Bon suppressing eye fate and promoting antennal fate. Our work has revealed that the Hippo pathway and Bon control cell fate decisions during Drosophila eye development at multiple levels.


title: Prevalence of mutualism in a simple model of microbial co-evolution

abstract: Evolutionary transitions among ecological interactions are widely known, although their detailed dynamics remain absent for most population models. Adaptive dynamics has been used to illustrate how the parameters of population models might shift through evolution, but within an ecological regime. Here we use adaptive dynamics combined with a generalised logistic model of population dynamics to show that transitions of ecological interactions might appear as a consequence of evolution. To this purpose we introduce a two-microbial toy model in which population parameters are determined by a bookkeeping of resources taken from (and excreted to) the environment, as well as from the byproducts of the other species. Despite its simplicity, this model exhibits all sorts of ecological transitions, some of which resemble those found in nature. Overall, the model shows a clear trend toward the emergence of mutualism.

category: ecology