phenopype: A phenotyping pipeline for Python O_LIDigital images are an intuitive way to capture, store and analyze organismal phenotypes. Many biologists are taking images to collect high-dimensional phenotypic information from specimens, which are key to investigating complex ecological, evolutionary and developmental phenomena, such as relationships between trait diversity and ecosystem function, multivariate natural selection, or developmental plasticity. As a consequence, images are being collected at ever increasing rates, but extraction of the contained phenotypic information still poses a veritable analytical bottleneck.
C_LIO_LIphenopype is a high throughput phenotyping pipeline for the programming language Python that aims at alleviating this bottleneck. The package facilitates immediate extraction of high dimensional phenotypic data from laboratory grade digital images with low levels of background noise and complexity. At the core, phenopype provides functionality for rapid signal processing based image preprocessing and segmentation, followed by trait extraction, visualization and numerous options for data export. The functionality is provided by wrapping low-level Python computer vision libraries (e.g. OpenCV) into accessible functions, facilitating their incorporation into scientific image analysis workflows. In addition, phenopype provides a project management ecosystem to further simplify rapid data collection and to increase reproducibility.
C_LIO_LIphenopype offers two different workflows that support users during different stages of scientific image analysis. The low-throughput workflow uses regular Python syntax and has greater flexibility at the cost of reproducibility, which is suitable for prototyping during the initial stage of a research project. The high-throughput workflow allows users to specify and store image-specific settings for analysis in human-readable YAML format, and then execute all functions in one step by means of an interactive parser. This approach facilitates rapid program-user interactions during batch processing, and greatly increases scientific reproducibility.
C_LIO_LIOverall, phenopype intends to make the features of powerful but technically involved low-level CV libraries available to biologists with little or no Python coding experience. Therefore, phenopype is aiming to augment, rather than replace the utility of existing Python CV libraries, allowing biologists to focus on rapid and reproducible data collection. Furthermore, image annotations produced by phenopype can be used as training data, thus presenting a stepping stone towards the application of deep learning architectures.
C_LI
A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea The importance of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae, Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including: (i) phylogenomics, (ii) the distribution of shared derived alleles, and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics, D-foil, and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggests that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species gene pools.
An automated workflow for multiplexed single-cell proteomics sample preparation at unprecedented sensitivity. The analysis of single-cell proteomes has recently become a viable complement to transcriptomics and genomics studies. Proteins are the primary driver of cellular functionality and mRNA levels are often an unreliable proxy of such. Therefore, the global analysis of the proteome is essential to study cellular identities. Multiplexed and label-free mass spectrometry-based approaches with single-cell resolution have lately attributed surprising heterogeneity to presumed homogenous cell populations. Even though specialized experimental designs and instrumentation have demonstrated remarkable advances, the efficient sample preparation of single cells still lag. Here, we introduce the proteoCHIP, a universal option for single-cell proteomics sample preparation at surprising sensitivity and throughput. The automated processing using a commercial system combining single-cell isolation and picoliter dispensing, the cellenONE(R), reduces final sample volumes to low nanoliters submerged in a hexadecane layer simultaneously eliminating error-prone manual sample handling and overcoming evaporation. The specialized proteoCHIP design allows direct injection of single cells via a standard autosampler resulting in around 1,500 protein groups per analytical run at remarkable reporter ion signal to noise while reducing or eliminating the carrier proteome. We identified close to 2,600 proteins across 170 multiplexed single cells from two highly similar human cell types. This dedicated loss-less workflow allows distinguishing in vitro co-differentiated cell types of self-organizing cardiac organoids based on indicative markers across 150 single cells. In-depth characterization revealed enhanced cellular motility of cardiac endothelial cells and sarcomere organization in cardiomyocytes. Our versatile and automated sample preparation has not only proven to be easily adaptable but is also sufficiently sensitive to drive biological applications of single-cell proteomics. FGF8-mediated signaling regulates tooth developmental pace during odontogenesis The developing human and mouse teeth constitute an ideal model system to study the regulatory mechanism underlying organ growth control due to the fact that their teeth share highly conserved and well-characterized developmental processes and their developmental tempo vary notably. In the current study, we manipulated heterogenous recombination between human and mouse dental tissues and demonstrate that the dental mesenchyme dominates the tooth developmental tempo and FGF8 could be a critical player during this developmental process. Forced activation of FGF8 signaling in the dental mesenchyme of mice promoted cell proliferation, prevented cell apoptosis via p38 and perhaps PI3K-Akt intracellular signaling, and impelled the transition of the cell cycle from G1- to S-phase in the tooth germ, resulting in the slowdown of the tooth developmental pace. Our results provide compelling evidence that extrinsic signals can profoundly affect tooth developmental tempo and the dental mesenchymal FGF8 could be a pivotal factor in controlling developmental pace in a non-cell-autonomous manner during mammalian odontogenesis. OncoLoop: A network-based precision cancer medicine framework At present, prioritizing cancer treatments at the individual patient level remains challenging, and performing co-clinical studies using patient-derived models in real-time is often not feasible. To circumvent these challenges, we introduce OncoLoop, a precision medicine framework to predict and validate drug sensitivity in human tumors and their pre-existing high-fidelity (cognate) model(s) by leveraging perturbational profiles of clinically-relevant oncology drugs. As proof-of-concept, we applied OncoLoop to prostate cancer (PCa) using a series of genetically-engineered mouse models (GEMMs) that recapitulate a broad spectrum of disease states, including castration-resistant, metastatic, and neuroendocrine prostate cancer. Interrogation of published cohorts using Master Regulator (MR) conservation analysis revealed that most patients were represented by at least one cognate GEMM-derived tumor (GEMM-DT). Drugs predicted to invert MR activity in patients and their cognate GEMM-DTs were successfully validated, including in two cognate allografts and one cognate patient-derived xenograft (PDX). OncoLoop is a highly generalizable framework that can be extended to other cancers and potentially other diseases.

Significance StatementOncoLoop is a transcriptomic-based experimental and computational framework that can support rapid-turnaround co-clinical studies to identify and validate drugs for individual patients, which can then be readily adapted to clinical practice. This framework should be applicable in many cancer contexts for which appropriate models and drug perturbation data are available. Mapping the genetic and environmental aetiology of autistic traits in Sweden and the United Kingdom BackgroundAutistic traits are influenced by both genetic and environmental factors, and are known to vary geographically in prevalence. But to what extent does their aetiology also vary from place to place?

MethodsWe applied a novel spatial approach to data on autistic traits from two large twin studies, the Child and Adolescent Twin Study in Sweden (CATSS; N=16,677, including 8,307 twin pairs) and the Twins Early Development Study in the UK (TEDS; N=11,594, including 5,796 twin pairs), to explore how the influence of nature and nurture on autistic traits varies from place to place.

ResultsWe present maps of gene- and environment-by geography interactions in Sweden and the United Kingdom (UK), showing geographical variation in both genetic and environmental influences across the two countries. In Sweden genetic influences appear higher in the far south and in a band running across the centre of the country. Environmental influences appear greatest in the south and north, with reduced environmental influence across the central band. In the UK genetic influences appear greater in the south, particularly in more central southern areas and the southeast, the Midlands and the north of England. Environmental influences appear greatest in the south and east of the UK, with less influence in the north and the west.

ConclusionsWe hope this systematic approach to identifying aetiological interactions will inspire research to examine a wider range of previously unknown environmental influences on the aetiology of autistic traits. By doing so, we will gain greater understanding of how these environments draw out or mask genetic predisposition and interact with other environmental influences in the development of autistic traits. 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. 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. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Action at a distance: Defects in division plane positioning in the root meristematic zone affect cell organization in the differentiation zone Cell division plane orientation is critical for plant and animal development and growth. TANGLED1 (TAN1) and AUXIN-INDUCED-IN-ROOT-CULTURES9 (AIR9) are division-site localized microtubule-binding proteins required for division plane positioning. tan1 and air9 Arabidopsis thaliana single mutants have minor or no noticeable phenotypes but the tan1 air9 double mutant has synthetic phenotypes including stunted growth, misoriented divisions, and aberrant cell-file rotation in the root differentiation zone. These data suggest that TAN1 plays a role in nondividing cells. To determine whether TAN1 is required in elongating and differentiating cells in the tan1 air9 double mutant, we limited its expression to actively dividing cells using the G2/M-specific promoter of the syntaxin KNOLLE (pKN:TAN1-YFP). Unexpectedly, in addition to rescuing division plane defects, pKN:TAN1-YFP rescued root growth and the root differentiation zone cell file rotation defects in the tan1 air9 double mutant. This suggests that defects that occur in the meristematic zone later affect the organization of elongating and differentiating cells.

Summary StatementExpression of TAN1 in the root meristematic zone rescues cell file rotation defects in tan1 air9 mutants, suggesting defects that occur in mitosis may influence organization of nondividing cells. NOMA: a high-throughput microchip for robust, sequential measurements of secretions from the same single-cells Despite recent advances in single-cell analysis technologies, lacking simple methods to keep the live single-cells traceable for longitudinal detection reliably poses a significant obstacle in single-cell secretion analysis. Here we reported the high-density NOMA (narrow-opening microwell array) microchip that realized the retention of[&ge;] 97% of trapped single cells in dedicated spatial locations during repetitive detection procedures, verified with both adherent and suspension cells by two researchers independently. We applied it in monitoring single-cell protein secretions sequentially from the same single cells, and we found the digital protein secretion patterns dominate the protein secretion. We also demonstrated the microchip for longitudinally tracking IL-8 and the CD81+EV secretions from the same single-cells over days, which revealed the presence of "super secretors" within the cell population be more persistent to secrete protein or extracellular vesicle for an extended period. The NOMA platform reported here is simple, robust, and easy to operate for tracking sequential measurements from the same single cells, representing a novel and informative tool to inspire new observations in biomedical research.

Table of Contents Graphic

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/427049v4_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@aa69c7org.highwire.dtl.DTLVardef@6bc79org.highwire.dtl.DTLVardef@3fcbbforg.highwire.dtl.DTLVardef@597f4b_HPS_FORMAT_FIGEXP  M_FIG C_FIG Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. Correlative all-optical quantification of mass density and mechanics of sub-cellular compartments with fluorescence specificity Quantitative measurements of physical parameters become increasingly important for understanding biological processes. Brillouin microscopy (BM) has recently emerged as one technique providing the 3D distribution of viscoelastic properties inside biological samples -- so far relying on the implicit assumption that refractive index (RI) and density can be neglected. Here, we present a novel method (FOB microscopy) combining BM with optical diffraction tomography and epi-fluorescence imaging for explicitly measuring the Brillouin shift, RI and absolute density with specificity to fuorescently labeled structures. We show that neglecting the RI and density might lead to erroneous conclusions. Investigating the nucleoplasm of wild-type HeLa cells, we find that it has lower density but higher longitudinal modulus than the cytoplasm. Thus, the longitudinal modulus is not merely sensitive to the water content of the sample -- a postulate vividly discussed in the field. We demonstrate the further utility of FOB on various biological systems including adipocytes and intracellular membraneless compartments. FOB microscopy can provide unexpected scientific discoveries and shed quantitative light on processes such as phase separation and transition inside living cells. Heritable changes in division speed accompany the diversification of singleT cell fate Rapid clonal expansion of antigen specific T cells is a fundamental feature of adaptive immune responses. It enables the outgrowth of an individual T cell into thousands of clonal descendants that diversify into short-lived effectors and long-lived memory cells. Clonal expansion is thought to be programmed upon priming of a single naive T cell and then executed by homogenously fast divisions of all of its descendants. However, the actual speed of cell divisions in such an emerging  T cell family has never been measured with single-cell resolution. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. This comprehensive mapping of T cell family trees identifies a short burst phase, in which division speed is homogenously fast and maintained independent of external cytokine availability or continued T cell receptor stimulation. Thereafter, however, division speed diversifies and model-based computational analysis using a novel Bayesian inference framework for tree-structured data reveals a segregation into heritably fast and slow dividing branches. This diversification of division speed is preceded already during the burst phase by variable expression of the interleukin-2 receptor alpha chain. Later it is accompanied by selective expression of memory marker CD62L in slower dividing branches. Taken together, these data demonstrate that T cell clonal expansion is structured into subsequent burst and diversification phases the latter of which coincides with specification of memory vs. effector fate.

SignificanceRapid clonal expansion of antigen-specific T cells is a fundamental feature of adaptive immune responses. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. Bayesian inference of tree-structured data reveals that clonal expansion is divided into a homogenously fast burst phase encompassing two to three divisions and a subsequent diversification phase during which T cells segregate into quickly dividing effector T cells and more slowly cycling memory precursors. Our work highlights cell cycle speed as a major heritable property that is regulated in parallel to key lineage decisions of activated T cells. An alternative UPF1 isoform drives conditional remodeling of nonsense-mediated mRNA decay The nonsense-mediated mRNA decay (NMD) pathway monitors translation termination to degrade transcripts with premature stop codons and regulate thousands of human genes. Here we show that an alternative mammalian-specific isoform of the core NMD factor UPF1, termed UPF1LL, enables condition-dependent remodeling of NMD specificity. Previous studies indicate that the extension of a conserved regulatory loop in the UPF1LL helicase core confers a decreased propensity to dissociate from RNA upon ATP hydrolysis relative to the major UPF1 isoform, designated UPF1 SL. Using biochemical and transcriptome-wide approaches, we find that UPF1LL overcomes the protective RNA binding proteins PTBP1 and hnRNP L to preferentially bind and down-regulate long 3UTRs normally shielded from NMD. Unexpectedly, UPF1LL supports induction of NMD on new populations of substrate mRNAs in response to activation of the integrated stress response and impaired translation efficiency. Thus, while canonical NMD is abolished by moderate translational repression, UPF1LL activity is enhanced, providing a mechanism to rapidly rewire NMD specificity in response to cellular stress. Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Modeling miRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Estimating the fundamental niche: accounting for the uneven availability of existing environmental conditions in the calibration area Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche. 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. Understanding transformation tolerant visual object representations in the human brain and convolutional neural networks Forming transformation-tolerant object representations is critical to high-level primate vision. Despite its significance, many details of tolerance in the human brain remain unknown. Likewise, despite the ability of convolutional neural networks (CNNs) to exhibit human-like object categorization performance, whether CNNs form tolerance similar to that of the human brain is unknown. Here we provide the first comprehensive documentation and comparison of three tolerance measures in the human brain and CNNs. We measured fMRI responses from human ventral visual areas to real-world objects across both Euclidean and non-Euclidean feature changes. In single fMRI voxels in higher visual areas, we observed robust object response rank-order preservation across feature changes. This is indicative of functional smoothness in tolerance at the fMRI meso-scale level that has never been reported before. At the voxel population level, we found highly consistent object representational structure across feature changes towards the end of ventral processing. Rank-order preservation, consistency, and a third tolerance measure, cross-decoding success (i.e., a linear classifiers ability to generalize performance across feature changes) showed an overall tight coupling. These tolerance measures were lower for Euclidean than non-Euclidean feature changes in lower visual areas, but increased over the course of ventral processing in most cases. These characteristics of tolerance, however, were absent in eight CNNs pretrained with ImageNet images with varying network architecture, depth, the presence/absence of recurrent processing, or whether a network was pretrained with the original or stylized ImageNet images that encouraged shape processing. Most notably, CNNs do not show increased representational consistency across feature changes at the higher layers. CNNs thus do not appear to develop the same kind of tolerance as the human brain over the course of visual processing.

Significant StatementPerceiving object identity among changes in non-identity features and forming transformation-tolerant object representations is essential to high-level primate vision. Here we provide a comprehensive documentation and comparison of three tolerance measures between the human brain and CNNs pretrained for object classification. While all three measures show increased tolerance in the human brain across four types of feature changes towards the end of ventral visual processing, CNNs fail to develop the same kind of tolerance with visual processing. EGFR-driven profibrotic signals in AKI-to-CKD transition are blocked by TNF inhibition: Opportunities for etanercept treatment Proximal-tubule-cell EGFR activation mediates tubule cell proliferation and repair early after kidney injury, while sustained EGFR activation causes kidney fibrosis. Inflammation is a key driver of AKI-to-CKD transition and fibrosis, but mechanisms of EGFR-driven profibrotic responses are not well understood. In a mouse model of AKI-to-CKD transition and CKD progression, we show that EGFR-inhibition significantly reduced kidney expression of many immunoregulatory molecules already by day two after injury, including the potent inflammatory cytokine tumor-necrosis-factor (TNF). Single nuclei RNA-sequencing analysis showed that macrophages were among the main early cellular sources of TNF in the injured kidney. In vitro, EGFR activation in macrophages increased macrophage TNF expression, while EGFR inhibition in vivo reduced kidney macrophage accumulation, as early as two days after injury. Thus, profibrotic EGFR signaling increases kidney TNF levels both directly and indirectly. TNF inhibition did not alter tubule EGFR activation and in contrast to EGFR inhibition did not reduce early macrophage accumulation in the kidney, suggesting that TNF does not promote early infiltration of immune cells in the kidney, but rather regulates profibrotic functions of kidney and/or immune cells. TNF inhibition with etanercept in vivo in AKI-injured mice downregulated a number of cytokines, including TNF itself. This cytokine downregulation overlapped with EGFR inhibition, and additional non-overlapping downregulated cytokines shared the same functions, predicting that TNF inhibition would prevent AKI-to-CKD transition like EGFR inhibition. Indeed, TNF-inhibition with etanercept reduced AKI-induced fibrosis to the same degree as the EGFR-inhibition, while the combination of both treatments showed no additive effect. In conclusion, our results identify TNF as a downstream effector of profibrotic EGFR activation and motivate the examination of TNF pathway inhibition in human AKI or CKD.

Translational statementProximal tubule epidermal-growth-factor (EGFR) activation in mice and likely also in humans drives inflammation in the kidney after acute-kidney-injury (AKI) and causes chronic-kidney-disease (CKD) with fibrosis in a process termed AKI-to-CKD transition. Recent retrospective data shows that patients treated with TNF inhibitors show decreased incidence and progression of CKD. Our current work shows that TNF inhibition in mice is equally as effective as EGFR inhibition in preventing AKI-to-CKD transition and fibrosis. Thus, our results have translational potential and may stimulate examination of short-term TNF inhibition in AKI to prevent AKI-to-CKD transition or possibly of longer-term TNF inhibition in CKD to prevent CKD progression. Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Scalable parametric encoding of multiple modalities AO_SCPLOWBSTRACTC_SCPLOWA flexible model is introduced which shares ideas with the Autoencoder, Canonical Correlation Analysis, Singular Value Decomposition, and Procrustes Analysis. It is proposed to find relevant maps to transform multiple datasets of various types from one modality to another. Here, the Generative Encoder is used to transform spatial gene expression from breast tissue, to the images of histology tissue measured with Spatial Transcriptomics. The model is directly interpretable given all parameters are linked to the data space. It is scalable on Big Data, training reasonably on several thousand RGB images of 100 by 100 pixels in under an hour, which equates to 30,000 pixel features per sample image. 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. Atypical integrative element with strand-biased circularization activity assists interspecies antimicrobial resistance gene transfer from Vibrio alfacsensis The exchange of antimicrobial resistance (AMR) genes between aquaculture and terrestrial microbial populations has emerged as a serious public health concern. However, the nature of the mobile genetic elements in marine bacteria is poorly documented. To gain insight into the genetic mechanisms underlying AMR gene transfer from marine bacteria, we mated a multidrug-resistant Vibrio alfacsensis strain with an Escherichia coli strain, and then determined the complete genome sequences of the donor and transconjugant strains. Sequence analysis revealed a multidrug resistance conjugative plasmid in the donor strain, which was integrated into the chromosome of the recipient. The plasmid backbone in the transconjugant chromosome was flanked by two copies of a 7.1 kb unclassifiable integrative element harboring a {beta}-lactamase gene. The 7.1 kb element (named SE-6945) and previously reported Tn6283 share four coding sequences, two of which encode the catalytic R-H-R-Y motif of tyrosine recombinases. Polymerase chain reaction and sequencing experiments revealed that SE-6945 generates a circular copy of its one specific strand without generating an empty donor site, in contrast to the movement of integron gene cassette or ICE/IMEs discovered to date. Copy numbers and location of SE-6945 in the chromosome affected antibiotic resistance levels of the transconjugants. SE-6945/Tn6283-related integrative elements are collectively termed SEs (strand-biased circularizing integrative elements). SEs were identified in the genomes of other Vibrio species. Overall, these results suggest that in addition to conjugative plasmids, SEs are involved in the spread of AMR genes among marine bacteria.