Fixed or random? On the reliability of mixed-effects models for a small number of levels in grouping variables O_LIBiological data are often intrinsically hierarchical. Due to their ability to account for such dependencies, mixed-effects models have become a common analysis technique in ecology and evolution. While many questions around their theoretical foundations and practical applications are solved, one fundamental question is still highly debated: When facing a low number of levels should we model a grouping (blocking, clustering) variable as a random or fixed effect? In such situation, the variance of the random effect is imprecise, but whether this affects the statistical properties of the population effect is unclear.
C_LIO_LIHere, we analyzed the consequences of including a grouping variable as fixed or random effect in the correctly specified and other possible miss-specified models (too simple or too complex models) for data with small number of levels (2 - 8). For all these options, we calculated type I error rates and power. Moreover, we show how these statistical properties change with the study design.
C_LIO_LIWe found that the model choice does not influence the statistical properties of the population effect when the effect is the same at all levels However, if an ecological effect differs among levels, using a random slope and intercept model, and switching to a fixed-effect model only in case of a singular fit, avoids overconfidence in the results. Additionally, power and type I error are strongly influenced by the number of and difference between levels.
C_LIO_LIWe conclude that inferring the correct random effect structure is of high importance to get correct statistical properties. When in doubt, we recommend starting with the simpler model and using model diagnostics to identify missing components. When having identified the correct structure, we encourage to start with a mixed-effects model independent of the number of levels and switch to a fixed-effect model only in case of a singular fit. With these recommendations, we allow for more informative choices about study design and data analysis and thus make ecological inference with mixed-effects models more robust for small number of levels.
C_LI
Rapid global phaseout of animal agriculture has the potential to stabilize greenhouse gas levels for 30 years and offset 68 percent of CO2 emissions this century Animal agriculture contributes significantly to global warming through ongoing emissions of the potent greenhouse gases methane and nitrous oxide, and displacement of biomass carbon on the land used to support livestock. However, because estimates of the magnitude of the effect of ending animal agriculture often focus on only one factor, the full potential benefit of a more radical change remains underappreciated. Here we quantify the full "climate opportunity cost" of current global livestock production, by modeling the combined, long-term effects of emission reductions and biomass recovery that would be unlocked by a phaseout of animal agriculture. We show that, even in the absence of any other emission reductions, persistent drops in atmospheric methane and nitrous oxide levels, and slower carbon dioxide accumulation, following a phaseout of livestock production would, through the end of the century, have the same cumulative effect on the warming potential of the atmosphere as a 25 gigaton per year reduction in anthropogenic CO2 emissions, providing half of the net emission reductions necessary to limit warming to 2{degrees}C. The magnitude and rapidity of these potential effects should place the reduction or elimination of animal agriculture at the forefront of strategies for averting disastrous climate change.

Significance StatementThe use of animals to produce food has a negative impact on the climate, but the benefits of a global switch to a plant based diet are underappreciated. We show that the global warming impact, through the rest of this century, of eliminating greenhouse gas emissions from livestock and allowing native ecosystems to regrow on the land currently used to house and feed livestock, would be equivalent to a 68% reduction in carbon dioxide emissions. We hope putting clearer numbers on the "climate opportunity cost" of our continued use of animals as food technology will help policymakers and the public properly prioritize dietary change as a climate defense strategy.

Declaration of Conflict of InterestPatrick Brown is the founder and CEO of Impossible Foods, a company developing alternatives to animals in food-production. Michael Eisen is an advisor to Impossible Foods. Both are shareholders in the company and thus stand to benefit financially from reduction of animal agriculture.
A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. The coiled coil of FtsLB is a structurally detuned element critical for modulating activation in bacterial cell division The FtsLB complex is a key regulator of bacterial cell division, existing in either an off or on state which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a "detuned" structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the "Y-model") in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model ("I-model"). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Liver Endothelium Promotes HER3-mediated Cell Survival in KRAS Wild-type and Mutant Colorectal Cancer Cells We previously showed that liver endothelial cells (ECs) secreted soluble factors in a paracrine fashion and activated human epidermal growth factor receptor 3 (HER3, also known as ERBB3) for promoting colorectal cancer (CRC) growth and chemoresistance. However, RAS proteins play a critical role in receptor tyrosine kinase signaling pathways, and KRAS mutations mediate CRC resistance to therapies targeting EGFR, another HER protein. Therefore, the role of KRAS mutation status in EC-induced HER3 activation and CRC survival was investigated as it has therapeutic implications. We used CRC cell lines and patient-derived xenografts harboring KRAS wild-type or mutant genes and demonstrated that liver EC-secreted factors promoted HER3-mediated CRC cell growth independent of KRAS mutation status. Also, blocking HER3 in CRC cells by siRNAs or a HER3 antibody seribantumab blocked EC-induced CRC survival. Our findings highlight the potential of utilizing HER3 targeted therapies for treating patients with mCRC independent of RAS mutational status. Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. Nontypeable Haemophilus influenzae infection impedes Pseudomonas aeruginosa colonization and persistence in mouse respiratory tract Patients with cystic fibrosis (CF) experience lifelong respiratory infections which are a significant cause of morbidity and mortality. These infections are polymicrobial in nature, and the predominant bacterial species undergo a predictable series of changes as patients age. Young patients have populations dominated by opportunists that are typically found within the microbiome of the human nasopharynx, such as nontypeable Haemophilus influenzae (NTHi); these are eventually supplanted and the population within the CF lung is later dominated by pathogens such as Pseudomonas aeruginosa (Pa). In this study, we investigated how initial colonization with NTHi impacts colonization and persistence of Pa in the respiratory tract.

