The dynamic range of voltage-dependent gap junction signaling is maintained by Ih-induced membrane potential depolarization Like their chemical counterparts, electrical synapses show complex dynamics such as rectification and voltage dependence that interact with other electrical processes in neurons. The consequences arising from these interactions for the electrical behavior of the synapse, and the dynamics they create, remain largely unexplored. Using a voltage-dependent electrical synapse between a descending modulatory projection neuron (MCN1) and a motor neuron (LG) in the crustacean stomatogastric ganglion, we find that the influence of the hyperpolarization-activated inward current (Ih) is critical to the function of the electrical synapse. When we blocked Ih with CsCl, the apparent voltage dependence of the electrical synapse shifted by 18.7 mV to more hyperpolarized voltages, placing the dynamic range of the electrical synapse outside of the range of voltages used by the LG motor neuron (-60.2 mV - -44.9 mV). With dual electrode current- and voltage-clamp recordings, we demonstrate that this voltage shift is not due to a change in the properties of the gap junction itself, but is a result of a sustained effect of Ih on the presynaptic MCN1 axon terminal membrane potential. Ih-induced depolarization of the axon terminal membrane potential increased the electrical postsynaptic potentials and currents. With Ih present, the axon terminal resting membrane potential depolarized, shifting the dynamic range of the electrical synapse towards the functional range of the motor neuron. We thus demonstrate that the function of an electrical synapse is critically influenced by a voltage-dependent ionic current (Ih).

New & NoteworthyElectrical synapses and voltage-gated ionic currents are often studied independently from one another, despite mounting evidence that their interactions can alter synaptic behavior. We show that the hyperpolarization-activated inward ionic current shifts the voltage dependence of an electrical synaptic transmission through its depolarizing effect on the membrane potential, enabling it to lie within the functional membrane potential range of a motor neuron. Thus, the electrical synapses function critically depends on the voltage-gated ionic current.
Thalamic control of sensory enhancement and sleep spindle properties in a biophysical model of thalamoreticular microcircuitry Thalamoreticular circuitry is known to play a key role in attention, cognition and the generation of sleep spindles, and is implicated in numerous brain disorders, but the cellular and synaptic mechanisms remain intractable. Therefore, we developed the first detailed computational model of mouse thalamus and thalamic reticular nucleus microcircuitry that captures morphological and biophysical properties of [~]14,000 neurons connected via [~]6M synapses, and recreates biological synaptic and gap junction connectivity. Simulations recapitulate multiple independent network-level experimental findings across different brain states, providing a novel unifying cellular and synaptic account of spontaneous and evoked activity in both wakefulness and sleep. Furthermore, we found that: 1.) inhibitory rebound produces frequency-selective enhancement of thalamic responses during wakefulness, in addition to its role in spindle generation; 2.) thalamic interactions generate the characteristic waxing and waning of spindle oscillations; and 3.) changes in thalamic excitability (e.g. due to neuromodulation) control spindle frequency and occurrence. The model is openly available and provides a new tool to interpret spindle oscillations and test hypotheses of thalamoreticular circuit function and dysfunction across different network states in health and disease.
Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation, however this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males, 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci, however these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16 % of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p-value = 4.6e-13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8,420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

SignificanceThe importance of uncovering biological sex differences and their translation to physiology has become increasingly evident. Using a large-scale meta-analysis of three cohorts, we perform the first comparison of genome-wide skeletal muscle DNA methylation between males and females, and identify thousands of genes that display sex-differential methylation. We then explore intrinsic biological factors that may be underlying the DNA methylation sex differences, such as fibre type proportions and sex hormones. Leveraging the GTEx database, we identify hundreds of genes with both sex-differential expression and DNA methylation in skeletal muscle. We further confirm the sex-biased genes with gene expression data from two cohorts included in the methylation meta-analysis. Our study integrates genomewide sex-biased DNA methylation and expression in skeletal muscle, shedding light on distinct sex differences in skeletal muscle. An NLR Integrated Domain toolkit to identify plant pathogen effector targets Plant resistance genes (or NLR "Nod-like Receptors") are known to contain atypical domains procuring them with a decoy capacity. Some of these integrated domains (or ID) allow the plant to lure the virulence determinants ("effectors") of pathogens and triggering a specific NLR immune reaction.

