Identification and analysis of novel variants in SARS-COV-2 genomes isolated from the Kingdom of Bahrain The challenges imposed by the ongoing outbreak of severe acute respiratory syndrome coronavirus-2 affects every aspect of our modern world, ranging from our health to our socio-economic needs. Our existence highly depends on the vaccines availability, which demands in-depth research of the available strains and their mutations. In this work, we have analyzed all the available SERS-CoV2 genomes isolated from the Kingdom of Bahrain in terms of their variance and origin analysis. We have predicted various known and unique mutations in the SERS-CoV2 isolated from Bahrain. The complexity of the phylogenetic tree and dot plot representation of the strains mentioned above with other isolates of Asia indicates the versatility and multiple origins of Bahrains SERS-CoV2 isolates. We have also identified two high impact spike mutations from these strains which increase the virulence of SARS-CoV2. Our research could have a high impact on vaccine development and distinguishes the source of SERS-CoV2 in the Kingdom of Bahrain.
YfiB: An Outer Membrane Protein Involved in the Virulence of Shigella flexneri The intracellular pathogen Shigella flexneri, which is the causative agent of bacillary dysentery, significantly influences the worldwide implication of diarrheal infections, consequentially causing about 1.1 million deaths each year. Due to a non-availability of an authorized vaccine and the upsurge of multidrug resistance amongst Shigella strains, there has been a huge demand for further genetic analyses which could help in the advancement of new/improved drugs, and finding vaccine candidates against the pathogen. Whilst many features about the invasion of colonic cells by Shigella have been identified, fundamental gaps in information concerning in what way the bacteria transit, survive, and control gene expression, remain. The present study aims to illustrate the role of yfiB gene in Shigella virulence, which is a part of the periplasmic YfiBNR tripartite signaling system. This system is involved in the regulation of cyclic-di-GMP levels inside the bacterial cells, which is a vital messenger molecule impacting varied cellular processes such as biofilm formation, cytotoxicity, motility, synthesis of exopolysaccharide, and other virulence mechanisms like adhesion and invasion of the bacteria. Through a combination of genetic, biochemical, and virulence assays, we show how knocking out the yfiB gene can disrupt the entire YfiBNR system and affect biofilm formation, bacterial invasion, host-surface attachment, and the overall virulence of Shigella. This study eventually improves our understanding of the in-vivo persistence and survival of Shigella, brings light to the c-di-GMP lead regulation of Shigella virulence, and provides a prospective new target to design anti-infective drugs and vaccines against S. flexneri and other bacterial pathogens.
NUCLEAR RETENTION OF LEPTIN/UPD2 REGULATES ORGANISMAL RESILIENCE TO NUTRIENT EXTREMES Adipokines released from the adipocytes function as a systemic adipometer; they impinge on neural circuits to signal nutrient status. On starvation, adipokines must be retained to signal energy deficit; else, it significantly reduces starvation survival. But how fat cells retain adipokines is unclear. Here, we demonstrate that Atg8, a cell-intrinsic autophagy factor, regulates the starvation-induced acute retention of the Leptin Drosophila ortholog Upd2. We show that on starvation, as a direct consequence of Atg8s lipidation, Upd2 accumulates in the nucleus. We illustrate that Upd2s nuclear retention is critical to fat mobilization and increased starvation resilience. Furthermore, nuclear Upd2 promotes the expression of a secreted innate immune gene signature. This hints at an unanticipated connection between adipokine nuclear retention and increased innate immunity. In conclusion, we propose that, during starvation, Atg8s role is not just limited to autophagy but is critical for withholding adipokines in the nucleus to promote starvation resilience.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=178 HEIGHT=200 SRC="FIGDIR/small/428913v5_ufig1.gif" ALT="Figure 1">
View larger version (33K):
org.highwire.dtl.DTLVardef@1826425org.highwire.dtl.DTLVardef@17cca52org.highwire.dtl.DTLVardef@1c70ac7org.highwire.dtl.DTLVardef@94a024_HPS_FORMAT_FIGEXP  M_FIG C_FIG O_LIIn fed state Upd2 requires Atg8 for nuclear exit and cytosolic localization.
C_LIO_LIAtg8s lipidation on starvation results in Upd2s nuclear accumulation.
C_LIO_LIUpd2 nuclear retention on starvation increases fat mobilization and post-starvation hunger.
C_LIO_LIOn starvation Upd2 nuclear retention increases expression of a secreted innate immune signature.
C_LI CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma. Beneficial effects of novel strains of Aureobasidium pullulans produced 1,3-1,6 β-glucans on non-esterified fatty acid levels in diabetic KKAy mice Obesity, metabolic syndrome, associated lipotoxicity and its cascade of events contribute to the majority of the burden related to non-communicable diseases globally. Preventive lifestyle changes aside, several beneficial effects have been reported in type II diabetes mellitus and dyslipidaemia patients with biological response modifier glucans (BRMG) produced as an exopolysaccharide by Aureobasidium pullulans. In this study, we compared two strains (AFO-202 and N-163) that produce beta glucans in alleviating lipotoxicity. This study was performed in obese diabetic mice model of KK-Ay mice, in four groups with six subjects in each group - Group 1: sacrificed on Day 0 for baseline values; Group 2: control (drinking water); Group 3: AFO-202 beta glucan--200 mg/kg/day; Group 4: N-163 beta glucan--300 mg/kg/day. The animals in groups 2-4 had the test solutions administered by gavage once daily for 28 consecutive days. Biochemical analyses were conducted of blood glucose, triglycerides, total cholesterol, LDL cholesterol, HDL cholesterol and non-esterified fatty acids (NEFA). Group 4 (N-163) had the lowest NEFA levels, as compared to the other groups, and marginally decreased triglyceride levels. The groups had no significant differences in blood glucose, HbA1c, triglycerides, or LDL and HDL cholesterol. N-163 produced by A. pullulans decreased NEFA in a diabetic mice model in 28 days. These results, although modest, warrant further in-depth research into lipotoxicity and associated inflammatory cascades in both healthy and disease affected subjects to develop novel strategies for prevention and management. Embryo development in Carica papaya Linn. Papaya (Carica papaya Linn.) is a tropical plant whose draft genome has been sequenced. Papaya produces large fruits rich in vitamins A and C and is an important cash crop in developing countries. Nonetheless, little is known about how the female gametophyte develops, how it is fertilized and how it develops into a mature seed containing an embryo and an endosperm. The Papaya female gametophyte displays a Polygonum-type architecture consisting of two synergid cells, an egg cell, a central cell, and three antipodal cells. Reports are available of the presumed existence of varieties in which cross fertilization is bypassed and autonomous development of embryos occurs (e.g., apomixis). In this study, we analyzed the development of female gametophytes in a commercial Hawaiian parental line and in the presumed apomictic Costa Rican line L1. Samples were collected before and after anthesis to compare the overall structure, size and transcriptional patterns of several genes that may be involved in egg and endosperm cell fate and proliferation. These genes were the putative papaya homologs of ARGONAUTE9 (AGO9), MEDEA (MEA), RETINOBLASTOMA RELATED-1 (RBR1), and SLOW WALKER-1 (SWA1). Our results suggest that its feasibly to identify structural features of Polygonum-type development, and that in bagged female flowers of line L1 we might have observed autonomous development of embryo-like structures. Possible downregulation of papaya homologs for AGO9, MEA, RBR1 and SWA1 was observed in embryo sacs from line L1 before and after anthesis, which may suggest a link between suspected apomixis and transcriptional downregulation of genes for RNA-directed DNA methylation, histone remodelers, and rRNA processing. Most notably, the large size of the papaya embryo sac suggests that it could be a cytological alternative to Arabidopsis thaliana for study. Significant variation in embryo sac size was observed between the varieties under study, suggesting wide differences in the genetic regulation of anatomical features. An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis 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. Pore dynamics and asymmetric cargo loading in an encapsulin nanocompartment Encapsulins are protein nanocompartments that house various cargo enzymes, including a family of decameric ferritin-like proteins. Here, we study a recombinant Haliangium ochraceum encapsulin:encapsulated ferritin complex using electron cryo-microscopy and hydrogen/deuterium exchange mass spectrometry to gain insight into the structural relationship between the encapsulin shell and its protein cargo. An asymmetric single particle reconstruction reveals four encapsulated ferritin decamers in a tetrahedral arrangement within the encapsulin nanocompartment. This leads to a symmetry mismatch between the protein cargo and the icosahedral encapsulin shell. The encapsulated ferritin decamers are offset from the interior face of the encapsulin shell. Using HDX-MS, we observed dynamic behavior of the major five-fold pore in the encapsulin shell and show the pore opening via the movement of the encapsulin A-domain. These data will accelerate efforts to engineer the encapsulation of heterologous cargo proteins and to alter the permeability of the encapsulin shell via pore modifications.

TeaserCryo-EM and HDX-MS analysis of an encapsulin nanocompartment shows that the pores at the five-fold icosahedral vertex of the shell are flexible. The multi-level phenotypic impact of synonymous substitutions: heterologous gene expression in human cells. Redundancy in the genetic code allows for differences in transcription and/or translation efficiency between sequences carrying synonymous polymorphisms, potentially leading to phenotypic changes. It is commonly admitted that the evolution codon usage bias ("CUB", the over-representation of certain codons in a genome, a gene or in positions along a gene) are driven by a combination of neutral and selective processes, but their relative contribution is a matter of debate, especially in mammals. Particularly, integrative studies quantifying the phenotypic impact of CUB at different molecular and cellular levels are lacking. Here we report a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells.

Six synonymous versions of the shble antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was assessed by means of: i) mRNA-to-DNA and protein-to-mRNA ratios for each shble version; ii) cellular fluorescence, using flow cytometry, as a proxy for single-cell level expression; and iii) real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness.

We show that differences in CUB strongly impact the molecular and the cellular phenotype: i) they result in large differences in mRNA and in protein levels, as well in mRNA-to-protein ratio; ii) they introduce splicing events not predicted by current algorithms; iii) they lead to reproducible phenotypic heterogeneity; iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression.

We interpret that CUB modulate mRNA availability and suitability for translation in human cells, leading to differences in protein levels and eventually eliciting phenotypic differences.

AUTHOR SUMMARYThe genetic code is redundant, with several codons encoding for the same amino acid. These synonymous codons are not used with equal frequencies. Instead, codon usage bias (CUB) varies between species, genes, and even positions along a gene. At each of these levels, CUB are shaped by the overall balance between mutational biases and selection forces. To shed light on the molecular mechanisms underlying molecular and organism phenotypes, integrative studies quantifying the phenotypic impact of CUB at different levels of biological integration are necessary. Here, we monitored the multiscale changes induced by six synonymous versions of an antibiotic resistance gene independently expressed in a human cell line. We show that: 1. both mRNA levels, protein levels, and their ratios are affected by CUB; 2. potential effects on mRNA (splicing and availability for translation) seem to constitute the main level of action; 3. cell fitness is severely impacted by a trade-off between the burden of heterologous expression and the benefits of antibiotic resistance. This integrative study provides new insights on translation regulation and the associated phenotypic impact in human cells, associated to CU. Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted. Interacting phenotypes and the coevolutionary process: Interspecific indirect genetic effects alter coevolutionary dynamics Coevolution occurs when species interact to influence one anothers fitness, resulting in reciprocal evolutionary change. In many coevolving lineages, trait expression in one species is modified by the genotypes and phenotypes of the other, forming feedback loops reminiscent of models of intraspecific social evolution. Here, we adapt the theory of within-species social evolution, characterized by indirect genetic effects and social selection imposed by interacting individuals, to the case of interspecific interactions. In a trait-based model, we derive general expressions for multivariate evolutionary change in two species and the expected between-species covariance in evolutionary change when selection varies across space. We show that reciprocal interspecific indirect genetic effects can dominate the coevolutionary process and drive patterns of correlated evolution beyond what is expected from direct selection alone. In extreme cases, interspecific indirect genetic effects can lead to coevolution when selection does not covary between species or even when one species lacks genetic variance. Moreover, our model indicates that interspecific indirect genetic effects may interact in complex ways with cross-species selection to determine the course of coevolution. Importantly, our model makes empirically testable predictions for how different forms of reciprocal interactions contribute to the coevolutionary process. Yap/Taz-Activated Tert-Expressing Acinar Cells Are Required for Pancreatic Regeneration The expression of TERT (telomerase reverse transcriptase) has been implicated in stem and progenitor cells, which are essential for tissue homeostasis and regeneration. However, the roles of TERT-expressing cells in the pancreas remain elusive. Employing genetically engineered Tert knock-in mouse model, herein, we located a rare population of Tert+ acinar cells. While Tert+ cells are quiescent in normal conditions, acinar cell injury leads to mitotic activation of Tert+ cells and subsequent generation of new acinar cells. Moreover, the genetic ablation of Tert+ cells impairs pancreatic regeneration. We further found that Yap/Taz activation is required for the expansion of Tert+ acinar cells. Our results identified Tert+ acinar cells as a distinct subset of acinar cells, which contributes to pancreatic regeneration via Yap/Taz activation. Nuclear alpha-synuclein is present in the human brain and is modified in dementia with Lewy bodies Dementia with Lewy bodies is pathologically defined by the cytoplasmic accumulation of alpha-synuclein within neuronal cells in the brain. Alpha-synuclein is predominately pre-synaptic, but has been reported present in various subcellular compartments in cell and animal models. In particular, nuclear alpha-synuclein is evident in-vitro and in disease models and has been associated with altered DNA integrity, gene transcription, nuclear homeostasis. However, owing to various factors, the presence of alpha-synuclein in the nuclei of human brain cells remains controversial, as does its role in synucleinopathies.

Here, we close this gap and provide a unique demonstration confirming the presence of nuclear alpha-synuclein in post-mortem brain tissue obtained from cases of dementia with Lewy bodies as well as from controls via immunohistochemistry, immunoblot, and label-free mass-spectrometry.

Discrete intra-nuclear alpha-synuclein puncta reactive against phosphorylated serine 129-alpha-synuclein and pan-alpha-synuclein antibodies were observed in cortical neurons and non-neuronal cells in fixed brain sections and in isolated nuclear preparations from Dementia with Lewy bodies cases and matched controls. Subsequent biochemical analysis of subcellular fractionated tissue confirmed alpha-synuclein as present in a nuclear fraction at levels ~ 10-fold lower than in the cytoplasm. Critically, however, an increase in monomeric nuclear alpha-synuclein phosphorylated as serine 129 was observed in cases of dementia with Lewy bodies alongside higher molecular weight pan- and phosphorylation reactive alpha-synuclein species, consistent with the formation of intranuclear phosphorylated alpha-synuclein oligomers. Furthermore, the presence of nuclear alpha-synuclein was confirmed via label free mass spectrometry, as 6 unique alpha-synuclein derived peptide sequences were identified in nuclear fractions (71.4% sequence coverage).

Collectively, our data confirm the presence of nuclear alpha-synuclein in human brain tissue and describe nuclear pathology associated with dementia with Lewy bodies. These findings address a major controversy in the synucleinopathy field by confirming the presence of nuclear alpha-synuclein in autoptic human brain tissue and, for the first time, identify that alpha-synuclein is aggregated into novel and potentially pathological assemblies in the nucleus as part of the disease process associated with dementia with Lewy bodies and thus may contribute to the disease phenotype. Oxytocin promotes social grooming in bonobos Oxytocin has attracted research attention due to its role in promoting social bonding. One notable recent hypothesis is the biobehavioral feedback loop, which posits that the oxytocin system has evolved to support the formation and maintenance of social bonds through a positive feedback loop, where oxytocin promotes social behaviours which then cause oxytocin release themselves. In the two Pan species, humans closest relatives, oxytocin is known to be released following key behaviours related to social bonding, such as social grooming in chimpanzees and female-female sexual behaviour in bonobos. However, no experimental evidence has demonstrated that oxytocin promotes such socio-positive behaviours. To test this, we administered nebulized oxytocin or saline placebo to a group of female bonobos and subsequently observed the change in their gross behavior during free interaction. We found that bonobos groomed other group members significantly more frequently in the oxytocin compared to placebo condition. Other behavioural measures did not largely differ between conditions, except for a nonsignificant trend for reduction in abnormal regurgitation/reingestion behaviour. Overall, we found that oxytocin promoted socio-positive interaction in bonobos, providing support for the biobehavioural feedback loop hypothesis of oxytocin in bonobo social evolution. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

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

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. A multi-dataset evaluation of frame censoring for motion correction in task-based fMRI Subject motion during fMRI can affect our ability to accurately measure signals of interest. In recent years, frame censoring--that is, statistically excluding motion-contaminated data within the general linear model using nuisance regressors--has appeared in several task-based fMRI studies as a mitigation strategy. However, there have been few systematic investigations quantifying its efficacy. In the present study, we compared the performance of frame censoring to several other common motion correction approaches for task-based fMRI using open data and reproducible workflows. We analyzed eight publicly-available datasets representing eleven distinct tasks in child, adolescent, and adult participants. Performance was quantified using maximum t-values in group analyses, and ROI-based mean activation and split-half reliability in single subjects. We compared frame censoring across several thresholds to the use of 6 and 24 canonical motion regressors, wavelet despiking, robust weighted least squares, and untrained ICA-based denoising, for a total of 240 separate analyses. Thresholds used to identify censored frames were based on both motion estimates (FD) and image intensity changes (DVARS). Relative to standard motion regressors, we found consistent improvements for modest amounts of frame censoring (e.g., 1-2% data loss), although these gains were frequently comparable to what could be achieved using other techniques. Importantly, no single approach consistently outperformed the others across all datasets and tasks. These findings suggest that the choice of a motion mitigation strategy depends on both the dataset and the outcome metric of interest. G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from {tau}RAMD Simulations There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of timescales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the {tau}-Random Acceleration Molecular Dynamics ({tau}RAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the {beta}2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

TABLE OF CONTENTS GRAPHIC

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=100 SRC="FIGDIR/small/449151v2_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@743dcdorg.highwire.dtl.DTLVardef@1857381org.highwire.dtl.DTLVardef@62359aorg.highwire.dtl.DTLVardef@1ad38a3_HPS_FORMAT_FIGEXP  M_FIG C_FIG A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Sex differences and risk factors for bleeding in Alagille syndrome Spontaneous lethal bleeds are major cause of death in the pediatric liver disease Alagille syndrome (ALGS), yet risk factors and screening methods have not been established. We performed a systematic review and identified significantly more female than male patients with idiopathic intracranial hemorrhage (10:1). We investigated bleeding and vasculature in patients and a mouse model for ALGS (Jag1Ndr/Ndr mice) and asked whether phenotypes identified in mice could be detected in patients non-invasively. Jag1Ndr/Ndr mice bled spontaneously, exhibiting a thin skull and vascular defects including artery-vein crossings, tortuous vessels, capillary breakdown and CADASIL-like sparse vascular smooth muscle cell coverage which was aggravated by hypertension. Retinographs from patients confirmed tortuous blood vessels and artery-vein crossings in ALGS. In conclusion, Jag1Ndr/Ndr mice could be used to develop interventions for vascular defects in ALGS, and retinography could provide a non-invasive method for vascular analysis in these pediatric patients. T Cell Receptor Beta Germline Variability is Revealed by Inference From Repertoire Data BackgroundT and B cell receptor (TCR, BCR) repertoires constitute the foundation of adaptive immunity. Adaptive immune receptor repertoire sequencing (AIRR-seq) is a common approach to study immune system dynamics. Understanding the genetic factors influencing the composition and dynamics of these repertoires is of major scientific and clinical importance. The chromosomal loci encoding for the variable regions of TCRs and BCRs are challenging to decipher due to repetitive elements and undocumented structural variants.

MethodsTo confront this challenge, AIRR-seq-based methods have recently been developed for B cells, enabling genotype and haplotype inference and discovery of undocumented alleles. However, this approach relies on complete coverage of the receptors variable regions, whereas most T cell studies sequence a small fraction of that region. Here, we adapted a B cell pipeline for undocumented alleles, genotype, and haplotype inference for full and partial TCR sequences. The pipeline also deals with gene assignment ambiguities, which is especially important in the analysis of data-sets of partial sequences.

ResultsWe identified 39 undocumented polymorphisms in T cell receptor Beta V (TRBV) and 31 undocumented 5 UTR sequences. A subset of these inferences was also observed using independent genomic approaches. We found that a single nucleotide polymorphism differentiating between the two documented T cell receptor Beta D2 (TRBD2) alleles is strongly associated with dramatic changes in the expressed repertoire.

ConclusionsWe reveal a rich picture of germline variability, and demonstrate how a single nucleotide polymorphism dramatically affects the composition of the whole repertoire. Our findings provide a basis for annotation of TCR repertoires for future basic and clinical studies. Classification, Using Field Potential Images, of Pacemaker Micro-coordination in the Mouse Intestine The flexible and sophisticated movement of the gastrointestinal (GI) tract implies the involvement of mechanisms other than enteric neural circuits, to coordinate excitation in microregions. We thus performed microimaging of pacemaker dynamics in the small intestine of mice since it contains typical network-forming pacemaker cells. A dialysis membrane-reinforced low-impedance microelectrode array (MEA) enabled field potentials over a wide frequency range to be stably measured in microregions. The pacemaker dynamics were classified into basic patterns despite large variations. In the developmental process, pacemaker activity was categorized as either an  expanding or a  migrating pattern that was initiated in or propagated to the MEA sensing area, respectively. The intercellular current of the volume conductor complicated the waveform of both activities. The existence of  expanding and  migrating patterns was attributable to duplicated pacemaker systems such as intracellular Ca2+ oscillation-activated and voltage-gated mechanisms. Additionally, from the spatio-temporal feature during the period of pacemaker events, the  bumpy/aberrant pattern was defined by aberrant, incoherent propagation, and associated with local impairment of excitability, while the  colliding/converging pattern involved the interaction of multiple activities in the MEA area. Interconversion between the four micro-coordination patterns occurred in the same microregion. 5-Hydroxytryptamine (5-HT) promoted  migrating activity, implying an improvement or restoration of spatial conductivity. These results agree well with the action of 5-HT to change GI movement toward propulsion. In conclusion, our MEA method of microimaging classification enables the quantitative assessment of spatio-temporal electric coordination underlying GI motility, suggesting its application to small model animals. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI 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. Impedance scaling for gold and platinum microelectrodes ObjectiveElectrical measurement of the activity of individual neurons is a primary goal for many invasive neural electrodes. Making these "single unit" measurements requires that we fabricate electrodes small enough so that only a few neurons contribute to the signal, but not so small that the impedance of the electrode creates overwhelming noise or signal attenuation. Thus, neuroelectrode design often must strike a balance between electrode size and electrode impedance, where the impedance is often assumed to scale linearly with electrode area.

Approach and main resultsHere we study how impedance scales with neural electrode area and find that the 1 kHz impedance of Pt electrodes (but not Au electrodes) transitions from scaling with area (r-2) to scaling with perimeter (r-1) when the electrode radius falls below 10 microns. This effect can be explained by the transition from planar to spherical diffusion behavior previously reported for electrochemical microelectrodes. Significance: These results provide important intuition for designing small, single unit recording electrodes. Specifically, for materials where the impedance is dominated by a pseudo-capacitance that is associated with a diffusion limited process, the total impedance will scale with perimeter rather than area when the electrode size becomes comparable with the diffusion layer thickness. For Pt electrodes this transition occurs around 10 um radius electrodes. At even lower frequencies (1 Hz) impedance approaches a constant. This transition to r-1 scaling implies that electrodes with a pseudo-capacitance can be made smaller than one might expect before thermal noise or voltage division limits the ability to acquire high-quality single-unit recordings. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Simpati: patient classifier identifies signature pathways based on similarity networks for the disease prediction BACKGROUNDPathway-based patient classification is a supervised learning task which implies a model learning pathways as features to predict the classes of patients. The counterpart of enrichment tools for the pathway analysis are fundamental methods for clinicians and biomedical scientists. They allow to find signature cellular functions which help to define and annotate a disease phenotype. They provide results which lead human experts to manually classify patients. It is a paradox that pathwaybased classifiers which natively resolve this objective are not strongly developed. They could simulate the human way of thinking, decipher hidden multivariate relationships between the deregulated pathways and the disease phenotype, and provide more information than a probability value. Instead, there are currently only two classifiers of such kind, they require a nontrivial hyperparameter tuning, are difficult to interpret and lack in providing new insights. There is the need of new classifiers which can provide novel perspectives about pathways, be easy to apply with different biological omics and produce new data enabling a further analysis of the patients.

RESULTSWe propose Simpati, an innovative and interpretable patient classifier based on pathway-specific patient similarity networks. The first classifier to adopt ad-hoc novel algorithms for such graph type. It standardizes the biological high-throughput dataset of patients profiles with a propagation algorithm that considers the interconnected nature of the cells molecules for inferring a new activity score. This allows Simpati to classify with dense, sparse, and non-homogenous omic data. Simpati organizes patients molecules in pathways represented by patient similarity networks for being interpretable, handling missing data and preserving the patient privacy. A network represents patients as nodes and a novel similarity measure determines how much every pair act co-ordinately in a pathway. Simpati detects signature biological processes based on how much the topological properties of the related networks separate the patient classes. In this step, it includes a new cohesive subgroup detection algorithm to handle patients not showing the same pathway activity as the other class members. An unknown patient is then classified by a unique recommender system which considers how much is similar to known patients and distant from being an outlier. Simpati outperforms previously published classifiers on five cancer datasets described with two biological omics, classifies well with sparse data, identifies more relevant pathways associated to the patients disease than the competitors and has the lowest computational requirements.

CONCLUSIONSimpati can serve as generic-purpose pathway-based classifier of patient classes. It provides signature pathways to unveil the altered biological mechanisms of a disease phenotype and to classify patients according to the learnt pathway-specific similarities. The signature condition and patient prediction can be deciphered considering the patient similarity networks which must reveal the members of a patient class more cohesive and similar than the non-members. Simpati divides the pathways in up and downinvolved. Upinvolved when the signaling cascades generated by the altered molecules of the disease patients impact stronger the pathway than the ones of the control class. We provide an R implementation, a graphical user interface and a visualization function for the patient similarity networks. The software is available at: https://github.com/LucaGiudice/Simpati