Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription.

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.

Allelic heterogeneity in NCF2 associated with systemic lupus erythematosus (SLE) susceptibility across four ethnic populations.

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.

Forward genetic screening identifies a small molecule that blocks Toxoplasma gondii growth by inhibiting both host- and parasite-encoded kinases.

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.

An enhancer element harboring variants associated with systemic lupus erythematosus engages the TNFAIP3 promoter to influence A20 expression.

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.

RBPJ mutations identified in two families affected by Adams-Oliver syndrome.

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.

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study.

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p(meta-Euro) = 2.08 × 10(-10)), transmembrane protein 39A (TMEM39A; rs1132200; p(meta-all) = 8.62 × 10(-9)), and 17q21 (rs1453560; p(meta-all) = 3.48 × 10(-10)) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10(-8) < p(meta-Euro) < 9.99 × 10(-5)) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.

High-throughput sequence analysis of Ciona intestinalis SL trans-spliced mRNAs: alternative expression modes and gene function correlates.

Pre-mRNA 5' spliced-leader (SL) trans-splicing occurs in some metazoan groups but not in others. Genome-wide characterization of the trans-spliced mRNA subpopulation has not yet been reported for any metazoan. We carried out a high-throughput analysis of the SL trans-spliced mRNA population of the ascidian tunicate Ciona intestinalis by 454 Life Sciences (Roche) pyrosequencing of SL-PCR-amplified random-primed reverse transcripts of tailbud embryo RNA. We obtained approximately 250,000 high-quality reads corresponding to 8790 genes, approximately 58% of the Ciona total gene number. The great depth of this data revealed new aspects of trans-splicing, including the existence of a significant class of "infrequently trans-spliced" genes, accounting for approximately 28% of represented genes, that generate largely non-trans-spliced mRNAs, but also produce trans-spliced mRNAs, in part through alternative promoter use. Thus, the conventional qualitative dichotomy of trans-spliced versus non-trans-spliced genes should be supplanted by a more accurate quantitative view recognizing frequently and infrequently trans-spliced gene categories. Our data include reads representing approximately 80% of Ciona frequently trans-spliced genes. Our analysis also revealed significant use of closely spaced alternative trans-splice acceptor sites which further underscores the mechanistic similarity of cis- and trans-splicing and indicates that the prevalence of +/-3-nt alternative splicing events at tandem acceptor sites, NAGNAG, is driven by spliceosomal mechanisms, and not nonsense-mediated decay, or selection at the protein level. The breadth of gene representation data enabled us to find new correlations between trans-splicing status and gene function, namely the overrepresentation in the frequently trans-spliced gene class of genes associated with plasma/endomembrane system, Ca(2+) homeostasis, and actin cytoskeleton.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

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.

Variants on the UBE2L3/YDJC Autoimmune Disease Risk Haplotype Increase UBE2L3 Expression by Modulating CCCTC-Binding Factor and YY1 Binding.

OBJECTIVE: Genetic variants spanning UBE2L3 are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme H7 (UbcH7), which facilitates activation of proinflammatory NF-κB signaling and susceptibility to autoimmune diseases. We undertook this study 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 with quantitative polymerase chain reaction (3C-qPCR), chromatin immunoprecipitation (ChIP)-qPCR, and small interfering RNA (siRNA) knockdown assays were performed on patient-derived Epstein-Barr virus-transformed B cells homozygous for the UBE2L3/YDJC nonrisk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. RESULTS: Of the 7 prioritized variants, 5 demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. High-throughput sequencing of chromosome conformation capture coupled with ChIP (HiChIP) and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and the downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3/YDJC risk haplotype, and correlated with the loss of CCCTC-binding factor (CTCF) and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the 2 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.

Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses.

OBJECTIVE: Sjö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. METHODS: Whole 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 (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays. RESULTS: Whole 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. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was 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. CONCLUSION: LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Ecogenomics: using massively parallel pyrosequencing to understand virus ecology.

Environmental samples have been analysed for viruses in metagenomic studies, but these studies have not linked individual viruses to their hosts. We designed a strategy to isolate double-stranded RNA, a hallmark of RNA virus infection, from individual plants and convert this to cDNA with a unique four nucleotide Tag at each end. Using 96 different Tags allowed us to pool samples and still retain the link to the original sample. We then analysed the sequence of pooled samples using massively parallel sequencing with Roche 454 pyrosequencing such that 384 samples could be assessed per picotiter plate. Using this method we have been able to analyse thousands of plants, and we have discovered several thousand new plant viruses, all linked to their specific plant hosts. Here we describe the method in detail, including the results and analysis for eight pools of samples. This technology will be extremely useful in understanding the full scope of plant virus biodiversity.

Role of Systemic Lupus Erythematosus Risk Variants With Opposing Functional Effects as a Driver of Hypomorphic Expression of TNIP1 and Other Genes Within a Three-Dimensional Chromatin Network.

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.

Non-cultivated plants of the Tallgrass Prairie Preserve of northeastern Oklahoma frequently contain virus-like sequences in particulate fractions.

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.

Cloning and characterization of small RNAs from Medicago truncatula reveals four novel legume-specific microRNA families.

MicroRNAs (miRNAs) and small-interfering RNAs (siRNAs) have emerged as important regulators of gene expression in higher eukaryotes. Recent studies indicate that genomes in higher plants encode lineage-specific and species-specific miRNAs in addition to the well-conserved miRNAs. Leguminous plants are grown throughout the world for food and forage production. To date the lack of genomic sequence data has prevented systematic examination of small RNAs in leguminous plants. Medicago truncatula, a diploid plant with a near-completely sequenced genome has recently emerged as an important model legume. We sequenced a small RNA library generated from M. truncatula to identify not only conserved miRNAs but also novel small RNAs, if any. Eight novel small RNAs were identified, of which four (miR1507, miR2118, miR2119 and miR2199) are annotated as legume-specific miRNAs because these are conserved in related legumes. Three novel transcripts encoding TIR-NBS-LRR proteins are validated as targets for one of the novel miRNA, miR2118. Small RNA sequence analysis coupled with the small RNA blot analysis, confirmed the expression of around 20 conserved miRNA families in M. truncatula. Fifteen transcripts have been validated as targets for conserved miRNAs. We also characterized Tas3-siRNA biogenesis in M. truncatula and validated three auxin response factor (ARF) transcripts that are targeted by tasiRNAs. These findings indicate that M. truncatula and possibly other related legumes have complex mechanisms of gene regulation involving specific and common small RNAs operating post-transcriptionally.

Microarray analysis of female- and larval-specific gene expression in the horn fly (Diptera: Muscidae).

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.

The genome of the versatile nitrogen fixer Azorhizobium caulinodans ORS571.

BACKGROUND: Biological nitrogen fixation is a prokaryotic process that plays an essential role in the global nitrogen cycle. Azorhizobium caulinodans ORS571 has the dual capacity to fix nitrogen both as free-living organism and in a symbiotic interaction with Sesbania rostrata. The host is a fast-growing, submergence-tolerant tropical legume on which A. caulinodans can efficiently induce nodule formation on the root system and on adventitious rootlets located on the stem. RESULTS: The 5.37-Mb genome consists of a single circular chromosome with an overall average GC of 67% and numerous islands with varying GC contents. Most nodulation functions as well as a putative type-IV secretion system are found in a distinct symbiosis region. The genome contains a plethora of regulatory and transporter genes and many functions possibly involved in contacting a host. It potentially encodes 4717 proteins of which 96.3% have homologs and 3.7% are unique for A. caulinodans. Phylogenetic analyses show that the diazotroph Xanthobacter autotrophicus is the closest relative among the sequenced genomes, but the synteny between both genomes is very poor. CONCLUSION: The genome analysis reveals that A. caulinodans is a diazotroph that acquired the capacity to nodulate most probably through horizontal gene transfer of a complex symbiosis island. The genome contains numerous genes that reflect a strong adaptive and metabolic potential. These combined features and the availability of the annotated genome make A. caulinodans an attractive organism to explore symbiotic biological nitrogen fixation beyond leguminous plants.

Evidence for novel viruses by analysis of nucleic acids in virus-like particle fractions from Ambrosia psilostachya.

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.

Enhancer histone-QTLs are enriched on autoimmune risk haplotypes and influence gene expression within chromatin networks.

Genetic variants can confer risk to complex genetic diseases by modulating gene expression through changes to the epigenome. To assess the degree to which genetic variants influence epigenome activity, we integrate epigenetic and genotypic data from lupus patient lymphoblastoid cell lines to identify variants that induce allelic imbalance in the magnitude of histone post-translational modifications, referred to herein as histone quantitative trait loci (hQTLs). We demonstrate that enhancer hQTLs are enriched on autoimmune disease risk haplotypes and disproportionately influence gene expression variability compared with non-hQTL variants in strong linkage disequilibrium. We show that the epigenome regulates HLA class II genes differently in individuals who carry HLA-DR3 or HLA-DR15 haplotypes, resulting in differential 3D chromatin conformation and gene expression. Finally, we identify significant expression QTL (eQTL) x hQTL interactions that reveal substructure within eQTL gene expression, suggesting potential implications for functional genomic studies that leverage eQTL data for subject selection and stratification.

Single-cell analysis of glandular T cell receptors in Sjögren's syndrome.

CD4+ T cells predominate in salivary gland (SG) inflammatory lesions in Sjögren'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.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity.

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.