Temporal composition of the cervicovaginal microbiome associates with hrHPV infection outcomes in a longitudinal study.

BACKGROUND: Persistent infections with high-risk human papillomavirus (hrHPV) can cause cervical squamous intraepithelial lesions (SIL) that may progress to cancer. The cervicovaginal microbiome (CVM) correlates with SIL, but the temporal composition of the CVM after hrHPV infections has not been fully clarified. METHODS: To determine the association between the CVM composition and infection outcome, we applied high-resolution microbiome profiling using the circular probe-based RNA sequencing technology on a longitudinal cohort of cervical smears obtained from 141 hrHPV DNA-positive women with normal cytology at first visit, of whom 51 were diagnosed by cytology with SIL six months later. RESULTS: Here we show that women with a microbial community characterized by low diversity and high Lactobacillus crispatus abundance at both visits exhibit low risk to SIL development, while women with a microbial community characterized by high diversity and Lactobacillus depletion at first visit have a higher risk of developing SIL. At the level of individual species, we observed that a high abundance for Gardnerella vaginalis and Atopobium vaginae at both visits associate with SIL outcomes. These species together with Dialister micraerophilus showed a moderate discriminatory power for hrHPV infection progression. CONCLUSIONS: Our results suggest that the CVM can potentially be used as a biomarker for cervical disease and SIL development after hrHPV infection diagnosis with implications on cervical cancer prevention strategies and treatment of SIL.

Effect of age and season on respiratory mucosal immune marker profiles.

BACKGROUND: The upper respiratory tract is continuously exposed to microorganisms and noxious elements, leading to local immune responses and the secretion of immune markers. While several studies describe immune marker profiles in respiratory mucosal samples in defined patient cohorts, mucosal immune profiles from the general population during the different seasons are lacking. Such baseline profiles are essential to understand the effect of various exposures to the mucosal immune system throughout life. OBJECTIVE: We sought to establish baseline local upper respiratory mucosal immune profiles in the general population and assess these profiles with regard to age, sex, seasonality, and basic health and lifestyle factors. METHODS: We measured the concentrations of 35 immune markers involved in a broad range of immunological processes at the mucosa in nasopharyngeal swab samples from 951 individuals, aged 0 to 86 years, from a nationwide study. RESULTS: Clustering analysis showed that immune marker profiles clearly reflected immunological functions, such as tissue regeneration and antiviral responses. Immune marker concentrations changed strongly with seasonality and age, with the most profound changes occurring in the first 25 years of life; they were also associated with sex, body mass index, smoking, mild symptoms of airway infection, and chronic asthma and hay fever. CONCLUSION: Immunological analyses of noninvasive mucosal samples provide insight into mucosal immune responses to microbial and noxious element exposure in the general population. These data provide a baseline for future studies on respiratory mucosal immune responses and for the development of mucosal immune-based diagnostics.

CD7 activation regulates cytotoxicity-driven pathology in systemic sclerosis, yielding a target for selective cell depletion.

OBJECTIVES: Cytotoxic T cells and natural killer (NK) cells are central effector cells in cancer and infections. Their effector response is regulated by activating and inhibitory receptors. The regulation of these cells in systemic autoimmune diseases such as systemic sclerosis (SSc) is less defined. METHODS: We conducted ex vivo analysis of affected skin and blood samples from 4 SSc patient cohorts (a total of 165 SSc vs 80 healthy individuals) using single-cell transcriptomics, flow cytometry and multiplex immunofluorescence staining. We further analysed the effects of costimulatory modulation in functional assays, and in a severely affected SSc patient who was treated on compassionate use with a novel anti-CD3/CD7 immunotoxin treatment. RESULTS: Here, we show that SSc-affected skin contains elevated numbers of proliferating T cells, cytotoxic T cells and NK cells. These cells selectively express the costimulatory molecule CD7 in association with cytotoxic, proinflammatory and profibrotic genes, especially in recent-onset and severe disease. We demonstrate that CD7 regulates the cytolytic activity of T cells and NK cells and that selective depletion of CD7+ cells prevents cytotoxic cell-induced fibroblast contraction and inhibits their profibrotic phenotype. Finally, anti-CD3/CD7 directed depletive treatment eliminated CD7+ skin cells and stabilised disease manifestations in a severely affected SSc patient. CONCLUSION: Together, the findings imply costimulatory molecules as key regulators of cytotoxicity-driven pathology in systemic autoimmune disease, yielding CD7 as a novel target for selective depletion of pathogenic cells.

Longitudinal analysis on the ecological dynamics of the cervicovaginal microbiome in hrHPV infection.

The cervicovaginal microbiome (CVM) is a dynamic continuous microenvironment that can be clustered in microbial community state types (CSTs) and is associated with women's cervical health. Lactobacillus-depleted communities particularly associate with an increased susceptibility for persistence of high-risk human papillomavirus (hrHPV) infections and progression of disease, but the long-term ecological dynamics of CSTs after hrHPV infection diagnosis remain poorly understood. To determine such dynamics, we examined the CVM of our longitudinal cohort of 141 women diagnosed with hrHPV infection at baseline with collected cervical smears at two timepoints six-months apart. Here we describe that the long-term microbiome dissimilarity has a positive correlation with microbial diversity at both visits and that women with high abundance and dominance for Lactobacillus iners at baseline exhibit more similar microbiome composition at second visit than women with Lactobacillus-depleted communities at baseline. We further show that the species Lactobacillus acidophilus and Megasphaera genomosp type 1 associate with CST changes between both visits. Lastly, we also observe that Gardnerella vaginalis is associated with the stability of Lactobacillus-depleted communities while L. iners is associated with the instability of Megasphaera genomosp type 1-dominated communities. Our data suggest dynamic patterns of cervicovaginal CSTs during hrHPV infection, which could be potentially used to develop microbiome-based therapies against infection progression towards disease.

In-depth insights into cervicovaginal microbial communities and hrHPV infections using high-resolution microbiome profiling.

The cervicovaginal microbiome (CVM) correlates with women's cervical health, and variations in its composition are associated with high-risk human papillomavirus (hrHPV) infection outcomes. Cervicovaginal microbes have been grouped into five community state types (CSTs) based on microbial community composition and abundance. However, studying the impact of CSTs in health and disease is challenging because the current sequencing technologies have limited confident discrimination between closely related and yet functionally different bacterial species. Circular probe-based RNA sequencing (ciRNAseq) achieves high-resolution microbiome profiling and therefore provides in-depth and unambiguous knowledge about the composition of the CVM. Based on ciRNAseq profiling of a large cohort of cervical smears (n = 541), we here define subgroups of CSTs I, III, and IV based on intra-CST differences with respect to abundances of Lactobacillus acidophilus (CSTs I-A vs. I-B and CSTs III-A vs. III-B), Lactobacillus iners (CSTs I-A vs. I-B and CSTs III-A vs. III-B), and Megasphaera genomosp type 1 (CSTs IV-A vs. IV-B). Our results further support the existence of subgroups of CST IV-C that are dominant for non-Lactobacillus species and have intermediate microbial diversity. We also show that CST V is associated with uninfected conditions, and CST IV-A associates with hrHPV-induced cervical disease. In conclusion, we characterized new subdivisions of cervicovaginal CSTs, which may further advance our understanding of women's cervical health and hrHPV-related progression to disease.

Assessing the Cervicovaginal Microbiota in the Context of hrHPV Infections: Temporal Dynamics and Therapeutic Strategies.

Cervical cancer is the third leading cause of female cancers globally, resulting in more than 300,000 deaths every year. The majority of all cervical cancers are caused by persistent infections with high-risk human papillomaviruses (hrHPV) that can progress to cancer via a series of premalignant lesions. Most women, however, clear this infection within a year, concomitant with disease regression. Both hrHPV clearance and disease regression have been associated with the composition of the cervicovaginal microenvironment, which is defined by the host immune system and the cervicovaginal microbiome (CVM). A healthy microbiome is generally characterized by a high abundance of Lactobacillus species, and a change in the composition may cause bacterial vaginosis (BV), which is associated with an increased susceptibility to persistent hrHPV infections and disease. In this review, the composition of the CVM is discussed, with emphasis on the possible causes that drive changes in the cervicovaginal microbiota in relation to hrHPV infections, disease progression, and disease regression. The literature search focused on the composition of the CVM and its correlation with hrHPV infections and neoplastic lesions as well as the current efforts to adjust the microbiome against adverse viral outcomes.

CG7630 is the Drosophila melanogaster homolog of the cytochrome c oxidase subunit COX7B.

The mitochondrial respiratory chain (MRC) is composed of four multiheteromeric enzyme complexes. According to the endosymbiotic origin of mitochondria, eukaryotic MRC derives from ancestral proteobacterial respiratory structures consisting of a minimal set of complexes formed by a few subunits associated with redox prosthetic groups. These enzymes, which are the "core" redox centers of respiration, acquired additional subunits, and increased their complexity throughout evolution. Cytochrome c oxidase (COX), the terminal component of MRC, has a highly interspecific heterogeneous composition. Mammalian COX consists of 14 different polypeptides, of which COX7B is considered the evolutionarily youngest subunit. We applied proteomic, biochemical, and genetic approaches to investigate the COX composition in the invertebrate model Drosophila melanogaster. We identified and characterized a novel subunit which is widely different in amino acid sequence, but similar in secondary and tertiary structures to COX7B, and provided evidence that this object is in fact replacing the latter subunit in virtually all protostome invertebrates. These results demonstrate that although individual structures may differ the composition of COX is functionally conserved between vertebrate and invertebrate species.

Targeted RNA next generation sequencing analysis of cervical smears can predict the presence of hrHPV-induced cervical lesions.

BACKGROUND: Because most cervical cancers are caused by high-risk human papillomaviruses (hrHPVs), cervical cancer prevention programs increasingly employ hrHPV testing as a primary test. The high sensitivity of HPV tests is accompanied by low specificity, resulting in high rates of overdiagnosis and overtreatment. Targeted circular probe-based RNA next generation sequencing (ciRNAseq) allows for the quantitative detection of RNAs of interest with high sequencing depth. Here, we examined the potential of ciRNAseq-testing on cervical scrapes to identify hrHPV-positive women at risk of having or developing high-grade cervical intraepithelial neoplasia (CIN). METHODS: We performed ciRNAseq on 610 cervical scrapes from the Dutch cervical cancer screening program to detect gene expression from 15 hrHPV genotypes and from 429 human genes. Differentially expressed hrHPV- and host genes in scrapes from women with outcome "no CIN" or "CIN2+" were identified and a model was built to distinguish these groups. RESULTS: Apart from increasing percentages of hrHPV oncogene expression from "no CIN" to high-grade cytology/histology, we identified genes involved in cell cycle regulation, tyrosine kinase signaling pathways, immune suppression, and DNA repair being expressed at significantly higher levels in scrapes with high-grade cytology and histology. Machine learning using random forest on all the expression data resulted in a model that detected 'no CIN' versus CIN2+ in an independent data set with sensitivity and specificity of respectively 85 ± 8% and 72 ± 13%. CONCLUSIONS: CiRNAseq on exfoliated cells in cervical scrapes measures hrHPV-(onco)gene expression and host gene expression in one single assay and in the process identifies HPV genotype. By combining these data and applying machine learning protocols, the risk of CIN can be calculated. Because ciRNAseq can be performed in high-throughput, making it cost-effective, it can be a promising screening technology to stratify women at risk of CIN2+. Further increasing specificity by model improvement in larger cohorts is warranted.

Novel high-resolution targeted sequencing of the cervicovaginal microbiome.

BACKGROUND: The cervicovaginal microbiome (CVM) plays a significant role in women's cervical health and disease. Microbial alterations at the species level and characteristic community state types (CST) have been associated with acquisition and persistence of high-risk human papillomavirus (hrHPV) infections that may result in progression of cervical lesions to malignancy. Current sequencing methods, especially most commonly used multiplex 16S rRNA gene sequencing, struggle to fully clarify these changes because they generally fail to provide sufficient taxonomic resolution to adequately perform species-level associative studies. To improve CVM species designation, we designed a novel sequencing tool targeting microbes at the species taxonomic rank and examined its potential for profiling the CVM. RESULTS: We introduce an accessible and practical circular probe-based RNA sequencing (CiRNAseq) technology with the potential to profile and quantify the CVM. In vitro and in silico validations demonstrate that CiRNAseq can distinctively detect species in a mock mixed microbial environment, with the output data reflecting its ability to estimate microbes' abundance. Moreover, compared to 16S rRNA gene sequencing, CiRNAseq provides equivalent results but with improved sequencing sensitivity. Analyses of a cohort of cervical smears from hrHPV-negative women versus hrHPV-positive women with high-grade cervical intraepithelial neoplasia confirmed known differences in CST occurring in the CVM of women with hrHPV-induced lesions. The technique also revealed variations in microbial diversity and abundance in the CVM of hrHPV-positive women when compared to hrHPV-negative women. CONCLUSIONS: CiRNAseq is a promising tool for studying the interplay between the CVM and hrHPV in cervical carcinogenesis. This technology could provide a better understanding of cervicovaginal CST and microbial species during health and disease, prompting the discovery of biomarkers, additional to hrHPV, that can help detect high-grade cervical lesions.

Novel insights from the Plasmodium falciparum sporozoite-specific proteome by probabilistic integration of 26 studies.

Plasmodium species, the causative agent of malaria, have a complex life cycle involving two hosts. The sporozoite life stage is characterized by an extended phase in the mosquito salivary glands followed by free movement and rapid invasion of hepatocytes in the human host. This transmission stage has been the subject of many transcriptomics and proteomics studies and is also targeted by the most advanced malaria vaccine. We applied Bayesian data integration to determine which proteins are not only present in sporozoites but are also specific to that stage. Transcriptomic and proteomic Plasmodium data sets from 26 studies were weighted for how representative they are for sporozoites, based on a carefully assembled gold standard for Plasmodium falciparum (Pf) proteins known to be present or absent during the sporozoite life stage. Of 5418 Pf genes for which expression data were available at the RNA level or at the protein level, 975 were identified as enriched in sporozoites and 90 specific to them. We show that Pf sporozoites are enriched for proteins involved in type II fatty acid synthesis in the apicoplast and GPI anchor synthesis, but otherwise appear metabolically relatively inactive in the salivary glands of mosquitos. Newly annotated hypothetical sporozoite-specific and sporozoite-enriched proteins highlight sporozoite-specific functions. They include PF3D7_0104100 that we identified to be homologous to the prominin family, which in human has been related to a quiescent state of cancer cells. We document high levels of genetic variability for sporozoite proteins, specifically for sporozoite-specific proteins that elicit antibodies in the human host. Nevertheless, we can identify nine relatively well-conserved sporozoite proteins that elicit antibodies and that together can serve as markers for previous exposure. Our understanding of sporozoite biology benefits from identifying key pathways that are enriched during this life stage. This work can guide studies of molecular mechanisms underlying sporozoite biology and potential well-conserved targets for marker and drug development.

Metabolic models predict bacterial passengers in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a complex multifactorial disease. Increasing evidence suggests that the microbiome is involved in different stages of CRC initiation and progression. Beyond specific pro-oncogenic mechanisms found in pathogens, metagenomic studies indicate the existence of a microbiome signature, where particular bacterial taxa are enriched in the metagenomes of CRC patients. Here, we investigate to what extent the abundance of bacterial taxa in CRC metagenomes can be explained by the growth advantage resulting from the presence of specific CRC metabolites in the tumor microenvironment. METHODS: We composed lists of metabolites and bacteria that are enriched on CRC samples by reviewing metabolomics experimental literature and integrating data from metagenomic case-control studies. We computationally evaluated the growth effect of CRC enriched metabolites on over 1500 genome-based metabolic models of human microbiome bacteria. We integrated the metabolomics data and the mechanistic models by using scores that quantify the response of bacterial biomass production to CRC-enriched metabolites and used these scores to rank bacteria as potential CRC passengers. RESULTS: We found that metabolic networks of bacteria that are significantly enriched in CRC metagenomic samples either depend on metabolites that are more abundant in CRC samples or specifically benefit from these metabolites for biomass production. This suggests that metabolic alterations in the cancer environment are a major component shaping the CRC microbiome. CONCLUSION: Here, we show with in sillico models that supplementing the intestinal environment with CRC metabolites specifically predicts the outgrowth of CRC-associated bacteria. We thus mechanistically explain why a range of CRC passenger bacteria are associated with CRC, enhancing our understanding of this disease. Our methods are applicable to other microbial communities, since it allows the systematic investigation of how shifts in the microbiome can be explained from changes in the metabolome.

Mapping actionable pathways and mutations in brain tumours using targeted RNA next generation sequencing.

Many biology-based precision drugs are available that neutralize aberrant molecular pathways in cancer. Molecular heterogeneity and the lack of reliable companion diagnostic biomarkers for many drugs makes targeted treatment of cancer inaccurate for many individuals. Identifying actionable hyperactive biological pathways in individual cancers may improve this situation.To achieve this we applied a novel targeted RNA next generation sequencing (t/RNA-NGS) technique to surgically obtained glioma tissues. The test combines mutation detection with analysis of biological pathway activities that are involved in tumour behavior in many cancer types (e.g. tyrosine kinase signaling, angiogenesis signaling, immune response, metabolism), via quantitative measurement of transcript levels and splice variants of hundreds of genes. We here present proof of concept that the technique, which uses molecular inversion probes, generates a histology-independent molecular diagnosis and identifies classifiers that are strongly associated with conventional histopathology diagnoses and even with patient prognosis. The test not only confirmed known glioma-associated molecular aberrations but also identified aberrant expression levels of actionable genes and mutations that have so far been considered not to be associated with glioma, opening up the possibility of drug repurposing for individual patients. Its cost-effectiveness makes t/RNA-NGS to an attractive instrument to aid oncologists in therapy decision making.

Isocitrate dehydrogenase 1-mutated human gliomas depend on lactate and glutamate to alleviate metabolic stress.

Diffuse gliomas often carry point mutations in isocitrate dehydrogenase ( IDH1mut), resulting in metabolic stress. Although IDHmut gliomas are difficult to culture in vitro, they thrive in the brain via diffuse infiltration, suggesting brain-specific tumor-stroma interactions that can compensate for IDH-1 deficits. To elucidate the metabolic adjustments in clinical IDHmut gliomas that contribute to their malignancy, we applied a recently developed method of targeted quantitative RNA next-generation sequencing to 66 clinical gliomas and relevant orthotopic glioma xenografts, with and without the endogenous IDH-1R132H mutation. Datasets were analyzed in R using Manhattan plots to calculate distance between expression profiles, Ward's method to perform unsupervised agglomerative clustering, and the Mann Whitney U test and Fisher's exact tests for supervised group analyses. The significance of transcriptome data was investigated by protein analysis, in situ enzymatic activity mapping, and in vivo magnetic resonance spectroscopy of orthotopic IDH1mut- and IDHwt-glioma xenografts. Gene set enrichment analyses of clinical IDH1mut gliomas strongly suggest a role for catabolism of lactate and the neurotransmitter glutamate, whereas, in IDHwt gliomas, processing of glucose and glutamine are the predominant metabolic pathways. Further evidence of the differential metabolic activity in these cancers comes from in situ enzymatic mapping studies and preclinical in vivo magnetic resonance spectroscopy imaging. Our data support an evolutionary model in which IDHmut glioma cells exist in symbiosis with supportive neuronal cells and astrocytes as suppliers of glutamate and lactate, possibly explaining the diffuse nature of these cancers. The dependency on glutamate and lactate opens the way for novel approaches in the treatment of IDHmut gliomas.-Lenting, K., Khurshed, M., Peeters, T. H., van den Heuvel, C. N. A. M., van Lith, S. A. M., de Bitter, T., Hendriks, W., Span, P. N., Molenaar, R. J., Botman, D., Verrijp, K., Heerschap, A., ter Laan, M., Kusters, B., van Ewijk, A., Huynen, M. A., van Noorden, C. J. F., Leenders, W. P. J. Isocitrate dehydrogenase 1-mutated human gliomas depend on lactate and glutamate to alleviate metabolic stress.

Genome-scale detection of positive selection in nine primates predicts human-virus evolutionary conflicts.

Hotspots of rapid genome evolution hold clues about human adaptation. We present a comparative analysis of nine whole-genome sequenced primates to identify high-confidence targets of positive selection. We find strong statistical evidence for positive selection in 331 protein-coding genes (3%), pinpointing 934 adaptively evolving codons (0.014%). Our new procedure is stringent and reveals substantial artefacts (20% of initial predictions) that have inflated previous estimates. The final 331 positively selected genes (PSG) are strongly enriched for innate and adaptive immunity, secreted and cell membrane proteins (e.g. pattern recognition, complement, cytokines, immune receptors, MHC, Siglecs). We also find evidence for positive selection in reproduction and chromosome segregation (e.g. centromere-associated CENPO, CENPT), apolipoproteins, smell/taste receptors and mitochondrial proteins. Focusing on the virus-host interaction, we retrieve most evolutionary conflicts known to influence antiviral activity (e.g. TRIM5, MAVS, SAMHD1, tetherin) and predict 70 novel cases through integration with virus-human interaction data. Protein structure analysis further identifies positive selection in the interaction interfaces between viruses and their cellular receptors (CD4-HIV; CD46-measles, adenoviruses; CD55-picornaviruses). Finally, primate PSG consistently show high sequence variation in human exomes, suggesting ongoing evolution. Our curated dataset of positive selection is a rich source for studying the genetics underlying human (antiviral) phenotypes. Procedures and data are available at https://github.com/robinvanderlee/positive-selection.

Interspecies differences in PTH-mediated PKA phosphorylation of the epithelial calcium channel TRPV5.

The epithelial calcium (Ca2+) channel TRPV5 (transient receptor potential vanilloid 5) is expressed in the distal convoluted tubule of the kidney and facilitates active Ca2+ reabsorption. This process is instrumental for the maintenance of Ca2+ homeostasis. Therefore, all aspects of TRPV5 function are tightly regulated by the calciotropic parathyroid hormone (PTH). Rabbit (rb)TRPV5 channel activity was shown to be stimulated upon PTH-mediated protein kinase A (PKA) phosphorylation. Since there is incomplete conservation of the PKA consensus motif (RR/QxT) across species, the aim of this study was to extend these findings to humans and characterize the expression and function of human (h)TRPV5. Functional differences between rbTRPV5 and hTRPV5 upon PTH stimulation were investigated using 45Ca2+ uptake assays, Fura-2 Ca2+ imaging, and cell surface biotinylation. While PTH treatment enhanced rbTRPV5 channel activity, it did not stimulate hTRPV5 activity. Mutation of the human RQxT motif into rabbit RRxT (hTRPV5 Q706R) partially restored the sensitivity to PTH. An ancestral sequence reconstruction of TRPV5 orthologues demonstrated that the change in the RRxT motif coincides with the creation of another putative PKA motif (RGAS to RRAS) in the amino terminus of hTRPV5. Interestingly, a constitutively phosphorylated hTRPV5 mutant (hTRPV5 S141D) displayed significantly decreased channel function, while its plasma membrane abundance was increased. Taken together, PTH-mediated stimulation of TRPV5, via PKA, is not conserved in humans. Our data suggest that PTH regulation of TRPV5 is altered in humans, an important observation for future studies that may add to new concepts on the role of PTH in renal Ca2+ handling.

Fifteen years of research on oral-facial-digital syndromes: from 1 to 16 causal genes.

Oral-facial-digital syndromes (OFDS) gather rare genetic disorders characterised by facial, oral and digital abnormalities associated with a wide range of additional features (polycystic kidney disease, cerebral malformations and several others) to delineate a growing list of OFDS subtypes. The most frequent, OFD type I, is caused by a heterozygous mutation in the OFD1 gene encoding a centrosomal protein. The wide clinical heterogeneity of OFDS suggests the involvement of other ciliary genes. For 15 years, we have aimed to identify the molecular bases of OFDS. This effort has been greatly helped by the recent development of whole-exome sequencing (WES). Here, we present all our published and unpublished results for WES in 24 cases with OFDS. We identified causal variants in five new genes (C2CD3, TMEM107, INTU, KIAA0753 and IFT57) and related the clinical spectrum of four genes in other ciliopathies (C5orf42, TMEM138, TMEM231 and WDPCP) to OFDS. Mutations were also detected in two genes previously implicated in OFDS. Functional studies revealed the involvement of centriole elongation, transition zone and intraflagellar transport defects in OFDS, thus characterising three ciliary protein modules: the complex KIAA0753-FOPNL-OFD1, a regulator of centriole elongation; the Meckel-Gruber syndrome module, a major component of the transition zone; and the CPLANE complex necessary for IFT-A assembly. OFDS now appear to be a distinct subgroup of ciliopathies with wide heterogeneity, which makes the initial classification obsolete. A clinical classification restricted to the three frequent/well-delineated subtypes could be proposed, and for patients who do not fit one of these three main subtypes, a further classification could be based on the genotype.

CTCF-mediated chromatin loops enclose inducible gene regulatory domains.

BACKGROUND: The CCTC-binding factor (CTCF) protein is involved in genome organization, including mediating three-dimensional chromatin interactions. Human patient lymphocytes with mutations in a single copy of the CTCF gene have reduced expression of enhancer-associated genes involved in response to stimuli. We hypothesize that CTCF interactions stabilize enhancer-promoter chromatin interaction domains, facilitating increased expression of genes in response to stimuli. Here we systematically investigate this model using computational analyses. RESULTS: We use CTCF ChIA-PET data from the ENCODE project to show that CTCF-associated chromatin loops have a tendency to enclose regions of enhancer-regulated stimulus responsive genes, insulating them from neighboring regions of constitutively expressed housekeeping genes. To facilitate cell type-specific CTCF loop identification, we develop an algorithm to predict CTCF loops from ChIP-seq data alone by exploiting the CTCF motif directionality in loop anchors. We apply this algorithm to a hundred ENCODE cell line datasets, confirming the universality of our observations as well as identifying a general distinction between primary and immortal cells in loop-enclosed gene content. Finally, we combine the existing evidence to propose a model for the formation of CTCF loops in which partner sites are brought together by chromatin template reeling through stationary RNA polymerases, consistent with the transcription factory hypothesis. CONCLUSIONS: We provide computational evidence that CTCF-mediated chromatin interactions enclose domains of stimulus responsive enhancer-regulated genes, insulating them from nearby housekeeping genes.

TMEM199 Deficiency Is a Disorder of Golgi Homeostasis Characterized by Elevated Aminotransferases, Alkaline Phosphatase, and Cholesterol and Abnormal Glycosylation.

Congenital disorders of glycosylation (CDGs) form a genetically and clinically heterogeneous group of diseases with aberrant protein glycosylation as a hallmark. A subgroup of CDGs can be attributed to disturbed Golgi homeostasis. However, identification of pathogenic variants is seriously complicated by the large number of proteins involved. As part of a strategy to identify human homologs of yeast proteins that are known to be involved in Golgi homeostasis, we identified uncharacterized transmembrane protein 199 (TMEM199, previously called C17orf32) as a human homolog of yeast V-ATPase assembly factor Vph2p (also known as Vma12p). Subsequently, we analyzed raw exome-sequencing data from families affected by genetically unsolved CDGs and identified four individuals with different mutations in TMEM199. The adolescent individuals presented with a mild phenotype of hepatic steatosis, elevated aminotransferases and alkaline phosphatase, and hypercholesterolemia, as well as low serum ceruloplasmin. Affected individuals showed abnormal N- and mucin-type O-glycosylation, and mass spectrometry indicated reduced incorporation of galactose and sialic acid, as seen in other Golgi homeostasis defects. Metabolic labeling of sialic acids in fibroblasts confirmed deficient Golgi glycosylation, which was restored by lentiviral transduction with wild-type TMEM199. V5-tagged TMEM199 localized with ERGIC and COPI markers in HeLa cells, and electron microscopy of a liver biopsy showed dilated organelles suggestive of the endoplasmic reticulum and Golgi apparatus. In conclusion, we have identified TMEM199 as a protein involved in Golgi homeostasis and show that TMEM199 deficiency results in a hepatic phenotype with abnormal glycosylation.

Evolution and structural organization of the mitochondrial contact site (MICOS) complex and the mitochondrial intermembrane space bridging (MIB) complex.

We have analyzed the distribution of mitochondrial contact site and cristae organizing system (MICOS) complex proteins and mitochondrial intermembrane space bridging complex (MIB) proteins over (sub)complexes and over species. The MICOS proteins are associated with the formation and maintenance of mitochondrial cristae. Indeed, the presence of MICOS genes in genomes correlates well with the presence of cristae: all cristae containing species have at least one MICOS gene and cristae-less species have none. Mic10 is the most widespread MICOS gene, while Mic60 appears be the oldest one, as it originates in the ancestors of mitochondria, the proteobacteria. In proteobacteria the gene occurs in clusters with genes involved in heme synthesis while the protein has been observed in intracellular membranes of the alphaproteobacterium Rhodobacter sphaeroides. In contrast, Mic23 and Mic27 appear to be the youngest MICOS proteins, as they only occur in opisthokonts. The remaining MICOS proteins, Mic10, Mic19, Mic25 and Mic12, the latter we show to be orthologous to human C19orf70/QIL1, trace back to the root of the eukaryotes. Of the remaining MIB proteins, also DNAJC11 shows a high correlation with the presence of cristae. In mitochondrial protein complexome profiles, the MIB complex occurs as a defined complex and as separate subcomplexes, potentially reflecting various assembly stages. We find three main forms of the complex: A) The MICOS complex, containing all the MICOS proteins, B) a membrane bridging subcomplex, containing in addition SAMM50, MTX2 and the previously uncharacterized MTX3, and C) the complete MIB complex containing in addition DNAJC11 and MTX1.

TMEM107 recruits ciliopathy proteins to subdomains of the ciliary transition zone and causes Joubert syndrome.

The transition zone (TZ) ciliary subcompartment is thought to control cilium composition and signalling by facilitating a protein diffusion barrier at the ciliary base. TZ defects cause ciliopathies such as Meckel-Gruber syndrome (MKS), nephronophthisis (NPHP) and Joubert syndrome (JBTS). However, the molecular composition and mechanisms underpinning TZ organization and barrier regulation are poorly understood. To uncover candidate TZ genes, we employed bioinformatics (coexpression and co-evolution) and identified TMEM107 as a TZ protein mutated in oral-facial-digital syndrome and JBTS patients. Mechanistic studies in Caenorhabditis elegans showed that TMEM-107 controls ciliary composition and functions redundantly with NPHP-4 to regulate cilium integrity, TZ docking and assembly of membrane to microtubule Y-link connectors. Furthermore, nematode TMEM-107 occupies an intermediate layer of the TZ-localized MKS module by organizing recruitment of the ciliopathy proteins MKS-1, TMEM-231 (JBTS20) and JBTS-14 (TMEM237). Finally, MKS module membrane proteins are immobile and super-resolution microscopy in worms and mammalian cells reveals periodic localizations within the TZ. This work expands the MKS module of ciliopathy-causing TZ proteins associated with diffusion barrier formation and provides insight into TZ subdomain architecture.