A self-supervised domain-general learning framework for human ventral stream representation Anterior regions of the ventral visual stream encode substantial information about object categories. Are top-down category-level forces critical for arriving at this representation, or can this representation be formed purely through domain-general learning of natural image structure? Here we present a fully self-supervised model which learns to represent individual images, rather than categories, such that views of the same image are embedded nearby in a low-dimensional feature space, distinctly from other recently encountered views. We find (i) category information implicitly emerges in the local similarity structure of this feature space, and (ii) these models learn hierarchical features which capture the structure of brain responses across the human ventral visual stream, on par with category-supervised models. These results provide computational support for a domain-general framework guiding the formation of visual representation, where the proximate goal is not explicitly about category information, but is instead to learn unique, compressed descriptions of the visual world.
Disruption of miR-18a alters proliferation, photoreceptor replacement kinetics, inflammatory signaling and microglia/macrophage numbers during retinal regeneration in zebrafish In mammals, photoreceptor loss causes permanent blindness, but in zebrafish (Danio rerio), photoreceptor loss reprograms Muller glia to function as stem cells, producing progenitors that regenerate photoreceptors. MicroRNAs (miRNAs) regulate CNS neurogenesis, but the roles of miRNAs in injury-induced neuronal regeneration are largely unknown. In the embryonic zebrafish retina, miR-18a regulates photoreceptor differentiation. The purpose of the current study was to determine, in zebrafish, the function of miR-18a during injury-induced photoreceptor regeneration. RT-qPCR, in-situ hybridization and immunohistochemistry showed that miR-18a expression increases throughout the retina between 1 and 5 days post-injury (dpi). To test miR-18a function during photoreceptor regeneration, we used homozygous miR-18a mutants (miR-18ami5012), and knocked down miR-18a with morpholino oligonucleotides. During photoreceptor regeneration, miR-18ami5012 retinas have fewer mature photoreceptors than WT at 7 and 10 dpi, but there is no difference at 14 dpi, indicating that photoreceptor regeneration is delayed. Labeling dividing cells with 5-bromo-2'-deoxyuridine (BrdU) showed that at 7 and 10 dpi, there are excess dividing progenitors in both mutants and morphants, indicating that miR-18a negatively regulates injury-induced proliferation. Tracing 5-ethynyl-2'-deoxyuridine (EdU) and BrdU-labeled cells showed that in miR-18ami5012 retinas excess progenitors migrate to other retinal layers in addition to the photoreceptor layer. Inflammation is critical for photoreceptor regeneration, and RT-qPCR showed that in miR-18ami5012 retinas, inflammatory gene expression and microglia activation are prolonged. Suppressing inflammation with dexamethasone rescues the miR-18ami5012 phenotype. Together, these data show that in the injured zebrafish retina, disruption of miR-18a alters proliferation, inflammation, the microglia/macrophage response, and the timing of photoreceptor regeneration.
A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. The role of Kabuki Syndrome genes KMT2D and KDM6A in development: Analysis in Human sequencing data and compared to mice and zebrafish KMT2D and KDM6A are epigenetic regulators that have been implicated in Kabuki Syndrome, a rare congenital birth defect with multiple tissue and organ abnormalities, including craniofacial and heart defects. Our previous study identified human families with mutations in the epigenetic modifiers KMT2D and KDM6A, which is implicated in 32% and 10% of Kabuki Syndrome patients respectively. To understand the connection to Kabuki syndrome patients, and the transcriptional targets of KMT2D and KDM6A in humans, we performed RNA sequencing (seq) of lymphoblastoid cells from Kabuki Syndrome patients carrying mutations in KMT2D and KDM6A. We identified 1995 significant changes in transcriptional targets for KMT2D and 1917 for KDM6A, as compared to control. When compared with RNA-seq datasets obtained from other mouse and zebrafish studies, our analysis revealed KMT2D mutations affect the expression of 76 orthologous genes across all three datasets. Similarliy, KDM6A afftects the expprssion of 7 orthologous genes across three datasets. Despite the differences in cell types, stages, and species in the comparison between the transcriptomic datasets, there are common gene expression changes associated with KMT2D and KDM6A mutations. qPCR on novel zebrafish mutants confirmed the differentially expression in KMT2D or KDM6A mutant backgrounds. Taken together, our results show that KMT2D and KDM6A regulate common and unique genes across humans, mice, and zebrafish for early craniofacial and cardiac development and this information contributes to the understanding of epigenetic dysregulation during development of Kabuki syndrome. A multi-state dynamic process confers mechano-adaptation to the bacterial flagellar motor Adaptation is a defining feature of living systems. The bacterial flagellar motor adapts to changes in the external mechanical environment by adding or removing torque-generating stator units. However, the molecular mechanism for mechanosensitive motor remodeling remains unclear. Here, we induced stator disassembly using electrorotation, followed by the time-dependent assembly of the individual stator units into the motor after electrorotation was terminated. From these experiments, we extracted detailed statistics of the dwell times that comprise the stochastic dynamics of the binding and unbinding of stator units. Dwell times reveal multiple timescales, indicating the existence of multiple binding states of the stator units. Based on these results, we propose a minimal model in which the stator unit can occupy four different states - two bound states with very different rates of unbinding, a diffusive unbound state, and a transiently detached state. Our minimal model quantitatively explains multiple features of the experimental data and allows us to determine the transition rates among all four states. Our experiments and modeling suggest a mechanism of mechano-adaptive remodeling of the bacterial flagellar motor in which torque generated by bound stator units controls their effective unbinding rate by modulating the transition between the two bound states. Furthermore, the binding rate of stator units with the motor has a non-monotonic dependence on the number of bound units, likely because of two counteracting effects of motor rotation on the binding process. Sensing Local Field Potentials with a Directional and Scalable Depth Electrode Array A variety of electrophysiology tools are available to the neurosurgeon for diagnosis, functional therapy, and neural prosthetics. However, no tool can currently address these three critical needs: (i) access to all cortical regions in a minimally invasive manner; (ii) recordings with microscale, mesoscale, and macroscale resolutions simultaneously; and (iii) access to spatially distant multiple brain regions that constitute distributed cognitive networks. We present a novel device for recording local field potentials (LFPs) with the form factor of a stereo-electroencephalographic electrode but combined with radially positioned microelectrodes and using the lead body to shield LFP sources, enabling directional sensitivity and scalability, referred to as the DISC array. As predicted by our electro-quasistatic models, DISC demonstrated significantly improved signal-to-noise ratio, directional sensitivity, and decoding accuracy from rat barrel cortex recordings during whisker stimulation. Critically, DISC demonstrated equivalent fidelity to conventional electrodes at the macroscale and uniquely, revealed stereoscopic information about current source density. Directional sensitivity of LFPs may significantly improve brain-computer interfaces and many diagnostic procedures, including epilepsy foci detection and deep brain targeting. SARS-CoV-2 B.1.617.2 Delta variant emergence, replication and immune evasion The SARS-CoV-2 B.1.617.2 (Delta) variant was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha). In vitro, B.1.617.2 is 6-fold less sensitive to serum neutralising antibodies from recovered individuals, and 8-fold less sensitive to vaccine-elicited antibodies as compared to wild type Wuhan-1 bearing D614G. Serum neutralising titres against B.1.617.2 were lower in ChAdOx-1 versus BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies against the receptor binding domain (RBD) and N-terminal domain (NTD), in particular to the clinically approved bamlavinimab and imdevimab monoclonal antibodies. B.1.617.2 demonstrated higher replication efficiency in both airway organoid and human airway epithelial systems as compared to B.1.1.7, associated with B.1.617.2 spike being in a predominantly cleaved state compared to B.1.1.7. Additionally we observed that B.1.617.2 had higher replication and spike mediated entry as compared to B.1.617.1, potentially explaining B.1.617.2 dominance. In an analysis of over 130 SARS-CoV-2 infected healthcare workers across three centres in India during a period of mixed lineage circulation, we observed substantially reduced ChAdOx-1 vaccine efficacy against B.1.617.2 relative to non-B.1.617.2. Compromised vaccine efficacy against the highly fit and immune evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era. Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. Entorhinal cortex epigenome-wide association study highlights four novel loci showing differential methylation in Alzheimer's disease BackgroundStudies on DNA methylation (DNAm) in Alzheimers disease (AD) have recently highlighted several genomic loci showing association with disease onset and progression.