Analysis of polymicrobial biofilms in vitro by confocal microscopy revealed that NTHi promoted greater levels of Pa biofilm volume and diffusion. However, sequential respiratory infection of mice with NTHi followed by Pa resulted in significantly lower Pa as compared to infection with Pa alone. Coinfected mice also had reduced airway tissue damage and lower levels of inflammatory cytokines as compared with Pa infected mice. Similar results were observed after instillation of heat-inactivated NTHi bacteria or purified NTHi lipooligosaccharide (LOS) endotoxin prior to Pa introduction. Based on these results, we conclude that NTHi significantly reduces susceptibility to subsequent Pa infection, most likely due to priming of host innate immunity rather than a direct competitive interaction between species. These findings have potential significance with regard to therapeutic management of early life infections in patients with CF. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. A pilot study on the home range and movement patterns of the Andean Fox Lycalopex culpaeus (Molina 1782) in Cotopaxi National Park, Ecuador This study reports movement patterns and home range estimates of an Andean fox (Lycalopex culpaeus) in Cotopaxi National Park in Ecuador, representing the first GPS-tagging of the species. The GPS functioned well during the 197-day tracking period. Home range sizes ranged between 4.9 - 8.1 km2, depending on the estimation method. Movement speeds averaged 0.17 km/hr at day vs. 0.87 km/hr. at night, and distance traveled averaged 0.23 km at day vs. 0.89 km at night. These preliminary results highlight the importance of collecting unbiased, high-quality data which enables an enhanced understanding on mammal behavior and human/animal interaction. Loss of SUMO-specific protease 2 causes isolated glucocorticoid deficiency by blocking adrenal cortex zonal transdifferentiation SUMOylation is a dynamic posttranslational modification, that provides fine-tuning of protein function involved in the cellular response to stress, differentiation, and tissue development. In the adrenal cortex, an emblematic endocrine organ that mediates adaptation to physiological demands, the SUMOylation gradient is inversely correlated with the gradient of cellular differentiation raising important questions about its role in functional zonation and the response to stress. Considering that SUMO-specific protease 2 (SENP2), a deSUMOylating enzyme, is upregulated by ACTH/PKA signalling within the zona Fasciculata (zF), we generated mice with adrenal-specific Senp2 loss to address these questions. Disruption of SENP2 activity in steroidogenic cells leads to specific hypoplasia of the zF, a blunted reponses to ACTH and isolated glucocorticoid deficiency. Mechanistically, overSUMOylation resulting from SENP2 loss shifts the balance between ACTH/PKA and WNT/{beta}-catenin signalling leading to repression of PKA activity and ectopic activation of {beta}-catenin. At the cellular level, this blocks transdifferentiation of {beta}-catenin-positive zona Glomerulosa cells into zF cells and sensitises them to premature apoptosis. Our findings indicate that the SUMO pathway is critical for adrenal homeostasis and stress responsiveness. 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. Essential Regression - a generalizable framework for inferring causal latent factors from multi-omic human datasets High-dimensional cellular and molecular profiling of human samples highlights the need for analytical approaches that can integrate multi-omic datasets to generate predictive biomarkers and prioritized causal inferences. Current methods are limited by high dimensionality of the combined datasets, the differences in their data distributions and their integration to infer causal relationships. Here we present Essential Regression (ER), an interpretable machine learning approach for high-dimensional multi-omic datasets, that addresses these problems by identifying latent factors and their likely cause-effect relationships with the system-wide outcome/properties of interest. ER is a novel data-distribution-free latent-factor regression model that integrates multi-omic datasets and identifies latent factors significantly associated with an outcome. ER outperforms a range of state-of-the-art methods in terms of prediction performance on simulated datasets. ER can be coupled with probabilistic graphical modeling thereby strengthening the causal inferences. ER generates novel cellular and molecular predictions, using multi-omic human systems immunology datasets, pertaining to immunosenescence and immune dysregulation. Exosomal TAR DNA binding protein 43 profile in canine model of amyotrophic lateral sclerosis: A preliminary study in developing blood-based biomarker for neurodegenerative diseases ObjectiveBlood-based biomarkers provide a crucial information in progress of neurodegenerative diseases with minimally invasive sampling method. Validated blood-based biomarker application in people with amyotrophic lateral sclerosis would derive numerous benefits. Canine degenerative myelopathy is a naturally occurring animal disease model to study the biology of human amyotrophic lateral sclerosis. Serum derived exosomes are potential carriers for cell-specific cargoes making them ideal venue to study biomarkers for a variety of diseases and biological processes. This study assessed the exosomal proteins that may be assigned as surrogate biomarker that may reflect biochemical changes in central nervous system.

MethodsExosomes were isolated from canine serum using commercial exosome isolation reagents. Exosomes target proteins contents were analysed by Western blotting method.

ResultsThe profiles of potential biomarker candidates in spinal cord homogenate and that of serum-derived exosomes were found elevated in dogs with degenerative myelopathy as compare to control subjects.

ConclusionsSerum-derived exosomal biomolecules can serve as surrogate biomarkers in neuro degenerative diseases.

Key MessagesO_LIA canine with degenerative myelopathy can serve as a model animal to study human amyotrophic lateral sclerosis.
C_LIO_LISerum-derived exosomes contains Transactive Response DNA Binding Protein 43 (TDP-43), potential biomarker candidate.
C_LIO_LIThe levels of spinal cord TDP-43 proteins and that of serum-derived exosomes exhibited a similar profiling. Therefore, serum derived exosomes may be used as a venue for establishing blood-based biomarkers for neurodegenerative diseases.
C_LI Kidney failure alters parathyroid Pin1 phosphorylation and parathyroid hormone mRNA binding proteins leading to secondary hyperparathyroidism Parathyroid hormone (PTH) regulates calcium metabolism and bone strength. Chronic kidney disease leads to secondary hyperparathyroidism (SHP) which increases morbidity and mortality. In experimental SHP, the increased PTH gene expression is due to enhanced PTH mRNA stability, mediated by changes in its interaction with stabilizing AUF1 and destabilizing KSRP. Pin1 isomerizes target proteins and leads to KSRP dephosphorylation. In SHP, Pin1 isomerase activity is decreased and phosphorylated KSRP fails to bind PTH mRNA, resulting in high PTH mRNA stability and levels. However, the up- and down-stream mechanisms by which kidney failure stimulates the parathyroid remain elusive. We now reveal a pathway where kidney failure induces parathyroid Pin1 phosphorylation, linking Pin1, KSRP and PTH mRNA stability as key players for the high PTH expression in SHP. We identified by mass-spectrometry, changes in rat parathyroid proteome and phosphoproteome profiles induced by impaired renal function, including KSRP phosphorylation at Pin1 target sites. Furthermore, both acute and chronic kidney failure led to parathyroid-specific Pin1 Ser16 and Ser71 phosphorylation, which disrupts Pin1 activity. Accordingly, pharmacologic Pin1 inhibition, that mimics the decreased Pin1 activity in SHP, increased PTH expression ex-vivo in parathyroid organ cultures and in transfected cells, through the PTH mRNA protein-interacting element and KSRP phosphorylation at potential Pin1-binding motifs. Therefore, kidney failure leads to loss of parathyroid Pin1 activity by inducing Pin1 phosphorylation. This predisposes parathyroids to increase PTH production through impaired PTH mRNA decay that is dependent on KSRP phosphorylation. Pin1 and KSRP phosphorylation and the Pin1-KSRP-PTH mRNA axis thus drive SHP. A memetic algorithm enables efficient local and global all-atom protein-protein docking with backbone and sidechain flexibility Protein complex formation is encoded by specific interactions at the atomic scale, but the computational cost of modeling proteins at this level often requires the use of simplified energy models and limited conformational flexibility. In particular, the use of all-atom energy functions, backbone and sidechain flexibility results in rugged energy landscapes that are difficult to explore. In this study we develop a protein-protein docking algorithm, EvoDOCK, that combine the strength of a differential evolution algorithm for efficient exploration of the global search space with the benefits of a local optimization method to refine detailed atomic interactions. EvoDOCK enabled accurate and fast local and global protein-protein docking using an all-atom energy function with side-chain flexibility. Comparison with a standard method built on Monte Carlo optimization demonstrated improved accuracy and with increases in computational speed of up to 35 times. The evolutionary algorithm also enabled efficient atomistic docking with backbone flexibility. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. 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. Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. Neutralization potential of Covishield vaccinated individuals sera against B.1.617.1 Covishield comprises the larger proportion in the vaccination program in India. Hence, it is of utmost importance to understand neutralizing capability of vaccine against the B.1.617.1 variant which is considered to responsible for surge of the cases in India. The neutralizing-antibody (NAb) titer against B.1.167.1 and prototype B.1 variant (D614G) was determined of the vaccine sera (4 weeks after second dose) of COVID-19 naive subjects (n=43) and COVID-19 recovered subjects (n=18). The results demonstrated that sera of COVID-19 recovered subjects (n=18) who received two doses of Covishield have higher NAb response compared to the COVID-19 naive with a significant difference (p<0.0001) in NAb titer against B.1 and B.1.617.1 In-spite of reduction in the neutralizing titer against B.1.617.1 variant; Covishield vaccine-induced antibodies are likely to be protective to limit the severity and mortality of the disease in the vaccinated individuals. 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. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Strong immobilization of C. elegans en masse on cultivation plates by thermoelectric cooling for high-throughput submicron-resolution microscopy Despite its profound impact on biology, most high-resolution in vivo microscopy approaches remain low throughput because current immobilization techniques require significant manual effort. We greatly accelerate imaging of the nematode Caenorhabditis elegans by implementing a simple cooling approach to easily immobilize entire populations directly on their cultivation plates. We optimize and characterize cooling immobilization. Counterintuitively, relatively warmer temperatures immobilize animals significantly more effectively than colder temperatures utilized in prior studies. This enhanced immobilization enables clear submicron-resolution fluorescence imaging, which is challenging to achieve with most current immobilization techniques. We demonstrate 64x magnification 3D imaging and timelapse recording of neurons in adults and embryos without motion blur. Compared to standard azide immobilization, cooling immobilization reduces the animal preparation and recovery time by >98%, significantly increasing experimental speed. By obviating individual animal manipulation, our approach could also empower automated imaging of large C. elegans populations within standard experimental setups and workflows. A quartz crystal microbalance method to quantify the size of hyaluronan and other glycosaminoglycans on surfaces Hyaluronan (HA) is a major component of peri- and extra-cellular matrices and plays important roles in many biological processes such as cell adhesion, proliferation and migration. The abundance, size distribution and presentation of HA dictate its biological effects and are also useful indicators of pathologies and disease progression. Methods to assess the molecular mass of free-floating HA and other glycosaminoglycans (GAGs) are well established. In many biological and technological settings, however, GAGs are displayed on surfaces, and methods to obtain the size of surface-attached GAGs are lacking. Here, we present a method to size HA that is end-attached to surfaces. The method is based on the quartz crystal microbalance with dissipation monitoring (QCM-D) and exploits that the softness and thickness of films of grafted HA increase with HA size. These two quantities are sensitively reflected by the ratio of the dissipation shift ({Delta}D) and the negative frequency shift (-{Delta}f) measured by QCM-D upon the formation of HA films. Using a series of size-defined HA preparations, ranging in size from [~]2 kDa tetrasaccharides to [~]1 MDa polysaccharides, we establish a monotonic yet non-linear standard curve of the {Delta}D/-{Delta}f ratio as a function of HA size, which reflects the distinct conformations adopted by grafted HA chains depending on their size and surface coverage. We demonstrate that the standard curve can be used to determine the mean size of HA, as well as other GAGs, such as chondroitin sulfate and heparan sulfate, of preparations of previously unknown size in the range from 1 to 500 kDa, with a resolution of better than 10%. For polydisperse samples, our analysis shows that the process of surface-grafting preferentially selects smaller GAG chains, and thus reduces the average size of GAGs that are immobilised on surfaces comparative to the original solution sample. Our results establish a quantitative method to size HA and other GAGs grafted on surfaces, and also highlight the importance of sizing GAGs directly on surfaces. The method should be useful for the development and quality control of GAG-based surface coatings in a wide range of research areas, from molecular interaction analysis to biomaterials coatings. Genetic coupling of life-history and aerobic performance in Atlantic salmon A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions (vgll3 and six6) makes Atlantic salmon (Salmo salar) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes -standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS) - and the life-history genomic regions and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Further, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile-stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve.