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.
T606-phosphorylation deprives the function of Kaiso as a transcription factor and cancer repressor It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer.
A meta-analysis of transcriptomic profiles from Huntington's disease patients points to a pathophysiological role of CDC42, NFY, DLX1 and PRMT3 Description of robust transcriptomic alterations in Huntingtons disease is essential to identify targets for biochemical studies and drug development. We analysed publicly available transcriptome data from the brain and blood of 220 HD patients and 241 healthy controls and identified 737 and 661 genes with robustly altered mRNA levels in the brain and blood of HD patients, respectively. In the brain, a subnetwork of 320 genes strongly correlated with HD and was enriched in transport-related genes. Bioinformatical analysis of this subnetwork highlighted CDC42, PAK1, YWHAH, NFY, DLX1, HMGN3, and PRMT3. Moreover, we found that CREB1 can regulate 78.0 % of genes whose mRNA levels correlated with HD in the blood of patients. Alterations in protein transport, metabolism, transcriptional regulation, and CDC42-mediated functions are likely central features of HD. Further our data substantiate the role of transcriptional regulators that have not been reported in the context of HD (e.g. DLX1, HMGN3 and PRMT3) and strongly suggest dysregulation of NFY and its target genes across tissues. A large proportion of the identified genes such as CDC42 were also altered in Parkinsons (PD) and Alzheimers disease (AD). The observed dysregulation of CDC42 and YWHAH in samples from HD, AD and PD patients indicates that those genes and their upstream regulators may be interesting therapeutic targets. Spatial organization of nuclear pores in Xenopus laevis oocytes Nuclear pores are protein assemblies inserted in the nuclear envelope of eukaryotic cells, acting as main gates for communication between nucleus and cytoplasm. So far, nuclear pores have been extensively studied to determine their structure and composition, yet their spatial organization and geometric arrangement on the nuclear surface are still poorly understood. Here, we analyze super-resolution images of the surface of Xenopus laevis oocyte nuclei during development, and characterize the arrangement of nuclear pores using tools commonly employed to study the atomic structural and topological features of soft matter. To interpret the experimental results, we hypothesize an effective interaction among nuclear pores and implemented it in extensive numerical simulations of octagonal clusters mimicking typical pore shapes. Thanks to our simple model, we find simulated spatial distributions of nuclear pores that are in excellent agreement with experiments, suggesting that an effective interaction among nuclear pores exists and could explain their geometrical arrangement. Furthermore, our results show that the statistical features of the geometric arrangement of nuclear pores do not depend on the type of pore-pore interaction, attractive or repulsive, but are mainly determined by the octagonal symmetry of each single pore. These results pave the way to further studies needed to determine the biological nature of pore-pore interactions. Implicit Laplacian of Enhanced Edge: An Unguided Algorithm for Accurate and Automated Quantitative Analysis of Cytoskeletal Images The eukaryotic cytoskeleton plays essential roles in cell signaling and trafficking, which is broadly associated with immunity and diseases of human and plants. To date, most analyses aiming at defining the temporal and spatial dynamics of the cytoskeleton have relied on qualitative and quantitative analysis of fluorescence images to describe cytoskeletal function. While state-of-the-art, these approaches have limitations: the diverse shape and brightness of the cytoskeleton cause considerable inaccuracy in both human-driven and automated approaches, and the widely adopted image projection process (3D to 2D) leads to bias and information loss. Here, we describe the development and application of Implicit Laplacian of Enhanced Edge (ILEE), an unguided approach that uses a 2D/3D-compatible local thresholding algorithm for the quantitative evaluation of cytoskeletal status and organization at high performance. Using ILEE, we constructed a Python library to enable automated cytoskeletal image analysis, providing numerous biologically-interpretable indices measuring the density, bundling, severing, branching, and directionality of the cytoskeleton. The data presented herein demonstrate that ILEE resolves the defects of classic cytoskeleton analysis approaches, enables the measurement of novel cytoskeletal features, and yields quantitatively descriptive data with superior accuracy, stability, and robustness. We released the ILEE algorithm as an open-source library and further developed a Google Colab interface as a community resource. The RAG1 Ubiquitin Ligase Domain Enhances Assembly and Selection of T Cell Receptor Genes to Restrain the Autoimmune Hazard of Generating T Cell Receptor Diversity RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte antigen receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. However, this process creates highly self-reactive cells that must be properly selected to suppress autoimmunity. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report that mice with a mutation that inactivates the RAG1 ubiquitin ligase in vitro exhibit modestly reduced thymic cellularity, decreased assembly and altered repertoires of T cell receptor (TCR) {beta} and  genes in thymocytes, and impaired thymocyte developmental transitions that require the assembly of TCR{beta} or  genes and signaling by their proteins. These RAG1 mutant mice also exhibit less efficient positive selection and superantigen-mediated negative selection of conventional {beta} T cells, 2) impaired differentiation of iNKT lineage {beta} T cells, and 3) CD4+ {beta} T cells with elevated autoimmune potential. Our findings demonstrate that the RAG1 ubiquitin ligase domain functions in vivo to stimulate the assembly and selection of TCR{beta} and TCR genes, thereby establishing replete diversity of {beta} TCRs and {beta} T cell lineages while restraining the inherent autoimmune hazard of generating diverse antigen specificities. Niche differences, not fitness differences, explain coexistence across ecological groups Understanding the drivers of species coexistence is an important objective in ecology. Yet, the multitude of methods to study coexistence hampers cross-community comparisons. Here, we standardized niche and fitness differences (i.e how species limit themselves compared to others and their competitive ability, respectively) across 1018 species pairs to investigate species coexistence across ecological groups and methodological settings (experimental setup, natural co-occurrence, population model used, and growth method). We find that, first, coexistence is driven by large niche differences, not by small fitness differences. Second, species group into clear clusters of coexisting and non-coexisting species along the niche axis. Finally, these clusters are not driven by ecological or methodological settings. This suggests differences between coexisting and non-coexisting communities transcending those measured in our empirical systems. Overall, our results show that species coexistence is mainly influenced by mechanisms acting on niche differences. Design of a radial multi-offset detection pattern for in vivo phase contrast imaging of the retinal ganglion cells in humans Previous work has shown that multi-offset detection in adaptive optics scanning laser ophthalmoscopy (AOSLO) can be used to image retinal ganglion cells (RGCs) in monkeys and humans. However, though images of RGCs in anesthetized monkeys with high light levels produced high contrast images of RGCs, images from humans failed to reach the same contrast due to several drawbacks in the previous dual-wavelength multi-offset approach. Our aim here was to design and build a multi-offset detection pattern for humans at safe light levels that could reveal the retinal ganglion cell layer neurons with a contrast, robustness and acquisition time approaching results only previously obtained in monkeys. Here, we present a new imaging system using only one light source, compared to the previous dual-wavelength used on monkeys. Our single-wavelength solution allows for increased light power and eliminates problematic chromatic aberrations. Then, we demonstrate that a radial multi-offset detection pattern with an offset distance of 8-10 Airy Disk Diameter (ADD) is optimal to detect photons multiply scattered in all directions from RGCs thereby enhancing their contrast. This new setup and image processing pipeline led to improved imaging of retinal ganglion cells using multi-offset imaging in AOSLO. A multi-dimensional integrative scoring framework for predicting functional variants in the human genome Attempts to identify and prioritize functional DNA elements in coding and noncoding regions, particularly through use of in silico functional annotation data,continue to increase in popularity. However, specific functional roles may vary widely from one variant to another, making it challenging to summarize different aspects of variant function. Here we propose Multi-dimensional Annotation Class Integrative Estimation (MACIE), an unsupervised multivariate mixed model framework capable of integrating annotations of diverse origin to assess multi-dimensional functional roles for both coding and noncoding variants. Unlike existing one-dimensional scoring methods, MACIE views variant functionality as a composite attribute encompassing multiple characteristics, and estimates the joint posterior functional probability vector of each genomic position, a quantity that offers richer and more interpretable information in the presence of multiple aspects of functionality. Applied to a variety of independent coding and non-coding datasets, MACIE demonstrates powerful and robust performance in discriminating between functional and non-functional variants. We also show an application of MACIE to fine-mapping using lipids GWAS summary statistics data from the European Network for Genetic and Genomic Epidemiology Consortium. Experimental evolution in TP53 deficient human gastric organoids recapitulates tumorigenesis The earliest events during human tumor initiation are poorly characterized but may hold clues as to how to detect and prevent malignancy. Here we model occult pre-neoplasia by engineering TP53 deficiency, a common and early event in gastric cancer, into human gastric organoids. By performing experimental evolution in multiple clonally derived cultures over two years we define causal relationships between this common initiating genetic lesion and resulting phenotypes. TP53 loss elicited progressive aneuploidy, including copy number alterations and complex structural variants that are common in gastric cancers and which follow preferred temporal orders. Longitudinal single cell sequencing of TP53 deficient gastric organoids similarly indicates progression towards malignant transcriptional programs. Moreover, lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programs repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a striking degree of phenotypic convergence in pre-malignant epithelial organoids, implying that the earliest stages of tumorigenesis may be predictable while illuminating evolutionary constraints and barriers to malignant transformation.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/487529v2_ufig1.gif" ALT="Figure 1">
View larger version (63K):
org.highwire.dtl.DTLVardef@122debborg.highwire.dtl.DTLVardef@4cd2daorg.highwire.dtl.DTLVardef@1a720d6org.highwire.dtl.DTLVardef@15ccaa7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Centimeter-scale quorum sensing dictates collective survival of differentiating embryonic stem cells Cells can help each other replicate by communicating with diffusible molecules. In cell cultures, molecules may diffuse within a cell colony or between adjacent or distant colonies. Determining which cell helps which cells replication is challenging. We developed a systematic approach, integrating modeling and experiments, for determining the length-scales of cell-cell communication (from microns to centimeters). With this approach, we discovered that differentiating murine ES cells, scattered across centimeters on a dish, communicate over millimeters to form one macroscopic entity that survives if and only if its centimeter-scale population-density is above a threshold value. Single-cell-level measurements, transcriptomics, and modeling revealed that this "macroscopic quorum sensing" arises from differentiating ES cells secreting and sensing survival-promoting FGF4 that diffuses over millimeters and activates YAP1-induced survival mechanisms. Through the same mechanism, a lone macroscopic, but not microscopic, colony survives differentiation. Our work rigorously establishes that in vitro ES-cell differentiation relies on macroscopic cooperation. Stochastic effects during the establishment of β-lactam resistant E. colimutants indicate conditions for collective resistance For antibiotic resistance to arise, new resistant mutants must establish in a bacterial population before they can spread via natural selection. Comprehending the stochastic factors that influence mutant establishment is crucial for a quantitative understanding of antibiotic resistance emergence. Here, we quantify the single-cell establishment probability of four Escherichia coli strains expressing {beta}-lactamase alleles with different activity against the antibiotic cefotaxime, as a function of antibiotic concentration in both unstructured (liquid) and structured (agar) environments. We show that concentrations well below the minimum inhibitory concentration (MIC) can substantially hamper establishment, particularly for highly resistant mutants. While the pattern of establishment suppression is comparable in both tested environments, we find greater variability in establishment probability on agar. Using a simple branching model, we investigate possible sources of this stochasticity, including environment-dependent lineage variability. Lastly, we use the single-cell establishment probability to predict each strains MIC in the absence of social interactions. We observe substantially higher measured than predicted MIC values, particularly for highly resistant strains, which indicates cooperative effects among resistant cells at large cell numbers, such as in standard MIC assays. 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. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. 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. 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. Updated vaccine protects from infection with SARS-CoV-2 variants, prevents transmission and is immunogenic against Omicron in hamsters Current first-generation COVID-19 vaccines are based on prototypic spike sequences from ancestral 2019 SARS-CoV-2 strains. However, the ongoing pandemic is fueled by variants of concern (VOC) that threaten to escape vaccine-mediated protection. Here we show in a stringent hamster model that immunization using prototypic spike expressed from a potent YF17D viral vector (1) provides vigorous protection against infection with ancestral virus (B lineage) and VOC Alpha (B.1.1.7), however, is insufficient to provide maximum protection against the Beta (B.1.351) variant. To improve vaccine efficacy, we created a revised vaccine candidate that carries an evolved spike antigen. Vaccination of hamsters with this updated vaccine candidate provides full protection against intranasal challenge with all four VOCs Alpha, Beta, Gamma (P.1) and Delta (B.1.617.2) resulting in complete elimination of infectious virus from the lungs and a marked improvement in lung pathology. Vaccinated hamsters did also no longer transmit the Delta variant to non-vaccinated sentinels. Hamsters immunized with our modified vaccine candidate also mounted marked neutralizing antibody responses against the recently emerged Omicron (B.1.1.529) variant, whereas the old vaccine employing prototypic spike failed to induce immunity to this antigenically distant virus. Overall, our data indicate that current first-generation COVID-19 vaccines need to be urgently updated to cover newly emerging VOCs to maintain vaccine efficacy and to impede virus spread at the community level.

