Copy number variation modifies expression profiles over developmental time

Abstract : A preliminary understanding of the phenotypic effect of copy number variation (CNV) of DNA segments is emerging. These rearrangements were shown to influence, in a somewhat dose-dependent manner, the expression of genes mapping within them. They were also shown to modify the expression of genes located on their Hanks, sometimes at great distance. Here, we demonstrate by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, 'we find that some brain- expressed genes mapping within CNVS appear to be under compensatory loops only at specific time-points, indicating that the effect of CNVS on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time-course expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but its timing of expression as well. Résume : Nous commençons à comprendre les effets phénotypiques liés aux séquences d'ADN qui changent de nombre de copies d'un individu a l'autre. Des travaux précédents ont montré que ces variante de nombre de copies (CNVS) avaient une influence sur l'expression non seulement des gènes se trouvant dans le réarrangement, mais aussi sur ceux se trouvant à une certaine distance. Le présent travail étudie ces effets à différents stades du développement de la souris allant d'un embryon de deux semaines à la souris adulte. Nous avons observé que certains gènes exprimés dans le cerveau semblent soumis à un contrôle plus strict a certaines étapes du développement suggèrent que l'effet des CNVs est modulé différemment au cours de la vie. Notre travail sur trois souches différentes de souris a permis de montrer que les gènes ayant un profil d'expression différent dans le temps entre souches sont enrichis en gènes se trouvant dans des CNVs. Ceci nous amène à penser que les CNVs ont, non seulement une influence sur le niveau d'expression des gènes, mais aussi sur les moments durant lesquels ils seront exprimés. Résumé pour un large public : De nombreuses maladies sont dues soit a un gain (on parle alors de duplication) soit à une perte de matériel génétique (il s'agit dune délétion). Bien que les recherches visant à identifier les mécanismes moléculaires liés à ces réarrangements de notre génome progressent continuellement, la plupart des causes des maladies génétiques restent à élucider. Certaines parties de notre génome sont présentes en un nombre de copies qui diffère d'un individu à l'autre sans pour autant provoquer une ou des maladies. Ces segments d'ADN qui varient en nombre sont appelés Copy Number Variant (CNVs). Ils couvrent environ 12% de notre matériel génétique. Des études menées sur différents modèles animaux ont montré que les CNVs avaient une influence aussi bien sur les gènes qui sont a l'intérieur des CNVs que sur ceux qui sont dans leur voisinage. Ce travail étudie l'effet des CNVs à travers différents stades du développement de la souris. Nous avons démontré que les segments d'ADN qui varient en nombre de copies ont des effets variables selon le stade auxquels ils sont mesurés. Ainsi, les CNVs ont non seulement un impact sur l'expression des gènes présents dans ces régions et dans leur voisinage, mais influencent également leurs profils d'expression au cours du temps.

Copy number variation modifies expression time courses.

A preliminary understanding into the phenotypic effect of DNA segment copy number variation (CNV) is emerging. These rearrangements were demonstrated to influence, in a somewhat dose-dependent manner, the expression of genes that map within them. They were also shown to modify the expression of genes located on their flanks and sometimes those at a great distance from their boundary. Here we demonstrate, by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, we find that some brain-expressed genes mapping within CNVs appear to be under compensatory loops only at specific time points, indicating that the effect of CNVs on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time courses of expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but also the timing of its expression.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Segmental copy number variation shapes tissue transcriptomes.

Copy number variation (CNV) is a key source of genetic diversity, but a comprehensive understanding of its phenotypic effect is only beginning to emerge. We have generated a CNV map in wild mice and classical inbred strains. Genome-wide expression data from six major organs show not only that expression of genes within CNVs tend to correlate with copy number changes, but also that CNVs influence the expression of genes in their vicinity, an effect that extends up to half a megabase. Genes within CNVs show lower expression and more specific spatial expression patterns than genes mapping elsewhere. Our analyses reveal differential constraint on copy number changes of genes expressed in different tissues. Dosage alterations of brain-expressed genes are less frequent than those of other genes and are buffered by tighter transcriptional regulation. Our study provides initial evidence that CNVs shape tissue transcriptomes on a global scale.

