Senataxin Mutation Reveals How R-Loops Promote Transcription by Blocking DNA Methylation at Gene Promoters.

R-loops are three-stranded nucleic acid structures found abundantly and yet often viewed as by-products of transcription. Studying cells from patients with a motor neuron disease (amyotrophic lateral sclerosis 4 [ALS4]) caused by a mutation in senataxin, we uncovered how R-loops promote transcription. In ALS4 patients, the senataxin mutation depletes R-loops with a consequent effect on gene expression. With fewer R-loops in ALS4 cells, the expression of BAMBI, a negative regulator of transforming growth factor ß (TGF-ß), is reduced; that then leads to the activation of the TGF-ß pathway. We uncovered that genome-wide R-loops influence promoter methylation of over 1,200 human genes. DNA methyl-transferase 1 favors binding to double-stranded DNA over R-loops. Thus, in forming R-loops, nascent RNA blocks DNA methylation and promotes further transcription. Hence, our results show that nucleic acid structures, in addition to sequences, influence the binding and activity of regulatory proteins.

cis Elements that Mediate RNA Polymerase II Pausing Regulate Human Gene Expression.

Aberrant gene expression underlies many human diseases. RNA polymerase II (Pol II) pausing is a key regulatory step in transcription. Here, we mapped the locations of RNA Pol II in normal human cells and found that RNA Pol II pauses in a consistent manner across individuals and cell types. At more than 1,000 genes including MYO1E and SESN2, RNA Pol II pauses at precise nucleotide locations. Characterization of these sites shows that RNA Pol II pauses at GC-rich regions that are marked by a sequence motif. Sixty-five percent of the pause sites are cytosines. By differential allelic gene expression analysis, we showed in our samples and a population dataset from the Genotype-Tissue Expression (GTEx) consortium that genes with more paused polymerase have lower expression levels. Furthermore, mutagenesis of the pause sites led to a significant increase in promoter activities. Thus, our data uncover that RNA Pol II pauses precisely at sites with distinct sequence features that in turn regulate gene expression.

Human proteins that interact with RNA/DNA hybrids.

RNA/DNA hybrids form when RNA hybridizes with its template DNA generating a three-stranded structure known as the R-loop. Knowledge of how they form and resolve, as well as their functional roles, is limited. Here, by pull-down assays followed by mass spectrometry, we identified 803 proteins that bind to RNA/DNA hybrids. Because these proteins were identified using in vitro assays, we confirmed that they bind to R-loops in vivo. They include proteins that are involved in a variety of functions, including most steps of RNA processing. The proteins are enriched for K homology (KH) and helicase domains. Among them, more than 300 proteins preferred binding to hybrids than double-stranded DNA. These proteins serve as starting points for mechanistic studies to elucidate what RNA/DNA hybrids regulate and how they are regulated.

The association of the UHRF1BP1 gene with systemic lupus erythematosus was replicated in a Han Chinese population from mainland China.

Single-nucleotide polymorphisms (SNPs) in the UHRF gene have been shown to be associated with systemic lupus erythematosus (SLE) in European and Hong Kong Chinese, but statistically significant evidence for association has not been found in a mainland Han Chinese population. Therefore, we selected SNP rs13205210 located in UHRF1BP1 as a candidate association from our previously published genome-wide association study (GWAS) data of SLE (1,047 cases and 1,205 controls from a mainland Han Chinese population) to explore the association between the UHRF1BP1 gene and SLE. We conducted a large-scale replication study in an additional independent sample of 3,509 cases and 8,246 controls from a mainland Han Chinese population. Real-time PCR was used to determine gene expression differences in peripheral blood mononuclear cells (PBMCs) from cases and controls. As a result, we replicated the association between the UHRF1BP1 gene and SLE (rs13205210, missense, Pmeta  = 2.26E-17, odds ratio = 1.41) by a meta-analysis of our previous GWAS and this replication study involving a total of 4,556 cases and 9,451 controls. The UHRF1BP1 mRNA expression level in PBMCs was significantly decreased in patients with SLE compared with that in healthy controls. SNP rs13205210 exhibited an expression quantitative trait loci effect on the UHRF1BP1 gene in PBMCs from patients. In conclusion, this study not only suggests that the UHRF1BP1 gene was associated with SLE in a mainland Han Chinese population, but also implied that it might be a common genetic factor contributing to SLE susceptibility in multiple populations.

RNA-DNA sequence differences in Saccharomyces cerevisiae.

Alterations of RNA sequences and structures, such as those from editing and alternative splicing, result in two or more RNA transcripts from a DNA template. It was thought that in yeast, RNA editing only occurs in tRNAs. Here, we found that Saccharomyces cerevisiae have all 12 types of RNA-DNA sequence differences (RDDs) in the mRNA. We showed these sequence differences are propagated to proteins, as we identified peptides encoded by the RNA sequences in addition to those by the DNA sequences at RDD sites. RDDs are significantly enriched at regions with R-loops. A screen of yeast mutants showed that RDD formation is affected by mutations in genes regulating R-loops. Loss-of-function mutations in ribonuclease H, senataxin, and topoisomerase I that resolve RNA-DNA hybrids lead to increases in RDD frequency. Our results demonstrate that RDD is a conserved process that diversifies transcriptomes and proteomes and provide a mechanistic link between R-loops and RDDs.

Bach2 regulates aberrant activation of B cell in systemic lupus erythematosus and can be negatively regulated by BCR-ABL/PI3K.

OBJECTIVE: This study was aimed to explore the effect of Bach2 on B cells in systemic lupus erythematosus (SLE), as well as the underlying mechanisms. METHODS: Expression of Bach2, phosphorylated-Bach2 (p-Bach2), Akt, p-Akt and BCR-ABL (p210) in B cells isolated from SLE patients and the healthy persons were assessed by Western blot. Immunofluorescence staining was performed to assess the localization of Bach2 in B cells. Enzyme-linked immunosorbent assay (ELISA) was employed to detect IgG produced by B cells. Cell counting kit-8 (CCK-8) and Annexin-V FITC/PI double staining assay were adopted to evaluate cell proliferation and apoptosis in B cells, respectively. RESULTS: Compared to the healthy controls, Bach2, p-Akt and p210 were significantly decreased, while nuclear translocation of Bach2, IgG, CD40 and CD86 obviously up-regulated in B cells from SLE patients. Bach2 significantly inhibited the proliferation, promoted apoptosis of B cells from SLE patients, whereas BCR-ABL dramatically reversed cell changes induced by Bach2. Besides, BCR-ABL also inhibited nuclear translocation of Bach2 in B cells from SLE patients. Further, LY294002 treatment had no effect on decreased expression of Bach2 induced by BCR-ABL, but significantly eliminated BCR-ABL-induced phosphorylation of Bach2 and restored reduced nuclear translocation of Bach2 induced by BCR-ABL in B cells from SLE. CONCLUSIONS: Bach2 may play a suppressive role in B cells from SLE, and BCR-ABL may inhibit the nuclear translocation of Bach2 via serine phosphorylation through the PI3K pathway.

Increased expression of PRKCB mRNA in peripheral blood mononuclear cells from patients with systemic lupus erythematosus.

The polymorphism of PRKCB has been proven to be associated with systemic lupus erythematosus (SLE) in our previous study. We aimed to investigate the relationship between expression of PRKCB mRNA and the Disease Activity Index (SLEDAI) and manifestations of SLE. Quantitative reverse transcription polymerase chain reaction (RT-PCR) was applied to examine the expression of PRKCB mRNA in peripheral blood mononuclear cells of 60 patients with SLE and 62 controls. The Sequenom MassArray System was used to detect genotype SNP rs16972959. The expression levels of PRKCB mRNA in SLE cases were significantly increased compared with those in healthy controls (P < 0.001). In addition, PRKCB mRNA expression levels were negatively correlated with the SLEDAI (P < 0.05, r = -0.322), with lower mRNA expression levels of PRKCB in patients found with higher SLEDAI, presence of a new rash (P < 0.01), and proteinuria (P < 0.05). No association evidence was observed between the genotype of the variant rs16972959 and PRKCB mRNA expression levels; however, SNP rs16972959 was found to be an expression quantitative trait loci for PRKCB with the SLE risk allele correlated with increased expression in naïve monocytes (PFDR  = 9.12 × 10-13 ) and stimulated monocytes (9.24 × 10-6  > PFDR  > 2.75 × 10-17 ). On the other hand, SNP rs16972959 of PRKCB was found to have suggestive significant associations with vasculitis (P = 0.00718) of SLE. These results indicated that expression of PRKCB mRNA may be correlated with the pathogenesis of SLE; however, more investigation is still needed.

Identification of ST3AGL4, MFHAS1, CSNK2A2 and CD226 as loci associated with systemic lupus erythematosus (SLE) and evaluation of SLE genetics in drug repositioning.

OBJECTIVES: Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with a strong genetic component in its pathogenesis. Through genome-wide association studies (GWAS), we recently identified 10 novel loci associated with SLE and uncovered a number of suggestive loci requiring further validation. This study aimed to validate those loci in independent cohorts and evaluate the role of SLE genetics in drug repositioning. METHODS: We conducted GWAS and replication studies involving 12 280 SLE cases and 18 828 controls, and performed fine-mapping analyses to identify likely causal variants within the newly identified loci. We further scanned drug target databases to evaluate the role of SLE genetics in drug repositioning. RESULTS: We identified three novel loci that surpassed genome-wide significance, including ST3AGL4 (rs13238909, pmeta=4.40E-08), MFHAS1 (rs2428, pmeta=1.17E-08) and CSNK2A2 (rs2731783, pmeta=1.08E-09). We also confirmed the association of CD226 locus with SLE (rs763361, pmeta=2.45E-08). Fine-mapping and functional analyses indicated that the putative causal variants in CSNK2A2 locus reside in an enhancer and are associated with expression of CSNK2A2 in B-lymphocytes, suggesting a potential mechanism of association. In addition, we demonstrated that SLE risk genes were more likely to be interacting proteins with targets of approved SLE drugs (OR=2.41, p=1.50E-03) which supports the role of genetic studies to repurpose drugs approved for other diseases for the treatment of SLE. CONCLUSION: This study identified three novel loci associated with SLE and demonstrated the role of SLE GWAS findings in drug repositioning.

Exome-wide association study identifies four novel loci for systemic lupus erythematosus in Han Chinese population.

OBJECTIVES: Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of considerable genetic predisposition. Genome-wide association studies have identified tens of common variants for SLE. However, the majority of them reside in non-coding sequences. The contributions of coding variants have not yet been systematically evaluated. METHODS: We performed a large-scale exome-wide study in 5004 SLE cases and 8179 healthy controls in a Han Chinese population using a custom exome array, and then genotyped 32 variants with suggestive evidence in an independent cohort of 13 246 samples. We further explored the regulatory effect of one novel non-coding single nucleotide polymorphism (SNP) in ex vivo experiments. RESULTS: We discovered four novel SLE gene regions (LCT, TPCN2, AHNAK2 and TNFRSF13B) encompassing three novel missense variants (XP_016859577.1:p.Asn1639Ser, XP_016859577.1:p.Val219Phe and XP_005267356.1:p.Thr4664Ala) and two non-coding variants (rs10750836 and rs4792801) with genome-wide significance (pmeta <5.00×10-8). These variants are enriched in several chromatin states of primary B cells. The novel intergenic variant rs10750836 exhibited an expression quantitative trait locus effect on the TPCN2 gene in immune cells. Clones containing this novel SNP exhibited gene promoter activity for TPCN2 (P=1.38×10-3) whose expression level was reduced significantly in patients with SLE (P<2.53×10-2) and was suggested to be further modulated by rs10750836 in CD19+ B cells (P=7.57×10-5) in ex vivo experiments. CONCLUSIONS: This study identified three novel coding variants and four new susceptibility gene regions for SLE. The results provide insights into the biological mechanism of SLE.

NFKB1 mediates Th1/Th17 activation in the pathogenesis of psoriasis.

This study was aimed to investigate whether NFKB1 participates in the pathogenesis of psoriasis by mediating Th1/Th17 cells. In this study, expression of NFKB1 was assessed in skin tissues from psoriasis patients and the healthy controls through Western blot and Immunohistochemistry. Enzyme-linked immunosorbent assay (ELISA) was used to analyze the serum levels of IFN-γ, IL-17 (IL-17A) and IL-17RA. The imiquimod-induced psoriasis mouse model was employed to examine the role of NFKB1 in psoriasis via the assessment of psoriasis area and severity index (PASI), including erythema, thickness and scales. The effects of NFKB1 on Th1/Th17 cells in were examined by flow cytometry. In vitro co-culture of Th1/Th17 cells isolated from different mice with HaCat cells was conducted to elucidate the effect of Th1/Th17 cells-mediated by NFKB1 on HaCat cells by MTT, wound healing and transwell invasion assay, respectively. The results showed that NF-κB p105/p50 expression in skin tissues was significantly increased in psoriasis (n = 21) compared to the healthy controls (n = 16), as well as levels of serum INF-γ and IL-17. Additionally, NF-κB p105/p50 expression in lesional skin tissues was much higher than that in non-lesional skin tissues of the same patients. In the psoriasis mouse model, NFKB1 overexpression significantly elevated the scores of erythema, thickness and scales. Besides, NFKB1 up-regulated the level of NF-κB p105/p50, INF-γ, T-bet, IL-17 and RORγt, as well as Th1/Th17 cells in skin tissues of psoriasis mice. Finally, in vitro assay confirmed that the activation of Th1 and Th17 mediated by NFKB1 in psoriasis promoted the proliferation, migration and invasion of keratinocytes. These findings suggest a critical role for NFKB1 in the regulation of Th1 and Th17 in psoriasis.

A Genetic Variant rs1020760at NFKB1 is Associated with Clinical Features of Psoriasis Vulgaris in a Han Chinese Population.

Psoriasis vulgaris is a chronic inflammatory skin disease associated with complex genetic susceptibility. Recently, we identified a single-nucleotide variant rs1020760 at NFKB1 significantly associated with psoriasis in a Han Chinese population in deep analysis of exome and targeted sequencing (P = 1.76 × 10(-8) ). To investigate the potential association between rs1020760 and phenotypes of psoriasis vulgaris, we performed a genotype-phenotype analysis. A total of 9946 cases and 9906 controls with detailed clinical and demographic information were involved in this study, while the genotype data of rs1020760 was available in the previous targeted sequencing study of psoriasis. Genotype-based association testing revealed the additive model might provide the best fit for rs1020760 (P = 5.44 × 10(-8) ). Case-only analysis showed that the distribution of allele G was significantly different between the cases with and without family history (Pallele = 4.07 × 10(-3) ,Pgenotype = 5.75 × 10(-3) ). The differences in allele and genotype frequencies were observed between all the subphenotypes and controls except for the genotype frequency of the late onset subgroup, while no difference was found in case-only analysis for the other two subphenotypes. Rs1020760 was preferentially associated with family history of psoriasis, implying that NFKB1 might not only play important roles in the development of psoriasis, but might also contribute to the special phenotypes of this disease.

Association analysis revealed one susceptibility locus for vitiligo with immune-related diseases in the Chinese Han population.

Generalized vitiligo is an autoimmune disease characterized by melanocyte loss, which results in patchy depigmentation of skin and hair, and is associated with an elevated risk of other immune-related diseases. However, there is no reported study on the associations between immune susceptibility polymorphisms and the risk of vitiligo with immune-related diseases. The aim of this study was to evaluate the potential influence of 10 single-nucleotide polymorphisms (SNPs) at 18q21.31 (rs10503019), 4p16.1 (rs11940117), 3q28 (rs1464510), 14q12 (rs2273844), 12q13.2 (rs2456973), 16q12.2 (rs3213758), 10q25.3 (rs4353229), 3q13.33 (rs59374417), and 10p15.1 (rs706779 and rs7090530) on vitiligo with immune-related diseases in the Chinese Han population. All SNPs were genotyped in 552 patients with vitiligo-associated immune-related diseases and 1656 controls using the Sequenom MassArray system. Data were analyzed with PLINK 1.07 software. The C allele of rs2456973 at 12q13.2 was observed to be significantly associated with vitiligo-associated immune-related diseases (autoimmune diseases and allergic diseases) (P = 0.0028, odds ratio (OR) = 1.27). In subphenotype analysis, the rs2456973 C allele was also significantly associated with early-onset vitiligo by comparing with controls (P = 0.0001) and in the case-only analysis (P = 0.0114). We confirmed that 12q13.2 was an important candidate locus for vitiligo with immune-related diseases (autoimmune diseases and allergic diseases) and affected disease phenotypes with early onset.

Exome sequencing identifies SLC17A9 pathogenic gene in two Chinese pedigrees with disseminated superficial actinic porokeratosis.

BACKGROUND: Disseminated superficial actinic porokeratosis (DSAP) is a rare autosomal dominant genodermatosis characterised by annular lesions that has an atrophic centre and a prominent peripheral ridge distributed on sun exposed area. It exhibits high heterogeneity, and five linkage loci have been reported. The mevalonate kinase (MVK) gene located on 12q24 has been confirmed as one of the disease-causing genes. But, the pathogenesis of a large part of DSAP remains unclear so far. METHODS: The recruited with DSAP carried no MVK coding mutations. Exome sequencing was performed in two affected and one unaffected individual in Family 1. Cosegregation of the candidate variants was tested in other family members. Sanger sequencing in 33 individuals with familial DSAP and 19 sporadic DSAP individuals was performed for validating the causative gene. RESULTS: An average of 1.35×10(5) variants were generated from exome data and 133 novel NS/SS/indels were identified as being shared by two affected individuals but absent in the unaffected individual. After functional prediction, 25 possible deleterious variants were identified. In Family 1, a missense variant c.932G>A (p.Arg311Gln) in exon 10 of SLC17A9 was observed in cosegregation with the phenotype; this amino acid substitution was located in a highly conserved major facilitator superfamily (MFS) domain in multiple mammalian. One additional missense variant c.25C>T (p.Arg9Cys) in exon 2 of SLC17A9 was found in Family 2. CONCLUSIONS: The result identified SLC17A9 as another pathogenic gene for DSAP, which suggests a correlation between the aberrant vesicular nucleotide transporter and the pathogenesis of DSAP.

Smoking behaviors before and after implementation of a smoke-free legislation in Guangzhou, China.

BACKGROUND: According to the partial smoke-free legislation implemented on 1 September 2010 in Guangzhou, China, smoke-free did not cover all indoor areas. Some places have a full smoking ban (100 % smoke-free), other places have a partial smoking ban, and homes have no ban. This study aimed to compare the smoking behaviors before and after implementation of a smoke-free legislation. METHOD: A repeated cross-sectional survey was conducted on smoking-related behaviors with a total of 4,900 respondents before, and 5,135 respondents after the legislation was instituted. For each wave of the survey, a three-stage stratified sampling process was used to obtain a representative sample. Pearson's Chi-square test was used to determine differences of smoking prevalence and quit ratio between the two samples. Logistic regression models were used to examine the associations of a smoke-free legislation with smoking behaviors. RESULTS: The overall daily smoking rate declined significantly from 20.8 % to 18.2 % (p < 0.05), especially among those aged 15-24 years. The quit ratios increased significantly (from 14.5 % to 17.9 %), but remained low among 15-44 year olds. The overall self-reported smoking behaviors in locations with a full smoking ban decreased significantly from 36.4 % to 24.3 % with the greater drops occurring in cultural venues, public transport vehicles, and government offices. Smoking in places with partial smoking bans remained high (89.6 % vs. 90.4 %), although a slight decrease was observed in some of these areas. The implementation of a smoke-free legislation did not lead to more smoking in homes (91.0 % vs 89.4 %), but smoking in homes remained high. CONCLUSIONS: These findings highlight the urgent need for a comprehensive smoke-free legislation covering all public places in Guangzhou, simultaneously educational interventions and campaigns promoting voluntary changes in home smoking need to occur.

Templates are available to model nearly all complexes of structurally characterized proteins.

Traditional approaches to protein-protein docking sample the binding modes with no regard to similar experimentally determined structures (templates) of protein-protein complexes. Emerging template-based docking approaches utilize such similar complexes to determine the docking predictions. The docking problem assumes the knowledge of the participating proteins' structures. Thus, it provides the possibility of aligning the structures of the proteins and the template complexes. The progress in the development of template-based docking and the vast experience in template-based modeling of individual proteins show that, generally, such approaches are more reliable than the free modeling. The key aspect of this modeling paradigm is the availability of the templates. The current common perception is that due to the difficulties in experimental structure determination of protein-protein complexes, the pool of docking templates is insignificant, and thus a broad application of template-based docking is possible only at some future time. The results of our large scale, systematic study show that, surprisingly, in spite of the limited number of protein-protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the known protein-protein interactions, provided the components themselves have a known structure or can be homology-built. About one-third of the templates are of good quality when they are compared to experimental structures in test sets extracted from the Protein Data Bank and would be useful starting points in modeling the complexes. This finding dramatically expands our ability to model protein interactions, and has far-reaching implications for the protein docking field in general.

MGAviewer: a desktop visualization tool for analysis of metagenomics alignment data.

SUMMARY: Numerous metagenomics projects have produced tremendous amounts of sequencing data. Aligning these sequences to reference genomes is an essential analysis in metagenomics studies. Large-scale alignment data call for intuitive and efficient visualization tool. However, current tools such as various genome browsers are highly specialized to handle intraspecies mapping results. They are not suitable for alignment data in metagenomics, which are often interspecies alignments. We have developed a web browser-based desktop application for interactively visualizing alignment data of metagenomic sequences. This viewer is easy to use on all computer systems with modern web browsers and requires no software installation. AVAILABILITY: http://weizhongli-lab.org/mgaviewer

Effectiveness of guselkumab in patients with moderate-to-severe plaque psoriasis: a retrospective study from China.

Data on guselkumab as treatment for moderate-to-severe plaque psoriasis, especially in different body regions, in China is limited. This study aimed to estimate the effectiveness of guselkumab in Chinese patients with moderate-to-severe plaque psoriasis, including effectiveness at different body regions. This multicentre, observational study retrospectively enrolled patients with moderate-to-severe plaque psoriasis. Effectiveness outcome was based on Psoriasis Area and Severity Index (PASI) response and improvement in Body Surface Area (BSA) and Dermatology Life Quality Index (DLQI). A total of 51 patients were included, with a median age of 44.00 (18.00, 74.00) years and median duration of psoriasis of 10.00 (0.50, 55.00) years. After 20 weeks of treatment, PASI response with 75% improvement from baseline (PASI 75) was reported in 96.1% of patients; 72.5% of patients achieved a DLQI score of 0-1 at week 20. The percentage of affected BSA was significantly decreased at week 4 (p<0.05), week 12 (p<0.001) and week 20 (p<0.001). PASI score significantly changed from baseline after four weeks (p<0.001), 12 weeks (p<0.001) and 20 weeks of treatment (p<0.001). DLQI score significantly increased at week 4 (p<0.001), week 12 (p<0.001) and week 20 (p<0.001). PASI 75 was achieved for the upper limbs in all cases and 100% PASI improvement (PASI 100) in 89.1%. The head and lower limbs were the areas least responsive to treatment, with PASI 100 achieved in only 68.6% and 70.6%, respectively. Guselkummab provided rapid and sustained PASI improvement, especially for the skin of the upper limbs and body trunk.

CD-HIT: accelerated for clustering the next-generation sequencing data.

SUMMARY: CD-HIT is a widely used program for clustering biological sequences to reduce sequence redundancy and improve the performance of other sequence analyses. In response to the rapid increase in the amount of sequencing data produced by the next-generation sequencing technologies, we have developed a new CD-HIT program accelerated with a novel parallelization strategy and some other techniques to allow efficient clustering of such datasets. Our tests demonstrated very good speedup derived from the parallelization for up to ∼24 cores and a quasi-linear speedup for up to ∼8 cores. The enhanced CD-HIT is capable of handling very large datasets in much shorter time than previous versions. AVAILABILITY: http://cd-hit.org. CONTACT: liwz@sdsc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.