Significance StatementSARS-CoV-2 keeps mutating rapidly, and the ongoing COVID-19 pandemic is fueled by new variants escaping immunity induced by current first-generation vaccines. There is hence an urgent need for universal vaccines that cover variants of concern (VOC). In this paper we show that an adapted version of our vaccine candidate YF-S0* provides full protection from infection, virus transmission and disease by VOCs Alpha, Beta, Gamma and Delta, and also results in markedly increased levels of neutralizing antibodies against recently emerged Omicron VOC in a stringent hamster model. Our findings underline the necessity to update COVID-19 vaccines to curb the pandemic, providing experimental proof on how to maintain vaccine efficacy in view of an evolving SARS-CoV-2 diversity. Modifying root/shoot ratios improves root water influxes in wheat under drought stress Drought intensity as experienced by plants depends upon soil moisture status and atmospheric variables such as temperature, radiation, and air vapour pressure deficit (VPD). Although the role of shoot architecture with these edaphic and atmospheric factors is well-characterized, the extent to which shoot and root dynamic interactions as a continuum are controlled by genotypic variation is less known. Here, we targeted these interactions using a wild emmer introgression line (IL20) with a distinct drought-induced shift in the shoot-to-root ratio and its drought-sensitive recurrent parent Svevo. Using a gravimetric platform, we show that IL20 maintained higher root water influx and gas exchange under drought stress, which supported a greater growth. Interestingly, the advantage of IL20 in root water influx and transpiration was expressed earlier during the daily diurnal cycle under lower VPD and therefore supported higher transpiration efficiency. Application of structural equation model indicates that under drought, VPD and radiation are antagonistic to transpiration rate, whereas the root water influx operates as feedback for the higher atmospheric responsiveness of leaves. Collectively, our results suggest that a drought-induced shift in root-to-shoot ratio can improve plant water uptake potential in a short preferable time window determined by both water and atmospheric parameters. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. A Positive Correlation between Bacterial GC Content and Growth Temperature BackgroundGC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate.

ResultsWith a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GCw) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GCw and growth temperature in bacteria. By including the GCw of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GCw and Topt appear especially after excluding the halophilic archaea whose GC contents might be are strongly shaped by intense UV radiation.

ConclusionsThis study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton. SARS-CoV-2 Beta Variant Substitutions Alter Spike Glycoprotein Receptor Binding Domain Structure and Stability The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic has visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations; and emerging variants of concern (VOC) are more infectious and reduce the efficacy of existing vaccines. While promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the RBD structure, stability, and ability to bind to ACE2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these VOC may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain. 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. Integrated annotation and analysis of genomic features reveal new types of functional elements and large-scale epigenetic phenomena in the developing zebrafish Zebrafish, a popular model for embryonic development and for modelling human diseases, has so far lacked a systematic functional annotation programme akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created the first central repository to store and process zebrafish developmental functional genomic data. Our Data Coordination Center (https://danio-code.zfin.org) combines a total of 1,802 sets of unpublished and reanalysed published genomics data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements in development, including novel classes with distinct features dependent on their activity in time and space. We delineated the distinction between regulatory elements active during zygotic genome activation and those active during organogenesis, identifying new aspects of how they relate to each other. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predict functional relationships between them beyond sequence similarity, extending the utility of zebrafish developmental genomics to mammals. F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. CBEA: Competitive balances for taxonomic enrichmentanalysis Research in human associated microbiomes often involves the analysis of taxonomic count tables generated via high-throughput sequencing. It is difficult to apply statistical tools as the data is high-dimensional, sparse, and compositional. An approachable way to alleviate high-dimensionality and sparsity is to aggregate variables into pre-defined sets. Set-based analysis is ubiquitous in the genomics literature, and has demonstrable impact in improving interpretability and power of downstream analysis. Unfortunately, there is a lack of sophisticated set-based analysis methods specific to microbiome taxonomic data, where current practice often employs abundance summation as a technique for aggregation. This approach prevents comparison across sets of different sizes, does not preserve inter-sample distances, and amplifies protocol bias. Here, we attempt to fill this gap with a new single sample taxon enrichment method that uses a novel log-ratio formulation based on the competitive null hypothesis commonly used in the enrichment analysis literature. Our approach, titled competitive balances for taxonomic enrichment analysis (CBEA), generates sample-specific enrichment scores as the scaled log ratio of the subcomposition defined by taxa within a set and the subcomposition defined by its complement. We provide sample-level significance testing by estimating an empirical null distribution of our test statistic with valid p-values. Herein we demonstrate using both real data applications and simulations that CBEA controls for type I error even under high sparsity and high inter-taxa correlation scenarios. Additionally, it provides informative scores that can be inputs to downstream analyses such as prediction tasks.

Author summaryThe study of human associated microbiomes relies on genomic surveys via high-throughput sequencing. However, microbiome taxonomic data is sparse and high dimensional which prevents the application of standard statistical techniques. One approach to address this problem is to perform analyses at the level of taxon sets. Set-based analysis has a long history in the genomics literature, with demonstrable impact in improving both power and interpretability. Unfortunately, there is limited interest in developing new set-based tools tailored for microbiome taxonomic data given its unique features compared to other  omics data types. We developed a new tool to generate taxon set enrichment scores at the sample level through a novel log-ratio formulation based on the competitive null hypothesis. Our scores can be used for statistical inference at both the sample and population levels, as well as inputs to other downstream analyses such as prediction models. We demonstrate the performance of our method against competing approaches across both real data analyses and simulation studies. A terrain treadmill to study animal locomotion through large obstacles A major challenge to understanding locomotion in complex 3-D terrain with large obstacles is to create tools for controlled, systematic lab experiments. Existing terrain arenas only allow observations at small spatiotemporal scales ([~]10 body length, [~]10 stride cycles). Here, we create a terrain treadmill to enable high-resolution observations of animal locomotion through large obstacles over large spatiotemporal scales. An animal moves through modular obstacles on an inner sphere, while a rigidly-attached, concentric, transparent outer sphere rotated with the opposite velocity via closed-loop feedback to keep the animal on top. During sustained locomotion, a discoid cockroach moved through pillar obstacles for 25 minutes ({approx} 2500 strides) over 67 m ({approx} 1500 body lengths), and was contained within a radius of 4 cm (0.9 body length) for 83% of the duration, even at speeds of up to 10 body length/s. The treadmill enabled observation of diverse locomotor behaviors and quantification of animal-obstacle interaction.

Summary StatementA device keeps a small animal on top of a sphere while it traverses obstacles, creating a "terrain treadmill" to study locomotion over a long time and distance at high-resolution.