In this work, our goal was to generate a library of known IDs that could be screened with plant pathogen effectors in order to identify putative new effector virulence targets and NLR-effector pairs.

We curated the IDs contained in NLRs from seven model and crop plant species. We cloned 52 IDs representing 31 distinct Pfam domains. This library was screened for interaction by yeast-two-hybrid with a set of 31 conserved Ralstonia solanacearum type III effectors. This screening and the further in planta interaction assay allowed us to identify three interactions, involving different IDs (kinase, DUF3542, WRKY) and two type III effectors (RipAE and PopP2).

PopP2 was found to physically interact with ID#85, an atypical WRKY domain integrated in the GmNLR-ID85 NLR protein from Soybean. Using a imaging method in living plant cells, we showed that PopP2 associates with ID#85 in the nucleus. But unlike the known WRKY-containing Arabidopsis RRS1-R NLR receptor, this newly identified soybean WRKY domain could not be acetylated by PopP2 and its atypical sequence (WRKYGKR) also probably renders it inefficient in plant immunity triggering.

This ID toolkit is available for screening with other plant pathogen effectors and should prove useful to discover new effectors targets and potentially engineer new plant resistance genes. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Interactions and information: Exploring task allocation in ant colonies using network analysis In animal societies, individuals may take on different roles to fulfil their own needs and the needs of their groups. Ant colonies display high levels of organisational complexity, with ants fulfilling different roles at different timescales (what is known as task allocation). Factors affecting task allocation can be at the individual level (e.g., physiology), or at the group level (e.g., the network of interactions). We focus on group level processes by exploring the relationship between interaction networks, task allocation and task switching using a previously published dataset (Mersch et al., 2013) tracking the behaviour of six Camponotus fellah colonies over 41 days. In our new analyses, our goal was to better explain the noisy process of task switching beyond simple age polyethism. First, we investigated the architecture of interaction networks using node (individual) level network measures and their relation to the individuals task - foraging, cleaning or nursing - and whether or not the ant switched tasks. We then explored how noisy information propagation was among ants, as a function of the colony composition (how many ants carried out which tasks), through the information-theoretic metric of Effective Information. Our results show that interaction history was tied to task allocation: ants who switched to a task are more likely to have interacted with other ants carrying out that task. The degree to which interactions related to task allocation, as well as the noise in those interactions, depended on which groups of ants were interacting. Overall, we showed that colony cohesion was stable even as ant-level network measures varied more for ants when they switched functional groups; thus, ant colonies maintained a high level of information flow as determined by network analysis, and ant functional groups played different roles in maintaining colony cohesion through varied information flows.

HighlightsO_LIWe analysed the interaction networks of six Camponotus fellah colonies
C_LIO_LIWe tested how centrality and information flow were tied to task switching
C_LIO_LINode-level network metrics and the information theoretic measure of Effective Information explained differences among functional groups
C_LIO_LIInteractions were correlated with task switching, but the strength of the correlation differed across functional groups
C_LI Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Protection induced by a human monoclonal antibody recognizing two different epitopes in a conserved region of streptococcal M proteins Group A streptococci have evolved multiple strategies to evade human antibodies, making it challenging to create effective vaccines or antibody treatments. Here, we have generated antibodies derived from the memory B cells of an individual who had successfully cleared a group A streptococcal infection. The antibodies bind with high affinity in the central region of the surface-bound M protein. Such antibodies are typically non-opsonic. However, one antibody could effectively promote vital immune functions, including phagocytosis and in vivo protection. Remarkably, this antibody primarily interacts through a bivalent dual-Fab cis mode, where the Fabs bind to two distinct epitopes in the M protein. The dual-Fab cis binding phenomenon is conserved across different groups of M types. In contrast, other antibodies binding with normal single-Fab mode to the same region can not bypass the M proteins virulent effects. A broadly binding, protective monoclonal antibody could be a candidate for anti-streptococcal therapy. Our findings highlight the concept of dual-Fab cis binding as a means to access conserved, and normally non-opsonic regions, for protective antibody targeting. Missing the sweet spot: one of the two N-glycans on human Gb3/CD77 synthase is expendable N-glycosylation is a ubiquitous posttranslational modification that may influence folding, subcellular localization, secretion, solubility and oligomerization of proteins. In this study, we examined the effects of N-glycans on the activity of human Gb3/CD77 synthase, which catalyzes the synthesis of glycosphingolipids with terminal Gal1[-&gt;]4Gal (Gb3 and the P1 antigen) and Gal1[-&gt;]4GalNAc disaccharides (the NOR antigen). The human Gb3/CD77 synthase contains two occupied N-glycosylation sites at positions N121 and N203. Intriguingly, we found that while the N-glycan at N203 is essential for activity and correct subcellular localization, the N-glycan at N121 is dispensable and its absence did not reduce, but, surprisingly, even increased the activity of the enzyme. The fully N-glycosylated human Gb3/CD77 synthase and its glycoform missing the N121 glycan correctly localized in the Golgi, whereas a glycoform without the N203 site partially mislocalized in the endoplasmic reticulum. A double mutein missing both N-glycans was inactive and accumulated in the endoplasmic reticulum. Our results suggest that the decreased specific activity of human Gb3/CD77 synthase glycovariants results from their improper subcellular localization and, to a smaller degree, a decrease in enzyme solubility. Taken together, our findings show that the two N-glycans of human Gb3/CD77 synthase have opposing effects on its properties, revealing a dual nature of N-glycosylation and potentially a novel regulatory mechanism controlling the biological activity of proteins. High performance single-cell gene regulatory network inference at scale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost Efficient and Economical Targeted Insertion in Plant Genomes via Protoplast Regeneration Versatile genome editing can be facilitated by the insertion of DNA sequences into specific locations. Current protocols involving clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins rely on low efficiency homology-directed repair or non-homologous end joining with modified double-stranded DNA oligonucleotides as donors. Our simple protocol eliminates the need for expensive equipment, chemical and enzymatic donor DNA modification or plasmid construction by using polyethylene glycol-calcium to deliver non-modified single-stranded DNA oligonucleotides and CRISPR-Cas9 ribonucleoprotein into protoplasts. Plants regenerated via edited protoplasts achieved targeted insertion frequencies of up to 50.0% in Nicotiana benthamiana and 13.6% in rapid cycling Brassica oleracea without antibiotic selection. Using a 60-nt donor containing 27 nt in each homologous arm, 6 of 22 regenerated N. benthamiana plants showed targeted insertions, and one contained a precise insertion of a 6-bp HindIII site. The inserted sequences were transmitted to the next generation and invite the possibility of future exploration of versatile genome editing by targeted DNA insertion in plants. Multiparametric profiling of cellular lipid metabolism in hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

MotivationWe have limited information on how cellular lipid uptake and processing differ between individuals and influence the development of metabolic diseases, such as hypercholesterolemia. Available assays are labor intensive, require skilled personnel and are difficult to scale to higher throughput, making it challenging to obtain systematic functional cell-based data from individuals. To overcome this problem, we established a scalable automated analysis pipeline enabling reliable quantification of multiple cellular readouts, including lipid uptake, storage and mobilization, from different white blood cell populations. This approach provides new personalized insights into the cellular basis of hypercholesterolemia and obesity.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development. Hot and cold spots of pest-induced US urban tree death, 2020-2050 O_LIUrban trees are important nature-based solutions for future wellbeing and livability but are at high risk of mortality from insect pests. In the United States (US), 82% of the population live in urban settings and this number is growing, making urban tree mortality a matter of concern for most of its population. Until now, the magnitudes and spatial distributions of risks were unknown.
C_LIO_LIHere, we combine new models of street tree populations in [~]30,000 US communities, species-specific spread predictions for 57 invasive insect species, and estimates of tree death due to insect exposure for 48 host tree genera.
C_LIO_LIWe estimate that 1.4 million street trees will be killed by invasive insects from 2020 through 2050, costing an annualized average of US$ 30M. However, these estimates hide substantial variation: 23% of urban centers will experience 95% of all insect-induced mortality. Further, 90% of all mortality will be due to emerald ash borer (Agrilus planipennis, EAB), which is expected to kill virtually all ash trees (Fraxinus spp.) in >6000 communities.
C_LIO_LIWe define an EAB high-impact zone spanning 902,500km2, largely within the southern and central US, within which we predict the death of 98.8% of all ash trees. "Mortality hotspot cities" include Milwaukee, WI; Chicago, IL; and New York, NY.
C_LIO_LIWe identify Asian wood borers of maple and oak trees as the highest risk future invaders, where a new establishment could cost US$ 4.9B over 30 years.
C_LIO_LIPolicy implications: To plan effective mitigation, managers must know which tree species in which communities will be at the greatest risk, as well as the highest-risk insects. We provide the first country-wide, spatial forecast of urban tree mortality due to invasive insect pests. This framework identifies dominant pest insects and spatial impact hotspots, which can provide the basis for spatial prioritization of spread control efforts such as quarantines and biological control release sites. Our results highlight the need for EAB early-detection efforts as far from current infestations as Seattle, WA. Further, these findings produce a list of biotic and spatiotemporal risk factors for future high-impact US urban forest insect pests.
C_LI Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

