Search

ObjectiveThe anti-α4β7 integrin antibody vedolizumab is administered at a fixed dose for the treatment of IBDs. This leads to a wide range of serum concentrations in patients and previous studies had suggested that highest exposure levels are associated with suboptimal clinical response. We aimed to determine the mechanisms underlying these non-linear exposure-efficacy characteristics of vedolizumab.DesignWe characterised over 500 samples from more than 300 subjects. We studied the binding of vedolizumab to T cells and investigated the functional consequences for dynamic adhesion, transmigration, gut homing and free binding sites in vivo. Employing single-cell RNA sequencing, we characterised α4β7 integrin-expressing T cell populations ‘resistant’ to vedolizumab and validated our findings in vitro and in samples from vedolizumab-treated patients with IBD. We also correlated our findings with a post-hoc analysis of the Gemini II and III studies.ResultsRegulatory T (TReg) cells exhibited a right-shifted vedolizumab binding profile compared with effector T (TEff) cells. Consistently, in a certain concentration range, the residual adhesion, transmigration, homing of and availability of functional α4β7 on TReg cells in vivo was higher than that of/on TEff cells. We identified a vedolizumab-‘resistant’ α4β7-expressing β1+PI16+ TReg cell subset with pronounced regulatory properties as the substrate for this effect. Our observations correlated with exposure-efficacy data from Gemini II and III trials.ConclusionCompletely blocking TEff cell trafficking with vedolizumab, while simultaneously permitting residual homing of powerful TReg cells in an optimal ‘therapeutic window’ based on target exposure levels might be a strategy to optimise treatment outcomes in patients with IBD.

Mycobacteria survive in macrophages despite triggering pattern recognition receptors and T cell–derived IFN-γ production. Mycobacterial cord factor trehalose-6,6-dimycolate (TDM) binds the C-type lectin receptor MINCLE and induces inflammatory gene expression. However, the impact of TDM on IFN-γ–induced macrophage activation is not known. In this study, we have investigated the cross-regulation of the mouse macrophage transcriptome by IFN-γ and by TDM or its synthetic analogue trehalose-6,6-dibehenate (TDB). As expected, IFN-γ induced genes involved in Ag presentation and antimicrobial defense. Transcriptional programs induced by TDM and TDB were highly similar but clearly distinct from the response to IFN-γ. The glycolipids enhanced expression of a subset of IFN-γ–induced genes associated with inflammation. In contrast, TDM/TDB exerted delayed inhibition of IFN-γ–induced genes, including pattern recognition receptors, MHC class II genes, and IFN-γ–induced GTPases, with antimicrobial function. TDM downregulated MHC class II cell surface expression and impaired T cell activation by peptide-pulsed macrophages. Inhibition of the IFN-γ–induced GTPase GBP1 occurred at the level of transcription by a partially MINCLE-dependent mechanism that may target IRF1 activity. Although activation of STAT1 was unaltered, deletion of Socs1 relieved inhibition of GBP1 expression by TDM. Nonnuclear Socs1 was sufficient for inhibition, suggesting a noncanonical, cytoplasmic mechanism. Taken together, unbiased analysis of transcriptional reprogramming revealed a significant degree of negative regulation of IFN-γ–induced Ag presentation and antimicrobial gene expression by the mycobacterial cord factor that may contribute to mycobacterial persistence.

Neuronal activity initiates transcriptional programs that shape long-term changes in plasticity. Although neuron subtypes differ in their plasticity response, most activity-dependent transcription factors (TFs) are broadly expressed across neuron subtypes and brain regions. Thus, how region- and neuronal subtype-specific plasticity are established on the transcriptional level remains poorly understood. We report that in young adult (i.e., 6–8 weeks old) mice, the developmental TF SOX11 is induced in neurons within 6 h either by electroconvulsive stimulation or by exploration of a novel environment. Strikingly, SOX11 induction was restricted to the dentate gyrus (DG) of the hippocampus. In the novel environment paradigm, SOX11 was observed in a subset of c-FOS expressing neurons (ca. 15%); whereas around 75% of SOX11+ DG granule neurons were c-FOS+, indicating that SOX11 was induced in an activity-dependent fashion in a subset of neurons. Environmental enrichment or virus-mediated overexpression of SOX11 enhanced the excitability of DG granule cells and downregulated the expression of different potassium channel subunits, whereas conditional Sox11/4 knock-out mice presented the opposite phenotype. We propose that Sox11 is regulated in an activity-dependent fashion, which is specific to the DG, and speculate that activity-dependent Sox11 expression may participate in the modulation of DG neuron plasticity.

Genome-wide analysis identifies 32 loci associated with breast cancer susceptibility, accounting for estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype(1-3). To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 x 10(-8)), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. We performed a genome-wide association study involving 10,146 kidney biopsy-diagnosed IgAN cases and 28,751 matched controls across 17 international cohorts. We defined 30 independent genome-wide significant risk loci jointly explaining 11% of disease risk. A total of 16 loci were novel, including TNFSF4, REL, CD28, CXCL8/PF4V1, LY86, LYN, ANXA3, TNFSF8/15, REEP3, ZMIZ1, RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The SNP-based heritability of IgAN was estimated at 23%. We observed a positive genetic correlation between IgAN and total serum IgA levels, allergy, tonsillectomy, and several infections, and a negative correlation with inflammatory bowel disease. All significant non-HLA loci shared with serum IgA levels had a concordant effect on the risk of IgAN. Moreover, IgAN loci were globally enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. The explained heritability was enriched in the regulatory elements of cells from the immune and hematopoietic systems and intestinal mucosa, providing support for the pathogenic role of extra-renal tissues. The polygenic risk of IgAN was associated with early disease onset, increased lifetime risk of kidney failure, as well as hematuria and several other traits in a phenome-wide association study of 590,515 individuals. In the comprehensive functional annotation analysis of candidate causal genes across genome-wide significant loci, we observed the convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Understanding seizure development requires an integrated knowledge of different scales of organization of epileptic networks. We developed a model of “epilepsy-in-a-dish” based on dissociated primary neuronal cells from neonatal rat hippocampus. We demonstrate how a single application of glutamate stimulated neurons to generate spontaneous synchronous spiking activity with further progression into spontaneous seizure-like events after a distinct latency period. By computational analysis, we compared the observed neuronal activity in vitro with intracranial electroencephalography (EEG) data recorded from epilepsy patients and identified strong similarities, including a related sequence of events with defined onset, progression, and termination. Next, a link between the neurophysiological changes with network composition and cellular structure down to molecular changes was established. Temporal development of epileptiform network activity correlated with increased neurite outgrowth and altered branching, increased ratio of glutamatergic over GABAergic synapses, and loss of calbindin-positive interneurons, as well as genome-wide alterations in DNA methylation. Differentially methylated genes were engaged in various cellular activities related to cellular structure, intracellular signaling, and regulation of gene expression. Our data provide evidence that a single short-term excess of glutamate is sufficient to induce a cascade of events covering different scales from molecule- to network-level, all of which jointly contribute to seizure development.

Recently, oligodendrocytes (Ol) have been attributed potential immunomodulatory effects. Yet, the exact mode of interaction with pathogenic CNS infiltrating lymphocytes remains unclear. Here, we attempt to dissect mechanisms of Ol modulation during neuroinflammation and characterize the interaction of Ol with pathogenic T cells. RNA expression analysis revealed an upregulation of immune-modulatory genes and adhesion molecules (AMs), ICAM-1 and VCAM-1, in Ol when isolated from mice undergoing experimental autoimmune encephalomyelitis (EAE). To explore whether AMs are involved in the interaction of Ol with infiltrating T cells, we performed co-culture studies on mature Ol and Th1 cells. Live cell imaging analysis showed direct interaction between both cell types. Eighty percentage of Th1 cells created contacts with Ol that lasted longer than 15 min, which may be regarded as physiologically relevant. Exposure of Ol to Th1 cells or their supernatant resulted in a significant extension of Ol processes, and upregulation of AMs as well as other immunomodulatory genes. Our observations indicate that blocking of oligodendroglial ICAM-1 can reduce the number of Th1 cells initially contacting the Ol. These results suggest that AMs may play a role in the interaction between Ol and Th1 cells. We identified Ol interacting with CD4+ cells in vivo in spinal cord tissue of EAE diseased mice indicating that our in vitro findings are of interest to further scientific research in this field. Further characterization and understanding of Ol interaction with infiltrating cells may lead to new therapeutic strategies enhancing Ol protection and remyelination potential. Oligodendrocytes regulate immune modulatory genes and adhesion molecules during autoimmune neuroinflammation Oligodendrocytes interact with Th1 cells in vitro in a physiologically relevant manner Adhesion molecules may be involved in Ol-Th1 cell interaction.

Preterm neonates are at a high risk for nephron loss under adverse clinical conditions. Renal damage potentially collides with postnatal nephrogenesis. Recent animal studies suggest that nephron loss within this vulnerable phase leads to renal damage later in life. Nephrogenic pathways are commonly reactivated after kidney injury supporting renal regeneration. We hypothesized that nephron loss during nephrogenesis affects renal development, which, in turn, impairs tissue repair after secondary injury. Neonates prior to 36 wk of gestation show an active nephrogenesis. In rats, nephrogenesis is ongoing until day 10 after birth. Mimicking the situation of severe nephron loss during nephrogenesis, male pups were uninephrectomized at day 1 of life (UNXd1). A second group of males was uninephrectomized at postnatal day 14 (UNXd14), after terminated nephrogenesis. Age-matched controls were sham operated. Three days after uninephrectomy transcriptional changes in the right kidney were analyzed by RNA-sequencing, followed by functional pathway analysis. In UNXd1, 1,182 genes were differentially regulated, but only 143 genes showed a regulation both in UNXd1 and UNXd14. The functional groups “renal development” and “kidney injury” were among the most differentially regulated groups and revealed distinctive alterations. Reduced expression of candidate genes concerning renal development ( Bmp7, Gdnf, Pdgf-B, Wt1) and injury ( nephrin, podocin, Tgf-β1) were detected. The downregulation of Bmp7 and Gdnf persisted until day 28. In UNXd14, Six2 was upregulated and Pax2 was downregulated. We conclude that nephron loss during nephrogenesis affects renal development and induces a specific regulation of genes that might hinder tissue repair after secondary kidney injury.

Th2 cells provide effector functions in type 2 immune responses to helminths and allergens. Despite knowledge about molecular mechanisms of Th2 cell differentiation, there is little information on Th2 cell heterogeneity and clonal distribution between organs. To address this, we performed combined single-cell transcriptome and TCR clonotype analysis on murine Th2 cells in mesenteric lymph nodes (MLN) and lung after infection withNippostrongylus brasiliensis(Nb) as a human hookworm infection model. We find organ-specific expression profiles, but also populations with conserved migration or effector/resident memory signatures that unexpectedly cluster with potentially regulatoryIl10posFoxp3negcells. A substantial MLN subpopulation with an interferon response signature suggests a role for interferon-signaling in Th2 differentiation or diversification. Further RNA-inferred developmental directions indicate proliferation as a hub for differentiation decisions. We also link long notedCxcr3expression in the Th2 compartment to a population ofIl4posNKT cells. Although the TCR repertoire is highly heterogeneous, we identified expanded clones and CDR3 motifs. Clonal relatedness between distant organs confirmed effective exchange of Th2 effector cells, although locally expanded clones dominated the response. These results provide new insights in Th2 cell subset diversity and clonal relatedness in distant organs.

Context CPE encodes carboxypeptidase E, an enzyme that converts proneuropeptides and propeptide hormones to bioactive forms. It is widely expressed in the endocrine and central nervous system. To date, 4 individuals from 2 families with core clinical features including morbid obesity, neurodevelopmental delay, and hypogonadotropic hypogonadism, harboring biallelic loss-of-function (LoF) CPE variants, have been reported. Objective We describe 4 affected individuals from 3 unrelated consanguineous families, 2 siblings of Syrian, 1 of Egyptian, and 1 of Pakistani descent, all harboring novel homozygous CPE LoF variants. Methods After excluding Prader-Willi syndrome (PWS), exome sequencing was performed in both Syrian siblings. The variants identified in the other 2 individuals were reported as research variants in a large-scale exome study and in the ClinVar database. Computational modeling of all possible missense alterations allowed assessing CPE tolerance to missense variants. Results All affected individuals were severely obese with neurodevelopmental delay and other endocrine anomalies. Three individuals from 2 families shared the same CPE homozygous truncating variant c.361C > T, p.(Arg121*), while the fourth carried the c.994del, p.(Ser333Alafs*22) variant. Comparison of clinical features with previously described cases and standardization according to the Human Phenotype Ontology terms indicated a recognizable clinical phenotype, which we termed Blakemore-Durmaz-Vasileiou (BDV) syndrome. Computational analysis indicated high conservation of CPE domains and intolerance to missense changes. Conclusion Biallelic truncating CPE variants are associated with BDV syndrome, a clinically recognizable monogenic recessive syndrome with childhood-onset obesity, neurodevelopmental delay, hypogonadotropic hypogonadism, and hypothyroidism. BDV syndrome resembles PWS. Our findings suggest missense variants may also be clinically relevant.

Mental disorders are highly comorbid and occur together with physical diseases, which are often considered to arise from separate pathogenic pathways. We observed in alcohol-dependent patients increased serum activity of neutral sphingomyelinase. A genetic association analysis in 456,693 volunteers found associations of haplotypes of SMPD3 coding for NSM-2 (NSM) with alcohol consumption, but also with affective state, and bone mineralisation. Functional analysis in mice showed that NSM controls alcohol consumption, affective behaviour, and their interaction by regulating hippocampal volume, cortical connectivity, and monoaminergic responses. Furthermore, NSM controlled bone–brain communication by enhancing osteocalcin signalling, which can independently supress alcohol consumption and reduce depressive behaviour. Altogether, we identified a single gene source for multiple pathways originating in the brain and bone, which interlink disorders of a mental–physical co-morbidity trias of alcohol abuse—depression/anxiety—bone disorder. Targeting NSM and osteocalcin signalling may, thus, provide a new systems approach in the treatment of a mental–physical co-morbidity trias. Deutsche Forschungsgemeinschaft (German Research Foundation) https://doi.org/10.13039/501100001659

Contains fulltext : 237689.pdf (Publisher’s version ) (Open Access) Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Transcription factor 4 (TCF4) is a crucial regulator of neurodevelopment and has been linked to the pathogenesis of autism, intellectual disability and schizophrenia. As a class I bHLH transcription factor (TF), it is assumed that TCF4 exerts its neurodevelopmental functions through dimerization with proneural class II bHLH TFs. Here, we aim to identify TF partners of TCF4 in the control of interhemispheric connectivity formation. Using a new bioinformatic strategy integrating TF expression levels and regulon activities from single cell RNA-sequencing data, we find evidence that TCF4 interacts with non-bHLH TFs and modulates their transcriptional activity in Satb2+ intercortical projection neurons. Notably, this network comprises regulators linked to the pathogenesis of neurodevelopmental disorders, e.g. FOXG1, SOX11 and BRG1. In support of the functional interaction of TCF4 with non-bHLH TFs, we find that TCF4 and SOX11 biochemically interact and cooperatively control commissure formation in vivo, and regulate the transcription of genes implicated in this process. In addition to identifying new candidate interactors of TCF4 in neurodevelopment, this study illustrates how scRNA-Seq data can be leveraged to predict TF networks in neurodevelopmental processes., Summary: Single-cell RNA sequencing identifies interactions of TCF4 with non-bHLH transcription factors linked to neurodevelopmental and neuropsychiatric disease in the regulation of interhemispheric projection neuron development.

OBJECTIVE 17q12 microdeletions containing HNF1B and intragenic variants within this gene are associated with variable developmental, endocrine, and renal anomalies, often already noted prenatally as hyperechogenic/cystic kidneys. Here, we describe prenatal and postnatal phenotypes of seven individuals with HNF1B aberrations and compare their clinical and genetic data to those of previous studies. METHODS Prenatal sequencing and postnatal chromosomal microarray analysis were performed in seven individuals with renal and/or neurodevelopmental phenotypes. We evaluated HNF1B-related clinical features from 82 studies and reclassified 192 reported intragenic HNF1B variants. RESULTS In a prenatal case, we identified a novel in-frame deletion p.(Gly239del) within the HNF1B DNA-binding domain, a mutational hot spot as demonstrated by spatial clustering analysis and high computational prediction scores. The six postnatally diagnosed individuals harbored 17q12 microdeletions. Literature screening revealed variable reporting of HNF1B-associated clinical traits. Overall, both mutation groups showed a high phenotypic heterogeneity. The reclassification of all previously reported intragenic HNF1B variants provided an up-to-date overview of the mutational spectrum. CONCLUSIONS We highlight the value of prenatal HNF1B screening in renal developmental diseases. Standardized clinical reporting and systematic classification of HNF1B variants are necessary for a more accurate risk quantification of prenatal and postnatal clinical features, improving genetic counseling and prenatal decision making.

Mutations in BCOR cause X-linked dominant and X-linked recessive forms of syndromic microphthalmia. By exome sequencing, we identified the recurrent BCOR mutation p.Pro85Leu in two brothers and their unaffected mother. While the older brother's phenotype completely fits the described phenotypic spectrum of X-linked recessive BCOR-associated Lenz microphthalmia syndrome, the younger brother showed developmental delay, microcephaly, and skeletal anomalies, but not the key feature of microphthalmia. In contrast to the previously published families, our findings demonstrate a large variability of BCOR-associated, syndromic phenotypes, indicating incomplete penetrance of p.Pro85Leu with regards to microphthalmia in males.

Background: Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hypervariable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants. Objectives: This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD. Methods: We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project. Results: The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably. Conclusions: Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction.

Simple Summary Breast cancer is the most common cancer in females worldwide. To date, many gene-environment interaction (GxE) studies have been conducted to better understand how genetic factors combine with environmental factors to influence risk. However, previous studies have not found or found only a few interactions by using SNPs which were discovered from genome-wide association studies and have been conducted, for the most part, within European populations. In this study, we focused on estrogen-related lifestyle factors that have been identified for breast cancer, including several well-established reproductive factors that are mediated by hormonal mechanisms. We aimed to examine whether there are any gene and environmental factor interactions related to estrogen exposure or metabolism using a candidate approach in Korean women. We found two interactions in this study, although they were not replicated in the independent large consortium data. These findings suggest specificity in Koreans for breast cancer risk. In this study we aim to examine gene-environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 x 10(-3)). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 x 10(-4)). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.

Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni < 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Maccari et al. identify a physiological population of highly proliferative CD38+CD45RA+, IL-10–producing TCRαβ+ T cells. They can be CD4+, CD8+, or double-negative and are controlled by FAS and CTLA4, while their survival is enhanced by mTOR and STAT3 signals., The identification and characterization of rare immune cell populations in humans can be facilitated by their growth advantage in the context of specific genetic diseases. Here, we use autoimmune lymphoproliferative syndrome to identify a population of FAS-controlled TCRαβ+ T cells. They include CD4+, CD8+, and double-negative T cells and can be defined by a CD38+CD45RA+T-BET− expression pattern. These unconventional T cells are present in healthy individuals, are generated before birth, are enriched in lymphoid tissue, and do not expand during acute viral infection. They are characterized by a unique molecular signature that is unambiguously different from other known T cell differentiation subsets and independent of CD4 or CD8 expression. Functionally, FAS-controlled T cells represent highly proliferative, noncytotoxic T cells with an IL-10 cytokine bias. Mechanistically, regulation of this physiological population is mediated by FAS and CTLA4 signaling, and its survival is enhanced by mTOR and STAT3 signals. Genetic alterations in these pathways result in expansion of FAS-controlled T cells, which can cause significant lymphoproliferative disease., Graphical Abstract

Dendritic cells (DCs) are professional antigen-presenting cells that induce and regulate adaptive immunity by presenting antigens to T cells. Due to their coordinative role in adaptive immune responses, DCs have been used as cell-based therapeutic vaccination against cancer. The capacity of DCs to induce a therapeutic immune response can be enhanced by re-wiring of cellular signalling pathways with microRNAs (miRNAs). Since the activation and maturation of DCs is controlled by an interconnected signalling network, we deploy an approach that combines RNA sequencing data and systems biology methods to delineate miRNA-based strategies that enhance DC-elicited immune responses.Through RNA sequencing of IKKβ-matured DCs that are currently being tested in a clinical trial on therapeutic anti-cancer vaccination, we identified 44 differentially expressed miRNAs. According to a network analysis, most of these miRNAs regulate targets that are linked to immune pathways, such as cytokine and interleukin signalling. We employed a network topology-oriented scoring model to rank the miRNAs, analysed their impact on immunogenic potency of DCs, and identified dozens of promising miRNA candidates with miR-15a and miR-16 as the top ones. The results of our analysis are incorporated in a database which constitutes a tool to identify DC-relevant miRNA-gene interactions with therapeutic potential (www.synmirapy.net/dc-optimization).

RNA polymerase II interacts with various other complexes and factors to ensure correct initiation, elongation, and termination of mRNA transcription. One of these proteins is SR-related CTD-associated factor 4 (SCAF4), which is important for correct usage of polyA sites for mRNA termination. Using exome sequencing and international matchmaking, we identified nine likely pathogenic germline variants in SCAF4 including two splice-site and seven truncating variants, all residing in the N-terminal two thirds of the protein. Eight of these variants occurred de novo, and one was inherited. Affected individuals demonstrated a variable neurodevelopmental disorder characterized by mild intellectual disability, seizures, behavioral abnormalities, and various skeletal and structural anomalies. Paired-end RNA sequencing on blood lymphocytes of SCAF4-deficient individuals revealed a broad deregulation of more than 9,000 genes and significant differential splicing of more than 2,900 genes, indicating an important role of SCAF4 in mRNA processing. Knockdown of the SCAF4 ortholog CG4266 in the model organism Drosophila melanogaster resulted in impaired locomotor function, learning, and short-term memory. Furthermore, we observed an increased number of active zones in larval neuromuscular junctions, representing large glutamatergic synapses. These observations indicate a role of CG4266 in nervous system development and function and support the implication of SCAF4 in neurodevelopmental phenotypes. In summary, our data show that heterozygous, likely gene-disrupting variants in SCAF4 are causative for a variable neurodevelopmental disorder associated with impaired mRNA processing.

We report a meta-analysis of breast, prostate, ovarian, and endometrial cancer genome-wide association data (effective sample size: 237,483 cases/317,006 controls). This identified 465 independent lead variants (P−8) across 192 genomic regions. Four lead variants were >1Mb from previously identified risk loci for the four cancers and an additional 23 lead variant-cancer associations were novel for one of the cancers. Bayesian models supported pleiotropic effects involving at least two cancers at 222/465 lead variants in 118/192 regions. Gene-level association analysis identified 13 shared susceptibility genes (P−6) in 13 regions not previously implicated in any of the four cancers and not uncovered by our variant-level meta-analysis. Several lead variants had opposite effects across cancers, including a cluster of such variants in the TP53 pathway. Fifty-four lead variants were associated with blood cell traits and suggested genetic overlaps with clonal hematopoiesis. Our study highlights the remarkable pervasiveness of pleiotropy across hormone-related cancers, further illuminating their shared genetic and mechanistic origins at variant- and gene-level resolution.

Intercortical connectivity is important for higher cognitive brain functions by providing the basis for integrating information from both hemispheres. We show that ablation of the neurodevelopmental disorder associated bHLH factor Tcf4 results in complete loss of forebrain commissural systems in mice. Applying a new bioinformatic strategy integrating transcription factor expression levels and regulon activities from single cell RNA-sequencing data predicted a TCF-4 interacting transcription factor network in intercortical projection neurons regulating commissure formation. This network comprises a number of regulators previously linked to the pathogenesis of intellectual disability, autism-spectrum disorders and schizophrenia, e.g. Foxg1, Sox11 and Brg1. Furthermore, we demonstrate that TCF-4 and SOX11 biochemically interact and cooperatively control commissure formation in vivo, and regulate the transcription of genes implied in this process. Our study provides a regulatory transcriptional network for the development of interhemispheric connectivity with potential pathophysiological relevance in neurodevelopmental disorders.

Familial, genome-wide association (GWAS), and sequencing studies and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWAS) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through novel application of a pleiotropy-guided conditional/conjunction false discovery rate approach for the first time in the setting of a TWAS. This identified 14 new candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 new candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were > 1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. Overlaying candidate causal risk variants identified by GWAS fine mapping onto expression prediction models for genes at known loci suggested that the association for 55% of these genes was driven by the underlying GWAS signal.SignificanceThe 22 new genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Arthritis typically involves recurrence and progressive worsening at specific predilection sites, but the checkpoints between remission and persistence remain unknown. Here, we defined the molecular and cellular mechanisms of this inflammation-mediated tissue priming. Re-exposure to inflammatory stimuli caused aggravated arthritis in rodent models. Tissue priming developed locally and independently of adaptive immunity. Repeatedly stimulated primed synovial fibroblasts (SFs) exhibited enhanced metabolic activity inducing functional changes with intensified migration, invasiveness and osteoclastogenesis. Meanwhile, human SF from patients with established arthritis displayed a similar primed phenotype. Transcriptomic and epigenomic analyses as well as genetic and pharmacological targeting demonstrated that inflammatory tissue priming relies on intracellular complement C3- and C3a receptor-activation and downstream mammalian target of rapamycin- and hypoxia-inducible factor 1α-mediated metabolic SF invigoration that prevents activation-induced senescence, enhances NLRP3 inflammasome activity, and in consequence sensitizes tissue for inflammation. Our study suggests possibilities for therapeutic intervention abrogating tissue priming without immunosuppression.

Although the bone marrow contains most hematopoietic activity during adulthood, hematopoietic stem and progenitor cells can be recovered from various extramedullary sites. Cells with hematopoietic progenitor properties have even been reported in the adult brain under steady-state conditions, but their nature and localization remain insufficiently defined. Here, we describe a heterogeneous population of myeloid progenitors in the leptomeninges of adult C57BL/6 mice. This cell pool included common myeloid, granulocyte/macrophage, and megakaryocyte/erythrocyte progenitors. Accordingly, it gave rise to all major myelo-erythroid lineages in clonogenic culture assays. Brain-associated progenitors persisted after tissue perfusion and were partially inaccessible to intravenous antibodies, suggesting their localization behind continuous blood vessel endothelium such as the blood-arachnoid barrier. Flt3Cre lineage tracing and bone marrow transplantation showed that the precursors were derived from adult hematopoietic stem cells and were most likely continuously replaced via cell trafficking. Importantly, their occurrence was tied to the immunologic state of the central nervous system (CNS) and was diminished in the context of neuroinflammation and ischemic stroke. Our findings confirm the presence of myeloid progenitors at the meningeal border of the brain and lay the foundation to unravel their possible functions in CNS surveillance and local immune cell production.

Generalized pustular psoriasis (GPP) is a severe multi-systemic inflammatory disease characterized by neutrophilic pustulosis and triggered by pro-inflammatory IL-36 cytokines in skin. While 19%–41% of affected individuals harbor bi-allelic mutations in IL36RN, the genetic cause is not known in most cases. To identify and characterize new pathways involved in the pathogenesis of GPP, we performed whole-exome sequencing in 31 individuals with GPP and demonstrated effects of mutations in MPO encoding the neutrophilic enzyme myeloperoxidase (MPO). We discovered eight MPO mutations resulting in MPO -deficiency in neutrophils and monocytes. MPO mutations, primarily those resulting in complete MPO deficiency, cumulatively associated with GPP (p = 1.85E−08; OR = 6.47). The number of mutant MPO alleles significantly differed between 82 affected individuals and >4,900 control subjects (p = 1.04E−09); this effect was stronger when including IL36RN mutations (1.48E−13) and correlated with a younger age of onset (p = 0.0018). The activity of four proteases, previously implicated as activating enzymes of IL-36 precursors, correlated with MPO deficiency. Phorbol-myristate-acetate-induced formation of neutrophil extracellular traps (NETs) was reduced in affected cells (p = 0.015), and phagocytosis assays in MPO-deficient mice and human cells revealed altered neutrophil function and impaired clearance of neutrophils by monocytes (efferocytosis) allowing prolonged neutrophil persistence in inflammatory skin. MPO mutations contribute significantly to GPP’s pathogenesis. We implicate MPO as an inflammatory modulator in humans that regulates protease activity and NET formation and modifies efferocytosis. Our findings indicate possible implications for the application of MPO inhibitors in cardiovascular diseases. MPO and affected pathways represent attractive targets for inducing resolution of inflammation in neutrophil-mediated skin diseases.

We evaluated genetics of hyperuricemia and gout, their interaction with kidney function and medication intake in chronic kidney disease (CKD) patients. Genome-wide association studies (GWAS) of urate and gout were performed in 4941 CKD patients in the German Chronic Kidney Disease (GCKD) study. Effect estimates of 26 known urate-associated population-based single nucleotide polymorphisms (SNPs) were examined. Interactions of urate-associated variants with urate-altering medications and clinical characteristics of gout were evaluated. Genome-wide significant associations with serum urate and gout were identified for known loci at SLC2A9 and ABCG2, but not for novel loci. Effects of the 26 known SNPs were of similar magnitude in CKD patients compared to population-based individuals, except for SNPs at ABCG2 that showed greater effects in CKD. Gene-medication interactions were not significant when accounting for multiple testing. Associations with gout in specific joints were significant for SLC2A9 rs12498742 in wrists and midfoot joints. Known genetic variants in SLC2A9 and ABCG2 were associated with urate and gout in a CKD cohort, with effect sizes for ABCG2 significantly greater in CKD compared to the general population. CKD patients are at high risk of gout due to reduced kidney function, diuretics intake and genetic predisposition, making treatment to target challenging.

Charcot-Marie-Tooth disease (CMT) represents a heterogeneous group of hereditary peripheral neuropathies. We previously reported a CMT locus on chromosome 19q13.3 segregating with the disease in a large Costa Rican family with axonal neuropathy and autosomal recessive pattern of inheritance (CMT2B2). We proposed a homozygous missense variant in the Mediator complex 25 (MED25) gene as causative of the disease. Nevertheless, the fact that no other CMT individuals with MED25 variants were reported to date led us to reevaluate the original family. Using exome sequencing, we now identified a homozygous nonsense variant (p.Gln517ter) in the last exon of an adjacent gene, the polynucleotide kinase 3'-phosphatase (PNKP) gene. It encodes a DNA repair protein recently associated with recessive ataxia with oculomotor apraxia type 4 (AOA4) and microcephaly, seizures, and developmental delay (MCSZ). Subsequently, five unrelated Costa Rican CMT2 subjects initially identified as being heterozygous for the same MED25 variant were found to be also compound heterozygote for PNKP. All were heterozygous for the same variant found homozygous in the large family and a second one previously associated with ataxia (p.Thr408del). Detailed clinical reassessment of the initial family and the new individuals revealed in all an adult-onset slowly progressive CMT2 associated with signs of cerebellar dysfunction such as slurred speech and oculomotor involvement, but neither microcephaly, seizures, nor developmental delay. We propose that PKNP variants are the major causative variant for the CMT2 phenotype in these individuals and that the milder clinical manifestation is due to an allelic effect.

High throughput sequencing has greatly advanced disease gene identification, especially in heterogeneous entities. Despite falling costs this is still an expensive and laborious technique, particularly when studying large cohorts. To address this problem we applied Exome Pool-Seq as an economic and fast screening technology in neurodevelopmental disorders (NDDs). Sequencing of 96 individuals can be performed in eight pools of 12 samples on less than one Illumina sequencer lane. In a pilot study with 96 cases we identified 27 variants, likely or possibly affecting function. Twenty five of these were identified in 923 established NDD genes (based on SysID database, status November 2016) (ACTB, AHDC1, ANKRD11, ATP6V1B2, ATRX, CASK, CHD8, GNAS, IFIH1, KCNQ2, KMT2A, KRAS, MAOA, MED12, MED13L, RIT1, SETD5, SIN3A, TCF4, TRAPPC11, TUBA1A, WAC, ZBTB18, ZMYND11), two in 543 (SysID) candidate genes (ZNF292, BPTF), and additionally a de novo loss-of-function variant in LRRC7, not previously implicated in NDDs. Most of them were confirmed to be de novo, but we also identified X-linked or autosomal-dominantly or autosomal-recessively inherited variants. With a detection rate of 28%, Exome Pool-Seq achieves comparable results to individual exome analyses but reduces costs by >85%. Compared with other large scale approaches using Molecular Inversion Probes (MIP) or gene panels, it allows flexible re-analysis of data. Exome Pool-Seq is thus well suited for large-scale, cost-efficient and flexible screening in characterized but heterogeneous entities like NDDs.

The role of AP-1 transcription factors in early B cell development and function is still incompletely characterized. Ubieta at al. describe the function of the Fra-2/AP-1 transcription factor as a regulator of Foxo1 and Irf4 expression in B cells. Fra-2 affects B cell proliferation and maintains their number in bone marrow., The role of AP-1 transcription factors in early B cell development and function is still incompletely characterized. Here we address the role of Fra-2 in B cell differentiation. Deletion of Fra-2 leads to impaired B cell proliferation in the bone marrow. In addition, IL-7–stimulated pro–B cell cultures revealed a reduced differentiation from large pre–B cells to small B cells and immature B cells. Gene profiling and chromatin immunoprecipitation sequencing analyses unraveled a transcriptional reduction of the transcription factors Foxo1, Irf4, Ikaros, and Aiolos in Fra-2–deficient B cells. Moreover, expression of IL7Rα and Rag 1/2, downstream targets of Irf4 and Foxo1, were also reduced in the absence of Fra-2. Pro–B cell proliferation and small pre–B cell differentiation were fully rescued by expression of Foxo1 and Irf4 in Fra-2–deficient pro–B cells. Hence, Fra-2 is a key upstream regulator of Foxo1 and Irf4 expression and influences proliferation and differentiation of B cells at multiple stages.

Patients with the Mayer–Rokitansky–Küster–Hauser syndrome (MRKH) have a congenital utero–vaginal cervical aplasia, but normal or hypoplastic adnexa and develop with normal female phenotype. Some reports mostly demonstrated regular steroid hormone levels in small MRKH cohorts including single MRKH patients with hyperandrogenemia and a clinical presentationof hirsutism and acne has also been shown. Genetically a correlation ofWNT4mutations with singular MRKH patients and hyperandrogenemia was noted. This study analyzed the hormone status of 215 MRKH patients by determining the levels of luteinizing hormone (LH), follicle-stimulating hormone (FSH), estradiol, 17-OH progesterone, testosterone, dehydroepiandrosterone sulfate (DHEAS), sex hormone–binding globulin (SHBG) and prolactin to determine the incidence of hyperandrogenemia and hyperprolactinemia in MRKH patients. Additional calculations and a ratio of free androgen index and biologically active testosterone revealed a hyperandrogenemia rate of 48.3%, hyperprolactinemia of 9.8% and combined hyperandrogenemia and hyperprolactinemia of 4.2% in MRKH patients. The rates of hirsutism, acne and especially polycystic ovary syndrome (PCOS) were in the normal range of the population and showed no correlation with hyperandrogenemia. A weekly hormone assessment over 30 days comparing 5 controls and 7 MRKH patients revealed high androgen and prolactin, but lower LH/FSH and SHBG levels with MRKH patients. The sequencing ofWNT4,WNT5A,WNT7AandWNT9Bdemonstrated no significant mutations correlating with hyperandrogenemia. Taken together, this study shows that over 52% of MRKH patients have hyperandrogenemia without clinical presentation and 14% hyperprolactinemia, which appeals for general hormone assessment and adjustments of MRKH patients.

Summary Plasmodesmata (PD) are microscopic pores connecting plant cells and enable cell-to-cell transport. Currently, little information is known about the molecular mechanisms regulating PD formation and development. To uncover components of PD development we made use of the 17 kDa movement protein (MP17) encoded by the Potato leafroll virus (PLRV). The protein is required for cell-to-cell movement of the virus and localises to complex PD. Forward genetic screening for Arabidopsis mutants with altered PD binding of MP17 revealed several mutant lines, while molecular genetics, biochemical and microscopic studies allowed further characterisation. Map-based cloning of one mutant revealed a point mutation in the choline transporter-like 1 (CHER1) protein, changing glycine247 into glutamate. Mutation in CHER1 resulted in a starch excess phenotype and stunted growth. Ultrastructure analysis of shoot apical meristems, developing and fully developed leaves showed reduced PD numbers and the absence of complex PD in fully developed leaves. This indicates that cher1 mutants are impaired in PD formation and development. Global lipid profiling revealed only slight modifications in the overall lipid composition, however, altered composition of PD-associated lipids cannot be ruled out. Thus, cher1 is devoid of complex PD in developed leaves and provides insights into the formation of complex PD at the molecular level.

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype. To identify novel loci, we performed a genome-wide association study (GWAS) including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) status and tumor grade. We identified 32 novel susceptibility loci (P-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate PIn-silico analyses showed these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 37.6% for triple-negative and 54.2% for luminal A-like disease. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

The kidneys integrate information from continuous systemic processes related to the absorption, distribution, metabolism and excretion (ADME) of metabolites. To identify underlying molecular mechanisms, we performed genome-wide association studies of the urinary concentrations of 1,172 metabolites among 1,627 patients with reduced kidney function. The 240 unique metabolite–locus associations (metabolite quantitative trait loci, mQTLs) that were identified and replicated highlight novel candidate substrates for transport proteins. The identified genes are enriched in ADME-relevant tissues and cell types, and they reveal novel candidates for biotransformation and detoxification reactions. Fine mapping of mQTLs and integration with single-cell gene expression permitted the prioritization of causal genes, functional variants and target cell types. The combination of mQTLs with genetic and health information from 450,000 UK Biobank participants illuminated metabolic mediators, and hence, novel urinary biomarkers of disease risk. This comprehensive resource of genetic targets and their substrates is informative for ADME processes in humans and is relevant to basic science, clinical medicine and pharmaceutical research. Genome-wide association analysis of 1,172 urinary metabolites identifies 240 metabolite–locus associations that when combined with UK Biobank data suggest novel metabolic mediators of disease and markers of disease risk.

Background The TUBA1A-associated tubulinopathy is clinically heterogeneous with brain malformations, microcephaly, developmental delay and epilepsy being the main clinical features. It is an autosomal dominant disorder mostly caused by de novo variants in TUBA1A. Results In three individuals with developmental delay we identified heterozygous de novo missense variants in TUBA1A using exome sequencing. While the c.1307G > A, p.(Gly436Asp) variant was novel, the two variants c.518C > T, p.(Pro173Leu) and c.641G > A, p.(Arg214His) were previously described. We compared the variable phenotype observed in these individuals with a carefully conducted review of the current literature and identified 166 individuals, 146 born and 20 fetuses with a TUBA1A variant. In 107 cases with available clinical information we standardized the reported phenotypes according to the Human Phenotype Ontology. The most commonly reported features were developmental delay (98%), anomalies of the corpus callosum (96%), microcephaly (76%) and lissencephaly (agyria-pachygyria) (70%), although reporting was incomplete in the different studies. We identified a total of 121 specific variants, including 15 recurrent ones. Missense variants cluster in the C-terminal region around the most commonly affected amino acid position Arg402 (13.3%). In a three-dimensional protein model, 38.6% of all disease-causing variants including those in the C-terminal region are predicted to affect the binding of microtubule-associated proteins or motor proteins. Genotype-phenotype analysis for recurrent variants showed an overrepresentation of certain clinical features. However, individuals with these variants are often reported in the same publication. Conclusions With 166 individuals, we present the most comprehensive phenotypic and genotypic standardized synopsis for clinical interpretation of TUBA1A variants. Despite this considerable number, a detailed genotype-phenotype characterization is limited by large inter-study variability in reporting. Electronic supplementary material The online version of this article (10.1186/s13023-019-1020-x) contains supplementary material, which is available to authorized users.

Graf, Norbert/0000-0002-2248-323X; Reis, Andre/0000-0002-6301-6363; Diets, Illja/0000-0001-7603-8898; Popp, Bernt/0000-0002-3679-1081; Thiel, Christian T/0000-0003-3817-7277; Roeleveld, Nel/0000-0002-3390-4466; Reis, AlessanRSS/0000-0001-8486-7469; Vasileiou, Georgia/0000-0002-1993-1134; Gessler, Manfred/0000-0002-7915-6045 WOS:000472571300008 PubMed ID: 30694527 Two percent of patients with Wilms tumors have a positive family history. In many of these cases the genetic cause remains unresolved. By applying germline exome sequencing in two families with two affected individuals with Wilms tumors, we identified truncating mutations in TRIM28. Subsequent mutational screening of germline and tumor DNA of 269 children affected by Wilms tumor was performed, and revealed seven additional individuals with germline truncating mutations, and one individual with a somatic truncating mutation in TRIM28. TRIM28 encodes a complex scaffold protein involved in many different processes, including gene silencing, DNA repair and maintenance of genomic integrity. Expression studies on mRNA and protein level showed reduction of TRIM28, confirming a loss-of-function effect of the mutations identified. The tumors showed an epithelial-type histology that stained negative for TRIM28 by immunohistochemistry. The tumors were bilateral in six patients, and 10/11 tumors are accompanied by perilobar nephrogenic rests. Exome sequencing on eight tumor DNA samples from six individuals showed loss-of-heterozygosity (LOH) of the TRIM28-locus by mitotic recombination in seven tumors, suggesting that TRIM28 functions as a tumor suppressor gene in Wilms tumor development. Additionally, the tumors showed very few mutations in known Wilms tumor driver genes, suggesting that loss of TRIM28 is the main driver of tumorigenesis. In conclusion, we identified heterozygous germline truncating mutations in TRIM28 in 11 children with mainly epithelial-type Wilms tumors, which become homozygous in tumor tissue. These data establish TRIM28 as a novel Wilms tumor predisposition gene, acting as a tumor suppressor gene by LOH. Dutch Cancer SocietyKWF Kankerbestrijding [KUN2012-5366]; KiKa Foundation [127] Grant sponsor: Dutch Cancer Society; Grant number: KUN2012-5366; Grant sponsor: The KiKa Foundation (project 127)

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439

Macrophages are considered to contribute to chronic inflammatory diseases such as rheumatoid arthritis1. However, both the exact origin and the role of macrophages in inflammatory joint disease remain unclear. Here we use fate-mapping approaches in conjunction with three-dimensional light-sheet fluorescence microscopy and single-cell RNA sequencing to perform a comprehensive spatiotemporal analysis of the composition, origin and differentiation of subsets of macrophages within healthy and inflamed joints, and study the roles of these macrophages during arthritis. We find that dynamic membrane-like structures, consisting of a distinct population of CX3CR1+ tissue-resident macrophages, form an internal immunological barrier at the synovial lining and physically seclude the joint. These barrier-forming macrophages display features that are otherwise typical of epithelial cells, and maintain their numbers through a pool of locally proliferating CX3CR1− mononuclear cells that are embedded into the synovial tissue. Unlike recruited monocyte-derived macrophages, which actively contribute to joint inflammation, these epithelial-like CX3CR1+ lining macrophages restrict the inflammatory reaction by providing a tight-junction-mediated shield for intra-articular structures. Our data reveal an unexpected functional diversification among synovial macrophages and have important implications for the general role of macrophages in health and disease. Analysis of macrophage subsets within joints reveals a population of CX3CR1+ tissue-resident macrophages that form a tight-junction-mediated barrier at the synovial lining, protecting the joint from the invasion of inflammatory cells.

Arginase 1 (Arg1), which converts l-arginine into ornithine and urea, exerts pleiotropic immunoregulatory effects. However, the function of Arg1 in inflammatory bowel disease (IBD) remains poorly characterized. Here, we found that Arg1 expression correlated with the degree of inflammation in intestinal tissues from IBD patients. In mice, Arg1 was upregulated in an IL-4/IL-13(–) and intestinal microbiota–dependent manner. Tie2-Cre Arg1(fl/fl) mice lacking Arg1 in hematopoietic and endothelial cells recovered faster from colitis than Arg1-expressing (Arg1(fl/fl)) littermates. This correlated with decreased vessel density, compositional changes in intestinal microbiota, diminished infiltration by myeloid cells, and an accumulation of intraluminal polyamines that promote epithelial healing. The proresolving effect of Arg1 deletion was reduced by an l-arginine–free diet, but rescued by simultaneous deletion of other l-arginine–metabolizing enzymes, such as Arg2 or Nos2, demonstrating that protection from colitis requires l-arginine. Fecal microbiota transfers from Tie2-Cre Arg1(fl/fl) mice into WT recipients ameliorated intestinal inflammation, while transfers from WT littermates into Arg1-deficient mice prevented an advanced recovery from colitis. Thus, an increased availability of l-arginine as well as altered intestinal microbiota and metabolic products accounts for the accelerated resolution from colitis in the absence of Arg1. Consequently, l-arginine metabolism may serve as a target for clinical intervention in IBD patients.

Biallelic variants in the AEBP1 gene cause a novel autosomal-recessive connective tissue disorder (CTD) reminiscent of Ehlers-Danlos Syndrome (EDS). The four previously reported individuals show considerable clinical variability. Unbiased high-throughput sequencing enables the rapid identification of additional cases for such rare entities. We identified the homozygous nonsense variant c.917dup, p.Tyr306* in AEBP1 using clinical exome sequencing in a female individual with previously unsolved CTD. Segregation testing confirmed homozygosity in the clinically affected brother and heterozygous carrier status in the healthy mother. Chromosomal microarray showed that the variant lies in a run of homozygosity, suggesting a common origin of this genomic segment. RT-PCR analysis in the mother revealed a monoallelic expression of the normal transcript supporting a nonsense-mediated mRNA decay and functional nullizygosity as disease mechanism. We describe two individuals from a fourth family with AEBP1-associated CTD. Our results further verify that autosomal-recessive inherited LOF variants in the AEBP1 gene cause clinical features of different EDS subtypes, but also of the marfanoid spectrum. As identification of further individuals is necessary to inform the clinical characterization, we stress the added value of exome sequencing for such rare diseases.

Genetic integrity of induced pluripotent stem cells (iPSCs) is essential for their validity as disease models and for potential therapeutic use. We describe the comprehensive analysis in the ForIPS consortium: an iPSC collection from donors with neurological diseases and healthy controls. Characterization included pluripotency confirmation, fingerprinting, conventional and molecular karyotyping in all lines. In the majority, somatic copy number variants (CNVs) were identified. A subset with available matched donor DNA was selected for comparative exome sequencing. We identified single nucleotide variants (SNVs) at different allelic frequencies in each clone with high variability in mutational load. Low frequencies of variants in parental fibroblasts highlight the importance of germline samples. Somatic variant number was independent from reprogramming, cell type and passage. Comparison with disease genes and prediction scores suggest biological relevance for some variants. We show that high-throughput sequencing has value beyond SNV detection and the requirement to individually evaluate each clone.

Whereas innate immune cells, such as NK and innate lymphoid cells (ILCs), have been characterized in different human tissues, knowledge on the thymic CD56-expressing cell subsets is limited. In this study, the rare subpopulations of thymic CD56+CD3− cells from samples of >100 patients have been successfully analyzed. The results revealed fundamental differences between thymic and peripheral blood (PB) CD56+CD3− cells. Thymic tissues lacked immunoregulatory CD56highCD16dim NK cells but showed two Eomes+CD56dim subsets on which common NK cell markers were significantly altered. CD56dimCD16high cells expressed high amounts of NKG2A, NKG2D, and CD27 with low CD57. Conversely, CD56dimCD16dim cells displayed high CD127 but low expression of KIR, NKG2D, and natural cytotoxicity receptors (NCRs). Thymic CD56+CD3− cells were able to gain cytotoxicity but were especially immunoregulatory cells, producing a broad range of cytokines. Finally, one population of thymic CD56+ cells resembled conventional NK cells, whereas the other represented a novel, noncanonical NK subset.

Although defects in intestinal barrier function are a key pathogenic factor in patients with inflammatory bowel diseases (IBDs), the molecular pathways driving disease-specific alterations of intestinal epithelial cells (IECs) are largely unknown. Here, we addressed this issue by characterizing the transcriptome of IECs from IBD patients using a genome-wide approach. We observed disease-specific alterations in IECs with markedly impaired Rho-A signaling in active IBD patients. Localization of epithelial Rho-A was shifted to the cytosol in IBDs, and inflammation was associated with suppressed Rho-A activation due to reduced expression of the Rho-A prenylation enzyme geranylgeranyltransferase-I (GGTase-l). Functionally, we found that mice with conditional loss of Rhoa or the gene encoding GGTase-1, Pggtlb, in IECs exhibit spontaneous chronic intestinal inflammation with accumulation of granulocytes and CD4(+) T cells. This phenotype was associated with cytoskeleton rearrangement and aberrant cell shedding, ultimately leading to loss of epithelial integrity and subsequent inflammation. These findings uncover deficient prenylation of Rho-A as a key player in the pathogenesis of IBDs. As therapeutic triggering of Rho-A signaling suppressed intestinal inflammation in mice with GGTase-l-deficient IECs, our findings suggest new avenues for treatment of epithelial injury and mucosal inflammation in IBD patients.

Height is a heritable and highly heterogeneous trait. Short stature affects 3% of the population and in most cases is genetic in origin. After excluding known causes, 67% of affected individuals remain without diagnosis. To identify novel candidate genes for short stature, we performed exome sequencing in 254 unrelated families with short stature of unknown cause and identified variants in 63 candidate genes in 92 (36%) independent families. Based on systematic characterization of variants and functional analysis including expression in chondrocytes, we classified 13 genes as strong candidates. Whereas variants in at least two families were detected for all 13 candidates, two genes had variants in 6 (UBR4) and 8 (LAMA5) families, respectively. To facilitate their characterization, we established a clustered network of 1025 known growth and short stature genes, which yielded 29 significantly enriched clusters, including skeletal system development, appendage development, metabolic processes, and ciliopathy. Eleven of the candidate genes mapped to 21 of these clusters, including CPZ, EDEM3, FBRS, IFT81, KCND1, PLXNA3, RASA3, SLC7A8, UBR4, USP45, and ZFHX3. Fifty additional growth-related candidates we identified await confirmation in other affected families. Our study identifies Mendelian forms of growth retardation as an important component of idiopathic short stature.

BackgroundThe TUBA1A-associated tubulinopathy is clinically heterogeneous with brain malformations, microcephaly, developmental delay and epilepsy being the main clinical features. It is an autosomal dominant disorder mostly caused by de novo variants in TUBA1A.ResultsIn three individuals with developmental delay we identified heterozygous de novo missense variants in TUBA1A using exome sequencing. While the c.1307G>A, p.(Gly436Asp) variant was novel, the two variants c.518C>T, p.(Pro173Leu) and c.641G>A, p.(Arg214His) were previously described. We compared the variable phenotype observed in these individuals with a carefully conducted review of the current literature and identified 166 individuals, 146 born and 20 fetuses with a TUBA1A variant. In 107 cases with available clinical information we standardized the reported phenotypes according to the Human Phenotype Ontology. The most commonly reported features were developmental delay (98%), anomalies of the corpus callosum (96%), microcephaly (76%) and lissencephaly (70%), although reporting was incomplete in the different studies. We identified a total of 121 distinct variants, including 15 recurrent ones. Missense variants cluster in the C-terminal region around the most commonly affected amino acid position Arg402 (13.3%). In a three-dimensional protein modelling, 38.6% of all disease causing variants including those in the C-terminal region are predicted to affect binding of microtubule-associated proteins or motor proteins. Genotype-phenotype analysis for recurrent variants showed an overrepresentation of certain clinical features. However, individuals with these variants are often reported in the same publication.ConclusionsWith 166 individuals, we present the most comprehensive phenotypic and genotypic standardized synopsis for clinical interpretation of TUBA1A variants. Despite this considerable number, a detailed genotype-phenotype characterization is limited by large inter-study variability in reporting.

Chromosomal rearrangements and specific aneuploidy patterns are initiating events and define subgroups in B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Here we analyzed 250 BCP-ALL cases and identified a novel subgroup (‘PAX5-plus’, n = 19) by distinct DNA methylation and gene expression profiles. All patients in this subgroup harbored mutations in the B-lineage transcription factor PAX5, with p.P80R as hotspot. Mutations either affected two independent codons, consistent with compound heterozygosity, or suffered LOH predominantly through chromosome 9p aberrations. These biallelic events resulted in disruption of PAX5 transcriptional programs regulating B-cell differentiation and tumor suppressor functions. Homozygous CDKN2A/B deletions and RAS-activating hotspot mutations were highly enriched as cooperating events in the genomic profile of PAX5-plus ALL. Together, this defined a specific pattern of triple alterations, exclusive to the novel subgroup. PAX5-plus ALL was observed in pediatric and adult patients. Although restricted by the limited sample size, a tendency for more favorable clinical outcome was observed, with 10 of 12 adult PAX5-plus patients achieving long-term survival. PAX5-plus represents the first BCP-ALL subgroup defined by sequence alterations in contrast to gross chromosomal events and exemplifies how deregulated differentiation (PAX5), impaired cell cycle control (CDKN2A/B) and sustained proliferative signaling (RAS) cooperatively drive leukemogenesis.

Mesenchymal stem cells (MSCs) represent key contributors to tissue homeostasis and promising therapeutics for hyperinflammatory conditions including graft-versus-host disease. Their immunomodulatory effects are controlled by microenvironmental signals. The MSCs' functional response towards inflammatory cues is known as MSC-"licensing" and includes indoleamine 2,3-dioxygenase (IDO) upregulation. MSCs use tryptophan-depleting IDO to suppress T-cells. Increasing evidence suggests that several functions are (co-)determined by the cells' metabolic commitment. MSCs are capable of both, high levels of glycolysis and of oxidative phosphorylation. Although several studies have addressed alterations of the immune regulatory phenotype elicited by inflammatory priming metabolic mechanisms controlling this process remain unknown. We demonstrate that inflammatory MSC-licensing causes metabolic shifts including enhanced glycolysis and increased fatty acid oxidation. Yet, only interfering with glycolysis impacts IDO upregulation and impedes T-cell-suppressivity. We identified the Janus kinase (JAK)/signal transducer and activator of transcription (STAT)1 pathway as a regulator of both glycolysis and IDO, and show that enhanced glucose turnover is linked to abundant STAT1 glycosylation. Inhibiting the responsible O-acetylglucosamine (O-GlcNAc) transferase abolishes STAT1 activity together with IDO upregulation. Our data suggest that STAT1-O-GlcNAcylation increases its stability towards degradation thus sustaining downstream effects. This pathway could represent a target for interventions aiming to enhance the MSCs' immunoregulatory potency.

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3′ gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09, 95% CI = 0.97–1.22, p = 0.15, N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13, 95% CI = 1.03–1.24, p = 0.01, N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10−28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.