Common genetic variants associated with urinary phthalate levels in children: A genome-wide study.

INTRODUCTION: Phthalates, or dieters of phthalic acid, are a ubiquitous type of plasticizer used in a variety of common consumer and industrial products. They act as endocrine disruptors and are associated with increased risk for several diseases. Once in the body, phthalates are metabolized through partially known mechanisms, involving phase I and phase II enzymes. OBJECTIVE: In this study we aimed to identify common single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) associated with the metabolism of phthalate compounds in children through genome-wide association studies (GWAS). METHODS: The study used data from 1,044 children with European ancestry from the Human Early Life Exposome (HELIX) cohort. Ten phthalate metabolites were assessed in a two-void pooled urine collected at the mean age of 8 years. Six ratios between secondary and primary phthalate metabolites were calculated. Genome-wide genotyping was done with the Infinium Global Screening Array (GSA) and imputation with the Haplotype Reference Consortium (HRC) panel. PennCNV was used to estimate copy number variants (CNVs) and CNVRanger to identify consensus regions. GWAS of SNPs and CNVs were conducted using PLINK and SNPassoc, respectively. Subsequently, functional annotation of suggestive SNPs (p-value < 1E-05) was done with the FUMA web-tool. RESULTS: We identified four genome-wide significant (p-value < 5E-08) loci at chromosome (chr) 3 (FECHP1 for oxo-MiNP_oh-MiNP ratio), chr6 (SLC17A1 for MECPP_MEHHP ratio), chr9 (RAPGEF1 for MBzP), and chr10 (CYP2C9 for MECPP_MEHHP ratio). Moreover, 115 additional loci were found at suggestive significance (p-value < 1E-05). Two CNVs located at chr11 (MRGPRX1 for oh-MiNP and SLC35F2 for MEP) were also identified. Functional annotation pointed to genes involved in phase I and phase II detoxification, molecular transfer across membranes, and renal excretion. CONCLUSION: Through genome-wide screenings we identified known and novel loci implicated in phthalate metabolism in children. Genes annotated to these loci participate in detoxification, transmembrane transfer, and renal excretion.

An encyclopedia of enhancer-gene regulatory interactions in the human genome.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

Treatment with MOG-DNA vaccines induces CD4+CD25+FoxP3+ regulatory T cells and up-regulates genes with neuroprotective functions in experimental autoimmune encephalomyelitis.

BACKGROUND: DNA vaccines represent promising therapeutic strategies in autoimmune disorders such as multiple sclerosis (MS). However, the precise mechanisms by which DNA vaccines induce immune regulation remain largely unknown. Here, we aimed to expand previous knowledge existing on the mechanisms of action of DNA vaccines in the animal model of MS, experimental autoimmune encephalomyelitis (EAE), by treating EAE mice with a DNA vaccine encoding the myelin oligodendrocyte glycoprotein (MOG), and exploring the therapeutic effects on the disease-induced inflammatory and neurodegenerative changes. METHODS: EAE was induced in C57BL6/J mice by immunization with MOG35₋55 peptide. Mice were intramuscularly treated with a MOG-DNA vaccine or vehicle in prophylactic and therapeutic approaches. Histological studies were performed in central nervous system (CNS) tissue. Cytokine production and regulatory T cell (Treg) quantification were achieved by flow cytometry. Gene expression patterns were determined using microarrays, and the main findings were validated by real-time PCR. RESULTS: MOG-DNA treatment reduced the clinical and histopathological signs of EAE when administered in both prophylactic and therapeutic settings. Suppression of clinical EAE was associated with dampening of antigen (Ag)-specific proinflammatory Th1 and Th17 immune responses and, interestingly, expansion of Treg in the periphery and upregulation in the CNS of genes encoding neurotrophic factors and proteins involved in remyelination. CONCLUSIONS: These results suggest for the first time that the beneficial effects of DNA vaccines in EAE are not limited to anti-inflammatory mechanisms, and DNA vaccines may also exert positive effects through hitherto unknown neuroprotective mechanisms.

Tyrosine kinase 2 variant influences T lymphocyte polarization and multiple sclerosis susceptibility.

The tyrosine kinase 2 variant rs34536443 has been established as a genetic risk factor for multiple sclerosis in a variety of populations. However, the functional effect of this variant on disease pathogenesis remains unclear. This study replicated the genetic association of tyrosine kinase 2 with multiple sclerosis in a cohort of 1366 French patients and 1802 controls. Furthermore, we assessed the functional consequences of this polymorphism on human T lymphocytes by comparing the reactivity and cytokine profile of T lymphocytes isolated from individuals expressing the protective TYK2(GC) genotype with the disease-associated TYK2(GG) genotype. Our results demonstrate that the protective C allele infers decreased tyrosine kinase 2 activity, and this reduction of activity is associated with a shift in the cytokine profile favouring the secretion of Th2 cytokines. These findings suggest that the rs34536443 variant effect on multiple sclerosis susceptibility might be mediated by deviating T lymphocyte differentiation toward a Th2 phenotype. This impact of tyrosine kinase 2 on effector differentiation is likely to be of wider importance because other autoimmune diseases also have been associated with polymorphisms within tyrosine kinase 2. The modulation of tyrosine kinase 2 activity might therefore represent a new therapeutic approach for the treatment of autoimmune diseases.

PLANdbAffy: probe-level annotation database for Affymetrix expression microarrays.

Standard Affymetrix technology evaluates gene expression by measuring the intensity of mRNA hybridization with a panel of the 25-mer oligonucleotide probes, and summarizing the probe signal intensities by a robust average method. However, in many cases, signal intensity of the probe does not correlate with gene expression. This could be due to the hybridization of the probe to a transcript of another gene, mapping of the probe to an intron, alternative splicing, single nucleotide polymorphisms and other reasons. We have developed a database, PLANdbAffy (available at http://affymetrix2.bioinf.fbb.msu.ru), that contains the results of the alignment of probe sequences from five Affymetrix expression microarrays to the human genome. We have determined the probes matching the transcript-coding regions in the correct orientation. For each such probe alignment region, we determined the mRNA and EST sequences that contain the probe sequence. In the textual part of the database interface we summarize the data on the sequences that cover the probe alignment region and SNPs that are located inside it. The graphical part of our database interface is implemented as custom tracks to the UCSC genome browser that allows one to utilize all the data that are offered by UCSC browser.

Positive selection in alternatively spliced exons of human genes.

Alternative splicing is a well-recognized mechanism of accelerated genome evolution. We have studied single-nucleotide polymorphisms and human-chimpanzee divergence in the exons of 6672 alternatively spliced human genes, with the aim of understanding the forces driving the evolution of alternatively spliced sequences. Here, we show that alternatively spliced exons and exon fragments (alternative exons) from minor isoforms experience lower selective pressure at the amino acid level, accompanied by selection against synonymous sequence variation. The results of the McDonald-Kreitman test suggest that alternatively spliced exons, unlike exons constitutively included in the mRNA, are also subject to positive selection, with up to 27% of amino acids fixed by positive selection.

Rodent-specific alternative exons are more frequent in rapidly evolving genes and in paralogs.

BACKGROUND: Alternative splicing is an important mechanism for generating functional and evolutionary diversity of proteins in eukaryotes. Here, we studied the frequency and functionality of recently gained, rodent-specific alternative exons. RESULTS: We projected the data about alternative splicing of mouse genes to the rat, human, and dog genomes, and identified exons conserved in the rat genome, but missing in more distant genomes. We estimated the frequency of rodent-specific exons while controlling for possible residual conservation of spurious exons. The frequency of rodent-specific exons is higher among predominantly skipped exons and exons disrupting the reading frame. Separation of all genes by the rate of sequence evolution and by gene families has demonstrated that rodent-specific cassette exons are more frequent in rapidly evolving genes and in rodent-specific paralogs. CONCLUSION: Thus we demonstrated that recently gained exons tend to occur in fast-evolving genes, and their inclusion rate tends to be lower than that of older exons. This agrees with the theory that gain of alternative exons is one of the major mechanisms of gene evolution.

Conserved and species-specific alternative splicing in mammalian genomes.

BACKGROUND: Alternative splicing has been shown to be one of the major evolutionary mechanisms for protein diversification and proteome expansion, since a considerable fraction of alternative splicing events appears to be species- or lineage-specific. However, most studies were restricted to the analysis of cassette exons in pairs of genomes and did not analyze functionality of the alternative variants. RESULTS: We analyzed conservation of human alternative splice sites and cassette exons in the mouse and dog genomes. Alternative exons, especially minor-isofom ones, were shown to be less conserved than constitutive exons. Frame-shifting alternatives in the protein-coding regions are less conserved than frame-preserving ones. Similarly, the conservation of alternative sites is highest for evenly used alternatives, and higher when the distance between the sites is divisible by three. The rate of alternative-exon and site loss in mouse is slightly higher than in dog, consistent with faster evolution of the former. The evolutionary dynamics of alternative sites was shown to be consistent with the model of random activation of cryptic sites. CONCLUSION: Consistent with other studies, our results show that minor cassette exons are less conserved than major-alternative and constitutive exons. However, our study provides evidence that this is caused not only by exon birth, but also lineage-specific loss of alternative exons and sites, and it depends on exon functionality.

Overlapping alternative donor splice sites in the human genome.

Over 50% of donor splice sites in the human genome have a potential alternative donor site at a distance of three to six nucleotides. Conservation of these potential sites is determined by the consensus requirements and by its exonic or intronic location. Several hundred pairs of overlapping sites are confirmed to be alternatively spliced as both sites in a pair are supported by a protein, by a full-length mRNA, or by expressed sequence tags (ESTs) from at least two independent clone libraries. Overlapping sites may clash with consensus requirements. Pairs with a site shift of four nucleotides are the most abundant, despite the frameshift in the protein-coding region that they introduce. The site usage in pairs is usually uneven, and the major site is more frequently conserved in other mammalian genomes. Overlapping alternative donor sites and acceptor sites may have different functional roles: alternative splicing of overlapping acceptor sites leads mainly to microvariations in protein sequences; whereas alternative donor sites often lead to frameshifts and thus either yield major differences in the protein sequence and structure, or generate nonsense-mediated decay-inducing mRNA isoforms likely involved in regulated unproductive splicing pathways.

Fast rate of evolution in alternatively spliced coding regions of mammalian genes.

BACKGROUND: At least half of mammalian genes are alternatively spliced. Alternative isoforms are often genome-specific and it has been suggested that alternative splicing is one of the major mechanisms for generating protein diversity in the course of evolution. Another way of looking at alternative splicing is to consider sequence evolution of constitutive and alternative regions of protein-coding genes. Indeed, it turns out that constitutive and alternative regions evolve in different ways. RESULTS: A set of 3029 orthologous pairs of human and mouse alternatively spliced genes was considered. The rate of nonsynonymous substitutions (dN), the rate of synonymous substitutions (dS), and their ratio (omega = dN/dS) appear to be significantly higher in alternatively spliced coding regions compared to constitutive regions. When N-terminal, internal and C-terminal alternatives are analysed separately, C-terminal alternatives appear to make the main contribution to the observed difference. The effects become even more pronounced in a subset of fast evolving genes. CONCLUSION: These results provide evidence of weaker purifying selection and/or stronger positive selection in alternative regions and thus one more confirmation of accelerated evolution in alternative regions. This study corroborates the theory that alternative splicing serves as a testing ground for molecular evolution.

No statistical support for correlation between the positions of protein interaction sites and alternatively spliced regions.

BACKGROUND: Alternative splicing is an efficient mechanism for increasing the variety of functions fulfilled by proteins in a living cell. It has been previously demonstrated that alternatively spliced regions often comprise functionally important and conserved sequence motifs. The objective of this work was to test the hypothesis that alternative splicing is correlated with contact regions of protein-protein interactions. RESULTS: Protein sequence spans involved in contacts with an interaction partner were delineated from atomic structures of transient interaction complexes and juxtaposed with the location of alternatively spliced regions detected by comparative genome analysis and spliced alignment. The total of 42 alternatively spliced isoforms were identified in 21 amino acid chains involved in biomolecular interactions. Using this limited dataset and a variety of sophisticated counting procedures we were not able to establish a statistically significant correlation between the positions of protein interaction sites and alternatively spliced regions. CONCLUSIONS: This finding contradicts a naïve hypothesis that alternatively spliced regions would correlate with points of contact. One possible explanation for that could be that all alternative splicing events change the spatial structure of the interacting domain to a sufficient degree to preclude interaction. This is indirectly supported by the observed lack of difference in the behaviour of relatively short regions affected by alternative splicing and cases when large portions of proteins are removed. More structural data on complexes of interacting proteins, including structures of alternative isoforms, are needed to test this conjecture.

Baseline gene expression signatures in monocytes from multiple sclerosis patients treated with interferon-beta.

BACKGROUND: A relatively large proportion of relapsing-remitting multiple sclerosis (RRMS) patients do not respond to interferon-beta (IFNb) treatment. In previous studies with peripheral blood mononuclear cells (PBMC), we identified a subgroup of IFNb non-responders that was characterized by a baseline over-expression of type I IFN inducible genes. Additional mechanistic experiments carried out in IFNb non-responders suggested a selective alteration of the type I IFN signaling pathway in the population of blood monocytes. Here, we aimed (i) to investigate whether the type I IFN signaling pathway is up-regulated in isolated monocytes from IFNb non-responders at baseline; and (ii) to search for additional biological pathways in this cell population that may be implicated in the response to IFNb treatment. METHODS: Twenty RRMS patients classified according to their clinical response to IFNb treatment and 10 healthy controls were included in the study. Monocytes were purified from PBMC obtained before treatment by cell sorting and the gene expression profiling was determined with oligonucleotide microarrays. RESULTS AND DISCUSSION: Purified monocytes from IFNb non-responders were characterized by an over-expression of type I IFN responsive genes, which confirms the type I IFN signature in monocytes suggested from previous studies. Other relevant signaling pathways that were up-regulated in IFNb non-responders were related with the mitochondrial function and processes such as protein synthesis and antigen presentation, and together with the type I IFN signaling pathway, may also be playing roles in the response to IFNb.

Search for specific biomarkers of IFNß bioactivity in patients with multiple sclerosis.

Myxovirus A (MxA), a protein encoded by the MX1 gene with antiviral activity, has proven to be a sensitive measure of IFNß bioactivity in multiple sclerosis (MS). However, the use of MxA as a biomarker of IFNß bioactivity has been criticized for the lack of evidence of its role on disease pathogenesis and the clinical response to IFNß. Here, we aimed to identify specific biomarkers of IFNß bioactivity in order to compare their gene expression induction by type I IFNs with the MxA, and to investigate their potential role in MS pathogenesis. Gene expression microarrays were performed in PBMC from MS patients who developed neutralizing antibodies (NAB) to IFNß at 12 and/or 24 months of treatment and patients who remained NAB negative. Nine genes followed patterns in gene expression over time similar to the MX1, which was considered the gold standard gene, and were selected for further experiments: IFI6, IFI27, IFI44L, IFIT1, HERC5, LY6E, RSAD2, SIGLEC1, and USP18. In vitro experiments in PBMC from healthy controls revealed specific induction of selected biomarkers by IFNß but not IFNγ, and several markers, in particular USP18 and HERC5, were shown to be significantly induced at lower IFNß concentrations and more selective than the MX1 as biomarkers of IFNß bioactivity. In addition, USP18 expression was deficient in MS patients compared with healthy controls (p = 0.0004). We propose specific biomarkers that may be considered in addition to the MxA to evaluate IFNß bioactivity, and to further explore their implication in MS pathogenesis.

Probe-Level Universal Search (PLUS) algorithm for gender differentiation in affymetrix datasets.

Affymetrix microarrays measure gene expression based on the intensity of hybridization of a panel of oligonucleotide probes (probe set) with mRNA. The signals from all probes within a probe set are converted into a single measure that represents the expression value of a gene. This step diminishes the number of independently measured parameters and eliminates from consideration individual "good-working" probes. We propose a new feature selection algorithm (Probe Level Universal Search or PLUS algorithm) for probe-level analysis of gene expression datasets. The algorithm evaluates the intensities of perfect-match Affymetrix probes individually and selects probes that allow one to distinguish two given classes of samples. The algorithm was used to differentiate the samples according to their gender ("gender differentiation"). The universal gender differentiating set of 3' Gene Affymetrix microarray probes was selected; the set consists of 38 probes from XIST gene of X-chromosome and 17 probes from five Y-chromosome genes: RPS4Y1, EIF1A, DDX3Y, JARID1D and USP9Y. The selection procedure based on the probes selected by PLUS algorithm differentiates the sex chromosome karyotype of the sample, reveals samples with incorrect gender labels and samples from patients with hereditary syndromes or cancer-associated chromosome abnormalities.

Low conservation of alternative splicing patterns in the human and mouse genomes.

Alternative splicing has recently emerged as a major mechanism of generating protein diversity in higher eukaryotes. We compared alternative splicing isoforms of 166 pairs of orthologous human and mouse genes. As the mRNA and EST libraries of human and mouse are not complete and thus cannot be compared directly, we instead analyzed whether known cassette exons or alternative splicing sites from one genome are conserved in the other genome. We demonstrate that about half of the analyzed genes have species-specific isoforms, and about a quarter of elementary alternatives are not conserved between the human and mouse genomes. The detailed results of this study are available at www.ig-msk.ru:8005/HMG_paper.