LEA motifs promote desiccation tolerance in vivo BackgroundCells and organisms typically cannot survive in the absence of water. However, there are some notable exceptions, including animals such as nematodes, tardigrades, rotifers, and some arthropods. One class of proteins known to play a role in desiccation resistance is the late embryogenesis abundant (LEA) proteins. These largely disordered proteins protect plants and animals from desiccation. A multitude of studies have characterized stress-protective capabilities of LEA proteins in vitro and in heterologous systems. However, the extent to which LEA proteins exhibit such functions in vivo, in their native contexts in animals, is unclear.

Furthermore, little is known about the distribution of LEA proteins in multicellular organisms or tissue-specific requirements in conferring stress protection.

ResultsTo study the endogenous function of an LEA protein in an animal, we created a true null mutant of C. elegans LEA-1, as well as endogenous fluorescent reporters of the protein. We confirmed that C. elegans lacking LEA-1 are sensitive to desiccation. LEA-1 mutant animals were also sensitive to heat and osmotic stress and were prone to protein aggregation. During desiccation, LEA-1 expression increased and became more widespread throughout the body. LEA-1 was required at high levels in body wall muscle for animals to survive desiccation and osmotic stress. We identified minimal motifs within C. elegans LEA-1 that are sufficient to increase desiccation survival of E. coli. To test whether such motifs are central to LEA-1s in vivo functions, we then replaced the sequence of lea-1 with these minimal motifs and found that C. elegans survived mild desiccation and osmotic stress at the same levels as worms with the full-length protein.

ConclusionsOur results provide insights into the endogenous functions and expression dynamics of an LEA protein in a multicellular animal. The results show that LEA-1 buffers animals from a broad range of stresses. Our identification of LEA motifs that can function in both bacteria and in a multicellular organism suggests the possibility of engineering LEA-1-derived peptides for optimized desiccation protection.
ConVarT: a search engine for matching human genetic variants with variants from non-human species The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely "matching variants" (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and C. elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease.
The Point Mutations in SARS-CoV-2 Variants Induce Long-Range Dynamical Perturbations in Neutralizing Antibodies Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. Identification of differentially distributed gene expression and distinct sets of cancer-related genes identified by changes in mean and variability There is increasing evidence that changes in the variability or overall distribution of gene expression are important both in normal biology and in diseases, particularly cancer. Genes whose expression differs in variability or distribution without a difference in mean are ignored by traditional differential expression-based analyses. Using a Bayesian hierarchical model that provides tests for both differential variability and differential distribution for bulk RNA-seq data, we report here an investigation into differential variability and distribution in cancer. Analysis of eight paired tumour-normal datasets from The Cancer Genome Atlas confirms that differential variability and distribution are able to identify cancer-related genes. We further demonstrate that differential variability identifies cancer-related genes that are missed by differential expression analysis, and that differential expression and differential variability identify functionally distinct sets of genes. These results suggest that differential variability analysis may provide insights into genetic aspects of cancer that would not be revealed by differential expression, and that differential distribution analysis may allow for more comprehensive identification of cancer-related genes than analyses based on changes in mean or variability alone. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis 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. An RNA-based feed-forward mechanism ensures motor switching in oskar mRNA transport Regulated recruitment and activity of motor proteins is essential for intracellular transport of cargoes, including ribonucleoprotein complexes (RNP). Here we show that orchestration of oskar RNP transport in the Drosophila germline relies on the interplay of two double-stranded RNA binding proteins, Staufen and the dynein adaptor Egalitarian (Egl). Following delivery of nurse cell-synthesized oskar mRNA into the oocyte by dynein, enhanced recruitment of Staufen to the RNPs results in dissociation of Egl and reduced minus end-directed transport without overtly affecting overall dynein levels on the cargo. This mechanism allows kinesin-1 to translocate the mRNA to its final destination at the posterior pole of the oocyte. Using in vitro motility assays with purified components, we find that Staufen directly antagonizes Egl-mediated activation of oskar mRNA transport by dynein. We additionally show that Egl associates with staufen mRNA in the nurse cells and mediates enrichment of the mRNA and its protein product in the ooplasm. Our observations identify a feed-forward mechanism, whereby dynein-dependent accumulation of staufen mRNA and protein in the oocyte enables motor switching during oskar mRNA transport by downregulating dynein activity. Semi-automated analysis of an optical ATP indicator in neurons SignificanceThe firefly enzyme luciferase has been used in a wide range of biological assays, including bioluminescence imaging of ATP. The biosensor Syn-ATP utilizes subcellular targeting of luciferase to nerve terminals for optical measurement of ATP in this compartment. Manual analysis of Syn-ATP signals is challenging due to signal heterogeneity and cellular motion in long imaging sessions. Here, we have leveraged machine learning tools to develop a method for analysis of bioluminescence images.

AimOur goal was to create a semi-automated pipeline for analysis of bioluminescence imaging to improve measurements of ATP content in nerve terminals.

ApproachWe developed an image analysis pipeline that applies machine learning toolkits to distinguish neurons from background signals, and excludes neural cell bodies, while also incorporating user input.

ResultsSide-by-side comparison of manual and semi-automated image analysis demonstrated that the latter improves precision and accuracy of ATP measurements.

ConclusionsOur method streamlines data analysis and reduces user-introduced bias, thus enhancing the reproducibility and reliability of quantitative ATP imaging in nerve terminals. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. A bispecific monomeric nanobody induces spike trimer dimers and neutralizes SARS-CoV-2 in vivo Antibodies binding to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike have therapeutic promise, but emerging variants show the potential for virus escape. This emphasizes the need for therapeutic molecules with distinct and novel neutralization mechanisms. Here we isolated a nanobody that interacts simultaneously with two RBDs from different spike trimers of SARS-CoV-2, rapidly inducing the formation of spike trimer-dimers leading to the loss of their ability to attach to the host cell receptor, ACE2. We show that this nanobody potently neutralizes SARS-CoV-2, including the B.1.351 variant, and cross-neutralizes SARS-CoV. Furthermore, we demonstrate the therapeutic potential of the nanobody against SARS-CoV-2 and the B.1.351 variant in a human ACE2 transgenic mouse model. This naturally elicited bispecific monomeric nanobody establishes a novel strategy for potent inactivation of viral antigens and represents a promising antiviral against emerging SARS-CoV-2 variants. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Robust and Remote Measurement of Heart Rate Based on A Surveillance Camera Arrhythmia is a marked symptom of many cardiovascular diseases (CVDs). The accurate and timely detection of heart rate can greatly reduce the harm of arrhythmia to people. However, it is still a challenge to robustly and remotely measure heart rate in daily life due to the changing environmental conditions during measurement, such as the varying light intensity, the movement of people, and the uncertain distance between the sensor to people. In this study, we propose a method to accurately measure human heart rate within a distance of 4.5 meters under different light intensities by simply using a surveillance camera. After a 20-second color video of a persons hand is captured by the camera, a method based on the Fast Fourier Transform (FFT) algorithm is designed to extract the blood volume pulse wave to calculate the heart rate. According to the comparison between the real heart rate and results measured by electrocardiography (ECG), the proposed method achieves an accuracy of 95.8% when the measurement is performed within a distance of 4.5 meters and 90% when within 5.0 meters. Our experiments show that when the illuminance varies between 100-1000 lux (lighting level indoor), we still get the correct results. Our experiments also demonstrate that the proposed method accurately obtain heart rate even when the light intensity is below 32 lux (300-500 lux in a workplace environment). The methods strong adaptability to changing environmental conditions makes it applicable to many scenarios, such as homes of the elderly, classrooms, and other public spaces. 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. Clofoctol inhibits SARS-CoV-2 replication and reduces lung pathology in mice Drug repurposing has the advantage of shortening regulatory preclinical development steps. Here, we screened a library of drug compounds, already registered in one or several geographical areas, to identify those exhibiting antiviral activity against SARS-CoV-2 with relevant potency. Of the 1,942 compounds tested, 21 exhibited a substantial antiviral activity in Vero-81 cells. Among them, clofoctol, an antibacterial drug used for the treatment of bacterial respiratory tract infections, was further investigated due to favorable safety profile and pharmacokinetic properties. Notably, the peak concentration of clofoctol that can be achieved in human lungs is more than 20 times higher than its IC50 measured against SARS-CoV-2 in human pulmonary cells. This compound inhibits SARS-CoV-2 at a post-entry step. Lastly, therapeutic treatment of human ACE2 receptor transgenic mice decreased viral load, reduced inflammatory gene expression and lowered pulmonary pathology. Altogether, these data strongly support clofoctol as a therapeutic candidate for the treatment of COVID-19 patients.

SummaryAntivirals targeting SARS-CoV-2 are sorely needed. In this study, we screened a library of drug compounds and identified clofoctol as an antiviral against SARS-CoV-2. We further demonstrated that, in vivo, this compound reduces inflammatory gene expression and lowers pulmonary pathology. Targeting SARS-CoV-2 infection through CAR-T like bispecific T cell engagers incorporating ACE2 Despite advances in antibody treatments and vaccines, COVID-19 caused by SARS-CoV-2 infection remains a major health problem resulting in excessive morbidity and mortality and the emergence of new variants has reduced the effectiveness of current vaccines. Here, as a proof-of-concept we engineered primary CD8 T cells to express SARS-CoV-2 Spike protein-specific CARs, using extracellular region of ACE2, and demonstrated their highly specific and potent cytotoxicity towards Spike-expressing target cells. To improve on this concept as a potential therapeutic, we developed a bispecific T cell engager combining ACE2 with an anti-CD3 scFv (ACE2-Bite) to target infected cells and the virus. Similar to CAR-T cell approach, ACE2-Bite endowed cytotoxic cells to selectively kill Spike-expressing targets. Furthermore, ACE2-Bite neutralized the pseudoviruses of SARS-CoV, SARS-CoV-2 wild-type and variants including Delta and Omicron, as a decoy protein. Remarkably, ACE2-Bite molecule showed a higher binding and neutralization affinity to Delta and Omicron variants compared to SARS-CoV-2 wild-type Spike proteins, suggesting the potential of this approach as a variant-proof, therapeutic strategy for future SARS-CoV-2 variants, employing both humoral and cellular arms of the adaptive immune response. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. MYB insufficiency disrupts proteostasis in hematopoietic stem cells leading to age-related neoplasia The Myb transcription factor plays critical roles in normal and malignant hematopoiesis. Acquired genetic dysregulation of Myb, which plays a central role in hematopoietic stem cell (HSC) gene regulation, is involved in the etiology of a number of leukemias. Also, inherited non-coding variants of the Myb gene are a factor in susceptibility to many hematological conditions, including myeloproliferative neoplasms (MPN), but the mechanisms by which variations in Myb levels predispose to disease, including age-dependency in disease occurrence, are completely unknown. Here, we address these key points by showing that Myb insufficiency in mice leads in later life to MPN, myelodysplasia, and leukemia, mirroring the age profile of equivalent human diseases. This age-dependence is intrinsic to HSC, involving progressive accumulation of subtle changes. Interestingly, and linking to previous studies showing the importance of proteostasis to the maintenance of normal HSC, we observed altered proteosomal activity in young Myb-insufficient mice and later elevated ribosome activity. We propose that these alterations collectively cause an imbalance in proteostasis, potentially creating a cellular milieu favoring disease initiation by driver mutations. Hydrophobic mismatch is a key factor in protein transport on the Tat pathway The twin-arginine translocation (Tat) pathway transports folded proteins across membranes in bacteria, thylakoid, plant mitochondria, and archaea. In most species, the active Tat machinery consists of three independent subunits, TatA, TatB and TatC. TatA and TatB from all bacterial species possess short transmembrane alpha-helices (TMHs), both of which are only fifteen residues long in E. coli. Such short TMHs cause a hydrophobic mismatch between Tat subunits and the membrane bilayer. Here, by modifying the length of the TMHs of E. coli TatA and TatB, we assess the functional importance of the hydrophobic mismatch in the Tat transport mechanism. Surprisingly, both TatA and TatB with as few as 11 residues in their respective TMHs are still able to insert into the membrane bilayer, albeit with a decline in membrane integrity. Three different assays, both qualitative and quantitative, were conducted to evaluate the Tat activity of the TMH length mutants. Our experiments indicate that the TMHs of TatA and TatB appear to be evolutionarily tuned to 15 amino acids, with activity dropping off with any modification of this length. We believe our study supports a model of Tat transport utilizing localized toroidal pores that form when the membrane bilayer is thinned to a critical threshold. In this context, the 15-residue length of the TatA and TatB TMHs can be seen as a compromise between the need for some hydrophobic mismatch to allow the membrane to reversibly reach the threshold thinness required for toroidal pore formation, and the permanently destabilizing effect of placing even shorter helices into these energy-transducing membranes. CDK9 and PP2A regulate RNA polymerase II transcription termination and coupled RNA maturation. CDK9 is a critical kinase required for the productive transcription of protein-coding genes by RNA polymerase II (pol II). As part of P-TEFb, CDK9 phosphorylates the carboxyl-terminal domain (CTD) of pol II and elongation factors, including SPT5, which allows pol II to elongate past the early elongation checkpoint (EEC) encountered soon after initiation. We show that, in addition to halting pol II at the EEC, loss of CDK9 activity causes premature termination of transcription across the last exon, loss of polyadenylation factors from chromatin, and loss of polyadenylation of nascent transcripts. Inhibition of the phosphatase PP2A abrogates the premature termination and loss of polyadenylation caused by CDK9 inhibition, indicating that this kinase/phosphatase pair regulates transcription elongation and RNA processing at the end of protein-coding genes. Our phosphoproteomic analyses after CDK9 inhibition, using either DRB or an analog-sensitive CDK9 cell line, confirm the splicing factor SF3B1 as an additional key target of this kinase. These results emphasize the important roles that CDK9 plays in coupling transcription elongation and termination to RNA maturation downstream of the EEC. Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer. Including imprecisely georeferenced specimens improves accuracy of species distribution models and estimates of niche breadth: Don't let the perfect be the enemy of the good Museum and herbarium specimen records are frequently used to assess species conservation status and responses to climate change. Typically, records of occurrence with imprecise geolocality information are discarded because they cannot be matched confidently to environmental conditions, and are thus expected to increase uncertainty in downstream analyses. However, using only those records that can be precisely georeferenced risks undersampling of species environmental and geographic distributions. Using simulated and real species, we compared the effect of discarding versus retaining imprecise records on the accuracy of ecological niche models (ENMs), and estimates of niche breadth and extent of occurrence. Imprecise records were assigned to locations or climate values using a conservative approach. For simulated species, including imprecisely-georeferenced records alongside precisely-georeferenced records improved the accuracy of ENMs projected to the present and the future. Improvements were especially notable for species with <~20 precise occurrences, but accuracy continued to improve and eventually matched the accuracy of ENMs using all records with no georeferencing error. Using only precise records underestimated loss in suitable habitat and overestimated the amount of suitable habitat in both the present and future. Including imprecise records also improved univariate and multivariate estimates of niche breadth, and of species extent of occurrence. An analysis of 44 species in the genus Asclepias (milkweeds) revealed similarly large differences between cases using only precise records and those using precise plus imprecise records. Although the inclusion of imprecisely-georeferenced occurrence records may not always be appropriate, careful consideration must be made regarding the trade-offs between fidelity of the match between occurrences and the environment, versus sample size and how well precisely-geolocated records sample species environmental and geographic distributions. The preponderance of imprecisely georeferences specimens in museums and herbaria could be gainfully employed to address known shortfalls in sampling of species distributions and climatic niche tolerances. MTD: a unique pipeline for host and meta-transcriptome joint and integrative analyses of RNA-seq data RNA-seq data contains not only host transcriptomes but also non-host information that comprises transcripts from active microbiota in the host cells. Therefore, joint and integrative analyses of both host and meta-transcriptome can reveal gene expression of microbial community in a given sample as well as the correlative and interactive dynamics of host response to the microbiome. However, there are no convenient tools that can systemically analyze host-microbiota interactions by simultaneously quantifying host and meta-transcriptome in the same sample at the tissue and the single-cell level, which poses a challenge for interested researchers with a limited expertise in bioinformatics. Here, we developed a software pipeline that can comprehensively and synergistically analyze and correlate the host and meta-transcriptome in a single sample using bulk and single-cell RNA-seq data. This pipeline, named MTD, can extensively identify and quantify microbiome, including viruses, bacteria, protozoa, fungi, plasmids, and vectors in the host cells and correlate the microbiome with the host transcriptome. MTD is easy to install and run, involving only a few lines of simple commands. It empowers researchers with unique genomics insights into host immune responses to microorganisms. Host phenology can drive the evolution of intermediate virulence strategies in some parasites AO_SCPLOWBSTRACTC_SCPLOWThe traditional mechanistic trade-offs resulting in a negative correlation between transmission and virulence are the foundation of nearly all current theory on the evolution of parasite virulence. Several ecological factors have been shown to modulate the optimal virulence strategies predicted from mechanistic trade-off models, but these ecological factors have not yet been shown to be sufficient to explain the intermediate virulence strategies observed in any natural system. The timing of seasonal activity, or phenology, is a common factor that influences the types and impact of many ecological interactions but is difficult to incorporate into virulence evolution studies. We develop a mathematical model of a disease system with seasonal host activity to study the evolutionary consequences of host phenology on the virulence of obligate-killer parasite. Results from this model demonstrated that seasonal host activity is sufficient to drive the evolution of intermediate parasite virulence in some types of natural disease systems, even when a traditional mechanistic trade-off between transmission and virulence is not assumed in the modeling framework. The optimal virulence strategy in these systems can be determined by both the duration of the host activity period as well as the variation in the host emergence timing. Parasites with low virulence strategies are favored in environments with long host activity periods and in environments in which hosts emerge synchronously. The results demonstrate that host phenology can be sufficient to select for intermediate optimal virulence strategies, providing an alternative mechanism to account for virulence evolution in some natural systems. Cell trajectory modeling identifies a primitive trophoblast state defined by BCAM enrichment In early placental development, progenitor cytotrophoblasts (CTBs) differentiate along one of two cellular trajectories: the villous or extravillous pathways. CTBs committed to the villous pathway fuse with neighboring CTBs to form the outer multinucleated syncytiotrophoblast (SCT), while CTBs committed to the extravillous pathway differentiate into invasive extravillous trophoblasts (EVT). Unfortunately, little is known about the processes controlling human CTB progenitor maintenance and differentiation. To address this, we established a single cell RNA sequencing (scRNA-seq) dataset from first trimester placentas to identify cell states important in trophoblast progenitor establishment, renewal, and differentiation. Multiple distinct trophoblast states were identified, representing progenitor CTBs, column CTBs, SCT precursors, and EVT. Lineage trajectory analysis identified a progenitor origin that was reproduced in human trophoblast stem cell organoids. Heightened expression of basal cell adhesion molecule (BCAM) defined this primitive state, where BCAM enrichment or gene silencing resulted in enhanced or diminished organoid growth. Together, this work describes at high-resolution trophoblast heterogeneity within the first trimester, resolves gene networks within human CTB progenitors, and identifies BCAM as a primitive progenitor marker and possible regulator.

Summary StatementLineage trajectory modeling identifies multiple human progenitor trophoblast states and defines trophoblast differentiation kinetics, where BCAM-expressing progenitors demonstrate enhanced regenerative ability. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space.