Linking the Human Gut Microbiome to Inflammatory Cytokine Production Capacity.

Gut microbial dysbioses are linked to aberrant immune responses, which are often accompanied by abnormal production of inflammatory cytokines. As part of the Human Functional Genomics Project (HFGP), we investigate how differences in composition and function of gut microbial communities may contribute to inter-individual variation in cytokine responses to microbial stimulations in healthy humans. We observe microbiome-cytokine interaction patterns that are stimulus specific, cytokine specific, and cytokine and stimulus specific. Validation of two predicted host-microbial interactions reveal that TNFα and IFNγ production are associated with specific microbial metabolic pathways: palmitoleic acid metabolism and tryptophan degradation to tryptophol. Besides providing a resource of predicted microbially derived mediators that influence immune phenotypes in response to common microorganisms, these data can help to define principles for understanding disease susceptibility. The three HFGP studies presented in this issue lay the groundwork for further studies aimed at understanding the interplay between microbial, genetic, and environmental factors in the regulation of the immune response in humans. PAPERCLIP.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation.

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

A rising tide of parasite transcriptomics propels pathogen biology.

Twenty years ago, the first transcriptome of the intraerythrocytic developmental cycle of the malaria parasite Plasmodium falciparum was published in PLOS Biology. Since then, transcriptomics studies have transformed the study of parasite biology.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles.

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis.

Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.

Cardelino: computational integration of somatic clonal substructure and single-cell transcriptomes.

Bulk and single-cell DNA sequencing has enabled reconstructing clonal substructures of somatic tissues from frequency and cooccurrence patterns of somatic variants. However, approaches to characterize phenotypic variations between clones are not established. Here we present cardelino (https://github.com/single-cell-genetics/cardelino), a computational method for inferring the clonal tree configuration and the clone of origin of individual cells assayed using single-cell RNA-seq (scRNA-seq). Cardelino flexibly integrates information from imperfect clonal trees inferred based on bulk exome-seq data, and sparse variant alleles expressed in scRNA-seq data. We apply cardelino to a published cancer dataset and to newly generated matched scRNA-seq and exome-seq data from 32 human dermal fibroblast lines, identifying hundreds of differentially expressed genes between cells from different somatic clones. These genes are frequently enriched for cell cycle and proliferation pathways, indicating a role for cell division genes in somatic evolution in healthy skin.

Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity.

Approximately 1.5 billion people worldwide are overweight or affected by obesity, and are at risk of developing type 2 diabetes, cardiovascular disease and related metabolic and inflammatory disturbances. Although the mechanisms linking adiposity to associated clinical conditions are poorly understood, recent studies suggest that adiposity may influence DNA methylation, a key regulator of gene expression and molecular phenotype. Here we use epigenome-wide association to show that body mass index (BMI; a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci with P < 1 × 10-7, range P = 9.2 × 10-8 to 6.0 × 10-46; n = 10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find that methylation loci are enriched for functional genomic features in multiple tissues (P < 0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P < 9.0 × 10-6, range P = 5.5 × 10-6 to 6.1 × 10-35, n = 1,785 samples). The methylation loci identify genes involved in lipid and lipoprotein metabolism, substrate transport and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future development of type 2 diabetes (relative risk per 1 standard deviation increase in methylation risk score: 2.3 (2.07-2.56); P = 1.1 × 10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type 2 diabetes and other adverse clinical consequences of obesity.

Contralateral hearing aid use in adult cochlear implant recipients: retrospective analysis of auditory outcomes.

OBJECTIVE: To investigate retrospectively the frequency of usage of bimodal stimulation among cochlear implant (CI) users, as well its clinical benefit relative to unilateral use. DESIGN: All subjects had been monitored with the clinical Minimal Outcome Measurements test battery. STUDY SAMPLES: 103 adults with bilateral postlingual profound sensorineural hearing loss and unilateral CI use were extracted from the local database. These were divided into two groups: those who only used a CI and those who used bimodal stimulation. RESULTS: The preoperative contralateral residual hearing in the bimodal group was significantly better than that of the CI-only group. In both groups, speech perception in quiet and in noise improved after CI, with no significant difference between postoperative unimodal conditions. For the bimodal group, an additional significant improvement was found for the bimodal condition compared to the unimodal. CONCLUSION: Given the observed auditory benefit of bimodal stimulation in comparison to unimodal stimulation and given the finding that degree of residual hearing is not correlated with bimodal benefits, it is recommended to encourage CI recipients to continue contralateral HA use after CI. As a result of expanding CI criteria worldwide, the population of bimodal users is expected to grow in the near future.

A Role for Basigin in Toxoplasma gondii Infection.

The role of specific host cell surface receptors during Toxoplasma gondii invasion of host cells is poorly defined. Here, we interrogated the role of the well-known malarial invasion receptor, basigin, in T. gondii infection of astrocytes. We found that primary astrocytes express two members of the BASIGIN (BSG) immunoglobulin family, basigin and embigin, but did not express neuroplastin. Antibody blockade of either basigin or embigin caused a significant reduction of parasite infectivity in astrocytes. The specific role of basigin during T. gondii invasion was further examined using a mouse astrocytic cell line (C8-D30), which exclusively expresses basigin. CRISPR-mediated deletion of basigin in C8-D30 cells resulted in decreased T. gondii infectivity. T. gondii replication and invasion efficiency were not altered by basigin deficiency, but parasite attachment to astrocytes was markedly reduced. We also conducted a proteomic screen to identify T. gondii proteins that interact with basigin. Toxoplasma-encoded cyclophilins, the protein 14-3-3, and protein disulfide isomerase (TgPDI) were among the putative basigin-ligands identified. Recombinant TgPDI produced in E. coli bound to basigin and pretreatment of tachyzoites with a PDI inhibitor decreased parasite attachment to host cells. Finally, mutagenesis of the active site cysteines of TgPDI abolished enzyme binding to basigin. Thus, basigin and its related immunoglobulin family members may represent host receptors that mediate attachment of T. gondii to diverse cell types.

The apical annuli of Toxoplasma gondii are composed of coiled-coil and signalling proteins embedded in the inner membrane complex sutures.

The apical annuli are among the most intriguing and understudied structures in the cytoskeleton of the apicomplexan parasite Toxoplasma gondii. We mapped the proteome of the annuli in Toxoplasma by reciprocal proximity biotinylation (BioID), and validated five apical annuli proteins (AAP1-5), Centrin2, and an apical annuli methyltransferase. Moreover, inner membrane complex (IMC) suture proteins connecting the alveolar vesicles were also detected and support annuli residence within the sutures. Super-resolution microscopy identified a concentric organisation comprising four rings with diameters ranging from 200 to 400 nm. The high prevalence of domain signatures shared with centrosomal proteins in the AAPs together with Centrin2 suggests that the annuli are related and/or derived from the centrosomes. Phylogenetic analysis revealed that the AAPs are conserved narrowly in coccidian, apicomplexan parasites that multiply by an internal budding mechanism. This suggests a role in replication, for example, to provide pores in the mother IMC permitting exchange of building blocks and waste products. However, presence of multiple signalling domains and proteins are suggestive of additional functions. Knockout of AAP4, the most conserved compound forming the largest ring-like structure, modestly decreased parasite fitness in vitro but had no significant impact on acute virulence in vivo. In conclusion, the apical annuli are composed of coiled-coil and signalling proteins assembled in a pore-like structure crossing the IMC barrier maintained during internal budding.

Lytic Cycle of Toxoplasma gondii: 15 Years Later.

Toxoplasmosis is the clinical and pathological consequence of acute infection with the obligate intracellular apicomplexan parasite Toxoplasma gondii. Symptoms result from tissue destruction that accompanies lytic parasite growth. This review updates current understanding of the host cell invasion, parasite replication, and eventual egress that constitute the lytic cycle, as well as the ways T. gondii manipulates host cells to ensure its survival. Since the publication of a previous iteration of this review 15 years ago, important advances have been made in our molecular understanding of parasite growth and mechanisms of host cell egress, and knowledge of the parasite's manipulation of the host has rapidly progressed. Here we cover molecular advances and current conceptual frameworks that include each of these topics, with an eye to what may be known 15 years from now.

Experimental and theoretical evidence for bilayer-by-bilayer surface melting of crystalline ice.

On the surface of water ice, a quasi-liquid layer (QLL) has been extensively reported at temperatures below its bulk melting point at 273 K. Approaching the bulk melting temperature from below, the thickness of the QLL is known to increase. To elucidate the precise temperature variation of the QLL, and its nature, we investigate the surface melting of hexagonal ice by combining noncontact, surface-specific vibrational sum frequency generation (SFG) spectroscopy and spectra calculated from molecular dynamics simulations. Using SFG, we probe the outermost water layers of distinct single crystalline ice faces at different temperatures. For the basal face, a stepwise, sudden weakening of the hydrogen-bonded structure of the outermost water layers occurs at 257 K. The spectral calculations from the molecular dynamics simulations reproduce the experimental findings; this allows us to interpret our experimental findings in terms of a stepwise change from one to two molten bilayers at the transition temperature.

Toxoplasma gondii chromosomal passenger complex is essential for the organization of a functional mitotic spindle: a prerequisite for productive endodyogeny.

The phylum Apicomplexa encompasses deadly pathogens such as malaria and Cryptosporidium. Apicomplexa cell division is mechanistically divergent from that of their mammalian host, potentially representing an attractive source of drug targets. Depending on the species, apicomplexan parasites can modulate the output of cell division, producing two to thousands of daughter cells at once. The inherent flexibility of their cell division mechanisms allows these parasites to adapt to different niches, facilitating their dissemination. Toxoplasma gondii tachyzoites divide using a unique form of cell division called endodyogeny. This process involves a single round of DNA replication, closed nuclear mitosis, and assembly of two daughter cells within a mother. In higher Eukaryotes, the four-subunit chromosomal passenger complex (CPC) (Aurora kinase B (ARKB)/INCENP/Borealin/Survivin) promotes chromosome bi-orientation by detaching incorrect kinetochore-microtubule attachments, playing an essential role in controlling cell division fidelity. Herein, we report the characterization of the Toxoplasma CPC (Aurora kinase 1 (Ark1)/INCENP1/INCENP2). We show that the CPC exhibits dynamic localization in a cell cycle-dependent manner. TgArk1 interacts with both TgINCENPs, with TgINCENP2 being essential for its translocation to the nucleus. While TgINCENP1 appears to be dispensable, interfering with TgArk1 or TgINCENP2 results in pronounced division and growth defects. Significant anti-cancer drug development efforts have focused on targeting human ARKB. Parasite treatment with low doses of hesperadin, a known inhibitor of human ARKB at higher concentrations, phenocopies the TgArk1 and TgINCENP2 mutants. Overall, our study provides new insights into the mechanisms underpinning cell cycle control in Apicomplexa, and highlights TgArk1 as potential drug target.

Comparative Adsorption of Acetone on Water and Ice Surfaces.

Small organic molecules on ice and water surfaces are ubiquitous in nature and play a crucial role in many environmentally relevant processes. Herein, we combine surface-specific vibrational spectroscopy and a controllable flow cell apparatus to investigate the molecular adsorption of acetone onto the basal plane of single-crystalline hexagonal ice with a large surface area. By comparing the adsorption of acetone on the ice/air and the water/air interface, we observed two different types of acetone adsorption, as apparent from the different responses of both the free O-H and the hydrogen-bonded network vibrations for ice and liquid water. Adsorption on ice occurs preferentially through interactions with the free OH group, while the interaction of acetone with the surface of liquid water appears less specific.

Interplay of host genetics and gut microbiota underlying the onset and clinical presentation of inflammatory bowel disease.

OBJECTIVE: Patients with IBD display substantial heterogeneity in clinical characteristics. We hypothesise that individual differences in the complex interaction of the host genome and the gut microbiota can explain the onset and the heterogeneous presentation of IBD. Therefore, we performed a case-control analysis of the gut microbiota, the host genome and the clinical phenotypes of IBD. DESIGN: Stool samples, peripheral blood and extensive phenotype data were collected from 313 patients with IBD and 582 truly healthy controls, selected from a population cohort. The gut microbiota composition was assessed by tag-sequencing the 16S rRNA gene. All participants were genotyped. We composed genetic risk scores from 11 functional genetic variants proven to be associated with IBD in genes that are directly involved in the bacterial handling in the gut: NOD2, CARD9, ATG16L1, IRGM and FUT2. RESULTS: Strikingly, we observed significant alterations of the gut microbiota of healthy individuals with a high genetic risk for IBD: the IBD genetic risk score was significantly associated with a decrease in the genus Roseburia in healthy controls (false discovery rate 0.017). Moreover, disease location was a major determinant of the gut microbiota: the gut microbiota of patients with colonic Crohn's disease (CD) is different from that of patients with ileal CD, with a decrease in alpha diversity associated to ileal disease (p=3.28×10-13). CONCLUSIONS: We show for the first time that genetic risk variants associated with IBD influence the gut microbiota in healthy individuals. Roseburia spp are acetate-to-butyrate converters, and a decrease has already been observed in patients with IBD.

Membrane skeletal association and post-translational allosteric regulation of Toxoplasma gondii GAPDH1.

When Toxoplasma gondii egresses from the host cell, glyceraldehyde-3-phosphate dehydrogenase 1 (GAPDH1), which is primary a glycolysis enzyme but actually a quintessential multifunctional protein, translocates to the unique cortical membrane skeleton. Here, we report the 2.25 Å resolution crystal structure of the GAPDH1 holoenzyme in a quaternary complex providing the basis for the molecular dissection of GAPDH1 structure-function relationships Knockdown of GAPDH1 expression and catalytic site disruption validate the essentiality of GAPDH1 in intracellular replication but we confirmed that glycolysis is not strictly essential. We identify, for the first time, S-loop phosphorylation as a novel, critical regulator of enzymatic activity that is consistent with the notion that the S-loop is critical for cofactor binding, allosteric activation and oligomerization. We show that neither enzymatic activity nor phosphorylation state correlate with the ability to translocate to the cortex. However, we demonstrate that association of GAPDH1 with the cortex is mediated by the N-terminus, likely palmitoylation. Overall, glycolysis and cortical translocation are functionally decoupled by post-translational modifications.

The emerging landscape of dynamic DNA methylation in early childhood.

BACKGROUND: DNA methylation has been found to associate with disease, aging and environmental exposure, but it is unknown how genome, environment and disease influence DNA methylation dynamics in childhood. RESULTS: By analysing 538 paired DNA blood samples from children at birth and at 4-5 years old and 726 paired samples from children at 4 and 8 years old from four European birth cohorts using the Illumina Infinium Human Methylation 450 k chip, we have identified 14,150 consistent age-differential methylation sites (a-DMSs) at epigenome-wide significance of p < 1.14 × 10-7. Genes with an increase in age-differential methylation were enriched in pathways related to 'development', and were more often located in bivalent transcription start site (TSS) regions, which can silence or activate expression of developmental genes. Genes with a decrease in age-differential methylation were involved in cell signalling, and enriched on H3K27ac, which can predict developmental state. Maternal smoking tended to decrease methylation levels at the identified da-DMSs. We also found 101 a-DMSs (0.71%) that were regulated by genetic variants using cis-differential Methylation Quantitative Trait Locus (cis-dMeQTL) mapping. Moreover, a-DMS-associated genes during early development were significantly more likely to be linked with disease. CONCLUSION: Our study provides new insights into the dynamic epigenetic landscape of the first 8 years of life.

A MORN1-associated HAD phosphatase in the basal complex is essential for Toxoplasma gondii daughter budding.

Apicomplexan parasites replicate by several budding mechanisms with two well-characterized examples being Toxoplasma endodyogeny and Plasmodium schizogony. Completion of budding requires the tapering of the nascent daughter buds toward the basal end, driven by contraction of the basal complex. This contraction is not executed by any of the known cell division associated contractile mechanisms and in order to reveal new components of the unusual basal complex we performed a yeast two-hybrid screen with its major scaffolding protein, TgMORN1. Here we report on a conserved protein with a haloacid dehalogenase (HAD) phosphatase domain, hereafter named HAD2a, identified by yeast two-hybrid. HAD2a has demonstrated enzyme-activity in vitro, localizes to the nascent daughter buds, and co-localizes with MORN1 to the basal complex during its contraction. Conditional knockout of HAD2a in Toxoplasma interferes with basal complex assembly, which leads to incomplete cytokinesis and conjoined daughters that ultimately results in disrupted proliferation. In Plasmodium, we further confirmed localization of the HAD2a ortholog to the basal complex toward the end of schizogony. In conclusion, our work highlights an essential role for this HAD phosphatase across apicomplexan budding and suggests a regulatory mechanism of differential phosphorylation on the structure and/or contractile function of the basal complex.

The conserved apicomplexan Aurora kinase TgArk3 is involved in endodyogeny, duplication rate and parasite virulence.

Aurora kinases are eukaryotic serine/threonine protein kinases that regulate key events associated with chromatin condensation, centrosome and spindle function and cytokinesis. Elucidating the roles of Aurora kinases in apicomplexan parasites is crucial to understand the cell cycle control during Plasmodium schizogony or Toxoplasma endodyogeny. Here, we report on the localization of two previously uncharacterized Toxoplasma Aurora-related kinases (Ark2 and Ark3) in tachyzoites and of the uncharacterized Ark3 orthologue in Plasmodium falciparum erythrocytic stages. In Toxoplasma gondii, we show that TgArk2 and TgArk3 concentrate at specific sub-cellular structures linked to parasite division: the mitotic spindle and intranuclear mitotic structures (TgArk2), and the outer core of the centrosome and the budding daughter cells cytoskeleton (TgArk3). By tagging the endogenous PfArk3 gene with the green fluorescent protein in live parasites, we show that PfArk3 protein expression peaks late in schizogony and localizes at the periphery of budding schizonts. Disruption of the TgArk2 gene reveals no essential function for tachyzoite propagation in vitro, which is surprising giving that the P. falciparum and P. berghei orthologues are essential for erythrocyte schizogony. In contrast, knock-down of TgArk3 protein results in pronounced defects in parasite division and a major growth deficiency. TgArk3-depleted parasites display several defects, such as reduced parasite growth rate, delayed egress and parasite duplication, defect in rosette formation, reduced parasite size and invasion efficiency and lack of virulence in mice. Our study provides new insights into cell cycle control in Toxoplasma and malaria parasites and highlights Aurora kinase 3 as potential drug target.