Ghost lineages highly influence the interpretation of introgression tests Most species are extinct; those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression and admixture, are therefore likely to involve "ghosts", i.e. extinct, unknown or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario.
Long-term translocation and hatchery release explain population genetic structure of a recreationally fished iconic species in Japan: A review and reanalysis Ayu (Plecoglossus altivelis altivelis) is an important freshwater fisheries resource and popular recreational fishing species in Japan, lives for only one year, and has a single breeding season. To supplement increased recreational fishing demand for this species, huge numbers of wild-born landlocked juveniles have been translocated from Lake Biwa into most Japanese rivers for more than 100 generations. Hatchery-reared juveniles born from captive-reared parents for more than 30 generations have also been extensively released. Previous studies have reported that landlocked and amphidromous forms of Ayu easily hybridise, but survival of landlocked larvae could be very low in seawater, leading to a general consensus among scientists, hatchery managers, and commercial and recreational fishers that the reproductive success of landlocked Ayu is very low, or even 0 in translocated rivers. Despite this, limited information exists regarding the reproductive success of landlocked Ayu in translocated rivers, and no study has evaluated the effects of translocation on population structure. Demonstrating that hybridisation occurs between the two forms is central to future management and conservation of this specie. To address this issue, a comprehensive literature search was undertaken, and three published genetic data sets are analysed. Analyses provide strong evidence for very high gene flow between populations, but population structure has been retained in several regions, and several populations are nested. Allele frequencies are similar in amphidromous and landlocked forms. Genetic diversity is homogeneous in amphidromous populations. Bayesian admixture analysis infers widespread introgression in Japanese rivers, with a mean introgression proportion of 24 {+/-} 8%. Maximum likelihood admixture graphs detect two migration events from Lake Biwa to anadromous populations. Analyses consistently indicate that hybridisation between translocated landlocked juveniles and native amphidromous Ayu occurs throughout Japanese rivers.
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. Highly efficient synthetic CRISPR RNA/Cas9-based mutagenesis for rapid cardiovascular phenotypic screening in F0 zebrafish The zebrafish is a valuable vertebrate model to study cardiovascular formation and function due to the facile visualization and rapid development of the circulatory system in its externally growing embryos. Despite having distinct advantages, zebrafish have paralogs of many important genes, making reverse genetics approaches inefficient since generating animals bearing multiple gene mutations requires substantial efforts. Here, we present a simple and robust synthetic CRISPR RNA/Cas9-based mutagenesis approach for generating biallelic F0 zebrafish knockouts. Using a dual-guide synthetic CRISPR RNA/Cas9 ribonucleoprotein (dgRNP) system, we compared the efficiency of biallelic gene disruptions following the injections of one, two, and three dgRNPs per gene into the cytoplasm or yolk. We show that simultaneous cytoplasmic injections of three distinct dgRNPs per gene into one-cell stage embryos resulted in the most efficient and consistent biallelic gene disruptions. Importantly, this triple dgRNP approach enables efficient inactivation of cell autonomous and cell non-autonomous gene function, likely due to the low mosaicism of biallelic disruptions. In support of this finding, we provide evidence that the F0 animals generated by this method fully phenocopied the endothelial and peri-vascular defects observed in corresponding stable mutant homozygotes. Moreover, this approach faithfully recapitulated the trunk vessel phenotypes resulting from the genetic interaction between two vegfr2 zebrafish paralogs. Mechanistically, investigation of genome editing and mRNA decay indicates that the combined mutagenic actions of three dgRNPs per gene lead to an increased probability of frameshift mutations, enabling efficient biallelic gene disruptions. Therefore, our approach offers a highly robust genetic platform to quickly assess novel and redundant gene function in F0 zebrafish. Simulation of foraging behavior using a decision-making agent with Bayesian and inverse Bayesian inference: Temporal correlations and power laws in displacement patterns It has been stated that in human migratory behavior, the step length series may have temporal correlation and that there is some relationship between this time dependency and the fact that the frequency distribution of step length follows a power-law distribution. Furthermore, the frequency of occurrence of the step length in some large marine organisms has been found to switch between power-law and exponential distributions, depending on the difficulty of prey acquisition. However, to date it has not been clarified how the aforementioned three phenomena arise: the positive correlation created in the step length series, the relation between the positive correlation of the step length series and the form of an individuals step length distribution, and the switching between power-law and exponential distributions depending on the abundance of prey. This study simulated foraging behavior using the Bayesian decision-making agent simultaneously performing both knowledge learning and knowledge-based inference to analyze how the aforementioned three phenomena arise. In the agent with learning and inference, past experiences were stored as hypotheses (knowledge) and they were used in current foraging behavior; at the same time, the hypothesis continued to be updated based on new experiences. The simulation results show that the agent with both learning and inference has a mechanism that simultaneously causes all the phenomena. Inter-chromosomal linkage disequilibrium and linked fitness cost loci associated with selection for herbicide resistance The adaptation of weedy plants to herbicide is both a significant problem in agriculture and a model for the study of rapid adaptation under regimes of strong selection. Despite recent advances in our understanding of simple genetic changes that lead to resistance, a significant gap remains in our knowledge of resistance controlled by many loci and the evolutionary factors that influence the maintenance of resistance over time. Here, we perform a multi-level analysis involving whole genome sequencing and assembly, resequencing and gene expression analysis to both uncover putative loci involved in nontarget herbicide resistance and to examine evolutionary forces underlying the maintenance of resistance in natural populations. We found loci involved in herbicide detoxification, stress sensing, and alterations in the shikimate acid pathway to be under selection, and confirmed that detoxification is responsible for glyphosate resistance using a functional assay. Furthermore, we found interchromosomal linkage disequilibrium (ILD), most likely associated with epistatic selection, to influence NTSR loci found on separate chromosomes thus potentially mediating resistance through generations. Additionally, by combining the selection screen, differential expression and LD analysis, we identified fitness cost loci that are strongly linked to resistance alleles, indicating the role of genetic hitchhiking in maintaining the cost. Overall, our work strongly suggests that NTSR glyphosate resistance in I. purpurea is conferred by multiple genes which are maintained through generations via ILD, and that the fitness cost associated with resistance in this species is a by-product of genetic-hitchhiking. Monitoring fish communities through environmental DNA metabarcoding in the fish pass system of the second largest hydropower plant in the world The Itaipu Hydroelectric Power Plant is the second largest in the world in power generation. The artificial barrier created by its dam imposes an obstacle for fish migration. Thus, in 2002, a fish pass system, named Piracema Channel, was built to allow fish to access areas upstream of the reservoir. We tested the potential of environmental DNA metabarcoding to monitor the impact of both the dam and associated fish pass system in the Parana River fish communities and to compare it with traditional monitoring methods. Using a fragment of the 12S gene, we characterized richness and community composition based on amplicon sequence variants, operational taxonomic units, and zero-radius OTUs. We combined GenBank and in-house data for taxonomic assignment. We found that different bioinformatics approaches showed similar results. Also, we found a decrease in fish diversity from 2019 to 2020 probably due to the recent extreme drought experienced in southeastern Brazil. The highest alpha diversity was recorded in the mouth of the fish pass system, located in a protected valley with the highest environmental heterogeneity. Despite the clear indication that the reference databases need to be continuously improved, our results demonstrate the analytical efficiency of the metabarcoding to monitor fish species. High-throughput analysis of DNA replication in single human cells reveals the complex nature of replication timing control DNA replication initiates from replication origins firing throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We developed an approach to profile DNA replication from whole-genome sequencing of thousands of single cells, which includes "in silico flow cytometry", a method for discriminating replicating and non-replicating cells. Using two microfluidic platforms, we analyzed up to 2,437 replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that most occur in confined genomic regions. While initiation order is remarkably similar across cells, we unexpectedly identify several subtypes of initiation regions in late-replicating regions. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. RESTORING GLUTAMATE RECEPTOR SIGNALING IN PANCREATIC ALPHA CELLS RESCUES GLUCAGON RESPONSES IN TYPE 1 DIABETES Glucagon secretion from pancreatic alpha cells is crucial to prevent hypoglycemia. For reasons still unknown, people with type 1 diabetes lose this glucoregulatory mechanism and are susceptible to dangerous hypoglycemia. Here we show that alpha cells in living pancreas slices from donors with type 1 diabetes failed to secrete glucagon in response to decreases in glucose concentration, thus mirroring the in vivo unresponsiveness to hypoglycemia. Glucagon content and responses to KCl depolarization were not affected, suggesting that alpha cells retained their secretory potential. By contrast, alpha cells had severely impaired signaling via glutamate receptors of the AMPA/kainate type. Under healthy conditions, activating these receptors was required to elicit full glucagon responses to decreases in glucose levels. In type 1 diabetes, reactivating residual glutamate receptor function with the positive allosteric modulators cyclothiazide and aniracetam restored glucagon secretion in response to hypoglycemia. These positive allosteric modulators are already approved to treat other conditions and could be repurposed to prevent hypoglycemia and improve management of diabetes. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Probing the rules of cell coordination in live tissues by interpretable machine learning based on graph neural networks Robustness in developing and homeostatic tissues is supported by various types of spatiotemporal cell-to-cell interactions. Although live imaging and cell tracking are powerful in providing direct evidence of cell coordination rules, extracting and comparing these rules across many tissues with potentially different length and timescales of coordination requires a versatile framework of analysis. Here we demonstrate that graph neural network (GNN) models are suited for this purpose, by showing how they can be applied to predict cell fate in tissues and utilized to infer the cell interactions governing the multicellular dynamics. Analyzing the live mammalian epidermis data, where spatiotemporal graphs constructed from cell tracks and cell contacts are given as inputs, GNN discovers distinct neighbor cell fate coordination rules that depend on the region of the body. This approach demonstrates how the GNN framework is powerful in inferring general cell interaction rules from live data without prior knowledge of the signaling involved. Relative contributions of sex hormones, sex chromosomes, and gonads to sex differences in tissue gene regulation Sex differences in physiology and disease in mammals result from the effects of three classes of factors that are inherently unequal in males and females: reversible (activational) effects of gonadal hormones, permanent (organizational) effects of gonadal hormones, and cell-autonomous effects of sex chromosomes, as well as genes driven by these classes of factors. Often, these factors act together to cause sex differences in specific phenotypes, but the relative contribution of each and the interactions among them remain unclear. Here, we used the Four Core Genotypes (FCG) mouse model with or without hormone replacement to distinguish the effects of each class of sex-biasing factors on transcriptome regulation in liver and adipose tissues. We found that the activational hormone levels have the strongest influence on gene expression, followed by the organizational gonadal sex effect and, lastly, sex chromosomal effect, along with interactions among the three factors. Tissue specificity was prominent, with a major impact of estradiol on adipose tissue gene regulation, and of testosterone on the liver transcriptome. The networks affected by the three sex-biasing factors include development, immunity and metabolism, and tissue-specific regulators were identified for these networks. Furthermore, the genes affected by individual sex-biasing factors and interactions among factors are associated with human disease traits such as coronary artery disease, diabetes, and inflammatory bowel disease. Our study offers a tissue-specific account of the individual and interactive contributions of major sex-biasing factors to gene regulation that have broad impact on systemic metabolic, endocrine, and immune functions. Transmembrane helices are an overlooked and evolutionarily conserved source of major histocompatibility complex class I and II epitopes Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most haplotypes of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHCI and MHC-II. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells respond more to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs therefore are less prone to escape mutations that enable pathogens to evade T cell responses. Brain deletion of Serf2 shifts amyloid conformation in a mouse model Loss of protein homeostasis accelerates ageing and contributes to age-related brain diseases like Parkinsons and Alzheimers diseases. Disease-specific proteins become aggregation-prone and form amyloid deposits that characterize disease, but the biological mechanisms driving their transition in cells and tissues are poorly understood. MOAG-4/SERF accelerates aggregation and toxicity of disease-related proteins in worm- and cell models. Whether MOAG-4/SERF affects amyloid pathology in brain, however, has remained unknown. Here we show that brain-specific knockout of SERF2 is viable and shifts aggregation in a mouse model for amyloid-beta aggregation. Without affecting amyloid-beta levels, knockout of SERF2 redirects its assembly into structurally and morphologically distinct amyloid fibrils and deposits. Our results show that SERF modifies the structural outcome of the aggregation process in mammalian brain. With structural polymorphisms recently being associated with disease severity and progression, our results may help to understand the molecular basis for structural diversity and find possibilities for structure-based interventions.

HighlightsO_LILoss of SERF2 slows embryonic development and causes perinatal lethality
C_LIO_LISERF2 affects proliferation in a cell-autonomous fashion
C_LIO_LIBrain-specific Serf2 knockout does not affect viability or A{beta} production
C_LIO_LIBrain deletion of Serf2 shifts the amyloid conformation of A{beta}
C_LI Trimeric Photosystem I facilitates energy transfer from phycobilisomes in Synechocystis sp. PCC 6803 In cyanobacteria, phycobilisomes serve as peripheral light-harvesting complexes of the two photosystems, extending their antenna size and the wavelength range of photons available for photosynthesis. The abundance of phycobilisomes, the number of phycobiliproteins they contain, and their light-harvesting function are dynamically adjusted in response to the physiological conditions. Phycobilisomes are also thought to be involved in state transitions that maintain the excitation balance between the two photosystems. Unlike its eukaryotic counterpart, PSI is trimeric in many cyanobacterial species and the physiological significance of this is not well understood. Here we compared the composition and light-harvesting function of phycobilisomes in cells of Synechocystis sp. PCC 6803, which has primarily trimeric PSI, and the{Delta} psaL mutant unable to form trimers. We also investigated a mutant additionally lacking the PsaJ and PsaF subunits of PSI. Both strains with monomeric PSI accumulated significantly more allophycocyanin per chlorophyll, indicating higher abundance of phycobilisomes. On the other hand, a higher phycocyanin:allophycocyanin ratio in WT suggests larger phycobilisomes or the presence of APC-less phycobilisomes (CpcL-type), that are not assembled in cells with monomeric PSI. Steady-state and time-resolved fluorescence spectroscopy at room temperature and 77 K revealed that PSII receives more energy from the phycobilisomes at the expense of PSI in cells with monomeric PSI, regardless of the presence of PsaF. Taken together, these results show that the oligomeric state of PSI has an impact on the excitation energy flow in Synechocystis.

One-sentence summaryCyanobacterial mutants with monomeric PSI show changes in the composition and abundance of phycobilisomes and in the excitation energy transfer to PSII and PSI. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. A somatosensory computation that unifies limbs and tools It is often claimed that tools are embodied by the user, but whether the brain actually repurposes its body-based computations to perform similar tasks with tools is not known. A fundamental computation for localizing touch on the body is trilateration. Here, the location of touch on a limb is computed by integrating estimates of the distance between sensory input and its boundaries (e.g., elbow and wrist of the forearm). As evidence of this computational mechanism, tactile localization on a limb is most precise near its boundaries and lowest in the middle. We show that the brain repurposes trilateration to localize touch on a tool. In a large sample of participants, we found that localizing touch on a tool produced the signature of trilateration, with highest precision close to the base and tip of the tool. A computational model of trilateration provided a good fit to the observed localization behavior. To further demonstrate the computational plausibility of repurposing trilateration, we implemented it in a three-layer neural network that was based on principles of probabilistic population coding. This network determined hit location in tool-centered coordinates by using a tools unique pattern of vibrations when contacting an object. Simulations demonstrated the expected signature of trilateration, in line with the behavioral patterns. Our results have important implications for how trilateration may be implemented by somatosensory neural populations. We conclude that trilateration is a fundamental spatial computation that unifies limbs and tools. De-noising Spatial Transcriptomics Data Based on Position and Image Information Spatial transcriptomics (ST) technologies provide gene expression close to or even superior to single-cell resolution while retaining the physical locations of sequencing and often also providing matched pathology images. However, the expression data captured by ST technologies suffer from high noise levels, as a result of the shallow coverage in each sequencing unit. The extra experimental steps for preserving the spatial locations of sequencing could result in even more severe noises, compared to regular single-cell RNA-sequencing (scRNA-seq). Fortunately, such noises could be largely removed by leveraging information from the physical locations of sequencing, and the tissue and cellular organization reflected by corresponding pathology images. In this work, we demonstrated the extensive levels of noise in ST data. We developed a mathematical model, named Sprod, to impute accurate ST gene expression based on latent space and graph learning of matched location and imaging data. We comprehensively validated Sprod and demonstrated its advantages over prior methods for removing drop-outs in scRNA-seq data. We further showed that, after adequate imputation by Sprod, differential expression analyses, pseudotime analyses, and cell-to-cell interaction inferences yield significantly more informative results. Overall, we envision denoising by Sprod to become a key first step to empower ST technologies for biomedical discoveries and innovations. Use of lignocellulosic residue from 2G ethanol production to enhance methane production through co-digestion This is a pioneer study evaluating the methane (CH4) production potential from residues of integrated 1st (vinasse and filter cake) and 2nd (deacetylation pretreatment liquor from straw) generation (1G2G) sugarcane biorefinery, providing a fully chemical characterization of them and their relation with the anaerobic digestion (AD) process. Small-scale assays provided fundamentals for basing the co-digestion optimization by assessing the optimal co-substrates synergistic conditions. Biochemical Methane Potential (BMP) tests showed co-digestion enhanced CH4 yield of isolated substrates, reaching up to 605 NmLCH4 gVS-1. The association of vinasse and deacetylation liquor as co-substrates increased the BMP by ~38% mostly by nutritionally benefiting the methanogenic activity. The kinetic analysis confirmed that the deacetylation liquor was the co-substrate responsible for improving the CH4 production in the co-digestion systems due to the highest CH4 conversion rate. The alkaline characteristic of the liquor (pH~12) also prevented alkalizing from being added to the co-digestion, an input that normally makes the process economically unfeasible to implement on an industrial scale due to the large quantities required for buffering the reactor. The filter cake had the lowest BMP (262 NmLCH4 gVS-1) and digestibility ([&le;] 40%), further limited by the required stirring to improve the mass transfer of biochemical reactions. The present study drives towards more sustainable use of vinasse, the most voluminous waste from the sugarcane industry, and lignin-rich residues derived from pre-treatment alkaline methods, aiming at an energy-efficient utilization, by at least 16% when compared to the traditional vinasse AD. The experimental and modeling elements from this work indicated the lignin-rich liquor is the main responsible for putting the co-digestion as a disruptive technological arrangement within the 1G2G sugarcane biorefineries, reinforcing the biogas production as the hub of the bioeconomy in the agroindustrial sector. Structural basis for the inhibition of the Bacillus subtilis c-di-AMP cyclase CdaA by the phosphoglucomutase GlmM Cyclic-di-adenosine monophosphate (c-di-AMP) is an important nucleotide signalling molecule, which plays a key role in osmotic regulation in bacteria. Cellular c-di-AMP levels are tightly regulated, as both high and low levels have a negative impact on bacterial growth. Here, we investigated how the activity of the main Bacillus subtilis c-di-AMP cyclase CdaA is regulated by the phosphoglucomutase GlmM. c-di-AMP is produced from two molecules of ATP by proteins containing a deadenylate cyclase (DAC) domain. CdaA is a membrane-linked cyclase with an N-terminal transmembrane domain followed by the cytoplasmic DAC domain. Here we show, using the soluble catalytic B. subtilis CdaACD domain and purified full-length GlmM or the GlmMF369 variant lacking the C-terminal flexible domain 4, that the cyclase and phosphoglucomutase form a stable complex in vitro and that GlmM is a potent cyclase inhibitor. We determined the crystal structure of the individual B. subtilis CdaACD and GlmM proteins, both of which form dimers in the structures, and of the CdaACD:GlmMF369 complex. In the complex structure, a CdaACD dimer is bound to a GlmMF369 dimer in such a manner that GlmM blocks the oligomerization of CdaACD and formation of active head-to-head cyclase oligomers, thus providing molecular details on how GlmM acts as cyclase inhibitor. The function of a key amino acid residue in CdaACD in complex formation was confirmed by mutagenesis analysis. As the amino acids at the CdaACD:GlmM interphase are conserved, we propose that the observed inhibition mechanism of CdaA by GlmM is conserved among Firmicutes. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. 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. H3K36 methylation and DNA-binding both promote Ioc4 recruitment and Isw1b remodeller function The Isw1b chromatin-remodelling complex is specifically recruited to gene bodies to help retain pre-existing histones during transcription by RNA polymerase II. Recruitment is dependent on H3K36 methylation and the Isw1b subunit Ioc4, which contains an N-terminal PWWP domain. Here, we present the crystal structure of the Ioc4-PWWP domain, including a detailed functional characterization of the domain on its own as well as in the context of full-length Ioc4 and the Isw1b remodeller. The Ioc4-PWWP domain preferentially binds H3K36me3-containing nucleosomes. Its ability to bind DNA is required for nucleosome binding. It is also furthered by the unique insertion motif present in Ioc4-PWWP. The ability to bind H3K36me3 and DNA promote the interaction of full-length Ioc4 with nucleosomes in vitro and they are necessary for its recruitment to gene bodies in vivo. Furthermore, a fully functional Ioc4-PWWP domain promotes efficient remodelling by Isw1b and the maintenance of ordered chromatin in vivo, thereby preventing the production of non-coding RNAs. Identifying and correcting multiple sources of misspecification in GWAS summary statistics for polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. A platform of patient-derived microtumors identifies individual treatment responses and therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment.