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.
oCEM: Automatic detection and analysis of overlapping co-expressed gene modules BackgroundWhen it comes to the co-expressed gene module detection, its typical challenges consist of overlap between identified modules and local co-expression in a subset of biological samples. The nature of module detection is the use of unsupervised clustering approaches and algorithms. Those methods are advanced undoubtedly, but the selection of a certain clustering method for sample- and gene-clustering tasks is separate, in which the latter task is often more complicated.

ResultsThis study presented an R-package, Overlapping CoExpressed gene Module (oCEM), armed with the decomposition methods to solve the challenges above. We also developed a novel auxiliary statistical approach to select the optimal number of principal components using a permutation procedure. We showed that oCEM outperformed state-of-the-art techniques in the ability to detect biologically relevant modules additionally.

ConclusionsoCEM helped non-technical users easily perform complicated statistical analyses and then gain robust results. oCEM and its applications, along with example data, were freely provided at https://github.com/huynguyen250896/oCEM.
Centriole distal-end proteins CP110 and Cep97 influence centriole cartwheel growth at the proximal-end Centrioles are composed of a central cartwheel tethered to nine-fold symmetric microtubule (MT) blades. The centriole cartwheel and MTs are thought to grow from opposite ends of these organelles, so it is unclear how they coordinate their assembly. We previously showed that an oscillation of Polo-like kinase 4 (Plk4) helps to initiate and time the growth of the cartwheel at the proximal end. Here, we show that CP110 and Cep97 form a complex close to the distal-end of the centriole MTs whose levels rise and fall as the new centriole MTs grow, entrained by the core Cdk/Cyclin oscillator that drives the nuclear divisions in these embryos. These CP110/Cep97 dynamics, however, do not appear to time the period of centriole MT growth directly. Instead, we find that changing the levels of CP110/Cep97 alters the Plk4 oscillation and the growth of the cartwheel at the proximal end. These findings reveal an unexpected crosstalk between factors normally concentrated at opposite ends of the growing centrioles, which may help to coordinate centriole growth. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Single amino-acid mutation in a Drosophila melanogaster ribosomal protein: an insight in uL11 transcriptional activity The ribosomal protein uL11 is located at the basis of the ribosome P-stalk and plays a paramount role in translational efficiency. In addition, no mutant for uL11 is available suggesting that this gene is haplo-insufficient as many other Ribosomal Protein Genes (RPGs). We have previously shown that overexpression of Drosophila melanogaster uL11 enhances the transcription of many RPGs and Ribosomal Biogenesis genes (RiBis) suggesting that uL11 might globally regulate the level of translation through its transcriptional activity. Moreover, uL11 trimethylated on lysine 3 (uL11K3me3) interacts with the chromodomain of the Enhancer of Polycomb and Trithorax Corto, and both proteins co- localize with RNA Polymerase II at many sites on polytene chromosomes. These data have led to the hypothesis that the N-terminal end of uL11, and more particularly the trimethylation of lysine 3, supports the extra-ribosomal activity of uL11 in transcription. To address this question, we mutated the lysine 3 codon using a CRISPR/Cas9 strategy and obtained several lysine 3 mutants. We describe here the first mutants of D. melanogaster uL11. Unexpectedly, the uL11K3A mutant, in which the lysine 3 codon is replaced by an alanine, displays a genuine Minute phenotype known to be characteristic of RPG deletions (longer development, low fertility, high lethality, thin and short bristles) whereas the uL11K3Y mutant, in which the lysine 3 codon is replaced by a tyrosine, is unaffected. In agreement, the rate of translation decreases in uL11K3A but not in uL11K3Y. Co-immunoprecipitation experiments show that the interaction between uL11 and the Corto chromodomain is impaired by both mutations. RNA- seq analysis from wing imaginal discs shows enrichment in the GO categories "glutathione metabolism" for up-regulated genes in both uL11K3A and uL11K3Y mutants and "regulation of transcription" for down-regulated genes in uL11K3A only. Analysis of the cis-regulatory sequences of these genes suggests that uL11 might regulate transcription of target genes in concert with the couple of transcription factors Mad/Med that mediate response to the Bone Morphogenetic Protein (BMP) signaling pathway. 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. Structural and Functional Analyses Explain Pea KAI2 Receptor Diversity and Reveal Stereoselective Catalysis During Signal Perception KAI2 are plant /{beta} hydrolase receptors, which perceive smoke-derived butenolide signals (karrikins) and putative endogenous, yet unidentified phytohormones (KAI2-ligands, KLs). The number of functional KAI2 receptors varies among plant species. It has been suggested that KAI2 gene duplication and sub-functionalization plays an adaptative role for diverse environments or ligand diversification by altering the receptor responsiveness to specific KLs. Legumes represent one of the largest families of flowering plants and contain many essential agronomic crops. Prior to legume diversification, KAI2 underwent duplication, resulting in KAI2A and KAI2B. Integrating plant genetics, ligand perception and enzymatic assays, and protein crystallography, we demonstrate that Pisum sativum KAI2A and KAI2B act as receptors and enzymes with divergent ligand stereoselectivity. KAI2B has a stronger affinity than KAI2A towards the KAI2-ligand (-)-GR24 and remarkably hydrolyses a broader range of substrates including the strigolactone-like isomer (+)-GR24. We determine the crystal structures of PsKAI2B in apo and butenolide-bound states. The biochemical and structural analyses as well as recorded mass spectra of KAI2s reveal a transient intermediate on the catalytic serine and a stable adduct on the catalytic histidine, further illuminating the role of KAI2 not only as receptors but also as bona fide enzymes. Our work uncovers the stereoselectivity of ligand perception and catalysis by evolutionarily diverged KAI2 receptors in KAR/KL signaling pathways and proposes adaptive sensitivity to KAR/KL and strigolactone phytohormones by KAI2B. 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. Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Atg11 promotes Kar9-mediated alignment of the mitotic spindle Asymmetric spindle pole body (SPB) inheritance requires a cascade of events that involve kinases, phosphatases and structural scaffold proteins including molecular motors and microtubule-associated proteins. Higher levels of an SPB component Spc72 and the spindle positioning factor Kar9 at the old SPB, which migrates to the daughter cell, ensure asymmetric SPB inheritance. In this study, by analyzing the protein interaction network of autophagy (Atg)-related proteins with those that constitute the chromosome segregation machinery, and growth inhibition kinetics in the presence of a microtubule poison, we identified Atg11 as a potential regulator of chromosome transmission from a cohort of 35 autophagy mutants in Saccharomyces cerevisiae. Cells lacking Atg11 did not show kinetochore defects but displayed a high rate of chromosome loss and delayed anaphase onset. Atg11 positively interacted with Kar9 and Kip2 and negatively with Dyn1 and Kar3 in mediating proper positioning and alignment of the mitotic spindle that was found to be under the surveillance of the spindle positioning checkpoint. Indeed, atg11{Delta} cells displayed an inverted SPB inheritance. Our experiments further revealed that via the Kar9 pathway, Atg11 promotes asymmetric localization of Spc72 and Kar9 on the old SPB. Atg11 physically interacted with Spc72 and transiently localized close to the old SPB during metaphase-to-anaphase progression. Thus, we uncover an autophagy-independent role of Atg11 in spindle alignment. This study emphasizes the importance of screens to identify factors mediating the complex and intricate crosstalk between biological processes fundamental to genomic integrity. Systematic monitoring of 2-Cys peroxiredoxin-derived redox signals unveiled its role in attenuating carbon assimilation rate Transmission of reductive cues from the photosynthetic electron transport chain to redox-regulated proteins plays a crucial role in activating chloroplast metabolism. However, deciphering the role of their counterbalanced oxidative signals is challenging due to monitoring difficulties. Here, we demonstrate the light-depended redox modification of chloroplast-targeted 2-Cys peroxiredoxins and introduce peroxiredoxin-based biosensors to monitor photosynthetically-derived oxidative signals. By employing a set of genetically encoded biosensors, we show the induction of oxidative signals under habitual light intensities and their inverse relationship with NADPH levels, unraveling the combined activity of reducing and oxidizing signals in fine-tuning chloroplast metabolism. A faster increase in carbon assimilation rates during photosynthesis induction phase was measured in plants deficient in 2-Cys peroxiredoxins compared to wild-type, suggesting the involvement of oxidative signals in attenuating photosynthesis under variable light environments. We suggest that oxidative signals measured by peroxiredoxin-based biosensors reflect the limitation to photosynthesis imposed by the redox regulatory system.

One-Sentence SummaryA genetically encoded biosensor unmasked the dominant role of photosynthetically-derived oxidative signals under habitual conditions. Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes BackgroundThe intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

ResultsThe mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with ten missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmFC gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

ConclusionsThe lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences. Decreased cell stiffness facilitates detachment and migration of breast cancer cells in 3D collagen matrices: An exploratory study BackgroundTumour-cell detachment is a critical early event in the metastatic cascade. However, the role of the cells mechanical properties in detachment and migration is not well understood. This exploratory study aimed to assess how intracellular stiffness changes these processes.

MethodsMDA-MB-231 cells were embedded as 10,000-cell spheroids in 2 and 4 mg/ml collagen matrices. Intracellular stiffness was assessed with mitochondria tracking microrheology of cells that migrated distances equivalent to four and six times the cell diameter (dC) from the spheroid and compared to cells at the spheroid surface (0dC), representing medium, high and no migration, respectively.

Results and discussionThe mitochondrial mean square displacement and intracellular stiffness decreased during detachment and migration for both collagen concentrations (i.e. rigidities). The mean square displacement of 4dC and 6dC cells was similar, whereas cell stiffness was lower for 4dC than for 6dC cells. With increasing matrix rigidity, the intracellular stiffness decreased for 0dC cells and did not change for 4dC and 6dC cells. It is proposed that decreased cell stiffness drives detachment and migration and increased matrix rigidity physically hinders migration, and cells need to become softer or remodel the environment to migrate. The independence of the stiffness of migrated cells from matrix rigidity suggests that cells remodel their environment through matrix proteins cleavage to migrate.

ConclusionsThe study revealed the collective effects of enhanced migratory conditions and increased matrix rigidity on the mechanical properties of the cells. The expression of matrix metalloproteinases and transforming growth factor {beta} and the role of cell volume on detachment and migration in matrices with varying pore sizes are proposed targets for further studies on metastatic cancer cells. Loss of glutamate transporter eaat2a leads to aberrant neuronal excitability, recurrent epileptic seizures and hypoactivity Astroglial excitatory amino acid transporter 2 (EAAT2, GLT-1, SLC1A2) regulates the duration and extent of neuronal excitation by removing glutamate from the synaptic cleft. Hence, an impairment in EAAT2 function could lead to an imbalanced neural network excitability. Here, we investigated the functional alterations of neuronal and astroglial networks associated with the loss of function in the astroglia predominant eaat2a gene in zebrafish. We observed that eaat2a-/- mutant zebrafish larvae display recurrent spontaneous and light-induced seizures in neurons and astroglia, which coincide with an abrupt increase in extracellular glutamate levels. In stark contrast to this hyperexcitability, basal neuronal and astroglial activity was surprisingly reduced in eaat2a-/- mutant animals, which manifested in decreased overall locomotion. Our results reveal an essential and mechanistic contribution of EAAT2a in balancing brain excitability, and its direct link to epileptic seizures. Two levels of selection of rhythmicity in gene expression: energy saving for rhythmic proteins and noise optimization for rhythmic transcripts Many genes have nycthemeral rhythms of expression, i.e. a 24-hours periodic variation, at either mRNA or protein level or both, and, most rhythmic genes are tissue-specific. Here, we investigate and discuss the evolutionary origins of rhythms in gene expression. Our results suggest that thythmicity of protein expression has been favored by selection for low cost. Trends are consistent in bacteria, plants and animals, and are also supported in tissue-specific patterns in mouse. Cost cannot explain rhythm at the RNA level, and we suggest instead it allows to periodically and drastically reduce expression noise. Noise control had strongest support in mouse, with limited power in other species. Genes under stronger purifying selection are rhythmically expressed at the mRNA level, probably because they are noise sensitive genes. We also suggest that mRNA rhythmicity allows to switch between optimal precision and higher stochasticity. Higher precision allows to maximize the robustness of gene expression when the function is most needed, while higher stochasticity allows to maintain oscillations and to exhibit diverse molecular phenotypes, i.e. "blind anticipation" of cells. The ability to alternate between these two states, enabled by rhythmicity at the mRNA level, might be adaptive in fluctuating environments. Finally, the adaptive role of rhythmic expression is also supported by rhythmic genes being highly expressed yet tissue-specific genes. This provides a good evolutionary explanation for the observation that nycthemeral rhythms are often tissue-specific. An Essential Role of UBXN3B in B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development. 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. Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2 Back and forth transmission of SARS-CoV-2 between humans and animals may lead to wild reservoirs of virus that can endanger efforts toward long-term control of COVID-19 in people, and protecting vulnerable animal populations that are particularly susceptible to lethal disease. Predicting high risk host species is key to targeting field surveillance and lab experiments that validate host zoonotic potential. A major bottleneck to predicting animal hosts is the small number of species with available molecular information about the structure of ACE2, a key cellular receptor required for viral cell entry. We overcome this bottleneck by combining species ecological and biological traits with 3D modeling of virus and host cell protein interactions using machine learning methods. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for over 5,000 mammals -- an order of magnitude more species than previously possible. The high accuracy predictions achieved by this approach are strongly corroborated by in vivo empirical studies. We identify numerous common mammal species whose predicted zoonotic capacity and close proximity to humans may further enhance the risk of spillover and spillback transmission of SARS-CoV-2. Our results reveal high priority areas of geographic overlap between global COVID-19 hotspots and potential new mammal hosts of SARS-CoV-2. With molecular sequence data available for only a small fraction of potential host species, predictive modeling integrating data across multiple biological scales offers a conceptual advance that may expand our predictive capacity for zoonotic viruses with similarly unknown and potentially broad host ranges. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. 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. IRI-CCE: A Platform for Evaluating CRISPR-based Base Editing Tools and Its Components Rapid assessment of CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated Cas protein)-based genome editing (GE) tools and their components are critical aspects for successful applications in different organisms. In many bacteria, double-stranded breaks (DSBs) generated by CRISPR/Cas tool generally cause cell death due to the lack of an efficient non-homologous end-joining pathway and restricts its use. CRISPR-based DSB-free base editors (BEs) have been applied for precise nucleotide editing in bacteria, which does not need to make DSBs. However, optimization of newer BE tools in bacteria is challenging owing to the toxic effects of BE reagents expressed using strong promoters. Improved variants of two main BEs capable of converting C-to-T (CBE) and A-to-G (ABE) have been recently developed but yet to be tested in bacteria. Here, we report a platform for in vivo rapid investigation of CRISPR-BE components in Escherichia coli (IRI-CCE) comprising different combinations of promoters/terminators. We demonstrate the use of IRI-CEE to characterize different variants of CBEs (PmCDA1, evoCDA1, APOBEC3A) and ABEs (ABE8e, ABE9e), exhibiting that each independent BE has its specific editing pattern for a given target site and promoter type. Additionally, the IRI-CCE platform offers a rapid way to screen functional gRNAs. In summary, CRISPR-BE components expressed by promoters of different strengths in the IRI-CCE allow an analysis of various BE tools, including cloned biopart modules and gRNAs. Identification of the three zinc-binding sites on Tau protein Tau protein has been extensively studied due to its key roles in microtubular cytoskeleton regulation and in the formation of aggregates found in some neurodegenerative diseases. Recently it has been shown that zinc is able to induce tau aggregation by interacting with several binding sites. However, the precise location of these sites and the molecular mechanism of zinc-induced aggregation remain unknown. Here we used Nuclear Magnetic Resonance (NMR) to identify zinc binding sites on hTau40 isoform. These experiments revealed three distinct zinc binding sites on tau, located in the N-terminal part (H14, H32, H94, and H121), the repeat region (H299, C322, H329 and H330) and the C-terminal part (H362, H374, H388 and H407). Further analysis enabled us to show that the C-terminal and the N-terminal sites are independent of each other. Using molecular simulations, we modeled the structure of each site in a complex with zinc. Given the clinical importance of zinc in tau aggregation, our findings pave the way for designing potential therapies for tauopathies.

HighlightsO_LIZinc is known to induce tau aggregation in neurodegenerative diseases
C_LIO_LIZinc binding locations and mechanism are not yet clear
C_LIO_LIUsing NMR we localized 3 zinc binding site on tau
C_LIO_LIBy molecular simulations, we proposed a modeled structure of each site
C_LIO_LIOur findings pave the way for designing potential therapies for tauopathies
C_LI Cargo surface fluidity reduces inter-motor mechanical interference, promotes load-sharing and enhances run-lengths in an ATP dependent manner In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of cytoskeletal filaments. How motors that typically mechanically interfere with each other, work as teams is unclear. Here we explored the possibility that purely physical mechanisms may potentially enhance teamwork, both at the single motor and cargo level. To explore these mechanisms, we developed a 3D dynamical simulation of cargo transport along microtubules by teams of canonically non-cooperative kinesin-1 motors. We accounted for cargo membrane fluidity by explicitly simulating the Brownian dynamics of motors on the cargo surface and considered both the load and ATP dependence of single motor functioning. We showed explicitly, for the first time, that surface fluidity leads to the reduction of negative mechanical interference between kinesins, enhancing load sharing thereby decreasing single motor off-rates and increasing processivity. Remarkably, we also showed that, independent of fluidity, increasing numbers of bound motors pull the cargo closer to the microtubule, increasing the on-rates of individual unbound motors, resulting in a cooperative increase in bound motor numbers that depends on 3D cargo geometry. At the cargo level, surface fluidity makes more motors available for binding, though this effect is significant only at low ATP or high motor density. Interestingly, we find that the fluidity induced reduction in mechanical interference dominates over the increased availability of motors for typical physiological and in vitro conditions, this allowing us to reconcile different experimental results in different regimes. Finally, we show that these effects can altogether result in enhanced mechanical efficiency and tunable cargo run-lengths for teams of molecular motors over physiological ranges of fluidity with implications for new experimental validation efforts, transport in vivo, artificial cargo design and multi-motor-driven mechanical processes in general.

Author summaryIn cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of protein filaments. How do intrinsically non-cooperative motors function as teams within cells? In this paper, we show, using computer simulations, that the fluid surfaces of cellular cargo reduce the mechanical interference between motors allowing better load sharing and decreasing their detachment rates, thereby increasing the distance over which they can carry cargo. We additionally found that motors pull cargo closer to the filament increasing the attachment rates of other unbound motors. These effects act synergistically with an increased availability of motors due to fluidity to further increase travel distances. Our simulation model accounted for the dynamics of the motors on the cargo surface and their load and ATP dependent kinetics, allowing us to connect single motor properties to overall transport. Our work on understanding how teamwork arises in mechanically coupled motors sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving mechanical efficiency and transport of artificial cargo powered by motor teams. Game of theories: an agent-based model of how variation in contest cost-benefits determines the evolution of contest resolution rules Game-theory based models are used to understand rules that animals use to settle contests over indivisible resources. However, the empirical literature of contests indicates controversial support to models, with some species supporting different models and other species showing no support to any model. Since strategies used to resolve contests may have different associated costs, it is possible that different conditions have determined the evolution of distinct assessment strategies used by animals. We used an agent-based model to explore the importance of the following conditions: resource availability, probability of reproduction with resource, and damage costs on evolution of assessment strategies. We used self- and mutual-assessment models as a heurist framework to build agents with different assessment strategies. In our model, agents competed for resources in scenarios with different combinations of resource availability, probability of reproduction with resource, and damage costs. We found that agents following the self-assessment with damage strategy were prevalent in scenarios with no probability of reproduction without the resource, independently of other variables. We also found that agents following the non-aggressive strategy occurred in all scenarios. However, agents using the non-aggressive strategy were prevalent only in scenarios with probability of reproduction with the resource. Finally, we observed that agents using mutual-assessment occurred only in a scenario with high risk of damage, low availability of resources, and with probability of reproduction without the resource. These results indicate that agents following the self-assessment with damage and non-aggressive strategies may be able to stay at most scenarios.