HOW CONCERNING IS A SARS-COV-2 VARIANT OF CONCERN? COMPUTATIONAL PREDICTIONS AND THE VARIANTS LABELING SYSTEM We herein report a study to evaluate the use of computational prediction of SARS-CoV-2 genetic variations in improving the current variants labeling system. First, we reviewed the basis of the system developed by the World Health Organization (WHO) for the labeling of SARS-CoV-2 genetic variants and the adaptations made to it by the United States Center of Diseases Control (CDC). We observed that the labeling system is based upon the virus major attributes. However, we found that the labeling criteria of the SARS-CoV-2 variants derived from these attributes are not accurately defined and are used differently by the two health management agencies. Consequently, discrepancies exist between the labels given by WHO and CDC to same variants. Our observations suggest that giving the VOC label to a new variant is premature and might not be appropriate. Therefore, we carried out a comparative computational study to predict the effects of the mutations on the virus structure and functions of five VOCs. By linking these data to the criteria used by WHO and the CDC for variant labeling, we ascertained that comparative computational predictions of the impact of genetic variations are a better ground for rapid and more accurate labelling of SARS-CoV-2 variants. We propose to label all emergent variants VUM or VBM and to carry out computational predictive studies and thorough variants comparison, upon which more appropriate and informative labels can be attributed. Furthermore, harmonization of the variants labeling system would be globally beneficial to communicate about and fight COVID-19 pandemic.
Identification and characterization of the WYL BrxR protein and its gene asseparable regulatory elements of a BREX phage restriction system Bacteriophage exclusion ( BREX) phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR ( BREX Regulator), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that the BrxR coding sequence is involved in cis regulation of anti-phage activity, while the BrxR protein plays an additional regulatory role, perhaps during horizontal transfer.
Key roles for kinesin-13 and kinesin-20 in malaria parasite proliferation, polarity, and transmission revealed by genome-wide functional analysis Kinesins are microtubule-based motors important in cell division, motility, polarity, and intracellular transport in many eukaryotes. However, they are poorly studied in the divergent eukaryotic pathogens-Plasmodium spp., the causative agents of malaria, which manifest atypical aspects of cell division and plasticity of morphology throughout the lifecycle in both mammalian and mosquito hosts. Here we describe a genome-wide screen of Plasmodium kinesins, revealing diverse subcellular locations and functions in spindle assembly, axoneme formation and cell morphology. Surprisingly, only kinesin-13 is essential for growth in the mammalian host while the other eight kinesins are required during the proliferative and invasive stages of parasite transmission through the mosquito vector. In-depth analyses of kinesin-13 and kinesin-20 revealed functions in microtubule dynamics during apical cell polarity formation, spindle assembly, and axoneme biogenesis. These findings help us to understand the importance of microtubule motors and may be exploited to discover new therapeutic interventions against malaria. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins - the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening. Coronavirus-associated molecular mimicry through homology to a SARS-CoV-2 peptide could be leading to susceptibility in patients with HLA-A*02:01 and HLA-A*24:02 serotypes This study aims to predict autoimmunity-related pathological mechanisms that possess risk for individuals with specific human leukocyte antigen (HLA) serotypes and shared by certain coronaviruses including SARS-CoV-2, based on homology to a SARS-CoV-2 peptide. With the given aim, 1-) coronavirus-associated sequences, which are homologous to the 15mer SARS-CoV-2 peptide CFLGYFCTCYFGLFC, are obtained. 2-) Human peptides that have at least 7 residue matches with those coronavirus sequences, and the SARS-CoV-2 15mer, are found. 3-) Epitope pairs, which are sourced by those aligned coronavirus and human sequences are identified. 4-) Epitope pairs that are predicted to bind strongly not only to the same HLA allele with each other but also to the same HLA allele as those of the respective alignment of the SARS-CoV-2 peptide are selected. Following are the identified proteins or peptides (with HLA-A*02:01 or HLA-A*24:02 epitopes), as described in 1-to-4: Immunoglobulin heavy chain junction regions, CRB1 isoform I precursor, slit homolog 2 protein, hCG1995581, hCG2028737, phospholipid phosphatase-related protein type 2. Among those, CRB1 isoform I precursor sequence with the predicted HLA-A*24:02 epitope aligns with the highest number of different sequences. Results imply autoimmunity risk in COVID-19 patients with HLA-A*02:01 and HLA-A*24:02 serotypes, through molecular mimicry, as a shared pathogenicity risk that can be prevalent upon getting infected with certain coronaviruses. These can pave way to improved risk groups assessment and autoimmunity treatment options, for COVID-19 and its associated diseases. Also, the approach in this study can be used to predict prospective pathologies of the transmissible variants in susceptible humans. Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann-fold proteins The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved {beta}{beta} motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-Adenosylmethionine (SAM), the oxidoreductases, depending on the family, bind nicotinamide (NAD, NADP) or flavin-based (FAD) cofactors. In this study, we show that despite its short length, the {beta}{beta} motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the {beta}{beta} motif. A benchmark on two independent test sets, one containing {beta}{beta} motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

