From smooth cortical gradients to discrete modules: A biologically plausible mechanism for the self-organization of modularity in grid cells The grid cell system is a paradigmatic example of the computational advantages of modular representations. Modular grid responses emerge along a strip of cortex with smooth gradients in several biophysical properties, within days of eye opening in juvenile animals. The mechanisms underlying this emergence are unknown. We show that multiple spatially and functionally discrete grid modules can spontaneously self-organize from lateral interactions that are continuously graded in scale, when combined with a fixed-scale interaction. We derive a comprehensive analytic theory to reveal that modularization through this mechanism is highly generic and robust. It predicts that each term in the sequence of grid period ratios should be distinct and furnishes values that are invariant to microscopic details about dynamics and connectivity, resulting in the most accurate match to data to date. Altogether, this work reveals novel self-organization mechanisms by which simple local interactions and smooth gradients may interact to produce macroscopic modular organization.
Analyses of infant vocal communication in mouse models of autism spectrum disorder using a novel unsupervised machine learning approach Infant crying is an innate communicative behavior that is frequently impaired in certain neurodevelopmental disorders (NDDs). Since advanced paternal age is a reported risk factor for NDDs in offspring, we evaluated the impact of a fathers age on early vocal development in C57BL/6J mice. We recorded and applied a unique combination of computational analyses to ultrasonic vocalizations (USVs) emitted by mouse pups sired by young and aged fathers. Our data showed that advanced paternal age reduced the number and duration of USVs, and altered the syllable composition in pups. Moreover, pups born to young fathers showed convergent vocal characteristics with a rich repertoire during postnatal development, while those born to aged fathers exhibited more divergent vocal patterns with limited repertoire. Principal component analysis in conjunction with clustering analysis demonstrated that pups from aged fathers deviated from typical trajectories of vocal development, which were considered as atypical individuals. Thus, our study indicates that advanced paternal age has a significant effect on offsprings early vocal development. It is suggested that the trajectories of vocal development could be a useful marker of the NDD-like phenotype associated with the advanced paternal age. In addition, our comprehensive computational analysis described here is an effective approach to characterize the altered individual diversity relevant to neurodevelopmental disorders.

One Sentence SummaryAdvanced paternal age affects vocal development in early postnatal mice, with more pups showing atypical developmental trajectories.
METAbolomics data Balancing with Over-sampling Algorithms (META-BOA): an online resource for addressing class imbalance MotivationClass imbalance, or unequal sample sizes between classes, is an increasing concern in machine learning for metabolomic and lipidomic data mining, which can result in overfitting for the over-represented class. Numerous methods have been developed for handling class imbalance, but they are not readily accessible to users with limited computational experience. Moreover, there is no resource that enables users to easily evaluate the effect of different over-sampling algorithms.

ResultsMETAbolomics data Balancing with Over-sampling Algorithms (META-BOA) is a web-based application that enables users to select between four different methods for class balancing, followed by data visualization and classification of the sample to observe the augmentation effects. META-BOA outputs a newly balanced dataset, generating additional samples in the minority class, according to the users choice of Synthetic Minority Over-sampling Technique (SMOTE), Borderline-SMOTE (BSMOTE), Adaptive Synthetic (ADASYN), or Random Over-Sampling Examples (ROSE). META-BOA further displays both principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) visualization of data pre- and post-over-sampling. Random forest classification is utilized to compare sample classification in both the original and balanced datasets, enabling users to select the most appropriate method for their analyses.

Availability and implementationMETA-BOA is available at https://complimet.ca/meta-boa.

