Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent one-year global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.
Identifying critical transcriptional targets of the MYC oncogene using a novel competitive precision genome editing (CGE) assay MYC is an oncogenic transcription factor that controls major pathways promoting cell growth and proliferation. MYC has been implicated in the regulation of large number of genes, but the exact target genes responsible for its proliferative effects are still not known. Here, we use a novel competitive genome editing (CGE) assay for studying the functional consequence of precise mutations of MYC binding sites on cell proliferation. The CGE method is based on precision genome editing, where a CRISPR/Cas9-induced DNA break is repaired using a template that either reconstitutes the original feature or introduces an altered sequence. Both types of repair templates harbor sequence tags that allow direct comparison between cells that carry original and mutant features and generate a large number of replicate cultures. The CGE method overcomes the limitations of CRISPR/Cas9-technology in analyzing the effect of genotype on phenotype, namely the difficulty of cutting DNA exactly at the intended site, and the decreased cell proliferation caused by the DNA cuts themselves. Importantly, it provides a powerful method for studying subtle effects elicited by mutation of individual transcription factor binding sites. We show here that E-box mutations at several MYC target gene promoters resulted in reduced cellular fitness, demonstrating a direct correlation between MYC-regulated cellular processes and MYC binding and identifying important transcriptional targets responsible for its functions.
Proinflammatory innate cytokines and metabolomic signatures predict T cell response in active COVID-19 The underlying factors contributing to the evolution of SARS-CoV-2-specific T cell responses during COVID-19 infection remain unidentified. To address this, we characterized innate and adaptive immune responses with metabolomic profiling longitudinally at three different time points (0-3, 7-9, and 14-16 days post-COVID-19 positivity) from young mildly symptomatic active COVID-19 patients infected with the ancestral virus in mid-2020. We observed that anti-RBD IgG and viral neutralization are significantly reduced against the Delta variant compared to the ancestral strain. In contrast, compared to the ancestral strain, T cell responses remain preserved against the delta and omicron variants. We determined innate immune responses during the early stage of active infection in response to TLR 3/7/8 mediated activation in PBMCs and serum metabolomic profiling. Correlation analysis indicated PBMCs-derived proinflammatory cytokines, IL-6, IL-1{beta}, and IL-23, and the abundance of plasma metabolites, pyroglutamate, itaconate, and azelaic acid, were predictive of a robust SARS-CoV-2-specific Th1 response at a later stage (two weeks after PCR positivity). These observations may contribute to designing effective vaccines and adjuvants that promote innate immune responses and metabolites to induce long-lasting anti-SARS-CoV-2 specific T cells response.

Key MessageO_LIEarly proinflammatory innate cytokine response is critical for a robust T cell response during COVID-19.
C_LIO_LIWe report a set of plasma metabolites that significantly correlate with the innate and adaptive immune responses.
C_LI

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=140 SRC="FIGDIR/small/483930v2_ufig1.gif" ALT="Figure 1">
View larger version (27K):
org.highwire.dtl.DTLVardef@1bdb656org.highwire.dtl.DTLVardef@dd23aorg.highwire.dtl.DTLVardef@164e404org.highwire.dtl.DTLVardef@b5f389_HPS_FORMAT_FIGEXP  M_FIG C_FIG Minimal requirements for a neuron to co-regulate many properties and the implications for ion channel correlations and robustness Neurons regulate their excitability by adjusting their ion channel levels. Degeneracy - achieving equivalent outcomes (excitability) using different solutions (channel combinations) - facilitates this regulation by enabling a disruptive change in one channel to be offset by compensatory changes in other channels. But neurons must co-regulate many properties. Pleiotropy - the impact of one channel on more than one property - complicates regulation because a compensatory ion channel change that restores one property to its target value often disrupts other properties. How then does a neuron simultaneously regulate multiple properties? Here we demonstrate that of the many channel combinations producing the target value for one property (the single-output solution set), few combinations produce the target value for other properties. Combinations producing the target value for two or more properties (the multi-output solution set) correspond to the intersection between single-output solution sets. Properties can be effectively co-regulated only if the number of adjustable channels (nin) exceeds the number of regulated properties (nout). Ion channel correlations emerge during homeostatic regulation when the dimensionality of solution space (nin - nout) is low. Even if each property can be regulated to its target value when considered in isolation, regulation as a whole fails if single-output solution sets do not intersect. Our results also highlight that ion channels must be co-adjusted with different ratios to regulate different properties, which suggests that each error signal drives modulatory changes independently, despite those changes ultimately affecting the same ion channels. Toxicity Profile of Artesunate in Rats and Dogs AO_SCPLOWBSTRACTC_SCPLOWO_ST_ABSOBJECTIVESC_ST_ABSThe objectives of these studies were to investigate the toxicity, safety and toxicokinetics of single and multiple doses of artesunate for injection in rats and dogs.

METHODSSprague-Dawley rats and Beagle dogs were treated intravenously or intramuscularly for 28 consecutive days with doses of up to 30 mg/kg artesunate, evaluating toxicity, kinetics, genotoxicity, and cardiovascular and central nervous safety parameters after single and 4-week repeated administrations. Furthermore, respiratory parameters were evaluated after a single intravenous administration in rats.

RESULTSArtesunate was well tolerated with no mortality and only minor effects on clinical pathology parameters. Following repeated intramuscular administration, local reactions at the injection site became evident.

CONCLUSIONSThe results obtained in these studies support the safe use of artesunate in humans.

Subchronic data for repeated intravenous and intramuscular artesunate in animals are available in the literature for a maximum duration of 14-days. Published results of genotoxicity studies in vitro in bacteria or cells of animal or human origin and in vivo in animals are equivocal, with a couple of negative results following oral administration and positive results in human cells.

With the studies performed in rats and dogs for a treatment period of 28-days including measures of safety pharmacology, toxicokinetic evaluation and in vivo genotoxicity data, a comparison of toxicity and exposure of intravenous versus intramuscular artesunate for a prolonged treatment period is possible, delivering important information, when a decision between both routes has to be made for use in patients. Results of the peripheral blood micronucleus assay in rats were positive, suggesting genotoxic potential for the intravenous route of administration of artesunate, what is, beside patients, likewise relevant for production and health care personnel. Shade-induced transcription of PIF-Direct-Target Genes precedes H3K4-trimethylation chromatin modification rises The phytochrome (phy)-PIF (Phytochrome Interacting Factor) sensory module perceives and transduces light signals to Direct-Target Genes (DTGs), which then drive the adaptational responses in plant growth and development, appropriate to the prevailing environment. These signals include the first exposure of etiolated seedlings to sunlight upon emergence from subterranean darkness, and the change in color of the light that is filtered through, or reflected from, neighboring vegetation ( shade). Previously, we identified three broad categories of rapidly signal-responsive genes: those repressed by light and conversely induced by shade; those repressed by light, but subsequently unresponsive to shade; and those responsive to shade only. Here, we investigate the potential role of epigenetic chromatin modifications in regulating these contrasting patterns of phy-PIF module-induced expression of DTGs. Using RNA-seq and ChlP-seq, time-resolved profiling of transcript and histone 3 lysine 4 trimethylation (H3K4me3) levels, respectively, we show that, whereas the initial dark-to-light transition triggers a rapid, apparently temporally-coincident decline of both parameters, the light-to-shade transition induces similarly rapid increases in transcript levels that precede increases in H3K4me3 levels. Together with other recent findings, these data raise the possibility that, rather than being causal in the shade-induced expression changes, H3K4me3 may function to buffer the rapidly fluctuating shade/light switching that is intrinsic to vegetational canopies under natural sunlight conditions. 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. Highly efficient synthetic CRISPR RNA/Cas9-based mutagenesis for rapid cardiovascular phenotypic screening in F0 zebrafish The zebrafish is a valuable vertebrate model to study cardiovascular formation and function due to the facile visualization and rapid development of the circulatory system in its externally growing embryos. Despite having distinct advantages, zebrafish have paralogs of many important genes, making reverse genetics approaches inefficient since generating animals bearing multiple gene mutations requires substantial efforts. Here, we present a simple and robust synthetic CRISPR RNA/Cas9-based mutagenesis approach for generating biallelic F0 zebrafish knockouts. Using a dual-guide synthetic CRISPR RNA/Cas9 ribonucleoprotein (dgRNP) system, we compared the efficiency of biallelic gene disruptions following the injections of one, two, and three dgRNPs per gene into the cytoplasm or yolk. We show that simultaneous cytoplasmic injections of three distinct dgRNPs per gene into one-cell stage embryos resulted in the most efficient and consistent biallelic gene disruptions. Importantly, this triple dgRNP approach enables efficient inactivation of cell autonomous and cell non-autonomous gene function, likely due to the low mosaicism of biallelic disruptions. In support of this finding, we provide evidence that the F0 animals generated by this method fully phenocopied the endothelial and peri-vascular defects observed in corresponding stable mutant homozygotes. Moreover, this approach faithfully recapitulated the trunk vessel phenotypes resulting from the genetic interaction between two vegfr2 zebrafish paralogs. Mechanistically, investigation of genome editing and mRNA decay indicates that the combined mutagenic actions of three dgRNPs per gene lead to an increased probability of frameshift mutations, enabling efficient biallelic gene disruptions. Therefore, our approach offers a highly robust genetic platform to quickly assess novel and redundant gene function in F0 zebrafish. Mitochondria preserve an autarkic one-carbon cycle to confer growth-independent cancer cell migration and metastasis Progression of primary cancer to metastatic disease is the most common cause of death in cancer patients with minimal treatment options available. Canonical drugs target mainly the proliferative capacity of cancer cells, which often leaves slow-proliferating, persistent cancer cells unaffected. Metabolic determinants that contribute to growth-independent functions supporting resistance and metastatic dissemination are still poorly understood.

In the present study, we revealed that antifolate treatment results in an uncoupled and autarkic mitochondrial one-carbon (1C) metabolism allowing sustained serine catabolism and formate overflow when cytosolic 1C metabolism is impaired. Interestingly, antifolate dependent growth-arrest did not correlate with decreased migration capacity. Therefore, using the antifolate Methotrexate as a tool compound allowed us to disentangle proliferation and migration to profile the metabolic phenotype of migrating (growth-arrested) cells. Supported by an increased NAD/NADH ratio, we observed increased serine de novo synthesis and increased serine catabolism to formate.

Consequently, inhibition of serine de novo synthesis using the competitive PHGDH-inhibitor BI-4916 or direct inhibition of mitochondrial 1C metabolism reduced cancer cell migration. Using an orthotopic breast cancer model, we show that sole inhibition of mitochondrial serine catabolism does not affect primary tumor growth but strongly inhibits pulmonary metastasis.

We conclude that mitochondrial 1C metabolism, despite being dispensable for proliferative capacities, confers an advantage to cancer cells by supporting their motility potential.

Our results improve our understanding of 1C metabolism and of metabolic determinants that support the process of cancer cell migration and metastasis. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. RNA polymerase II is recruited to DNA double-strand breaks for dilncRNA transcription in Drosophila DNA double-strand breaks are among the most toxic lesions that can occur in a genome and their faithful repair is thus of great importance. Recent findings have uncovered local transcription that initiates at the break and forms a non-coding transcript, called damage-induced long non-coding RNA (dilncRNA), which helps to coordinate the DNA transactions necessary for repair. We provide nascent RNA sequencing-based evidence that RNA polymerase II transcribes the dilncRNA in Drosophila and that this is more efficient for DNA breaks in an intron-containing gene, consistent with the higher damage-induced siRNA levels downstream of an intron. The spliceosome thus stimulates recruitment of RNA polymerase II to the break, rather than merely promoting the annealing of sense and antisense RNA to form the siRNA precursor. In contrast, RNA polymerase III nascent RNA libraries did not contain reads corresponding to the cleaved loci and selective inhibition of RNA polymerase III did not reduce the yield of damage-induced siRNAs. Finally, the damage-induced siRNA density was unchanged downstream of a T8 sequence, which terminates RNA polymerase III transcription. We thus found no evidence for a participation of RNA polymerase III in dilncRNA transcription in cultured Drosophila cells. 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. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Experimental characterization of in silico red-shift predicted iLOVL470T/Q489K and iLOVV392K/F410V/A426S mutants iLOV is a flavin mononucleotide-binding fluorescent protein used for in vivo cellular imaging similar to the green fluorescent protein. To expand the range of applications of iLOV, spectrally tuned red-shifted variants are desirable to have reduced phototoxicity and better tissue penetration. In this report, we experimentally tested two iLOV mutants, iLOVL470T/Q489K and iLOVV392K/F410V/A426S, which were previously computationally proposed by Khrenova et al. (DOI: 10.1021/acs.jpcb.7b07533) to have red-shifted excitation and emission spectra. We found that mutants bearing the V392K mutation lost the ability to bind FMN. While iLOVL470T/Q489K is about 20% brighter compared to WT iLOV in vitro, it exhibits a blue shift in contrast to QM/MM predictions. Additionally, both mutants are expressed at low levels and have undetectable fluorescence in living cells, which prevents their utilization in imaging applications. Our results suggest that rational fluorescent protein design efforts can benefit from in silico protein stability and ligand affinity calculations. NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action. Deficits in seizure threshold and other behaviors in adult mice without gross neuroanatomic injury after late gestation transient prenatal hypoxia Intrauterine hypoxia is a common cause of brain injury in children resulting in a broad spectrum of long-term neurodevelopmental sequela, including life-long disabilities that can occur even in the absence of severe neuroanatomic damage. Postnatal hypoxia-ischemia rodent models are commonly used to understand the effects of ischemia and transient hypoxia on the developing brain. Postnatal models, however, have some limitations. First, they do not test the impact of placental pathologies on outcomes from hypoxia. Second, they primarily recapitulate severe injury because they provoke substantial cell death, which is not seen in children with mild hypoxic injury. Lastly, they do not model preterm hypoxic injury. Prenatal models of hypoxia in mice may allow us to address some of these limitations to expand our understanding of developmental brain injury. The published rodent models of prenatal hypoxia employ multiple days of hypoxic exposure or complicated surgical procedures, making these models challenging to perform consistently in mice. Furthermore, large animal models suggest that transient prenatal hypoxia without ischemia is sufficient to lead to significant functional impairment to the developing brain. However, these large animal studies are resource-intensive and not readily amenable to mechanistic molecular studies. Therefore, here we characterized the effect of late gestation (embryonic day 17.5) transient prenatal hypoxia (5% inspired oxygen) on long-term anatomical and neurodevelopmental outcomes in mice. Late gestation transient prenatal hypoxia increased hypoxia-inducible factor 1 alpha protein levels (a marker of hypoxic exposure) in the fetal brain. Hypoxia exposure predisposed animals to decreased weight at postnatal day 2, which normalized by day 8. However, hypoxia did not affect gestational age at birth, litter size at birth, or pup survival. No differences in fetal brain cell death or long-term gray or white matter changes resulted from hypoxia. Animals exposed to prenatal hypoxia did have several long-term functional consequences, including sex-dichotomous changes. Hypoxia exposure was associated with a decreased seizure threshold and abnormalities in hindlimb strength and repetitive behaviors in males and females. Males exposed to hypoxia had increased anxiety-related deficits, whereas females had deficits in social interaction. Neither sex developed any motor or visual learning deficits. This study demonstrates that late gestation transient prenatal hypoxia in mice is a simple, clinically relevant paradigm for studying putative environmental and genetic modulators of the long-term effects of hypoxia on the developing brain. Population genetics reveals divergent lineages and ongoing hybridization in a declining migratory fish species complex Deciphering the effects of historical and recent demographic processes responsible for the spatial patterns of genetic diversity and structure is a key objective in evolutionary and conservation biology. Using population genetic analyses, we investigated the demographic history, the contemporary genetic diversity and structure, and the occurrence of hybridization and introgression of two species of anadromous fish with contrasting life history strategies and which have undergone recent demographic declines, the allis shad (Alosa alosa) and the twaite shad (Alosa fallax). We genotyped 706 individuals from 20 rivers and 5 sites at sea in Southern Europe at thirteen microsatellite markers. Genetic structure between populations was lower for the nearly semelparous species A. alosa, which disperses greater distances compared to the iteroparous species, A. fallax. Individuals caught at sea were assigned at the river level for A. fallax and at the region level for A. alosa. Using an approximate Bayesian computation framework, we inferred that the most likely long term historical divergence scenario between both species and lineages involved historical separation followed by secondary contact accompanied by strong population size decline. Accordingly, we found evidence for contemporary hybridization and bidirectional introgression due to gene flow between both species and lineages. Moreover, our results support the existence of at least one distinct species in the Mediterrannean sea: A. agone in Golfe du Lion area, and another divergent lineage in Corsica. Overall, our results shed light on the interplay between historical and recent demographic processes and life history strategies in shaping population genetic diversity and structure of closely related species. The recent demographic decline of these species populations and their hybridization should be carefully considered while implementing conservation programs. Muscle Activation Profile During Perturbed Walking is Modulated According to Body State During an unexpected loss of balance, avoiding a fall requires people to readjust their footing rapidly and effectively. We investigated the activation patterns of the ankle and knee muscles, and muscle fiber type recruitment resulting from unannounced, medio-lateral (i.e., right/left) horizontal-surface walking perturbations in twenty healthy adults (27.00{+/-}2.79 years, 10 females). Surface electromyography (sEMG) total spectral power for specific frequency bands (40-60Hz, 60-150Hz, 150-250Hz, 250-400Hz and 400-1000Hz), from tibialis anterior (TA) and vastus lateralis (VL) muscles were analyzed. Compared to non-perturbed walking, we found a significant increase in the total spectral power of lower-extremity muscles during the first 3 seconds after perturbation. When two feet were on the ground in time of perturbation we found a different muscle fiber type recruitment pattern between VL and TA muscles. This was not significant for perturbations implemented when one foot was on the ground. Our findings suggest that muscle operating frequency is modulated in real time to fit body state and somatosensory input in the context of functional goal requirements such as a rapid change of footing in response to unexpected loss of balance in single and double-support phases of gait.

New & NoteworthyTo study muscle spectral profiles in response to loss of balance, we investigated the dynamics of muscle spectral power changes, across different frequency bands after unannounced mechanical perturbations during walking. We showed increased activation of high-frequency motor units of the lower-limb muscles, subside 3 seconds after perturbation. Differences in power increase of specific frequency bands suggest that muscle activation is modulated in real time to fit body state in the context of functional goal requirements. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. Individual-based multiple-unit dissimilarity: novel indices and null model for assessing temporal variability in community composition Beta-diversity was originally defined spatially, i.e., as variation in community composition among sites in a region. However, the concept of beta-diversity has since been expanded to temporal contexts. This is referred to as "temporal beta-diversity", and most approaches are simply an extension of spatial beta-diversity. The persistence and turnover of individuals over time is a unique feature of temporal beta-diversity. Nakadai (2020) introduced the "individual-based beta-diversity" concept, and provided novel indices to evaluate individual turnover and compositional shift by comparing individual turnover between two periods at a given site. However, the proposed individual-based indices are applicable only to pairwise dissimilarity, not to multiple-temporal (or more generally, multiple-unit) dissimilarity. Here, individual-based beta-diversity indices are extended to multiple-unit cases. In addition, a novel type of random permutation criterion related to these multiple-unit indices for detecting patterns of individual persistence is introduced in the present study. To demonstrate the usage the properties of these indices compared to average pairwise measures, I applied them to a dataset for a permanent 50-ha forest dynamics plot on Barro Colorado Island in Panama. Information regarding "individuals" is generally missing from community ecology and biodiversity studies of temporal dynamics. In this context, the methods proposed here are expected to be useful for addressing a wide range of research questions regarding temporal changes in biodiversity, especially studies using traditional individual-tracked forest monitoring data. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. MODEL PREDICTIVE CONTROL OF CANCER CELLULAR DYNAMICS: A NEW STRATEGY FOR THERAPY DESIGN Recent advancements in Cybergenetics have led to the development of new computational and experimental platforms that enable to robustly steer cellular dynamics by applying external feedback control. Such technologies have never been applied to regulate intracellular dynamics of cancer cells. Here, we show in silico that adaptive model predictive control (MPC) can effectively be used to steer signalling dynamics in Non-Small Cell Lung Cancer (NSCLC) cells to resemble those of wild-type cells, and to support the design of combination therapies. Our optimisation-based control algorithm enables tailoring the cost function to force the controller to alternate different drugs and/or reduce drug exposure, minimising both drug-induced toxicity and resistance to treatment. Our results pave the way for new cybergenetics experiments in cancer cells for drug combination therapy design. α-Phenylalanyl tRNA synthetase competes with Notch signaling through its N-terminal domain The alpha subunit of the cytoplasmic Phenylalanyl tRNA synthetase (-PheRS, FARSA in humans) displays cell growth and proliferation activities and its elevated levels can induce cell fate changes and tumor-like phenotypes that do neither dependent on the canonical function of charging tRNAPhe with phenylalanine nor on stimulating general translation. In intestinal stem cells of Drosophila midguts, -PheRS levels are naturally slightly elevated and human FARSA mRNA levels are elevated in multiple cancers. In the Drosophila midgut model, elevated -PheRS levels caused the accumulation of many additional proliferating cells resembling intestinal stem cells (ISCs) and enteroblasts (EBs). This phenotype resembles partially the tumor-like phenotype described as Notch RNAi phenotype for the same cells. Genetic interactions between -PheRS and Notch suggest that their activities neutralize each other and that elevated -PheRS levels attenuate Notch signaling whether Notch induces differentiation into enterocytes, neuroblast stem cell proliferation, or transcription of a Notch reporter. These non-canonical functions all map to the N-terminal part of -PheRS which accumulates naturally in the intestine. This truncated version of -PheRS (-S) also localizes to nuclei and displays weak sequence similarity to the Notch intracellular domain (NICD), suggesting that -S might compete with the NICD for binding to a common target. Supporting this hypothesis, the tryptophan (W) residue reported to be key for the interaction between the NICD and the Su(H) BTD domain is not only conserved in -PheRS and -S, but also essential for attenuating Notch signaling.

Author SummaryAminoacyl tRNA synthetases charge tRNAs with their cognate amino acid to ensure proper decoding of the genetic code during translation. Independent of its aminoacylation function, the alpha subunit of Drosophila cytoplasmic Phenylalanyl tRNA synthetase (-PheRS, FARSA in humans) has an additional activity that promotes growth and proliferation (Ho et al., 2021). Here we describe that elevated -PheRS levels also induce cell fate changes and tumorous phenotypes in Drosophila midguts. Excessive proliferating cells with stem and progenitor cell characteristics accumulate and the composition of the terminally differentiated cells changes, too. This phenotype together with observed genetic interactions between -PheRS and Notch levels show that -PheRS counteracts Notch signaling in many different tissues and developmental stages. This novel activity of -PheRS maps to its N-terminal part, which is naturally produced. The fragment contains a DNA binding domain, translocates into nuclei, and displays essential similarities to a Notch domain that binds to the downstream transcription factor. This suggests that it might be competing with Notch for binding to a common target. Not only because Notch plays important roles in many tumors, but also because FARSA mRNA levels are considerably upregulated in many tumors, this novel activity deserves more attention for cancer research. A protease-activatable luminescent biosensor andreporter cell line for authentic SARS-CoV-2 infection Efforts to define serological correlates of protection against COVID-19 have been hampered by the lack of a simple, scalable, standardised assay for SARS-CoV-2 infection and antibody neutralisation. Plaque assays remain the gold standard, but are impractical for high-throughput screening. In this study, we show that expression of viral proteases may be used to quantitate infected cells. Our assays exploit the cleavage of specific oligopeptide linkers, leading to the activation of cell-based optical biosensors. First, we characterise these biosensors using recombinant SARS-CoV-2 proteases. Next, we confirm their ability to detect viral protease expression during replication of authentic virus. Finally, we generate reporter cells stably expressing an optimised luciferase-based biosensor, enabling viral infection to be measured within 24 h in a 96- or 384-well plate format, including variants of concern. We have therefore developed a luminescent SARS-CoV-2 reporter cell line, and demonstrated its utility for the relative quantitation of infectious virus and titration of neutralising antibodies.

Author summaryTechniques for measuring infection with SARS-CoV-2 in the laboratory are laborious and time-consuming, and different laboratories use different approaches. There is therefore no generally agreed way to quantitate neutralising antibodies against SARS-CoV-2, which block infection with the virus and protect people from COVID-19. In this study, we describe a new way to measure SARS-CoV-2 infection, which is much simpler and faster than existing methods. It relies on the production of a specific protease enzyme by the virus, which is able to cleave and activate an engineered protein biosensor in infected cells. This biosensor emits light in the presence of viral infection, and the amount of light released is used as a readout for the amount of infectious SARS-CoV-2 present. The signal is very sensitive, so the number of infected cells required is very small, and the method can be scaled-up to test many samples at once. In particular, we demonstrate how it can be used to detect different variants of SARS-CoV-2, and quantitate neutralising antibodies against these viruses. High-throughput variant detection using a color-mixing strategy Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single qPCR tube is not possible due to a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Here, we experimentally validate our novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously. The color-mixing strategy relies on a simple intra-tube assay that can probe for 15 variants as part of an inter-tube assay that can probe for an exponentially increased number of variants. The color-mixing strategy is achieved using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wildtype or variant templates. We used the color-mixing strategy to probe for 21 pathogenic mutations in thalassemia and to distinguish between heterozygous and homozygous variants in 6 tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, we used the same strategy to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in rpoB gene, plus five isoniazid-resistance mutations in inhA and katG genes.