Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Trypanosoma cruzi strains isolated from human, vector, and animal reservoir in the same endemic region in Mexico and typed as T. cruzi I, discrete typing unit 1 exhibit considerable biological diversity

In this study, three strains of Trypanosoma cruzi were isolated at the same time and in the same endemic region in Mexico from a human patient with chronic chagasic cardiomyopathy (RyC-H); vector (Triatoma barberi) (RyC-V); and rodent reservoir (Peromyscus peromyscus) (RyC-R). The three strains were characterized by multilocus enzyme electrophoresis, random amplified polymorphic DNA, and by pathological profiles in experimental animals (biodemes). Based on the analysis of genetic markers the three parasite strains were typed as belonging to T. cruzi I major group, discrete typing unit 1. The pathological profile of RyC-H and RyC-V strains indicated medium virulence and low mortality and, accordingly, the strains should be considered as belonging to biodeme Type III. On the other hand, the parasites from RyC-R strain induced more severe inflammatory processes and high mortality (> 40%) and were considered as belonging to biodeme Type II. The relationship between genotypes and biological characteristics in T. cruzi strains is still debated and not clearly understood. An expert committee recommended in 1999 that Biodeme Type III would correspond to T. cruzi I group, whereas Biodeme Type II, to T. cruzi II group. Our findings suggest that, at least for Mexican isolates, this correlation does not stand and that biological characteristics such as pathogenicity and virulence could be determined by factors different from those identified in the genotypic characterization

Schizodeme analysis of Trypanosoma cruzi Colombian strain clones isolated from the acute phase of murine infection

Colombian strain of Trypanosoma cruzi, biodeme Type III (T. cruzi I), has been cloned by micromanipulation at two phases of the acute infection: early (10 days ) and advanced (30 days). Twelve clones were obtained therefrom. Characterization by their biological and biochemical behavior showed an identity among the several clones and their parental strain, albeit with different degrees of virulence. Molecular characterization of the kinetoplast DNA (kDNA) after amplification by polymerase chain reaction revealed identical profiles of the bands from the kDNA minicircle by the analysis of restriction fragment lenght polymorphism for the isolated clones, their parental strain, and to the clones isolated at two different phases of the infection. Results suggest the predominance of a "principal clone", in the composition of the Colombian strain, responsible for the biological and biochemical behavior. However, no relationship was detected between the molecular profile of kDNA and the degree of virulence presented by the several clones.

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Biological aspects of crosses between Triatoma maculata (Erichson, 1848) and Triatoma pseudomaculata Corrêa & Espínola, 1964 (Hemiptera: Reduviidae)

Due to morphological similarities between Triatoma maculata and T. pseudomaculata, which comprise the "maculate complex", both had been regarded as the same species until 1964. Considering that the studies on triatomine hybridization permit hypotheses formulation concerning origin and divergence of species, enabling a quantitative analysis of taxonomic relationships between species, the present investigation was aimed at broadening further understanding related to the capacity of hybrid production by determining the degree of reproductive isolation between T. maculata and T. pseudomaculata. Our results have demonstrated that T. maculata and T. pseudomaculata showed no differences regarding reproduction patterns and they are able to cross, generating infertile hybrids.

Functional and transcriptional induction of aquaporin-1 gene by hypoxia; analysis of promoter and role of Hif-1α.

Aquaporin-1 (AQP1) is a water channel that is highly expressed in tissues with rapid O(2) transport. It has been reported that this protein contributes to gas permeation (CO(2), NO and O(2)) through the plasma membrane. We show that hypoxia increases Aqp1 mRNA and protein levels in tissues, namely mouse brain and lung, and in cultured cells, the 9L glioma cell line. Stopped-flow light-scattering experiments confirmed an increase in the water permeability of 9L cells exposed to hypoxia, supporting the view that hypoxic Aqp1 up-regulation has a functional role. To investigate the molecular mechanisms underlying this regulatory process, transcriptional regulation was studied by transient transfections of mouse endothelial cells with a 1297 bp 5' proximal Aqp1 promoter-luciferase construct. Incubation in hypoxia produced a dose- and time-dependent induction of luciferase activity that was also obtained after treatments with hypoxia mimetics (DMOG and CoCl(2)) and by overexpressing stabilized mutated forms of HIF-1α. Single mutations or full deletions of the three putative HIF binding domains present in the Aqp1 promoter partially reduced its responsiveness to hypoxia, and transfection with Hif-1α siRNA decreased the in vitro hypoxia induction of Aqp1 mRNA and protein levels. Our results indicate that HIF-1α participates in the hypoxic induction of AQP1. However, we also demonstrate that the activation of Aqp1 promoter by hypoxia is complex and multifactorial and suggest that besides HIF-1α other transcription factors might contribute to this regulatory process. These data provide a conceptual framework to support future research on the involvement of AQP1 in a range of pathophysiological conditions, including edema, tumor growth, and respiratory diseases.

Clinical pharmacogenomic testing of KRAS, BRAF and EGFR mutations by high resolution melting analysis and ultra-deep pyrosequencing

BACKGROUND: Epidermal growth factor receptor (EGFR) and its downstream factors KRAS and BRAF are mutated in several types of cancer, affecting the clinical response to EGFR inhibitors. Mutations in the EGFR kinase domain predict sensitivity to the tyrosine kinase inhibitors gefitinib and erlotinib in lung adenocarcinoma, while activating point mutations in KRAS and BRAF confer resistance to the anti-EGFR monoclonal antibody cetuximab in colorectal cancer. The development of new generation methods for systematic mutation screening of these genes will allow more appropriate therapeutic choices. METHODS: We describe a high resolution melting (HRM) assay for mutation detection in EGFR exons 19-21, KRAS codon 12/13 and BRAF V600 using formalin-fixed paraffin-embedded samples. Somatic variation of KRAS exon 2 was also analysed by massively parallel pyrosequencing of amplicons with the GS Junior 454 platform. RESULTS: We tested 120 routine diagnostic specimens from patients with colorectal or lung cancer. Mutations in KRAS, BRAF and EGFR were observed in 41.9%, 13.0% and 11.1% of the overall samples, respectively, being mutually exclusive. For KRAS, six types of substitutions were detected (17 G12D, 9 G13D, 7 G12C, 2 G12A, 2 G12V, 2 G12S), while V600E accounted for all the BRAF activating mutations. Regarding EGFR, two cases showed exon 19 deletions (delE746-A750 and delE746-T751insA) and another two substitutions in exon 21 (one showed L858R with the resistance mutation T590M in exon 20, and the other had P848L mutation). Consistent with earlier reports, our results show that KRAS and BRAF mutation frequencies in colorectal cancer were 44.3% and 13.0%, respectively, while EGFR mutations were detected in 11.1% of the lung cancer specimens. Ultra-deep amplicon pyrosequencing successfully validated the HRM results and allowed detection and quantitation of KRAS somatic mutations. CONCLUSIONS: HRM is a rapid and sensitive method for moderate-throughput cost-effective screening of oncogene mutations in clinical samples. Rather than Sanger sequence validation, next-generation sequencing technology results in more accurate quantitative results in somatic variation and can be achieved at a higher throughput scale.

Sensitive headspace gas chromatography analysis of free and conjugated 1-methoxy-2-propanol in urine

Glycol ethers still continue to be a workplace hazard due to their important use on an industrial scale. Currently, chronic occupational exposures to low levels of xenobiotics become increasingly relevant. Thus, sensitive analytical methods for detecting biomarkers of exposure are of interest in the field of occupational exposure assessment. 1-Methoxy-2-propanol (1M2P) is one of the dominant glycol ethers and the unmetabolized urinary fraction has been identified to be a good biological indicator of exposure. An existing analytical method including a solid-phase extraction and derivatization before GC/FID analysis is available but presents some disadvantages. We present here an alternative method for the determination of urinary 1M2P based on the headspace gas chromatography technique. We determined the 1M2P values by the direct headspace method for 47 samples that had previously been assayed by the solid-phase extraction and derivatization gas chromatography procedure. An inter-method comparison based on a Bland-Altman analysis showed that both techniques can be used interchangeably. The alternative method showed a tenfold lower limit of detection (0.1 mg/L) as well as good accuracy and precision which were determined by several urinary 1M2P analyses carried out on a series of urine samples obtained from a human volunteer study. The within- and between-run precisions were generally about 10%, which corresponds to the usual injection variability. We observed that the differences between the results obtained with both methods are not clinically relevant in comparison to the current biological exposure index of urinary 1M2P. Accordingly, the headspace gas chromatography technique turned out to be a more sensitive, accurate, and simple method for the determination of urinary 1M2P.[Authors]

Robust regression in biological assay: application to the evaluation of alternative experimental techniques.

Robust Huber type regression and testing of linear hypotheses are adapted to statistical analysis of parallel line and slope ratio assays. They are applied in the evaluation of results of several experiments carried out in order to compare and validate alternatives to animal experimentation based on embryo and cell cultures. Computational procedures necessary for the application of robust methods of analysis used the conversational statistical package ROBSYS. Special commands for the analysis of parallel line and slope ratio assays have been added to ROBSYS.

Demethylation analysis of the FOXP3 locus shows quantitative defects of regulatory T cells in IPEX-like syndrome.

Immune dysregulation, Polyendocrinopathy, Enteropathy X-linked (IPEX) syndrome is a unique example of primary immunodeficiency characterized by autoimmune manifestations due to defective regulatory T (Treg) cells, in the presence of FOXP3 mutations. However, autoimmune symptoms phenotypically resembling IPEX often occur in the absence of detectable FOXP3 mutations. The cause of this "IPEX-like" syndrome presently remains unclear. To investigate whether a defect in Treg cells sustains the immunological dysregulation in IPEX-like patients, we measured the amount of peripheral Treg cells within the CD3(+) T cells by analysing demethylation of the Treg cell-Specific-Demethylated-Region (TSDR) in the FOXP3 locus and demethylation of the T cell-Specific-Demethylated-Region (TLSDR) in the CD3 locus, highly specific markers for stable Treg cells and overall T cells, respectively. TSDR demethylation analysis, alone or normalized for the total T cells, showed that the amount of peripheral Treg cells in a cohort of IPEX-like patients was significantly reduced, as compared to both healthy subjects and unrelated disease controls. This reduction could not be displayed by flow cytometric analysis, showing highly variable percentages of FOXP3(+) and CD25(+)FOXP3(+) T cells. These data provide evidence that a quantitative defect of Treg cells could be considered a common biological hallmark of IPEX-like syndrome. Since Treg cell suppressive function was not impaired, we propose that this reduction per se could sustain autoimmunity.

The biological in vitro effect and selectivity of aromatic dicationic compounds on Trypanosoma cruzi

Trypanosoma cruzi is a parasite that causes Chagas disease, which affects millions of individuals in endemic areas of Latin America. One hundred years after the discovery of Chagas disease, it is still considered a neglected illness because the available drugs are unsatisfactory. Aromatic compounds represent an important class of DNA minor groove-binding ligands that exhibit potent antimicrobial activity. This study focused on the in vitro activity of 10 aromatic dicationic compounds against bloodstream trypomastigotes and intracellular forms of T. cruzi. Our data demonstrated that these compounds display trypanocidal effects against both forms of the parasite and that seven out of the 10 compounds presented higher anti-parasitic activity against intracellular parasites compared with the bloodstream forms. Additional assays to determine the potential toxicity to mammalian cells showed that the majority of the dicationic compounds did not considerably decrease cellular viability. Fluorescent microscopy analysis demonstrated that although all compounds were localised to a greater extent within the kinetoplast than the nucleus, no correlation could be found between compound activity and kDNA accumulation. The present results stimulate further investigations of this class of compounds for the rational design of new chemotherapeutic agents for Chagas disease.

Correlation of biological serum markers with the degree of hepatic fibrosis and necroinflammatory activity in hepatitis C and schistosomiasis patients

Liver biopsy is the gold-standard method to stage fibrosis; however, it is an invasive procedure and is potentially dangerous. The main objective of this study was to evaluate biological markers, such as cytokines IL-13, IFN-γ, TNF-α and TGF-β, platelets, bilirubins (Bil), alanine aminotransferase (ALT) and aspartate aminotransferase (AST), total proteins, γ-glutamil transferase (γ-GT) and alkaline phosphatase (AP), that could be used to predict the severity of hepatic fibrosis in schistosomiasis and hepatitis C (HC) as isolated diseases or co-infections. The following patient groups were selected: HC (n = 39), HC/hepatosplenic schistosomiasis (HSS) (n = 19), HSS (n = 22) and a control group (n = 13). ANOVA and ROC curves were used for statistical analysis. P < 0.05 was considered significant. With HC patients we showed that TNF-α (p = 0.020) and AP (p = 0.005) could differentiate mild and severe fibrosis. With regard to necroinflammatory activity, AST (p = 0.002), γ-GT (p = 0.034) and AP (p = 0.001) were the best markers to differentiate mild and severe activity. In HC + HSS patients, total Bil (p = 0.008) was capable of differentiating between mild and severe fibrosis. In conclusion, our study was able to suggest biological markers that are non-invasive candidates to evaluate fibrosis and necroinflammatory activity in HC and HC + HSS.

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.

Construction and analysis of a new conditional ferritin H knockout mouse strain/

Résumé Le fer joue un rôle important dans la plupart des fonctions biologiques mais sa présence excessive provoque la production de molécules réactives d'oxygène (ROS) qui peuvent contribuer à diverses maladies. La protéine de stockage du fer, la ferritine H, capte l'excès en fer et le stocke sous forme non-toxique, ce qui empêche des dommages potentiels. La délétion de la ferritine H dans des souris knock-out a été essayée antérieurement, mais ces souris mouraient au stade précoce du développement embryonnaire. Pour étudier l'importance du fer, et en particulier son stockage dans la ferritine, et pour pouvoir mieux comprendre les fonctions de la ferritine H, nous avons créé un modèle de souris knock-out conditionnelles de la ferritine H, selon le système classique de Cre-LoxP. Le premier exon et la région du promoteur du gène de la ferritine H ont été entourés de sites loxP. La mortalité embryonnaire provoquée par la délétion constitutive du gène de la ferritine H a été confirmée en croisant nos souris avec des souris exprimant nestin-Cre1. En croisant nos souris avec des souris transgéniques Mx-Cre, nous avons observé que l'induction de Cre par injection de polyI-polyC provoque la délétion presque complète de la ferritine H dans le foie (> 99%) et la rate (> 88%). Ces tissus ont également perdu une grande partie de leur réserve de fer. Cette observation apporte pour la première fois la preuve in vivo que la ferritine H est indispensable pour le stockage du fer, que les fonctions de la ferritine H et de la ferritine L ne sont pas équivalentes, et que la ferritine L ne peut pas assumer seule la fonction de stockage du fer. Dans le foie des souris knock-out, l'expression de l'ARN messager de l'hepcidine a été induite après 10 jours. En même temps, l'expression de l'ARN messager des gènes codant pour des protéines de l'absorption de fer (DMT1, ferroportin, Dcytb1 et hephaestin) a été réprimée mais dans le duodénum seulement. L'expression d'hepcidine est inversément corrélée avec celle des gènes liés à l'absorption de fer. Cette observation corrobore des études antérieures. Mais, en plus, elle montre également que cette répression se produit seulement dans l'intestin. Nous pouvons ainsi tirer la conclusion suivante : ou bien l'hepcidine a un récepteur spécifique dans le duodénum ou bien les gènes liés à l'absorption de fer dans le duodénum ont un facteur spécifique de transcription sensible à l'hepcidine. Aucune répression de DMT1 et de ferroportin n'a été observée dans les macrophages de la rate après l'induction d'hepcidine. La délétion de ferritine H a entraîné une augmentation du taux de mortalité des cellules hépatiques, ainsi que des altérations dans l'architecture normale du tissu de la rate. Vu par l'immunohistologie, le nombre de lymphocytes B et T était réduit dans la rate, tendant à démontrer que la ferritine H et l'homéostase du fer jouent un rôle dans l'immunité. En conclusion, le modèle de souris knock-out conditionnelles de la ferritine H nous fournit un outil précieux pour l'étude in vivo du rôle joué par la ferritine dans l'homéostase du fer, dans les dommages créés par les ROS, ainsi que dans l'apoptose et l'immunité. Summary Iron plays an important role in most biological functions. However, excess of iron results in production of reactive oxygen species (ROS) which could substantially contribute to pathology of various diseases. Ferritin H scavenges excess of iron and stores it in non-toxic form and potentially prevents the damage. Fenitin H targeting in mice has been attempted before, however, straight knockout was lethal in early embryonic stage. To study the role of iron and its storage protein ferritin and to further elucidate ferritin H functions, we aimed at creating a conditional ferritin H knockout mouse model by classical Cre-LoxP system. First exon along with promoter region of the ferritin H gene was foxed. Embryonic lethality of the constitutive ferritin H deletion was confirmed by crossing the foxed mice with mice expressing nestin Cre-1 as transgene. Almost complete deletion was observed in liver (> 99%) and spleen (>88%) upon induction of Cre by injecting polyI-polyC in Fth Lox/Lox; MxCre mice. These tissues also lost substantial fraction of their iron stores. This provides first in vivo evidence that ferritin H is required for iron storage, ferritin H and L functions are not redundant and that ferritin L cannot perform iron storage function alone. Hepcidin mRNA expression was induced after 10 days in the livers of deleted mice and, simultaneously, mRNA expression of iron absorption related genes (DMT 1, ferroportin, Dcytb1 and hephaestin) was repressed in duodenum only. Hepcidin expression is inversely correlated with that of duodenal iron absorption related genes. This is in agreement with previous studies. However, we also show that this repression happens only in intestine. This leads to the conclusion that either hepcidin has a specific receptor in duodenum or the iron absorption related genes have duodenum specific transcription factor that is responsive to hepcidin. No repression of DMT1 and ferroportin was observed in spleen macrophages upon hepcidin induction. Ferritin H deletion showed increased cell death in liver and disruption of normal architecture of spleen. B lymphocytes were reduced in spleen on immunohistology which point towards a role of ferritin H and iron homeostasis in immunity. In conclusion, ferritin H conditional knockout mouse model provides us with an invaluable tool to study the in vivo role of ferritin H in iron homeostasis, ROS mediated damage, apoptosis and immunity.