Supplementary InformationSupplementary material is available at Bioinformatics online. From fitting the average to fitting the individual: A cautionary tale for mathematicalmodelers An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Author summaryAs we enter the era of healthcare where personalized medicine becomes a more common approach to treating cancer patients, harnessing the power of mathematical models will only become more essential. Using a preclinical dataset on cancer immunotherapy, we explore the challenges and limitations that arise when trying to move from fitting and making predictions for the population-level average, to fitting and making predictions for an individual. We find that the standard of approach of picking a single fitting methodology and a single cost function may end up having limited predictive value when applied to individual data. We also show how having a small amount of the right additional experimental data, and establishing the robustness of average treatment response, can help improve confidence in personalized model predictions. Emergence of a novel mobile RND-type efflux pump gene cluster, tmexC3D2-toprJ1b, in Pseudomonas species Tigecycline exhibits promising activity against multidrug-resistant Gram-negative bacteria (MDR-GNB), whose infections are difficult to treat with antimicrobials. However, mobile tigecycline resistance genes, such as tmexCD-toprJ, have emerged in Enterobacterales isolated in China, Vietnam, and possibly other countries in the world. In this study, we investigated tigecycline-nonsusceptible GNB in Japan. Eight tigecycline- and carbapenem-nonsusceptible isolates of Pseudomonas alcaligenes were obtained from sewage water from a medical institution in Japan in 2020. Whole genome analysis of all P. alcaligenes isolates was performed using short-read sequencing, and the isolate KAM426 was further analyzed using long-read sequencing. For important antimicrobial resistance genes, analysis of surrounding structures and comparison with similar sequences in the public genome database were performed. We identified a novel hybrid type of tmexCD-toprJ gene cluster, tmexC3D2-toprJ1b consisting of tmexC3, tmexC2, and toprJ1b, in phylogenetically clonal isolates of P. alcaligenes. The complete genome sequence of KAM426 revealed that this isolate co-harbors tmexC3D2-toprJ1b and two copies of the carbapenemase gene blaIMP-1 on the chromosome. tmexC3D2-toprJ1b in KAM426 was flanked by the IS5/IS1182 family transposase gene, suggesting that the gene cluster was acquired by horizontal gene transfer (HGT). tmexC3D2-toprJ1b seemed to have spread to other Pseudomonas species such as Pseudomonas aeruginosa via HGT mediated by mobile gene elements such as a plasmid. This study identified tmexCD-toprJ-like tigecycline resistance genes in Japan for the first time and suggests that diverse tmexCD-toprJ-like gene clusters, including tmexC3D2-toprJ1b, have spread among MDR-GNB worldwide. Further epidemiological genomic studies in clinical and environmental settings are needed. Identification and characterization of the WYL BrxR protein and its gene asseparable regulatory elements of a BREX phage restriction system Bacteriophage exclusion ( BREX) phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR ( BREX Regulator), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that the BrxR coding sequence is involved in cis regulation of anti-phage activity, while the BrxR protein plays an additional regulatory role, perhaps during horizontal transfer. Muscle Activation Profile During Perturbed Walking is Modulated According to Body State During an unexpected loss of balance, avoiding a fall requires people to readjust their footing rapidly and effectively. We investigated the activation patterns of the ankle and knee muscles, and muscle fiber type recruitment resulting from unannounced, medio-lateral (i.e., right/left) horizontal-surface walking perturbations in twenty healthy adults (27.00{+/-}2.79 years, 10 females). Surface electromyography (sEMG) total spectral power for specific frequency bands (40-60Hz, 60-150Hz, 150-250Hz, 250-400Hz and 400-1000Hz), from tibialis anterior (TA) and vastus lateralis (VL) muscles were analyzed. Compared to non-perturbed walking, we found a significant increase in the total spectral power of lower-extremity muscles during the first 3 seconds after perturbation. When two feet were on the ground in time of perturbation we found a different muscle fiber type recruitment pattern between VL and TA muscles. This was not significant for perturbations implemented when one foot was on the ground. Our findings suggest that muscle operating frequency is modulated in real time to fit body state and somatosensory input in the context of functional goal requirements such as a rapid change of footing in response to unexpected loss of balance in single and double-support phases of gait.

