Functional characterization of E- and P-cadherin in invasive breast cancer cells

BACKGROUND Alterations in the cadherin-catenin adhesion complexes are involved in tumor initiation, progression and metastasis. However, the functional implication of distinct cadherin types in breast cancer biology is still poorly understood. METHODS To compare the functional role of E-cadherin and P-cadherin in invasive breast cancer, we stably transfected these molecules into the MDA-MB-231 cell line, and investigated their effects on motility, invasion and gene expression regulation. RESULTS Expression of either E- and P-cadherin significantly increased cell aggregation and induced a switch from fibroblastic to epithelial morphology. Although expression of these cadherins did not completely reverse the mesenchymal phenotype of MDA-MB-231 cells, both E- and P-cadherin decreased fibroblast-like migration and invasion through extracellular matrix in a similar way. Moreover, microarray gene expression analysis of MDA-MB-231 cells after expression of E- and P-cadherins revealed that these molecules can activate signaling pathways leading to significant changes in gene expression. Although the expression patterns induced by E- and P-cadherin showed more similarities than differences, 40 genes were differentially modified by the expression of either cadherin type. CONCLUSION E- and P-cadherin have similar functional consequences on the phenotype and invasive behavior of MDA-MB-231 cells. Moreover, we demonstrate for the first time that these cadherins can induce both common and specific gene expression programs on invasive breast cancer cells. Importantly, these identified genes are potential targets for future studies on the functional consequences of altered cadherin expression in human breast cancer.

Analysis of CXCL12 3 ' UTR G > A polymorphism in colorectal cancer

Colorectal cancer is one of the most prevalent cancers in developed countries. However, the genetic factors influencing its appearance remain far from being fully characterized. Recently, a G>A functional transition mapping the 3' untranslated region of the CXCL12 gene (rs1801157) has been found to be under-represented among rectal cancer patients when compared to colon cancer patients from a Swedish series. Here we present the results from an independent analysis of CXCL12 rs1801157 in a larger CRC series of Spanish origin in order to analyse the robustness of this association within a different European population. No significant difference was observed between controls and colon or rectal cancer patients. We were also unable to find a correlation between rs1801157 and different prognostic markers such as metastasis development or disease-free survival time. The epidemiologic data involving CXCL12 rs1801157 in colorectal cancer risk are discussed.

Molecular evidence for a functional ecdysone signaling system in Brugia malayi.

BACKGROUND: Filarial nematodes, including Brugia malayi, the causative agent of lymphatic filariasis, undergo molting in both arthropod and mammalian hosts to complete their life cycles. An understanding of how these parasites cross developmental checkpoints may reveal potential targets for intervention. Pharmacological evidence suggests that ecdysteroids play a role in parasitic nematode molting and fertility although their specific function remains unknown. In insects, ecdysone triggers molting through the activation of the ecdysone receptor: a heterodimer of EcR (ecdysone receptor) and USP (Ultraspiracle). METHODS AND FINDINGS: We report the cloning and characterization of a B. malayi EcR homologue (Bma-EcR). Bma-EcR dimerizes with insect and nematode USP/RXRs and binds to DNA encoding a canonical ecdysone response element (EcRE). In support of the existence of an active ecdysone receptor in Brugia we also cloned a Brugia rxr (retinoid X receptor) homolog (Bma-RXR) and demonstrate that Bma-EcR and Bma-RXR interact to form an active heterodimer using a mammalian two-hybrid activation assay. The Bma-EcR ligand-binding domain (LBD) exhibits ligand-dependent transactivation via a GAL4 fusion protein combined with a chimeric RXR in mammalian cells treated with Ponasterone-A or a synthetic ecdysone agonist. Furthermore, we demonstrate specific up-regulation of reporter gene activity in transgenic B. malayi embryos transfected with a luciferase construct controlled by an EcRE engineered in a B. malayi promoter, in the presence of 20-hydroxy-ecdysone. CONCLUSIONS: Our study identifies and characterizes the two components (Bma-EcR and Bma-RXR) necessary for constituting a functional ecdysteroid receptor in B. malayi. Importantly, the ligand binding domain of BmaEcR is shown to be capable of responding to ecdysteroid ligands, and conversely, ecdysteroids can activate transcription of genes downstream of an EcRE in live B. malayi embryos. These results together confirm that an ecdysone signaling system operates in B. malayi and strongly suggest that Bma-EcR plays a central role in it. Furthermore, our study proposes that existing compounds targeting the insect ecdysone signaling pathway should be considered as potential pharmacological agents against filarial parasites.

Discovery and fine mapping of serum protein loci through transethnic meta-analysis.

Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p < 5 × 10(-8)) for serum albumin (HPN-SCN1B, GCKR-FNDC4, SERPINF2-WDR81, TNFRSF11A-ZCCHC2, FRMD5-WDR76, and RPS11-FCGRT, in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein (TNFRS13B, 6q21.3, and ELL2, in up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Anatomic and functional outcome after 23-gauge vitrectomy, peeling, and intravitreal triamcinolone for idiopathic macular epiretinal membrane.

PURPOSE: To report both the functional and anatomic outcome and safety profile of 23-gauge pars plana vitrectomy combined with membrane peeling and intravitreal injection of triamcinolone acetonide in eyes with idiopathic macular epiretinal membranes. METHODS: Retrospective study of 39 consecutive patients who underwent 23-gauge transconjunctival sutureless vitrectomy, membrane peeling, and intravitreal triamcinolone acetonide injection for an idiopathic macular epiretinal membrane between February 2007 and February 2008. Minimum follow-up was 6 months. RESULTS: Thirty-nine eyes of 39 patients were included in the study. The mean follow-up was 7 +/- 2.2 months (range, 6-15 months). Twenty-two eyes (56%) were pseudophakic and 17 (44%) were phakic at the time of surgery. Five of the phakic eyes (29.4%) had worsening of cataracts during the follow-up period. Mean preoperative intraocular pressure was 14 +/- 3.5 mmHg. At the final follow-up, mean intraocular pressure was 14.5 +/- 2.7 mmHg, which did not differ significantly from the intraocular pressure at baseline (P = 0.14, two-tailed t-test). Five (13%) patients needed topical antiglaucoma treatment. Mean preoperative best-corrected visual acuity (BCVA) was 0.28 decimal equivalent (20/71 Snellen equivalent; logarithm of the minimum angle of resolution 0.54 +/- 0.2, range: 1.0-0.2) and improved significantly (P < 0.0001, two-tailed t-test) to a mean of 0.6 decimal equivalent (20/33 Snellen equivalent; logarithm of the minimum angle of resolution 0.22 +/- 0.16, range: 0.6-0) at the final follow-up. The BCVA improved by a mean of 3.2 +/- 2.1 lines (range: 0-8). Twenty-nine patients (74%) demonstrated a gain of > or =3 lines. Mean central macular thickness was 456 +/- 77 microm (mean +/- SD) at baseline, which was significantly reduced at the final follow-up to 327 +/- 79 microm (mean +/- SD; P < 0.0001, two-tailed t-test). Average central macular thickness reduction was 131 +/- 77 microm (mean +/- SD; range: 36-380 microm). A subgroup analysis of 15 selected cases, which had central macular thickness and BCVA measurements after the first postoperative week, demonstrated that 84% of the total final reduction in central macular thickness and 84% of the total final improvement in BCVA occurred already during the first postoperative week. CONCLUSION: Twenty-three-gauge sutureless transconjunctival vitrectomy is a safe and effective technique for the treatment of idiopathic macular epiretinal membranes. The concomitant administration of intravitreal triamcinolone acetonide after pars plana vitrectomy may speed up and improve the anatomic and functional outcome.

Structural components in functional data

Analyzing functional data often leads to finding common factors, for which functional principal component analysis proves to be a useful tool to summarize and characterize the random variation in a function space. The representation in terms of eigenfunctions is optimal in the sense of L-2 approximation. However, the eigenfunctions are not always directed towards an interesting and interpretable direction in the context of functional data and thus could obscure the underlying structure. To overcome such difficulty, an alternative to functional principal component analysis is proposed that produces directed components which may be more informative and easier to interpret. These structural components are similar to principal components, but are adapted to situations in which the domain of the function may be decomposed into disjoint intervals such that there is effectively independence between intervals and positive correlation within intervals. The approach is demonstrated with synthetic examples as well as real data. Properties for special cases are also studied.

Analysis of Centrosome duplication in "Caenorhabditis elegans" : the role of "sas" genes

Abstract: The centrosome is the major microtubule organizing center (MTOC) of most animal cells. As such, it is essential for a number of processes, including polarized secretion or bipolar spindle assembly. Hence, centrosome number needs to be controlled precisely in coordination with DNA replication. Cells early in the cell cycle contain one centrosome that duplicates during S-phase to give rise to two centrosomes that organize a bipolar spindle during mitosis. A failure in this process is likely to engage the spindle assembly checkpoint and threaten genome stability. Despite its importance for normal and uncontrolled proliferation the mechanisms underlying centrosome duplication are still unclear. The Caenorhabditis elegans embryo is well suited to study the mechanisms of centrosome duplication. It allows for the analysis of cellular processes with high temporal and spatial resolution. Gene identification and inactivation techniques are very powerful and a wide set of mutant and transgenic strains facilitates analysis. My thesis project consisted of characterizing three sas-genes: sas-4, sas-5 and sas-¬6. Embryos lacking these genes fail to form a bipolar spindle, hence their name (spindle assembly). I established that sas-4(RNAi) and sas-6(RNAi) embryos do not form daughter centrioles and thus do not duplicate their centrosomes. Furthermore, I showed that both proteins localize to the cytoplasm and are strikingly enriched at centrioles throughout the cell cycle. By performing fluorescent recovery after photobleaching (FRAP) experiments and differentially labeling centrioles, I established that both proteins are recruited to centrioles once per cell cycle when daughter centrioles form. In contrast, SAS-5, PLK-1 and SPD-2 shuttle permanently between the cytoplasm and centrioles. By showing that SAS-5 and SAS-6 interact in vivo, I established a functional relationship between the proteins. Testing the putative human homologue of SAS-6 (HsSAS-6) and a distant relative of SAS-4 (CPAP), I was able to show that these proteins are required for centrosome duplication in human cells. In addition I found that overexpression of GFP¬HsSAS-6 leads to formation of extra centrosomes. In conclusion, we identified and gained important insights into proteins required for centrosome duplication in C. elegans and in human cells. Thus, our work contributes to further elucidate an important step of cell division in normal and malignant tissues. Eventually, this may allow for the development of novel diagnostic or therapeutic reagents to treat cancer patients. Résumé: Le centrosome est le principal centre organisateur des microtubules dans les cellules animales. De ce fait, il est essentiel pour un certain nombre de processus, comme l'adressage polarisé ou la mise en place d'un fuseau bipolaire. Le nombre de centrosome doit être contrôlé de façon précise et en coordination avec la réplication de l'ADN. Au début du cycle cellulaire, les cellules n'ont qu'un seul centrosome qui se duplique au cours de la phase S pour donner naissance à deux centrosomes qui forment le fuseau bipolaire pendant la mitose. Des défauts dans ce processus déclencheront probablement le "checkpoint" d'assemblage du fuseau et menaceront la stabilité du génome. Malgré leurs importances pour la prolifération normale ou incontrôlée des cellules, les mécanismes gouvernant la duplication des centrosomes restent obscures. L'embryon de Caenorhabditis elegans est bien adapté pour étudier les mécanismes de duplication des centrosomes. Il permet l'analyse des processus cellulaires avec une haute résolution spatiale et temporelle. L'identification des gènes et les techniques d'inactivation sont très puissantes et de larges collections de mutants et de lignées transgéniques facilitent les analyses. Mon projet de thèse a consisté à caractérisé trois gènes: sas-4, sas-5 et sas-6. Les embryons ne possédant pas ces gènes ne forment pas de fuseaux bipolaires, d'où leur nom (spindle assembly). J'ai établi que les embryons sas-4(RNAi) et sas-6(RNAi) ne forment pas de centrioles fils, et donc ne dupliquent pas leur centrosome. De plus, j'ai montré que les deux protéines sont localisées dans le cytoplasme et sont étonnamment enrichies aux centrioles tout le long du cycle cellulaire. En réalisant des expériences de FRAP (fluorscence recovery after photobleaching) et en marquant différentiellement les centrioles, j'ai établi que ces deux protéines sont recrutées une fois par cycle cellulaire aux centrioles, au moment de la duplication. Au contraire, SAS-5, PLK-1 et SPD-2 oscillent en permanence entre le cytoplasme et les centrioles. En montrant que SAS-5 et SAS-6 interagissent in vivo, j'ai établi une relation fonctionnelle entre les deux protéines. En testant les homologues humains putatifs de SAS-6 (HsSAS-6) et de SAS-4 (CPAP), j'ai été capable de montrer que ces protéines étaient aussi requises pour la duplication des centrosomes dans les cellules humaines. De plus, j'ai montré que la surexpression de GFP-HsSAS-6 entrainait la formation de centrosomes surnuméraires. En conclusion, nous avons identifié et progressé dans la compréhension de protéines requises pour la duplication des centrosomes chez C. elegans et dans les cellules humaines. Ainsi, notre travail contribue à mieux élucider une étape importante du la division cellulaire dans les cellules normales et malignes. A terme, ceci devrait aider au développement de nouveaux diagnostics ou de traitements thérapeuthiques pour soigner les malades du cancer.

Analysis of a boundary protein delimiting chromatin domains in yeast

SUMMARY The expression state of a eukaryotic gene depends in part on its location in the chromosome. This position effect results from the organization of eukaryotic genomes into discrete functional domains, defined by local differences in chromatin structure. The expression of genes within each domain appears to be defined and maintained by the concerted action of regulatory elements such as promoters, enhancers, silencers and locus control regions. Individual domains may be bordered by boundary elements that separate regions of permissive and silent chromatin. When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres towards more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one of these transcription factors, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino-acid substitutions that similarly inhibited the boundary and histone-binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains. RESUME L'expression des gènes eucaryotes dépend en partie de leur localisation sur les chromosomes. Cet effet de position résulte de l'organisation des génomes eucaryotes en domaines fonctionnels, définis par des changements locaux au niveau de la structure de la chromatine. Dans chacun de ces domaines, l'expression des gènes est définie et maintenue par l'action concertée de différents éléments régulateurs tels que les promoteurs, les amplificateurs, les silenceurs et les locus control régions. Ces domaines peuvent être entourés par des éléments barrière, séparant les régions de chromatine répressive des régions permissive pour l'expression des gènes. Lorsqu'ils se situent à proximité d'éléments chromosomiques comme les telomères, les gènes peuvent être réprimés de manière épigénétique. Chez la levure, cette répression est établie par la propagation des protéines SIR depuis les télomères vers les régions centromériques. Certains facteurs de transcription peuvent empêcher la répression d'un gène, lorsqu'ils sont placés entre ce gène et le télomère. Nous avons étudié un de ces facteurs, CTF-1, qui a la particularité de lier directement l'histone H3. La délétion de certaines parties de CTF-1 a permis de déterminer que la région responsable de l'activité barrière correspond au domaine d'interaction avec H3. Plusieurs mutations points effectuées dans ce domaine inhibent à la fois l'activité barrière et la capacité de lier H3. Des expériences d'immuno-précipitation de la chromatine indiquent que la protéine barrière CTF-1 prévient efficacement la propagation des protéines SIR et sépare des domaines contenant des histones hypo-acétylées de ceux constitués d'histones hyper-acétylées. Ces résultats suggèrent que CTF-1 interagit directement avec l'histone H3 pour empêcher la propagation de la chromatine répressive, délimitant ainsi des domaines de chromatine permissive et des domaines de chromatine silencieuse.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Functional Characterization of a Novel Frameshift Mutation in the C-terminus of the Nav1.5 Channel Underlying a Brugada Syndrome with Variable Expression in a Spanish Family.

INTRODUCTION We functionally analyzed a frameshift mutation in the SCN5A gene encoding cardiac Na(+) channels (Nav1.5) found in a proband with repeated episodes of ventricular fibrillation who presented bradycardia and paroxysmal atrial fibrillation. Seven relatives also carry the mutation and showed a Brugada syndrome with an incomplete and variable expression. The mutation (p.D1816VfsX7) resulted in a severe truncation (201 residues) of the Nav1.5 C-terminus. METHODS AND RESULTS Wild-type (WT) and mutated Nav1.5 channels together with hNavβ1 were expressed in CHO cells and currents were recorded at room temperature using the whole-cell patch-clamp. Expression of p.D1816VfsX7 alone resulted in a marked reduction (≈90%) in peak Na(+) current density compared with WT channels. Peak current density generated by p.D1816VfsX7+WT was ≈50% of that generated by WT channels. p.D1816VfsX7 positively shifted activation and inactivation curves, leading to a significant reduction of the window current. The mutation accelerated current activation and reactivation kinetics and increased the fraction of channels developing slow inactivation with prolonged depolarizations. However, late INa was not modified by the mutation. p.D1816VfsX7 produced a marked reduction of channel trafficking toward the membrane that was not restored by decreasing incubation temperature during cell culture or by incubation with 300 μM mexiletine and 5 mM 4-phenylbutirate. CONCLUSION Despite a severe truncation of the C-terminus, the resulting mutated channels generate currents, albeit with reduced amplitude and altered biophysical properties, confirming the key role of the C-terminal domain in the expression and function of the cardiac Na(+) channel.

Lack of association between the protein tyrosine phosphatase non-receptor type 22 R263Q and R620W functional genetic variants and endogenous non-anterior uveitis.

OBJECTIVE Endogenous uveitis is a major cause of visual loss mediated by the immune system. The protein tyrosine phosphatase non-receptor type 22 (PTPN22) gene encodes a lymphoid-specific phosphatase that plays a key role in T-cell receptor (TCR) signaling. Two independent functional missense single nucleotide polymorphisms (SNPs) located within the PTPN22 gene (R263Q and R620W) have been associated with different autoimmune disorders. We aimed to analyze for the first time the influence of these PTPN22 genetic variants on endogenous non-anterior uveitis susceptibility. METHODS We performed a case-control study of 217 patients with endogenous non-anterior uveitis and 718 healthy controls from a Spanish population. The PTPN22 polymorphisms (rs33996649 and rs2476601) were genotyped using TaqMan allelic discrimination assays. The allele, genotype, carriers, and allelic combination frequencies were compared between cases and controls with χ(2) analysis or Fisher's exact test. RESULTS Our results showed no influence of the studied SNPs in the global susceptibility analysis (rs33996649: allelic P- value=0.92, odds ratio=0.97, 95% confidence interval=0.54-1.75; rs2476601: allelic P- value=0.86, odds ratio=1.04, 95% confidence interval=0.68-1.59). Similarly, the allelic combination analysis did not provide additional information. CONCLUSIONS Our results suggest that the studied polymorphisms of the PTPN22 gene do not play an important role in the pathophysiology of endogenous non-anterior uveitis.

Structural and functional brain changes in middle-aged type 2 diabetic patients: a cross-sectional study.

BACKGROUND Type 2 diabetes mellitus (T2DM) is an emerging risk factor for cognitive impairment. Whether this impairment is a direct effect of this metabolic disorder on brain function, a consequence of vascular disease, or both, remains unknown. Structural and functional neuroimaging studies in patients with T2DM could help to elucidate this question. OBJECTIVE We designed a cross-sectional study comparing 25 T2DM patients with 25 age- and gender-matched healthy control participants. Clinical information, APOE genotype, lipid and glucose analysis, structural cerebral magnetic resonance imaging including voxel-based morphometry, and F-18 fluorodeoxyglucose positron emission tomography were obtained in all subjects. METHODS Gray matter densities and metabolic differences between groups were analyzed using statistical parametric mapping. In addition to comparing the neuroimaging profiles of both groups, we correlated neuroimaging findings with HbA1c levels, duration of T2DM, and insulin resistance measurement (HOMA-IR) in the diabetic patients group. Results: Patients with T2DM presented reduced gray matter densities and reduced cerebral glucose metabolism in several fronto-temporal brain regions after controlling for various vascular risk factors. Furthermore, within the T2DM group, longer disease duration, and higher HbA1c levels and HOMA-IR were associated with lower gray matter density and reduced cerebral glucose metabolism in fronto-temporal regions. CONCLUSION In agreement with previous reports, our findings indicate that T2DM leads to structural and metabolic abnormalities in fronto-temporal areas. Furthermore, they suggest that these abnormalities are not entirely explained by the role of T2DM as a cardiovascular risk factor.

Two functional variants of IRF5 influence the development of macular edema in patients with non-anterior uveitis.

OBJECTIVE Interferon (IFN) signaling plays a crucial role in autoimmunity. Genetic variation in interferon regulatory factor 5 (IRF5), a major regulator of the type I interferon induction, has been associated with risk of developing several autoimmune diseases. In the current study we aimed to evaluate whether three sets of correlated IRF5 genetic variants, independently associated with SLE and with different functional roles, are involved in uveitis susceptibility and its clinical subphenotypes. METHODS Three IRF5 polymorphisms, rs2004640, rs2070197 and rs10954213, representative of each group, were genotyped using TaqMan® allelic discrimination assays in a total of 263 non-anterior uveitis patients and 724 healthy controls of Spanish origin. RESULTS A clear association between two of the three analyzed genetic variants, rs2004640 and rs10954213, and the absence of macular edema was observed in the case/control analysis (P FDR =5.07E-03, OR=1.48, CI 95%=1.14-1.92 and P FDR =3.37E-03, OR=1.54, CI 95%=1.19-2.01, respectively). Consistently, the subphenotype analysis accordingly with the presence/absence of this clinical condition also reached statistical significance (rs2004640: P=0.037, OR=0.69, CI 95%=0.48-0.98; rs10954213: P=0.030, OR=0.67, CI 95%=0.47-0.96), thus suggesting that both IRF5 genetic variants are specifically associated with the lack of macular edema in uveitis patients. CONCLUSION Our results clearly showed for the first time that two functional genetic variants of IRF5 may play a role in the development of macular edema in non-anterior uveitis patients. Identifying genetic markers for macular edema could lead to the possibility of developing novel treatments or preventive therapies.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.