Embracing imperfection: machine-assisted invertebrate classification in difficult datasets O_LIDespite growing concerns over the health of global invertebrate diversity, terrestrial invertebrate monitoring efforts remain poorly geographically distributed. Machine-assisted classification has been proposed as a potential solution to quickly gather large amounts of data; however, previous studies have often used unrealistic or idealized datasets to train their models.
C_LIO_LIIn this study, we describe a practical methodology for including machine learning in ecological data acquisition pipelines. Here we train and test machine learning algorithms to classify over 56,000 bulk terrestrial invertebrate specimens from morphometric data and contextual metadata. All vouchered specimens were collected in pitfall traps by the National Ecological Observatory Network (NEON) at 27 locations across the United States in 2016. Specimens were photographed, and morphometric data was extracted as feature vectors using ImageJ. Issues stemming from inconsistent taxonomic label specificity were resolved by making classifications at the lowest identified taxonomic level (LITL). Taxa with too few specimens to be included in the training dataset were classified by the model using zero-shot classification.
C_LIO_LIWhen classifying specimens that were known and seen by our models, we reached an accuracy of 72.7% using extreme gradient boosting (XGBoost) at the LITL. Models that were trained without contextual metadata underperformed models with contextual metadata by an average of 7.2%. We also classified invertebrate taxa that were unknown to the model using zero-shot classification, with an accuracy of 39.4%, resulting in an overall accuracy of 71.5% across the entire NEON dataset.
C_LIO_LIThe general methodology outlined here represents a realistic application of machine learning as a tool for ecological studies. Hierarchical and LITL classifications allow for flexible taxonomic specificity at the input and output layers. These methods also help address the  long tail problem of underrepresented taxa missed by machine learning models. Finally, we encourage researchers to consider more than just morphometric data when training their models, as we have shown that the inclusion of contextual metadata can provide significant improvements to accuracy.
C_LI
Assessing the impact of grassland management on landscape multifunctionality Land-use intensification has contrasting effects on different ecosystem services, often leading to land-use conflicts. While multiple studies have demonstrated how landscape-scale strategies can minimise the trade-off between agricultural production and biodiversity conservation, little is known about which land-use strategies maximise the landscape-level supply of multiple ecosystem services (landscape multifunctionality), a common goal of stakeholder communities.