New & NoteworthyTo study muscle spectral profiles in response to loss of balance, we investigated the dynamics of muscle spectral power changes, across different frequency bands after unannounced mechanical perturbations during walking. We showed increased activation of high-frequency motor units of the lower-limb muscles, subside 3 seconds after perturbation. Differences in power increase of specific frequency bands suggest that muscle activation is modulated in real time to fit body state in the context of functional goal requirements. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Early postnatal activation of the hypoxia pathway disrupts β-cell function Hypoxic insults in the perinatal period can lead to persistent hyperinsulinism and profound hypoglycemia in newborns. We studied the impact of the hypoxia-inducible factor 1A (HIF1A) pathway on postnatal {beta}-cell function. Rat pups were treated daily between postnatal day (P)7 to P10 with adaptaquin (AQ), an inhibitor of prolyl hydroxylases, which stabilizes HIF1A. AQ-treated pups were hypoglycemic and had higher plasma insulin concentrations. Their islets had a decreased glucose threshold for insulin secretion, indicative of a delay in {beta}-cell postnatal functional maturation. Histology analyses revealed that AQ-treated pups had increased pancreatic insulin-positive area but no changes in the number of islets or number of {beta}-cells per islet, suggesting larger average {beta}-cell size. AQ-treated rat pups had decreased expression of cell cycle genes and decreased numbers of proliferating {beta}-cells. In conclusion, pharmacologic activation of the HIF1A pathway in the early postnatal period leads to hyperinsulinism, due to the persistence of a low glucose threshold for insulin secretion, and to decreased early postnatal {beta}-cell proliferation, suggesting it can impact adult {beta}-cell mass and diabetes risk. Native American Ancestry and Pigmentation Allele Contributions to Skin Color in a Caribbean Population Interest in the genetic basis of variation in skin pigmentation in Native American populations led us to seek indigenous populations of the Western Hemisphere with African and minimal European admixture to study the effect of Native American ancestry on skin color. Admixture analysis from DNA collected from 458 individuals in the Kalinago territory of the Commonwealth of Dominica showed shared ancestry with East Asians at K=3 and 55% Native American, 32% African, and 11% European ancestry at K=6, the highest Native American ancestry of Caribbean populations. Skin pigmentation was 20 to 80 melanin units, averaging 46. Three albino individuals were homozygous for multi-nucleotide polymorphism OCA2NW273KV of African origin, whose population allele frequency was 0.03 and single allele effect size was -8 melanin units. Hypopigmenting allele frequencies for SLC24A5A111T and SLC45A2L374F were 0.14 and 0.05, whose single allele effect sizes were -6 and -3, respectively. Skin color plots of individuals lacking known hypopigmenting alleles suggests that Native American Ancestry reduced pigmentation by more than 20 melanin units (low and high estimates 21.8 and 28.5). Shared ancestry with East Asians at K=3 suggests potential sharing of one or more pigmentation alleles. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Conformity leading self-organised rescue from a suboptimal decision-making bias Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Cellular feedback to organotelluranes displays modulation of antioxidant proteins gene expression Organotelluranes RT3 and RT4 are thiol reagents that induce mitochondrial transition pore (MTP) opening in a sensitive and insensitive manner to cyclosporin A. Although RT3 and RT4 promote glutathione depletion, paradoxically, they are also an efficient antioxidant for membrane lipids. These compounds antagonistic effects elicited the challenging question of how the gene expression of antioxidant enzymes would respond to treatment with these compounds. The influence of RT3 and RT4 on antioxidant enzyme expression was investigated in cultured aortic smooth muscle cells (ASMC). RT3 and RT4 promoted disruption of ionic calcium homeostasis, mitochondrial transmembrane potential ({Delta}{Psi}), and cell death in a dose-dependent manner. The cell death mechanisms responded qualitatively to the increase of the organotellurane concentration and changed from apoptosis to necrosis. RT3 and RT4 increased the expression of thioredoxin significantly. RT3 also increased the expression of glutaredoxin and glutathione peroxidase, slightly the catalase expression without significant effects on SOD expression. The results are consistent with GSH and protein thiol depletion and discussed based on the cell toxicity mechanism exhibited by these compounds.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/425411v4_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@52bb82org.highwire.dtl.DTLVardef@19476d9org.highwire.dtl.DTLVardef@6aa99org.highwire.dtl.DTLVardef@1aa605e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Computational modeling and quantitative cell physiology reveal central parameters for the brassinosteroid-regulated cell growth of the Arabidopsis root Brassinosteroids (BR) are key hormonal regulators of plant development. However, whereas the individual components of BR perception and signaling are well characterized experimentally, the question of how they can act and whether they are sufficient to carry out the critical function of cellular elongation remains open. Here, we combined computational modeling with quantitative cell physiology to understand the dynamics of the plasma membrane (PM)-localized BR response pathway during the initiation of cell elongation in the epidermis of the Arabidopsis root tip. The model, consisting of ordinary differential equations, comprises the BR induced hyperpolarization of the PM, the acidification of the apoplast and subsequent cell wall swelling. We demonstrated that the competence of the root epidermal cells for the BR response predominantly depends on the amount and activity of H+-ATPases in the PM. The model further predicted that an influx of cations is required to compensate for the shift of positive charges caused by the apoplastic acidification. A potassium channel was subsequently identified and experimentally characterized, fulfilling this function. Thus, we established the landscape of components and parameters capable of triggering and guiding cellular elongation through the fast response to BRs, a central process in plant development. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Omicron booster in ancestral strain vaccinated mice augments protective immunities against both the Delta and Omicron variants A booster vaccination is called for constraining the evolving epidemic of SARS-CoV-2. However, the necessity of a new COVID-19 vaccine is currently unclear. To compare the effect of an Omicron-matched S DNA vaccine and an ancestral S DNA vaccine in boosting cross-reactive immunities, we firstly immunized mice with two-dose of a DNA vaccine encoding the spike protein of the ancestral Wuhan strain. Then the mice were boosted with DNA vaccines encoding spike proteins of either the Wuhan strain or the Omicron variant. Specific antibody and T cell responses were measured at 4 weeks post boost. Our data showed that the Omicron-matched vaccine efficiently boosted RBD binding antibody and neutralizing antibody responses against both the Delta and the Omicron variants. Of note, antibody responses against the Omicron variant elicited by the Omicron-matched vaccine were much stronger than those induced by the ancestral S DNA vaccine. Meanwhile, CD8+ T cell responses against both the ancestral Wuhan strain and the Omicron strain also tended to be higher in mice boosted by the Omicron-matched vaccine than those in mice boosted with the ancestral S DNA vaccine, albeit no significant difference was observed. Our findings suggest that an Omicron-matched vaccine is preferred for boosting cross-reactive immunities. A liquid +TIP-network regulates microtubule dynamics through tubulin condensation Regulation of microtubule dynamics is essential for diverse cellular functions, and proteins that bind to dynamic microtubule ends can regulate network dynamics. Here we show, that two conserved microtubule end-binding proteins, CLIP-170 and EB3, undergo phase separation and form dense liquid-networks. When CLIP-170 and EB3 act together the multivalency of the network increases, which synergistically increases the amount of protein in the dense phase. In vitro and in cells these liquid networks can condense tubulin. In vitro in the presence of microtubules, EB3/CLIP-170 phase separation can co-condense tubulin all along the microtubule. At this condition microtubule growth speed increases up to two-fold and depolymerization events are strongly reduced, compared to conditions with phase separation deficient networks. Our data show that phase separated EB3/CLIP-170 networks impact microtubule growth dynamics beyond direct protein-microtubule interactions. Bioinformatic and mechanistic analysis of the palmerolide PKS-NRPS biosynthetic pathway from the microbiome of an Antarctic ascidian Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the [~]75 kb trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicohabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression. 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. SARS-Cov-2 Spike binding to ACE2 is stronger and longer ranged due to glycan interaction Highly detailed steered Molecular Dynamics simulations are performed on differently glycosylated receptor binding domains of the SARS-CoV-2 spike protein. The binding strength and the binding range increases with glycosylation. The interaction energy rises very quickly with pulling the proteins apart and only slowly drops at larger distances. We see a catch slip type behavior where interactions during pulling break and are taken over by new interactions forming. The dominant interaction mode are hydrogen bonds but Lennard-Jones and electrostatic interactions are relevant as well.

Statement of SignificanceGlycosylation of the receptor binding domain of the Spike protein of SARS-CoV-2 as well as the ACE2 receptor leads to stronger and longer ranged binding interactions between the proteins. Particularly, at shorter distances the interactions are between residues of the proteins themselves whereas at larger distances these interactions are mediated by the glycans.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=164 SRC="FIGDIR/small/452507v2_ufig1.gif" ALT="Figure 1">
View larger version (32K):
org.highwire.dtl.DTLVardef@7b9021org.highwire.dtl.DTLVardef@338aadorg.highwire.dtl.DTLVardef@1d5921eorg.highwire.dtl.DTLVardef@6c27e7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. 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. Single cell profiling of functionally cured Chronic Hepatitis B patients reveals the emergence of activated innate and an altered adaptive immune response in the intra-hepatic environment Hepatitis B surface antigen (HBsAg) loss or functional cure (FC), is considered the desirable therapeutic outcome for chronic hepatitis B (CHB) patients. However, the immuno-pathological biomarkers and underlying mechanisms remain unclear. Here we present a comprehensive single cell-transcriptomic atlas together with immune-phenotyping of disease-associated cell states (DACS) isolated from intra-hepatic tissue and matched PBMCs of either CHB or FC patients. We find that the intra-hepatic environment displays specific cell identities and molecular signatures that are distinct from PBMCs. FC is associated with emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes (CD4-CTLs), and an activated innate response represented by liver-resident natural killer (LR-NK) cells. Overall, these findings provide novel insights into immuno-pathological cell states associated with FC that could serve as prognostic biomarkers. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. '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