METHODSWe combined mice models with kidney biopsies and renal epithelial cells sampled from human urine to characterize the renal inflammatory network of nephronophthisis (NPH), the most frequent renal ciliopathy in children.

RESULTSIn human, mutations in cilia genes involved in NPH enhance urine excretion of the chemokine CCL2, causing abnormal macrophage recruitment in kidney tissues from NPH patients. Differing from ADPKD, inactivating Ccl2 specifically in mouse tubular cells does not rescue the NPH phenotype, suggesting that other inflammatory mediators are involved. Using transcriptional data from 2 NPH models, we identify a set of pro-inflammatory cytokines upregulated in this disease, independently of CCL2. The majority of detectable transcripts from this set are specifically upregulated in kidney cells from NPH patients. In line with the function of these cytokines, NPH kidneys show disproportionate neutrophils and T cells infiltrates compared to healthy subject or hypertensive and diabetic chronic kidney disease patients.

CONCLUSIONSThis study reveals that inflammation is a central aspect in human NPH and delineates a specific set of inflammatory mediators that regulates immune cell recruitment in human NPH.

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. Alteration of gut microbiota with a broad-spectrum antibiotic does not impair maternal care in the European earwig. AO_SCPLOWBSTRACTC_SCPLOWThe microbes residing within the gut of an animal host often increase their own fitness by modifying their hosts physiological, reproductive, and behavioural functions. Whereas recent studies suggest that they may also shape host sociality and therefore have critical effects on animal social evolution, the impact of the gut microbiota on maternal care remains unexplored. This is surprising, as this behaviour is widespread among animals, often determines the fitness of both juveniles and parents, and is essential in the evolution of complex animal societies. Here, we tested whether life-long alterations of the gut microbiota with rifampicin - a broad-spectrum antibiotic - impair pre- and post-hatching maternal care in the European earwig. Our results first confirm that rifampicin altered the mothers gut microbial communities and indicate that the composition of the gut microbiota differs before and after egg care. Contrary to our predictions, however, the rifampicin-induced alterations of the gut microbiota did not modify pre- or post-hatching care. Independent of maternal care, rifampicin increased the females feces production and resulted in lighter eggs and juveniles. By contrast, rifampicin altered none of the other 21 physiological, reproductive and longevity traits measured over the 300 days of a females lifetime. Overall, these findings reveal that altering the gut microbiota with a large spectrum antibiotic such as rifampicin does not necessarily affect host sociality. They also emphasize that not all animals have evolved a co-dependence with their microbiota and call for caution when generalizing the central role of gut microbes in host biology. Modeling the SARS-CoV-2 nsp1--5'-UTR complex via extended ensemble simulations Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2s own RNA evades nsp1s host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results showed that the stem loop structure of the SARS-CoV-2 RNA 5-untranslated region (SL1) is recognized by both nsp1s globular region and intrinsically disordered region. The recognition presumably enables selective translation of viral RNAs. Cluster analysis of the binding mode and detailed analysis of the binding poses revealed several residues involved in the SL1 recognition mechanism. The simulation results imply that the nsp1 C-terminal helices are lifted from the 40S ribosome upon the binding of SL1 to nsp1, unblocking translation of the viral RNA. Contrasting resistance and resilience to light variation of the coupled oxic and anoxic componentsof an experimental microbial ecosystem Understanding how microbial communities of aquatic ecosystems respond to environmental change remains a critical challenge in microbial ecology. In this study, we used phototrophic oxic-anoxic micro-ecosystems to understand how the functioning and diversity of aerobic and anaerobic lake analog communities is affected by light deprivation. Continuous measurements were performed to describe oxygen dynamics (mean/min/max/amplitude) and time-series of full-length 16S rRNA sequencing were used to quantify changes in alpha- and beta-diversity. In the top oxic layer, oxygen concentration decreased significantly under light deprivation, but showed resilience in mean, minimum and maximum after light conditions were restored. Only the amplitude of diurnal fluctuations in oxygen concentrations did not recover fully, and instead tended to remain lower in treated ecosystems. Alpha-diversity of the top oxic layer communities showed a delayed increase after light conditions were restored, and was not resilient. In contrast, alpha-diversity of the anoxic bottom layer communities increased due to the stressor, but was resilient. Community composition changed significantly during light deprivation, showed resilience in the oxic layer and lack of resilience in the anoxic layer. Alpha-diversity and the amplitude of oxygen within and among treatments were strongly correlated, suggesting that higher diversity could lead to less variable oxygen concentrations, or vice versa. Our experiment showed that light deprivation induces multifaceted responses of community function and structure, hence focusing on a single stability component could potentially be misleading. MMMVI: Detecting SARS-CoV-2 Variants of Concern in Metagenomic Wastewater Samples MotivationSARS-CoV-2 is the causative agent of the COVID-19 pandemic. Variants of Concern (VOCs) and Variants of Interest (VOIs) are lineages that represent a greater risk to public health, and can be differentiated from the wildtype virus based on unique profiles of signature mutations. Metagenomic sequence analysis of wastewater represents an emerging form of surveillance that can capture early signal for these variants in a community prior to detection through public health testing or genomic surveillance activities. However, because multiple viral genomes are likely to be present in a metagenomic sample, additional analytical scrutiny of the sequencing reads beyond variant calling is needed to increase confidence in diagnostic determinations.

ResultsWhere multiple signature mutations are present on a given read, they can be used as enhanced biomarkers to confirm the presence of a VOC/VOI in the sample. We have developed MMMVI, a tool to aggregate and report on the likely presence of a VOC/VOI in a set of metagenomic reads based on the detection of reads bearing multiple signature mutations.

AvailabilityMMMVI is implemented in Python, and is available under the MIT licence from https://github.com/dorbarker/voc-identify/

Contactdillon.barker@canada.ca Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. Single-cell Transcriptomic Profiling Unveils Cardiac Cell-type Specific Response to Maternal Hyperglycemia Underlying the Risk of Congenital Heart Defects Congenital heart disease (CHD) is the most prevalent structural malformations of the heart affecting [~]1% of live births. To date, both damaging genetic variations and adverse environmental exposure such as maternal diabetes have been found to cause CHD. Clinical studies show [~]fivefold higher risk of CHD in the offspring of mothers with pregestational diabetes. Maternal pregestational diabetes affects the gene regulatory networks key to proper cardiac development in the fetus. However, the cell-type specificity of these gene regulatory responses to maternal diabetes and their association with the observed cardiac defects in the fetuses remains unknown. To uncover the transcriptional responses to maternal diabetes in the early embryonic heart, we used an established murine model of pregestational diabetes. In this model, we have previously demonstrated an increased incidence of CHD. Here, we show maternal hyperglycemia (matHG) elicits diverse cellular responses during heart development by single-cell RNA-sequencing in embryonic hearts exposed to control and matHG environment. Through differential gene-expression and pseudotime trajectory analyses of this data, we identified changes in lineage specifying transcription factors, predominantly affecting Isl1+ second heart field progenitors and Tnnt2+ cardiomyocytes with matHG. Using in vivo cell-lineage tracing studies, we confirmed that matHG exposure leads to impaired second heart field-derived cardiomyocyte differentiation. Finally, this work identifies matHG-mediated transcriptional determinants in cardiac cell lineages elevate CHD risk and show perturbations in Isl1-dependent gene-regulatory network (Isl1-GRN) affect cardiomyocyte differentiation. Functional analysis of this GRN in cardiac progenitor cells will provide further mechanistic insights into matHG-induced severity of CHD associated with diabetic pregnancies.