Evolutionary fate of retroposed gene copies in the human genome.

Given that retroposed copies of genes are presumed to lack the regulatory elements required for their expression, retroposition has long been considered a mechanism without functional relevance. However, through an in silico assay for transcriptional activity, we identify here >1,000 transcribed retrocopies in the human genome, of which at least approximately 120 have evolved into bona fide genes. Among these, approximately 50 retrogenes have evolved functions in testes, more than half of which were recruited as functional autosomal counterparts of X-linked genes during spermatogenesis. Generally, retrogenes emerge "out of the testis," because they are often initially transcribed in testis and later evolve stronger and sometimes more diverse spatial expression patterns. We find a significant excess of transcribed retrocopies close to other genes or within introns, suggesting that retrocopies can exploit the regulatory elements and/or open chromatin of neighboring genes to become transcribed. In direct support of this hypothesis, we identify 36 retrocopy-host gene fusions, including primate-specific chimeric genes. Strikingly, 27 intergenic retrogenes have acquired untranslated exons de novo during evolution to achieve high expression levels. Notably, our screen for highly transcribed retrocopies also uncovered a retrogene linked to a human recessive disorder, gelatinous drop-like corneal dystrophy, a form of blindness. These functional implications for retroposition notwithstanding, we find that the insertion of retrocopies into genes is generally deleterious, because it may interfere with the transcription of host genes. Our results demonstrate that natural selection has been fundamental in shaping the retrocopy repertoire of the human genome.

Emergent geomorphic-vegetation interactions on a subalpine alluvial fan

Following perturbation, an ecosystem (flora, fauna, soil) should evolve as a function of time at a rate conditioned by external variables (relief, climate, geology). More recently, biogeomorphologists have focused upon the notion of co-development of geomorphic processes with ecosystems over very short through to very long (evolutionary) timescales. Alpine environments have been a particular focus of this co-development. However, work in this field has tended to adopt a simplified view of the relationship between perturbation and succession, including: how the landform and ecosystem itself conditions the impact of a perturbation to create a complex spatial response impact; and how perturbations are not simply ecosystem destroyers but can be a significant source of ecosystem resources. What this means is that at the within landform scale, there may well be a complex and dynamic topographic and sedimentological template that co-develops with soil, flora and fauna. Here, we present and test a conceptual model of this template for a subalpine alluvial fan. We combine detailed floristic inventory with soil inventory, determination of edaphic variables and analysis of historical aerial imagery. Spatial variation in the probability of perturbation of sites on the fan surface was associated with down fan variability in the across-fan distribution of fan ages, fan surface channel characteristics and fan surface sedimentology. Floristic survey confirmed that these edaphic factors distinguished site floristic richness and plant communities up until the point that the soil-vegetation system was sufficiently developed to sustain plant communities regardless of edaphic conditions. Thus, the primary explanatory variable was the estimated age of each site, which could be tied back into perturbation history and its spatial expression due to the geometry of the fan: distinct plant communities were emergent both across fan and down fan, a distribution maintained by the way in which the fan dissipates potentially perturbing events.

Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding.

A new metabolite profiling approach combined with an ultrarapid sample preparation procedure was used to study the temporal and spatial dynamics of the wound-induced accumulation of jasmonic acid (JA) and its oxygenated derivatives in Arabidopsis thaliana. In addition to well known jasmonates, including hydroxyjasmonates (HOJAs), jasmonoyl-isoleucine (JA-Ile), and its 12-hydroxy derivative (12-HOJA-Ile), a new wound-induced dicarboxyjasmonate, 12-carboxyjasmonoyl-l-isoleucine (12-HOOCJA-Ile) was discovered. HOJAs and 12-HOOCJA-Ile were enriched in the midveins of wounded leaves, strongly differentiating them from the other jasmonate metabolites studied. The polarity of these oxylipins at physiological pH correlated with their appearance in midveins. When the time points of accumulation of different jasmonates were determined, JA levels were found to increase within 2-5 min of wounding. Remarkably, these changes occurred throughout the plant and were not restricted to wounded leaves. The speed of the stimulus leading to JA accumulation in leaves distal to a wound is at least 3 cm/min. The data give new insights into the spatial and temporal accumulation of jasmonates and have implications in the understanding of long-distance wound signaling in plants.

Expression and developmental regulation of Ehk-1, a neuronal Elk-like receptor tyrosine kinase in brain.

Protein tyrosine kinases are pivotal in central nervous tissue development and maintenance. Here we focus on the expression of Ehk-1, a novel Elk-related receptor tyrosine kinase. Ehk-1 gene expression is observed in the developing and adult central nervous system and is highly regulated throughout development at both the messenger RNA and protein levels. Three messenger RNA transcripts of 8.5, 5.9 and 5.1 kb are detectable in the rat brain and a variety of splice possibilities have been identified. However, a major protein species of around M(r) 120,000 predominates throughout development. Ehk-1 messenger RNA and protein levels are highest in the first postnatal week. By in situ messenger RNA hybridization the gene is expressed by all neurons of the adult brain, but mostly in the hippocampus, cerebral cortex and large neurons of the deep cerebellar nuclei, as well as the Purkinje and granular cells of the cerebellum. At earlier stages of development, transcripts are most prominent in the periventricular germinal layers of the brain. Immunohistochemistry reveals a pronounced membrane associated protein expression in immature neurons. In the adult animal, peak reactivity was found in the neuropil with sparing of most perikarya. The spatial and temporal pattern of ehk-1 gene expression suggests a role in both the development and maintenance of differentiated neurons of the central nervous system.

The effect of translocation-induced nuclear reorganization on gene expression.

Translocations are known to affect the expression of genes at the breakpoints and, in the case of unbalanced translocations, alter the gene copy number. However, a comprehensive understanding of the functional impact of this class of variation is lacking. Here, we have studied the effect of balanced chromosomal rearrangements on gene expression by comparing the transcriptomes of cell lines from controls and individuals with the t(11;22)(q23;q11) translocation. The number of differentially expressed transcripts between translocation-carrying and control cohorts is significantly higher than that observed between control samples alone, suggesting that balanced rearrangements have a greater effect on gene expression than normal variation. Many of the affected genes are located along the length of the derived chromosome 11. We show that this chromosome is concomitantly altered in its spatial organization, occupying a more central position in the nucleus than its nonrearranged counterpart. Derivative 22-mapping chromosome 22 genes, on the other hand, remain in their usual environment. Our results are consistent with recent studies that experimentally altered nuclear organization, and indicated that nuclear position plays a functional role in regulating the expression of some genes in mammalian cells. Our study suggests that chromosomal translocations can result in hitherto unforeseen, large-scale changes in gene expression that are the consequence of alterations in normal chromosome territory positioning. This has consequences for the patterns of gene expression change seen during tumorigenesis-associated genome instability and during the karyotype changes that lead to speciation.

SOST expression is restricted to the great arteries during embryonic and neonatal cardiovascular development.

Spatial-temporal regulation of bone morphogenetic protein (BMP) and Wnt activity is essential for normal cardiovascular development, and altered activity of these growth factors causes maldevelopment of the cardiac outflow tract and great arteries. In the present study, we show that SOST, a Dan family member reported to antagonize BMP and Wnt activity, is expressed within the medial vessel wall of the great arteries containing smooth muscle cells. The ascending aorta, aortic arch, brachiocephalic artery, common carotids, and pulmonary trunk were all associated with SOST expressing smooth muscle cells, while the heart itself, including the valves, and more distal arteries, that is, pulmonary arteries, subclavian arteries, and descending aorta, were negative. SOST was expressed from embryonic day 15.5 up to the neonatal period. SOST expression, however, did not correspond with inhibition of Smad-dependent BMP activity or beta-catenin-dependent Wnt activity in the great arteries. Activity of both signaling pathways was already down-regulated before induction of SOST expression.

Spatial, temporal and subcellular localization of islet-brain 1 (IB1), a homologue of JIP-1, in mouse brain.

Islet-brain 1 (IB1) was recently identified as a DNA-binding protein of the GLUT2 gene promoter. The mouse IB1 is the rat and human homologue of the Jun-interacting protein 1 (JIP-1) which has been recognized as a key player in the regulation of c-Jun amino-terminal kinase (JNK) mitogen-activated protein kinase (MAPK) pathways. JIP-1 is involved in the control of apoptosis and may play a role in brain development and aging. Here, IB1 was studied in adult and developing mouse brain tissue by in situ hybridization, Northern and Western blot analysis at cellular and subcellular levels, as well as by immunocytochemistry in brain sections and cell cultures. IB1 expression was localized in the synaptic regions of the olfactory bulb, retina, cerebral and cerebellar cortex and hippocampus in the adult mouse brain. IB1 was also detected in a restricted number of axons, as in the mossy fibres from dentate gyrus in the hippocampus, and was found in soma, dendrites and axons of cerebellar Purkinje cells. After birth, IB1 expression peaks at postnatal day 15. IB1 was located in axonal and dendritic growth cones in primary telencephalon cells. By biochemical and subcellular fractionation of neuronal cells, IB1 was detected both in the cytosolic and membrane fractions. Taken together with previous data, the restricted neuronal expression of IB1 in developing and adult brain and its prominent localization in synapses suggest that the protein may be critical for cell signalling in developing and mature nerve terminals.

Formation of a spatial memory trace : a behavioral, cellular and molecular study

THESIS ABSTRACTThis thesis project was aimed at studying the molecular mechanisms underlying learning and memory formation, in particular as they relate to the metabolic coupling between astrocytes and neurons. For that, changes in the metabolic activity of different mice brain regions after 1 or 9 days of training in an eight-arm radial maze were assessed by (14C) 2-deoxyglucose (2DG) autoradiography. Significant differences in the areas engaged during the behavioral task at day 1 (when animals are confronted for the first time to the learning task) and at day 9 (when animals are highly performing) have been identified. These areas include the hippocampus, the fornix, the parietal cortex, the laterodorsal thalamic nucleus and the mammillary bodies at day 1 ; and the anterior cingulate, the retrosplenial cortex and the dorsal striatum at day 9. Two of these cerebral regions (those presenting the greatest changes at day 1 and day 9: the hippocampus and the retrosplenial cortex, respectively) were microdissected by laser capture microscopy and selected genes related to neuron-glia metabolic coupling, glucose metabolism and synaptic plasticity were analyzed by RT-PCR. 2DG and gene expression analysis were performed at three different times: 1) immediately after the end of the behavioral paradigm, 2) 45 minutes and 3) 6 hours after training. The main goal of this study was the identification of the metabolic adaptations following the learning task. Gene expression results demonstrate that the learning task profoundly modulates the pattern of gene expression in time, meaning that these two cerebral regions with high 2DG signal (hippocampus and retrosplenial cortex) have adapted their metabolic molecular machinery in consequence. Almost all studied genes show a higher expression in the hippocampus at day 1 compared to day 9, while an increased expression was found in the retrosplenial cortex at day 9. We can observe these molecular adaptations with a short delay of 45 minutes after the end of the task. However, 6 hours after training a high gene expression was found at day 9 (compared to day 1) in both regions, suggesting that only one day of training is not sufficient to detect transcriptional modifications several hours after the task. Thus, gene expression data match 2DG results indicating a transfer of information in time (from day 1 to day 9) and in space (from the hippocampus to the retrosplenial cortex), and this at a cellular and a molecular level. Moreover, learning seems to modify the neuron-glia metabolic coupling, since several genes involved in this coupling are induced. These results also suggest a role of glia in neuronal plasticity.RESUME DU TRAVAIL DE THESECe projet de thèse a eu pour but l'étude des mécanismes moléculaires qui sont impliqués dans l'apprentissage et la mémoire et, en particulier, à les mettre en rapport avec le couplage métabolique existant entre les astrocytes et les neurones. Pour cela, des changements de l'activité métabolique dans différentes régions du cerveau des souris après 1 ou 9 jours d'entraînement dans un labyrinthe radial à huit-bras ont été évalués par autoradiographie au 2-désoxyglucose (2DG). Des différences significatives dans les régions engagées pendant la tâche comportementale au jour 1 (quand les animaux sont confrontés pour la première fois à la tâche) et au jour 9 (quand les animaux ont déjà appris) ont été identifiés. Ces régions incluent, au jour 1, l'hippocampe, le fornix, le cortex pariétal, le noyau thalamic laterodorsal et les corps mamillaires; et, au jour 9, le cingulaire antérieur, le cortex retrosplenial et le striatum dorsal. Deux de ces régions cérébrales (celles présentant les plus grands changements à jour 1 et à jour 9: l'hippocampe et le cortex retrosplenial, respectivement) ont été découpées par microdissection au laser et quelques gènes liés au couplage métabolique neurone-glie, au métabolisme du glucose et à la plasticité synaptique ont été analysées par RT-PCR. L'étude 2DG et l'analyse de l'expression de gènes ont été exécutés à trois temps différents: 1) juste après entraînement, 2) 45 minutes et 3) 6 heures après la fin de la tâche. L'objectif principal de cette étude était l'identification des adaptations métaboliques suivant la tâche d'apprentissage. Les résultats de l'expression de gènes démontrent que la tâche d'apprentissage module profondément le profile d'expression des gènes dans le temps, signifiant que ces deux régions cérébrales avec un signal 2DG élevé (l'hippocampe et le cortex retrosplenial) ont adapté leurs « machines moléculaires » en conséquence. Presque tous les gènes étudiés montrent une expression plus élevée dans l'hippocampe au jour 1 comparé au jour 9, alors qu'une expression accrue a été trouvée dans le cortex retrosplenial au jour 9. Nous pouvons observer ces adaptations moléculaires avec un retard court de 45 minutes après la fin de la tâche. Cependant, 6 heures après l'entraînement, une expression de gènes élevée a été trouvée au jour 9 (comparé à jour 1) dans les deux régions, suggérant que seulement un jour d'entraînement ne suffit pas pour détecter des modifications transcriptionelles plusieurs heures après la tâche. Ainsi, les données d'expression de gènes corroborent les résultats 2DG indiquant un transfert d'information dans le temps (de jour 1 à jour 9) et dans l'espace (de l'hippocampe au cortex retrosplenial), et ceci à un niveau cellulaire et moléculaire. D'ailleurs, la tâche d'apprentissage semble modifier le couplage métabolique neurone-glie, puisque de nombreux gènes impliqués dans ce couplage sont induits. Ces observations suggèrent un rôle important de la glie dans les mécanismes de plasticité du système nerveux.

Temporal quantification of MAPK induced expression in single yeast cells.

The quantification of gene expression at the single cell level uncovers novel regulatory mechanisms obscured in measurements performed at the population level. Two methods based on microscopy and flow cytometry are presented to demonstrate how such data can be acquired. The expression of a fluorescent reporter induced upon activation of the high osmolarity glycerol MAPK pathway in yeast is used as an example. The specific advantages of each method are highlighted. Flow cytometry measures a large number of cells (10,000) and provides a direct measure of the dynamics of protein expression independent of the slow maturation kinetics of the fluorescent protein. Imaging of living cells by microscopy is by contrast limited to the measurement of the matured form of the reporter in fewer cells. However, the data sets generated by this technique can be extremely rich thanks to the combinations of multiple reporters and to the spatial and temporal information obtained from individual cells. The combination of these two measurement methods can deliver new insights on the regulation of protein expression by signaling pathways.

Tissue-specific FLAGELLIN-SENSING 2 (FLS2) expression in roots restores immune responses in Arabidopsis fls2 mutants.

The flagellin receptor of Arabidopsis, At-FLAGELLIN SENSING 2 (FLS2), has become a model for mechanistic and functional studies on plant immune receptors. Responses to flagellin or its active epitope flagellin 22 (flg22) have been extensively studied in Arabidopsis leaves. However, the perception of microbe-associated molecular patterns (MAMPs) and the immune responses in roots are poorly understood. Here, we show that isolated root tissue is able to induce pattern-triggered immunity (PTI) responses upon flg22 perception, in contrast to elf18 (the active epitope of elongation factor thermo unstable (EF-Tu)). Making use of fls2 mutant plants and tissue-specific promoters, we generated transgenic Arabidopsis lines expressing FLS2 only in certain root tissues. This allowed us to study the spatial requirements for flg22 responses in the root. Remarkably, the intensity of the immune responses did not always correlate with the expression level of the FLS2 receptor, but depended on the expressing tissue, supporting the idea that MAMP perception and sensitivity in different tissues contribute to a proper balance of defense responses according to the expected exposure to elicitors. In summary, we conclude that each investigated root tissue is able to perceive flg22 if FLS2 is present and that tissue identity is a major element of MAMP perception in roots.

Revisiting spatial distribution and biochemical composition of calcium-containing crystals in human osteoarthritic articular cartilage.

INTRODUCTION: Calcium-containing (CaC) crystals, including basic calcium phosphate (BCP) and calcium pyrophosphate dihydrate (CPP), are associated with destructive forms of osteoarthritis (OA). We assessed their distribution and biochemical and morphologic features in human knee OA cartilage. METHODS: We prospectively included 20 patients who underwent total knee replacement (TKR) for primary OA. CaC crystal characterization and identification involved Fourier-transform infra-red spectrometry and scanning electron microscopy of 8 to 10 cartilage zones of each knee, including medial and lateral femoral condyles and tibial plateaux and the intercondyle zone. Differential expression of genes involved in the mineralization process between cartilage with and without calcification was assessed in samples from 8 different patients by RT-PCR. Immunohistochemistry and histology studies were performed in 6 different patients. RESULTS: Mean (SEM) age and body mass index of patients at the time of TKR was 74.6 (1.7) years and 28.1 (1.6) kg/m², respectively. Preoperative X-rays showed joint calcifications (chondrocalcinosis) in 4 cases only. The medial femoro-tibial compartment was the most severely affected in all cases, and mean (SEM) Kellgren-Lawrence score was 3.8 (0.1). All 20 OA cartilages showed CaC crystals. The mineral content represented 7.7% (8.1%) of the cartilage weight. All patients showed BCP crystals, which were associated with CPP crystals for 8 joints. CaC crystals were present in all knee joint compartments and in a mean of 4.6 (1.7) of the 8 studied areas. Crystal content was similar between superficial and deep layers and between medial and femoral compartments. BCP samples showed spherical structures, typical of biological apatite, and CPP samples showed rod-shaped or cubic structures. The expression of several genes involved in mineralization, including human homolog of progressive ankylosis, plasma-cell-membrane glycoprotein 1 and tissue-nonspecific alkaline phosphatase, was upregulated in OA chondrocytes isolated from CaC crystal-containing cartilages. CONCLUSIONS: CaC crystal deposition is a widespread phenomenon in human OA articular cartilage involving the entire knee cartilage including macroscopically normal and less weight-bearing zones. Cartilage calcification is associated with altered expression of genes involved in the mineralisation process.