Single-cell landscape of idiopathic multicentric Castleman disease in identical twins.

ABSTRACT: Idiopathic multicentric Castleman disease (iMCD) is a rare cytokine-driven disorder characterized by systemic inflammation, generalized lymphadenopathy, and organ dysfunction. Here, we present an unusual occurrence of iMCD in identical twins and examined the immune milieu within the affected lymphoid organs and the host circulation using multiomic high-dimensional profiling. Using spatial enhanced resolution omics sequencing (Stereo-seq) transcriptomic profiling, we performed unsupervised spatially constrained clustering to identify different anatomic structures, mapping the follicles and interfollicular regions. After a cell segmentation approach, interleukin 6 (IL-6) pathway genes significantly colocalized with endothelial cells and fibroblastic reticular cells, confirming observations using a single-cell sequencing approach (10× Chromium). Furthermore, single-cell sequencing of peripheral blood mononuclear cells revealed an "inflammatory" peripheral monocytosis enriched for the expression of S100A family genes in both twins. In summary, we provided evidence of the putative cell-of-origin of IL-6 signals in iMCD and described a distinct monocytic host immune response phenotype through a unique identical twin model.

Predicting synthetic lethal interactions using heterogeneous data sources.

MOTIVATION: A synthetic lethal (SL) interaction is a relationship between two functional entities where the loss of either one of the entities is viable but the loss of both entities is lethal to the cell. Such pairs can be used as drug targets in targeted anticancer therapies, and so, many methods have been developed to identify potential candidate SL pairs. However, these methods use only a subset of available data from multiple platforms, at genomic, epigenomic and transcriptomic levels; and hence are limited in their ability to learn from complex associations in heterogeneous data sources. RESULTS: In this article, we develop techniques that can seamlessly integrate multiple heterogeneous data sources to predict SL interactions. Our approach obtains latent representations by collective matrix factorization-based techniques, which in turn are used for prediction through matrix completion. Our experiments, on a variety of biological datasets, illustrate the efficacy and versatility of our approach, that outperforms state-of-the-art methods for predicting SL interactions and can be used with heterogeneous data sources with minimal feature engineering. AVAILABILITY AND IMPLEMENTATION: Software available at https://github.com/lianyh. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MultiDCoX: Multi-factor analysis of differential co-expression.

BACKGROUND: Differential co-expression (DCX) signifies change in degree of co-expression of a set of genes among different biological conditions. It has been used to identify differential co-expression networks or interactomes. Many algorithms have been developed for single-factor differential co-expression analysis and applied in a variety of studies. However, in many studies, the samples are characterized by multiple factors such as genetic markers, clinical variables and treatments. No algorithm or methodology is available for multi-factor analysis of differential co-expression. RESULTS: We developed a novel formulation and a computationally efficient greedy search algorithm called MultiDCoX to perform multi-factor differential co-expression analysis. Simulated data analysis demonstrates that the algorithm can effectively elicit differentially co-expressed (DCX) gene sets and quantify the influence of each factor on co-expression. MultiDCoX analysis of a breast cancer dataset identified interesting biologically meaningful differentially co-expressed (DCX) gene sets along with genetic and clinical factors that influenced the respective differential co-expression. CONCLUSIONS: MultiDCoX is a space and time efficient procedure to identify differentially co-expressed gene sets and successfully identify influence of individual factors on differential co-expression.

Analysis of non-synonymous-coding variants of Parkinson's disease-related pathogenic and susceptibility genes in East Asian populations.

To evaluate the contribution of non-synonymous-coding variants of known familial and genome-wide association studies (GWAS)-linked genes for Parkinson's disease (PD) to PD risk in the East Asian population, we sequenced all the coding exons of 39 PD-related disease genes and evaluated the accumulation of rare non-synonymous-coding variants in 375 early-onset PD cases and 399 controls. We also genotyped 782 non-synonymous-coding variants of these genes in 710 late-onset PD cases and 9046 population controls. Significant enrichment of LRRK2 variants was observed in both early- and late-onset PD (odds ratio = 1.58; 95% confidence interval = 1.29-1.93; P = 8.05 × 10(-6)). Moderate enrichment was also observed in FGF20, MCCC1, GBA and ITGA8. Half of the rare variants anticipated to cause loss of function of these genes were present in healthy controls. Overall, non-synonymous-coding variants of known familial and GWAS-linked genes appear to make a limited contribution to PD risk, suggesting that clinical sequencing of these genes will provide limited information for risk prediction and molecular diagnosis.

Linking a genome-wide association study signal to a LRRK2 coding variant in Parkinson's disease.

BACKGROUND: Genome-wide association studies have identified several loci associated with Parkinson's disease (PD). Whole-exome sequencing detects rare coding variants, but their links with PD genome-wide association study loci are unknown. Our objective was to investigate whether nonsynonymous variants in LRRK2 can explain associations at the PD-associated locus tagged by rs1994090. METHODS: We sequenced all coding exons of LRRK2 in 453 East Asian samples and evaluated linkage disequilibrium between each nonsynonymous variant and rs1994090. We then tested selected variants and haplotypes for association with PD in 13,581 East Asian samples. RESULTS: Of all the nonsynonymous variants, only p.Gly2385Arg was in moderate linkage disequilibrium with rs1994090 and was observed on haplotypes tagged by the rs1994090-C risk allele. Conditional analyses showed that associations at these 2 variants are not independent. CONCLUSIONS: LRRK2 p.Gly2385Arg can explain most if not all of the PD association at rs1994090 in East Asians, but other nonsynonymous variants are independent. © 2015 International Parkinson and Movement Disorder Society.

Targeted next-generation sequencing to diagnose disorders of HDL cholesterol.

A low level of HDL cholesterol (HDL-C) is a common clinical scenario and an important marker for increased cardiovascular risk. Many patients with very low or very high HDL-C have a rare mutation in one of several genes, but identification of the molecular abnormality in patients with extreme HDL-C is rarely performed in clinical practice. We investigated the accuracy and diagnostic yield of a targeted next-generation sequencing (NGS) assay for extreme levels of HDL-C. We developed a targeted NGS panel to capture the exons, intron/exon boundaries, and untranslated regions of 26 genes with highly penetrant effects on plasma lipid levels. We sequenced 141 patients with extreme HDL-C levels and prioritized variants in accordance with medical genetics guidelines. We identified 35 pathogenic and probably pathogenic variants in HDL genes, including 21 novel variants, and performed functional validation on a subset of these. Overall, a molecular diagnosis was established in 35.9% of patients with low HDL-C and 5.2% with high HDL-C, and all prioritized variants identified by NGS were confirmed by Sanger sequencing. Our results suggest that a molecular diagnosis can be identified in a substantial proportion of patients with low HDL-C using targeted NGS.

An association study of TOLL and CARD with leprosy susceptibility in Chinese population.

Previous genome-wide association studies (GWASs) identified multiple susceptibility loci that have highlighted the important role of TLR (Toll-like receptor) and CARD (caspase recruitment domain) genes in leprosy. A large three-stage candidate gene-based association study of 30 TLR and 47 CARD genes was performed in the leprosy samples of Chinese Han. Of 4363 SNPs investigated, eight SNPs showed suggestive association (P < 0.01) in our previously published GWAS datasets (Stage 1). Of the eight SNPs, rs2735591 and rs4889841 showed significant association (P < 0.001) in an independent series of 1504 cases and 1502 controls (Stage 2), but only rs2735591 (next to BCL10) showed significant association in the second independent series of 938 cases and 5827 controls (Stage 3). Rs2735591 showed consistent association across the three stages (P > 0.05 for heterogeneity test), significant association in the combined validation samples (Pcorrected = 5.54 × 10(-4) after correction for 4363 SNPs tested) and genome-wide significance in the whole GWAS and validation samples (P = 1.03 × 10(-9), OR = 1.24). In addition, we demonstrated the lower expression of BCL10 in leprosy lesions than normal skins and a significant gene connection between BCL10 and the eight previously identified leprosy loci that are associated with NFκB, a major regulator of downstream inflammatory responses, which provides further biological evidence for the association. We have discovered a novel susceptibility locus on 1p22, which implicates BCL10 as a new susceptibility gene for leprosy. Our finding highlights the important role of both innate and adaptive immune responses in leprosy.

Impaired development of neural-crest cell-derived organs and intellectual disability caused by MED13L haploinsufficiency.

MED13L is a component subunit of the Mediator complex, an important regulator of transcription that is highly conserved across eukaryotes. Here, we report MED13L disruption in a translocation t(12;19) breakpoint of a patient with Pierre-Robin syndrome, moderate intellectual disability, craniofacial anomalies, and muscular defects. The phenotype is similar to previously described patients with MED13L haploinsufficiency. Knockdown of MED13L orthologue in zebrafish, med13b, showed early defective migration of cranial neural crest cells (NCCs) that contributed to cartilage structure deformities in the later stage, recapitulating craniofacial anomalies seen in human patients. Notably, we observed abnormal distribution of developing neurons in different brain regions of med13b morphant embryos, which could be rescued upon introduction of full-length human MED13L mRNA. To compare with mammalian system, we suppressed MED13L expression by short-hairpin RNA in ES-derived human neural progenitors, and differentiated them into neurons. Transcriptome analysis revealed differential expression of components of Wnt and FGF signaling pathways in MED13L-deficient neurons. Our finding provides a novel insight into the mechanism of overlapping phenotypic outcome targeting NCCs derivatives organs in patients with MED13L haploinsufficiency, and emphasizes a clinically recognizable syndromic phenotype in these patients.

ASTER: A Method to Predict Clinically Relevant Synthetic Lethal Genetic Interactions.

A Synthetic Lethal (SL) interaction is a functional relationship between two genes or functional entities where the loss of either entity is viable but the loss of both is lethal. Such pairs can be used to develop targeted anticancer therapies with fewer side effects and reduced overtreatment. However, finding clinically relevant SL interactions remains challenging. Leveraging unified gene expression data of both disease-free and cancerous samples, we design a new technique based on statistical hypothesis testing, called ASTER, to identify SL pairs. We empirically find that the patterns of mutually exclusivity ASTER finds using genomic and transcriptomic data provides a strong signal of synthetic lethality. For large-scale multiple hypothesis testing, we develop an extension called ASTER++ that can utilize additional input gene features within the hypothesis testing framework. Our computational and functional experiments demonstrate the efficacy of ASTER in identifying SL pairs with potential therapeutic benefits.

An Algorithm to Mine Therapeutic Motifs for Cancer From Networks of Genetic Interactions.

Study of pairwise genetic interactions, such as mutually exclusive mutations, has led to understanding of underlying mechanisms in cancer. Investigation of various combinatorial motifs within networks of such interactions can lead to deeper insights into its mutational landscape and inform therapy development. One such motif called the Between-Pathway Model (BPM) represents redundant or compensatory pathways that can be therapeutically exploited. Finding such BPM motifs is challenging since most formulations require solving variants of the NP-complete maximum weight bipartite subgraph problem. In this paper we design an algorithm based on Integer Linear Programming (ILP) to solve this problem. In our experiments, our approach outperforms the best previous method to mine BPM motifs. Further, our ILP-based approach allows us to easily model additional application-specific constraints. We illustrate this advantage through a new application of BPM motifs that can potentially aid in finding combination therapies to combat cancer.

Integrative Analysis Identifies Key Molecular Signatures Underlying Neurodevelopmental Deficits in Fragile X Syndrome.

BACKGROUND: Fragile X syndrome (FXS) is a neurodevelopmental disorder caused by epigenetic silencing of FMR1 and loss of FMRP expression. Efforts to understand the molecular underpinnings of the disease have been largely performed in rodent or nonisogenic settings. A detailed examination of the impact of FMRP loss on cellular processes and neuronal properties in the context of isogenic human neurons remains lacking. METHODS: Using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 to introduce indels in exon 3 of FMR1, we generated an isogenic human pluripotent stem cell model of FXS that shows complete loss of FMRP expression. We generated neuronal cultures and performed genome-wide transcriptome and proteome profiling followed by functional validation of key dysregulated processes. We further analyzed neurodevelopmental and neuronal properties, including neurite length and neuronal activity, using multielectrode arrays and patch clamp electrophysiology. RESULTS: We showed that the transcriptome and proteome profiles of isogenic FMRP-deficient neurons demonstrate perturbations in synaptic transmission, neuron differentiation, cell proliferation and ion transmembrane transporter activity pathways, and autism spectrum disorder-associated gene sets. We uncovered key deficits in FMRP-deficient cells demonstrating abnormal neural rosette formation and neural progenitor cell proliferation. We further showed that FMRP-deficient neurons exhibit a number of additional phenotypic abnormalities, including neurite outgrowth and branching deficits and impaired electrophysiological network activity. These FMRP-deficient related impairments have also been validated in additional FXS patient-derived human-induced pluripotent stem cell neural cells. CONCLUSIONS: Using isogenic human pluripotent stem cells as a model to investigate the pathophysiology of FXS in human neurons, we reveal key neural abnormalities arising from the loss of FMRP.

Evaluation of novel Parkinson's disease candidate genes in the Chinese population.

Recent whole-exome sequencing studies in European patients with Parkinson's disease (PD) have identified potential risk variants across 33 novel PD candidate genes. We aim to determine if these reported candidate genes are similarly implicated in Asians by assessing common, rare, and novel nonsynonymous coding variants by sequencing all 33 genes in 198 Chinese samples and genotyping coding variants in an independent set of 9756 Chinese samples. We carried out further targeted sequencing of CD36 in an additional 576 Chinese and Korean samples. We found that only 8 of 43 reported risk variants were polymorphic in our Chinese samples. We identified several heterozygotes for rare loss-of-function mutations, including the reported CD36 p.Gln74Ter variant, in both cases and controls. We also observed 2 potential compound heterozygotes among PD cases for rare loss-of-function mutations in CD36 and SSPO. The other reported variants were common in East Asians and not associated with PD, completely absent, or only found in controls. Therefore, the 33 reported candidate genes and associated variants are unlikely to confer significant PD risk in the East Asian population.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

A Synonymous Variant in IL10RA Affects RNA Splicing in Paediatric Patients with Refractory Inflammatory Bowel Disease.

Interleukin-10 receptor [IL10R] mutations are associated with severe childhood inflammatory bowel disease [IBD]. Two unrelated patients who died of very early-onset severe IBD and sepsis were identified as harbouring the same compound heterozygous mutations in IL10RA [p.R101W; p.T179T]. A third patient was found to be homozygous for p.T179T. The missense change of p.R101W has been reported. The synonymous change of p.T179T, with a minor allele frequency of 0.035% in the population, was novel. The p.T179T mutation was located before the 5' splice donor site, leading to exon skipping and out-of-frame fusion of exons 3 and 5, causing altered STAT3 phosphorylation in IL10-induced peripheral blood mononuclear cells. The patient developed colitis at 6 years of age, the oldest reported age of onset among patients with IL10RA mutations, and did not suffer from perianal disease. We report three paediatric patients with a rare, synonymous p.T179T variant causing a splicing error in IL10RA.

Discovery of six new susceptibility loci and analysis of pleiotropic effects in leprosy.

Genome-wide association studies (GWAS) have led to the discovery of several susceptibility loci for leprosy with robust evidence, providing biological insight into the role of host genetic factors in mycobacterial infection. However, the identified loci only partially explain disease heritability, and additional genetic risk factors remain to be discovered. We performed a 3-stage GWAS of leprosy in the Chinese population using 8,313 cases and 16,017 controls. Besides confirming all previously published loci, we discovered six new susceptibility loci, and further gene prioritization analysis of these loci implicated BATF3, CCDC88B and CIITA-SOCS1 as new susceptibility genes for leprosy. A systematic evaluation of pleiotropic effects demonstrated a high tendency for leprosy susceptibility loci to show association with autoimmunity and inflammatory diseases. Further analysis suggests that molecular sensing of infection might have a similar pathogenic role across these diseases, whereas immune responses have discordant roles in infectious and inflammatory diseases.

Discovery of a novel genetic susceptibility locus on X chromosome for systemic lupus erythematosus.

INTRODUCTION: Systemic lupus erythematosus (SLE) is an autoimmune connective tissue disease affecting predominantly females. To discover additional genetic risk variants for SLE on the X chromosome, we performed a follow-up study of our previously published genome-wide association study (GWAS) data set in this study. METHODS: Twelve single nucleotide polymorphisms (SNPs) within novel or unpublished loci with P-value < 1.00 × 10(-02) were selected for genotype with a total of 2,442 cases and 2,798 controls(including 1,156 cases and 2,330 controls from central China, 1,012 cases and 335 controls from southern China and 274 cases and 133 controls from northern China) using Sequenom Massarry system. Associaton analyses were performed using logistic regression with sample region as a covariate through PLINK 1.07 software. RESULTS: Combined analysis in discovery and central validation dataset discovered a novel locus rs5914778 within LINC01420 associated with SLE at genome-wide significance (P = 1.00 × 10(-08); odds ratio (OR) = 1.32). We also confirmed rs5914778 in the southern Chinese sample cohort (P = 5.31 × 10(-05); OR = 1.51), and meta-analysis of the samples from the discovery, central and southern validations regions provided robust evidence for the association of rs5914778 (P = 5.26 × 10(-12); OR = 1.35). However, this SNP did not show association with SLE in the northern sample (P = 0.33). Further analysis represent the association of northern was significantly heterogeneous compared to central and southern respectively. CONCLUSIONS: Our study increases the number of established susceptibility loci for SLE in Han Chinese population and has further demonstrated the important role of X-linked genetic risk variants in the pathogenesis of SLE in Chinese Han population.

Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility.

Psoriasis is a common inflammatory skin disease with complex genetics and different degrees of prevalence across ethnic populations. Here we present the largest trans-ethnic genome-wide meta-analysis (GWMA) of psoriasis in 15,369 cases and 19,517 controls of Caucasian and Chinese ancestries. We identify four novel associations at LOC144817, COG6, RUNX1 and TP63, as well as three novel secondary associations within IFIH1 and IL12B. Fine-mapping analysis of MHC region demonstrates an important role for all three HLA class I genes and a complex and heterogeneous pattern of HLA associations between Caucasian and Chinese populations. Further, trans-ethnic comparison suggests population-specific effect or allelic heterogeneity for 11 loci. These population-specific effects contribute significantly to the ethnic diversity of psoriasis prevalence. This study not only provides novel biological insights into the involvement of immune and keratinocyte development mechanism, but also demonstrates a complex and heterogeneous genetic architecture of psoriasis susceptibility across ethnic populations.

Combined linkage and family-based association analysis improves candidate gene detection in Genetic Analysis Workshop 18 simulation data.

Because the genotype-phenotype correlation information is investigated differently by linkage and association analyses, various efforts have been made to model linkage and association jointly. However, joint modeling methods are usually computationally intensive; hence they cannot currently accommodate large pedigrees with dense markers. This article proposes a simple method to combine the linkage and association evidence with the aim of improving the detection power of disease susceptibility genes. Our detection power comparisons show that the combined linkage-association p values can improve remarkably the causal gene detection power in Genetic Analysis Workshop 18 simulation data.

DNAJ mutations are rare in Chinese Parkinson's disease patients and controls.

Mutations in DNAJC13, DNAJC6 and DNAJC5 have been implicated in Parkinson's disease (PD). To determine if rare coding variants in these genes play a role in PD risk in the Chinese population, we sequenced all coding exons of the three genes in 99 early-onset PD cases and 99 controls, and genotyped 8 missense variants in another 711 PD cases and 539 controls. Besides two common missense variants that did not show association with PD, the remaining missense variants were extremely rare (<0.5%), found in healthy population controls and did not show enrichment in PD cases. Our results suggest that missense mutations in DNAJC13, DNAJC5 and DNAJC6 do not play a major role in PD in the Chinese population.

Analysis of POFUT1 gene mutation in a Chinese family with Dowling-Degos disease.

Dowling-Degos disease (DDD) is an autosomal dominant genodermatosis characterized by reticular pigmented anomaly mainly affecting flexures. Though KRT5 has been identified to be the causal gene of DDD, the heterogeneity of this disease was displayed: for example, POFUT1 and POGLUT1 were recently identified and confirmed to be additional pathogenic genes of DDD. To identify other DDD causative genes, we performed genome-wide linkage and exome sequencing analyses in a multiplex Chinese DDD family, in which the KRT5 mutation was absent. Only a novel 1-bp deletion (c.246+5delG) in POFUT1 was found. No other novel mutation or this deletion was detected in POFUT1 in a second DDD family and a sporadic DDD case by Sanger Sequencing. The result shows the genetic-heterogeneity and complexity of DDD and will contribute to the further understanding of DDD genotype/phenotype correlations and to the pathogenesis of this disease.