Sequential analysis of global gene expression profiles in immature and in vitro matured bovine oocytes: potential molecular markers of oocyte maturation

Background: Without intensive selection, the majority of bovine oocytes submitted to in vitro embryo production (IVP) fail to develop to the blastocyst stage. This is attributed partly to their maturation status and competences. Using the Affymetrix GeneChip Bovine Genome Array, global mRNA expression analysis of immature (GV) and in vitro matured (IVM) bovine oocytes was carried out to characterize the transcriptome of bovine oocytes and then use a variety of approaches to determine whether the observed transcriptional changes during IVM was real or an artifact of the techniques used during analysis. Results: 8489 transcripts were detected across the two oocyte groups, of which similar to 25.0% (2117 transcripts) were differentially expressed (p < 0.001); corresponding to 589 over-expressed and 1528 under-expressed transcripts in the IVM oocytes compared to their immature counterparts. Over expression of transcripts by IVM oocytes is particularly interesting, therefore, a variety of approaches were employed to determine whether the observed transcriptional changes during IVM were real or an artifact of the techniques used during analysis, including the analysis of transcript abundance in oocytes in vitro matured in the presence of a-amanitin. Subsets of the differentially expressed genes were also validated by quantitative real-time PCR (qPCR) and the gene expression data was classified according to gene ontology and pathway enrichment. Numerous cell cycle linked (CDC2, CDK5, CDK8, HSPA2, MAPK14, TXNL4B), molecular transport (STX5, STX17, SEC22A, SEC22B), and differentiation (NACA) related genes were found to be among the several over-expressed transcripts in GV oocytes compared to the matured counterparts, while ANXA1, PLAU, STC1and LUM were among the over-expressed genes after oocyte maturation. Conclusion: Using sequential experiments, we have shown and confirmed transcriptional changes during oocyte maturation. This dataset provides a unique reference resource for studies concerned with the molecular mechanisms controlling oocyte meiotic maturation in cattle, addresses the existing conflicting issue of transcription during meiotic maturation and contributes to the global goal of improving assisted reproductive technology.

Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia

BACKGROUND: It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS: The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS: The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS: This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Novel markers of the human follicle-associated epithelium identified by genomic profiling and microdissection.

BACKGROUND & AIMS: Regulation of gene expression in the follicle-associated epithelium (FAE) over Peyer's patches is largely unknown. CCL20, a chemokine that recruits immature dendritic cells, is one of the few FAE-specific markers described so far. Lymphotoxin beta (LTalpha1beta2) expressed on the membrane of immune cells triggers CCL20 expression in enterocytes. In this study, we measured expression profiles of LTalpha1beta2-treated intestinal epithelial cells and selected CCL20 -coregulated genes to identify new FAE markers. METHODS: Genomic profiles of T84 and Caco-2 cell lines treated with either LTalpha1beta2, flagellin, or tumor necrosis factor alpha were measured using the Affymetrix GeneChip U133A. Clustering analysis was used to select CCL20 -coregulated genes, and laser dissection microscopy and real-time polymerase chain reaction on human biopsy specimens was used to assess the expression of the selected markers. RESULTS: Applying a 2-way analysis of variance, we identified regulated genes upon the different treatments. A subset of genes involved in inflammation and related to the nuclear factor kappaB pathway was coregulated with CCL20 . Among these genes, the antiapoptotic factor TNFAIP3 was highly expressed in the FAE. CCL23 , which was not coregulated in vitro with CCL20 , was also specifically expressed in the FAE. CONCLUSIONS: We have identified 2 novel human FAE specifically expressed genes. Most of the CCL20 -coregulated genes did not show FAE-specific expression, suggesting that other signaling pathways are critical to modulate FAE-specific gene expression.

Étude transcriptionnelle d'une souche pathogène aviaire de Escherichia coli (APEC) et son mutant Pst (phosphate specific transport)

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Profilage Moléculaire de la Scoliose Idiopathique de l’Adolescent

La scoliose idiopathique de l’adolescent (SIA) est la déformation la plus fréquente en orthopédie pédiatrique. C'est une maladie qui génère des déformations complexes du rachis, du thorax et du bassin affectant plus sévèrement et majoritairement les filles à la puberté. Malgré de nombreuses années de recherches approfondies dans le domaine de la SIA, la cause n'a toujours pas été résolue. OBJECTIFS : Il a été démontré qu’il existe une dysfonction de la signalisation de la mélatonine par les protéines G inhibitrices (Gi) dans les ostéoblastes isolés de patients atteints de SIA. Ceci a conduit à une stratification des patients en trois groupes fonctionnels. Le but de ce projet de maîtrise est d’établir les profils d’expression moléculaire correspondant à chacun des groupes fonctionnels. Les objectifs spécifiques sont d’identifier les gènes responsables de l’initiation de la SIA et du défaut différentiel observé dans la signalisation des récepteurs couplés aux protéines Gi. MÉTHODES : Des ARNs ont été préparés à partir d’ostéoblastes isolés de patients atteints de SIA sélectionnés pour chaque groupe fonctionnel et comparés aux ARNs obtenus d’ostéoblastes de sujets sains. En plus, des cellules sanguines (PBMCs) isolées d’une paire de jumelles monozygotes discordantes pour la scoliose ont été étudiées pour comparer leur profil d’expression génétique. Nous avons utilisé des biopuces à ADN contenant un ensemble de 54,000 sondes permettant l’analyse du niveau d’expression de plus de 47,000 transcrits représentant 30,000 gènes humains bien caractérisés (Affymetrix, GeneChip® Human gene 1.0 ST array). Les gènes retenus ont par la suite été validés par qPCR sur un plus grand nombre de patients afin de tester la spécificité des profils d’expression pour chaque groupe de patients à partir des cellules osseuses dérivées lors de la chirurgie. RÉSULTATS: Le profilage moléculaire proposé permettra d’établir les fondements moléculaires de la SIA selon le test fonctionnel développé par le Dr Moreau et son équipe et d’identifier de nouveaux gènes associés à la SIA. Ce projet a permis de mettre en évidence des gènes possiblement impliqués dans le défaut de signalisation associé aux protéines Gi communs aux trois groupes, ainsi que des gènes spécifiques à certains groupes, pouvant contribuer au développement de certaines co-morbidités et/ou au risque d’aggravation des déformations rachidiennes. L’étude préliminaire des jumelles monozygotes discordantes pour la scoliose a mis en évidence un nombre limité de gènes possiblement associés au défaut de signalisation observé chez la jumelle scoliotique, gènes dont l’expression pourrait être influencée par des modifications d’ordre épigénétique liées à l’environnement. CONCLUSION: Les données obtenues par des approches de transcriptomiques dans le cadre de ce projet soutiennent notre méthode de stratification fonctionnelle des patients SIA et ont conduit à l’identification de nouveaux gènes associés à la SIA. Ces résultats jettent un éclairage nouveau sur la SIA et contribuent à une meilleure compréhension des mécanismes et facteurs impliqués dans l’étiologie de la SIA. À cet égard, nos résultats permettront éventuellement d’identifier des cibles thérapeutiques potentielles et des traitements personnalisés compte tenu des profils moléculaires distincts observés pour chaque groupe de patients.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

A genome wide association study of mathematical ability reveals an association at chromosome 3q29, a locus associated with autism and learning difficulties: a preliminary study

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10-5, 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10-6). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10-4). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.

A pooled genome-wide association study of Asperger Syndrome

Asperger Syndrome (AS) is a neurodevelopmental condition characterized by impairments in social interaction and communication, alongside the presence of unusually repetitive, restricted interests and stereotyped behaviour. Individuals with AS have no delay in cognitive and language development. It is a subset of Autism Spectrum Conditions (ASC), which are highly heritable and has a population prevalence of approximately 1%. Few studies have investigated the genetic basis of AS. To address this gap in the literature, we performed a genome-wide pooled DNA association study to identify candidate loci in 612 individuals (294 cases and 318 controls) of Caucasian ancestry, using the Affymetrix GeneChip Human Mapping version 6.0 array. We identified 11 SNPs that had a p-value below 1x10-5. These SNPs were independently genotyped in the same sample. Three of the SNPs (rs1268055, rs7785891 and rs2782448) were nominally significant, though none remained significant after Bonferroni correction. Two of our top three SNPs (rs7785891 and rs2782448) lie in loci previously implicated in ASC. However, investigation of the three SNPs in the ASC genome-wide association dataset from the Psychiatric Genomics Consortium indicated that these three SNPs were not significantly associated with ASC. The effect sizes of the variants were modest, indicating that our study was not sufficiently powered to identify causal variants with precision.

Search for Genomic Alterations in Monozygotic Twins Discordant for Cleft Lip and/or Palate

Phenotypically discordant monozygotic twins offer the possibility of gene discovery through delineation of molecular abnormalities in one member of the twin pair. One proposed mechanism of discordance is postzygotically occurring genomic alterations resulting from mitotic recombination and other somatic changes. Detection of altered genomic fragments can reveal candidate gene loci that can be verified through additional analyses. We investigated this hypothesis using array comparative genomic hybridization; the 50K and 250K Affymetrix GeneChip (R) SNP arrays and an Illumina custom array consisting of 1,536 SNPs, to scan for genomic alterations in a sample of monozygotic twin pairs with discordant cleft lip and/or palate phenotypes. Paired analysis for deletions, amplifications and loss of heterozygosity, along with sequence verification of SNPs with discordant genotype calls did not reveal any genomic discordance between twin pairs in lymphocyte DNA samples. Our results demonstrate that postzygotic genomic alterations are not a common cause of monozygotic twin discordance for isolated cleft lip and/or palate. However, rare or balanced genomic alterations, tissue-specific events and small aberrations beyond the detection level of our experimental approach cannot be ruled out. The stability of genomes we observed in our study samples also suggests that detection of discordant events in other monozygotic twin pairs would be remarkable and of potential disease significance.

Genetische und funktionelle Analyse von Vesikelmembranprotein-Mutanten in Maus und Caenorhabditis elegans

Tetraspan vesicle membrane proteins (TVPs) sind konservierte, ubiquitär vorkommende Membranproteine synaptischer Vesikel und zytoplasmatischer Transportvesikel. Bei Säugetieren lassen sie sich in die Physine, Gyrine und SCAMPs (secretory carrier-associated membrane proteins) unterteilen, die im Nematoden C. elegans jeweils nur durch ein einzelnes Polypeptid vertreten sind (Synaptophysin-1 [SPH-1], Synaptogyrin-1 [SNG-1] und SCAMP-1 [SCM-1]). Obwohl den TVPs eine Beteiligung bei der Regulation des Vesikelzyklus zugesprochen wurde, sind Synaptophysin-1-Knockout-Mäuse und vollständig TVP-defiziente Würmer gesund und weisen nur geringgradige Veränderungen auf. In dieser Arbeit sollten daher zum einen genomweite komparative Transkriptomanalysen durchgeführt werden, um mögliche Kompensationsmechanismen in der Maus und C. elegans zu finden, zum anderen sollten mit Hilfe pharmakologischer Stressassays und genetischer Verfahren Schwachstellen und Redundanzen identifiziert werden. Erstaunlicherweise konnten durch Affymetrix GeneChip-Analysen der RNA in der Retina von Synaptophysin-1-/--Mäusen keine differenziell exprimierten Gene gefunden werden. Bei der Untersuchung der C. elegans-TVP-Dreifachmutante wurden hingegen 17 Gene mit erhöhter und 3 mit erniedrigter Transkription identifiziert. Die Befunde für 12 hochregulierte Gene wurden durch quantitative Real-Time RT-PCR bestätigt. Das am stärksten hochregulierte Gen arf-1.1 kodiert für eine GTPase, die vermutlich an der Regulation der Vesikelbildung beteiligt ist. Von den ebenso identifizierten Genen cdr-2, cdr-4 und pgp-9 ist bekannt, dass sie in Stresssituationen, z. B. in Gegenwart von Cadmium, verstärkt transkribiert werden. ugt-62 und ugt-19 kodieren für Glucuronosyltransferasen. Für arf-1.1, cdr-2, ugt-62 sowie für das Gen T16G1.6, das für eine coiled-coil-Domäne kodiert, wurden im Folgenden fluoreszierende Promoterkonstrukte hergestellt, um Koexpressionsmuster mit TVPs zu bestimmen. Es stellte sich heraus, dass alle vier Promoterkonstrukte im Darm zusammen mit SPH-1 und SCM-1 im Darm transkribiert werden. Mit fluoreszierenden Translationschimären konnte weiterhin gezeigt werden, dass ARF-1.1 und CDR-2 mit den Darm-spezifischen TVPs im apikalen Bereich der Darmzellen kolokalisieren. Um mehr über die Funktion von TVPs im Vesikelzyklus zu erfahren, wurden pharmakologische und genetische Analysen von Würmern durchgeführt, in denen die Expression des Neuronen-spezifischen SNG-1 verändert ist. Deletion oder Überexpression führte zu einer Resistenz gegenüber dem Acetylcholinesterase-Inhibitor Aldicarb und zu erhöhter Empfindlichkeit gegenüber dem GABA-Rezeptor-Antagonisten Pentylentetrazol. Auf genetischer Ebene zeigte sich, dass sng-1 synthetisch mit den Genen für Synaptotagmin-1, Endophilin A sowie Synaptojanin wirkt. Die beobachteten Effekte weisen auf alternative Funktionen in der synaptischen Übertragung hin und unterstützen zugleich die Hypothese, dass SNG-1 im synaptischen Vesikelzyklus eine wichtige Funktion erfüllt, die möglicherweise einem noch unbekannten redundanten Kompartiment-spezifischen Signalweg der synaptischen Transmission zuzuordnen ist.

