Cytokinin oxidase/dehydrogenase family genes exhibit functional divergence and overlap in rice growth and development Cytokinins play key roles in plant growth and development; hence, cytokinin biosynthesis and degradation have been extensively studied. Cytokinin oxidase/dehydrogenases (CKXs) are a group of enzymes that regulate oxidative cleavage to maintain cytokinin homeostasis. In rice, 11 OsCKX genes have been identified to date; however, most of their functions remain unknown. Here, we comprehensively analyzed the expression patterns and functions of OsCKX genes. Using CRISPR/Cas9 technology, we constructed mutants of all OsCKX genes to determine the functions of OsCKXs in rice development. The results revealed that the single osckx and higher-order osckx4 osckx9 mutant lines showed functional overlap and subfunctionalization. Notably, the osckx1 osckx2 and osckx4 osckx9 double mutants displayed contrasting phenotypic changes in tiller number and panicle size compared to the wild type. Moreover, we identified several genes with significantly altered expression in osckx4 and osckx9 single and double mutant plants. Many differentially expressed genes were found to be associated with auxin and cytokinin pathways. Additionally, the cytokinins in osckx4 osckx9 mutants were increased compared to the wild type. Overall, our findings provide new insights into the functions of OsCKX genes in rice growth and may be used as a foundation for future studies aimed at improving rice yield. Key words: Cytokinin, expression pattern, OsCKX, panicle, phenotype, rice, tiller

HighlightThe osckx4 osckx9 double mutant had a significantly greater number of tillers, whereas the osckx1 osckx2 double mutant showed the opposite phenotypic change, compared to the wild type
S-acylation stabilizes ligand-induced receptor kinase complex formation during plant pattern-triggered immune signalling. Plant receptor kinases are key transducers of physical extracellular stimuli, such as the presence of beneficial or pathogenic microbes or secreted signalling molecules. Receptor kinases are regulated by numerous post-translational modifications. Here, using the bacterial flagellin perceiving receptor kinase FLS2, we show that S-acylation at an evolutionarily conserved cysteine is crucial for function. S-acylation involves the addition of long-chain fatty acids to cysteine residues within proteins, altering their biophysical properties and behaviour within the membrane environment. We observe S-acylation of FLS2 at C-terminal kinase domain cysteine residues within minutes of treatment with flg22 ligand, and in a BAK1 co-receptor dependent manner. We demonstrate that S-acylation is essential for FLS2-mediated immune signalling, including anti-bacterial immunity. Similarly, mutating the corresponding conserved cysteine residue in the immune receptor kinase EFR supressed elf18 mediated signalling. Biochemical analysis of unstimulated and activated FLS2 containing complexes using detergents and native membrane DIBMA nanodiscs indicates that S-acylation assists the stabilization of activated receptor kinase complexes withing the membrane environment to increase signalling efficiency.
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. Reduced metabolic activity in established Gardnerella spp. biofilms contributes to protection from the bactericidal effects of metronidazole Gardnerella spp. can form biofilm and it has been suggested that failure of antibiotic treatment of bacterial vaginosis and recurrent vaginosis are linked to this property, however no specific mechanisms have been identified. Here, we performed broth microdilution assays to measure the minimum inhibitory concentration (MIC) of metronidazole for thirty-five Gardnerella isolates in one medium in which Gardnerella spp. grow primarily as biofilm and another in which growth is primarily planktonic. The MIC of isolates in the two conditions were highly correlated (R2= 0.69, p <0.001). No significant reduction in viable cell count of 7/9 isolates was observed when established biofilms were exposed to metronidazole at levels double the MIC. Subsequent measurements of metabolic activity in established biofilms formed by a larger collection of 31 isolates showed reduced metabolic activity following treatment with 128 g/ml of metronidazole relative to untreated controls in 27/31 cases. The amount of biofilm produced by Gardnerella isolates was not enhanced by metronidazole exposure, even at sub-MIC levels. Our results suggest that reduced metabolic activity in Gardnerella spp. growing in established biofilms may be a mechanism of protection from the bactericidal effects of metronidazole. Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. Gene expression underlying parenting and being parented reflect canalization, not plasticity, in response to environmental stress Parenting buffers offspring from hostile environments, but it is not clear how or if the genes that underlie parenting change their expression under environmental stress. We recently demonstrated that for the subsocial carrion beetle, Nicrophorus orbicollis, temperature during parenting does not affect parenting phenotypes. Here, we ask if transcriptional changes associated with parenting are likewise robust to environmental stress. The absence of a transcriptional response for parenting under stress would suggest that the genetic  programs for parenting and being parented are canalized. Conversely, a robust transcriptional response would suggest that plasticity of underlying gene expression is critical for maintaining behavioral stability, and that these mechanisms provide a potential target for selection in the face of environmental change. We test these alternatives by characterizing gene expression of parents and offspring with and without parent-offspring interactions under a benign and a stressful temperature. We found that parent-offspring interactions elicit distinct transcriptional responses of parents and larvae irrespective of temperature. We further detected robust changes of gene expression in beetles breeding at 24{degrees} C compared to 20{degrees} C irrespective of family interaction. However, no strong interaction between parent-offspring interaction and temperature was detected for either parents or larvae. We therefore conclude that canalization, not plasticity of gene expression, most likely explains the absence of behavioral plasticity under thermal stress. This result suggests that species may not have the genetic variation needed to respond to all environmental change, especially for complex phenotypes.

IMPACT STATEMENTParenting is thought to evolve to buffer offspring from harsh or inhospitable conditions. However, once evolved, does parenting continue to buffer against new environmental challenges such as might occur under global warming? One way organisms achieve different phenotypes without evolving different genes for each trait is by changing timing or extent of gene expression. Changes in gene expression influencing physiological traits such as hypoxia have been shown to facilitate adaptation to changing environments. It is less clear if changes in gene expression are also important for behavioral plasticity in complex social traits, such as parental care. In this study, we ask whether families of the subsocial burying beetle, Nicrophorus orbicollis, show altered expression of parent or offspring genes that may help compensate for reduced offspring success under global warming. To test the extent that gene expression underlies behavioral plasticity in parenting and offspring responses to parenting, we characterize the gene expression profiles of whole families (males, females, and larvae) with and without parent-offspring interactions under two thermal conditions: 20{degrees} C and 24{degrees} C. We found that heat stress predominantly induced downregulation of genes, consistent with a constrained acclimation response. However, when temperature was considered in conjunction with parent-offspring interactions, genetic plasticity was extremely limited. That is, the genes underlying parent-offspring interactions did not show strong signatures of differential expression under more stressful thermal conditions. Given that we previously found no effects of temperature on parental behavior, we suggest that constrained, not plastic gene expression, underpins stability of parent-offspring interactions and limits the buffering capacity for family interactions under environmental stress. Insights from this study deepen our understanding of the interplay between genetic and phenotypic plasticity in complex behavioral traits and how traits such as parenting may limit responses to novel environmental challenges. Beyond linear regression: mapping models in cognitive neuroscience should align with research goals Many cognitive neuroscience studies use large feature sets to predict and interpret brain activity patterns. Feature sets take many forms, from human stimulus annotations to representations in deep neural networks. Of crucial importance in all these studies is the mapping model, which defines the space of possible relationships between features and neural data. Until recently, most encoding and decoding studies have used linear mapping models. Increased availability of large datasets and computing resources has allowed some researchers to employ more flexible nonlinear mapping models instead; however, the question of whether nonlinear mapping models can yield meaningful scientific insights remains debated. Here, we discuss the choice of a mapping model in the context of three overarching desiderata: predictive accuracy, interpretability, and biological plausibility. We show that, contrary to popular intuition, these desiderata do not map cleanly onto the linear/nonlinear divide; instead, each desideratum can refer to multiple research goals, each of which imposes its own constraints on the mapping model. Moreover, we argue that, instead of categorically treating the mapping models as linear or nonlinear, we should instead aim to estimate the complexity of these models. We show that, in many cases, complexity provides a more accurate reflection of restrictions imposed by various research goals. Finally, we outline several complexity metrics that can be used to effectively evaluate mapping models. Enhancing CDK4/6 inhibitor therapy for medulloblastoma using nanoparticle delivery and scRNA-seq-guided combination with sapanisertib CDK4/6 inhibitors hold promise for brain tumor treatment, but efficacy has been limited by recurrence in both preclinical models and clinical trials. To address recurrence, we tested a nanoparticle formulation of the CDK4/6 inhibitor palbociclib (POx-palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma. We then analyzed medulloblastomas in mice receiving palbociclib treatment, and compared the efficacy of combining palbociclib with specific inhibitors suggested by our analysis. POx-Palbo showed reduced toxicity compared to conventional palbociclib, was tolerable in parenteral administration, improved CNS pharmacokinetics, and extended survival of mice with medulloblastoma. Recurrence, however, remained problematic as fractions of tumor cells proliferated during therapy. ScRNA-seq identified a gene expression pattern unique to proliferating medulloblastoma cells in POx-Palbo-treated mice, marked by up-regulation of the glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes. Reduced mTORC1 signaling, suggested by ribosomal suppression in POx-Palbo-treated tumors was confirmed by decreased 4EBP1 phosphorylation (p4EBP1). Further reducing mTORC1 activity by combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects, with increased suppression of both pRB and p4EBP1, and prolonged mouse survival compared to either agent alone. In contrast, targeting cell cycle progression by combining POx-Palbo with the SHH-pathway inhibitor vismodegib, or with the replication-targeting agents gemcitabine or etoposide, failed to enhance efficacy. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment, and identify POx-palbo plus sapanisertib as effective combinatorial therapy for SHH medulloblastoma. This combination may be appropriate for testing in patients with recurrence, who need new options. The nucleus bypasses obstacles by deforming like a drop with surface tension mediated by lamin A/C Migrating cells must deform their stiff cell nucleus to move through pores and fibers in tissue. Lamin A/C is known to hinder cell migration by limiting nuclear deformation and passage through confining channels, but its role in nuclear deformation and passage through fibrous environments is less clear. We studied cell and nuclear migration through discrete, closely spaced, slender obstacles which mimic the mechanical properties of collagen fibers. Nuclei bypassed slender obstacles while preserving their overall morphology by deforming around them with deep local invaginations of little resisting force. The obstacles did not impede the nuclear trajectory or cause a rupture of the nuclear envelope. Nuclei likewise deformed around single collagen fibers in cells migrating in 3D collagen gels. In contrast to its limiting role in nuclear passage through confining channels, lamin A/C facilitated nuclear deformation and passage through fibrous environments; nuclei in lamin-null (Lmna-/-) cells lost their overall morphology and became entangled on the obstacles. Analogous to surface tension-mediated deformation of a liquid drop, lamin A/C imparts a surface tension on the nucleus that allows nuclear invaginations with little mechanical resistance, preventing nuclear entanglement and allowing nuclear passage through fibrous environments. Implicating Gene and Cell Networks Responsible for Differential COVID-19 Host Responses via an Interactive Single Cell Web Portal Numerous studies have provided single-cell transcriptome profiles of host responses to SARS-CoV-2 infection. Critically lacking however is a datamine that allows users to compare and explore cell profiles to gain insights and develop new hypotheses. To accomplish this, we harmonized datasets from COVID-19 and other control condition blood, bronchoalveolar lavage, and tissue samples, and derived a compendium of gene signature modules per cell type, subtype, clinical condition, and compartment. We demonstrate approaches to probe these via a new interactive web portal (http://toppcell.cchmc.org/ COVID-19). As examples, we develop three hypotheses: (1) a multicellular signaling cascade among alternatively differentiated monocyte-derived macrophages whose tasks include T cell recruitment and activation; (2) novel platelet subtypes with drastically modulated expression of genes responsible for adhesion, coagulation and thrombosis; and (3) a multilineage cell activator network able to drive extrafollicular B maturation via an ensemble of genes strongly associated with risk for developing post-viral autoimmunity. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. N-VEGF, the autoregulatory arm of VEGF-A Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify changes in metabolic flux, and to quantify contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in liver.

HighlightsO_LIWe created OMELET to infer metabolic flux and its regulation from multi-omic data.
C_LIO_LIGluconeogenic and pyruvate cycle fluxes increased in fasting ob/ob mice.
C_LIO_LITranscripts increases mediated the increase in gluconeogenic fluxes in ob/ob mice.
C_LIO_LIIncreases in transcripts and substrates enhanced pyruvate cycle flux in ob/ob mice.
C_LI Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Glucose dynamics during ozone exposure measured using radiotelemetry: Stress drivers BackgroundStress-related neurobehavioral and metabolic disorders are associated with altered circulating adrenal-derived hormones and hyperglycemia. Temporal assessment of glucose and these hormones is critical for insights on an individuals health.

ObjectivesHere we used implantable-telemetry in rats to assess real-time changes in circulating glucose during and after exposure to the air pollutant ozone, and link responses to circulating neuroendocrine stress and metabolic hormones. We also compared rodent glucose and corticosterone (cortisol in humans) responses to humans exposed to ozone.

MethodsFirst, using a cross-over design, we monitored glucose levels during single or repeated ozone exposures (0.0, 0.2, 0.4 and 0.8-ppm) and non-exposure periods in male Wistar-Kyoto-rats implanted with glucose-telemeters. A second cohort of un-implanted rats was exposed to ozone (0.0, 0.4 or 0.8-ppm) for 30-min, 1-hour, 2-hour, or 4-hour with hormones measured immediately after exposure. Then we assessed glucose metabolism in sham and adrenalectomized rats with or without pharmacological interventions of adrenergic and glucocorticoid receptors. Finally, we assessed glucose and cortisol in serum samples from a clinical study involving exposure of human volunteers to air or 0.3 ppm ozone.

ResultsOzone (0.8-ppm) caused hyperglycemia and hypothermia beginning 90-min into exposure, with reversal of effects 4-6 hours post-exposure. Glucose monitoring during four daily 4-hour ozone exposures revealed duration of hyperglycemia, adaptation, and diurnal variations. Ozone-induced hyperglycemia was preceded by increased adrenocorticotropic hormone, corticosterone, and epinephrine, but depletion of thyroid-stimulating, prolactin, and luteinizing hormones. Hyperglycemia was inhibited in rats that were adrenalectomized and/or treated with glucocorticoid inhibitor. The depletion of cortisol was dampened in humans exposed to ozone during intermittent exercise.

DiscussionWe demonstrate for the first time the temporality of neuroendocrine-stress-mediated biological sequalae responsible for ozone-induced metabolic dysfunction as exposure occurs. Real-time glucose monitoring with stress hormones assessment may be useful in identifying interactions among pollutants and stress-related illnesses. Towards a cumulative science of vocal markers of autism: a cross-linguistic meta-analysis-based investigation of acoustic markers in American and Danish autistic children Acoustic atypicalities in speech production are argued to be potential markers of clinical features in Autism Spectrum Disorder (ASD). A recent meta-analysis highlighted shortcomings in the field, in particular small sample sizes and study heterogeneity (Fusaroli et al., 2017). We showcase a cumulative (i.e., explicitly building on previous studies both conceptually and statistically) yet self-correcting (i.e., critically assessing the impact of cumulative statistical techniques) approach to prosody in ASD to overcome these issues.

We relied on the recommendations contained in the meta-analysis to build and analyze a cross-linguistic corpus of multiple speech productions in 77 autistic and 72 neurotypical children and adolescents (>1000 recordings in Danish and US English). We used meta-analytically informed and skeptical priors, with informed priors leading to more generalizable inference. We replicated findings of a minimal cross-linguistically reliable distinctive acoustic profile for ASD (higher pitch and longer pauses) with moderate effect sizes. We identified novel reliable differences between the two groups for normalized amplitude quotient, maxima dispersion quotient, and creakiness. However, the differences were small, and there is likely no one acoustic profile characterizing all autistic individuals. We identified reliable relations of acoustic features with individual differences (age, gender), and clinical features (speech rate and ADOS sub-scores).

Besides cumulatively building our understanding of acoustic atypicalities in ASD, the study shows how to use systematic reviews and meta-analyses to guide the design and analysis of follow-up studies. We indicate future directions: larger and more diverse cross-linguistic datasets, focus on heterogeneity, self-critical cumulative approaches and open science.

Lay SummaryAutistic individuals are reported to speak in distinctive ways. Distinctive vocal production can affect social interactions and social development and could represent a noninvasive way to support the assessment of ASD. We systematically checked whether acoustic atypicalities highlighted in previous articles could be actually found across multiple recordings and two languages. We find a minimal acoustic profile of ASD: higher pitch, longer pauses, increased hoarseness and creakiness of the voice. However, there is much individual variability (by age, sex, language, and clinical characteristics). This suggests that the search for one common "autistic voice" might be naive and more fine-grained approaches are needed. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Embedded 3D printing in self-healing annealable composites for precise patterning of functionally mature human neural constructs Human in vitro models of neural tissue with controllable cellular identity, tunable microenvironment, and defined spatial arrangement are needed to facilitate studies of brain development and disease. Towards this end, embedded printing in jammed microgel supports (i.e., granular gels) holds great promise as it allows precise and programmable patterning of extremely soft and compliant tissue constructs. However, in contrast to the vast material landscape available for bulk hydrogels, granular printing support formulations are restricted to a handful of materials without the ability for facile adjustment of biofunctional properties of the cellular microenvironment. Therefore, there has been a need for novel materials that take advantage of versatile biomimicry of bulk hydrogels while providing high-fidelity support for embedded printing akin to granular gels. To address this need, we present a modular platform for bioengineering of neuronal networks via direct embedded 3D printing of human stem cells inside Self-Healing Annealable Particle-Extracellular matrix (SHAPE) composites. SHAPE composites consist of soft microgels immersed in viscous extracellular-matrix solution to enable precise freeform patterning of human stem cells and consequent generation and long-term maintenance of mature subtype-specific neurons that extend projections within the volume of the annealed support. The developed approach further allows multi-ink deposition, live spatial and temporal monitoring of oxygen levels, as well as creation of vascular channels. Due to its modularity, SHAPE biomanufacturing toolbox not only offers a solution for functional modeling of mechanically sensitive neural constructs, but also has potential to be applied to a wide range of biomaterials with different crosslinking mechanisms to model tissues and diseases where recapitulation of complex architectural features and topological cues is essential. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. 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. Colicin E1 opens its hinge to plug TolC The double membrane architecture of Gram-negative bacteria forms a barrier that is effectively impermeable to extracellular threats. Bacteriocin proteins evolved to exploit the accessible, surface-exposed proteins embedded in the outer membrane to deliver cytotoxic cargo. Colicin E1 is a bacteriocin produced by, and lethal to, Escherichia coli that hijacks the outer membrane proteins TolC and BtuB to enter the cell. Here we capture the colicin E1 translocation domain inside its membrane receptor, TolC, by high-resolution cryoEM, the first reported structure of a bacteriocin bound to TolC. Colicin E1 binds stably to TolC as an open hinge through the TolC pore--an architectural rearrangement from colicin E1s unbound conformation. This binding is stable in live E. coli cells as indicated by single-molecule fluorescence microscopy. Finally, colicin E1 fragments binding to TolC plugs the channel, inhibiting its native efflux function as an antibiotic efflux pump and heightening susceptibility to three antibiotic classes. In addition to demonstrating that these protein fragments are useful starting points for developing novel antibiotic potentiators, this method could be expanded to other colicins to inhibit other outer membrane protein functions. Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Loss of Resf1 reduces the efficiency of embryonic stem cell self-renewal and germline entry Retroelement silencing factor 1 (RESF1) interacts with the key regulators of mouse embryonic stem cells (ESCs) OCT4 and NANOG, and its absence results in sterility of mice. However, the function of RESF1 in ESCs and germ line specification is poorly understood. In this study, we used Resf1 knockout cell lines to determine the requirements of RESF1 for ESCs self-renewal and for in vitro specification of ESCs into primordial germ cell-like cells (PGCLCs). We found that deletion of Resf1 in ESCs cultured in serum and LIF reduces self-renewal potential whereas episomal expression of RESF1 has a modest positive effect on ESC self-renewal. In addition, RESF1 is not required for the capacity of NANOG and its downstream target ESRRB to drive self-renewal in the absence of LIF. However, Resf1 deletion reduces efficiency of PGCLC differentiation in vitro. These results identify Resf1 as a novel player in the regulation of pluripotent stem cells and germ cell specification. Estimating the fundamental niche: accounting for the uneven availability of existing environmental conditions in the calibration area Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche.