Diverse selective regimes shape genetic diversity at ADAR genes and at their coding targets.

A-to-I RNA editing operated by ADAR enzymes is extremely common in mammals. Several editing events in coding regions have pivotal physiological roles and affect protein sequence (recoding events) or function. We analyzed the evolutionary history of the 3 ADAR family genes and of their coding targets. Evolutionary analysis indicated that ADAR evolved adaptively in primates, with the strongest selection in the unique N-terminal domain of the interferon-inducible isoform. Positively selected residues in the human lineage were also detected in the ADAR deaminase domain and in the RNA binding domains of ADARB1 and ADARB2. During the recent history of human populations distinct variants in the 3 genes increased in frequency as a result of local selective pressures. Most selected variants are located within regulatory regions and some are in linkage disequilibrium with eQTLs in monocytes. Finally, analysis of conservation scores of coding editing sites indicated that editing events are counter-selected within regions that are poorly tolerant to change. Nevertheless, a minority of recoding events occurs at highly conserved positions and possibly represents the functional fraction. These events are enriched in pathways related to HIV-1 infection and to epidermis/hair development. Thus, both ADAR genes and their targets evolved under variable selective regimes, including purifying and positive selection. Pressures related to immune response likely represented major drivers of evolution for ADAR genes. As for their coding targets, we suggest that most editing events are slightly deleterious, although a minority may be beneficial and contribute to antiviral response and skin homeostasis.

Data on cardiac lncRNA STX18-AS1 expression in developing human hearts and function during in vitro hESC-cardiomyocyte differentiation.

This article presents data concerning STX18-AS1, a long noncoding RNA gene identified from a Genome-wide association study of Atrial Septal Defect (ASD). The data describes its expression patterns in human tissues and functions in regulating cardiomyocyte differentiation in vitro. STX18-AS1 is a lncRNA with a higher abundance in developing tissues, including hearts. Its transcription distribution within the embryonic hearts during key heart septation stages supports STX18-AS1's association with risk SNPs for ASD. The CRISPR stem cell pool in which STX18-AS1 was knocked down, showed reduced CM differentiation efficiency and lower expression of key cardiac transcriptional factors. This indicated its regulative role in supporting the lineage specification from cardiac mesoderm into cardiac progenitors and cardiomyocytes. These data can benefit the understanding of human embryonic heart developmental biology, and the time-course changes of cardiac transcriptional factors during in vitro cardiomyocyte differentiation from human embryonic stem cells.

Mendelian randomization analysis identified genes potentially pleiotropically associated with periodontitis.

OBJECTIVE: To prioritize genes that were pleiotropically or potentially causally associated with periodontitis. METHODS: We applied the summary data-based Mendelian randomization (SMR) method integrating genome-wide association study (GWAS) for periodontitis and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with periodontitis. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analysis were done for participants of European and East Asian ancestries, separately. RESULTS: We identified multiple genes showing pleiotropic association with periodontitis in participants of European ancestry and participants of East Asian ancestry. PDCD2 (corresponding probe: ILMN_1758915) was the top hit showing pleotropic association with periodontitis in the participants of European ancestry using CAGE eQTL data, and BX093763 (corresponding probe: ILMN_1899903) and AC104135.3 (corresponding probe: ENSG00000204792.2) were the top hits in the participants of East Asian ancestry using CAGE eQTL data and GTEx eQTL data, respectively. CONCLUSION: We identified multiple genes that may be involved in the pathogenesis of periodontitis in participants of European ancestry and participants of East Asian ancestry. Our findings provided important leads to a better understanding of the mechanisms underlying periodontitis and revealed potential therapeutic targets for the effective treatment of periodontitis.

Dissecting the contribution of single nucleotide polymorphisms in CCR9 and CCL25 genomic regions to the celiac disease phenotype.

PURPOSE AND OBJECTIVES: Given their role in homing immune cells to the intestine, CC motif chemokine receptor 9 (CCR9) and its specific ligand CC motif chemokine ligand 25 (CCL25) are interesting candidate genes for celiac disease. These genes are located in regions previously shown to be associated with or linked to celiac disease, but no investigations on their association with various celiac disease phenotypes have so far been conducted. Here we studied such associations of both genotyped and imputed single nucleotide polymorphisms (SNPs) with either regulatory function or exonic location of the CCR9 and CCL25 loci. RESULTS: Exploiting a carefully phenotyped cohort of 625 celiac disease patients and 1817 non-celiac controls, we identified that multiple SNPs with predicted regulatory function (RegulomeDB score ≤3a and/or eQTL effect) located between 100 kB upstream and downstream of CCR9 and CCL25 are associated with celiac disease and/or selected phenotypes. Of the genotyped SNPs in the CCR9 loci, rs213360 with an eQTL effect on CCR9 expression in blood was associated with celiac disease and all investigated phenotypes except high HLA risk. Rs1545985 with an eQTL on CCR9 expression and rs7652331 and rs12493471, both with RegulomeDB score ≤3a, were all associated with gastrointestinal symptoms and malabsorption and the latter additionally with anemia. The genotyped CCL25 SNPs rs952444 and rs882951, with RegulomeDB scores 1d and 1f respectively and eQTL effect on CCL25 expression in small intestine, were associated with gastrointestinal symptoms and malabsorption. The CCL25 SNP rs2303165 identified in sequencing followed by imputation was associated with partial villous atrophy. However, the association did not pass the permutation based multiple testing correction (PEMP2 > 0.05). CONCLUSIONS: We conclude that SNPs in the region of CCR9 and CCL25 with predicted functional effect or exonic localization likely contribute only modestly to various celiac disease phenotypes.

RSAT variation-tools: An accessible and flexible framework to predict the impact of regulatory variants on transcription factor binding.

Gene regulatory regions contain short and degenerated DNA binding sites recognized by transcription factors (TFBS). When TFBS harbor SNPs, the DNA binding site may be affected, thereby altering the transcriptional regulation of the target genes. Such regulatory SNPs have been implicated as causal variants in Genome-Wide Association Study (GWAS) studies. In this study, we describe improved versions of the programs Variation-tools designed to predict regulatory variants, and present four case studies to illustrate their usage and applications. In brief, Variation-tools facilitate i) obtaining variation information, ii) interconversion of variation file formats, iii) retrieval of sequences surrounding variants, and iv) calculating the change on predicted transcription factor affinity scores between alleles, using motif scanning approaches. Notably, the tools support the analysis of haplotypes. The tools are included within the well-maintained suite Regulatory Sequence Analysis Tools (RSAT, http://rsat.eu), and accessible through a web interface that currently enables analysis of five metazoa and ten plant genomes. Variation-tools can also be used in command-line with any locally-installed Ensembl genome. Users can input personal collections of variants and motifs, providing flexibility in the analysis.

A Gene Module-Based eQTL Analysis Prioritizing Disease Genes and Pathways in Kidney Cancer.

Clear cell renal cell carcinoma (ccRCC) is the most common and most aggressive form of renal cell cancer (RCC). The incidence of RCC has increased steadily in recent years. The pathogenesis of renal cell cancer remains poorly understood. Many of the tumor suppressor genes, oncogenes, and dysregulated pathways in ccRCC need to be revealed for improvement of the overall clinical outlook of the disease. Here, we developed a systems biology approach to prioritize the somatic mutated genes that lead to dysregulation of pathways in ccRCC. The method integrated multi-layer information to infer causative mutations and disease genes. First, we identified differential gene modules in ccRCC by coupling transcriptome and protein-protein interactions. Each of these modules consisted of interacting genes that were involved in similar biological processes and their combined expression alterations were significantly associated with disease type. Then, subsequent gene module-based eQTL analysis revealed somatic mutated genes that had driven the expression alterations of differential gene modules. Our study yielded a list of candidate disease genes, including several known ccRCC causative genes such as BAP1 and PBRM1, as well as novel genes such as NOD2, RRM1, CSRNP1, SLC4A2, TTLL1 and CNTN1. The differential gene modules and their driver genes revealed by our study provided a new perspective for understanding the molecular mechanisms underlying the disease. Moreover, we validated the results in independent ccRCC patient datasets. Our study provided a new method for prioritizing disease genes and pathways.