Key pointsO_LIThe Rossmann fold encompasses a multitude of diverse enzymes involved in most of the essential cellular pathways
C_LIO_LIProteins belonging to the Rossmann fold co-evolved with their nucleoside-based cofactors and require them for the functioning
C_LIO_LIManipulating the cofactor specificity is an important step in the process of enzyme engineering
C_LIO_LIWe developed an end-to-end pipeline for the prediction and design of the cofactor specificity of the Rossmann fold proteins
C_LIO_LIOwing to the utilization of deep learning approaches the pipeline achieved nearly perfect accuracy
C_LI Diurnal activity budget and conservation implication of Hippopotamus amphibius at Bui National Park, Ghana Animals apportion time for their daily behavioral activities. We studied the activity budget of Hippopotamus amphibius at the Black Volta River in the Bui National Park, Ghana. We performed instantaneous scan sample using ground count survey between August 2020 to July 2021. We observed that H. amphibius spent most of their day time resting (54.75 %), followed by feeding (22.93 %), walking (19.2 %) and touching (3.12 %). We found significant difference in the time spent between the daily activities (H = 41.67; p < 0.0001). Except activities involving touching, we found no significant seasonal differences in feeding, resting and walking activities by H. amphibius. Understanding the behavioral activities of H. amphibius will enhance management and conservation of the animals. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Differential regulation of degradation and immune pathways underlies adaptation of the ectosymbiotic nematode Laxus oneistus to oxic-anoxic interfaces Eukaryotes may experience oxygen deprivation under both physiological and pathological conditions. Because oxygen shortage leads to a reduction in cellular energy production, all eukaryotes studied so far conserve energy by suppressing their metabolism. However, the molecular physiology of animals that naturally and repeatedly experience anoxia is underexplored. One such animal is the marine nematode Laxus oneistus. It thrives, invariably coated by its sulfur-oxidizing symbiont Candidatus Thiosymbion oneisti, in anoxic sulfidic or hypoxic sand. Here, transcriptomics and proteomics showed that, whether in anoxia or not, L. oneistus mostly expressed genes involved in ubiquitination, energy generation, oxidative stress response, immune response, development, and translation. Importantly, ubiquitination genes were also upregulated when the nematode was subjected to anoxic sulfidic conditions, together with genes involved in autophagy, detoxification and ribosome biogenesis. We hypothesize that these degradation pathways were induced to recycle damaged cellular components (mitochondria) and misfolded proteins into nutrients. Remarkably, when L. oneistus was subjected to anoxic sulfidic conditions, lectin and mucin genes were also upregulated, potentially to promote the attachment of its thiotrophic symbiont. Furthermore, the nematode appeared to survive oxygen deprivation by using an alternative electron carrier (rhodoquinone) and acceptor (fumarate), to rewire the electron transfer chain. On the other hand, under hypoxia, genes involved in costly processes (e.g., amino acid biosynthesis, development, feeding, mating) were upregulated, together with the worms Toll- like innate immunity pathway and several immune effectors (e.g., Bacterial Permeability Increasing proteins, fungicides).

In conclusion, we hypothesize that, in anoxic sulfidic sand, L. oneistus upregulates degradation processes, rewires oxidative phosphorylation and by reinforces its coat of bacterial sulfur-oxidizers. In upper sand layers, instead, it appears to produce broad-range antimicrobials and to exploit oxygen for biosynthesis and development. RNA degradation of blood reveals a time diagnostic window in a controlled environment Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. RNAi screen of RING/U-box domain ubiquitin ligases identifies critical regulators of tissue regeneration in planarians Regenerative processes depend on the interpretation of signals to coordinate cell behaviors. The role of ubiquitin-mediated signaling is known to be important in many cellular and biological contexts, but its role in regeneration is not well understood. To investigate how ubiquitylation impacts tissue regeneration in vivo, we are studying planarians that are capable of regenerating after nearly any injury using a population of stem cells. Here we used RNAi to screen RING/U-box E3 ubiquitin ligases that are highly expressed in planarian stem cells and stem cell progeny. RNAi screening identified nine genes with functions in regeneration, including the spliceosomal factor prpf19 and histone modifier rnf2; based on their known roles in developmental processes, we further investigated these two genes. We found that prpf19 was required for animal survival but not for stem cell maintenance, suggesting a role in promoting cell differentiation. Because RNF2 is the catalytic subunit of the Polycomb Repressive Complex 1 (PRC1), we also examined other putative members of this complex (CBX and PHC). We observed a striking phenotype of regional tissue misspecification in cbx and phc RNAi planarians. To identify genes regulated by PRC1, we performed RNA-seq after knocking down rnf2 or phc. Although these proteins are predicted to function in the same complex, we found that the set of genes differentially expressed in rnf2 versus phc RNAi were largely non-overlapping. Using in situ hybridization, we showed that rnf2 regulates gene expression levels within a tissue type, whereas phc is necessary for the spatial restriction of gene expression, findings consistent with their respective in vivo phenotypes. This work not only uncovered roles for RING/U-box E3 ligases in stem cell regulation and regeneration, but also identified differential gene targets for two putative PRC1 factors required for maintaining cell-type-specific gene expression in planarians. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. Low photorespiratory capacity is sufficient to block the induction of photosynthetic capacity in high light The induction of high photosynthetic capacity in high light (HL) is a common response among many herbaceous dicot plants, however, the signals that control this response remain largely unknown. Here, multiple independent lines of evidence are presented in support of the conclusion that low photorespiratory capacity acts a negative signal to limit photosynthetic capacity acclimation in HL in Arabidopsis thaliana. Using a panel of natural accessions, primary nitrogen (N) assimilation and photorespiration rates early after a shift to growth in HL, as well as activities for key enzymes in these pathways, were shown to positively correlate with the magnitude of the subsequent induction of photosynthetic capacity, which occurred several days later. Time-resolved metabolomic data during acclimation to HL were collected using a strongly acclimating ecotype and a weakly acclimating ecotype, revealing in greater detail the differences in N assimilation, photorespiration, and triose-phosphate utilization pathways underlying efficient photosynthetic capacity acclimation. When shifted into HL growth conditions under non-photorespiratory conditions, weakly acclimating ecotypes and even photorespiratory mutants gained the ability to strongly induce high photosynthetic capacity in HL. Thus, a negative, photorespiration-dependent signal early in the HL shift appears to block photosynthetic capacity acclimation in accessions with low photorespiratory capacity, whereas accessions with high photorespiratory capacity are licensed to increase photosynthetic capacity. Rapid encoding of task regularities in the human hippocampus guides sensorimotor timing The brain encodes the statistical regularities of the environment in a task-specific yet flexible and generalizable format. Here, we seek to understand this process by converging two parallel lines of research, one centered on sensorimotor timing, and the other on cognitive mapping in the hippocampal system. By combining functional magnetic resonance imaging (fMRI) with a fastpaced time-to-contact (TTC) estimation task, we found that the hippocampus signaled behavioral feedback received in each trial as well as performance improvements across trials along with reward-processing regions. Critically, it signaled performance improvements independent from the tested intervals, and its activity accounted for the trial-wise regression-to-the-mean biases in TTC estimation. This suggests that the hippocampus supports the rapid encoding of temporal context even on short time scales in a behavior-dependent manner. Our results emphasize the central role of the hippocampus in statistical learning and position it at the core of a brain-wide network updating sensorimotor representations in real time for flexible behavior. Parallel Characterization of cis-Regulatory Elements for Multiple Genes Using CRISPRpath Current pooled CRISPR screens for cis-regulatory elements (CREs) can only accommodate one gene based on its expression level. Here, we describe CRISPRpath, a scalable screening strategy for parallelly characterizing CREs of genes linked to the same biological pathway and converging phenotypes. We demonstrate the ability of CRISPRpath for simultaneously identifying functional enhancers of six genes in the 6-thioguanine-induced DNA mismatch repair pathway using both CRISPR interference (CRISPRi) and CRISPR nuclease (CRISPRn) approaches. 60% of the identified enhancers are known promoters with distinct epigenomic features compared to other active promoters, including increased chromatin accessibility and interactivity. Furthermore, by imposing different levels of selection pressure, CRISPRpath can distinguish enhancers exerting strong impact on gene expression from those exerting weak impact. Our results offer a nuanced view of cis-regulation and demonstrate that CRISPRpath can be leveraged for understanding the complex gene regulatory program beyond transcriptional output at scale. Competition and environmental gradients structure assemblages of Canidae across the planet The phylogenetic information of assemblages carries the signature of ecological and evolutionary processes that assembled these communities. Identifying the mechanisms that shape communities is not simple, as they can vary spatially. Here, we investigated how the phylogenetic structure of Canidae across the globe is affected by the environment and competition. We first identified phylogenetically clustered and overdispersed assemblages of canids over the planet taking into account regional pools of species in South and North America, Eurasia, and Africa. Then, we applied Structural Equation Models in these communities in order to identify the effect of temperature, vegetation cover, human impact, and body size dissimilarity on the spatial distribution of the phylogenetic information of canids. Only southern South America and the Middle East had strong phylogenetic clustered patterns, while the rest of the planet predominantly showed overdispersed assemblages. Body size dissimilarity and vegetation cover were the most important variables to explain the spatial patterns of phylogenetic structure in clustered and overdispersed communities, respectively. Canidae community composition across the world presents significant patterns of clustering and overdispersion, which vary following mainly the environmental gradient, suggesting habitat filtering as the main force acting on Canidae assemblages. Trypanosomal variant surface glycoprotein expression in human African trypanosomiasis patients Trypanosoma brucei gambiense, an extracellular protozoan parasite, is the primary causative agent of human African trypanosomiasis. T. b. gambiense is endemic to West and Central Africa, where it is transmitted by the bite of infected tsetse flies. In the bloodstream of an infected host, the parasite evades antibody recognition by altering the Variant Surface Glycoprotein (VSG) that forms a dense coat on its cell surface through a process known as antigenic variation. Each VSG has a variable N-terminal domain exposed to the host and a less variable C-terminal domain that is partially hidden from host antibodies. Our lab developed VSG-seq, a targeted RNA-seq method, to study VSG gene expression in T. brucei. Studies using VSG-seq to characterize antigenic variation in a mouse model have revealed marked diversity in VSG expression within parasite populations, but this finding has not yet been validated in a natural human infection. Here, we used VSG-seq to analyze VSGs expressed in the blood of 12 patients infected with T. b. gambiense. The number of VSGs identified per patient ranged from one to 14; notably, two VSGs were shared by more than one patient. Analysis of expressed VSG N-terminal domain types revealed that 82% of expressed VSGs encoded a type B N-terminus, a bias not seen in datasets from other T. brucei subspecies. C-terminal types in T. b. gambiense infection were also restricted. These results demonstrate a bias either in the underlying VSG repertoire of T. b. gambiense or in this subspecies expression of VSGs during infection. This work demonstrates the feasibility of using VSG-seq to study antigenic variation in human infections and highlights the importance of understanding VSG repertoires in the field.

Author SummaryHuman African Trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their dense surface coat of Variant Surface Glycoprotein (VSG). Despite the important role of VSGs in prolonging infection, VSG expression during natural human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat and improve how trypanosomiasis is diagnosed in humans. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that parasites in both natural and experimental infections express diverse sets of VSGs. 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 linkage of methionine addiction, overmethylation of histone H3 lysines and malignancy demonstrated when cancer cells revert to methionine-independence Methionine addiction is a fundamental and general hallmark of cancer and is an area of current intense interest. Methionine addiction results from the overuse of methionine by cancer cells for excess transmethylation reactions. In order to identify excess transmethylation reactions in cancer and further understand the basis of methionine addiction, we compared the histone H3 lysine-methylation status and malignancy between methionine-addicted cancer cells and their methionine-independent revertants which have regained the ability to grow on low levels of methionine or independently of exogenous methionine. The levels of trimethylated histone H3 lysine marks were reduced in methionine-independent revertants compared to parental cancer cells in vitro. Tumorigenicity and experimental metastatic potential in nude mice were also highly reduced in the methionine-independent revertants compared to the parental cells. Our present results demonstrate that overmethylation of histone H3 lysines is linked with methionine addiction of cancer and to malignancy which suggests a possible causal relationship. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. Ecological interactions in breast cancer: Cell facilitation promotes growth and survival under drug pressure The interplay of positive and negative interactions between drug sensitive and resistant cells determines the efficacy of treatment in heterogeneous cancer cell populations. We developed isogenic estrogen receptor positive (ER+) breast cancer cell lineages sensitive and resistant to ribociclib-induced CDK4/6 inhibition to investigate cell interactions under selective pressure. During treatment with ribociclib alone or in combination with endocrine therapies, sensitive cells survive and proliferate better when grown together with resistant cells than when grown in monoculture. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays detected estrone, a less active metabolite, in medium taken from each lineage but detected highly active estradiol only from resistant cells. Because sensitive cells are more dependent on estrogen signaling than resistant cells, they benefit from the increase in local estradiol from resistant cells, leading to facilitated growth when grown under treatment in coculture. We developed a suite of mathematical models to quantify the strength of competition and facilitation between cell lines and evaluate the effects of treatments that modify these interactions. These analyses support the conclusion that resistant cell lineages facilitate growth of sensitive lineages during therapy through the increased production and sharing of local estradiol and themselves suffer from increased competition and suppression as a result. These models show that blocking facilitation has the potential to control both resistant and sensitive cell populations and inhibit the emergence of a refractory population. This combination of experiments and models shows that interactions within heterogeneous cancer cell populations can combine positive and negative interactions to shape the dynamics of resistance. Priority conservation areas and a global population estimate for the Critically Endangered Philippine Eagle derived from modelled range metrics using remote sensing habitat characteristics Many range-restricted taxa are currently experiencing population declines yet lack fundamental information regarding distribution and population size. Establishing baseline estimates for both these key biological parameters is however critical for directing conservation planning for at-risk range-restricted species. The International Union for the Conservation of Nature (IUCN) Red List uses three range metrics that define species distributions and inform extinction risk assessments: extent of occurrence (EOO), area of occupancy (AOO) and area of habitat (AOH). However, calculating all three metrics using standard IUCN approaches relies on a geographically representative sample of locations, which for rare species is often spatially biased. Here, we apply model-based interpolation using Species Distribution Models (SDMs), correlating occurrences with remote-sensing covariates, to calculate IUCN range metrics, protected area coverage and a global population estimate for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi). Our final range wide continuous SDM had high predictive accuracy (Continuous Boyce Index = 0.927) and when converted to a binary model estimated an AOH = 23,185 km2, a maximum EOO = 605,759 km2, a minimum EOO = 272,272 km2, with an AOO = 53,867 km2. Based on inferred habitat from the AOH metric, we estimate a global population of 318 breeding pairs (range: 258-362 pairs), or 636 mature individuals, across the Philippine Eagle global range. Protected areas covered 34 % of AOH, 15 % less than the target representation, with the continuous model identifying key habitat as priority conservation areas. We demonstrate that even when occurrences are geographically biased, robust habitat models can be built that enable quantification of baseline IUCN range metrics, protected area coverage, and a population size estimate. In the absence of adequate location data for many rare and threatened taxa, our method is a promising spatial modelling tool with widespread applications, in particular for island endemics facing high extinction risk. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed.