Highly variable biological effects of statins on cancer, non-cancer, and stem cells in vitro.

Statins, the drugs used for the treatment of hypercholesterolemia, have come into the spotlight not only as chemoadjuvants, but also as potential stem cell modulators in the context of regenerative therapy. In our study, we compared the in vitro effects of all clinically used statins on the viability of human pancreatic cancer (MiaPaCa-2) cells, non-cancerous human embryonic kidney (HEK 293) cells and adipose-derived mesenchymal stem cells (ADMSC). Additionally, the effect of statins on viability of MiaPaCa-2 and ADMSC cells spheroids was tested. Furthermore, we performed a microarray analysis on ADMSCs treated with individual statins (12 µM) and compared the importance of the effects of statins on gene expression between stem cells and pancreatic cancer cells. Concentrations of statins that significantly affected cancer cells viability (< 40 µM) did not affect stem cells viability after 24 h. Moreover, statins that didn´t affect viability of cancer cells grown in a monolayer, induce the disintegration of cancer cell spheroids. The effect of statins on gene expression was significantly less pronounced in stem cells compared to pancreatic cancer cells. In conclusion, the low efficacy of statins on non-tumor and stem cells at concentrations sufficient for cancer cells growth inhibition, support their applicability in chemoadjuvant tumor therapy.

Astrocyte-like subpopulation of NG2 glia in the adult mouse cortex exhibits characteristics of neural progenitor cells.

Glial cells expressing neuron-glial antigen 2 (NG2), also known as oligodendrocyte progenitor cells (OPCs), play a critical role in maintaining brain health. However, their ability to differentiate after ischemic injury is poorly understood. The aim of this study was to investigate the properties and functions of NG2 glia in the ischemic brain. Using transgenic mice, we selectively labeled NG2-expressing cells and their progeny in both healthy brain and after focal cerebral ischemia (FCI). Using single-cell RNA sequencing, we classified the labeled glial cells into five distinct subpopulations based on their gene expression patterns. Additionally, we examined the membrane properties of these cells using the patch-clamp technique. Of the identified subpopulations, three were identified as OPCs, whereas the fourth subpopulation had characteristics indicative of cells likely to develop into oligodendrocytes. The fifth subpopulation of NG2 glia showed astrocytic markers and had similarities to neural progenitor cells. Interestingly, this subpopulation was present in both healthy and post-ischemic tissue; however, its gene expression profile changed after ischemia, with increased numbers of genes related to neurogenesis. Immunohistochemical analysis confirmed the temporal expression of neurogenic genes and showed an increased presence of NG2 cells positive for Purkinje cell protein-4 at the periphery of the ischemic lesion 12 days after FCI, as well as NeuN-positive NG2 cells 28 and 60 days after injury. These results suggest the potential development of neuron-like cells arising from NG2 glia in the ischemic tissue. Our study provides insights into the plasticity of NG2 glia and their capacity for neurogenesis after stroke.

Terminal differentiation of villus tip enterocytes is governed by distinct Tgfß superfamily members.

The protective and absorptive functions of the intestinal epithelium rely on differentiated enterocytes in the villi. The differentiation of enterocytes is orchestrated by sub-epithelial mesenchymal cells producing distinct ligands along the villus axis, in particular Bmps and Tgfß. Here, we show that individual Bmp ligands and Tgfß drive distinct enterocytic programs specific to villus zonation. Bmp4 is expressed from the centre to the upper part of the villus and activates preferentially genes connected to lipid uptake and metabolism. In contrast, Bmp2 is produced by villus tip mesenchymal cells and it influences the adhesive properties of villus tip epithelial cells and the expression of immunomodulators. Additionally, Tgfß induces epithelial gene expression programs similar to those triggered by Bmp2. Bmp2-driven villus tip program is activated by a canonical Bmp receptor type I/Smad-dependent mechanism. Finally, we establish an organoid cultivation system that enriches villus tip enterocytes and thereby better mimics the cellular composition of the intestinal epithelium. Our data suggest that not only a Bmp gradient but also the activity of individual Bmp drives specific enterocytic programs.

Scdrake: a reproducible and scalable pipeline for scRNA-seq data analysis.

Motivation: While the workflow for primary analysis of single-cell RNA-seq (scRNA-seq) data is well established, the secondary analysis of the feature-barcode matrix is usually done by custom scripts. There is no fully automated pipeline in the R statistical environment, which would follow the current best programming practices and requirements for reproducibility. Results: We have developed scdrake, a fully automated workflow for secondary analysis of scRNA-seq data, which is fully implemented in the R language and built within the drake framework. The pipeline includes quality control, cell and gene filtering, normalization, detection of highly variable genes, dimensionality reduction, clustering, cell type annotation, detection of marker genes, differential expression analysis and integration of multiple samples. The pipeline is reproducible and scalable, has an efficient execution, provides easy extendability and access to intermediate results and outputs rich HTML reports. Scdrake is distributed as a Docker image, which provides a straightforward setup and enhances reproducibility. Availability and implementation: The source code and documentation are available under the MIT license at https://github.com/bioinfocz/scdrake and https://bioinfocz.github.io/scdrake, respectively. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Retinitis pigmentosa-associated mutations in mouse Prpf8 cause misexpression of circRNAs and degeneration of cerebellar granule cells.

A subset of patients with retinitis pigmentosa (RP) carry mutations in several spliceosomal components including the PRPF8 protein. Here, we established two alleles of murine Prpf8 that genocopy or mimic aberrant PRPF8 found in RP patients-the substitution p.Tyr2334Asn and an extended protein variant p.Glu2331ValfsX15. Homozygous mice expressing the aberrant Prpf8 variants developed within the first 2 mo progressive atrophy of the cerebellum because of extensive granule cell loss, whereas other cerebellar cells remained unaffected. We further show that a subset of circRNAs were deregulated in the cerebellum of both Prpf8-RP mouse strains. To identify potential risk factors that sensitize the cerebellum for Prpf8 mutations, we monitored the expression of several splicing proteins during the first 8 wk. We observed down-regulation of all selected splicing proteins in the WT cerebellum, which coincided with neurodegeneration onset. The decrease in splicing protein expression was further pronounced in mouse strains expressing mutated Prpf8. Collectively, we propose a model where physiological reduction in spliceosomal components during postnatal tissue maturation sensitizes cells to the expression of aberrant Prpf8 and the subsequent deregulation of circRNAs triggers neuronal death.

IL-17-driven induction of Paneth cell antimicrobial functions protects the host from microbiota dysbiosis and inflammation in the ileum.

Interleukin (IL)-17 protects epithelial barriers by inducing the secretion of antimicrobial peptides. However, the effect of IL-17 on Paneth cells (PCs), the major producers of antimicrobial peptides in the small intestine, is unclear. Here, we show that the targeted ablation of the IL-17 receptor (IL-17R) in PCs disrupts their antimicrobial functions and decreases the frequency of ileal PCs. These changes become more pronounced after colonization with IL-17 inducing segmented filamentous bacteria. Mice with PCs that lack IL-17R show an increased inflammatory transcriptional profile in the ileum along with the severity of experimentally induced ileitis. These changes are associated with a decrease in the diversity of gut microbiota that induces a severe ileum pathology upon transfer to genetically susceptible mice, which can be prevented by the systemic administration of IL-17a/f in microbiota recipients. In an exploratory analysis of a small cohort of pediatric patients with Crohn's disease, we have found that a portion of these patients exhibits a low number of lysozyme-expressing ileal PCs and a high ileitis severity score, resembling the phenotype of mice with IL-17R-deficient PCs. Our study identifies IL-17R-dependent signaling in PCs as an important mechanism that maintains ileal homeostasis through the prevention of dysbiosis.

Gut microbiota variation between climatic zones and due to migration strategy in passerine birds.

Introduction: Decreasing biotic diversity with increasing latitude is an almost universal macroecological pattern documented for a broad range of taxa, however, there have been few studies focused on changes in gut microbiota (GM) across climatic zones. Methods: Using 16S rRNA amplicon profiling, we analyzed GM variation between temperate (Czechia) and tropical (Cameroon) populations of 99 passerine bird species and assessed GM similarity of temperate species migrating to tropical regions with that of residents/short-distance migrants and tropical residents. Our study also considered the possible influence of diet on GM. Results: We observed no consistent GM diversity differences between tropical and temperate species. In the tropics, GM composition varied substantially between dry and rainy seasons and only a few taxa exhibited consistent differential abundance between tropical and temperate zones, irrespective of migration behavior and seasonal GM changes. During the breeding season, trans-Saharan migrant GM diverged little from species not overwintering in the tropics and did not show higher similarity to tropical passerines than temperate residents/short-distance migrants. Interestingly, GM of two temperate-breeding trans-Saharan migrants sampled in the tropical zone matched that of tropical residents and converged with other temperate species during the breeding season. Diet had a slight effect on GM composition of tropical species, but no effect on GM of temperate hosts. Discussion: Consequently, our results demonstrate extensive passerine GM plasticity, the dominant role of environmental factors in its composition and limited effect of diet.

Two complementary approaches for efficient isolation of Sertoli cells for transcriptomic analysis.

Sertoli cells (SCs) are the only somatic cells that reside in seminiferous tubules of testis. They directly interact with and support the development of germ cells, thus have an indispensable role in the process of spermatogenesis. SCs first appear in a proliferative state and then, with the initiation of the first wave of spermatogenesis, progress to a mature "nurturing" state which supports lifelong continuous sperm production. During this development, the SC transcriptome must adapt rapidly as obstacles in SC maturation often result in deficiencies in male fertility. Due to its importance in spermatogenesis, a reliable, rapid, and precise method for the isolation of high purity, viable and unadulterated SC has been largely missing. We have developed an improved method for the preparation of a testicular single cell suspension comprised of two alternative protocols to separate SCs from the rest of the testicular cells by FACS. The first sorting scheme is based on their co-expression of surface specific markers, FSHr and Occludin-1, while the second focuses on the co-staining of SCs with FSHr-specific antibody and Hoechst 33342, which discriminates DNA content of testicular cells. The entire procedure can be completed in less than 3 h which permits the analysis of the development-related transcriptional profile of these cells. Notably, our comparative study showed that this method resulted in a SC transcriptome that is largely comparable to SCs which were briskly isolated due to their cell-specific expression of fluorescent protein. Interestingly, we also show that SCs sorted as FSHr+Occludin+ cells contained a tangible portion of transcripts from all types of testicular germ cells. Sorting of SCs according to their 2C DNA content significantly reduced the presence of these transcripts, thus seems to be the most suitable approach for accurate determination of the SC transcriptome. We believe that these novel approaches for the isolation of SCs will assist researchers in the elucidation of their function as well as their role in spermatogenesis and disorders related to male infertility.

TROP2 Represents a Negative Prognostic Factor in Colorectal Adenocarcinoma and Its Expression Is Associated with Features of Epithelial-Mesenchymal Transition and Invasiveness.

Trophoblastic cell surface antigen 2 (TROP2) is a membrane glycoprotein overexpressed in many solid tumors with a poor prognosis, including intestinal neoplasms. In our study, we show that TROP2 is expressed in preneoplastic lesions, and its expression is maintained in most colorectal cancers (CRC). High TROP2 positivity correlated with lymph node metastases and poor tumor differentiation and was a negative prognostic factor. To investigate the role of TROP2 in intestinal tumors, we analyzed two mouse models with conditional disruption of the adenomatous polyposis coli (Apc) tumor-suppressor gene, human adenocarcinoma samples, patient-derived organoids, and TROP2-deficient tumor cells. We found that Trop2 is produced early after Apc inactivation and its expression is associated with the transcription of genes involved in epithelial-mesenchymal transition, the regulation of migration, invasiveness, and extracellular matrix remodeling. A functionally similar group of genes was also enriched in TROP2-positive cells from human CRC samples. To decipher the driving mechanism of TROP2 expression, we analyzed its promoter. In human cells, this promoter was activated by ß-catenin and additionally by the Yes1-associated transcriptional regulator (YAP). The regulation of TROP2 expression by active YAP was verified by YAP knockdown in CRC cells. Our results suggest a possible link between aberrantly activated Wnt/ß-catenin signaling, YAP, and TROP2 expression.

Role of LGR5-positive mesenchymal cells in craniofacial development.

Leucine Rich Repeat Containing G Protein-Coupled Receptor 5 (LGR5), a Wnt pathway member, has been previously recognised as a stem cell marker in numerous epithelial tissues. In this study, we used Lgr5-EGFP-CreERT2 mice to analyse the distribution of LGR5-positive cells during craniofacial development. LGR5 expressing cells were primarily located in the mesenchyme adjacent to the craniofacial epithelial structures undergoing folding, such as the nasopharyngeal duct, lingual groove, and vomeronasal organ. To follow the fate of LGR5-positive cells, we performed lineage tracing using an inducible Cre knock-in allele in combination with Rosa26-tdTomato reporter mice. The slight expansion of LGR5-positive cells was found around the vomeronasal organ, in the nasal cavity, and around the epithelium in the lingual groove. However, most LGR5 expressing cells remained in their original location, possibly supporting their signalling function for adjacent epithelium rather than exerting their role as progenitor cells for the craniofacial structures. Moreover, Lgr5 knockout mice displayed distinct defects in LGR5-positive areas, especially in the reduction of the nasopharyngeal duct, the alteration of the palatal shelves shape, abnormal epithelial folding in the lingual groove area, and the disruption of salivary gland development. The latter defect manifested as an atypical number and localisation of the glandular ducts. The gene expression of several Wnt pathway members (Rspo1-3, Axin2) was altered in Lgr5-deficient animals. However, the difference was not found in sorted EGFP-positive cells obtained from Lgr5 +/+ and Lgr5 -/- animals. Expression profiling of LGR5-positive cells revealed the expression of several markers of mesenchymal cells, antagonists, as well as agonists, of Wnt signalling, and molecules associated with the basal membrane. Therefore, LGR5-positive cells in the craniofacial area represent a very specific population of mesenchymal cells adjacent to the epithelium undergoing folding or groove formation. Our results indicate a possible novel role of LGR5 in the regulation of morphogenetic processes during the formation of complex epithelial structures in the craniofacial areas, a role which is not related to the stem cell properties of LGR5-positive cells as was previously defined for various epithelial tissues.

BRAT1 links Integrator and defective RNA processing with neurodegeneration.

Mutations in BRAT1, encoding BRCA1-associated ATM activator 1, have been associated with neurodevelopmental and neurodegenerative disorders characterized by heterogeneous phenotypes with varying levels of clinical severity. However, the underlying molecular mechanisms of disease pathology remain poorly understood. Here, we show that BRAT1 tightly interacts with INTS9/INTS11 subunits of the Integrator complex that processes 3' ends of various noncoding RNAs and pre-mRNAs. We find that Integrator functions are disrupted by BRAT1 deletion. In particular, defects in BRAT1 impede proper 3' end processing of UsnRNAs and snoRNAs, replication-dependent histone pre-mRNA processing, and alter the expression of protein-coding genes. Importantly, impairments in Integrator function are also evident in patient-derived cells from BRAT1 related neurological disease. Collectively, our data suggest that defects in BRAT1 interfere with proper Integrator functions, leading to incorrect expression of RNAs and proteins, resulting in neurodegeneration.

Within-community variation of interspecific divergence patterns in passerine gut microbiota.

Gut microbiota (GM) often exhibit variation between different host species and co-divergence with hosts' phylogeny. Identifying these patterns is a key for understanding the mechanisms that shaped symbiosis between GM and its hosts. Therefore, both GM-host species specificity and GM-host co-divergence have been investigated by numerous studies. However, most of them neglected a possibility that different groups of bacteria within GM can vary in the tightness of their association with the host. Consequently, unlike most of these studies, we aimed to directly address how the strength of GM-host species specificity and GM-host co-divergence vary across different GM clades. We decomposed GM communities of 52 passerine species (394 individuals), characterized by 16S rRNA amplicon sequence variant (ASV) profiles, into monophyletic Binned Taxonomic units (BTUs). Subsequently, we analyzed strength of host species specificity and correlation with host phylogeny separately for resulting BTUs. We found that most BTUs exhibited significant host-species specificity in their composition. Notably, BTUs exhibiting high host-species specificity comprised bacterial taxa known to impact host's physiology and immune system. However, BTUs rarely displayed significant co-divergence with host phylogeny, suggesting that passerine GM evolution is not shaped primarily through a shared evolutionary history between the host and its gut microbes.

Retinitis pigmentosa-linked mutation in DHX38 modulates its splicing activity.

Retinitis pigmentosa (RP) is a hereditary disease affecting tens of thousands of people world-wide. Here we analyzed the effect of an amino acid substitution in the RNA helicase DHX38 (Prp16) causing RP. DHX38 has been proposed as the helicase important for the 2nd step of splicing. We showed that DHX38 associates with key splicing factors involved in both splicing steps but did not find any evidence that the RP mutations changes DHX38 interaction profile with the spliceosome. We further downregulated DHX38 and monitored changes in splicing. We observed only minor perturbations of general splicing but detected modulation of ~70 alternative splicing events. Next, we probed DHX38 function in splicing of retina specific genes and found that FSCN2 splicing is dependent on DHX38. In addition, RHO splicing was inhibited specifically by expression of DHX38 RP variant. Finally, we showed that overexpression of DHX38 promotes usage of canonical as well as cryptic 5' splice sites in HBB splicing reporter. Together, our data show that DHX38 is a splicing factor that promotes splicing of cryptic splice sites and regulate alternative splicing. We further provide evidence that the RP-linked substitution G332D modulates DHX38 splicing activity.

Sodium Accumulation and Blood Capillary Rarefaction in the Skin Predispose Spontaneously Hypertensive Rats to Salt Sensitive Hypertension.

Recent studies in humans and rats suggested that increased Na+ storage in the skin without parallel water retention may predispose to salt-sensitive hypertension. In the current studies, we compared tissue Na+ storage in salt sensitive spontaneously hypertensive rats (SHR) versus salt resistant normotensive Brown Norway (BN-Lx) rats. After salt loading (10 days drinking 1% NaCl solution), the SHR showed significant parallel increase in Na+-to-water as well as (Na++K+)-to-water ratios suggesting increased storage of osmotically inactive Na+ in the skin while no significant changes in skin electrolyte concentrations were observed in BN-Lx rats. SHR rats after salt treatment exhibited a nonsignificant decrease in skin blood capillary number (rarefaction) while BN-Lx rats showed significantly increased skin blood capillary density. Analysis of dermal gene expression profiles in BN-Lx rats after salt treatment showed significant up-regulation of genes involved in angiogenesis and proliferation of endothelial cells contrary to the SHR. Since the skin harbors most of the body's resistance vessels it is possible that blood capillary rarefaction may lead to increased peripheral resistance and salt sensitivity in the SHR.

Expression of Interferons Lambda 3 and 4 Induces Identical Response in Human Liver Cell Lines Depending Exclusively on Canonical Signaling.

Lambda interferons mediate antiviral immunity by inducing interferon-stimulated genes (ISGs) in epithelial tissues. A common variant rs368234815TT/∆G creating functional gene from an IFNL4 pseudogene is associated with the expression of major ISGs in the liver but impaired clearance of hepatitis C. To explain this, we compared Halo-tagged and non-tagged IFNL3 and IFNL4 signaling in liver-derived cell lines. Transfection with non-tagged IFNL3, non-tagged IFNL4 and Halo-tagged IFNL4 led to a similar degree of JAK-STAT activation and ISG induction; however, the response to transfection with Halo-tagged IFNL3 was lower and delayed. Transfection with non-tagged IFNL3 or IFNL4 induced no transcriptome change in the cells lacking either IL10R2 or IFNLR1 receptor subunits. Cytosolic overexpression of signal peptide-lacking IFNL3 or IFNL4 in wild type cells did not interfere with JAK-STAT signaling triggered by interferons in the medium. Finally, expression profile changes induced by transfection with non-tagged IFNL3 and IFNL4 were highly similar. These data do not support the hypothesis about IFNL4-specific non-canonical signaling and point out that functional studies conducted with tagged interferons should be interpreted with caution.

Phenotypic and Clonal Stability of Antigen-Inexperienced Memory-like T Cells across the Genetic Background, Hygienic Status, and Aging.

Ag-inexperienced memory-like T (AIMT) cells are functionally unique T cells, representing one of the two largest subsets of murine CD8+ T cells. However, differences between laboratory inbred strains, insufficient data from germ-free mice, a complete lack of data from feral mice, and an unclear relationship between AIMT cells formation during aging represent major barriers for better understanding of their biology. We performed a thorough characterization of AIMT cells from mice of different genetic background, age, and hygienic status by flow cytometry and multiomics approaches, including analyses of gene expression, TCR repertoire, and microbial colonization. Our data showed that AIMT cells are steadily present in mice, independent of their genetic background and hygienic status. Despite differences in their gene expression profiles, young and aged AIMT cells originate from identical clones. We identified that CD122 discriminates two major subsets of AIMT cells in a strain-independent manner. Whereas thymic CD122LOW AIMT cells (innate memory) prevail only in young animals with high thymic IL-4 production, peripheral CD122HIGH AIMT cells (virtual memory) dominate in aged mice. Cohousing with feral mice changed the bacterial colonization of laboratory strains but had only minimal effects on the CD8+ T cell compartment, including AIMT cells.

Depletion of Retinal Dopaminergic Activity in a Mouse Model of Rod Dysfunction Exacerbates Experimental Autoimmune Uveoretinitis: A Role for the Gateway Reflex.

The gateway reflex is a mechanism by which neural inputs regulate chemokine expression at endothelial cell barriers, thereby establishing gateways for the invasion of autoreactive T cells into barrier-protected tissues. In this study, we hypothesized that rod photoreceptor dysfunction causes remodeling of retinal neural activity, which influences the blood-retinal barrier and the development of retinal inflammation. We evaluated this hypothesis using Gnat1rd17 mice, a model of night blindness with late-onset rod-cone dystrophy, and experimental autoimmune uveoretinitis (EAU). Retinal remodeling and its effect on EAU development were investigated by transcriptome profiling, target identification, and functional validation. We showed that Gnat1rd17 mice primarily underwent alterations in their retinal dopaminergic system, triggering the development of an exacerbated EAU, which was counteracted by dopamine replacement with L-DOPA administered either systemically or locally. Remarkably, dopamine acted on retinal endothelial cells to inhibit NF-κB and STAT3 activity and the expression of downstream target genes such as chemokines involved in T cell recruitment. These results suggest that rod-mediated dopamine release functions in a gateway reflex manner in the homeostatic control of immune cell entry into the retina, and the loss of retinal dopaminergic activity in conditions associated with rod dysfunction increases the susceptibility to autoimmune uveitis.

Pathogenic ARH3 mutations result in ADP-ribose chromatin scars during DNA strand break repair.

Neurodegeneration is a common hallmark of individuals with hereditary defects in DNA single-strand break repair; a process regulated by poly(ADP-ribose) metabolism. Recently, mutations in the ARH3 (ADPRHL2) hydrolase that removes ADP-ribose from proteins have been associated with neurodegenerative disease. Here, we show that ARH3-mutated patient cells accumulate mono(ADP-ribose) scars on core histones that are a molecular memory of recently repaired DNA single-strand breaks. We demonstrate that the ADP-ribose chromatin scars result in reduced endogenous levels of important chromatin modifications such as H3K9 acetylation, and that ARH3 patient cells exhibit measurable levels of deregulated transcription. Moreover, we show that the mono(ADP-ribose) scars are lost from the chromatin of ARH3-defective cells in the prolonged presence of PARP inhibition, and concomitantly that chromatin acetylation is restored to normal. Collectively, these data indicate that ARH3 can act as an eraser of ADP-ribose chromatin scars at sites of PARP activity during DNA single-strand break repair.

Toll-like receptor signaling in thymic epithelium controls monocyte-derived dendritic cell recruitment and Treg generation.

The development of thymic regulatory T cells (Treg) is mediated by Aire-regulated self-antigen presentation on medullary thymic epithelial cells (mTECs) and dendritic cells (DCs), but the cooperation between these cells is still poorly understood. Here we show that signaling through Toll-like receptors (TLR) expressed on mTECs regulates the production of specific chemokines and other genes associated with post-Aire mTEC development. Using single-cell RNA-sequencing, we identify a new thymic CD14+Sirpα+ population of monocyte-derived dendritic cells (CD14+moDC) that are enriched in the thymic medulla and effectively acquire mTEC-derived antigens in response to the above chemokines. Consistently, the cellularity of CD14+moDC is diminished in mice with MyD88-deficient TECs, in which the frequency and functionality of thymic CD25+Foxp3+ Tregs are decreased, leading to aggravated mouse experimental colitis. Thus, our findings describe a TLR-dependent function of mTECs for the recruitment of CD14+moDC, the generation of Tregs, and thereby the establishment of central tolerance.

Codiversification of gastrointestinal microbiota and phylogeny in passerines is not explained by ecological divergence.

Vertebrate gut microbiota (GM) is comprised of a taxonomically diverse consortium of symbiotic and commensal microorganisms that have a pronounced effect on host physiology, immune system function and health status. Despite much research on interactions between hosts and their GM, the factors affecting inter- and intraspecific GM variation in wild populations are still poorly known. We analysed data on faecal microbiota composition in 51 passerine species (319 individuals) using Illumina MiSeq sequencing of bacterial 16S rRNA (V3-V4 variable region). Despite pronounced interindividual variation, GM composition exhibited significant differences at the interspecific level, accounting for approximately 20%-30% of total GM variation. We also observed a significant correlation between GM composition divergence and host's phylogenetic divergence, with strength of correlation higher than that of GM vs. ecological or life history traits and geographic variation. The effect of host's phylogeny on GM composition was significant, even after statistical control for these confounding factors. Hence, our data do not support codiversification of GM and passerine phylogeny solely as a by-product of their ecological divergence. Furthermore, our findings do not support that GM vs. host's phylogeny codiversification is driven primarily through trans-generational GM transfer as the GM vs. phylogeny correlation does not increase with higher sequence similarity used when delimiting operational taxonomic units. Instead, we hypothesize that the GM vs. phylogeny correlation may arise as a consequence of interspecific divergence of genes that directly or indirectly modulate composition of GM.