Optimal strategies for learning multi-ancestry polygenic scores vary across traits Polygenic scores (PGSs) are individual-level measures that aggregate the genome-wide genetic predisposition to a given trait. As PGS have predominantly been developed using European-ancestry samples, trait prediction using such European ancestry-derived PGS is less accurate in non-European ancestry individuals. Although there has been recent progress in combining multiple PGS trained on distinct populations, the problem of how to maximize performance given a multiple-ancestry cohort is largely unexplored. Here, we investigate the effect of sample size and ancestry composition on PGS performance for fifteen traits in UK Biobank. For some traits, PGS estimated using a relatively small African-ancestry training set outperformed, on an African-ancestry test set, PGS estimated using a much larger European-ancestry only training set. We observe similar, but not identical, results when considering other minority-ancestry groups within UK Biobank. Our results emphasise the importance of targeted data collection from underrepresented groups in order to address existing disparities in PGS performance.
Chromosomal rearrangements with stable repertoires of genes and transposable elements in an invasive forest-pathogenic fungus AO_SCPLOWBSTRACTC_SCPLOWChromosomal rearrangements have been largely described among eukaryotes, and may have important consequences on evolution of species. High genome plasticity has been often reported in Fungi, which may explain their apparent ability to quickly adapt to new environments. Cryphonectria parasitica, causing the Chestnut blight disease, is an invasive fungal pathogen species associated with several recent host shifts during its successive introductions from Asia to North America and Europe. Previous cytological karyotyping and genomic studies suggested several chromosomal rearrangements which remains to be described in detail for this species. A serious limitation for valid genome comparisons is the access to robust genome assemblies that usually contain genomic regions of low complexity. We present a new de novo whole-genome assembly obtained from a new method of DNA extraction and long-reads sequencing Nanopore technology obtained from a Japanese isolate sampled in the native area of the species. The comparison with a recently published reference genome revealed stable gene and transposable elements (TEs) repertoires. We also showed that the C. parasitica genome is lowly compartmentalized, with a poor association between TEs and genes, such as those potentially involved in host interactions (i.e., genes coding for small secreted proteins or for secondary metabolites). This genome comparison, however, detected several large chromosomal rearrangements that may have important consequences in gene regulations and sexual mating in this invasive species. This study opens the way for more comparisons of high-quality assembled genomes, and questions the role of structural variations in the invasive success of this fungal pathogen species.
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. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Cellular Fitness Phenotypes of Cancer Target Genes from Oncobiology to Cancer Therapeutics To define the growing significance of cellular targets of cancer drugs, we examined the fitness dependency of cellular targets or effectors of cancer drug targets across human cancer cells from 19 cancer types. We observed that the deletion of 35 out of 47 cellular mediators or targets of oncology drugs did not result in the expected loss of cell fitness in appropriate cancer types for which drugs targeting or utilizing these molecules were approved. Additionally, our analysis recognized 43 cellular targets as fitness genes in several cancer types in which these drugs were not approved, and thus, providing clues repurposing approved oncology drugs in cancer types. For example, we found the widespread upregulation and fitness dependency of the components of the mevalonate and purine biosynthesis pathways (currently targeted by bisphosphonates, statins, and pemetrexed in certain cancers) and an association between the overexpression of these targets and reduction in the overall survival duration of patients with breast and other hard-to-treat cancers, for which such drugs are not approved. In brief, the present analysis raised cautions about off-target and undesirable effects of certain oncology drugs in a subset of cancers where the in-tended cellular effectors of drug might not be fitness genes and offers a potential rationale for repurposing certain approved oncology drugs for targeted therapeutics in additional cancer types. Host cell membrane capture by the SARS CoV-2 spike protein fusion intermediate Cell entry by SARS-CoV-2 is accomplished by the S2 subunit of the spike S protein on the virion surface by capture of the host cell membrane and fusion with the viral envelope. Capture and fusion require the prefusion S2 to transit to its potent, fusogenic form, the fusion intermediate (FI). However, the FI structure is unknown, detailed computational models of the FI are unavailable, and the mechanisms and timing of membrane capture and fusion are not established. Here, we constructed a full-length model of the CoV-2 FI by extrapolating from known CoV-2 pre- and postfusion structures. In atomistic and coarse-grained molecular dynamics simulations the FI was remarkably flexible and executed large bending and extensional fluctuations due to three hinges in the C-terminal base. Simulations suggested a host cell membrane capture time of [~] 2 ms. Isolated fusion peptide simulations identified an N-terminal helix that directed and maintained binding to the membrane but grossly underestimated the binding time, showing that the fusion peptide environment is radically altered when attached to its host fusion protein. The large configurational fluctuations of the FI generated a substantial exploration volume that aided capture of the target membrane, and may set the waiting time for fluctuation-triggered refolding of the FI that draws the viral envelope and host cell membrane together for fusion. These results describe the FI as a machinery designed for efficient membrane capture and suggest novel potential drug targets.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=117 SRC="FIGDIR/small/439051v3_ufig1.gif" ALT="Figure 1">
View larger version (29K):
org.highwire.dtl.DTLVardef@1374b59org.highwire.dtl.DTLVardef@133d21forg.highwire.dtl.DTLVardef@127f2ddorg.highwire.dtl.DTLVardef@1b9ab8f_HPS_FORMAT_FIGEXP  M_FIG C_FIG Venomics and antivenomics of Indian spectacled cobra (Naja naja) from the Western Ghats Venom proteome profiling of Naja naja from the Western Ghats region in Kerala was achieved through SDS-PAGE and RP-HPLC followed by Q-TOF LC-MS/MS analysis, incorporating PEAKS and Novor assisted de novo sequencing methodologies. A total of 115 proteins distributed across 17 different enzymatic and non-enzymatic venom protein families were identified through conventional and 39 peptides through homology-driven proteomics approaches. Fourteen peptides derived through de novo complements the Mascot data indicating the importance of homology-driven approaches in improving protein sequence information. Among the protein families identified, glutathione peroxidase and endonuclease were reported for the first time in the Indian cobra venom. Immunological cross-reactivity assessed using Indian polyvalent antivenoms suggested that VINS showed better EC50 (2.48 g/mL) value than that of PSAV (6.04 g/mL) and Virchow (6.03 g/mL) antivenoms. Western blotting experiments indicated that all the antivenoms elicited poor binding specificities, especially towards low molecular mass proteins. Second-generation antivenomics studies revealed that VINS antivenom was less efficient to detect many low molecular mass proteins such as three-finger toxins and Kunitz-type serine protease Inhibitors. Taken together, the present study enabled a large-scale characterization of the venom proteome of Naja naja from the Western Ghats and emphasized the need for developing more efficient antivenoms.

HighlightsO_LIProteomics of cobra venom resulted in the identification of 115 proteins representing 17 snake venom protein families.
C_LIO_LIDe novo approaches exclusively yielded 39 peptides harbouring multiple amino acid mutations.
C_LIO_LIGlutathione peroxidase and endonuclease were identified for the first time in Indian cobra venom.
C_LIO_LIIndian polyvalent antivenoms showed varying cross-reactivity towards cobra venom.
C_LIO_LIVINS antivenom was less efficient to detect many low molecular mass proteins (< 20 kDa).
C_LI XAV-19, a swine glyco-humanized polyclonal antibody against SARS-CoV-2 Spike receptor-binding domain, targets multiple epitopes and broadly neutralizes variants. Amino acid substitutions and deletions in Spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants can reduce the effectiveness of monoclonal antibodies (mAbs). In contrast, heterologous polyclonal antibodies raised against S protein, through the recognition of multiple target epitopes, have the potential to maintain neutralization capacities. XAV-19 is a swine glyco-humanized polyclonal neutralizing antibody raised against the receptor binding domain (RBD) of the Wuhan-Hu-1 Spike protein of SARS-CoV-2. XAV-19 target epitopes were found distributed all over the RBD and particularly cover the receptor binding motives (RBM), in direct contact sites with the Angiotensin Converting Enzyme-2 (ACE-2). Therefore, in Spike/ACE2 interaction assays, XAV-19 showed potent neutralization capacities of the original Wuhan Spike and of the United Kingdom (Alpha/B.1.1.7) and South African (Beta/B.1.351) variants. These results were confirmed by cytopathogenic assays using Vero E6 and live virus variants including the Brazil (Gamma/P.1) and the Indian (Delta/B.1.617.2) variants. In a selective pressure study with the Beta strain on Vero E6 cells conducted over 1 month, no mutation was associated with addition of increasing doses XAV-19. The potential to reduce viral load in lungs was confirmed in a human ACE2 transduced mouse model. XAV-19 is currently evaluated in patients hospitalized for COVID-19-induced moderate pneumonia in a phase 2a-2b (NCT04453384) where safety was already demonstrated and in an ongoing 2/3 trial (NCT04928430) to evaluate the efficacy and safety of XAV-19 in patients with moderate-to-severe COVID-19. Owing to its polyclonal nature and its glyco-humanization, XAV-19 may provide a novel safe and effective therapeutic tool to mitigate the severity of coronavirus disease 2019 (Covid-19) including the different variants of concern identified so far. Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we developed an inference index (C-score) for Cell-specific Expression- BUffering (CEBU), whereby functional buffering is mediated via expression of buffering genes in specific cells and tissues in humans. By computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq, we report that C- score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results reveal CEBU as a buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

Summary blurbWe unveil a genome-wide functional buffering mechanism, termed Cell-specific Expression Buffering (CEBU), whereby gene expression contributes to functional buffering in specific cell types and tissues. We link CEBU to genetic interactions, tissue homeostasis and cancer robustness. Structural and molecular determinants for the interaction of ExbB from Serratia marcescens and HasB, a TonB paralog ExbB and ExbD are cytoplasmic membrane proteins that associate with TonB to convey the energy of the proton-motive force to outer membrane receptors in Gram-negative bacteria for iron uptake. The opportunistic pathogen Serratia marcescens (Sm) possesses both TonB and a heme-specific TonB paralog, HasB. ExbBSm has a long periplasmic extension absent in other bacteria such as E. coli (Ec). Long ExbBs are found in several genera of Alphaproteobacteria, most often in correlation with a hasB gene. We investigated specificity determinants of ExbBSm and HasB. We determined the cryo-EM structures of ExbBSm and of the ExbB-ExbDSm complex from S. marcescens. ExbBSm alone is a stable pentamer, and its complex includes two ExbD monomers. We showed that ExbBSm extension interacts with HasB and is involved in heme acquisition and we identified key residues in the membrane domain of ExbBSm and ExbBEc, essential for function and likely involved in the interaction with TonB/HasB. Our results shed light on the new class of inner membrane energy machinery formed by ExbB,ExbD and HasB. Interacting phenotypes and the coevolutionary process: Interspecific indirect genetic effects alter coevolutionary dynamics Coevolution occurs when species interact to influence one anothers fitness, resulting in reciprocal evolutionary change. In many coevolving lineages, trait expression in one species is modified by the genotypes and phenotypes of the other, forming feedback loops reminiscent of models of intraspecific social evolution. Here, we adapt the theory of within-species social evolution, characterized by indirect genetic effects and social selection imposed by interacting individuals, to the case of interspecific interactions. In a trait-based model, we derive general expressions for multivariate evolutionary change in two species and the expected between-species covariance in evolutionary change when selection varies across space. We show that reciprocal interspecific indirect genetic effects can dominate the coevolutionary process and drive patterns of correlated evolution beyond what is expected from direct selection alone. In extreme cases, interspecific indirect genetic effects can lead to coevolution when selection does not covary between species or even when one species lacks genetic variance. Moreover, our model indicates that interspecific indirect genetic effects may interact in complex ways with cross-species selection to determine the course of coevolution. Importantly, our model makes empirically testable predictions for how different forms of reciprocal interactions contribute to the coevolutionary process. Novel MOG analogues to explore the MCT2 pharmacophore, α-ketoglutarate biology and cellular effects of N-oxalylglycine -ketoglutarate (KG) is a central metabolic node with a broad influence on cellular physiology. The KG analogue N-oxalylglycine (NOG) and its membrane-permeable pro-drug derivative dimethyl-oxalylglycine (DMOG) have been extensively used as tools to study prolyl hydroxylases (PHDs) and other KG-dependent processes. In cell culture media, DMOG is rapidly converted to MOG, which enters cells through monocarboxylate transporter MCT2, leading to intracellular NOG concentrations that are sufficiently high to inhibit glutaminolysis enzymes and cause cytotoxicity. Therefore, the degree of (D)MOG instability together with MCT2 expression levels determine the intracellular targets NOG engages with and, ultimately, its effects on cell viability. Here we designed and characterised a series of MOG analogues with the aims of improving stability and exploring the functional requirements for interaction with MCT2, a relatively understudied member of the SLC16 family. We report MOG analogues that maintain ability to enter cells via MCT2, and identify compounds that do not inhibit glutaminolysis or cause cytotoxicity but can still inhibit PHDs. We use these analogues to show that glutaminolysis-induced activation of mTORC1 can be uncoupled from PHD activity. Therefore, these new compounds can help deconvolute cellular effects that result from the polypharmacological action of NOG. Protective host-dependent antagonism among 	Pseudomonas in the Arabidopsis phyllosphere The plant microbiome is a rich biotic environment, comprising numerous taxa. The community structure of these colonizers is constrained by multiple factors, including host-microbe and microbe-microbe interactions, as well as the interplay between the two. While much can be learned from pairwise relationships between individual hosts and microbes, or individual microbes with themselves, the ensemble of interrelations between the host and microbial consortia may lead to different outcomes that are not easily predicted from the individual interactions. Their study can thus provide new insights into the complex relationship between plants and microbes. Of particular importance is how strain-specific such plant-microbe-microbe interactions are, and how they eventually affect plant health. Here, we test strain-level interactions in the phyllosphere between groups of co-existing commensal and pathogenic Pseudomonas among each other and with A. thaliana, by employing synthetic communities of genome-barcoded isolates. We found that commensal Pseudomonas prompted a host response leading to a selective inhibition of a specific pathogenic lineage, resulting in plant protection. The extent of plant protection, however, was dependent on plant genotype, indicating that these effects are host-mediated. There were similar genotype-specific effects on the microbe side, as we could pinpoint an individual Pseudomonas isolate as the predominant cause for this differential interaction. Collectively, our work highlights how within-species genetic differences on both the host and microbe side can have profound effects on host-microbe-microbe dynamics. The paradigm that we have established provides a platform for the study of host-dependent microbe-microbe competition and cooperation in the A. thaliana-Pseudomonas system. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Manduca sexta experience high parasitoid pressures in the field but minor fitness costs of consuming plant secondary compounds O_LIPlant-herbivore co-evolutionary interactions have led to a range of plant defenses that minimize insect damage and a suite of counter-adaptations that allow herbivores to feed on defended plants. Consuming plant secondary compounds results in herbivore growth and developmental costs but can have beneficial effects such as deterrence or harm of parasitoid enemies. Therefore, the role of secondary compounds on herbivore fitness must be considered in the context of the abundance and level of harm from natural enemies and the costs herbivores incur feeding on plant secondary compounds.
C_LIO_LIIn this study, I combined field measurements of Cotesia congregata wasp parasitism pressure with detailed measurements of the costs of plant secondary compounds across developmental stages in the herbivore host, Manduca sexta.
C_LIO_LII show that C. congregata parasitoids exert large negative selective pressures, killing 31-57% of M. sexta larvae in the field. Manduca sexta developed fastest during instars most at risk for parasitoid oviposition but growth was slowed by consumption of plant secondary compounds. The negative effects of consuming plant secondary compounds as larvae influenced adult size traits but there were no immune, survival, or fecundity costs.
C_LIO_LIThese results suggest that developmental costs experienced by M. sexta herbivores consuming defensive compounds are minor in comparison to the strong negative survival pressures from abundant parasitoid enemies.
C_LI Prediction of infectivity of SARS-CoV2: Mathematical model with analysis of docking simulation for spike proteins and angiotensin-converting enzyme 2 Variants of a coronavirus (SARS-CoV-2) have been spreading in a global pandemic. Improved understanding of the infectivity of future new variants is important so that effective countermeasures against them can be quickly undertaken. In our research reported here, we aimed to predict the infectivity of SARS-CoV-2 by using a mathematical model with molecular simulation analysis, and we used phylogenetic analysis to determine the evolutionary distance of the spike protein gene (S gene) of SARS-CoV-2. We subjected the six variants and the wild type of spike protein and human angiotensin-converting enzyme 2 (ACE2) to molecular docking simulation analyses to understand the binding affinity of spike protein and ACE2. We then utilized regression analysis of the correlation coefficient of the mathematical model and the infectivity of SARS-CoV-2 to predict infectivity. The evolutionary distance of the S gene correlated with the infectivity of SARS-CoV-2 variants. The coefficient of the mathematical model obtained with results of molecular docking simulation also correlated with the infectivity of SARS-CoV-2 variants. These results suggest that the data from the docking simulation for the receptor binding domain of variant spike proteins and human ACE2 were valuable for prediction of SARS-CoV-2 infectivity. In addition, we developed a mathematical model for prediction of SARS-CoV-2 variant infectivity by using binding affinity obtained via molecular docking and the evolutionary distance of the S gene. Astrocyte-secreted glypican-4 drives APOE4-dependent tau pathology The aggregation of tau proteins into insoluble filaments and the spread of these filaments across brain regions are the major drivers of neurodegeneration in tauopathies, including in Alzheimers disease (AD). Apolipoprotein E4 (APOE4), a crucial genetic risk factor for late-onset AD, has been shown to exacerbate tau pathology in mouse models. However, the exact mechanisms through which APOE4 induces tau pathology remains unknown. Here, we report that the astrocyte-secreted protein glypican-4 (GPC-4), which we identify as a novel binding partner of APOE4, drives tau pathology. We discovered that GPC-4 preferentially interacts with APOE4 in comparison to APOE2, considered to be a protective allele to AD, and post-mortem APOE4-carrying AD brains highly express GPC-4 in neurotoxic astrocytes. The astrocyte-secreted GPC-4 induced both tau accumulation and propagation in vitro. CRISPR/dCas9 mediated activation of GPC-4 in a tauopathy mouse model robustly induced tau pathology. Using in vitro Tau FRET-biosensor cells, human iPSCs-derived astrocytes and an in vivo mouse model, we found that APOE4-induced tau pathology was greatly diminished in the absence of GPC-4. We further show that APOE4-mediated surface trafficking of APOE receptor low-density lipoprotein receptor-related protein 1 (LRP1) through GPC-4 can be a gateway to tau spreading. Together, our data comprehensively demonstrate that APOE4-induced tau pathologies are directly mediated by GPC-4. Loss of SUMO-specific protease 2 causes isolated glucocorticoid deficiency by blocking adrenal cortex zonal transdifferentiation SUMOylation is a dynamic posttranslational modification, that provides fine-tuning of protein function involved in the cellular response to stress, differentiation, and tissue development. In the adrenal cortex, an emblematic endocrine organ that mediates adaptation to physiological demands, the SUMOylation gradient is inversely correlated with the gradient of cellular differentiation raising important questions about its role in functional zonation and the response to stress. Considering that SUMO-specific protease 2 (SENP2), a deSUMOylating enzyme, is upregulated by ACTH/PKA signalling within the zona Fasciculata (zF), we generated mice with adrenal-specific Senp2 loss to address these questions. Disruption of SENP2 activity in steroidogenic cells leads to specific hypoplasia of the zF, a blunted reponses to ACTH and isolated glucocorticoid deficiency. Mechanistically, overSUMOylation resulting from SENP2 loss shifts the balance between ACTH/PKA and WNT/{beta}-catenin signalling leading to repression of PKA activity and ectopic activation of {beta}-catenin. At the cellular level, this blocks transdifferentiation of {beta}-catenin-positive zona Glomerulosa cells into zF cells and sensitises them to premature apoptosis. Our findings indicate that the SUMO pathway is critical for adrenal homeostasis and stress responsiveness. Cas9-mediated gene disruption in tetraploid Giardia intestinalis CRISPR/Cas9 system is an extremely powerful technique that is extensively used for various genome modifications in different organisms including parasitic protists. Giardia intestinalis, a protist parasite infecting about 280 million people around the world each year, has been eluding the routine use of CRISPR/Cas9 for generating knock-out cell lines due to the presence of four copies of each gene in its two nuclei. Apart from single exception employing rather laborious Cre/loxP system, no full knock-out cell line has been established yet. In this work, we show the ability of in-vitro assembled CRISPR/Cas9 components to successfully edit the genome of G. intestinalis. We further established a cell line stably expressing Cas9 in both G. intestinalis nuclei. Subsequent introduction of a template for homologous recombination containing the transcription units for the resistance marker and gRNA resulted in the removal of all gene copies at once for three independent experimental genes, mem, cwp1 and mlf1. The method was also applicable for the incomplete disruption of an essential gene, as documented by markedly decreased expression of tom40. Finally, testing the efficiency of Cas9-induced recombination revealed that homologous arms as short as 150 bp can be sufficient to establish a full knock-out cell line in G. intestinalis. 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. 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. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Liver ductal organoids reconstruct intrahepatic biliary trees in decellularized liver grafts Three-dimensional scaffolds decellularized from native organs are a promising technique to establish engineered liver grafts and overcome the current shortage of donor organs. However, limited sources of bile duct cells and inappropriate cell distribution in bioengineered liver grafts have hindered their practical application. Organoid technology is anticipated to be an excellent tool for the advancement of regenerative medicine. In the present study, we reconstructed intrahepatic bile ducts in a rat decellularized liver graft by recellularization with liver ductal organoids. Using an ex vivo perfusion culture system, we demonstrated the biliary characteristics of repopulated mouse liver organoids, which maintained bile duct markers and reconstructed biliary tree-like networks with luminal structures. We also established a method for the co-recellularization with engineered bile ducts and primary hepatocytes, revealing the appropriate cell distribution to mimic the native liver. We then utilized this model in human organoids to demonstrate the reconstructed bile ducts. Our results show that liver ductal organoids are a potential cell source for bile ducts from bioengineered liver grafts using three-dimensional scaffolds. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development.