The information bottleneck method for genome-wide association studies

In population studies, most current methods focus on identifying one outcome-related SNP at a time by testing for differences of genotype frequencies between disease and healthy groups or among different population groups. However, testing a great number of SNPs simultaneously has a problem of multiple testing and will give false-positive results. Although, this problem can be effectively dealt with through several approaches such as Bonferroni correction, permutation testing and false discovery rates, patterns of the joint effects by several genes, each with weak effect, might not be able to be determined. With the availability of high-throughput genotyping technology, searching for multiple scattered SNPs over the whole genome and modeling their joint effect on the target variable has become possible. Exhaustive search of all SNP subsets is computationally infeasible for millions of SNPs in a genome-wide study. Several effective feature selection methods combined with classification functions have been proposed to search for an optimal SNP subset among big data sets where the number of feature SNPs far exceeds the number of observations. ^ In this study, we take two steps to achieve the goal. First we selected 1000 SNPs through an effective filter method and then we performed a feature selection wrapped around a classifier to identify an optimal SNP subset for predicting disease. And also we developed a novel classification method-sequential information bottleneck method wrapped inside different search algorithms to identify an optimal subset of SNPs for classifying the outcome variable. This new method was compared with the classical linear discriminant analysis in terms of classification performance. Finally, we performed chi-square test to look at the relationship between each SNP and disease from another point of view. ^ In general, our results show that filtering features using harmononic mean of sensitivity and specificity(HMSS) through linear discriminant analysis (LDA) is better than using LDA training accuracy or mutual information in our study. Our results also demonstrate that exhaustive search of a small subset with one SNP, two SNPs or 3 SNP subset based on best 100 composite 2-SNPs can find an optimal subset and further inclusion of more SNPs through heuristic algorithm doesn't always increase the performance of SNP subsets. Although sequential forward floating selection can be applied to prevent from the nesting effect of forward selection, it does not always out-perform the latter due to overfitting from observing more complex subset states. ^ Our results also indicate that HMSS as a criterion to evaluate the classification ability of a function can be used in imbalanced data without modifying the original dataset as against classification accuracy. Our four studies suggest that Sequential Information Bottleneck(sIB), a new unsupervised technique, can be adopted to predict the outcome and its ability to detect the target status is superior to the traditional LDA in the study. ^ From our results we can see that the best test probability-HMSS for predicting CVD, stroke,CAD and psoriasis through sIB is 0.59406, 0.641815, 0.645315 and 0.678658, respectively. In terms of group prediction accuracy, the highest test accuracy of sIB for diagnosing a normal status among controls can reach 0.708999, 0.863216, 0.639918 and 0.850275 respectively in the four studies if the test accuracy among cases is required to be not less than 0.4. On the other hand, the highest test accuracy of sIB for diagnosing a disease among cases can reach 0.748644, 0.789916, 0.705701 and 0.749436 respectively in the four studies if the test accuracy among controls is required to be at least 0.4. ^ A further genome-wide association study through Chi square test shows that there are no significant SNPs detected at the cut-off level 9.09451E-08 in the Framingham heart study of CVD. Study results in WTCCC can only detect two significant SNPs that are associated with CAD. In the genome-wide study of psoriasis most of top 20 SNP markers with impressive classification accuracy are also significantly associated with the disease through chi-square test at the cut-off value 1.11E-07. ^ Although our classification methods can achieve high accuracy in the study, complete descriptions of those classification results(95% confidence interval or statistical test of differences) require more cost-effective methods or efficient computing system, both of which can't be accomplished currently in our genome-wide study. We should also note that the purpose of this study is to identify subsets of SNPs with high prediction ability and those SNPs with good discriminant power are not necessary to be causal markers for the disease.^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Involvement of {\it p53\/} in PUVA-induced skin cancers

A combination of psoralen and ultraviolet-A radiation, commonly referred to as "PUVA," is widely used in the treatment of psoriasis. However, PUVA treatment increases the risk of developing skin cancer in psoriasis patients and induces skin cancer in mice. It is, however unknown whether the increased incidence of skin cancer in PUVA treated psoriasis patients is due to the carcinogenic effects of PUVA therapy or due to an indirect effect such as immunosuppression, which can permit the growth of tumors induced by UVB radiation. In this study, we used the p53 tumor suppressor gene as a molecular marker to determine whether PUVA-induced mouse skin cancers contain unique mutations in p53 that are different from UV-induced mutations, and if so, determine whether skin cancers from PUVA treated patients have PUVA-type or UV-type p53 mutations. Since the DNA lesions induced by PUVA are quite different from those induced by UV, we hypothesize that p53 mutations induced by PUVA may also be different from those induced by UV.^ Analysis of PUVA-induced murine skin cancers for p53 mutations revealed that 14 of 15 (93%) missense mutations detected in these cancers were localized at 5$\sp\prime$-TA/5$\sp\prime$-TAT sites, potential sites of psoralen photoadditions. Mutations at these sequences are exceedingly rare in UV-induced murine skin cancers. In addition, PUVA-induced murine skin cancers did not contain UV signature (C $\to$ T or CC $\to$ TT transitions) mutations in p53. These results suggest that PUVA induces unique mutations in p53 that can be distinguished from those induced by UV.^ Next we determined whether SCCs arising in PUVA treated psoriasis patients have PUVA-type or UV-type p53 mutations. The results indicated that 16 of 25 (64%) missense p53 mutations detected in SCCs from PUVA treated patients were located at 5$\sp\prime$-TG, 5$\sp\prime$-TA and 5$\sp\prime$-TT sites, putative sites of psoralen photobinding. Interestingly, about 32% of p53 mutations detected in SCCs from PUVA treated patients had the UV signature. Taken together these results suggest that both PUVA and UVB play a role in the development of SCCs in psoriasis patients undergoing PUVA therapy. ^

Acantosis Nigricans

Fil: Parra, Viviana Gladys. Hospital Luis Lagomaggiore (Mendoza, Argentina). Servicio de Dermatología

Papel de la interleuquina 23 en la diferenciación hacia Th17

The incidence of inflammatory and autoimmune diseases has increased among developed countries in the past 30 years, creating a demand for the development of effective and economic therapies for these diseases. Interleukin 23 (IL-23) is a pro-inflammatory cytokine whose increased production has been shown to play a key role in the establishment and maintenance of inflammatory and autoimmune diseases in different murine models such as inflammatory bowel disease, psoriasis and experimental autoimmune encephalomyelitis. More importantly, increased levels of IL-23 have been found in biopsies from patients with Crohn’s disease and ulcerative colitis, and psoriasis. The pathological consequences of excessive IL-23 signalling have been linked to its ability to promote the production of interleukin 17 (IL-17), particularly in the subpopulation of CD4 T cells Th17. However, the precise molecular mechanisms by which IL-23 sustains the Th17 response and induces pathogenic effector functions in these cells remain largely unknown. The global objective of the experiments carried out in this work was to determine the effect of IL-23 on the proliferation, survival and IL-17 and interferon gamma (IFN-ɣ) production in Th17 cells. These experiments have shown that IL-23 does not promote proliferation or survival of in vitro generated Th17 cells, and that there is no difference in the production of IL -17 in the absence or presence of IL -23. The IL-23 receptor, like other cytokine receptors, lacks intrinsic enzymatic activity. Instead, IL-23 receptor associates with members of the Janus tyrosine kinase family (Jaks). Cytokine binding to a Jak-associated receptor triggers the activation of the Signal Transducers and Activators of Transcription (STAT) family of transcription factors. Previous work indicated that the IL-23 receptor complex is associated with the tyrosine kinases Jak2 and Tyk2 that promote STAT3 phosphorylation. Subsequent studies showed that IL23 activation of STAT3 induces the expression of the transcription factor RORγt, which is crucial for IL-17 production. This work has explored the IL-23 signalling cascade, determining the optimal conditions for STAT3 activation and demonstrating the activation of other transcription factors such as STAT4, STAT5 and STAT1 that contribute to IL-23-mediated signalling pathways.

An antagonistic vascular endothelial growth factor (VEGF) variant inhibits VEGF-stimulated receptor autophosphorylation and proliferation of human endothelial cells

Vascular endothelial growth factor (VEGF) is a potent mitogen with a unique specificity for endothelial cells and a key mediator of aberrant endothelial cell proliferation and vascular permeability in a variety of human pathological situations, such as tumor angiogenesis, diabetic retinopathy, rheumatoid arthritis, or psoriasis. VEGF is a symmetric homodimeric molecule with two receptor binding interfaces lying on each pole of the molecule. Herein we report on the construction and recombinant expression of an asymmetric heterodimeric VEGF variant with an intact receptor binding interface at one pole and a mutant receptor binding interface at the second pole of the dimer. This VEGF variant binds to VEGF receptors but fails to induce receptor activation. In competition experiments, the heterodimeric VEGF variant antagonizes VEGF-stimulated receptor autophosphorylation and proliferation of endothelial cells. A 15-fold excess of the heterodimer was sufficient to inhibit VEGF-stimulated endothelial cell proliferation by 50%, and a 100-fold excess resulted in an almost complete inhibition. By using a rational approach that is based on the structure of VEGF, we have shown the feasibility to construct a VEGF variant that acts as an VEGF antagonist.

A 12R-lipoxygenase in human skin: Mechanistic evidence, molecular cloning, and expression

A recognized feature of psoriasis and other proliferative dermatoses is accumulation in the skin of the unusual arachidonic acid metabolite, 12R-hydroxyeicosatetraenoic acid (12R-HETE). This hydroxy fatty acid is opposite in chirality to the product of the well-known 12S-lipoxygenase and heretofore in mammals is known only as a product of cytochrome P450s. Here we provide mechanistic evidence for a lipoxygenase route to 12R-HETE in human psoriatic tissue and describe a 12R-lipoxygenase that can account for the biosynthesis. Initially we demonstrated retention of the C-12 deuterium of octadeuterated arachidonic acid in its conversion to 12R-HETE in incubations of psoriatic scales, indicating the end product is not formed by isomerization from 12S-H(P)ETE via the 12-keto derivative. Secondly, analysis of product formed from [10R-3H] and [10S-3H]-labeled arachidonic acids revealed that 12R-HETE synthesis is associated with stereospecific removal of the pro-R hydrogen from the 10-carbon of arachidonate. This result is compatible with 12R-lipoxygenase-catalyzed formation of 12R-HETE and not with a P450-catalyzed route to 12R-HETE in psoriatic scales. We cloned a lipoxygenase from human keratinocytes; the cDNA and deduced amino acid sequences share ≤50% identity to other human lipoxygenases. This enzyme, when expressed in Hela cells, oxygenates arachidonic acid to 12-HPETE, >98% 12R in configuration. The 12R-lipoxygenase cDNA is detectable by PCR in psoriatic scales and as a 2.5-kilobase mRNA by Northern analysis of keratinocytes. Identification of this enzyme extends the known distribution of R-lipoxygenases to humans and presents an additional target for potential therapeutic interventions in psoriasis.

Identification and characterization of a retinoid-induced class II tumor suppressor/growth regulatory gene

Retinoids, synthetic and natural analogs of retinoic acid, exhibit potent growth inhibitory and cell differentiation activities that account for their beneficial effects in treating hyperproliferative diseases such as psoriasis, actinic keratosis, and certain neoplasias. Tazarotene is a synthetic retinoid that is used in the clinic for the treatment of psoriasis. To better understand the mechanism of retinoid action in the treatment of hyperproliferative diseases, we used a long-range differential display–PCR to isolate retinoid-responsive genes from primary human keratinocytes. We have identified a cDNA, tazarotene-induced gene 3 (TIG3; Retinoic Acid Receptor Responder 3) showing significant homology to the class II tumor suppressor gene, H-rev 107. Tazarotene treatment increases TIG3 expression in primary human keratinocytes and in vivo in psoriatic lesions. Increased TIG3 expression is correlated with decreased proliferation. TIG3 is expressed in a number of tissues, and expression is reduced in cancer cell lines and some primary tumors. In breast cancer cell lines, retinoid-dependent TIG3 induction is observed in lines that are growth suppressed by retinoids but not in nonresponsive lines. Transient over-expression of TIG3 in T47D or Chinese hamster ovary cells inhibits colony expansion. Finally, studies in 293 cells expressing TIG3 linked to an inducible promoter demonstrated decreased proliferation with increased TIG3 levels. These studies suggest that TIG3 may be a growth regulator that mediates some of the growth suppressive effects of retinoids.

Signature p53 mutation at DNA cross-linking sites in 8-methoxypsoralen and ultraviolet A (PUVA)-induced murine skin cancers.

A combination of psoralen and ultraviolet A radiation (PUVA) is widely used in the treatment of psoriasis. However, PUVA treatment increases the risk of developing skin cancer in psoriasis patients and induces skin cancer in mice. Since the DNA damage induced by PUVA is quite different from that induced by UV, we investigated whether PUVA-induced mouse skin cancers display carcinogen-specific mutations in the p53 tumor suppressor gene. The results indicated that 10 of 13 (77%) PUVA-induced skin tumors contained missense mutations predominantly at exons 6 and 7. In contrast, tumor-adjacent, PUVA-exposed skin from tumor-bearing animals did not exhibit p53 mutation in exons 4-8. Interestingly, about 40% of all mutations in PUVA-induced skin tumors occurred at 5'-TA sites, and an equal number of mutations occurred at one base flanking 5'TA or 5'-TAT sites. Since PUVA induces DNA cross-links exclusively at these sites and since UV "signature" mutations were rarely detected in PUVA-induced skin cancers, we can conclude that PUVA acts as a carcinogen by inducing unique PUVA signature mutations in p53. This finding may have implications for identifying the etiology of skin cancer in psoriasis patients who have undergone PUVA therapy.

A polygenic mouse model of psoriasiform skin disease in CD18-deficient mice.

Previously, a hypomorphic mutation in CD18 was generated by gene targeting, with homozygous mice displaying increased circulating neutrophil counts, defects in the response to chemically induced peritonitis, and delays in transplantation rejection. When this mutation was backcrossed onto the PL/J inbred strain, virtually all homozygous mice developed a chronic inflammatory skin disease with a mean age of onset of 11 weeks after birth. The disease was characterized by erythema, hair loss, and the development of scales and crusts. The histopathology revealed hyperplasia of the epidermis, subcorneal microabscesses, orthohyperkeratosis, parakeratosis, and lymphocyte exocytosis, which are features in common with human psoriasis and other hyperproliferative inflammatory skin disorders. Repetitive cultures failed to demonstrate bacterial or fungal organisms potentially involved in the pathogenesis of this disease, and the dermatitis resolved rapidly after subcutaneous administration of dexamethasone. Homozygous mutant mice on a (PL/J x C57BL/6J)F1 background did not develop the disease and backcross experiments suggest that a small number of genes (perhaps as few as one), in addition to CD18, determine susceptibility to the disorder. This phenotype provides a model for inflammatory skin disorders, may have general relevance to polygenic human inflammatory diseases, and should help to identify genes that interact with the beta2 integrins in inflammatory processes.

Biologic therapies and non-melanoma skin cancer rates in three chronic immune-mediated diseases

Trabalho Final do Curso de Mestrado Integrado em Medicina, Faculdade de Medicina, Universidade de Lisboa, 2014

Transdermal delivery of a tetrapeptide: Evaluation of passive diffusion

Skin penetration of the tetrapeptide Ac-Ala-Ala-Pro-Val-NH2 was assessed. This peptide sequence fits the P-P-1 subsites of elastase and inhibits human neutrophil elastase competitively. Consequently this peptide may be therapeutically useful in a variety of inflammatory disorders, including psoriasis. in which elevated levels of human neutrophil elastase have been reported. Peptide penetration was assessed across whole human skin, whole skin with the stratum corneum removed by tape stripping and epidermis, which had been removed from the dermis by heat separation. The influence of 75% aqueous ethanol as a potential penetration enhancer of the tetrapeptide across epidermis was also assessed. The tetrapeptide did not penetrate whole human skin or epidermis, even under the influence of 75% aqueous ethanol. However, when the stratum corneum was removed tetrapeptide flux of 73.39 mug cm(-2) h(-1) was achieved. The study demonstrates that the stratum corneum is the main barrier to tetrapeptide skin penetration and must be overcome if therapeutically relevant amounts of tetrapeptide are to be delivered to the skin.

Psoriatic arthritis: epidemiology, clinical features, course, and outcome

Psoriatic arthritis (PsA) has been defined as a unique inflammatory arthritis associated with psoriasis. Its exact prevalence is unknown, but estimates vary from 0.3% to 1% of the population. The clinical features described initially are recognised by most experienced clinicians, although they are most distinct in early disease. Initially, PsA typically presents as an oligoarticular and mild disease. However, with time PsA becomes polyarticular, and it is a severe disease in at least 20% of patients. Patients with PsA who present with polyarticular disease are at risk for disease progression. In addition to progression of clinical and radiological damage, health related quality of life is reduced among patients with PsA. It important to note that patients included in recent drug trials resemble patients followed prospectively in a clinic.

Epidermal penetration of a therapeutic peptide by lipid conjugation; Stereo-selective peptide availability of a topical diastereomeric lipopeptide

In inflammatory disorders (e.g. psoriasis), local concentrations of human neutrophil elastase (HNE), also known as polymorphonuclear leukocyte elastase (HLE), possibly overwhelm its natural inhibitors leading to extracellular matrix degradation. Elevated levels of HNE have been reported in a variety of inflammatory disorders, including psoriasis. Peptidic HNE inhibitors have a common hydrophobic sequence (Ala-Ala-Pro-Val). This peptide sequence inhibits HNE competitively; however the stratum corneum presents an effective barrier to the delivery of this tetrapeptide across the skin. The current work investigates the delivery of the modified peptide whereby the tetrapeptide was lipidated to enhance its ability to penetrate the stratum The tetrapeptide Was Coupled to a racaemic mixture of a short chain lipoamino acid (LAA) resulting in two diastereomers of the lipoamino acid-modified tetrapeptide. The penetration of the lipopeptide mixture was assessed across human epidermis in vitro. The percentage of applied dose penetrating to the receptor over 8 h following administration was 2.53% for the D-LAA conjugate and 1.47% for the L-diastereomer, compared to 0% for the peptide. The D-diastereomer appears to be relatively stable but the L-diastereomer appears to degrade releasing possibly the tetrapeptide and peptide fragment(s). Therefore the results clearly indicate that coupling the tetrapeptide to a short chain LAA enhances its delivery across human epidermis.

Complement factor 5a as a therapeutic target

The complement system is an innate immune defense mechanism that protects the host from infection and injury. Complement activation results in the formation of anaphylatoxins, including the biologically active protein C5a. This anaphylatoxin is a potent chemotactic agent for immune and inflammatory cells and induces cell activation. In situations of excessive or uncontrolled complement activation, the overproduction of C5a can cause deleterious effects to the host, and this process is implicated in the pathogenesis of numerous immunoinflammatory disease states, including rheumatoid arthritis, psoriasis, inflammatory bowel disease, ischemia-reperfusion injuries and others. The presence of C5a in a wide variety of condition's has prompted many groups to examine the potential of inhibiting this complement activation product, with the aim of controlling these diseases and reducing the pathologic process. However, to date there is no clinically available specific C5a inhibitor and development of this new drug class is still in a relatively early stage, although limited phase I and phase II human clinical trials have been undertaken in the last few years with selected agents. In this review, examination of the current evidence supporting a specific role of C5a in selected disease states and an overview of potential therapeutic C5a inhibitors will enable the critical evaluation of the potential for C5a as a therapeutic target.