Comparing correspondence analysis and the log-ratio alternative for representing categorical data

We compare correspondance análisis to the logratio approach based on compositional data. We also compare correspondance análisis and an alternative approach using Hellinger distance, for representing categorical data in a contingency table. We propose a coefficient which globally measures the similarity between these approaches. This coefficient can be decomposed into several components, one component for each principal dimension, indicating the contribution of the dimensions to the difference between the two representations. These three methods of representation can produce quite similar results. One illustrative example is given

Analysis of the cytolytic T lymphocyte response of melanoma patients to the naturally HLA-A*0201-associated tyrosinase peptide 368-376.

The human tyrosinase gene codes for two distinct antigens that are recognized by HLA-A*0201-restricted CTLs. For one of them, tyrosinase peptide 368-376, the sequence identified by mass spectrometry in melanoma cell eluates differs from the gene-encoded sequence as a result of posttranslational modification of amino acid residue 370 (asparagine to aspartic acid). Here, we used fluorescent tetrameric complexes ("tetramers") of HLA-A*0201 and tyrosinase peptide 368-376 (YMDGTMSQV) to characterize the CD8+ T-cell response to this antigen in lymphoid cell populations from HLA-A2 melanoma patients. Taking advantage of the presence of significant numbers of tetramer-positive CD8+ T cells in tumor-infiltrated lymph node cells from a melanoma patient, we derived polyclonal and monoclonal tyrosinase peptide 368-376-specific CTLs by tetramer-guided flow cytometric sorting. These CTLs efficiently and specifically lysed HLA-A*0201- and tyrosinase-positive melanoma cells. As assessed with tyrosinase peptide variants, the fine antigen specificity of the CTLs was quite diverse at the clonal level. Flow cytometric analysis of PBMCs stained with tetramers showed that tyrosinase peptide 368-376-specific CD8+ T cells were hardly detectable in peripheral blood of melanoma patients. However, significant numbers of such cells were detected after short-term stimulation of CD8+ lymphocytes with tyrosinase peptide 368-376 in 6 of 10 HLA-A2 melanoma patients. Taken together, these findings emphasize the significant contribution of the natural tyrosinase peptide 368-376 to the antigenic specificities recognized by the tumor-reactive CTLs that may develop in HLA-A2 melanoma patients.

A meta-analysis of the association of fracture risk and body mass index in women.

Several recent studies suggest that obesity may be a risk factor for fracture. The aim of this study was to investigate the association between body mass index (BMI) and future fracture risk at different skeletal sites. In prospective cohorts from more than 25 countries, baseline data on BMI were available in 398,610 women with an average age of 63 (range, 20-105) years and follow up of 2.2 million person-years during which 30,280 osteoporotic fractures (6457 hip fractures) occurred. Femoral neck BMD was measured in 108,267 of these women. Obesity (BMI ≥ 30 kg/m(2) ) was present in 22%. A majority of osteoporotic fractures (81%) and hip fractures (87%) arose in non-obese women. Compared to a BMI of 25 kg/m(2) , the hazard ratio (HR) for osteoporotic fracture at a BMI of 35 kg/m(2) was 0.87 (95% confidence interval [CI], 0.85-0.90). When adjusted for bone mineral density (BMD), however, the same comparison showed that the HR for osteoporotic fracture was increased (HR, 1.16; 95% CI, 1.09-1.23). Low BMI is a risk factor for hip and all osteoporotic fracture, but is a protective factor for lower leg fracture, whereas high BMI is a risk factor for upper arm (humerus and elbow) fracture. When adjusted for BMD, low BMI remained a risk factor for hip fracture but was protective for osteoporotic fracture, tibia and fibula fracture, distal forearm fracture, and upper arm fracture. When adjusted for BMD, high BMI remained a risk factor for upper arm fracture but was also a risk factor for all osteoporotic fractures. The association between BMI and fracture risk is complex, differs across skeletal sites, and is modified by the interaction between BMI and BMD. At a population level, high BMI remains a protective factor for most sites of fragility fracture. The contribution of increasing population rates of obesity to apparent decreases in fracture rates should be explored. © 2014 American Society for Bone and Mineral Research.

The T-381C SNP in BNP gene may be modestly associated with type 2 diabetes: an updated meta-analysis in 49 279 subjects.

