@article{CaliskanCaliskanRasbachetal.2023, author = {Caliskan, Aylin and Caliskan, Deniz and Rasbach, Lauritz and Yu, Weimeng and Dandekar, Thomas and Breitenbach, Tim}, title = {Optimized cell type signatures revealed from single-cell data by combining principal feature analysis, mutual information, and machine learning}, series = {Computational and Structural Biotechnology Journal}, volume = {21}, journal = {Computational and Structural Biotechnology Journal}, issn = {2001-0370}, doi = {10.1016/j.csbj.2023.06.002}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-349989}, pages = {3293-3314}, year = {2023}, abstract = {Machine learning techniques are excellent to analyze expression data from single cells. These techniques impact all fields ranging from cell annotation and clustering to signature identification. The presented framework evaluates gene selection sets how far they optimally separate defined phenotypes or cell groups. This innovation overcomes the present limitation to objectively and correctly identify a small gene set of high information content regarding separating phenotypes for which corresponding code scripts are provided. The small but meaningful subset of the original genes (or feature space) facilitates human interpretability of the differences of the phenotypes including those found by machine learning results and may even turn correlations between genes and phenotypes into a causal explanation. For the feature selection task, the principal feature analysis is utilized which reduces redundant information while selecting genes that carry the information for separating the phenotypes. In this context, the presented framework shows explainability of unsupervised learning as it reveals cell-type specific signatures. Apart from a Seurat preprocessing tool and the PFA script, the pipeline uses mutual information to balance accuracy and size of the gene set if desired. A validation part to evaluate the gene selection for their information content regarding the separation of the phenotypes is provided as well, binary and multiclass classification of 3 or 4 groups are studied. Results from different single-cell data are presented. In each, only about ten out of more than 30000 genes are identified as carrying the relevant information. The code is provided in a GitHub repository at https://github.com/AC-PHD/Seurat_PFA_pipeline.}, language = {en} } @article{EngstlerBeneke2023, author = {Engstler, Markus and Beneke, Tom}, title = {Gene editing and scalable functional genomic screening in Leishmania species using the CRISPR/Cas9 cytosine base editor toolbox LeishBASEedit}, series = {eLife}, volume = {12}, journal = {eLife}, doi = {10.7554/eLife.85605}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-350002}, year = {2023}, abstract = {CRISPR/Cas9 gene editing has revolutionised loss-of-function experiments in Leishmania, the causative agent of leishmaniasis. As Leishmania lack a functional non-homologous DNA end joining pathway however, obtaining null mutants typically requires additional donor DNA, selection of drug resistance-associated edits or time-consuming isolation of clones. Genome-wide loss-of-function screens across different conditions and across multiple Leishmania species are therefore unfeasible at present. Here, we report a CRISPR/Cas9 cytosine base editor (CBE) toolbox that overcomes these limitations. We employed CBEs in Leishmania to introduce STOP codons by converting cytosine into thymine and created http://www.leishbaseedit.net/ for CBE primer design in kinetoplastids. Through reporter assays and by targeting single- and multi-copy genes in L. mexicana, L. major, L. donovani, and L. infantum, we demonstrate how this tool can efficiently generate functional null mutants by expressing just one single-guide RNA, reaching up to 100\% editing rate in non-clonal populations. We then generated a Leishmania-optimised CBE and successfully targeted an essential gene in a plasmid library delivered loss-of-function screen in L. mexicana. Since our method does not require DNA double-strand breaks, homologous recombination, donor DNA, or isolation of clones, we believe that this enables for the first time functional genetic screens in Leishmania via delivery of plasmid libraries.}, language = {en} } @article{SalihogluSrivastavaLiangetal.2023, author = {Salihoglu, Rana and Srivastava, Mugdha and Liang, Chunguang and Schilling, Klaus and Szalay, Aladar and Bencurova, Elena and Dandekar, Thomas}, title = {PRO-Simat: Protein network simulation and design tool}, series = {Computational and Structural Biotechnology Journal}, volume = {21}, journal = {Computational and Structural Biotechnology Journal}, issn = {2001-0370}, doi = {10.1016/j.csbj.2023.04.023}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-350034}, pages = {2767-2779}, year = {2023}, abstract = {PRO-Simat is a simulation tool for analysing protein interaction networks, their dynamic change and pathway engineering. It provides GO enrichment, KEGG pathway analyses, and network visualisation from an integrated database of more than 8 million protein-protein interactions across 32 model organisms and the human proteome. We integrated dynamical network simulation using the Jimena framework, which quickly and efficiently simulates Boolean genetic regulatory networks. It enables simulation outputs with in-depth analysis of the type, strength, duration and pathway of the protein interactions on the website. Furthermore, the user can efficiently edit and analyse the effect of network modifications and engineering experiments. In case studies, applications of PRO-Simat are demonstrated: (i) understanding mutually exclusive differentiation pathways in Bacillus subtilis, (ii) making Vaccinia virus oncolytic by switching on its viral replication mainly in cancer cells and triggering cancer cell apoptosis and (iii) optogenetic control of nucleotide processing protein networks to operate DNA storage. Multilevel communication between components is critical for efficient network switching, as demonstrated by a general census on prokaryotic and eukaryotic networks and comparing design with synthetic networks using PRO-Simat. The tool is available at https://prosimat.heinzelab.de/ as a web-based query server.}, language = {en} } @article{ConradKehlMuelleretal.2023, author = {Conrad, David and Kehl, Alexandra and M{\"u}ller, Tobias and Klopfleisch, Robert and Aupperle-Lellbach, Heike}, title = {Immunohistochemical and molecular genetic analysis of canine digital mast cell tumours}, series = {Animals}, volume = {13}, journal = {Animals}, number = {10}, issn = {2076-2615}, doi = {10.3390/ani13101694}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-319199}, year = {2023}, abstract = {Grading, immunohistochemistry and c-kit mutation status are criteria for assessing the prognosis and therapeutic options of canine cutaneous mast cell tumours (MCTs). As a subset, canine digital MCTs have rarely been explored in this context. Therefore, in this retrospective study, 68 paraffin-embedded canine digital MCTs were analysed, and histological grading was assessed according to Patnaik and Kiupel. The immunohistochemical markers KIT and Ki67 were used, as well as polymerase chain reaction (PCR) for mutational screening in c-kit exons 8, 9, 11 and 14. Patnaik grading resulted in 22.1\% grade I, 67.6\% grade II and 10.3\% grade III tumours. Some 86.8\% of the digital MCTs were Kiupel low-grade. Aberrant KIT staining patterns II and III were found in 58.8\%, and a count of more than 23 Ki67-positive cells in 52.3\% of the cases. Both parameters were significantly associated with an internal tandem duplication (ITD) in c-kit exon 11 (12.7\%). French Bulldogs, which tend to form well-differentiated cutaneous MCTs, had a higher proportion of digital high-grade MCTs and ITD in c-kit exon 11 compared with mongrels. Due to its retrospective nature, this study did not allow for an analysis of survival data. Nevertheless, it may contribute to the targeted characterisation of digital MCTs.}, language = {en} } @phdthesis{Nirchal2024, author = {Nirchal, Naveen Kumar}, title = {Mechanistische Regulierung des gastro{\"o}sophagealen {\"U}bergangs und die Rolle der Retins{\"a}ure bei der Entwicklung des Barrett-{\"O}sophagus}, doi = {10.25972/OPUS-31155}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-311556}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2024}, abstract = {Der gastro{\"o}sophageale {\"U}bergang (GEJ), der die Region abgrenzt, in der der distale {\"O}sophagus auf die proximale Magenregion trifft, ist bekannt f{\"u}r die Entwicklung pathologischer Zust{\"a}nde, wie Metaplasie und Adenokarzinom des {\"O}sophagus (EAC). Es ist wichtig, die Mechanismen der Entwicklungsstadien zu verstehen, die zu EAC f{\"u}hren, da die Inzidenzrate von EAC in den letzten 4 Jahrzehnten um das 7-fache gestiegen ist und die Gesamt{\"u}berlebensrate von 5 Jahren 18,4 \% betr{\"a}gt. In den meisten F{\"a}llenwird die Diagnose im fortgeschrittenen Stadium ohne vorherige Symptome erstellt. Der Hauptvorl{\"a}ufer f{\"u}r die Entwicklung von EAC ist eine pr{\"a}maligne Vorstufe namens Barrett-{\"O}sophagus (BE). BE ist der metaplastische Zustand, bei dem das mehrschichtige Plattenepithel des nativen {\"O}sophagus durch ein spezialisiertes einschichtiges S{\"a}ulenepithel ersetzt wird, das die molekularen Eigenschaften des Magen- sowie des Darmepithels aufweist. Zu den wichtigsten Risikofaktoren f{\"u}r die Entwicklung von BE geh{\"o}ren die chronische gastro{\"o}sophageale Refluxkrankheit (GERD), eine ver{\"a}nderte Mikrobiota und ver{\"a}nderte Retins{\"a}ure-Signalwege (RA). Es ist unklar, welche Zelle der Ursprung f{\"u}r BE ist, da es keine eindeutigen Beweisen f{\"u}r den Prozess der BE-Initiation gibt. In dieser Arbeit habe ich untersucht, wie die GEJ-Hom{\"o}ostase in gesundem Gewebe durch stammzellregulatorische Morphogene aufrechterhalten wird, welche Rolle der Vitamin-A (RA-Signal{\"u}bertragung) spieltund wie ihre Ver{\"a}nderung zur BE-Entwicklung beitr{\"a}gt. Im ersten Teil meiner Dissertation habe ich anhand von Einzelmolek{\"u}l-RNA in situ-Hybridisierung und Immunhistochemie eindeutig das Vorhandensein von zwei Arten von Epithelzellen nachweisen k{\"o}nnen, dem Plattenepithel in der Speiser{\"o}hre und dem S{\"a}ulenepithel imMagenbereich des GEJ. Mittels Abstammungsanalysen im Mausmodell konnte ich zeigen, dass die Epithelzellen des {\"O}sophagus und des Magens von zwei verschiedenen epithelialen Stammzelllinien imGEJ abstammen. Die Grenze zwischen Plattenepithel und S{\"a}ulenepithelzellen im SCJ des GEJ wirddurch gegens{\"a}tzliche Wnt-Mikroumgebungen streng reguliert. Plattenepithelstammzellen des {\"O}sophagus werden durch das Wnt-hemmende Mikroumgebungssignal aufrechterhalten, w{\"a}hrend Magens{\"a}ulenepithelzellen durch das Wnt-aktivierende Signal aus dem Stromakompartiment erhalten werden. Ich habe die in vivo Erhaltung der Epithelstammzellen des GEJ mit Hilfe eines in vitro Epithel-3D-Organoidkulturmodells rekonstruiert. Das Wachstum und die Vermehrung von Magens{\"a}ulenepithel-Organoiden h{\"a}ngen von Wnt-Wachstumsfaktoren ab, w{\"a}hrend das Wachstum von Plattenepithel-Organoiden von Wnt-defizienten Kulturbedingungen abh{\"a}ngt. Dar{\"u}ber hinaus zeigte die Einzelzell-RNA-Sequenzanalyse (scRNA-seq) der aus Organoiden gewonnenenEpithelzellen, dass der nicht-kanonische Wnt/ planar cell polarity (PCP) Signalweg an der Regulierung der Plattenepithelzellen beteiligt ist. Im Gegensatz dazu werden s{\"a}ulenf{\"o}rmige Magenepithelzellen durch den kanonischen Wnt/beta-Catenin- und den nicht-kanonischen Wnt/Ca2+-Weg reguliert. Meine Daten zeigen, dass die SCJ-Epithelzellen, die am GEJ verschmelzen, durch entgegengesetzte stromale Wnt-Faktoren und unterschiedliche Wnt-Weg-Signalee in den Epithelzellen reguliert werden. Im zweiten Teil der Dissertation untersuchte ich die Rolle der bioaktiven Vitamin A Verbindung RA auf {\"O}sophagus- und Magenepithelstammzellen. Die In-vitro-Behandlung von epithelialen Organoiden der Speiser{\"o}hre und des Magens mitRA oder seinem pharmakologischen Inhibitors BMS 493 zeigte, dass jeder Zelltyp unterschiedlich reguliert wurde. Ich beobachtete, dass eine verst{\"a}rkte RA die Differenzierung von Stammzellen und den Verlust der Schichtung f{\"o}rderte, w{\"a}hrend die RA-Hemmung zu einer verst{\"a}rkten Stammzellbildung und Regeneration im mehrschichtigen Epithel der Speiser{\"o}hre f{\"u}hrte. Im Gegensatz zur Speiser{\"o}hre ist der RA-Signalweg in Magen-Organoiden aktiv, und die Hemmung von RA hat ein reduziertes Wachstum von Magen-Organoiden. Globale transkriptomische Daten und scRNA-seq-Daten zeigten, dass derRA-Signalweg einen Ruheph{\"a}notyp in den {\"O}sophaguszellen induziert. Dagegen f{\"u}hrt das Fehlen von RA in Magenepithelzellen zur Expression von Genen, die mit BE assoziiert sind. Daher isteine r{\"a}umlich definierte Regulation der Wnt- und Retins{\"a}ure-Signalgebung amGEJ entscheidend f{\"u}r eine gesunde Hom{\"o}ostase, und ihre St{\"o}rung f{\"u}hrt zur Entwicklung von Krankheiten.}, subject = {Retinoes{\"a}ure}, language = {en} } @phdthesis{Weisert2024, author = {Weisert, Nadine}, title = {Characterization of telomere-associated proteins in \(Trypanosoma\) \(brucei\)}, doi = {10.25972/OPUS-35273}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-352732}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2024}, abstract = {The unicellular pathogen Trypanosoma brucei is the causative agent of African trypanosomiasis, an endemic disease prevalent in sub-Saharan Africa. Trypanosoma brucei alternates between a mammalian host and the tsetse fly vector. The extracellular parasite survives in the mammalian bloodstream by periodically exchanging their ˈvariant surface glycoproteinˈ (VSG) coat to evade the host immune response. This antigenic variation is achieved through monoallelic expression of one VSG variant from subtelomeric ˈbloodstream form expression sitesˈ (BES) at a given timepoint. During the differentiation from the bloodstream form (BSF) to the procyclic form (PCF) in the tsetse fly midgut, the stage specific surface protein is transcriptionally silenced and replaced by procyclins. Due to their subtelomeric localization on the chromosomes, VSG transcription and silencing is partly regulated by homologues of the mammalian telomere complex such as TbTRF, TbTIF2 and TbRAP1 as well as by ˈtelomere-associated proteinsˈ (TelAPs) like TelAP1. To gain more insights into transcription regulation of VSG genes, the identification and characterization of other TelAPs is critical and has not yet been achieved. In a previous study, two biochemical approaches were used to identify other novel TelAPs. By using ˈco-immunoprecipitationˈ (co-IP) to enrich possible interaction partners of TbTRF and by affinity chromatography using telomeric repeat oligonucleotides, a listing of TelAP candidates has been conducted. With this approach TelAP1 was identified as a novel component of the telomere complex, involved in the kinetics of transcriptional BES silencing during BSF to PCF differentiation. To gain further insights into the telomere complex composition, other previously enriched proteins were characterized through a screening process using RNA interference to deplete potential candidates. VSG expression profile changes and overall proteomic changes after depletion were analyzed by mass spectrometry. With this method, one can gain insights into the functions of the proteins and their involvement in VSG expression site regulation. To validate the interaction of proteins enriched by co-IP with TbTRF and TelAP1 and to identify novel interaction proteins, I performed reciprocal affinity purifications of the four most promising candidates (TelAP2, TelAP3, PPL2 and PolIE) and additionally confirmed colocalization of two candidates with TbTRF via immunofluorescence (TelAP2, TelAP3). TelAP3 colocalizes with TbTRF and potentially interacts with TbTRF, TbTIF2, TelAP1 and TelAP2, as well as with two translesion polymerases PPL2 and PolIE in BSF. PPL2 and PolIE seem to be in close contact to each other at the telomeric ends and fulfill different roles as only PolIE is involved in VSG regulation while PPL2 is not. TelAP2 was previously characterized to be associated with telomeres by partially colocalizing with TbTRF and cells show a VSG derepression phenotype when the protein was depleted. Here I show that TelAP2 interacts with the telomere-binding proteins TbTRF and TbTIF2 as well as with the telomere-associated protein TelAP1 in BSF and that TelAP2 depletion results in a loss of TelAP1 colocalization with TbTRF in BSF. In conclusion, this study demonstrates that characterizing potential TelAPs is effective in gaining insights into the telomeric complex's composition and its role in VSG regulation in Trypanosoma brucei. Understanding these interactions could potentially lead to new therapeutic targets for combatting African trypanosomiasis.}, subject = {Telomer }, language = {en} } @article{AmatobiOzbekUnalSchaebleretal.2023, author = {Amatobi, Kelechi M. and Ozbek-Unal, Ayten Gizem and Sch{\"a}bler, Stefan and Deppisch, Peter and Helfrich-F{\"o}rster, Charlotte and Mueller, Martin J. and Wegener, Christian and Fekete, Agnes}, title = {The circadian clock is required for rhythmic lipid transport in Drosophila in interaction with diet and photic condition}, series = {Journal of Lipid Research}, volume = {64}, journal = {Journal of Lipid Research}, number = {10}, doi = {10.1016/j.jlr.2023.100417}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-349961}, pages = {100417}, year = {2023}, abstract = {Modern lifestyle is often at odds with endogenously driven rhythmicity, which can lead to circadian disruption and metabolic syndrome. One signature for circadian disruption is a reduced or altered metabolite cycling in the circulating tissue reflecting the current metabolic status. Drosophila is a well-established model in chronobiology, but day-time dependent variations of transport metabolites in the fly circulation are poorly characterized. Here, we sampled fly hemolymph throughout the day and analyzed diacylglycerols (DGs), phosphoethanolamines (PEs) and phosphocholines (PCs) using LC-MS. In wild-type flies kept on sugar-only medium under a light-dark cycle, all transport lipid species showed a synchronized bimodal oscillation pattern with maxima at the beginning and end of the light phase which were impaired in period01 clock mutants. In wild-type flies under constant dark conditions, the oscillation became monophasic with a maximum in the middle of the subjective day. In strong support of clock-driven oscillations, levels of the targeted lipids peaked once in the middle of the light phase under time-restricted feeding independent of the time of food intake. When wild-type flies were reared on full standard medium, the rhythmic alterations of hemolymph lipid levels were greatly attenuated. Our data suggest that the circadian clock aligns daily oscillations of DGs, PEs, and PCs in the hemolymph to the anabolic siesta phase, with a strong influence of light on phase and modality.}, language = {en} } @article{SchuhmannScheiner2023, author = {Schuhmann, Antonia and Scheiner, Ricarda}, title = {A combination of the frequent fungicides boscalid and dimoxystrobin with the neonicotinoid acetamiprid in field-realistic concentrations does not affect sucrose responsiveness and learning behavior of honeybees}, series = {Ecotoxicology and Environmental Safety}, volume = {256}, journal = {Ecotoxicology and Environmental Safety}, doi = {10.1016/j.ecoenv.2023.114850}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-350047}, year = {2023}, abstract = {The increasing loss of pollinators over the last decades has become more and more evident. Intensive use of plant protection products is one key factor contributing to this decline. Especially the mixture of different plant protection products can pose an increased risk for pollinators as synergistic effects may occur. In this study we investigated the effect of the fungicide Cantus® Gold (boscalid/dimoxystrobin), the neonicotinoid insecticide Mospilan® (acetamiprid) and their mixture on honeybees. Since both plant protection products are frequently applied sequentially to the same plants (e.g. oilseed rape), their combination is a realistic scenario for honeybees. We investigated the mortality, the sucrose responsiveness and the differential olfactory learning performance of honeybees under controlled conditions in the laboratory to reduce environmental noise. Intact sucrose responsiveness and learning performance are of pivotal importance for the survival of individual honeybees as well as for the functioning of the entire colony. Treatment with two sublethal and field relevant concentrations of each plant protection product did not lead to any significant effects on these behaviors but affected the mortality rate. However, our study cannot exclude possible negative sublethal effects of these substances in higher concentrations. In addition, the honeybee seems to be quite robust when it comes to effects of plant protection products, while wild bees might be more sensitive. Highlights • Mix of SBI fungicides and neonicotinoids can lead to synergistic effects for bees. • Combination of non-SBI fungicide and neonicotinoid in field-realistic doses tested. • Synergistic effect on mortality of honeybees. • No effects on sucrose responsiveness and learning performance of honeybees. • Synergistic effects by other pesticide mixtures or on wild bees cannot be excluded.}, language = {en} } @article{CoelhoAlvesMonteiroetal.2019, author = {Coelho, Luis Pedro and Alves, Renato and Monteiro, Paulo and Huerta-Cepas, Jaime and Freitas, Ana Teresa and Bork, Peer}, title = {NG-meta-profiler: fast processing of metagenomes using NGLess, a domain-specific language}, series = {Microbiome}, volume = {7}, journal = {Microbiome}, number = {84}, doi = {10.1186/s40168-019-0684-8}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-223161}, year = {2019}, abstract = {Background Shotgun metagenomes contain a sample of all the genomic material in an environment, allowing for the characterization of a microbial community. In order to understand these communities, bioinformatics methods are crucial. A common first step in processing metagenomes is to compute abundance estimates of different taxonomic or functional groups from the raw sequencing data. Given the breadth of the field, computational solutions need to be flexible and extensible, enabling the combination of different tools into a larger pipeline. Results We present NGLess and NG-meta-profiler. NGLess is a domain specific language for describing next-generation sequence processing pipelines. It was developed with the goal of enabling user-friendly computational reproducibility. It provides built-in support for many common operations on sequencing data and is extensible with external tools with configuration files. Using this framework, we developed NG-meta-profiler, a fast profiler for metagenomes which performs sequence preprocessing, mapping to bundled databases, filtering of the mapping results, and profiling (taxonomic and functional). It is significantly faster than either MOCAT2 or htseq-count and (as it builds on NGLess) its results are perfectly reproducible. Conclusions NG-meta-profiler is a high-performance solution for metagenomics processing built on NGLess. It can be used as-is to execute standard analyses or serve as the starting point for customization in a perfectly reproducible fashion. NGLess and NG-meta-profiler are open source software (under the liberal MIT license) and can be downloaded from https://ngless.embl.de or installed through bioconda.}, language = {en} } @article{CoelhoKultimaCosteaetal.2018, author = {Coelho, Luis Pedro and Kultima, Jens Roat and Costea, Paul Igor and Fournier, Coralie and Pan, Yuanlong and Czarnecki-Maulden, Gail and Hayward, Matthew Robert and Forslund, Sofia K. and Schmidt, Thomas Sebastian Benedikt and Descombes, Patrick and Jackson, Janet R. and Li, Qinghong and Bork, Peer}, title = {Similarity of the dog and human gut microbiomes in gene content and response to diet}, series = {Microbiome}, volume = {6}, journal = {Microbiome}, doi = {10.1186/s40168-018-0450-3}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-223177}, year = {2018}, abstract = {Background Gut microbes influence their hosts in many ways, in particular by modulating the impact of diet. These effects have been studied most extensively in humans and mice. In this work, we used whole genome metagenomics to investigate the relationship between the gut metagenomes of dogs, humans, mice, and pigs. Results We present a dog gut microbiome gene catalog containing 1,247,405 genes (based on 129 metagenomes and a total of 1.9 terabasepairs of sequencing data). Based on this catalog and taxonomic abundance profiling, we show that the dog microbiome is closer to the human microbiome than the microbiome of either pigs or mice. To investigate this similarity in terms of response to dietary changes, we report on a randomized intervention with two diets (high-protein/low-carbohydrate vs. lower protein/higher carbohydrate). We show that diet has a large and reproducible effect on the dog microbiome, independent of breed or sex. Moreover, the responses were in agreement with those observed in previous human studies. Conclusions We conclude that findings in dogs may be predictive of human microbiome results. In particular, a novel finding is that overweight or obese dogs experience larger compositional shifts than lean dogs in response to a high-protein diet.}, language = {en} }