MethodsHere, we conducted an epigenome-wide association study (EWAS) using DNAm profiles in entorhinal cortex (EC) from 149 AD patients and control brains and combined these with two previously published EC datasets by meta-analysis (total n=337).

ResultsWe identified 12 cytosine-phosphate-guanine (CpG) sites showing epigenome-wide significant association with either case-control status or Braaks tau-staging. Four of these CpGs, located in proximity to CNFN/LIPE, TENT5A, PALD1/PRF1, and DIRAS1, represent novel findings. Integrating DNAm levels with RNA sequencing-based mRNA expression data generated in the same individuals showed significant DNAm-mRNA correlations for 6 of the 12 significant CpGs. Lastly, by calculating rates of epigenetic age acceleration using two recently proposed "epigenetic clock" estimators we found a significant association with accelerated epigenetic aging in AD patients vs. controls.

ConclusionIn summary, our study represents the hitherto most comprehensive EWAS in AD using EC and highlights several novel differentially methylated loci with potential effects on gene expression. Structural basis for cell-type specific evolution of viral fitness by SARS-CoV-2 As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants which is of particular concern due to their potential for increased transmissibility and pathology. In addition to this entrenched variant diversity in circulation, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriS{Delta} variant, originally identified as a viral subpopulation by passaging SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike glycoprotein an eight amino-acid deletion encompassing the furin recognition motif and S1/S2 cleavage site. Here, we elucidate the structure, function and molecular dynamics of this variant spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Moreover, our study reveals long-range allosteric communication between functional regions within the spike that differ in wild-type and deletion variant. Our results support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host. A Wnt/β-catenin signaling driven adipocyte population required for beiging and glucose homeostasis Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Intron-mediated induction of phenotypic heterogeneity Introns are universally present in the nuclear genomes of eukaryotes1. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein (RP) genes differing primarily in their introns2-4. Despite recent findings on the role of RP introns under stress and starvation5-7, understanding the contribution of introns to ribosome regulation remains challenging. Here, combining isogrowth profiling8 with single-cell protein measurements9, we found that introns can mediate inducible phenotypic heterogeneity conferring a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, mediated by an intron in the 5 untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low Rps22B protein levels resulted in prolonged survival under sustained starvation, while high Rps22B levels enabled cells to grow faster after transient starvation. Further, yeast growing at high sugar concentrations - similar to those in ripe grapes - exhibit bimodal Rps22B expression when approaching stationary phase. Differential intron-mediated regulation of RP genes thus provides a way to diversify the population when starvation looms in natural environments. Our findings reveal a new role for introns in inducing phenotypic heterogeneity in changing environments and suggest that duplicated RP genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness10. Shedding Light on the Circadian Clock of the Threespine Stickleback The circadian clock is an internal timekeeping system shared by most organisms, and knowledge about its functional importance and evolution in natural environments is still needed. Here, we investigated the circadian clock of wild-caught threespine sticklebacks (Gasterosteus aculeatus) at the behavioural and molecular levels. While their behaviour, ecology, and evolution are well studied, information on their circadian rhythms are scarce. We quantified the daily locomotor activity rhythm under a light-dark cycle (LD) and under constant darkness (DD). Under LD, all fish exhibited significant daily rhythmicity, while under DD, only 18% of individuals remained rhythmic. This interindividual variation suggests that the circadian clock controls activity only in certain individuals. Moreover, under LD, some fish were almost exclusively nocturnal, while others were active around the clock. Furthermore, the most nocturnal fish were also the least active. These results suggest that light masks activity (i.e. suppresses activity without entraining the internal clock) more strongly in some individuals than others. Finally, we quantified the expression of five clock genes in the brain of sticklebacks under DD using qPCR. We did not detect circadian rhythmicity, which could either indicate that the clock molecular oscillator is highly light-dependent, or that there was an oscillation but that we were unable to detect it. Overall, our study suggests that a strong circadian control on behavioural rhythms may not necessarily be advantageous in a natural population of sticklebacks and that the daily phase of activity varies greatly between individuals because of a differential masking effect of light.

Summary statementWe found that in wild-caught threespine sticklebacks, the circadian clock does not control locomotor activity in most, but not all, individuals. Sticklebacks are mostly nocturnal, although interindividual variation exists. 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. Modelling interference between vectors of non-persistently transmitted plant viruses to identify effective control strategies Aphids are the primary vector of plant viruses. Transient aphids, which probe several plants per day, are considered to be the principal vectors of non-persistently transmitted (NPT) viruses. However, resident aphids, which can complete their life cycle on a single host and are affected by agronomic practices, can transmit NPT viruses as well. Moreover, they can interfere both directly and indirectly with transient aphids, eventually shaping plant disease dynamics. By mean of an epidemiological model, originally accounting for ecological principles and agronomic practices, we explore the consequences of fertilization and irrigation, pesticide deployment and roguing of infected plants on the spread of viral diseases in crops. Our results indicate that the spread of NPT viruses can be i) both reduced or increased by fertilization and irrigation, depending on whether the interference is direct or indirect; ii) counter-intuitively increased by pesticide application and iii) reduced by roguing infected plants. We show that a better understanding of vectors interactions would enhance our understanding of disease transmission, supporting the development of disease management strategies.

Author summaryA range of both experimental and theoretical studies show that the behaviour and population dynamics of insects depend strongly upon interactions with other insect species. These interactions have the potential to greatly affect the dynamics of insect-vectored plant disease, as transmission of viruses is intimately dependent on the local density of vectors, as well as how they select and move between potential host plants. Surprisingly, the effects of interaction between vector species on epidemics remains little studied and even worse understood, probably because experimentation is costly and difficult. Here, we present a model which permits to investigate the effect of interaction between a virus, two vector species, and the host plant on the spread of viral disease in crops. In this study, our model is used to explore the consequences of common agronomic practices on epidemics. Our study highlights the importance of exploring vectors interactions to enhance the understanding of disease transmission, supporting the development of disease management strategies. TCAB1 prevents nucleolar accumulation of the telomerase RNA to promote telomerase assembly Localization of a wide variety of RNAs to non-membrane bound cellular compartments such as nucleoli, Cajal bodies, and stress-granules is critical for their function and stability. The molecular mechanisms that underly the recruitment and exclusion of specific RNAs from these phase-separated organelles is poorly understood. Telomerase is a ribonucleoprotein (RNP), that is composed of the reverse transcriptase protein TERT, the telomerase RNA (TR), and several auxiliary proteins that associate with TR, including TCAB1. Here we show that, that in the absence of TCAB1, TR is sequestered in the nucleolus, while TERT localizes to the nucleoplasm and is excluded from the nucleolus, which prevents telomerase assembly. Thus, nuclear compartmentalization by the non-membrane bound nucleolus counteracts telomerase assembly and TCAB1 is required to exclude the telomerase RNA from the nucleolus. Our work provides general insight into the mechanism and functional consequences of RNA recruitment to organelles formed by phase-separation and proposes a new model explaining the critical role of TCAB1 in telomerase function. Omicron BA.2 specifically evades broad sarbecovirus neutralizing antibodies Omicron sub-lineage BA.2 has rapidly surged globally, accounting for over 60% of recent SARS-CoV-2 infections. Newly acquired RBD mutations and high transmission advantage over BA.1 urge the investigation of BA.2s immune evasion capability. Here, we show that BA.2 causes strong neutralization resistance, comparable to BA.1, in vaccinated individuals plasma. However, BA.2 displays more severe antibody evasion in BA.1 convalescents, and most prominently, in vaccinated SARS convalescents plasma, suggesting a substantial antigenicity difference between BA.2 and BA.1. To specify, we determined the escaping mutation profiles1,2 of 714 SARS-CoV-2 RBD neutralizing antibodies, including 241 broad sarbecovirus neutralizing antibodies isolated from SARS convalescents, and measured their neutralization efficacy against BA.1, BA.1.1, BA.2. Importantly, BA.2 specifically induces large-scale escape of BA.1/BA.1.1-effective broad sarbecovirus neutralizing antibodies via novel mutations T376A, D405N, and R408S. These sites were highly conserved across sarbecoviruses, suggesting that Omicron BA.2 arose from immune pressure selection instead of zoonotic spillover. Moreover, BA.2 reduces the efficacy of S309 (Sotrovimab)3,4 and broad sarbecovirus neutralizing antibodies targeting the similar epitope region, including BD55-5840. Structural comparisons of BD55-5840 in complexes with BA.1 and BA.2 spike suggest that BA.2 could hinder antibody binding through S371F-induced N343-glycan displacement. Intriguingly, the absence of G446S mutation in BA.2 enabled a proportion of 440-449 linear epitope targeting antibodies to retain neutralizing efficacy, including COV2-2130 (Cilgavimab)5. Together, we showed that BA.2 exhibits distinct antigenicity compared to BA.1 and provided a comprehensive profile of SARS-CoV-2 antibody escaping mutations. Our study offers critical insights into the humoral immune evading mechanism of current and future variants. Gene structure-based homology search identifies highly divergent putative effector gene family Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors. Sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall- forming aphid, Tetraneura ulmi, and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of Hormaphis cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. 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. 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. sgcocaller and comapr: personalised haplotype assembly and comparative crossover map analysis using single-gamete sequencing data Profiling gametes of an individual enables the construction of personalised haplotypes and meiotic crossover landscapes, now achievable at larger scale than ever through the availability of high-throughput single-cell sequencing technologies. However, high-throughput single-gamete data commonly have low depth of coverage per gamete, which challenges existing gametebased haplotype phasing methods. In addition, haplotyping a large number of single gametes from high-throughput singlecell DNA sequencing data and constructing meiotic crossover profiles using existing methods requires intensive processing. Here, we introduce efficient software tools for the essential tasks of generating personalised haplotypes and calling crossovers in gametes from single-gamete DNA sequencing data (sgcocaller), and constructing, visualising, and comparing individualised crossover landscapes from single gametes (comapr). With additional data pre-possessing, the tools can also be applied to bulk-sequenced samples. We demonstrate that sgcocaller is able to generate impeccable phasing results for high-coverage datasets, on which it is more accurate and stable than existing methods, and also performs well on low-coverage single-gamete sequencing datasets for which current methods fail. Our tools achieve highly accurate results with user-friendly installation, comprehensive documentation, efficient computation times and minimal memory usage. 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. 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. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS.