Fast event-related mapping of population fingertip tuning properties in human sensorimotor cortex at 7T fMRI studies that investigate somatotopic tactile representations in the human cortex typically use either block or phase-encoded stimulation designs. Event-related (ER) designs allow for more flexible and unpredictable stimulation sequences than the other methods, but they are less efficient. Here we compared an efficiency-optimized fast ER design (2.8s average intertrial interval, ITI) to a conventional slow ER design (8s average ITI) for mapping voxelwise fingertip tactile tuning properties in the sensorimotor cortex of 6 participants at 7 Tesla. The fast ER design yielded more reliable responses compared to the slow ER design, but with otherwise similar tuning properties. Concatenating the fast and slow ER data, we demonstrate in each individual brain the existence of two separate somatotopically-organized representations of the fingertips, one in the primary somatosensory cortex (S1) on the post-central gyrus, and the other shared across the motor and pre-motor cortices on the pre-central gyrus. In both S1 and motor representations, fingertip selectivity decreased progressively, from narrowly-tuned Brodmann areas 3b and 4a respectively, towards associative parietal and frontal regions that responded equally to all fingertips, suggesting increasing information integration along these two pathways. In addition, fingertip selectivity in S1 decreased from the cortical representation of the thumb to that of the pinky.

Significance StatementSensory and motor cortices in the human brain contain map-like representations of the body in which adjacent brain regions respond to adjacent body parts. The properties of these somatotopic maps provide important insight into how tactile and motor information is processed by the brain. Here, we describe an efficient mapping method using functional magnetic resonance imaging (fMRI) to measure somatotopic maps and their tuning properties. We used a fast event-related sequence to map the five fingers of the left hand in six human participants, and show that this method is superior to a conventional, slower event-related design. Furthermore, we confirm previously-identified tuning properties of fingertip representations in somatosensory cortex, and reveal a hitherto unknown fingertip map in the motor cortex.
Cholecystokinin 1 Receptor (Cck1R) Normalizes mTORC1 signaling and is Protective to Purkinje cells in SCA Mice Spinocerebellar Ataxias (SCAs) are a group of genetic diseases characterized by progressive ataxia and neurodegeneration, often in cerebellar Purkinje neurons. A SCA1 mouse model, Pcp2-ATXN1[30Q]D776, has severe ataxia in absence of progressive Purkinje neuron degeneration and death. Previous RNA-seq analyses identified cerebellar up-regulation of the peptide hormone Cholecystokinin (Cck) in Pcp2-ATXN1[30Q]D776 mice. Importantly, absence of Cck1 receptor (Cck1R) in Pcp2-ATXN1[30Q]D776 mice confers a progressive disease with Purkinje neuron death. A Cck1R agonist, A71623 administered to Pcp2-ATXN1[30Q]D776;Cck-/- and Pcp2-AXTN1[82Q] mice dampened Purkinje neuron pathology and associated deficits in motor performance. In addition, A71623 administration improved motor performance of Pcp2-ATXN2[127Q] SCA2 mice. Moreover, the Cck1R agonist A71623 corrected mTORC1 signaling and improved expression of calbindin in cerebella of AXTN1[82Q] and ATXN2[127Q] mice. These results indicate that manipulation of the Cck-Cck1R pathway is a potential therapeutic target for treatment of diseases involving Purkinje neuron degeneration.
A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. Effect of Adolescent Stress on Adult Morphine-Induced Behavioral Sensitization is Dependent Upon Genetic Background 1Deaths related to opioid use have skyrocketed in the United States, leading to a public health epidemic. Research has shown that both biological (genes) and environmental (stress) precursors are linked to opioid use. In particular, stress during adolescence - a critical period of frontal lobe development - influences the likelihood of abusing drugs. However, little is known about the biological mechanisms through which adolescent stress leads to long-term risk of opioid use, or whether genetic background moderates this response. Male and female C57BL/6J and BALB/cJ mice were exposed to chronic variable social stress (CVSS) or control conditions throughout adolescence and then tested for morphine locomotor sensitization or morphine consumption in adulthood. To examine possible mechanisms that underlie stress-induced changes in morphine behaviors, we assessed physiological changes in response to acute stress exposure and prefrontal cortex miRNA gene expression. Adolescent stress did not influence morphine sensitization or consumption in BALB/cJ animals, and there was limited evidence of stress effects in female C57BL/6J mice. In contrast, male C57BL/6J mice exposed to adolescent CVSS had blunted morphine sensitization compared to control animals; no differences were observed in the acute locomotor response to morphine administration or morphine consumption. Physiologically, C57BL/6J mice exposed to CVSS had an attenuated corticosterone recovery following an acute stressor and downregulation of twelve miRNA in the prefrontal cortex compared to control mice. The specificity of the effects for C57BL/6J versus BALB/cJ mice provides evidence of a gene by environmental interaction influencing opioid behaviors. Long-term differences in stress reactivity or miRNA expression suggests two possible biological mechanisms to evaluate in future research. Temporally robust occupancy frequency distributions in riverine ecosystems explained by local biodiversity regulation The distribution of species site occupancy--the occupancy frequency distribution (OFD)--is a fundamental pattern in spatial ecology. Despite decades of research, the mechanisms responsible for the shape of the distribution remain incompletely understood. Here, we simultaneously examine both spatial and temporal patterns of site occupancy in communities of macroinverterbrates, macrophytes and diatoms. We find that the shape of the OFD at the catchment scale is typically time invariant despite substantial turnover in composition and species proportional occupancy over the 30 year time series. This can been seen as evidence in favour of an ecological mechanism for the shape of the OFD driven by local regulation of species richness and other biodiversity properties. We show that a simple, locally saturated patch occupancy model explains both the shape of the OFD and the temporal turnover rate in English rivers. This model is mathematically analogous to Hubbells neutral theory however concerns patches occupied by a given species rather its population size. Our analysis supports the view that metacommunities exist in a dynamic steady state arising from local ecological constraints and that this intrinsic regulation can drives the emergence of a wide variety of macroecological patterns. We also show that interspecific differences are unlikely to play a dominant role in driving large scale ecosystem structure. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. 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. A chromosome-level genome assembly of the European Beech (Fagus sylvatica) reveals anomalies for organelle DNA integration, repeat content and distribution of SNPs The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 Mb and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb of the genome in which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome- wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate. TNFa-induced LDL cholesterol accumulation involve elevated LDLR cell surface levels and SR-B1 downregulation in human arterial endothelial cellsan arterial endothelial cells Excess lipid droplets are frequently observed in arterial endothelial cells at sites of advanced atherosclerotic plaques. Here, the role of tumor necrosis factor alpha (TNF) in modulating low density lipoprotein (LDL) content in confluent primary human aortic endothelial cells (pHAECs) was investigated. TNF promoted up to 2 folds increase in cellular cholesterol, which was resistant to ACAT inhibition. The cholesterol increase was associated with increased 125I-LDL surface binding. Using the non-hydrolysable label, Dil, TNF could induce a massive increase in Dil-LDL by over 200 folds. The elevated intracellular Dil-LDL was blocked with excess unlabeled LDL and PCSK9, but not oxidized LDL (oxLDL), receptor associated protein (RAP), or apolipoprotein (apoE) depletion. Moreover, TNF-induced increase of LDL-derived lipids was elevated through lysosome inhibition. Using specific LDLR antibody, the Dil-LDL accumulation was reduced by over 99%. Effects of TNF included LDLR cell surface increase by 138%, and very large increases in ICAM-1 total and surface proteins, respectively. In contrast, that of scavenger receptor B1 (SR-B1) was reduced. Additionally, LDLR antibody bound rapidly in TNF-treated cells by about 30 folds, inducing a migrating shift in the LDLR protein. The effect of TNF on Dil-LDL accumulation was inhibited by the antioxidant tetramethythiourea (TMTU) dose-dependently, but not by inhibitors against NF-{kappa}B, stress kinases, ASK1, JNK, p38, or apoptosis caspases. Grown on transwell inserts, TNF did not enhance apical to basolateral LDL cholesterol or Dil release. It is concluded that TNF promotes LDLR functions through combined increase at the cell surface and SR-B1 downregulation. E-Cadherin mediated AMIS localisation Individual cells within de novo polarising tubes and cavities must integrate their forming apical domains into a centralised apical membrane initiation site (AMIS). This is necessary to enable organised lumen formation within multi-cellular tissue. Despite the well documented importance of cell division in localising the AMIS, we have found a division-independent mechanism of AMIS localisation that relies instead on Cadherin-mediated cell-cell adhesion. Our study of de novo polarising mouse embryonic stem cells (mESCs) cultured in 3D suggest that cell-cell adhesion directs the localisation of apical proteins such as PAR-6 to a centralised AMIS. Unexpectedly, we also found that mESC cell clusters lacking functional E-Cadherin were still able to form a lumen-like cavity in the absence of AMIS localisation and did so at a later stage of development via a  closure mechanism, instead of via hollowing. This work suggests that there are two, interrelated mechanisms of apical polarity localisation: cell adhesion and cell division. Alignment of these mechanisms in space allows for redundancy in the system and ensures the development of a coherent epithelial structure within a growing organ. De novo identification of maximally deregulated subnetworks based on multi-omics data with DeRegNet BackgroundWith a growing amount of (multi-)omics data being available, the extraction of knowledge from these datasets is still a difficult problem. Classical enrichment-style analyses require predefined pathways or gene sets that are tested for significant deregulation to assess whether the pathway is functionally involved in the biological process under study. De novo identification of these pathways can reduce the bias inherent in predefined pathways or gene sets. At the same time, the definition and efficient identification of these pathways de novo from large biological networks is a challenging problem.

ResultsWe present a novel algorithm, DeRegNet, for the identification of maximally deregulated subnetworks on directed graphs based on deregulation scores derived from (multi-)omics data. DeRegNet can be interpreted as maximum likelihood estimation given a certain probabilistic model for de-novo subgraph identification. We use fractional integer programming to solve the resulting combinatorial optimization problem. We can show that the approach outperforms related algorithms on simulated data with known ground truths. On a publicly available liver cancer dataset we can show that DeRegNet can identify biologically meaningful subgraphs suitable for patient stratification. DeRegNet is freely available as open-source software.

ConclusionThe proposed algorithmic framework and its available implementation can serve as a valuable heuristic hypothesis generation tool contextualizing omics data within biomolecular networks. An energy-dissipative state of the major antenna complex of plants Plants and algae are faced with a conundrum: harvesting sufficient light to drive their metabolic needs while dissipating light in excess to prevent photodamage, a process known as non-photochemical quenching. A slowly relaxing form of energy dissipation, termed qH, is critical for plants survival under abiotic stress. Here, we tested whether we could isolate photosynthetic subcomplexes (from plants in which qH was induced) that would remain in an energy-dissipative state. Interestingly chlorophyll fluorescence lifetimes were decreased by qH in isolated major trimeric antenna complexes, providing a natively quenched complex with physiological relevance to natural conditions. Next, we monitored the changes in thylakoid pigment, protein or lipid content of antenna with active or inactive qH, and no evident differences were detected. Finally, we investigated whether specific antenna subunits of the major antenna were required for qH but found it insensitive to trimer composition. Because qH can occur in the absence of specific xanthophylls, and no changes in pigments were detected, we propose that the energy-dissipative state reported here may stem from chlorophyll-chlorophyll excitonic interaction. PAK1 and PAK2 in cell metabolism regulation P21-activated kinases (PAKs) regulate processes associated with cytoskeletal rearrangements, such as cell division, adhesion, and migration. The possible regulatory role of PAKs in cell metabolism has not been well explored, but increasing evidence suggests that a cell metabolic phenotype is related to cell interactions with the microenvironment. We analyzed the impact of PAK inhibition by small molecule inhibitors, small interfering RNA, or gene knockout on the rates of mitochondrial respiration and aerobic glycolysis. Pharmacological inhibition of PAK group I by IPA-3 induced a strong decrease in metabolic rates in human adherent cancer cell lines, leukemia/lymphoma cell lines, and primary leukemia cells. The immediate effect of FRAX597, which inhibits PAK kinase activity, was moderate, indicating that PAK nonkinase functions are essential for cell metabolism. Selective downregulation or deletion of PAK2 was associated with a shift toward oxidative phosphorylation. Alternatively, PAK1 knockout resulted in increased glycolysis. However, the overall metabolic capacity was not substantially reduced by PAK1 or PAK2 deletion, possibly due to partial redundancy in PAK1/PAK2 regulatory roles or to activation of other compensatory mechanisms. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Clonal Dynamics of Cancer-Initiating Cells in a Mouse Genetic Mosaic Model of Serous Ovarian Cancer High-grade serous ovarian cancer (HGSOC) is the most lethal gynecologic malignancy. Early detection remains the most promising strategy to decrease mortality and requires a deep understanding of the cancer initiation process. Recent breakthroughs identified the fallopian tube (FT) Pax8+ cells as a cell-of-origin for HGSOC. However, the progressive trajectory from an initiated Pax8+ cell toward malignancy remains elusive. Here, we use a mouse genetic mosaic system called MADM (Mosaic Analysis with Double Markers), which generates sporadic mutant cells labeled with GFP, to trace the clonal dynamics of HGSOC development. Surprisingly, we observed dichotomous clonal expansion kinetics of individual cancer-initiating Pax8+ cells; while the vast majority stalled immediately, a few significantly expanded. Clonal analysis of wildtype cells showed concordant dichotomy, suggesting that a small fraction of Pax8+ cells intrinsically carry higher expansion potential. Furthermore, although all wildtype clones stalled shortly after development, some mutant clones continued to progress, indicating that oncogenic mutations prolong clonal expansion beyond the cell-intrinsic potential. Finally, intra-clonal fate mapping revealed that progressing mutant clones have a higher propensity to produce Pax8+ cells, linking imbalanced cell differentiation to disease progression. In summary, our studies reveal that the earliest clonal expansion during HGSOC initiation emerges when oncogenic mutations exaggerate the intrinsic expansion potential of a rare Pax8+ subset and bias its differentiation potential, warranting future investigation as a readout for early detection and cancer prevention.

SIGNIFICANCE STATEMENTEffective early detection is the most promising approach to manage cancer, but it requires a thorough understanding of cancer initiation and early progression. Serous ovarian cancer has been found to originate from Pax8+ cells in the fallopian tube. To assess the clonal dynamics of cancer-initiating cells, we established a mouse genetic mosaic model that enables clonal analysis of mutant Pax8+ cells in the fallopian tube. We uncovered that the earliest clonal event in ovarian cancer development emerges from a rare subset of stem/progenitor-like cells with intrinsic expansion potential, and that oncogenic mutations exaggerate and prolong this potential to promote cancer initiation. Our findings pinpoint a precise target for future early-detection and cancer prevention studies.

HIGHLIGHTSO_LIOvarian cancer initiation revealed by a mouse genetic mosaic model.
C_LIO_LIDichotomous expansion capacity of initiated cells at the clonal level.
C_LIO_LISynergies between cell-intrinsic properties and onco-mutations fuels clonal expansion.
C_LIO_LIImbalanced differentiation within highly expansive clones leads to Pax8 clusters.
C_LI Deep autoencoder for interpretable tissue-adaptive deconvolution and cell-type-specific gene analysis Single-cell RNA-sequencing (RNA-seq) has become a powerful tool to study biologically significant characteristics at explicitly high resolution. However, its application on emerging data is currently limited by its intrinsic techniques. Though many methods have been proposed to analyze bulk data using single-cell profile as a reference, they are limited on the interpretability, processing speed, and data size requirement. Here, we introduce Tissue-AdaPtive autoEncoder (TAPE), a deep learning method connecting bulk RNA-seq and single-cell RNA-seq, to achieve precise prediction in short time. By constructing an interpretable decoder and training under a unique scheme, TAPE can predict cell-type fractions and cell-type-specific gene expression tissue-adaptively. Compared with existing methods on several benchmarking datasets, TAPE is more accurate (up to 40% improvement on the real bulk data) and faster. Thus, it is sensitive enough to provide biologically meaningful predictions. For example, only TAPE can predict the tendency of increasing monocytes-to-lymphocytes ratio in COVID-19 patients from mild to serious symptoms, of which estimated indices are consistent with laboratory data. More importantly, through the analysis of clinical data, TAPE shows its ability to predict cell-type-specific gene expression profiles with biological significance. Combining with single-sample gene set enrichment analysis (ssGSEA), TAPE also provides valuable clues to investigate the immune response in different virus-infected patients. We believe that TAPE will enable and accelerate the precise analysis of high-throughput clinical data in a wide range. 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. Using Obviously-Related Instrumental Variables to Increase the Predictive Power of Polygenic Scores The conventional way of boosting the predictive power of polygenic scores is to increase the GWAS sample size by meta-analyzing GWAS results of multiple cohorts. In this paper, we challenge this convention. Through simulations, we show that Instrumental Variable (IV) regression using two polygenic scores constructed from independent GWAS summary statistics outperforms the typical Ordinary Least Squares (OLS) model employing a single meta-analysis based polygenic score in terms of bias, root mean squared error, and statistical power. We verify the empirical validity of the simulations by predicting educational attainment (EA) and height in a sample of siblings from the UK Biobank. We show that IV regression between-families approaches the SNP-based heritability, and improves the predictive power of polygenic scores by 12% (height) to 22% (EA). Furthermore, IV regression within-families provides the tightest lower bound for the direct genetic effect, increasing the lower bound for EA from 2.0% to 3.4%, and for height from 28.9% to 37.7%. Engineered bacteria detect tumor DNA Advances in bacterial engineering have catalysed the development of living cell diagnostics and therapeutics1-3, including microbes that respond to diseases such as gut inflammation4, intestinal bleeding5, pathogens6 and hypoxic tumors7. Bacteria can easily access the entire gastrointestinal tract via oral administration8, and they can produce outputs that can be noninvasively measured in stool4 or urine7. Cellular memory, such as bistable switches4,9,10 or genomic rearrangement11, has been used to allow bacteria to store information over time. However, living biosensors have not yet been engineered to detect specific DNA sequences or mutations from outside the cell. Here, we engineer naturally competent Acinetobacter baylyi to detect donor DNA from the genomes of colorectal cancer (CRC) cells and organoids. We characterize the functionality of the biosensors in vitro with co-culture assays and then validate in vivo with sensor bacteria delivered orally or rectally into mice injected with orthotopic donor CRC organoids. We observe horizontal gene transfer from the tumor to the sensor bacteria in vivo, allowing their detection in stool. The sensor bacteria achieved 100% discrimination between mice with and without CRC using both delivery methods. Our findings establish a framework for biosensing applications that require the detection of mutations or organisms within environments that are difficult to sample. In addition, the platform can be readily expanded to include in situ production and delivery of therapeutic payloads at the detection site. A novel method for producing functionalized vesicles that efficiently deliver oligonucleotides in vitro in cancer cells and in vivo in mice Nano-based delivery systems have enhanced our ability to administer and target drugs and macromolecules to their targets. Oligonucleotide drugs have great therapeutic potential but often have off-target effects and stability issues. Therefore, they are often encapsulated in vesicles with targeting ligands such as antibodies (Ab) to deliver their cargo. Herein, we describe a novel, scalable and straightforward approach to producing functionalized vesicles called the "Functionalized Lipid Insertion Method," which differs from the older "Detergent-Dialysis Method." The latter method required excess detergent and extensive dialysis over many hours to produce the functionalized vesicles. With our method, only the functionalized lipid is detergent-solubilized during the engineering of the vesicle. The approach reduces the dialysis time, keeps the vesicles intact while orienting the targeting moieties of the functionalized lipid toward the outside of the vesicle. Pilot in vitro and in vivo experiments was performed to show the feasibility of our method. Dynamic light scattering (DLS) experiments suggested that the original vesicular structure was relatively unperturbed, and the functionalized lipid was inserted externally. Our approach efficiently delivered oligonucleotides and affected the function of liver cancer HepG2 cells. Furthermore, functionalized vesicles achieved targeted delivery of oligonucleotides in mice without inducing a significant innate immune response. The industrial and therapeutic significance and implications of functionalized vesicles produced by our method are also discussed. Additional experiments and analyses are recommended to bring out the full potential of this molecular delivery technology. Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex Corona Virus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) has caused a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha variant, B.1.1.7 lineage, reported from the United Kingdom (UK), is of great significance in that it contains several mutations that increase its infection and transmission rates as evident from clinical reports. Specifically, the N501Y mutation in the SARS-CoV-2 spike S1 receptor binding domain (S1-RBD) has been shown to possess an increased affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased affinity using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing drugs targeting the ACE2-S1-RBD interaction.

SignificanceThe emergence of the new SARS-CoV-2 lineage in the UK in December 2020 has further aggravated the COVID-19 pandemic due to an increased ability of the variant to infect human hosts, likely due to mutations in the viral S1 spike protein including the N501Y S1-RBD mutation that is located at the interface of S1-RBD and ACE2, the host cell receptor for SARS-CoV-2. Given its location at the interface, N501Y S1-RBD mutation can therefore potentially alter the interfacial interaction. Multiple, all-atom, explicit solvent MD simulations of the ACE2-S1-RBD complex carried here indicated a more stable interaction between the N501Y mutant S1-RBD and ACE2 through stabilizing interfacial interactions of residues at one end of the interface that are either sequentially or physically near the mutation site. These mechanistic details will aid in better understanding the mechanism by which the alpha variant has increased infectivity as well as in designing better therapeutics including ACE2-S1 spike protein inhibitors that will, in turn, help thwarting the current and future pandemic.

HighlightsO_LIN501 in the wild type SARS-CoV-2 S1-RBD forms unsustained hydrogen bonds with residues in the ACE2, namely Y41 and K353
C_LIO_LIY501 in the N501Y mutant SARS-CoV-2 S1-RBD is not capable of forming substantial hydrogen bonds with ACE2 within the time span of the current simulation
C_LIO_LIEvidence from analyzing the simulation results suggests that Y501 of S1-RBD could form other types of non-covalent interactions with ACE2, such as van der Waals interactions
C_LIO_LIN501Y S1-RBD mutation stabilizes the position of interfacial residues neighboring to the mutation site, as well as other non-interfacial residues that are distant from the mutation site
C_LIO_LIThese altered dynamics results in more stable interaction of S1-RBD with ACE2 which could be the main reason underlying the reported enhanced affinity of S1-RBD in the SARS-CoV-2 alpha variant (UK B.1.1.7 lineage) to ACE2
C_LI REGEN-COV protects against viral escape in preclinical and human studies Monoclonal antibodies against SARS-CoV-2 are a clinically validated therapeutic option against COVID-19. As rapidly emerging virus mutants are becoming the next major concern in the fight against the global pandemic, it is imperative that these therapeutic treatments provide coverage against circulating variants and do not contribute to development of treatment emergent resistance. To this end, we investigated the sequence diversity of the spike protein and monitored emergence of minor virus variants in SARS-COV-2 isolates found in COVID-19 patients or identified from preclinical in vitro and in vivo studies. This study demonstrates that a combination of non-competing antibodies, REGEN-COV, not only provides full coverage against current variants of concern/interest but also protects against emergence of new such variants and their potential seeding into the population in a clinical setting. A bacterial extracellular vesicle-based intranasal vaccine against SARS-CoV-2 protects against disease and elicits neutralizing antibodies to wild-type virus and Delta variant Several vaccines have been introduced to combat the coronavirus infectious disease-2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Current SARS-CoV-2 vaccines include mRNA-containing lipid nanoparticles or adenoviral vectors that encode the SARS-CoV-2 Spike (S) protein of SARS-CoV-2, inactivated virus, or protein subunits. Despite growing success in worldwide vaccination efforts, additional capabilities may be needed in the future to address issues such as stability and storage requirements, need for vaccine boosters, desirability of different routes of administration, and emergence of SARS-CoV-2 variants such as the Delta variant. Here, we present a novel, well-characterized SARS-CoV-2 vaccine candidate based on extracellular vesicles (EVs) of Salmonella typhimurium that are decorated with the mammalian cell culture-derived Spike receptor-binding domain (RBD). RBD-conjugated outer membrane vesicles (RBD-OMVs) were used to immunize the golden Syrian hamster (Mesocricetus auratus) model of COVID-19. Intranasal immunization resulted in high titers of blood anti-RBD IgG as well as detectable mucosal responses. Neutralizing antibody activity against wild-type and Delta variants was evident in all vaccinated subjects. Upon challenge with live virus, hamsters immunized with RBD-OMV, but not animals immunized with unconjugated OMVs or a vehicle control, avoided body mass loss, had lower virus titers in bronchoalveolar lavage fluid, and experienced less severe lung pathology. Our results emphasize the value and versatility of OMV-based vaccine approaches. 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. Prenatal environmental conditions underlie alternative reproductive tactics that drive the formation of a mixed-kin cooperative society Although animal societies often evolve due to limited natal dispersal that results in kin clustering and facilitates cooperation among relatives, many species form cooperative groups with low and variable kin structure. Groups in such mixed-kin societies often comprise residents and immigrants of the same sex that compete for breeding opportunities. To understand how such mixed-kin societies form despite their potential for social conflict, we investigated the environmental causes and subsequent fitness consequences of dispersal decisions in male cooperatively breeding superb starlings (Lamprotornis superbus) that live in a climatically unpredictable environment. We show that the two alternative reproductive tactics--natal dispersal or philopatry--exhibit reproductive tradeoffs resulting in equal lifetime inclusive fitness. The tactic an individual adopts is governed by the environment its parents experience prior to laying rather than the environment it experiences during its juvenile stage. When individuals adopt the tactic not predicted by early life environmental conditions, their fitness is reduced. Ultimately, climate-driven oscillating selection may help stabilize mixed-kin animal societies, despite their reduced kin structure and potential for social conflict.