Search

Objective Genetic variants spanning the ubiquitin-conjugating enzyme E2 L3 (UBE2L3) gene are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme, UbcH7, that facilitates activation of proinflammatory NF-κB signaling, and susceptibility to autoimmune diseases. This study aims to delineate how genetic variants carried on the UBE2L3-YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression. Methods We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture (3C)-qPCR, ChIP-qPCR, and siRNA knockdown assays were performed on patient-derived EBV-transformed B cells homozygous for the UBE2L3-YDJC non-risk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. Results Five of the seven prioritized variants demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. HiChIP and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3-YDJC risk haplotype, and correlated with the loss of CTCF binding and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the two promoters and reduced UBE2L3 expression. Conclusion The UBE2L3-YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

ObjectiveSjögren’s disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised.MethodsWhole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo−); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2fold change ≥2 or ≤0.5; padjLINC01871was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targetedLINC01871deletion (LINC01871−/−) and in vitro stimulation assays.ResultsWhole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene,RSAD2, in SjDRo+(r≥0.65 or ≤−0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines.LINC01871was upregulated in all SjD cases. RNA-seq and pathway analyses ofLINC01871−/−cells implicated roles in cytotoxic function, differentiation and IFNγ induction.LINC01871was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells.ConclusionLINC01871influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. AlteredLINC01871expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Friedreich ataxia (FRDA) is typically caused by homozygosity for an expanded GAA triplet-repeat in intron 1 of the FXN gene, which results in transcriptional deficiency via epigenetic silencing. Most patients are homozygous for alleles containing > 500 triplets, but a subset (~20%) have at least one expanded allele with 500 triplets, a significantly higher prevalence of unmethylated epialleles (median = 9.8%) was observed in patients with at least one allele containing 20%) and later onset (>15 years). The higher prevalence in mild FRDA of somatic FXN epialleles devoid of DNA methylation is consistent with variegated epigenetic silencing mediated by expanded triplet-repeats. The proportion of unsilenced somatic FXN genes is an unrecognized phenotypic determinant in FRDA and has implications for the deployment of effective therapies.

Woodpeckers are found in nearly every part of the world and have been important for studies of biogeography, phylogeography, and macroecology. Woodpecker hybrid zones are often studied to understand the dynamics of introgression between bird species. Notably, woodpeckers are gaining attention for their enriched levels of transposable elements (TEs) relative to most other birds. This enrichment of TEs may have substantial effects on molecular evolution. However, comparative studies of woodpecker genomes are hindered by the fact that no high-contiguity genome exists for any woodpecker species. Using hybrid assembly methods combining long-read Oxford Nanopore and short-read Illumina sequencing data, we generated a highly contiguous genome assembly for the Golden-fronted Woodpecker (Melanerpes aurifrons). The final assembly is 1.31 Gb and comprises 441 contigs plus a full mitochondrial genome. Half of the assembly is represented by 28 contigs (contig L50), each of these contigs is at least 16 Mb in size (contig N50). High recovery (92.6%) of bird-specific BUSCO genes suggests our assembly is both relatively complete and relatively accurate. Over a quarter (25.8%) of the genome consists of repetitive elements, with 287 Mb (21.9%) of those elements assignable to the CR1 superfamily of transposable elements, the highest proportion of CR1 repeats reported for any bird genome to date. Our assembly should improve comparative studies of molecular evolution and genomics in woodpeckers and allies. Additionally, the sequencing and bioinformatic resources used to generate this assembly were relatively low-cost and should provide a direction for development of high-quality genomes for studies of animal biodiversity.

OBJECTIVE Genetic variants in the region of tumor necrosis factor-induced protein 3-interacting protein 1 (TNIP1) are associated with autoimmune disease and reduced TNIP1 gene expression. The aim of this study was to define the functional genetic mechanisms driving TNIP1 hypomorphic expression imparted by the systemic lupus erythematosus-associated TNIP1 H1 risk haplotype. METHODS Dual luciferase expression and electrophoretic mobility shift assays were used to evaluate the allelic effects of 11 risk variants on enhancer function and nuclear protein binding in immune cell line models (Epstein-Barr virus [EBV]-transformed human B cells, Jurkat cells, and THP-1 cells), left in a resting state or stimulated with phorbol 12-myristate 13-acetate/ionomycin. HiChIP was used to define the regulatory 3-dimensional (3-D) chromatin network of the TNIP1 haplotype by detecting in situ long-range DNA contacts associated with H3K27ac-marked chromatin in EBV B cells. Then, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to determine the expression of genes within the 3-D chromatin network. RESULTS Bioinformatics analyses of 50 single-nucleotide polymorphisms on the TNIP1 H1 risk haplotype identified 11 non-protein-coding variants with a high likelihood of influencing TNIP1 gene expression. Eight variants in EBV B cells, 5 in THP-1 cells, and 2 in Jurkat cells exhibited various allelic effects on enhancer activation, resulting in a cumulative suppressive effect on TNIP1 expression (net effect of risk variants -7.14 fold, -6.80 fold, and -2.44 fold, respectively; n > 3). Specifically, in EBV B cells, only 2 variants (rs10057690 and rs13180950) exhibited allele-specific loss of both enhancer activity and nuclear protein binding (each P < 0.01 relative to nonrisk alleles). In contrast, the rs10036748 risk allele reduced binding affinities of the transcriptional repressors basic helix-loop-helix family member 40/differentially expressed in chondrocytes 1 (bHLHe40/DEC1) (P < 0.05 relative to nonrisk alleles) and CREB-1 (P not significant) in EBV B cells, resulting in a gain of enhancer activity (P < 0.05). HiChIP and qRT-PCR analyses revealed that overall transcriptional repression of the TNIP1 haplotype extended to the neighboring genes DCTN4 and GMA2, both of which also showed decreased expression in the presence of the TNIP1 risk haplotype (P < 0.001 and P < 0.01, respectively, relative to the nonrisk haplotype); notably, it was found that these genes share a 3-D chromatin network. CONCLUSION Hypomorphic TNIP1 expression results from the combined concordant and opposing effects of multiple risk variants carried on the TNIP1 risk haplotype, with the strongest regulatory effect in B lymphoid lineage cells. Furthermore, the TNIP1 risk haplotype effect extends to neighboring genes within a shared chromatin network.

Glioblastoma is an aggressive brain tumour found in adults, and the therapeutic approaches available have not significantly increased patient survival. Recently, we discovered that ELTD1, an angiogenic biomarker, is highly expressed in human gliomas. Polyclonal anti‐ELTD1 treatments were effective in glioma pre‐clinical models, however, pAb binding is potentially promiscuous. Therefore, the aim of this study was to determine the effects of an optimized monoclonal anti‐ELTD1 treatment in G55 xenograft glioma models. MRI was used to assess the effects of the treatments on animal survival, tumour volumes, perfusion rates and binding specificity. Immunohistochemistry and histology were conducted to confirm and characterize microvessel density and Notch1 levels, and to locate the molecular probes. RNA‐sequencing was used to analyse the effects of the mAb treatment. Our monoclonal anti‐ELTD1 treatment significantly increased animal survival, reduced tumour volumes, normalized the vasculature and showed higher binding specificity within the tumour compared with both control‐ and polyclonal‐treated mice. Notch1 positivity staining and RNA‐seq results suggested that ELTD1 has the ability to interact with and interrupt Notch1 signalling. Although little is known about ELTD1, particularly about its ligand and pathways, our data suggest that our monoclonal anti‐ELTD1 antibody is a promising anti‐angiogenic therapeutic in glioblastomas.

Treatment of glioblastoma (GBM) remains a challenge using conventional chemotherapy, such as temozolomide (TMZ), and is often ineffective as a result of drug resistance. We have assessed a novel nitrone-based agent, OKN-007, and found it to be effective in decreasing tumor volumes and increasing survival in orthotopic GBM xenografts by decreasing cell proliferation and angiogenesis and increasing apoptosis. In this study, we assessed combining OKN-007 with TMZ in vivo in a human G55 GBM orthotopic xenograft model and in vitro in TMZ-resistant and TMZ-sensitive human GBM cell lines. For the in vivo studies, magnetic resonance imaging was used to assess tumor growth and vascular alterations. Percent animal survival was also determined. For the in vitro studies, cell growth, IC50 values, RNA-seq, RT-PCR, and ELISA were used to assess growth inhibition, possible mechanism-of actions (MOAs) associated with combined OKN-007 + TMZ versus TMZ alone, and gene and protein expression levels, respectively. Microarray analysis of OKN-007–treated rat F98 glioma tumors was also carried out to determine possible MOAs of OKN-007 in glioma-bearing animals either treated or not treated with OKN-007. OKN-007 seems to elicit its effect on GBM tumors via inhibition of tumorigenic TGF-β1, which affects the extracellular matrix. When combined with TMZ, OKN-007 significantly increases percent survival, decreases tumor volumes, and normalizes tumor blood vasculature in vivo compared to untreated tumors and seems to affect TMZ-resistant GBM cells possibly via IDO-1, SUMO2, and PFN1 in vitro. Combined OKN-007 + TMZ may be a potentially potent treatment strategy for GBM patients.

Genomic variation underlying major depressive disorder (MDD) likely involves the interaction and regulation of multiple genes in a network. Data-driven co-expression network module inference has the potential to account for variation within regulatory networks, reduce the dimensionality of RNA-Seq data, and detect significant gene-expression modules associated with depression severity. We performed an RNA-Seq gene co-expression network analysis of mRNA data obtained from the peripheral blood mononuclear cells of unmedicated MDD (n = 78) and healthy control (n = 79) subjects. Across the combined MDD and HC groups, we assigned genes into modules using hierarchical clustering with a dynamic tree cut method and projected the expression data onto a lower-dimensional module space by computing the single-sample gene set enrichment score of each module. We tested the single-sample scores of each module for association with levels of depression severity measured by the Montgomery-Åsberg Depression Scale (MADRS). Independent of MDD status, we identified 23 gene modules from the co-expression network. Two modules were significantly associated with the MADRS score after multiple comparison adjustment (adjusted p = 0.009, 0.028 at 0.05 FDR threshold), and one of these modules replicated in a previous RNA-Seq study of MDD (p = 0.03). The two MADRS-associated modules contain genes previously implicated in mood disorders and show enrichment of apoptosis and B cell receptor signaling. The genes in these modules show a correlation between network centrality and univariate association with depression, suggesting that intramodular hub genes are more likely to be related to MDD compared to other genes in a module.

Age-related phenotypic changes of cerebromicrovascular endothelial cells lead to dysregulation of cerebral blood flow and blood-brain barrier disruption, promoting the pathogenesis of vascular cognitive impairment (VCI). In recent years, endothelial cell senescence has emerged as a potential mechanism contributing to microvascular pathologies opening the avenue to the therapeutic exploitation of senolytic drugs in preclinical studies. However, difficulties with the detection of senescent endothelial cells in wild type mouse models of aging hinder the assessment of the efficiency of senolytic treatments. To detect senescent endothelial cells in the aging mouse brain, we analyzed 4233 cells in fractions enriched for cerebromicrovascular endothelial cells and other cells associated with the neurovascular unit obtained from young (3-month-old) and aged (28-month-old) C57BL/6 mice. We define 13 transcriptomic cell types by deep, single-cell RNA sequencing. We match transcriptomic signatures of cellular senescence to endothelial cells identified on the basis of their gene expression profile. Our study demonstrates that with advanced aging, there is an increased ratio of senescent endothelial cells (~ 10%) in the mouse cerebral microcirculation. We propose that our single-cell RNA sequencing–based method can be adapted to study the effect of aging on senescence in various brain cell types as well as to evaluate the efficiency of various senolytic regimens in multiple tissues.

BackgroundWoodpeckers are found in nearly every part of the world, absent only from Antarctica, Australasia, and Madagascar. Woodpeckers have been important for studies of biogeography, phylogeography, and macroecology. Woodpeckers hybrid zones are often studied to understand the dynamics of introgression between bird species. Notably, woodpeckers are gaining attention for their enriched levels of transposable elements (TEs) relative to most other birds. This enrichment of TEs may have substantial effects on woodpecker molecular evolution. The Golden-fronted Woodpecker (Melanerpes aurifrons) is a member of the largest radiation of New World woodpeckers. However, comparative studies of woodpecker genomes are hindered by the fact that no high-contiguity genome exists for any woodpecker species.FindingsUsing hybrid assembly methods that combine long-read Oxford Nanopore and short-read Illumina sequencing data, we generated a highly contiguous genome assembly for the Golden-fronted Woodpecker. The final assembly is 1.31 Gb and comprises 441 contigs plus a full mitochondrial genome. Half of the assembly is represented by 28 contigs (contig N50), each of these contigs is at least 16 Mb in size (contig L50). High recovery (92.6%) of bird-specific BUSCO genes suggests our assembly is both relatively complete and relatively accurate. Accuracy is also demonstrated by the recovery of a putatively error-free mitochondrial genome. Over a quarter (25.8%) of the genome consists of repetitive elements, with 287 Mb (21.9%) of those elements assignable to the CR1 superfamily of transposable elements, the highest proportion of CR1 repeats reported for any bird genome to date.ConclusionOur assembly provides a useful tool for comparative studies of molecular evolution and genomics in woodpeckers and allies, a group emerging as important for studies on the role that TEs may play in avian evolution. Additionally, the sequencing and bioinformatic resources used to generate this assembly were relatively low-cost and should provide a direction for the development of high-quality genomes for future studies of animal biodiversity.

Previously published studies revealed that variation in expression of the DNA-binding protein ARID3a in B lymphocytes from patients with systemic lupus erythematosus (SLE) correlated with levels of disease activity ("Disease activity in systemic lupus erythematosus correlates with expression of the transcription factor AT-rich-interactive domain 3A" (J.M. Ward, K. Rose, C. Montgomery, I. Adrianto, J.A. James, J.T. Merrill et al., 2014) [1]). The data presented here compare DNA methylation patterns from SLE peripheral blood mononuclear cells obtained from samples with high numbers of ARID3a expressing B cells (ARID3a(H)) versus SLE samples with normal numbers of ARID3a(+) B cells (ARID3a(N)). The methylation data is available at the gene expression omnibus (GEO) repository, "Gene Expression Omnibus: NCBI gene expression and hybridization array data repository" (R. Edgar, M. Domrachev, A.E. Lash, 2002) [2]. Isolated B cells from SLE ARID3a(H) and ARID3a(N) B samples were also evaluated via qRT-PCR for Type I interferon (IFN) signature and pathway gene expression levels by qRT-PCR. Similarly, healthy control B cells and B cells stimulated to express ARID3a with the TLR agonist, CpG, were also compared via qRT-PCR. Primers designed to detect 6 IFNa subtype mRNAs were tested in 4 IFNa, Epstein-Barr Virus-transformed B cell lines ("Reduced interferon-alpha production by Epstein-Barr virus transformed B-lymphoblastoid cell lines and lectin-stimulated lymphocytes in congenital dyserythropoietic anemia type I" (S.H. Wickramasinghe, R. Hasan, J. Smythe, 1997) [3]). The data in this article support the publication, "Human effector B lymphocytes express ARID3a and secrete interferon alpha" (J.M. Ward, M.L. Ratliff, M.G. Dozmorov, G. Wiley, J.M. Guthridge, P.M. Gaffney, J.A. James, C.F. Webb, 2016) [4].

Background Cardiomyopathy may occur due to viral infections or drug induced heart damage. Cardiac myosin released from damaged heart has been shown to be a damage associated molecular pattern which binds to TLR2 or TLR8 and can act as an adjuvant to induce a strong autoimmune response against the heart. The result is autoimmunity against the heart which can lead to apoptosis, fibrosis and heart failure. Purpose Immune biomarkers of the early stages of heart failure are needed to identify individuals who develop progressive heart failure, do not recover their ejection fraction and may be candidates for immunotherapies. Methods Forty-one patients with myocarditis and heart failure Results Autoantibodies against human cardiac myosin and the beta-adrenergic receptor were significantly elevated in our cohort and functionally acted on cardiomyocytes to activate protein kinase A. Concomitantly, a Th17+ immunophenotype was significantly elevated in blood as well as in cardiac biopsies. CD4+IL17+ T cells (p=0.0008) and Th17-promoting cytokines TGF beta (p1year were identified using Reactome which revealed significant (FDR = 1.52E-13) enrichment of neutrophil degranulation (48 genes). Conclusion Our study illustrates a strong Th17 signature in more severe heart failure early in disease with elevated anti-cardiac myosin autoantibodies in non-recovery of left ventricular function. We observed a strong correlation with Th17-related neutrophil degranulation pathways in later disease, which may be biomarkers of fibrosis progression and disease severity in patients with heart failure. Cardiomyopathy with a Th17 signature might be treated with preventive immunomodulatory therapies such as anti-IL17A. Acknowledgement/Funding National Heart Lung and Blood Institute, Bethesda, MD, USA

Background Previous studies suggest that autoantibodies against cardiac myosin lead to dilated cardiomyopathy (DCM). Anti-cardiac myosin antibodies cross-react with the beta adrenergic receptor (βAR) and signal cAMP-dependent protein kinase A (PKA) in cardiomyocytes leading to apoptosis, fibrosis, dilated cardiomyopathy and arrhythmias. Purpose To determine if cross-reactive anti-cardiac myosin/anti-βAR autoantibodies which signal cardiomyocytes through PKA might play a role to establish DCM by promoting remodeling, apoptosis, and fibrosis. Methods Forty-one adults with DCM were enrolled Results Anti-HCM autoantibodies including autoantibody responses against 32 overlapping synthetic peptides of the S2 fragment of HCM were significantly elevated in patients whose ejection fraction did not improve over 1-year compared to those with improved ejection fraction. The human mAb confirmed our results with HCM, βAR, specific HCM peptides, and PKA signaling. RNA sequencing revealed differentially expressed genes in serum/mAb-treated cardiomyocytes compared to genes identified after RNA sequencing of peripheral blood of patients (n=10) with DCM for >1 year from onset. A primary heart cell line (H9c2-ATCC) treated with myocarditis/DCM patient sera or human mAb revealed differentially expressed genes associated with cardiac hypertrophy and heart failure, and included inflammasome component NLRP3 and complement factor H. Ingenuity Pathway Analyses revealed 27, 7, and 1 differentially expressed genes related to apoptosis, fibrosis, and hypoxia, respectively. Gene expression of CASZ1, a transcription factor important in protection against DCM, was strongly correlated with PKA signaling (r=0.89). The KDM6B gene for lysine demethylase associated with hypoxia and apoptosis pathways and was shared between cardiomyocyte and peripheral blood analysis of DCM patients. Overall, 5 genes were shared in heart failure vs in vitro Ab-treated cardiomyocyte RNA sequencing analysis: CYP4F3, KDM6B, MBOAT7, SMAP2, and DDIT4, which affects phosphorylation of mTOR to promote autophagy and cell death, cardiac hypertrophy and dysfunction. Conclusions Significantly higher responses to cardiac myosin in patients with DCM were related to lack of left ventricular function improvement and to differential expression of genes promoting apoptosis, fibrosis and disease severity. These studies identify autoantibody-directed gene signaling as a potential novel therapeutic target in DCM. Acknowledgement/Funding National Heart, Lung, and Blood Institute, Bethesda, MD, USA

Background:Sjögren’s syndrome (SS) is a chronic, heterogenous autoimmune disease characterized by inflammatory destruction of the exocrine glands. Long non-coding RNAs (lncRNAs) have emerged as a functionally diverse class of non-protein coding RNA (ncRNA) with increasing implications in interferon signaling, immune cell regulation, and autoimmune disease pathology. The potential role of lncRNAs in SS pathogenesis is unknown.Objectives:To identify and characterize candidate lncRNAs with potential relevance to SS pathology.Methods:RNA-seq was used on whole blood from SS patients (n=30 antibody negative (Ro-); n=27 antibody positive (Ro+)) and healthy controls (HC, n=27) to identify differentially expressed (DE) lncRNAs (log2 fold change (FC) ≥ 2 or ≤ 0.5; padjIn vitrotime course experiments in HSB2 T cell lymphoblasts stimulated with PMA/Ionomycin (PMA/I) or type I interferon (IFN) were used to assess biological relevance.LINC01871function was further investigated by RNA-seq on a single cell clone of HSB2 with confirmed CRISPR-targetedLINC01871deletion (LINC01871-/-).Results:We identified a total of 1054 unique DE ncRNAs between Ro+, Ro-and/or a combined analysis relative to HC; of these, 45 (1 long intergenic ncRNA (lincRNA), 1 antisense, 43 pseudogenes) were overexpressed in all 3 SS subsets. To begin investigating the function of the previously undescribed lincRNA,LINC01871(SSRo-: FC=2.85; padj=1.1x10-4), we performed a correlation analysis of the SSRo-transcriptome, which found several co-expressed protein coding RNAs involved in immune regulation (THEMIS,TBX21,IL10RA,IL2RB,among many others). Similarly, Ingenuity Pathway Analysis of the SS transcriptome compared to HC, as well as several gene ontology enrichment analyses of publicly available RNA expression correlation databases, identified shared immune-related pathways including cytotoxic T cell, natural killer cell, and T cell regulation. To further study the role ofLINC01871in cytotoxic T cells, we used qRT-PCR to resolve the effects of PMA/I or type I IFN stimulation onLINC01871expression in the T lymphoblastoid HSB2 cells.LINC01871expression was downregulated after PMA/I stimulation, but unchanged with type I IFN stimulation. To explore the regulatory function ofLINC01871in T cells, we targetedLINC01871in HSB2 cells using CRISPR. To this end, we generated a single cellLINC01871-/-clone with no RNA expression by qPCR and confirmed homozygous deletion using DNA sequencing. RNA-seq analysis ofLINC01871-/-compared to unmodified HSB2 cells identified 1166 DE transcripts. Pathway analyses clustered the DE transcripts into similar immune regulatory, cytotoxic and T cell pathways identified in SSRo-whole blood RNA-seq and publicly available RNA-seq databases. Further, several prominent T cell regulatory transcripts that exhibited correlated upregulation withLINC01871in SSRo-whole blood RNA-seq also demonstrated downregulation afterLINC01871deletion:CD109(FC=-9.7; padj=5.3x10-16),IL22(FC=-8.1; padj=7.6x10-11),PDCD1(FC=-6.2; padj=1.1x10-6),THEMIS(FC=-3.8; padj=2.7x10-165) andTBX21(FC=-2.1; padj=3.3x10-25).Conclusion:LncRNAs are emerging as important regulators of immune function with increasing evidence of autoimmune disease relevance. Here, we leveraged RNA-seq, extensive bioinformatic data, and CRISPR technology to identify and functionally characterizeLINC01871as a potential mediator of the dysregulated T cell inflammatory response pathways implicated in SS pathogenesis.Disclosure of Interests:Michelle L Joachims: None declared, Bhuwan Khatri: None declared, Kandice L Tessneer: None declared, Anna M Stolarczyk: None declared, Graham B Wiley: None declared, Astrid Rasmussen Speakers bureau: Novartis, ThermoFischer, Joel Guthridge Grant/research support from: Xencor, Bristol Myers Squibb, DXterity, Judith A. James Grant/research support from: Progentec Diagnostics, Inc, Consultant of: Abbvie, Novartis, Jannsen, R Hal Scofield Grant/research support from: Pfizer, Kathy L Sivils: None declared, Indra Adrianto: None declared, Christopher Lessard: None declared

Genomic variation underlying major depressive disorder (MDD) likely involves the interaction and regulation of multiple genes in a network. Data-driven co-expression network module inference has the potential to account for variation within regulatory networks, reduce the dimensionality of RNA-Seq data, and detect significant gene-expression modules associated with depression severity. We performed an RNA-Seq gene co-expression network analysis of mRNA data obtained from the peripheral blood mononuclear cells of unmedicated MDD (n = 78) and healthy control (n = 79) subjects. Across the combined MDD and HC groups, we assigned genes into modules using hierarchical clustering with a dynamic tree cut method and projected the expression data onto a lower-dimensional module space by computing the single-sample gene set enrichment score of each module. We tested the single-sample scores of each module for association with levels of depression severity measured by the Montgomery-Åsberg Depression Scale (MADRS). Independent of MDD status, we identified 23 gene modules from the co-expression network. Two modules were significantly associated with the MADRS score after multiple comparison adjustment (adjusted p = 0.009, 0.028 at 0.05 FDR threshold), and one of these modules replicated in a previous RNA-Seq study of MDD (p = 0.03). The two MADRS-associated modules contain genes previously implicated in mood disorders and show enrichment of apoptosis and B cell receptor signaling. The genes in these modules show a correlation between network centrality and univariate association with depression, suggesting that intramodular hub genes are more likely to be related to MDD compared to other genes in a module.

Site-1 protease (S1P), encoded by MBTPS1, is a serine protease in the Golgi. S1P regulates lipogenesis, endoplasmic reticulum (ER) function, and lysosome biogenesis in mice and in cultured cells. However, how S1P differentially regulates these diverse functions in humans has been unclear. In addition, no human disease with S1P deficiency has been identified. Here, we report a pediatric patient with an amorphic and a severely hypomorphic mutation in MBTPS1. The unique combination of these mutations results in a frequency of functional MBTPS1 transcripts of approximately 1%, a finding that is associated with skeletal dysplasia and elevated blood lysosomal enzymes. We found that the residually expressed S1P is sufficient for lipid homeostasis but not for ER and lysosomal functions, especially in chondrocytes. The defective S1P function specifically impairs activation of the ER stress transducer BBF2H7, leading to ER retention of collagen in chondrocytes. S1P deficiency also causes abnormal secretion of lysosomal enzymes due to partial impairment of mannose-6-phosphate–dependent delivery to lysosomes. Collectively, these abnormalities lead to apoptosis of chondrocytes and lysosomal enzyme–mediated degradation of the bone matrix. Correction of an MBTPS1 variant or reduction of ER stress mitigated collagen-trafficking defects. These results define a new congenital human skeletal disorder and, more importantly, reveal that S1P is particularly required for skeletal development in humans. Our findings may also lead to new therapies for other genetic skeletal diseases, as ER dysfunction is common in these disorders.

βIV spectrin links ankyrinG (AnkG) and clustered ion channels at axon initial segments (AISs) and nodes of Ranvier to the axonal cytoskeleton. Here, we report bi-allelic pathogenic SPTBN4 variants (three homozygous and two compound heterozygous) that cause a severe neurological syndrome that includes congenital hypotonia, intellectual disability, and motor axonal and auditory neuropathy. We introduced these variants into βIV spectrin, expressed these in neurons, and found that 5/7 were loss-of-function variants disrupting AIS localization or abolishing phosphoinositide binding. Nerve biopsies from an individual with a loss-of-function variant had reduced nodal Na+ channels and no nodal KCNQ2 K+ channels. Modeling the disease in mice revealed that although ankyrinR (AnkR) and βI spectrin can cluster Na+ channels and partially compensate for the loss of AnkG and βIV spectrin at nodes of Ranvier, AnkR and βI spectrin cannot cluster KCNQ2- and KCNQ3-subunit-containing K+ channels. Our findings define a class of spectrinopathies and reveal the molecular pathologies causing nervous-system dysfunction.

Autoimmune diseases develop through the exposure of genetically susceptible hosts to environmental triggers. Advances in high-throughput genotyping and sequencing technologies coupled with comprehensive databases of human genetic variation and the assembly of large cohorts of case and control subjects have led to substantial progress in defining the genetic risk factors that underlie autoimmune diseases (1). The workhorse statistical methodology has been the genome-wide association study (GWAS). Studies using the GWAS approach have convincingly and reproducibly identified ∼700 genomic regions in 183 published studies of autoimmune diseases at the level of genome-wide statistical significance (P < 10−8) (www.genome.gov/gwastudies). The majority of these genetic associations are near genes that map to critical immunoregulatory pathways, illuminating genetic effects that are shared across multiple autoimmune diseases and other genetic effects that are restricted to only a few (1). GWAS studies detect most causal variants indirectly, by leveraging linkage disequilibrium (LD) throughout the human genome. Classically defined, LD is the nonrandom association of two or more loci, resulting in segments of the genome being inherited as haplotype “blocks.” Knowledge of the allele at one variant predicts with high likelihood the alleles at the other variants on the same haplotype block (2). Though LD makes locus discovery by GWAS very efficient, because only one or two variants per haplotype block need to be genotyped in order to detect association, the high correlation of variants on associated haplotypes confounds the ability of genetic association methods to distinguish causal from noncausal variants. The UBE2L3 locus associated with multiple autoimmune diseases, including systemic lupus erythematosus (SLE) (3), illustrates the problem vividly, where 34 SLE-associated variants are located within a 67-kb haplotype, but we assume that only a few are likely to be causal (3). With a situation such as this, which variants are we to choose for functional screening? No specific guidelines exist on how to answer this question. In general, the approach to overcoming the LD problem is to first comprehensively understand the genetic architecture at a given locus. Doing this in multiple ethnic populations when possible, adds significant power to our ability to discern causal variants by allowing the comparison of haplotypes across populations. Comprehensive characterization of a locus may include the following activities: 1) locus enrichment—capture for analysis all available genetic variation present on the risk haplotype, 2) locus refinement—winnow down the associated variants within a locus to a prioritized list for functional testing, and 3) functional testing—identify allele-specific differences in biologic function that support the variant’s role(s) in causality. In this review, we discuss each of these steps and describe in more detail the available molecular methods that can provide the functional evidence required to assign causality to variants associated with autoimmune diseases.

The reniform nematode (RN), Rotylenchulus reniformis, is an agriculturally important pest with a broad host range that results in a large economic impact in tropical, subtropical and in warm temperate zones. In an initial effort to understand the transcriptome and gene expression in RN, we present EST results that reveal numerous putative parasitism-related genes some of which play roles in plant cell wall modification. The characterized contigs included 8362 (40.6%) matches to unique proteins. Coding contigs predicted were 10,656 (51.7%) or 3079 (14.9%), that was similar to those identified in Brugia malayi and Caenorhabditis elegans as reference organisms respectively. Specific transcripts studied in more detail include putative plant parasitism genes, prominent among them were several plant cell wall modification genes. Contigs matching 14 parasitism genes found in sedentary endoparasitic nematodes included expansins, hexosaminidase, glycosyl hydrolases family, 14-3-3 protein, xylanases, glutathione peroxidase, pectate lyase, β-1,4-endoglucanase, major sperm protein, aminopeptidase, c-type lectin, chitin synthase, FMR famide-like peptide, and calreticulin. These genes function in suppression of host defenses and development of feeding sites.

Previously, we determined that enhanced disease activity in patients with systemic lupus erythematosus (SLE) was associated with dramatic increases in numbers of B lymphocytes expressing the transcription factor ARID3a. Our data now indicate ARID3a is important for interferon alpha (IFNa) expression and show a strong association between ARID3a expression and transcription of genes associated with lupus IFN signatures. Furthermore, both ARID3a and IFNa production were elicited in healthy control B cells upon stimulation with the TLR 9 agonist, CpG. Importantly, secretion of IFNa from ARID3a+ healthy B lymphocytes stimulated increased IFNa production in plasmacytoid dendritic cells. These data identify ARID3a+ B cells as a novel type of effector B cell, and link ARID3a expression in B lymphocytes to IFN-associated inflammatory responses in SLE.

CD4+ T cells predominate in salivary gland (SG) inflammatory lesions in Sjogren's syndrome (SS). However, their antigen specificity, degree of clonal expansion, and relationship to clinical disease features remain unknown. We used multiplex reverse-transcriptase PCR to amplify paired T cell receptor α (TCRα) and β transcripts of single CD4+CD45RA- T cells from SG and peripheral blood (PB) of 10 individuals with primary SS, 9 of whom shared the HLA DR3/DQ2 risk haplotype. TCRα and β sequences were obtained from a median of 91 SG and 107 PB cells per subject. The degree of clonal expansion and frequency of cells expressing two productively rearranged α genes were increased in SG versus PB. Expanded clones from SG exhibited complementary-determining region 3 (CDR3) sequence similarity both within and among subjects, suggesting antigenic selection and shared antigen recognition. CDR3 similarities were shared among expanded clones from individuals discordant for canonical Ro and La autoantibodies, suggesting recognition of alternative SG antigen(s). The extent of SG clonal expansion correlated with reduced saliva production and increased SG fibrosis, linking expanded SG T cells with glandular dysfunction. Knowledge of paired TCRα and β sequences enables further work toward identification of target antigens and development of novel therapies.

The human immune system defends the body against microbes and other threats. However, if this process goes wrong the immune system can attack the body’s own healthy cells, which can lead to serious autoimmune diseases. Systemic lupus erythematosus (SLE) is an autoimmune disease in which immune cells often attack internal organs – including the kidneys, nervous system and heart. Over the past decade, multiple genes have been linked with an increased risk of SLE. However, it is largely unknown how the sequences of these genes differ between individuals with SLE and healthy individuals, and the precise changes that lead to an increased risk of SLE are also not clear. Now, Raj, Rai et al. have determined the genetic sequences of over 700 people with SLE and over 500 healthy individuals and looked for differences that influence susceptibility to the disease. The vast majority of differences were discovered in stretches of DNA that regulate the expression of nearby genes, rather than in DNA that encodes the structures of proteins. Notably, extensive differences were found in a region of the human genome that regulates the production of proteins called Human Leukocyte Antigen class II molecules; which are known to play a critical role in activating the immune system. Raj, Rai et al. found that slight changes to the regulatory DNA sequences resulted in an overabundance of these proteins, which led to a hyperactive immune system that is strongly associated with SLE. Future studies could now ask if the changes to the regulatory DNA sequences highlighted by Raj, Rai et al. increase susceptibility to other autoimmune disorders as well. It may also be possible to use the increased understanding of how the immune system is regulated to develop new ways to minimize the rejection of organ transplants.

Through exome resequencing, we identified two unique mutations in recombination signal binding protein for immunoglobulin kappa J (RBPJ) in two independent families affected by Adams-Oliver syndrome (AOS), a rare multiple-malformation disorder consisting primarily of aplasia cutis congenita of the vertex scalp and transverse terminal limb defects. These identified mutations link RBPJ, the primary transcriptional regulator for the Notch pathway, with AOS, a human genetic disorder. Functional assays confirmed impaired DNA binding of mutated RBPJ, placing it among other notch-pathway proteins altered in human genetic syndromes.

Viruses are most frequently discovered because they cause disease. To expand knowledge of plant-associated viruses beyond these narrow constraints, non-cultivated plants of the Tallgrass Prairie of the United States were systematically surveyed for evidence of viruses. This report discusses putative viruses of the family Secoviridae identified by the survey. Sequence analysis suggests the presence of at least six viruses in the study site, including Bean pod mottle virus, Maize chlorotic dwarf virus, three previously undescribed viruses within the subfamily Comovirinae and one unclassifiable virus.

Chantier qualité GA; International audience; Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation1. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species2. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ~94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa’s genomic toolbox.

Viruses are most frequently discovered because they cause disease in organisms of importance to humans. To expand knowledge of plant-associated viruses beyond these narrow constraints, non-cultivated plants of the Tallgrass Prairie Preserve, Osage County, Oklahoma, USA were systematically surveyed for evidence of the presence of viruses. This report discusses viruses of the family Tombusviridae putatively identified by the survey. Evidence of two carmoviruses, a tombusvirus, a panicovirus and an unclassifiable tombusvirid was found. The complete genome sequence was obtained for putative TGP carmovirus 1 from the legume Lespedeza procumbens, and the virus was detected in several other plant species including the fern Pellaea atropurpurea. Phylogenetic analysis of the sequence and partial sequence of a related virus supported strongly the placement of these viruses in the genus Carmovirus. Polymorphisms in the sequences suggested existence of two populations of TGP carmovirus 1 in the study area and year-to-year variations in infection by TGP carmovirus 3.

A striking evolutionary variation in eukaryotic gene expression mechanisms is the presence or absence in diverse organismal groups of a form of RNA splicing—pre-mRNA spliced leader (SL) trans-splicing (Davis 1996; Nilsen 2001; Hastings 2005). SL trans-splicing is closely related to conventional RNA splicing, or cis-splicing. The same nucleotide sequence features define donor and acceptor sites, and both processes occur in spliceosomes and involve formation of a 5′,3′,2′ branchpoint upstream of the acceptor site (Agabian 1990; Nilsen 1993). Whereas cis-splicing joins paired donor and acceptor sites within a single RNA molecule, in SL trans-splicing the donor exon is the 5′-segment of a specialized small 5′-capped RNA molecule—the SL RNA—and the target is an unpaired acceptor site near the 5′-end of a pre-mRNA molecule. Transfer of the SL RNA 5′-segment, the SL sequence, to the pre-mRNA molecule occurs with loss of the pre-mRNA's initial 5′-segment upstream of the acceptor site—a segment termed the “outron” (Conrad et al. 1991). In organisms that carry out SL trans-splicing, many different pre-mRNAs are trans-spliced with the same SL RNA species, with the result that many mature mRNA species share a common 5′-end sequence. Apart from being short, from 16 to ∼50 nucleotides (nt), the SL sequences of diverse organisms are not similar (Nilsen 2001). The best understood function of SL trans-splicing is to resolve polycistronic operon transcripts into individual 5′-capped mRNAs by trans-splicing to unpaired acceptor sites adjacent to downstream cistron open reading frames (Clayton 2002; Blumenthal and Gleason 2003). However, it is likely that SL trans-splicing has additional unknown functions because in trans-splicing metazoa the majority of trans-spliced genes are not present in operons. A variety of possible functions for SL trans-splicing in monocistronic genes have been hypothesized but not firmly established, including direct effects of the leader itself on mRNA stability or translation, and an indirect effect, i.e., the removal of potentially deleterious elements within the outron (5′ untranslated region [UTR] “sanitization”) (Hastings 2005). Due to its sporadic phylogenetic distribution, it is not clear whether SL trans-splicing is an ancestral eukaryotic feature that has been lost in several lineages or whether it was absent from the ancestor and arose independently in several lineages (Nilsen 2001; Stover and Steele 2001). In the deuterostome division of the metazoa, we discovered SL trans-splicing in the chordate tunicate ascidian Ciona intestinalis (Vandenberghe et al. 2001) and it has since been found in several other tunicate species (Yuasa et al. 2002; Ganot et al. 2004). Tunicates are of particular evolutionary interest because, as chordates, they are related to vertebrate ancestors (Dehal et al. 2002). Moreover, because vertebrates are among the groups that almost certainly do not carry out SL trans-splicing (Nilsen 2001; Hastings 2005), it follows that either SL trans-splicing has been lost in the vertebrate lineage or was invented in the tunicate lineage after the tunicate/vertebrate divergence. Thus, in-depth knowledge of trans-splicing in the tunicates may provide valuable insight into genome evolution in the chordates and the evolutionary dynamics of RNA splicing. Among the metazoan organisms that carry out SL trans-splicing, a significant fraction of genes are conventionally expressed, i.e., are not trans-spliced (e.g., in Ciona, ∼50% of genes are apparently not trans-spliced; Satou et al. 2006). It is not known why some monocistronic genes are trans-spliced and others are not. Genome-wide knowledge of the trans-splicing status of individual genes would provide a basis for elucidating those aspects of gene structure or function that may affect the “choice” of SL trans-splicing versus conventional gene expression. Although there have been several studies of individual genes or of small samples of the trans-spliced or non-trans-spliced gene sets (Satou et al. 2006; Sierro et al. 2009), there has not yet been a genome-wide identification of trans-spliced and non-trans-spliced gene populations in any organism. Because the trans-spliced and non-trans-spliced gene subpopulations represent a significant fraction of the genome, generating a comprehensive overview has been beyond the reach of conventional sequencing methodologies. However, recently developed high-throughput methods have made it feasible to study genome-wide processes through DNA sequencing. We have developed and employed an approach to high-throughput characterization of the trans-spliced mRNA population of the ascidan Ciona intestinalis based on 454 Life Sciences (Roche) pyrosequencing. In order to sample a wide range of expressed genes we used the whole organism at a stage of active development and differentiation, the tailbud embryo. Our results identify the majority of trans-spliced genes in Ciona and precisely localize their trans-splice acceptor sites. Analysis of this extensive data set provides new insight into the splicing mechanism and into the nature of the trans-spliced and non-trans-spliced gene sets. In addition to improved understanding of trans-splicing, our results also provide a wealth of specific genetic information on Ciona, a key model organism for genomic and developmental genetics studies relevant to vertebrate early development and evolution (Satou et al. 2005; Imai et al. 2006; Munro et al. 2006).

We used an expressed sequence tag and 454 pyrosequencing approach to initiate a study of the genome of the screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae). Two normalized cDNA libraries were constructed from RNA isolated from embryos and second instar larvae from the Panama 95 strain. Approximately 5,400 clones from each library were sequenced from both the 5' and 3' directions using the Sanger method. In addition, double-stranded cDNA was prepared from random-primed polyA RNA purified from embryos, second-instar larvae, adult males, and adult females. These four cDNA samples were used for 454 pyrosequencing that produced approximately 300,000 independent sequences. Sequences were assembled into a database of assembled contigs and singletons and used to search public protein databases and annotate the sequences. The full database consists of 6,076 contigs and 58,221 singletons assembled from both the traditional expressed sequence tag (EST) and 454 sequences. Annotation of the data led to the identification of several gene coding regions with possible roles in sex determination in the screwworm. This database will facilitate the design of microarray and other experiments to study screwworm gene expression on a larger scale than previously possible.

The diversity of viruses associated with non-cultivated plants was assessed from plant samples collected in the Tallgrass Prairie Preserve of northeastern Oklahoma, USA. The samples were processed to determine the sequences of nucleic acids extracted from the virus-like particle fraction of plant homogenates. Sequences from 95 specimens of 52 plant species included those of probable origin from the genomes of plants (including retroelements), bacteria, fungi, other organisms, and viruses. Virus-like sequences were identified in sequences from 25% of the specimens, coming from 19% of the plant species. Evidence of a member of the genus Tymovirus was found in 16 specimens of 6 plant species, making it the most predominant virus associated with the sampled plants. There was evidence of the presence of more than one virus in each of six specimens.

The horn fly, Haematobia irritans L., is an obligate blood-feeding parasite of cattle, and control of this pest is a continuing problem because the fly is becoming resistant to pesticides. Dominant conditional lethal gene systems are being studied as population control technologies against agricultural pests. One of the components of these systems is a female-specific gene promoter that drives expression of a lethality-inducing gene. To identify candidate genes to supply this promoter, microarrays were designed from a horn fly expressed sequence tag (EST) database and probed to identify female-specific and larval-specific gene expression. Analysis of dye swap experiments found 432 and 417 transcripts whose expression levels were higher or lower in adult female flies, respectively, compared with adult male flies. Additionally, 419 and 871 transcripts were identified whose expression levels were higher or lower in first-instar larvae compared with adult flies, respectively. Three transcripts were expressed more highly in adult females flies compared with adult males and also higher in the first-instar larval lifestage compared with adult flies. One of these transcripts, a putative nanos ortholog, has a high female-to-male expression ratio, a moderate expression level in first-instar larvae, and has been well characterized in Drosophila. melanogaster (Meigen). In conclusion, we used microarray technology, verified by reverse transcriptase-polymerase chain reaction and massively parallel pyrosequencing, to study life stage- and sex-specific gene expression in the horn fly and identified three gene candidates for detailed evaluation as a gene promoter source for the development of a female-specific conditional lethality system.

To test the hypothesis that many viruses remain to be discovered in plants, a procedure was developed to sequence nucleic acids cloned randomly from virus-like particle fractions of plant homogenates. As a test of the efficiency of the procedure we targeted Ambrosia psilostachya, western ragweed, plants growing at the Tallgrass Prairie Preserve of northeastern Oklahoma. Amplifiable nucleic acid was found in the fractions from six of twelve specimens and sequences were characterized from four of them. Evidence was obtained for the presence of viruses belonging to two families (Caulimoviridae, Flexiviridae). Multiple viral species were found in two of the four specimens and their level within the isolated nucleic acid population varied from less than 1-37%. None of the sequences were derived from reported sequences of known viruses. Thus, the analysis of nucleic acid from virus-like particles is a useful tool to expand our knowledge of the universe of viruses to non-cultivated species.

The ascomycete Epichloë festucae is a model endophyte that 1) switches between mutualistic and antagonistic states, 2) is seed transmissible, 3) has a sexual state amenable to genetic analysis, and 4) is rich in bioprotective alkaloids. This fungus grows systemically and intercellularly throughout the life of its host plant. On each reproductive tiller the fungus either infects benignly and transmits clonally in seeds, or produces its sexual state (stroma) and chokes inflorescence development. The E. festucae genome was estimated at 29 Mb in six chromosomes. The genome sequence was assembled from cloned insert end reads (4.2 x coverage) and preassembled pyrosequencing reads (454-sequencing: 20 x raw, 1.7 x assembled), giving 3967 supercontigs, of which 1004 were larger than 2 kb and covered 92% of the genome. Gene prediction with FGENESH identified ~10,000 putative genes. We also sequenced 25,000 ESTs from each of two normalised libraries — one of choked inflorescences, the other of benignly infected inflorescences — yielding 5077 E. festucae unigenes, annotated by BLAST and InterPro. Sequence data and annotations are stored in a database for visualisation and inspection with the GBrowse browser. The genomic sequences can be queried by BLAST at http://www.genome.ou.edu/blast/ ef_blastall.html. Keywords: bioinformatics, DNA sequence, Epichloë festucae, expressed sequence tags, Festucae pratensis, fungal genomics, Lolium pratense

Antiviral defenses are inappropriately activated in systemic lupus erythematosus (SLE) and association between SLE and the antiviral helicase gene, IFIH1, is well established. We sought to extend the previously reported association of pathogenic soluble mediators and autoantibodies with mouse Mda5 to its human ortholog, IFIH1. To better understand the role this gene plays in human lupus, we assessed association of IFIH1 variants with soluble mediators and autoantibodies in 357 European-American SLE patients, first-degree relatives, and unrelated, unaffected healthy controls. Association between each of 135 genotyped SNPs in IFIH1 and four lupus-associated plasma mediators, IL-6, TNF-α, IFN-β, and IP-10, were investigated via linear regression. No significant associations were found to SNPs orthologous to those identified in exon 13 of the mouse. However, outside of this region there were significant associations between IL-6 and rs76162067 (p = 0.008), as well as IP-10 and rs79711023 (p = 0.003), located in a region of IFIH1 previously shown to directly influence MDA-5 mediated IP-10 and IL-6 secretion. SLE patients and FDRs carrying the minor allele for rs79711023 demonstrated lower levels of IP-10, while only FDRs carrying the minor allele for rs76162067 demonstrated an increased level of IL-6. This would suggest that the change in IP-10 is genotypically driven, while the change in IL-6 may be reflective of SLE transition status. These data suggest that IFIH1 may contribute to SLE pathogenesis via altered inflammatory mechanisms.

The simultaneous targeting of host and pathogen processes represents an untapped approach for the treatment of intracellular infections. Hypoxia-inducible factor-1 (HIF-1) is a host cell transcription factor that is activated by and required for the growth of the intracellular protozoan parasite Toxoplasma gondii at physiological oxygen levels. Parasite activation of HIF-1 is blocked by inhibiting the family of closely related Activin-Like Kinase (ALK) host cell receptors ALK4, ALK5, and ALK7, which was determined in part by use of an ALK4,5,7 inhibitor named SB505124. Besides inhibiting HIF-1 activation, SB505124 also potently blocks parasite replication under normoxic conditions. To determine whether SB505124 inhibition of parasite growth was exclusively due to inhibition of ALK4,5,7 or because the drug inhibited a second kinase, SB505124-resistant parasites were isolated by chemical mutagenesis. Whole-genome sequencing of these mutants revealed mutations in the Toxoplasma MAP kinase, TgMAPK1. Allelic replacement of mutant TgMAPK1 alleles into wild-type parasites was sufficient to confer SB505124 resistance. SB505124 independently impacts TgMAPK1 and ALK4,5,7 signaling since drug resistant parasites could not activate HIF-1 in the presence of SB505124 or grow in HIF-1 deficient cells. In addition, TgMAPK1 kinase activity is inhibited by SB505124. Finally, mice treated with SB505124 had significantly lower tissue burdens following Toxoplasma infection. These data therefore identify SB505124 as a novel small molecule inhibitor that acts by inhibiting two distinct targets, host HIF-1 and TgMAPK1., Author Summary Understanding how a compound blocks growth of an intracellular pathogen is important not only for developing these compounds into drugs that can be prescribed to patients, but also because these data will likely provide novel insight into the biology of these pathogens. Forward genetic screens are one established approach towards defining these mechanisms. But performing these screens with intracellular parasites has been limited not only because of technical limitations but also because the compounds may have off-target effects in either the host or parasite. Here, we report the first compound that kills a pathogen by simultaneously inhibiting distinct host- and parasite-encoded targets. Because developing drug resistance simultaneously to two targets is less likely, this work may highlight a new approach to antimicrobial drug discovery.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Signaling through the store-operated Ca(2+) release-activated Ca(2+) (CRAC) channel regulates critical cellular functions, including gene expression, cell growth and differentiation, and Ca(2+) homeostasis. Loss-of-function mutations in the CRAC channel pore-forming protein ORAI1 or the Ca(2+) sensing protein stromal interaction molecule 1 (STIM1) result in severe immune dysfunction and nonprogressive myopathy. Here, we identify gain-of-function mutations in the cytoplasmic domain of STIM1 (p.R304W) associated with thrombocytopenia, bleeding diathesis, miosis, and tubular myopathy in patients with Stormorken syndrome, and in ORAI1 (p.P245L), associated with a Stormorken-like syndrome of congenital miosis and tubular aggregate myopathy but without hematological abnormalities. Heterologous expression of STIM1 p.R304W results in constitutive activation of the CRAC channel in vitro, and spontaneous bleeding accompanied by reduced numbers of thrombocytes in zebrafish embryos, recapitulating key aspects of Stormorken syndrome. p.P245L in ORAI1 does not make a constitutively active CRAC channel, but suppresses the slow Ca(2+)-dependent inactivation of the CRAC channel, thus also functioning as a gain-of-function mutation. These data expand our understanding of the phenotypic spectrum of dysregulated CRAC channel signaling, advance our knowledge of the molecular function of the CRAC channel, and suggest new therapies aiming at attenuating store-operated Ca(2+) entry in the treatment of patients with Stormorken syndrome and related pathologic conditions.

Genotyping variants in the human genome has proven to be an efficient method to identify genetic associations with phenotypes. The distribution of variants within families or populations can facilitate identification of the genetic factors of disease. Illumina's panel of genotyping BeadChips allows investigators to genotype thousands or millions of single nucleotide polymorphisms (SNPs) or to analyze other genomic variants, such as copy number, across a large number of DNA samples. These SNPs can be spread throughout the genome or targeted in specific regions in order to maximize potential discovery. The Infinium assay has been optimized to yield high-quality, accurate results quickly. With proper setup, a single technician can process from a few hundred to over a thousand DNA samples per week, depending on the type of array. This assay guides users through every step, starting with genomic DNA and ending with the scanning of the array. Using propriety reagents, samples are amplified, fragmented, precipitated, resuspended, hybridized to the chip, extended by a single base, stained, and scanned on either an iScan or Hi Scan high-resolution optical imaging system. One overnight step is required to amplify the DNA. The DNA is denatured and isothermally amplified by whole-genome amplification; therefore, no PCR is required. Samples are hybridized to the arrays during a second overnight step. By the third day, the samples are ready to be scanned and analyzed. Amplified DNA may be stockpiled in large quantities, allowing bead arrays to be processed every day of the week, thereby maximizing throughput.

Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10(-9)), and rs13306575 in HS and KR (P(HS) = 7.04 × 10(-7), P(KR) = 3.30 × 10(-3)). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10(-7)), implying that SLE predisposing variants were tagged. Significant SNP-SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance ('missing heritability') of complex diseases like SLE.

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.

Functional characterization of causal variants present on risk haplotypes identified through genome-wide association studies (GWAS) is a primary objective of human genetics. In this report, we evaluate the function of a pair of tandem polymorphic dinucleotides, 42 kb downstream of the promoter of TNFAIP3, (rs148314165, rs200820567, collectively referred to as TT>A) recently nominated as causal variants responsible for genetic association of systemic lupus erythematosus (SLE) with tumor necrosis factor alpha inducible protein 3 (TNFAIP3). TNFAIP3 encodes the ubiquitin-editing enzyme, A20, a key negative regulator of NF-κB signaling. A20 expression is reduced in subjects carrying the TT>A risk alleles; however, the underlying functional mechanism by which this occurs is unclear. We used a combination of electrophoretic mobility shift assays (EMSA), mass spectrometry (MS), reporter assays, chromatin immunoprecipitation-PCR (ChIP-PCR) and chromosome conformation capture (3C) EBV transformed lymphoblastoid cell lines (LCL) from individuals carrying risk and non-risk TNFAIP3 haplotypes to characterize the effect of TT>A on A20 expression. Our results demonstrate that the TT>A variants reside in an enhancer element that binds NF-κB and SATB1 enabling physical interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB to the TT>A risk alleles or knockdown of SATB1 expression by shRNA, inhibits the looping interaction resulting in reduced A20 expression. Together, these data reveal a novel mechanism of TNFAIP3 transcriptional regulation and establish the functional basis by which the TT>A risk variants attenuate A20 expression through inefficient delivery of NF-κB to the TNFAIP3 promoter. These results provide critical functional evidence supporting a direct causal role for TT>A in the genetic predisposition to SLE., Author Summary A key objective of human genetics is the identification and characterization of variants responsible for association with complex diseases. A pair of single nucleotide polymorphisms (rs148314165, rs200820567) 42 kb downstream from the promoter of TNFAIP3, have been proposed as the variants responsible for association with systemic lupus erythematosus based on comprehensive genetic and bioinformatic analyses. TNFAIP3 encodes for the ubiquitin-editing enzyme, A20, which plays a central role in maintaining immune system homeostasis through restriction of NF-κB signaling. Cells that carry this risk haplotype express low levels of TNFAIP3 compared to cells carrying the nonrisk haplotype. How the risk alleles of rs148314165 and rs200820567 might influence low TNFAIP3 expression is unknown. In this paper, we demonstrate that these variants reside in an enhancer element that binds NF-κB and SATB1 enabling the interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB directly to the risk alleles or shRNA-mediated knockdown of SATB1 inhibits interaction of the enhancer with the TNFAIP3 promoter resulting in reduced A20 expression. These results clarify the functional mechanism by which rs148314165 and rs200820567 attenuate A20 expression and support a causal role for these variants in the predisposition to autoimmune disease.

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Objectives: Homozygous C1q deficiency is an extremely rare condition and strongly associated with systemic lupus erythematosus. To assess and characterize C1q deficiency in an African-American lupus pedigree, C1q genomic region was evaluated in the lupus cases and family members. Methods: Genomic DNA from patient was obtained and C1q A, B and C gene cluster was sequenced using next generation sequencing method. The identified mutation was further confirmed by direct Sanger sequencing method in the patient and all blood relatives. C1q levels in serum were measured using sandwich ELISA method. Results: In an African-American patient with lupus and C1q deficiency, we identified and confirmed a novel homozygote start codon mutation in C1qA gene that changes amino acid methionine to arginine at position 1. The Met1Arg mutation prevents protein translation (Met1Arg). Mutation analyses of the patient’s family members also revealed the Met1Arg homozygote mutation in her deceased brother who also had lupus with absence of total complement activity consistent with a recessive pattern of inheritance. Conclusion: The identification of new mutation in C1qA gene that disrupts the start codon (ATG to AGG (Met1Arg)) has not been reported previously and it expands the knowledge and importance of the C1q gene in the pathogenesis of lupus especially in the high-risk African-American population.

Pyrosequencing analysis of 16S rRNA genes was used to examine impacts of elevated CO 2 (eCO 2 ) on soil microbial communities from 12 replicates each from ambient CO 2 (aCO 2 ) and eCO 2 settings. The results suggest that the soil microbial community composition and structure significantly altered under conditions of eCO 2 , which was closely associated with soil and plant properties.

Systemic lupus erythematosus (SLE) is an autoimmune disease with diverse clinical manifestations characterized by the development of pathogenic autoantibodies manifesting in inflammation of target organs such as the kidneys, skin and joints. Genome-wide association studies have identified genetic variants in the UBE2L3 region that are associated with SLE in subjects of European and Asian ancestry. UBE2L3 encodes an ubiquitin-conjugating enzyme, UBCH7, involved in cell proliferation and immune function. In this study, we sought to further characterize the genetic association in the region of UBE2L3 and use molecular methods to determine the functional effect of the risk haplotype. We identified significant associations between variants in the region of UBE2L3 and SLE in individuals of European and Asian ancestry that exceeded a Bonferroni-corrected threshold (P

Objective Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-?B signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-?B pathway. Methods We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. Results We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. Conclusion Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis. Copyright © 2012 by the American College of Rheumatology.