Synaptic targets of photoreceptors specialized to detect color and skylight polarization in Drosophila Color and polarization provide complementary information about the world and are detected by specialized photoreceptors. However, the downstream neural circuits that process these distinct modalities are incompletely understood in any animal. Using electron microscopy, we have systematically reconstructed the synaptic targets of the photoreceptors specialized to detect color and skylight polarization in Drosophila, and we have used light microscopy to confirm many of our findings. We identified known and novel downstream targets that are selective for different wavelengths or polarized light, and followed their projections to other areas in the optic lobes and the central brain. Our results revealed many synapses along the photoreceptor axons between brain regions, new pathways in the optic lobes, and spatially segregated projections to central brain regions. Strikingly, photoreceptors in the polarization-sensitive dorsal rim area target fewer cell types, and lack strong connections to the lobula, a neuropil involved in color processing. Our reconstruction identifies shared wiring and modality-specific specializations for color and polarization vision, and provides a comprehensive view of the first steps of the pathways processing color and polarized light inputs.
Meningeal dendritic cells-astrocytes interactions elevate the kynurenine metabolic pathway to sustain neuropathic pain Neuropathic pain is triggered by injury to the somatosensory system, and is one of the most important types of chronic pain. Nevertheless, critical pathophysiological mechanisms that maintain neuropathic pain are poorly understood. Here, we show that neuropathic pain is abrogated when the kynurenine metabolic pathway (KYNPATH) initiated by the enzyme indoleamine 2,3-dioxygenase (IDO) is ablated pharmacologically or genetically. Mechanistically, it was found that IDO upregulation in dendritic cells that accumulate in the dorsal root leptomeninges led to increased levels of kynurenine (Kyn) in the spinal cord, where Kyn is metabolized by astrocytes-expressed kynurenine-3-monooxygenase into a pro-nociceptive metabolite 3-hydroxykynurenine. In conclusion, these data reveal a novel role for KYNPATH as an important factor maintaining neuropathic pain during neuroimmune-glia cells interactions. This novel paradigm offers potential new targets for drug development against this type of chronic pain.
Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Recurrent exon-deleting activating mutations in AHR act as drivers of urinary tract cancer Bladder cancer has a high recurrence rate and low survival of advanced stage patients. Few genetic drivers of bladder cancer have thus far been identified. We performed in-depth structural variant analysis on whole-genome sequencing data of 206 metastasized urinary tract cancers. In ~10% of the patients, we identified recurrent in-frame deletions of exons 8 and 9 in the aryl hydrocarbon receptor gene (AHR{Delta}e8-9), which codes for a ligand-activated transcription factor. Pan-cancer analyses show that AHR{Delta}e8-9 is highly specific to urinary tract cancer and mutually exclusive with other bladder cancer drivers. In AHR{Delta}e8-9 the ligand-binding domain is disrupted and we show that this results in ligand-independent AHR-pathway activation. In bladder organoids, AHR{Delta}e8-9 induces a transformed phenotype that is characterized by upregulation of AHR target genes, downregulation of differentiation markers and upregulation of genes associated with stemness and urothelial cancer. Furthermore, AHR{Delta}e8-9 expression results in anchorage independent growth of bladder organoids, indicating tumorigenic potential. DNA-binding deficient AHR{Delta}e8-9 fails to induce transformation, suggesting a role for AHR target genes in the acquisition of the oncogenic phenotype. In conclusion, we show that AHR{Delta}e8-9 is a novel driver of urinary tract cancer and that the AHR pathway could be an interesting therapeutic target. Cell trajectory modeling identifies a primitive trophoblast state defined by BCAM enrichment In early placental development, progenitor cytotrophoblasts (CTBs) differentiate along one of two cellular trajectories: the villous or extravillous pathways. CTBs committed to the villous pathway fuse with neighboring CTBs to form the outer multinucleated syncytiotrophoblast (SCT), while CTBs committed to the extravillous pathway differentiate into invasive extravillous trophoblasts (EVT). Unfortunately, little is known about the processes controlling human CTB progenitor maintenance and differentiation. To address this, we established a single cell RNA sequencing (scRNA-seq) dataset from first trimester placentas to identify cell states important in trophoblast progenitor establishment, renewal, and differentiation. Multiple distinct trophoblast states were identified, representing progenitor CTBs, column CTBs, SCT precursors, and EVT. Lineage trajectory analysis identified a progenitor origin that was reproduced in human trophoblast stem cell organoids. Heightened expression of basal cell adhesion molecule (BCAM) defined this primitive state, where BCAM enrichment or gene silencing resulted in enhanced or diminished organoid growth. Together, this work describes at high-resolution trophoblast heterogeneity within the first trimester, resolves gene networks within human CTB progenitors, and identifies BCAM as a primitive progenitor marker and possible regulator.

Summary StatementLineage trajectory modeling identifies multiple human progenitor trophoblast states and defines trophoblast differentiation kinetics, where BCAM-expressing progenitors demonstrate enhanced regenerative ability. Omicron-specific mRNA vaccination alone and as a heterologous booster against SARS-CoV-2 The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has high transmissibility and recently swept the globe. Due to the extensive number of mutations, this variant has high level of immune evasion, which drastically reduced the efficacy of existing antibodies and vaccines. Thus, it is important to test an Omicron-specific vaccine, evaluate its immune response against Omicron and other variants, and compare its immunogenicity as boosters with existing vaccine designed against the reference wildtype virus (WT). Here, we generated an Omicron-specific lipid nanoparticle (LNP) mRNA vaccine candidate, and tested its activity in animals, both alone and as a heterologous booster to existing WT mRNA vaccine. Our Omicron-specific LNP-mRNA vaccine elicited strong and specific antibody response in vaccination-naive mice. Mice that received two-dose WT LNP-mRNA, the one mimicking the commonly used Pfizer/Moderna mRNA vaccine, showed a >40-fold reduction in neutralization potency against Omicron variant than that against WT two weeks post second dose, which further reduced to background level >3 months post second dose. As a booster shot for two-dose WT mRNA vaccinated mice, a single dose of either a homologous booster with WT LNP-mRNA or a heterologous booster with Omicron LNP-mRNA restored the waning antibody response against Omicron, with over 40-fold increase at two weeks post injection as compared to right before booster. Interestingly, the heterologous Omicron LNP-mRNA booster elicited neutralizing titers 10-20 fold higher than the homologous WT booster against the Omicron variant, with comparable titers against the Delta variant. All three types of vaccination, including Omicron mRNA alone, WT mRNA homologous booster, and Omicron heterologous booster, elicited broad binding antibody responses against SARS-CoV-2 WA-1, Beta, and Delta variants, as well as other Betacoronavirus species such as SARS-CoV, but not Middle East respiratory syndrome coronavirus (MERS-CoV). These data provided direct proof-of-concept assessments of an Omicron-specific mRNA vaccination in vivo, both alone and as a heterologous booster to the existing widely-used WT mRNA vaccine form. Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry Aiming to elucidate the grounds of system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (proteome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide (BPM) and 4-(N-succinimidylcarboxy)benzophenone (BPS). Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2. Of particular interest are the interactions of CYP2E1 with the latter two enzymes, which are involved in the synthesis and disposal of vasoactive and pro-inflammatory eicosanoids. We also detected the conjugates of CYP2E1 with UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1 (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes. Hepatic MIR20B promotes nonalcoholic fatty liver disease by suppressing PPARA BackgroundNon-alcoholic fatty liver disease (NAFLD) is characterized by excessive lipid accumulation and imbalances in lipid metabolism in the liver. Although nuclear receptors (NRs) play a crucial role in hepatic lipid metabolism, the underlying mechanisms of NR regulation in NAFLD remain largely unclear.

MethodsUsing network analysis and RNA-seq to determine the correlation between NRs and microRNA in human NAFLD patients, we revealed that MIR20B specifically targets PPARA. MIR20B mimic and anti-MIR20B were administered to human HepG2 and Huh-7 cells and mouse primary hepatocytes as well as high fat diet (HFD)- or methionine-deficient diet (MCD)-fed mice to verify the specific function of MIR20B in NAFLD. We tested the inhibition of the therapeutic effect of a PPAR agonist, fenofibrate, by MIR20B and the synergic effect of combination of fenofibrate with anti-MIR20B.

ResultsWe revealed that MIR20B specifically targets PPARA through miRNA regulatory network analysis of nuclear receptor genes in NAFLD. The expression of MIR20B was upregulated in free fatty acid (FA)- treated hepatocytes and the livers of both obesity-induced mice and NAFLD patients. Overexpression of MIR20B significantly increased hepatic lipid accumulation and triglyceride levels. Furthermore, MIR20B significantly reduced FA oxidation and mitochondrial biogenesis by targeting PPARA. In MIR20B-introduced mice, the effect of fenofibrate to ameliorate hepatic steatosis was significantly suppressed. Finally, inhibition of MIR20B significantly increased FA oxidation and uptake, resulting in improved insulin sensitivity and a decrease in NAFLD progression. Moreover, combination of fenofibrate and anti-MIR20B exhibited the synergic effect on improvement of NAFLD in MCD-fed mice.

ConclusionsTaken together, our results demonstrate that the novel MIR20B targets PPARA, plays a significant role in hepatic lipid metabolism, and present an opportunity for the development of novel therapeutics for NAFLD.

FundingThis research was funded by Korea Mouse Phenotyping Project (2016M3A9D5A01952411), the National Research Foundation of Korea (NRF) grant funded by the Korea government (2020R1F1A1061267, 2018R1A5A1024340, NRF-2021R1I1A2041463, 2020R1I1A1A01074940), and the Future-leading Project Research Fund (1.210034.01) of UNIST. Small extracellular vesicles from young mice prevent frailty, improve healthspan and decrease epigenetic age in old mice. Aging is associated with an increased risk of frailty, disability, comorbidities, institutionalization, falls, fractures, hospitalization, and mortality. Searching for strategies to delay the degenerative changes associated with aging and frailty is interesting. We treated old animals intravenously with small extracellular vesicles (sEVs) derived from adipose mesenchymal stem cells (ADSCs) of young animals, and we found an improvement of several functional parameters usually altered with aging, such as motor coordination, grip strength, fatigue resistance, fur regeneration, and renal function. Frailty index analysis showed that 40% of old control mice were frail, whereas none of the old ADSCs-sEVs treated mice were. Molecular and structural benefits in muscle and kidney accompanied this functional improvement. ADSCs-sEVs induced pro-regenerative effects and a decrease in oxidative stress, inflammation, and senescence markers. Moreover, predicted epigenetic age was lower in tissues of old mice treated with ADSCs-sEVs and their metabolome changed to a youth-like pattern. Finally, we gained some insight into the miRNAs contained in sEVs that might be, at least in part, responsible for the effects observed. We propose that young sEVs treatment can be beneficial against frailty and therefore can promote healthy aging. Three-Dimensional Loop Extrusion Loop extrusion convincingly describes how certain Structural Maintenance of Chromosome (SMC) proteins mediate the formation of large DNA loops. Yet, most of the existing computational models cannot reconcile recent in vitro observations showing that condensins can traverse each other, bypass large roadblocks and perform steps longer than its own size. To fill this gap, we propose a three-dimensional (3D) "trans-grabbing" model for loop extrusion which not only reproduces the experimental features of loop extrusion by one SMC complex, but also predicts the formation of so-called "Z-loops" via the interaction of two or more SMCs extruding along the same DNA substrate. By performing Molecular Dynamics simulations of this model we discover that the experimentally observed asymmetry in the different types of Z-loops is a natural consequence of the DNA tethering in vitro. Intriguingly, our model predicts this bias to disappear in absence of tethering and a third type of Z-loop, which has not yet been identified in experiments, to appear. Our model naturally explains road-block bypassing and the appearance of steps larger than the SMC size as a consequence of non-contiguous DNA grabbing. Finally, it is the first to our knowledge to address how Z-loops and bypassing might occur in a way that is broadly consistent with existing cis-only 1D loop extrusion models. 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. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Low heritability and high phenotypic plasticity of salivary cortisol in response to environmental heterogeneity in a wild pinniped Individuals are unique in how they interact with and respond to their environment. Correspondingly, unpredictable challenges or environmental stressors often produce an individualized response of the hypothalamic-pituitary-adrenal axis and its downstream effector cortisol. We used a fully crossed, repeated measures design to investigate the factors shaping individual variation in baseline cortisol in Antarctic fur seal pups and their mothers. Saliva samples were collected from focal individuals at two breeding colonies, one with low and the other with high density, during two consecutive years of contrasting food availability. Mothers and pups were sampled concurrently at birth and shortly before weaning, while pups were additionally sampled every 20 days. We found that heritability was low for baseline cortisol, while within-individual repeatability and among-individual variability were high. A substantial proportion of the variation in baseline cortisol could be explained in pups and mothers by a combination of intrinsic and extrinsic factors including sex, weight, day, season, and colony of birth. Our findings provide detailed insights into the individualization of endocrine phenotypes and their genetic and environmental drivers in a wild pinniped. Furthermore, the strong associations between cortisol and life history traits that we report in fur seals could have important implications for understanding the population dynamics of species impacted by environmental change. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice LinearTurboFold: Linear-Time Global Prediction of Conserved Structures for RNA Homologs with Applications to SARS-CoV-2 The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in SARS-CoV-2 genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length, and are thus infeasible for coronaviruses, which possess the longest genomes ([~]30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurbo-Folds purely in silico prediction not only is close to experimentally-guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5 and 3 UTRs ([~]29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies novel conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, siRNAs, CRISPR-Cas13 guide RNAs and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies, and will be a useful tool in fighting the current and future pandemics.

Significance StatementConserved RNA structures are critical for designing diagnostic and therapeutic tools for many diseases including COVID-19. However, existing algorithms are much too slow to model the global structures of full-length RNA viral genomes. We present LinearTurboFold, a linear-time algorithm that is orders of magnitude faster, making it the first method to simultaneously fold and align whole genomes of SARS-CoV-2 variants, the longest known RNA virus ([~]30 kilobases). Our work enables unprecedented global structural analysis and captures long-range interactions that are out of reach for existing algorithms but crucial for RNA functions. LinearTurboFold is a general technique for full-length genome studies and can help fight the current and future pandemics. 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. Choice of pre-processing pipeline influences clustering quality of scRNA-seq datasets BackgroundSingle-cell RNA sequencing (scRNA-seq) has quickly become one of the most dominant techniques in modern transcriptome assessment. In particular, 10X Genomics Chromium system, with its high throughput approach, turn key and thorough user guide made this cutting-edge technique accessible to many laboratories using diverse animal models. However, standard pre-processing, including the alignment and cell filtering pipelines might not be ideal for every organism or tissue. Here we applied an alternative strategy, based on the pseudoaligner kallisto, on twenty-two publicly available single cell sequencing datasets from a wide range of tissues of eight organisms and compared the results with the standard 10X Genomics Cell Ranger pipeline.

ResultsIn most of the tested samples, kallisto produced higher sequencing read alignment rates and total gene detection rates in comparison to Cell Ranger. Although datasets processed with Cell Ranger had higher cell counts, outside of human and mouse datasets, these additional cells were routinely of low quality, containing low gene detection rates. Thorough downstream analysis of one kallisto processed dataset, obtained from the zebrafish pineal gland, revealed clearer clustering, allowing the identification of an additional photoreceptor cell type that previously went undetected. The finding of the new cluster suggests that the photoreceptive pineal gland is essentially a bi-chromatic tissue containing both green and red cone-like photoreceptors and implies that the alignment and pre-processing pipeline can affect the discovery of biologically-relevant cell types.

ConclusionWhile Cell Ranger favors higher cell numbers, using kallisto results in datasets with higher median gene detection per cell. We could demonstrate that cell type identification was not hampered by the lower cell count, but in fact improved as a result of the high gene detection rate and the more stringent filtering. Depending on the acquired dataset, it can be beneficial to favor high quality cells and accept a lower cell count, leading to an improved classification of cell types. Distribution and molecular evolution of the anti-CRISPR family AcrIF7 Anti-CRISPRs are proteins capable of blocking CRISPR-Cas systems and typically encoded in mobile genetic elements. Since their discovery, numerous anti-CRISPR families have been identified. However, little is known about the distribution and sequence diversity of members within a family, nor how these traits influence the anti-CRISPRs function and evolution. Here we use AcrIF7 to explore the dissemination and molecular evolution of an anti-CRISPR family. We uncovered five sub-clusters and prevalent anti-CRISPR variants within the group. Remarkably, AcrIF7 homologs display high similarity despite their broad geographical, ecological and temporal distribution. Although mainly associated with Pseudomonas aeruginosa, AcrIF7 was identified in distinct genetic backgrounds indicating horizontal dissemination, primarily by phages. Using mutagenesis, we recreated variation observed in databases but also extended the sequence diversity of the group. Characterisation of the variants identified residues key for the anti-CRISPR function and other contributing to its mutational tolerance. Moreover, molecular docking revealed that variants with affected function lose key interactions with its CRISPR-Cas target. Analysis of publicly available data and the generated variants suggests that the dominant AcrIF7 variant corresponds to the minimal and optimal anti-CRISPR selected in the family. Our study provides a blueprint to investigate the molecular evolution of anti-CRISPR families. Ecological Differentiation Among Globally Distributed Lineages of the Rice Blast Fungus Pyricularia oryzae Many invasive fungal species coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. Here, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus (Pyricularia oryzae) covering a broad geographical range. We showed that the four lineages of P. oryzae were found in areas with different prevailing environmental conditions and types of rice grown, indicating niche separation. Pathogenicity tests revealed that specialization to rice subspecies contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong barriers to gene flow between these lineages. Our results demonstrate that the spread of a pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct invasive lineages. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Optical regulation of endogenous RhoA reveals selection of cellular responses by signal amplitude How protein signaling networks respond to different input strengths is an important but poorly understood problem in cell biology. For example, the small GTPase RhoA regulates both focal adhesion (FA) growth or disassembly, but whether RhoA serves as a switch selecting between cellular outcomes, or if outcomes are simply modulated by additional factors in the cell, is not clear. Here, we develop a photoswitchable RhoA guanine exchange factor, psRhoGEF, to precisely control endogenous RhoA activity. We also develop a FRET-based biosensor to allow visualization of RhoA activity together with psRhoGEF control. Using these new optical tools, we discover that low levels of RhoA activation preferentially induce FA disassembly in a Src-dependent manner, while high levels induce both FA growth and disassembly in a ROCK-dependent manner. Thus, rheostatic control of RhoA activation with photoswitchable RhoGEF reveals that cells can use signal amplitude to produce multiple responses to a single biochemical signal. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. The impact of housing conditions on mesenchymal stromal/stem cell populations differ between adipose tissue and skeletal muscle BackgroundIn pigs, the ratio between lean mass and fat mass in the carcass determines production efficiency and is strongly influenced by the number and size of cells in tissues. During growth, the increase in the number of cells results from the recruitment of different populations of multipotent mesenchymal stromal/stem cells (MSCs) residing in the tissues. We hypothesized that the impact of hygiene of housing conditions during growth on the proportions of MSCs in adipose tissue and skeletal muscle may differ between pigs with different residual feed intake (RFI), a measure of feed efficiency.

MethodsAt the age of 12 weeks, Large White pigs from two lines divergently selected for low and high RFI were housed in two contrasting hygiene conditions (good vs poor). After six weeks, pigs were slaughtered (n = 30; 5-9/group). Samples of subcutaneous adipose tissue and longissimus skeletal muscle were collected, and cells from the stromal vascular fraction (SVF), which includes mesenchymal stromal/stem cell populations, were isolated from each tissue. Adipose and muscle cell populations from the SVF were phenotyped by flow cytometry using antibodies that targeted different cell surface markers (CD45 to separate hematopoietic cells from MSCs; CD34, CD38, CD56 and CD140a to identify MSC populations with adipogenic and/or myogenic potential).

ResultsAdipose tissue and muscle shared some common MSC populations although MSC diversity was higher in muscle than in adipose tissue. In muscle, the CD45-CD56+CD34-CD140a+ and CD45-CD56+CD34+CD140a+ cell populations were abundant. Of these two cell populations, only the proportions of CD45-CD56+CD34+CD140a+ cells increased (P < 0.05) in pigs housed in poor hygiene as compared with pigs in good hygiene conditions. For the CD45-CD56-CD34- cell population, present in low proportion, there was an interaction between hygiene condition and genetic line (P < 0.05) with a decrease in low RFI pigs housed in poor hygiene conditions. In adipose tissue, the two abundant MSC populations were CD45-CD56-CD34- and CD45-CD56+CD34-. The proportion of CD45-CD56-CD34- cells increased (P < 0.05) whereas the proportion of CD45-CD56+CD34- tended to decrease (P < 0.1) in pigs housed in poor conditions. This study shows that the proportions of some MSC populations were affected by hygiene of housing conditions in a tissue-dependent manner in pigs of both RFI lines. It suggests that MSCs may play a significant role in adipose tissue and skeletal muscle homeostasis and may influence later growth and body composition in growing animals. Network Inference from Perturbation Time Course Data Networks underlie much of biology from subcellular to ecological scales. Yet, understanding what experimental data are needed and how to use them for unambiguously identifying the structure of even small networks remains a broad challenge. Here, we integrate a dynamic least squares framework into established modular response analysis (DL-MRA), that specifies sufficient experimental perturbation time course data to robustly infer arbitrary two and three node networks. DL-MRA considers important network properties that current methods often struggle to capture: (i) edge sign and directionality; (ii) cycles with feedback or feedforward loops including self-regulation; (iii) dynamic network behavior; (iv) edges external to the network; and (v) robust performance with experimental noise. We evaluate the performance of and the extent to which the approach applies to cell state transition networks, intracellular signaling networks, and gene regulatory networks. Although signaling networks are often an application of network reconstruction methods, the results suggest that only under quite restricted conditions can they be robustly inferred. For gene regulatory networks, the results suggest that incomplete knockdown is often more informative than full knockout perturbation, which may change experimental strategies for gene regulatory network reconstruction. Overall, the results give a rational basis to experimental data requirements for network reconstruction and can be applied to any such problem where perturbation time course experiments are possible.