Exome variant prioritization in a large cohort of hearing-impaired individuals indicates IKZF2 to be associated with non-syndromic hearing loss and guides future research of unsolved cases.

Although more than 140 genes have been associated with non-syndromic hereditary hearing loss (HL), at least half of the cases remain unexplained in medical genetic testing. One reason is that pathogenic variants are located in 'novel' deafness genes. A variant prioritization approach was used to identify novel (candidate) genes for HL. Exome-wide sequencing data were assessed for subjects with presumed hereditary HL that remained unexplained in medical genetic testing by gene-panel analysis. Cases in group AD had presumed autosomal dominantly inherited HL (n = 124), and in group AR, presumed autosomal recessive HL (n = 337). Variants in known and candidate deafness genes were prioritized based on allele frequencies and predicted effects. Selected variants were tested for their co-segregation with HL. Two cases were solved by variants in recently identified deafness genes (ABHD12, TRRAP). Variant prioritization also revealed potentially causative variants in candidate genes associated with recessive and X-linked HL. Importantly, missense variants in IKZF2 were found to co-segregate with dominantly inherited non-syndromic HL in three families. These variants specifically affected Zn2+-coordinating cysteine or histidine residues of the zinc finger motifs 2 and 3 of the encoded protein Helios. This finding indicates a complex genotype-phenotype correlation for IKZF2 defects, as this gene was previously associated with non-syndromic dysfunction of the immune system and ICHAD syndrome, including HL. The designed strategy for variant prioritization revealed that IKZF2 variants can underlie non-syndromic HL. The large number of candidate genes for HL and variants therein stress the importance of inclusion of family members for variant prioritization.

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.

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.

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.

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.

RNA-based high-risk HPV genotyping and identification of high-risk HPV transcriptional activity in cervical tissues.

Nearly all cervical cancers are initiated by a persistent infection with one of the high-risk human papillomaviruses (high-risk HPV). High-risk HPV DNA testing is highly sensitive but cannot distinguish between active, productive infections and dormant infections or merely deposited virus. A solution for this shortcoming may be the detection of transcriptional activity of viral oncogenes instead of mere presence of high-risk HPVs. In this study, fresh-frozen cervical tissues (n = 22) were subjected to high-risk HPV DNA detection using the line probe assay and to targeted RNA next-generation sequencing using single-molecule molecular inversion probes. Targeted RNA sequencing was applied for (1) RNA-based genotyping of high-risk HPV, giving information on specific HPV-subtype (2) discrimination of E2, E6, and E7 transcripts and (3) discovery of possible non-HPV cancer biomarkers. Data were analyzed using computational biology. Targeted RNA sequencing enabled reliable genotyping of high-risk HPV subtypes and allowed quantitative detection of E2, E6, and E7 viral gene expression, thereby discriminating cervical lesions from normal cervical tissues. Moreover, targeted RNA sequencing identified possible cervical cancer biomarkers other than high-risk HPV. Interestingly, targeted RNA sequencing also provided high-quality transcription profiles from cervical scrape samples, even after 1 week of dry storage or storage in Preservcyt fixative. This proof of concept study shows that targeted RNA sequencing can be used for high-risk HPV genotyping and simultaneous detection of high-risk HPV gene activity. Future studies are warranted to investigate the potential of targeted RNA sequencing for risk assessment for the development of cervical lesions, based on molecular analysis of cervical scrapes.

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.

Molecular Profiling of Druggable Targets in Clear Cell Renal Cell Carcinoma Through Targeted RNA Sequencing.

Clear cell renal cell carcinoma (ccRCC) comprises more than 80% of all renal cancers and when metastasized leads to a 5-year survival rate of only 10%. The high rate of therapy failure and resistance development calls for reliable methods that provide information on the actionable biological pathways and predict optimal treatment protocols for individual patients. We here applied targeted RNA sequencing (t/RNA-NGS) using single molecule Molecular Inversion Probes on tumor nephrectomy samples of five ccRCC patients, comparing tumor with healthy kidney tissues. Transcriptome profiling focused on expression of genes with involvement in ccRCC biology that can be targeted with clinically available drugs. Results confirm high expression of vascular endothelial growth factor-A (VEGF-A) in tumor tissue relative to healthy-appearing kidney, in line with the angiogenic nature of ccRCC. PDGFRα and KIT, targets of the multi-kinase inhibitor sunitinib which is one of the current choices of first-line drug in metastasized ccRCC patients, were expressed at relatively low levels in tumor tissues, whereas significantly increased in normal kidney. Of all measured druggable tyrosine kinases, MET, AXL, or EGFR were expressed at higher levels in tumors than in normal kidney tissues, although intertumor differences were observed. Using cancer cell lines we show that t/RNA-NGS gene expression profiles can be used to predict in vitro sensitivity to targeted drugs. In conclusion, t/RNA-NGS analysis may provide insights into the (druggable) molecular make-up of individual renal cancers, and may guide personalized therapy of renal cell cancers.

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.

The Eukaryotic-Specific ISD11 Is a Complex-Orphan Protein with Ability to Bind the Prokaryotic IscS.

The eukaryotic protein Isd11 is a chaperone that binds and stabilizes the central component of the essential metabolic pathway responsible for formation of iron-sulfur clusters in mitochondria, the desulfurase Nfs1. Little is known about the exact role of Isd11. Here, we show that human Isd11 (ISD11) is a helical protein which exists in solution as an equilibrium between monomer, dimeric and tetrameric species when in the absence of human Nfs1 (NFS1). We also show that, surprisingly, recombinant ISD11 expressed in E. coli co-purifies with the bacterial orthologue of NFS1, IscS. Binding is weak but specific suggesting that, despite the absence of Isd11 sequences in bacteria, there is enough conservation between the two desulfurases to retain a similar mode of interaction. This knowledge may inform us on the conservation of the mode of binding of Isd11 to the desulfurase. We used evolutionary evidence to suggest Isd11 residues involved in the interaction.