A recent study reported an association between the brain natriuretic peptide (BNP) promoter T-381C polymorphism (rs198389) and protection against type 2 diabetes (T2D). As replication in several studies is mandatory to confirm genetic results, we analyzed the T-381C polymorphism in seven independent case-control cohorts and in 291 T2D-enriched pedigrees totalling 39 557 subjects of European origin. A meta-analysis of the seven case-control studies (n = 39 040) showed a nominal protective effect [odds ratio (OR) = 0.86 (0.79-0.94), P = 0.0006] of the CC genotype on T2D risk, consistent with the previous study. By combining all available data (n = 49 279), we further confirmed a modest contribution of the BNP T-381C polymorphism for protection against T2D [OR = 0.86 (0.80-0.92), P = 1.4 x 10(-5)]. Potential confounders such as gender, age, obesity status or family history were tested in 4335 T2D and 4179 normoglycemic subjects and they had no influence on T2D risk. This study provides further evidence of a modest contribution of the BNP T-381C polymorphism in protection against T2D and illustrates the difficulty of unambiguously proving modest-sized associations even with large sample sizes.

DRB1*03:01 haplotypes: differential contribution to multiple sclerosis risk and specific association with the presence of intrathecal IgM bands.

BACKGROUND Multiple sclerosis (MS) is a multifactorial disease with a genetic basis. The strongest associations with the disease lie in the Human Leukocyte Antigen (HLA) region. However, except for the DRB1*15:01 allele, the main risk factor associated to MS so far, no consistent effect has been described for any other variant. One example is HLA-DRB1*03:01, with a heterogeneous effect across populations and studies. We postulate that those discrepancies could be due to differences in the diverse haplotypes bearing that allele. Thus, we aimed at studying the association of DRB1*03:01 with MS susceptibility considering this allele globally and stratified by haplotypes. We also evaluated the association with the presence of oligoclonal IgM bands against myelin lipids (OCMB) in cerebrospinal fluid. METHODS Genotyping of HLA-B, -DRB1 and -DQA1 was performed in 1068 MS patients and 624 ethnically matched healthy controls. One hundred and thirty-nine MS patients were classified according to the presence (M+, 58 patients)/absence (M-, 81 patients) of OCMB. Comparisons between groups (MS patients vs. controls and M+ vs. M-) were performed with the chi-square test or the Fisher exact test. RESULTS Association of DRB1*03:01 with MS susceptibility was observed but with different haplotypic contribution, being the ancestral haplotype (AH) 18.2 the one causing the highest risk. Comparisons between M+, M- and controls showed that the AH 18.2 was affecting only M+ individuals, conferring a risk similar to that caused by DRB1*15:01. CONCLUSIONS The diverse DRB1*03:01-containing haplotypes contribute with different risk to MS susceptibility. The AH 18.2 causes the highest risk and affects only to individuals showing OCMB.

Single cell analysis reveals similar functional competence of dominant and nondominant CD8 T-cell clonotypes.

Immune protection from infectious diseases and cancer is mediated by individual T cells of different clonal origin. Their functions are tightly regulated but not yet fully characterized. Understanding the contribution of each T cell will improve the prediction of immune protection based on laboratory assessment of T-cell responses. Here we developed techniques for simultaneous molecular and functional assessment of single CD8 T cells directly ex vivo. We studied two groups of patients with melanoma after vaccination with two closely related tumor antigenic peptides. Vaccination induced T cells with strong memory and effector functions, as found in virtually all T cells of the first patient group, and fractions of T cells in the second group. Interestingly, high functionality was not restricted to dominant clonotypes. Rather, dominant and nondominant clonotypes acquired equal functional competence. In parallel, this was also found for EBV- and CMV-specific T cells. Thus, the nondominant clonotypes may contribute similarly to immunity as their dominant counterparts.

An integrated genetic and functional analysis of the role of type II transmembrane serine proteases (TMPRSSs) in hearing loss.

Building on our discovery that mutations in the transmembrane serine protease, TMPRSS3, cause nonsyndromic deafness, we have investigated the contribution of other TMPRSS family members to the auditory function. To identify which of the 16 known TMPRSS genes had a strong likelihood of involvement in hearing function, three types of biological evidence were examined: 1) expression in inner ear tissues; 2) location in a genomic interval that contains a yet unidentified gene for deafness; and 3) evaluation of hearing status of any available Tmprss knockout mouse strains. This analysis demonstrated that, besides TMPRSS3, another TMPRSS gene was essential for hearing and, indeed, mice deficient for Hepsin (Hpn) also known as Tmprss1 exhibited profound hearing loss. In addition, TMPRSS2, TMPRSS5, and CORIN, also named TMPRSS10, showed strong likelihood of involvement based on their inner ear expression and mapping position within deafness loci PKSR7, DFNB24, and DFNB25, respectively. These four TMPRSS genes were then screened for mutations in affected members of the DFNB24 and DFNB25 deafness families, and in a cohort of 362 sporadic deaf cases. This large mutation screen revealed numerous novel sequence variations including three potential pathogenic mutations in the TMPRSS5 gene. The mutant forms of TMPRSS5 showed reduced or absent proteolytic activity. Subsequently, TMPRSS genes with evidence of involvement in deafness were further characterized, and their sites of expression were determined. Tmprss1, 3, and 5 proteins were detected in spiral ganglion neurons. Tmprss3 was also present in the organ of Corti. TMPRSS1 and 3 proteins appeared stably anchored to the endoplasmic reticulum membranes, whereas TMPRSS5 was also detected at the plasma membrane. Collectively, these results provide evidence that TMPRSS1 and TMPRSS3 play and TMPRSS5 may play important and specific roles in hearing.

Contribution to Asterisk Open Source Project

With this final master thesis we are going to contribute to the Asterisk open source project. Asterisk is an open source project that started with the main objective of develop an IP telephony platform, completely based on Software (so not hardware dependent) and under an open license like GPL. This project was started on 1999 by the software engineer Mark Spencer at Digium. The main motivation of that open source project was that the telecommunications sector is lack of open solutions, and most of the available solutions are based on proprietary standards, which are close and not compatible between them. Behind the Asterisk project there is a company, Digum, which is the project leading since the project was originated in its laboratories. This company has some of its employees fully dedicated to contribute to the Asterisk project, and also provide the whole infrastructure required by the open source project. But the business of Digium isn't based on licensing of products due to the open source nature of Asterisk, but it's based on offering services around Asteriskand designing and selling some hardware components to be used with Asterisk. The Asterisk project has grown up a lot since its birth, offering in its latest versions advanced functionalities for managing calls and compatibility with some hardware that previously was exclusive of proprietary solutions. Due to that, Asterisk is becoming a serious alternative to all these proprietaries solutions because it has reached a level of maturity that makes it very stable. In addition, as it is open source, it can be fully customized to a givenrequirement, which could be impossible with the proprietaries solutions. Due to the bigness that is reaching the project, every day there are more companies which develop value added software for telephony platforms, that are seriously evaluating the option of make their software fully compatible withAsterisk platforms. All these factors make Asterisk being a consolidated project but in constant evolution, trying to offer all those functionalities offered by proprietaries solutions. This final master thesis will be divided mainly in two blocks totally complementaries. In the first block we will analyze Asterisk as an open source project and Asterisk as a telephony platform (PBX). As a result of this analysis we will generate a document, written in English because it is Asterisk project's official language, which could be used by future contributors as an starting point on joining Asterisk. On the second block we will proceed with a development contribution to the Asterisk project. We will have several options in the form that we do the contribution, such as solving bugs, developing new functionalities or start an Asterisk satellite project. The type of contribution will depend on the needs of the project on that moment.

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Comparative analysis of cancer genes in the human and chimpanzee genomes

Background: Cancer is a major medical problem in modern societies. However, the incidence of this disease in non-human primates is very low. To study whether genetic differences between human and chimpanzee could contribute to their distinct cancer susceptibility, we have examined in the chimpanzee genome the orthologous genes of a set of 333 human cancer genes. Results: This analysis has revealed that all examined human cancer genes are present in chimpanzee, contain intact open reading frames and show a high degree of conservation between both species. However, detailed analysis of this set of genes has shown some differences in genes of special relevance for human cancer. Thus, the chimpanzee gene encoding p53 contains a Pro residue at codon 72, while this codon is polymorphic in humans and can code for Arg or Pro, generating isoforms with different ability to induce apoptosis or interact with p73. Moreover, sequencing of the BRCA1 gene has shown an 8 Kb deletion in the chimpanzee sequence that prematurely truncates the co-regulated NBR2 gene. Conclusion: These data suggest that small differences in cancer genes, as those found in tumor suppressor genes, might influence the differences in cancer susceptibility between human and chimpanzee. Nevertheless, further analysis will be required to determine the exact contribution of the genetic changes identified in this study to the different cancer incidence in non-human primates.

Androgenetic alopecia: identification of four genetic risk loci and evidence for the contribution of WNT signaling to its etiology.

The pathogenesis of androgenetic alopecia (AGA, male-pattern baldness) is driven by androgens, and genetic predisposition is the major prerequisite. Candidate gene and genome-wide association studies have reported that single-nucleotide polymorphisms (SNPs) at eight different genomic loci are associated with AGA development. However, a significant fraction of the overall heritable risk still awaits identification. Furthermore, the understanding of the pathophysiology of AGA is incomplete, and each newly associated locus may provide novel insights into contributing biological pathways. The aim of this study was to identify unknown AGA risk loci by replicating SNPs at the 12 genomic loci that showed suggestive association (5 × 10(-8)<P<10(-5)) with AGA in a recent meta-analysis. We analyzed a replication set comprising 2,759 cases and 2,661 controls of European descent to confirm the association with AGA at these loci. Combined analysis of the replication and the meta-analysis data identified four genome-wide significant risk loci for AGA on chromosomes 2q35, 3q25.1, 5q33.3, and 12p12.1. The strongest association signal was obtained for rs7349332 (P=3.55 × 10(-15)) on chr2q35, which is located intronically in WNT10A. Expression studies in human hair follicle tissue suggest that WNT10A has a functional role in AGA etiology. Thus, our study provides genetic evidence supporting an involvement of WNT signaling in AGA development.

Morphostructural analysis of an alpine debris flows catchment: implication for debris supply

Rock slope instabilities are implicitly linked to the supply of sediment and debris recharging channels prone to debris flow. Hence, the incorporation of bedrock structure and terrain morphology can be relevant in the analysis of sediment budget and debris flow hazard assessment. Here, the mode of debris production of the Manival catchment (northern French Alps) is documented by the study of its morphostructural aspects extracted from high resolution DEM. Terrain implication in the process of debris supply is evaluated by: a) A systematic classification of the major morphological units based on the slope gradient that enables a spatial analysis of zones of debris production and deposition. b) A detailed structural analysis performed on DEM in order to identify potential unstable slopes. c) An analysis of the gullies orientation that informs in term of structural control of the sources zones. d) Localisation of high density joints sets that document about whether sources of continuous debris production are controlled by the structural setting of the catchment. These DEM-based indicators can be used as proxies for assessing the influences of the current topography and enable to quantify a degree of susceptibility to mass wasting and hillslope erosion activity. This present contribution suggests some directions for characterizing sediment flux dynamic in small alpine catchment.

The contributions of rare objects in correspondence analysis

Correspondence analysis, when used to visualize relationships in a table of counts(for example, abundance data in ecology), has been frequently criticized as being too sensitiveto objects (for example, species) that occur with very low frequency or in very few samples. Inthis statistical report we show that this criticism is generally unfounded. We demonstrate this inseveral data sets by calculating the actual contributions of rare objects to the results ofcorrespondence analysis and canonical correspondence analysis, both to the determination ofthe principal axes and to the chi-square distance. It is a fact that rare objects are oftenpositioned as outliers in correspondence analysis maps, which gives the impression that theyare highly influential, but their low weight offsets their distant positions and reduces their effecton the results. An alternative scaling of the correspondence analysis solution, the contributionbiplot, is proposed as a way of mapping the results in order to avoid the problem of outlying andlow contributing rare objects.

Contribution biplots

In order to interpret the biplot it is necessary to know which points usually variables are the ones that are important contributors to the solution, and this information is available separately as part of the biplot s numerical results. We propose a new scaling of the display, called the contribution biplot, which incorporates this diagnostic directly into the graphical display, showing visually the important contributors and thus facilitating the biplot interpretation and often simplifying the graphical representation considerably. The contribution biplot can be applied to a wide variety of analyses such as correspondence analysis, principal component analysis, log-ratio analysis and the graphical results of a discriminant analysis/MANOVA, in fact to any method based on the singular-value decomposition. In the contribution biplot one set of points, usually the rows of the data matrix, optimally represent the spatial positions of the cases or sample units, according to some distance measure that usually incorporates some form of standardization unless all data are comparable in scale. The other set of points, usually the columns, is represented by vectors that are related to their contributions to the low-dimensional solution. A fringe benefit is that usually only one common scale for row and column points is needed on the principal axes, thus avoiding the problem of enlarging or contracting the scale of one set of points to make the biplot legible. Furthermore, this version of the biplot also solves the problem in correspondence analysis of low-frequency categories that are located on the periphery of the map, giving the false impression that they are important, when they are in fact contributing minimally to the solution.

Contribution à la cytotaxonomie des Soricidés (Mammalia, Insectivora) de l'Afrique occidentale.

Chromosomes analysis of 28 shrews belonging to the genus Crocidura from West Africa allowed the detection of five new karyotypes: C. wimmeri (2n = 50, NF = 84), C. cf. gracilipes (56/86), C. cf. nimbae (46/68), and C. cf. planiceps (44/72). Among the individuals identified morphologically as C. poensis (52/70) a distinct species was recognized on the basis of its chromosome complement, C. nigeriae (50/76). The karyotype of C. occidentalis was confirmed with the study of another subspecies, C. occidentalis manni (50/66). An identical chromosome set characterizes C. odorata giffardi and provides evidence for its relationship with C. occidentalis.