Revisiting regulatory decoherence and phenotypic integration: accounting for temporal bias in co-expression analyses Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets.
An estimate of the deepest branches of the tree of life from ancient vertically-evolving genes Core gene phylogenies provide a window into early evolution, but different gene sets and analytical methods have yielded substantially different views of the tree of life. Trees inferred from a small set of universal core genes have typically supported a long branch separating the archaeal and bacterial domains. By contrast, recent analyses of a broader set of non-ribosomal genes have suggested that Archaea may be less divergent from Bacteria, and that estimates of inter-domain distance are inflated due to accelerated evolution of ribosomal proteins along the inter-domain branch. Resolving this debate is key to determining the diversity of the archaeal and bacterial domains, the shape of the tree of life, and our understanding of the early course of cellular evolution. Here, we investigate the evolutionary history of the marker genes key to the debate. We show that estimates of a reduced Archaea-Bacteria (AB) branch length result from inter-domain gene transfers and hidden paralogy in the expanded marker gene set. By contrast, analysis of a broad range of manually curated marker gene datasets from an evenly sampled set of 700 Archaea and Bacteria reveal that current methods likely underestimate the AB branch length due to substitutional saturation and poor model fit; that the best-performing phylogenetic markers tend to support longer inter-domain branch lengths; and that the AB branch lengths of ribosomal and non-ribosomal marker genes are statistically indistinguishable. Furthermore, our phylogeny inferred from the 27 highest-ranked marker genes recovers a clade of DPANN at the base of the Archaea, and places CPR within Bacteria as the sister group to the Chloroflexota.
Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Loss of Stk11 in basolateral amygdalar projection neurons impairs taste aversion learning by altering the temporal pattern of taste response plasticity in gustatory cortex Gustatory cortex (GC) responds to tastes on the tongue with dynamic ensemble activity that represents first the presence, then the identity, and finally the hedonic value (palatability) of tastes. This final state of the taste response is uniquely altered by conditioned taste aversion (CTA) - a powerful one-trial learning paradigm in which a taste becomes aversive after association with gastric malaise - a process requiring coordination between GC and basolateral amygdala (BLA). Previously, we found that one key requirement for learning in this circuit is expression of the serine/threonine kinase 11 (Stk11) gene (a tumor suppression gene which has only recently been associated with learning; Levitan et al., 2020). When Stk11 is knocked out in BLA projection neurons (BLApn), CTA learning fails to occur. Here we have examined how learning-related response plasticity in GC taste responses is impacted by the knockout of Stk11 in BLApn. Contrary to the commonly held assumption that a lack of learning means a lack of such plasticity, but consistent with the fact that Stk11KO has been shown to increase the excitability of BLApn, our data reveal that the knockout of Stk11 in BLApn does not eliminate plasticity; rather, it shifts the impact of CTA training on GC taste responses to an earlier, learning-inappropriate epoch. Even naive taste representations are altered--specifically, the pattern of similarities and differences among the different taste responses are rendered abnormal by Stk11 KO, and these relationships fail to change with training. Finally, the latency of behavior-related dynamic ensemble features of the GC taste response, which is also abnormal in naive KO mice, is rendered disorganized by CTA. Together, these results suggest that Stk11 plays a role in governing the coordination of GC activity by BLA, and demonstrate that alterations in the function of BLApn caused by Stk11KO inhibit learning not by inhibiting plasticity but by changing its temporal properties. Isotopic signatures induced by upwelling reveal regional fish populations in Lake Tanganyika Lake Tanganyikas pelagic fish sustain the second largest inland fishery in Africa and are under pressure from heavy fishing and global warming related increases in stratification. Only little is known about whether basin-scale hydrodynamics - including a more stratified north and an upwelling-driven south - induce ecological and genetic differences among populations of highly mobile, pelagic fish inhabiting these different areas. Here, we examine whether the basin-scale dynamics leave distinct isotopic imprints in the pelagic fish of Lake Tanganyika, which may reveal differences in habitat, diet, or lipid content. We conducted two lake-wide campaigns during different seasons and collected physical, nutrient, chlorophyll, phytoplankton and zooplankton data. Additionally, we analyzed the pelagic fish - the clupeids Stolothrissa tanganicae, Limnothrissa miodon and four Lates species - for their isotopic and elemental carbon (C) and nitrogen (N) compositions. The {delta}13C values were significantly higher in the productive south after the upwelling/mixing period across all trophic levels, implying that the fish have regional foraging grounds, and thus record these latitudinal isotope gradients. By combining our isotope data with genetics, we demonstrate that the fish form regional populations on a seasonal to multiannual time scale. Based on {delta}15N and C:N ratios, we found no strong evidence for varying diets or lipid contents between those regional populations.

Additional analyses revealed that isotopic variations between specimens from the same location are not linked to genetic differences. We suggest that the development of basinscale ecological differences in response to the prevailing hydrodynamic regimes may be inhibited by lake-wide gene flow on the long term. Our findings show that the pelagic fish species are genetically adapted to the whole lake, but they form regional populations on short time scales. This implies that sustainable management strategies may adopt basin-scale fishing quotas. Adaptive Wild Introgression Mediate Temporal Wheat Root Water Influxes under Drought Stress through Plant-Atmospheric Interrelations Drought intensity as experienced by plants depends upon soil moisture status and atmospheric variables such as temperature, radiation, and air vapour pressure deficit (VPD). Although the role of shoot architecture with these edaphic and atmospheric factors is well-characterized, the extent to which shoot and root dynamic interactions as a continuum are controlled by genotypic variation is less known. Here, we targeted these interactions using a wild emmer introgression line (IL20) with a distinct drought-induced shift in the shoot-to-root ratio and its drought-sensitive recurrent parent Svevo. Using a gravimetric platform, we show that IL20 maintained higher root water influx and gas exchange under drought stress, which supported a greater growth. Interestingly, the advantage of IL20 in root water influx and transpiration was expressed earlier during the daily diurnal cycle under lower VPD and therefore supported higher transpiration efficiency. Application of structural equation model indicates that under drought, VPD and radiation are antagonistic to transpiration rate, whereas the root water influx operates as feedback for the higher atmospheric responsiveness of leaves. Collectively, our results suggest that a drought-induced shift in root-to-shoot ratio can improve plant water uptake potential in a short preferable time window determined by both water and atmospheric parameters. Leaf Aqueous Extract of Manilkara hexandra influenced Glucose Metabolism in fish: An observation with Labeo rohita Fingerlings as model organism Manilkara hexandra (Roxb.) Dubard (common name Khirni, family: Sapotaceae) is an evergreen traditionally used medicinal plant. The present study aimed to see the effects of leaf aqueous extract of M. hexandra (LAEMH) on digestive and glucose metabolic enzyme action and serum glucose level in rohu, Labeo rohita, fingerlings. Experimental fish were fed a basal diet (Group I), and diets supplemented with LAEMH at 300 mg (Group II) and 600 mg (Group III) kg-1body weight. A significant reduction in serum glucose was noticed in the treated groups when compared with the control group. The reduced amylase activity was noticed in the treated fingerlings in a dose-dependent manner. However, lipase and protease activities didnt differ significantly among the experimental groups. Reduced serum glucose level might be correlated with the decline in the activity of digestive amylase in fish. Group II (3.35{+/-}0.19 U) and Group III (3.49{+/-}0.13 U) were recorded with reduced glucose-6-phosphatase activities than the control group (4.4{+/-}0.39 U). Moreover, the study revealed a decline in fructose-1,6-phosphatase activities in the treated groups in comparison to the control group. The decline in the activities of the metabolic enzymes might be associated with the non-availability of glucose owing to reduced activity of the digestive amylase in the treated groups. In conclusion, the present study established the hypoglycaemic effect of the leaf aqueous extract of M. hexandra in a fish model. Effect of SARS-CoV-2 spike mutations on animal ACE2 usage and in vitro neutralization sensitivity The emergence of SARS-CoV-2 variants poses greater challenges to the control of COVID-19 pandemic. Here, we parallelly investigated three important characteristics of seven SARS-CoV-2 variants, including two mink-associated variants, the B.1.617.1 variant, and the four WHO-designated variants of concerns (B.1.1.7, B.1.351, P.1, and B.1.617.2). We first investigated the ability of these variants to bind and use animal ACE2 orthologs as entry receptor. We found that, in contrast to a prototype variant, the B.1.1.7, B.1.351, and P.1 variants had significantly enhanced affinities to cattle, pig, and mouse ACE2 proteins, suggesting increased susceptibility of these species to these SARS-CoV-2 variants. We then evaluated in vitro neutralization sensitivity of these variants to four monoclonal antibodies in clinical use. We observed that all the variants were partially or completely resistant against at least one of the four tested antibodies, with B.1.351 and P.1 showing significant resistance to three of them. As ACE2-Ig is a broad-spectrum anti-SARS-CoV-2 drug candidate, we then evaluated in vitro neutralization sensitivity of these variants to eight ACE2-Ig constructs previously described in three different studies. All the SARS-CoV-2 variants were efficiently neutralized by these ACE2-Ig constructs. Interestingly, compared to the prototype variant, most tested variants including the variants of concern B.1.1.7, B.1.351, P.1, and B.1.617.2 showed significantly increased (up to [~]15-fold) neutralization sensitivity to ACE2-Ig constructs that are not heavily mutated in the spike-binding interface of the soluble ACE2 domain, suggesting that SARS-CoV-2 evolves toward better utilizing ACE2, and that ACE2-Ig is an attractive drug candidate for coping with SARS-CoV-2 mutations. Persistent Immune and Clotting Dysfunction Detected in Saliva and Blood Plasma after COVID-19 A growing number of studies indicate that coronavirus disease 2019 (COVID-19) is associated with inflammatory sequelae, but molecular signatures governing the normal vs. pathologic convalescence process have not been well-delineated. We characterized global immune and proteome responses in matched plasma and saliva samples obtained from COVID-19 patients collected between 4-6 weeks after initial clinical symptoms resolved. Convalescent subjects showed robust IgA and IgG responses and positive antibody correlations between matched saliva and plasma samples. However, global shotgun proteomics revealed persistent inflammatory patterns in convalescent samples including dysfunction of salivary innate immune cells and clotting factors in plasma (e.g., fibrinogen and antithrombin), with positive correlations to acute COVID-19 disease severity. Saliva samples were characterized by higher concentrations of IgA, and proteomics showed altered pathways that correlated positively with IgA levels. Our study positions saliva as a viable fluid to monitor immunity beyond plasma to document COVID-19 immune, inflammatory, and coagulation-related sequelae. Learning and motor inhibitory control in crows and domestic chickens Cognitive abilities allow animals to navigate through complex, fluctuating environments. In the present study, we tested the performance of a captive group of eight crows, Corvus corone, and ten domestic chickens, Gallus gallus domesticus, in the cylinder task, as a test of motor inhibitory control, and reversal learning as a measure of learning ability and behavioural flexibility. Four crows and nine chickens completed the cylinder task, eight crows and six chickens completed the reversal learning experiment. Crows performed better in the cylinder task compared to chickens. In the reversal learning experiment, species did not significantly differ in the number of trials until learning criterion was reached. Performance in the reversal learning experiment did not correlate with performance in the cylinder task in chickens. Our results suggest crows to possess better motor inhibitory control compared to chickens. In contrast, learning performance in a reversal learning task did not differ between the species, indicating similar levels of behavioural flexibility. Interestingly, we describe notable individual differences in performance. We stress the importance not only to compare cognitive performance between species but also between individuals of the same species when investigating the evolution of cognitive skills. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. Interpretable deep learning model reveals subsequences of various functions for long non-coding RNA identification Long non-coding RNAs (lncRNAs) play crucial roles in many biological processes and are implicated in several diseases. With the next-generation sequencing technologies, substantial un-annotated transcripts have been discovered. Classifying unannotated transcripts using biological experiments is more time-consuming and expensive than computational approaches. Several tools for identifying long non-coding RNAs are available. These tools, however, did not explain which features in their tools contributed to the prediction results. Here, we present Xlnc1DCNN, a tool for distinguishing long non-coding RNAs (lncRNAs) from protein-coding transcripts (PCTs) using a one-dimensional convolutional neural network with prediction explanations. The evaluation results of the human test set showed that Xlnc1DCNN outperformed other state-of-the-art tools in terms of accuracy and F1-score. The explanation results revealed that lncRNA transcripts were mainly identified as sequences with no conserved regions or with a region of transmembrane helix while protein-coding transcripts were mostly classified by conserved protein domains or families. The explanation results also conveyed the probably inconsistent annotations among the public databases, lncRNA transcripts which contain protein domains or families, as well as protein-coding transcripts which are nonsense-mediated decay or processed transcripts. Xlnc1DCNN is freely available at https://github.com/cucpbioinfo/Xlnc1DCNN. Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. MTG16 (CBFA2T3) regulates colonic epithelial differentiation, colitis, and tumorigenesis by repressing E protein transcription factors Aberrant epithelial differentiation and regeneration contribute to colon pathologies including inflammatory bowel disease (IBD) and colitis-associated cancer (CAC). MTG16 (CBFA2T3) is a transcriptional corepressor expressed in the colonic epithelium. MTG16 deficiency in mice exacerbates colitis and increases tumor burden in CAC, though the underlying mechanisms remain unclear. Here, we identified MTG16 as a central mediator of epithelial differentiation, promoting goblet and restraining enteroendocrine cell development in homeostasis and enabling regeneration following dextran sulfate sodium (DSS)-induced colitis. Transcriptomic analyses implicated increased E box-binding transcription factor (E protein) activity in MTG16-deficient colon crypts. Using a novel mouse model with a point mutation that disrupts MTG16:E protein complex formation (Mtg16P209T), we established that MTG16 exerts control over colonic epithelial differentiation and regeneration by repressing E protein-mediated transcription. Mimicking murine colitis, MTG16 expression was increased in biopsies from patients with active IBD compared to unaffected controls. Finally, uncoupling MTG16:E protein interactions only partially phenocopied the enhanced tumorigenicity of Mtg16-/- colon in the azoxymethane(AOM)/DSS-induced model of CAC, indicating that MTG16 protects from tumorigenesis through additional mechanisms. Collectively, our results demonstrate that MTG16, via its repression of E protein targets, is a key regulator of cell fate decisions during colon homeostasis, colitis, and cancer.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/467178v3_ufig1.gif" ALT="Figure 1">
View larger version (30K):
org.highwire.dtl.DTLVardef@21948borg.highwire.dtl.DTLVardef@9f60eeorg.highwire.dtl.DTLVardef@194b9fdorg.highwire.dtl.DTLVardef@13eb1a7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Intragenomic rearrangements of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses. Specifically, we report numerous genomic insertions of 5-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain 5-UTR-derived inserted sequences. In limited instances, there is mounting evidence that these insertions alter the fundamental biological properties of mutant viruses. Intragenomic rearrangements add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

SignificanceUnderstanding mechanisms of variation in coronaviruses is vital to control of their associated diseases. Beyond point mutations, insertions, deletions and recombination, we here describe for the first time intragenomic rearrangements and their relevance to changes in transmissibility, immune escape and/or virulence documented during the SARS-CoV-2 pandemic. Full-length Genome of a Ogataea polymorpha strain CBS4732 ura3Δ reveals large duplicated segments in subtelomeric regions BackgroundCurrently, methylotrophic yeasts (e.g., Pichia pastoris, Ogataea polymorpha, and Candida boindii) are subjects of intense genomics studies in basic research and industrial applications. In the genus Ogataea, most research is focused on three basic O. polymorpha strains--CBS4732, NCYC495, and DL-1. However, the relationship between CBS4732, NCYC495, and DL-1 remains unclear, as the genomic differences between them have not be exactly determined without their high-quality complete genomes. As a nutritionally deficient mutant derived from CBS4732, the O. polymorpha strain CBS4732 ura3{Delta} (named HU-11) is being used for high-yield production of several important proteins or peptides. HU-11 has the same reference genome as CBS4732 (noted as HU-11/CBS4732), because the only genomic difference between them is a 5-bp insertion.

ResultsIn the present study, we have assembled the full-length genome of O. polymorpha HU-11/CBS4732 using high-depth PacBio and Illumina data. Long terminal repeat (LTR) retrotransposons, rDNA, 5 and 3 telomeric, subtelomeric, low complexity and other repeat regions were curated to improve the genome quality. Particularly, we detected large duplicated segments (LDSs) in the subtelomeric regions and exactly determined all the structural variations (SVs) between CBS4732 and NCYC495.

New findings mainly include(1) the genomic differences between HU-11/CBS4732 and NCYC495 include single nucleotide polymorphisms, small insertions and deletions, and only three SVs; (2) six genes were incorporated into CBS4732 from Cyberlindnera jadinii by horizontal gene transfer and may bring HU-11/CBS4732 new biological functions or physiological properties; (3) many recombination events may have occurred on chromosome 4 and 5 of CBS4732 and NCYC495 ancestors and two large segments were acquired by CBS4732 and NCYC495 from chromosome 6 and C. jadinii during recombination, respectively; and (4) the genome expansion in methylotrophic yeasts is mainly driven by large segment duplication in subtelomeric regions.

ConclusionsThe present study preliminarily revealed the complex relationship between CBS4732, NCYC495, and DL-1. The new findings provide new opportunities for in-depth understanding of genome evolution in methylotrophic yeasts and lay the foundations for the industrial applications of O. polymorpha CBS4732, NCYC495, DL-1, and their derivative strains. The full-length genome of the O. polymorpha strain HU-11/CBS4732 should be included into the NCBI RefSeq database for future studies of Ogataea spp.. Ion Permeation, Selectivity, and Electronic Polarization in Fluoride Channels Fluoride channels (Fluc) export toxic F- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein we have applied constant-pH molecular dynamics and free energy sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F-, by traversing through a non-bonded network. The extraordinary F- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared to the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F- - Fluc interactions.

O_TEXTBOXSIGNIFICANCE

A comprehensive atomistic-level computational study is presented of the mechanism of fluoride permeation through fluoride channels. The mechanism of fluoride permeation of the F- anion is established and the microscopic determinants of F- selectivity revealed. The essential nature of electronic polarization during F- permeation is also demonstrated through the computational modeling.

C_TEXTBOX HIF-dependent expression of creatine kinase brain isoform (CKB) promotes breast cancer metastasis, whereas cyclocreatine therapy impairs invasion and improves the efficacy of conventional chemotherapies The oxygen-responsive Hypoxia Inducible Factor (HIF)-1 promotes several steps of the metastatic cascade. A hypoxic gene signature is enriched in triple negative breast cancers (TNBCs) and correlates with poor patient survival. Since inhibiting the HIF transcription factors with small molecules is challenging, we sought to identify genes downstream of HIF-1 that could be targeted to block invasion and metastasis. Creatine kinase brain isoform (CKB) was identified as a highly differentially expressed gene in a screen of HIF-1 wild type and knockout mammary tumor cells derived from a transgenic model of metastatic breast cancer. CKB is a cytosolic enzyme that reversibly catalyzes the phosphorylation of creatine, generating phosphocreatine (PCr) in the forward reaction, and regenerating ATP in the reverse reaction. Creatine kinase activity is inhibited by the creatine analog cyclocreatine (cCr). Loss and gain of function genetic approaches were used in combination with cCr therapy to define the contribution of CKB expression or creatine kinase activity to cell proliferation, migration, invasion, and metastasis in ER-negative breast cancers. CKB was necessary for cell invasion in vitro and strongly promoted tumor growth and metastasis in vivo. Similarly, cyclocreatine therapy repressed cell migration, cell invasion, formation of invadopodia, and lung metastasis. Moreover, in common TNBC cell line models, the addition of cCr to conventional cytotoxic chemotherapy agents was either additive or synergistic to repress tumor cell growth. ConVarT: a search engine for orthologous variants and functional inference of human genetic variants The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely "matching variants" (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and C. elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease. Define protein variant functions with high-complexity mutagenesis libraries and enhanced mutation detection software ASMv1.0 Open reading frame (ORF) variant libraries have advanced our ability to query the functions of a large number of variants of a protein simultaneously in a single experiment. Variant libraries targeting full-length ORFs typically consists of all possible single-amino-acid substitutions and a stop codon at each amino-acid position. Because a variant differs from the template ORF by a single codon variation, variant quantification presents the most profound challenge to this technology. Efforts such as dividing a library into sub-libraries for direct sequencing, or tag-directed subassembly are practical only for small ORFs. Our approach, however, features an enhanced variant-calling tool that processes massive shotgun sequencing data, allowing one to screen a single variant library for ORFs sized up to [~]3600 bases. The variant-calling tool detects variants reliably, and also presents variant annotations in datafiles enabling analyses that have reshaped our strategies governing library design, screen deconvolution, sequencing and sequencing data analysis. Snrpb, the gene mutated in CCMS, is haploinsufficient and required in neural crest cells for proper splicing of developmentally important genes essential for craniofacial morphogenesis Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for Cerebrocostomandibular Syndrome (CCMS). We show that Snrpb heterozygous embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon-skipping and intron-retention in association with increased 5 splice site strength. We found increased exon-skipping in negative regulators of the P53-pathway, increased levels of nuclear P53, P53-target genes. However, removing TrP53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon-skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53-activity and craniofacial-specific genes both contribute to craniofacial malformations. 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. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. 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. 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. Muscle activation strategies of the vastus lateralis according to sex AimDespite men exhibiting greater muscle strength and fatigibility than women, it remains unclear if there are sex-based differences in muscle recruitment strategies e.g. motor unit (MU) recruitment and modulation of firing rate (FR) at normalised forces and during progressive increases in force.

MethodsTwenty-nine healthy male and thirty-one healthy female participants (18-35 years) were studied. Intramuscular electromyography was used to record individual motor unit potentials (MUPs) and near fibre MUPs from the vastus lateralis (VL) during 10% and 25% maximum isometric voluntary contractions (MVC), and spike-triggered averaging was used to obtain motor unit number estimates (MUNE) of the VL. Multilevel mixed-effects linear regression models were used to investigate the effects of sex at each contraction level.

ResultsMen exhibited greater muscle strength (p<0.001) and size (p<0.001) than women, with no difference in force steadiness at 10% or 25% MVC. Women had smaller MUs and higher FR at 10% MVC (both p<0.02), similar to that at 25% MVC in MU size (p=0.062) and FR (p=0.031). However, both sexes showed similar increases in MU size and FR when moving from low-to mid-level contractions. There were no sex differences in any near fibre MUP parameters or in MUNE.

ConclusionIn the vastus lateralis, women produce muscle force via different neuromuscular recruitment strategies to men which is characterised by smaller MUs discharging at higher rates. However, similar strategies are employed to increase force production from low to moderate contractions. These findings of similar proportional increases between sexes support the use of mixed sex cohorts in studies of this nature.

Key pointsO_LIIncreases in muscle force production are mediated by motor unit (MU) recruitment, and MU firing rate (FR).
C_LIO_LIWomen are underrepresented in studies of human neuromuscular research and markedly differ to men in a number of aspects of neuromuscular function, yet little is known of the recruitment strategies of each.
C_LIO_LIHere we demonstrate men and women have similar vastus lateralis MU number estimates, yet women recruit smaller MUs with higher FR than men at normalised contraction levels. However, increases in force are achieved via similar trajectories of MU recruitment and MU FR in men and women.
C_LIO_LIAlthough men and women exhibit divergent neuromuscular recruitment strategies to achieve normalised forces, increases in force are achived similarly and support the inclusion of mixed sex cohorts in studies of this nature.
C_LI