We combine comprehensive data collected from 150 German grassland sites with a simulation approach to identify landscape compositions, with differing proportions of low-, medium-, and high-intensity grasslands, that minimise trade-offs between the six main grassland ecosystem services prioritised by local stakeholders: biodiversity conservation, aesthetic value, productivity, carbon storage, foraging, and regional identity. Results are made accessible through an online tool that provides information on which compositions best meet any combination of user-defined priorities (https://neyret.shinyapps.io/landscape_composition_for_multifunctionality/).

Results show that an optimal landscape composition can be identified for any pattern of ecosystem service priorities. However, multifunctionality was similar and low for all landscape compositions in cases where there are strong trade-offs between services (e.g. aesthetic value and fodder production), where many services were prioritised, and where drivers other than land use played an important role. We also found that if moderate service levels are deemed acceptable, then strategies in which both high and low intensity grasslands are present can deliver landscape multifunctionality. The tool presented can aid informed decision-making by predicting the impact of future changes in landscape composition, and by allowing for the relative roles of stakeholder priorities and biophysical trade-offs to be understood by scientists and practitioners alike.

HighlightsO_LIAn online tool identifies optimal landscape compositions for desired ecosystem services
C_LIO_LIWhen the desired services are synergic, the optimum is their common best landscape composition
C_LIO_LIWhen the desired services trade-off, a mix of grassland intensity is most multifunctional
C_LIO_LISuch tools could support decision-making processes and aid conflict resolution
C_LI

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=111 SRC="FIGDIR/small/208199v5_ufig1.gif" ALT="Figure 1">
View larger version (27K):
org.highwire.dtl.DTLVardef@3273f6org.highwire.dtl.DTLVardef@5b3a5aorg.highwire.dtl.DTLVardef@153f526org.highwire.dtl.DTLVardef@103cbea_HPS_FORMAT_FIGEXP  M_FIG C_FIG
A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. Targeting CISH enhances natural cytotoxicity receptor signaling and reduces NK cell exhaustion to improve solid tumor immunity Cytokine inducible SH2-containing protein (CISH) is a natural killer (NK) cell negative regulator of cytokine signaling pathway. To further understand CISH functions in NK cells, we developed a conditional Cish-deficient mouse model in NK cells (Cishfl/flNcr1Ki/+). We detected no developmental or homeostatic difference in NK cells.

However, global gene expression of Cishfl/flNcr1Ki/+ NK cells compared to Cish+/+Ncr1Ki/+ NK cells revealed upregulation of pathways and genes associated with NK cell cycling and activation. We show that CISH does not only regulate interleukin-15 (IL-15) signaling pathways but also natural cytotoxicity receptors (NCR) pathways. Indeed, CISH protein expression level increases upon NCR triggering. Primed Cishfl/flNcr1Ki/+ NK cells display increased activation upon NCR stimulation. Cishfl/flNcr1Ki/+ NK cells display lower activation thresholds and Cishfl/flNcr1Ki/+ mice are more resistant to tumor metastasis. Remarkably, we found that Cishfl/flNcr1Ki/+ mice were also more resistant to primary breast cancer growth in addition to superior control of spontaneous tumor metastasis. CISH deletion favors NK cell accumulation to the primary tumor, optimizes NK cell killing properties and decreases TIGIT immune checkpoint receptor expression, limiting NK cell exhaustion. Finally, we argue that specifically enhancing NK cell function is sufficient to boost anti-tumor response to both primary and secondary tumor models. Using CRISPRi, we then targeted CISH in human NK-92 or primary NK cells. According to the results in our mouse model, CISH deletion favors NCR signaling and anti-tumor functions in human NK cells. Our results validate CISH as an emerging therapeutic target to enhance NK cell immunotherapy. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. 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. Genome-wide association study suggests an independent genetic basis of zinc and cadmium concentrations in fresh sweet corn kernels Despite being one of the most consumed vegetables in the United States, the elemental profile of sweet corn (Zea mays L.) is limited in its dietary contributions. To address this through genetic improvement, a genome-wide association study was conducted for the concentrations of 15 elements in fresh kernels of a sweet corn association panel. In concordance with mapping results from mature maize kernels, we detected a probable pleiotropic association of zinc and iron concentrations with nicotianamine synthase5 (nas5), which purportedly encodes an enzyme involved in synthesis of the metal chelator nicotianamine. Additionally, a pervasive association signal was identified for cadmium concentration within a recombination suppressed region on chromosome 2. The likely causal gene underlying this signal was heavy metal ATPase3 (hma3), whose counterpart in rice, OsHMA3, mediates vacuolar sequestration of cadmium and zinc in roots, whereby regulating zinc homeostasis and cadmium accumulation in grains. In our association panel, hma3 associated with cadmium but not zinc accumulation in fresh kernels. This finding implies that selection for low cadmium will not affect zinc levels in fresh kernels. Although less resolved association signals were detected for boron, nickel, and calcium, all 15 elements were shown to have moderate predictive abilities via whole-genome prediction. Collectively, these results help enhance our genomics-assisted breeding efforts centered on improving the elemental profile of fresh sweet corn kernels. Deep Heterogeneous Dilation of LSTM for Transient-phase Gesture Prediction through High-density Electromyography: Towards Application in Neurorobotics Deep networks have been recently proposed to estimate motor intention using conventional bipolar surface electromyography (sEMG) signals for myoelectric control of neurorobots. In this regard, Deepnets are generally challenged by long training times (affecting practicality and calibration), complex model architectures (affecting the predictability of the outcomes), and a large number of trainable parameters (increasing the need for big data). Capitalizing on our recent work on homogeneous temporal dilation in a Recurrent Neural Network (RNN) model, this paper proposes, for the first time, heterogeneous temporal dilation in an LSTM model and applies that to high-density surface electromyography (HD-sEMG), allowing for the decoding of dynamic temporal dependencies with tunable temporal foci. In this paper, a 128-channel HD-sEMG signal space is considered due to the potential for enhancing the spatiotemporal resolution of human-robot interfaces. Accordingly, this paper addresses a challenging motor intention decoding problem of neurorobots, namely, transient intention identification. Our approach uses only the dynamic and transient phase of gesture movements when the signals are not stabilized or plateaued, which can significantly enhance the temporal resolution of human-robot interfaces. This would eventually enhance seamless real-time implementations. Additionally, this paper introduces the concept of "dilation foci" to modulate the modeling of temporal variation in transient phases. In this work a high number (e.g., 65) of gestures is included, which adds to the complexity and significance of the understudied problem. Our results show state-of-the-art performance for gesture prediction in terms of accuracy, training time, and model convergence. 3D Reconstructions of Mouse Skeletal Muscle Reveal a Decrease in the MICOS Complex and Altered Mitochondrial Networks Background: Skeletal muscle gradually loses mass, strength, endurance, and oxidative capacity during aging. Studies of bioenergetics and protein turnover show that mitochondria mediate this decline in function. Mitochondria cristae are essential for the production of ATP. While mitochondrial aging is associated with endoplasmic reticulum stress, fragmented mitochondria, and decreased mitochondrial capacity, the genes associated with morphological changes in mitochondria during aging still requires further elucidation Further, it is not completely understood how 3D mitochondrial networks and the specialization of mitochondria are altered during aging. Methods: We measured changes in mitochondrial morphology and mitochondrial connectivity during the aging of the mouse gastrocnemius muscle through serial block facing-scanning electron microscopy and 3D reconstruction. Nanotunnels were also measured through 3D reconstruction. CRISPR/Cas9 KO was performed to examine changes in mitochondria upon loss of the mitochondrial contact site and cristae organizing system (MICOS) complex and optic atrophy 1 (OPA-1). Metabolomics were used to find key metabolite pathways changed upon off of the MICOS complex. Results: We found changes in mitochondrial network configuration, nanotunneling, size, shape, number, contact sites, cristae organizing system (MICOS) dynamics and gene expression in skeletal muscle during aging. We also found an association of OPA-1 and the MICOS complex in the gastrocnemius with mitochondrial aging. Further, the loss of the MICOS complex was linked with decreased oxidative capacity, and altered mitochondrial metabolism. Conclusions: To our knowledge, studies of mitochondrial changes in skeletal muscle throughout aging remains limited. 3D reconstructions of nanotunnels elucidated novel patterns of mitochondrial aging in skeletal muscle. MICOS proteins decreased with age and mitochondrial morphology was similar between aged skeletal muscle and that of young mice with MICOS protein loss. In tandem, our data suggest a relationship between the MICOS complex and aging, which could be potentially linked to disease states with additional 3D reconstruction. 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. Yersinia remodels epigenetic histone modifications in human macrophages Various pathogens systematically reprogram gene expression in macrophages, but the underlying mechanisms are largely unknown. We investigated whether the enteropathogen Yersinia enterocolitica alters chromatin states to reprogram gene expression in primary human macrophages. Genome-wide chromatin immunoprecipitation (ChIP) seq analyses showed that pathogen-associated molecular patterns (PAMPs) induced up- or down-regulation of histone modifications (HMod) at approximately 14500 loci in promoters and enhancers. Effectors of Y. enterocolitica reorganized about half of these dynamic HMod, with the effector YopP being responsible for about half of these modulatory activities. The reorganized HMod were associated with genes involved in immune response and metabolism. Remarkably, the altered HMod also associated with 61 % of all 534 known Rho GTPase pathway genes, revealing a new level in Rho GTPase regulation and a new aspect of bacterial pathogenicity. Changes in HMod were associated to varying degrees with corresponding gene expression, e. g. depending on chromatin localization and cooperation of the HMod. In summary, infection with Y. enterocolitica remodels HMod in human macrophages to modulate key gene expression programs of the innate immune response.

Author SummaryHuman pathogenic bacteria can affect epigenetic histone modifications to modulate gene expression in host cells. However, a systems biology analysis of this bacterial virulence mechanism in immune cells has not been performed. Here we analyzed genome-wide epigenetic histone modifications and associated gene expression changes in primary human macrophages infected with enteropathogenic Yersinia enterocolitica. We demonstrate that Yersinia virulence factors extensively modulate histone modifications and associated gene expression triggered by the pathogen-associated molecular patterns (PAMPs) of the bacteria. The epigenetically modulated genes are involved in several key pathways of the macrophage immune response, including the Rho GTPase pathway, revealing a novel level of Rho GTPase regulation by a bacterial pathogen. Overall, our findings provide an in-depth view of epigenetic and gene expression changes during host-pathogen interaction and might have further implications for understanding of the innate immune memory in macrophages. Characterising ontogeny of numerosity discrimination in zebrafish reveals the use of multiple, numerical and non- numerical mechanisms. A sense of non-symbolic numerical magnitudes is widespread in the animal kingdom and has been documented in adult zebrafish. Here we investigated the ontogeny of this ability using a group size preference task in juvenile zebrafish. Fish showed group size preference from 21 days post fertilization (dpf) and reliably chose the larger group when presented with discriminations of between 1 vs. 3, 2 vs. 5 and 2 vs. 3 conspecifics but not 2 vs. 4 conspecifics. When the ratio between the number of conspecifics in each group was maintained at 1:2, fish could discriminate between 1 vs. 2 individuals and 3 vs. 6, but again, not when given a choice between 2 vs. 4 individuals. These findings are in agreement with studies in other species suggesting the systems involved in quantity representation do not operate separately from other cognitive mechanisms. Rather they suggest quantity processing in fish may be the result of an interplay between attentional, cognitive and memory-related mechanisms as in humans and other animals. Our results emphasise the potential of the use of zebrafish to explore the genetic and neural processes underlying the ontogeny and function of number cognition. "Neighbourhood watch" model: embryonic epiblast cells assess positional information in relation to their neighbours In many developing and regenerating systems, tissue pattern is established through gradients of informative morphogens, but we know little about how cells interpret these. Using experimental manipulation of early chick embryos including misexpression of an inducer (VG1 or ACTIVIN) and an inhibitor (BMP4), we test two alternative models for their ability to explain how the site of primitive streak formation is positioned relative to the rest of the embryo. In one model, cells read morphogen concentrations cell-autonomously. In the other, cells sense changes in morphogen status relative to their neighbourhood. We find that only the latter model can account for the experimental results, including some counter-intuitive predictions. This mechanism (which we name "neighbourhood watch" model) illuminates the classic "French Flag Problem" and how positional information is interpreted by a sheet of cells in a large developing system.

Summary statementIn a large developing system, the chick embryo before gastrulation, cells interpret gradients of positional signals relative to their neighbours to position the primitive streak, establishing bilateral symmetry. S-acylation stabilizes ligand-induced receptor kinase complex formation during plant pattern-triggered immune signalling. Plant receptor kinases are key transducers of physical extracellular stimuli, such as the presence of beneficial or pathogenic microbes or secreted signalling molecules. Receptor kinases are regulated by numerous post-translational modifications. Here, using the bacterial flagellin perceiving receptor kinase FLS2, we show that S-acylation at an evolutionarily conserved cysteine is crucial for function. S-acylation involves the addition of long-chain fatty acids to cysteine residues within proteins, altering their biophysical properties and behaviour within the membrane environment. We observe S-acylation of FLS2 at C-terminal kinase domain cysteine residues within minutes of treatment with flg22 ligand, and in a BAK1 co-receptor dependent manner. We demonstrate that S-acylation is essential for FLS2-mediated immune signalling, including anti-bacterial immunity. Similarly, mutating the corresponding conserved cysteine residue in the immune receptor kinase EFR supressed elf18 mediated signalling. Biochemical analysis of unstimulated and activated FLS2 containing complexes using detergents and native membrane DIBMA nanodiscs indicates that S-acylation assists the stabilization of activated receptor kinase complexes withing the membrane environment to increase signalling efficiency. Modeling microRNA-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. Aldh2 is a lineage-specific metabolic gatekeeper in melanocyte stem cells Melanocyte stem cells (McSCs) in zebrafish serve as an on-demand source of melanocytes during growth and regeneration, but metabolic programs associated with their activation and regenerative processes are not well known. Here, using live imaging coupled with scRNA-sequencing, we discovered that quiescent McSCs during regeneration activate a dormant embryonic neural crest transcriptional program followed by an aldehyde dehydrogenase (Aldh) 2 metabolic switch to generate progeny. Unexpectedly, while ALDH2 is well known for its aldehyde clearing mechanisms we find that in regenerating McSCs, Aldh2 activity is required to generate formate - the one-carbon (1C) building block for nucleotide biosynthesis - through formaldehyde metabolism. Consequently, we find that disrupting the 1C cycle with low-doses of methotrexate caused melanocyte regeneration defects. In the absence of Aldh2, we find that purines (but not pyrimidines) are the metabolic end product sufficient for activated McSCs to generate progeny. Together, our work reveals McSCs undergo a two-step cell state transition during regeneration, and that the reaction products of Aldh2 enzymes have tissue-specific stem cell functions that meet metabolic demands in regeneration.

SUMMARY STATEMENTIn melanocyte regeneration, quiescent McSCs respond by re-expressing a neural crest identity, followed by an Aldh2-dependent metabolic switch to generate progeny. Exploring functional protein covariation across single cells using nPOP Many biological processes, such as the cell division cycle, are reflected in protein covariation across single cells. This covariation can be quantified and interpreted by single-cell proteomics with sufficiently high throughput and accuracy. Toward this goal, we developed the nano-ProteOmic sample Preparation (nPOP) method for single-cell proteomics. nPOP uses piezo acoustic dispensing to isolate individual cells in 300 picoliter volumes and performs all subsequent preparation steps in small droplets on a fluorocarbon-coated slide. This design enables simultaneous sample preparation of thousands of single cells, including lysing, digesting, and labeling individual cells in volumes below 20 nl. We used nPOP to prepare 1,888 single cells and 128 negative controls in a single batch. Their analysis enabled quantifying the covariation between thousands of proteins and cell-cycle protein markers. Many protein sets covaried with the cell cycle similarly across all cell types and states, reflecting cell-type independent cell cycle functions. However, the cell cycle covariation of other protein sets differed markedly between cell types, even within subpopulation of melanoma cells expressing markers for drug-resistance priming. The cells expressing these markers accumulated in the G1 phase of the cell cycle and exhibited different covariation of enzymes catabolizing glucose. These results demonstrate that protein covariation across single cells may reveal functionally concerted biological differences between closely related cell states.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=59 SRC="FIGDIR/small/441211v3_ufig1.gif" ALT="Figure 1">
View larger version (53K):
org.highwire.dtl.DTLVardef@1e23b1org.highwire.dtl.DTLVardef@b93faeorg.highwire.dtl.DTLVardef@3c00forg.highwire.dtl.DTLVardef@362b41_HPS_FORMAT_FIGEXP  M_FIG C_FIG Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity Non-alcoholic fatty liver disease is a continuum of disorders among which non-alcoholic steatohepatitis (NASH) is particularly associated with a negative prognosis. Hepatocyte lipotoxicity is one of the main pathogenic factors of liver fibrosis and NASH. However, the molecular mechanisms regulating this process are poorly understood. The main aim of this study was to dissect transcriptional mechanisms regulated by lipotoxicity in hepatocytes. We achieved this aim by combining transcriptomic, proteomic and chromatin accessibility analyses from human liver and mouse hepatocytes. This integrative approach revealed several transcription factor networks deregulated by NASH and lipotoxicity. To validate these predictions, genetic deletion of the transcription factors MAFK and TCF4 was performed, resulting in hepatocytes that were better protected against saturated fatty acid oversupply. MAFK- and TCF4-regulated gene expression profiles suggest a mitigating effect against cell stress, while promoting cell survival and growth. Moreover, in the context of lipotoxicity, some MAFK and TCF4 target genes were to the corresponding differentially regulated transcripts in human liver fibrosis. Collectively, our findings comprehensively profile the transcriptional response to lipotoxicity in hepatocytes, revealing new molecular insights and providing a valuable resource for future endeavours to tackle the molecular mechanisms of NASH. Laminar Mechanisms of Saccadic Suppression in Primate Visual Cortex Saccades are a ubiquitous and crucial component of our visual system, allowing for the efficient deployment of the fovea and its accompanying neural resources. Initiation of a saccade is known to cause saccadic suppression, a temporary reduction in visual sensitivity1, 2 and visual cortical firing rates3-6. While saccadic suppression has been well characterized at the level of perception and single neurons, relatively little is known about the visual cortical networks governing this phenomenon. Here we examine the effects of saccadic suppression on distinct neural subpopulations within visual area V4. We find cortical layer- and cell type-specific differences in the magnitude and timing of peri-saccadic modulation. Neurons in the input layer show changes in firing rate and inter-neuronal correlations prior to saccade onset, indicating that this layer receives information about impending saccades. Putative inhibitory interneurons in the input layer elevate their firing rate during saccades, suggesting they play a role in suppressing the activity of other cortical subpopulations. A computational model of this circuit recapitulates our empirical observations and demonstrates that an input layer-targeting pathway can initiate saccadic suppression by enhancing local inhibitory activity. Collectively, our results provide a mechanistic understanding of how eye movement signaling interacts with cortical circuitry to enforce visual stability. Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. 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. GUIdEStaR (G-quadruplex, uORF, IRES, Epigenetics, Small RNA, Repeats), the integrated matadatabase: transcript-indexed binary information creation for chaining with neural network analysis GUIdEStaR integrates existing databases of various types of G-quadruplex, upstream Open Reading Frame (uORF), Internal Ribosome Entry Site (IRES), methylation to RNA and histone protein, small RNA, and repeats. GUIdEStaR consists of approx. 40,000 genes and 320,000 transcripts. An mRNA transcript is divided into 5 regions (5UTR, 3UTR, exon, intron, and biological region) where each region contains presence-absence data of 169 different types of elements. Recently, artificial intelligence (AI) based analysis of sequencing data has been gaining popularity in the area of bioinformatics. GUIdEStaR generates datasets that can be used as inputs to AI methods. At the GUIdEStaR homepage, users submit gene symbols by clicking a "Send" button, and shortly result files in CSV format are available for download at the result website. Users have an option to send the result files to their email addresses. Additionally, the entire database and the example Java codes are also freely available for download. Here, we demonstrate the database usage with three neural network classification studies-1) small RNA study for classifying transcription factor (TF) genes into either one of TF mediated by small RNA originated from SARS-CoV-2 or by human microRNA (miRNA), 2) cell membrane receptor study for classifying receptor genes as either with virus interaction or without one, and 3) nonsense mediated mRNA decay (NMD) study for classifying cell membrane and nuclear receptors as either NMD target or non-target. GUIdEStaR is available for access to the easy-to-use web-based database at www.guidestar.kr and for download at https://sourceforge.net/projects/guidestar. Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. Vapor exposure to Δ9-tetrahydrocannabinol (THC) slows locomotion of the Maine Lobster (Homarus americanus) RationaleDespite a long history of use in synaptic physiology, the lobster has been a neglected model for behavioral pharmacology. A restaurateur proposed that exposing lobster to cannabis smoke reduces anxiety and pain during the cooking process. It is unknown if lobster gill respiration in air would result in significant {Delta}9-tetrahydrocannabinol (THC) uptake and whether this would have any detectable behavioral effects.

ObjectiveThe primary goal was to determine tissue THC levels in the lobster after exposure to THC vapor. Secondary goals were to determine if THC vapor altered locomotor behavior or nociception.

MethodsTissue samples were collected (including muscle, brain and hemolymph) from Homarus americanus (N=3 per group) following 30 or 60 minutes of exposure to vapor generated by an e-cigarette device using THC (100 mg/mL in a propylene glycol vehicle). Separate experiments assessed locomotor behavior and hot water nociceptive responses following THC vapor exposure.

ResultsTHC vapor produced duration-related THC levels in all tissues examined. Locomotor activity was decreased (distance, speed, time-mobile) by 30 min inhalation of THC. Lobsters exhibit a temperature-dependent withdrawal response to immersion of tail, antennae or claws in warm water; this is novel evidence of thermal nociception for this species. THC exposure for 60 minutes had only marginal effect on nociception under the conditions assessed.

ConclusionsVapor exposure of lobsters, using an e-cigarette based model, produces dose-dependent THC levels in all tissues and reduces locomotor activity. Hot water nociception was temperature dependent, but only minimal anti-nociceptive effect of THC exposure was confirmed. An individualized causal framework for learning intercellular communicationnetworks that define microenvironments of individual tumors Cells within a tumor microenvironment (TME) dynamically communicate and influence each others cellular states through an intercellular communication network (ICN). In cancers, intercellular communications underlie immune evasion mechanisms of individual tumors. We developed an individualized causal analysis framework for discovering tumor-specific ICNs. Using head and neck squamous cell carcinoma (HNSCC) tumors as a testbed, we first mined single-cell RNA-sequencing data to discover gene expression modules (GEMs) that reflect the states of transcriptomic processes within tumor and stromal single cells. By deconvoluting bulk transcriptomes of HNSCC tumors profiled by The Cancer Genome Atlas (TCGA), we estimated the activation states of these transcriptomic processes in individual tumors. Finally, we applied individualized causal network learning to discover an ICN within each tumor. Our results show that cellular states of cells in TMEs are coordinated through ICNs that enable multi-way communications among epithelial, fibroblast, endothelial, and immune cells. Further analyses of individual ICNs revealed structural patterns that were shared across subsets of tumors, leading to the discovery of 4 different subtypes of networks that underlie disparate TMEs of HNSCC. Patients with distinct TMEs exhibited significantly different clinical outcomes. Our results show that the capability of estimating individual ICNs reveals heterogeneity of ICNs and sheds light on the importance of intercellular communication in impacting disease development and progression.