Examination of FGFRL1 as a candidate gene for diaphragmatic defects at chromosome 4p16.3 shows that Fgfrl1 null mice have reduced expression of Tpm3, sarcomere genes and Lrtm1 in the diaphragm

Fgfrl1 (also known as Fgfr5; OMIM 605830) homozygous null mice have thin, amuscular diaphragms and die at birth because of diaphragm hypoplasia. FGFRL1 is located at 4p16.3, and this chromosome region can be deleted in patients with congenital diaphragmatic hernia (CDH). We examined FGFRL1 as a candidate gene for the diaphragmatic defects associated with 4p16.3 deletions and re-sequenced this gene in 54 patients with CDH. We confirmed six known coding single nucleotide polymorphisms (SNPs): c.209G > A (p.Pro20Pro), c.977G > A (p.Pro276Pro), c.1040T > C (p.Asp297Asp), c.1234C > A (p.Pro362Gln), c.1420G > T (p.Arg424Leu), and c.1540C > T (p.Pro464Leu), but we did not identify any gene mutations. We genotyped additional CDH patients for four of these six SNPs, including the three non-synonymous SNPs, to make a total of 200 chromosomes, and found that the allele frequency for the four SNPs, did not differ significantly between patients and normal controls (p > or = 0.05). We then used Affymetrix Genechip Mouse Gene 1.0 ST arrays and found eight genes with significantly reduced expression levels in the diaphragms of Fgfrl1 homozygous null mice when compared with wildtype mice-Tpm3, Fgfrl1 (p = 0.004), Myl2, Lrtm1, Myh4, Myl3, Myh7 and Hephl1. Lrtm1 is closely related to Slit3, a protein associated with herniation of the central tendon of the diaphragm in mice. The Slit proteins are known to regulate axon branching and cell migration, and inhibition of Slit3 reduces cell motility and decreases the expression of Rac and Cdc42, two genes that are essential for myoblast fusion. Further studies to determine if Lrtm1 has a similar function to Slit3 and if reduced Fgfrl1 expression can cause diaphragm hypoplasia through a mechanism involving decreased myoblast motility and/or myoblast fusion, seem indicated.

Life stage differences in mammary gland gene expression profile in non-human primates

Breast cancer (BC) is the most common malignancy of women in the developed world. To better understand its pathogenesis, knowledge of normal breast development is crucial, as BC is the result of disregulation of physiologic processes. The aim of this study was to investigate the impact of reproductive life stages on the transcriptional profile of the mammary gland in a primate model. Comparative transcriptomic analyses were carried out using breast tissues from 28 female cynomolgus macaques (Macaca fascicularis) at the following life stages: prepubertal (n = 5), adolescent (n = 4), adult luteal (n = 5), pregnant (n = 6), lactating (n = 3), and postmenopausal (n = 5). Mammary gland RNA was hybridized to Affymetrix GeneChip(®) Rhesus Macaque Genome Arrays. Differential gene expression was analyzed using ANOVA and cluster analysis. Hierarchical cluster analysis revealed distinct separation of life stage groups. More than 2,225 differentially expressed mRNAs were identified. Gene families or pathways that changed across life stages included those related to estrogen and androgen (ESR1, PGR, TFF1, GREB1, AR, 17HSDB2, 17HSDB7, STS, HSD11B1, AKR1C4), prolactin (PRLR, ELF5, STAT5, CSN1S1), insulin-like growth factor signaling (IGF1, IGFBP1, IGFBP5), extracellular matrix (POSTN, TGFB1, COL5A2, COL12A1, FOXC1, LAMC1, PDGFRA, TGFB2), and differentiation (CD24, CD29, CD44, CD61, ALDH1, BRCA1, FOXA1, POSTN, DICER1, LIG4, KLF4, NOTCH2, RIF1, BMPR1A, TGFB2). Pregnancy and lactation displayed distinct patterns of gene expression. ESR1 and IGF1 were significantly higher in the adolescent compared to the adult animals, whereas differentiation pathways were overrepresented in adult animals and pregnancy-associated life stages. Few individual genes were distinctly different in postmenopausal animals. Our data demonstrate characteristic patterns of gene expression during breast development. Several of the pathways activated during pubertal development have been implicated in cancer development and metastasis, supporting the idea that other developmental markers may have application as biomarkers for BC.

Stochastic Models Based on Molecular Hybridization Theory for Short Oligonucleotide Microarrays

High density oligonucleotide expression arrays are a widely used tool for the measurement of gene expression on a large scale. Affymetrix GeneChip arrays appear to dominate this market. These arrays use short oligonucleotides to probe for genes in an RNA sample. Due to optical noise, non-specific hybridization, probe-specific effects, and measurement error, ad-hoc measures of expression, that summarize probe intensities, can lead to imprecise and inaccurate results. Various researchers have demonstrated that expression measures based on simple statistical models can provide great improvements over the ad-hoc procedure offered by Affymetrix. Recently, physical models based on molecular hybridization theory, have been proposed as useful tools for prediction of, for example, non-specific hybridization. These physical models show great potential in terms of improving existing expression measures. In this paper we demonstrate that the system producing the measured intensities is too complex to be fully described with these relatively simple physical models and we propose empirically motivated stochastic models that compliment the above mentioned molecular hybridization theory to provide a comprehensive description of the data. We discuss how the proposed model can be used to obtain improved measures of expression useful for the data analysts.

Experimental heart transplantation: effect of cyclosporine on expression and activity of metzincins

Metzincins, such as matrix metalloproteases (MMP), and extracellular matrix (ECM) proteins are differentially regulated in inflammation. We hypothesised that metzincins are also dysregulated in experimental acute cardiac allograft rejection. We investigated the Dark Agouti-to-Lewis (DA-to-Lew) rat model of acute cardiac allograft rejection. Cyclosporine (CsA) (7.5 mg/kg/d) was given from transplantation to sacrifice (day +5). At that time, mRNA levels were analysed by Affymetrix genechip and quantitative reverse transcription polymerase chain reaction (qRTPCR). MMP protein and activities were analysed by immunohistology, fluorometry, zymography and Western blots. In untreated rejected DA allografts, mRNA levels of MMP-2/-7/-9/-/12-/14, a disintegrin and metalloprotease (ADAM)-17, tissue inhibitor of metalloprotease (TIMP)-1/-3 were increased, whereas MMP-11/-16/-24 and TIMP-2/-4 were lowered compared to native DA hearts. With respect to these untreated allografts, CsA lowered mRNA levels of MMP-7, TIMP-1/-3 (TIMP-2/-4 remained relatively low) and ADAM17, but augmented mRNA levels of MMP-11/-16/-23 and of many ECM genes. Immunohistology showed increased staining of MMP-2 in acute rejection (AR). Overall MMP activity was augmented in both transplanted groups, but CsA reduced MMP-9 activity and MMP-14 production. Taken together, MMP and TIMP were upregulated during acute AR. CsA ameliorated histology of rejection but showed potential pro-fibrotic effects. Thus, MMP and TIMP may play a role in acute cardiac allograft rejection, and beneficial modification of the MMP-ECM balance requires interventions beyond CsA.