The Gly482Ser genotype at the PPARGC1A gene and elevated blood pressure: a meta-analysis involving 13,949 individuals.

The protein encoded by the PPARGC1A gene is expressed at high levels in metabolically active tissues and is involved in the control of oxidative stress via reactive oxygen species detoxification. Several recent reports suggest that the PPARGC1A Gly482Ser (rs8192678) missense polymorphism may relate inversely with blood pressure. We used conventional meta-analysis methods to assess the association between Gly482Ser and systolic (SBP) or diastolic blood pressures (DBP) or hypertension in 13,949 individuals from 17 studies, of which 6,042 were previously unpublished observations. The studies comprised cohorts of white European, Asian, and American Indian adults, and adolescents from South America. Stratified analyses were conducted to control for population stratification. Pooled genotype frequencies were 0.47 (Gly482Gly), 0.42 (Gly482Ser), and 0.11 (Ser482Ser). We found no evidence of association between Gly482Ser and SBP [Gly482Gly: mean = 131.0 mmHg, 95% confidence interval (CI) = 130.5-131.5 mmHg; Gly482Ser mean = 133.1 mmHg, 95% CI = 132.6-133.6 mmHg; Ser482Ser: mean = 133.5 mmHg, 95% CI = 132.5-134.5 mmHg; P = 0.409] or DBP (Gly482Gly: mean = 80.3 mmHg, 95% CI = 80.0-80.6 mmHg; Gly482Ser mean = 81.5 mmHg, 95% CI = 81.2-81.8 mmHg; Ser482Ser: mean = 82.1 mmHg, 95% CI = 81.5-82.7 mmHg; P = 0.651). Contrary to previous reports, we did not observe significant effect modification by sex (SBP, P = 0.966; DBP, P = 0.715). We were also unable to confirm the previously reported association between the Ser482 allele and hypertension [odds ratio: 0.97, 95% CI = 0.87-1.08, P = 0.585]. These results were materially unchanged when analyses were focused on whites only. However, statistical evidence of gene-age interaction was apparent for DBP [Gly482Gly: 73.5 (72.8, 74.2), Gly482Ser: 77.0 (76.2, 77.8), Ser482Ser: 79.1 (77.4, 80.9), P = 4.20 x 10(-12)] and SBP [Gly482Gly: 121.4 (120.4, 122.5), Gly482Ser: 125.9 (124.6, 127.1), Ser482Ser: 129.2 (126.5, 131.9), P = 7.20 x 10(-12)] in individuals <50 yr (n = 2,511); these genetic effects were absent in those older than 50 yr (n = 5,088) (SBP, P = 0.41; DBP, P = 0.51). Our findings suggest that the PPARGC1A Ser482 allele may be associated with higher blood pressure, but this is only apparent in younger adults.

Absence of Association of Metabolic Syndrome with PPARGC1A, PPARG and UCP1 Gene Polymorphisms in Asian Indians

The objective of this study was to evaluate the association of PPARG coactivator1 alpha (PPARGC1A), peroxisome proliferator activated receptor gamma (PPARG), and uncoupling protein1 (UCP1) gene polymorphisms with the metabolic syndrome (MS) in an Asian Indian population. Nine common polymorphisms were genotyped via polymerase chain reaction restriction fragment length polymorphism and direct sequencing in 950 normal glucose-tolerant subjects and 550 type 2 diabetic subjects, chosen randomly from the Chennai Urban Rural Epidemiological Study, an ongoing population based study in Southern India. Among the 9 polymorphisms examined, only the Thr394Thr variant of the PPARGC1A gene was significantly associated with diabetes and obesity. The genotype frequency of GA of Thr394Thr variant was 16% (138/887) in the nonMS group and 22% (136/613) in the MS group, and this genotype frequency was significantly higher with MS both in males (p = 0.01) and females (p = 0.05), compared to the without-MS group. Logistic regression analysis revealed that the odds ratio for MS for the susceptible genotype GA of Thr394Thr was 1.411 [95% CI: 1.03-1.84, p = 0.012]. In the multiple logistic regression analysis, however, there was no association of this polymorphism as an independent factor with MS. Hence, the study shows that the polymorphisms in the PPARGC1A, PPARG and UCP1 genes are not associated with MS in Asian Indians.

Habitual energy expenditure modifies the association between NOS3 gene polymorphisms and blood pressue

BACKGROUND: The endothelial nitric-oxide synthase (NOS3) gene encodes the enzyme (eNOS) that synthesizes the molecule nitric oxide, which facilitates endothelium-dependent vasodilation in response to physical activity. Thus, energy expenditure may modify the association between the genetic variation at NOS3 and blood pressure. METHODS: To test this hypothesis, we genotyped 11 NOS3 polymorphisms, capturing all common variations, in 726 men and women from the Medical Research Council (MRC) Ely Study (age (mean +/- s.d.): 55 +/- 10 years, body mass index: 26.4 +/- 4.1 kg/m(2)). Habitual/non-resting energy expenditure (NREE) was assessed via individually calibrated heart rate monitoring over 4 days. RESULTS: The intronic variant, IVS25+15 [G-->A], was significantly associated with blood pressure; GG homozygotes had significantly lower levels of diastolic blood pressure (DBP) (-2.8 mm Hg; P = 0.016) and systolic blood pressure (SBP) (-1.9 mm Hg; P = 0.018) than A-allele carriers. The interaction between NREE and IVS25+15 was also significant for both DBP (P = 0.006) and SBP (P = 0.026), in such a way that the effect of the GG-genotype on blood pressure was stronger in individuals with higher NREE (DBP: -4.9 mm Hg, P = 0.02. SBP: -3.8 mm Hg, P= 0.03 for the third tertile). Similar results were observed when the outcome was dichotomously defined as hypertension. CONCLUSIONS: In summary, the NOS3 IVS25+15 is directly associated with blood pressure and hypertension in white Europeans. However, the associations are most evident in the individuals with the highest NREE. These results need further replication and have to be ideally tested in a trial before being informative for targeted disease prevention. Eventually, the selection of individuals for lifestyle intervention programs could be guided by knowledge of genotype.

Thr54 allele carriers of the Ala54Thr variant of FABP2 gene have associations with metabolic syndrome and hypertriglyceridemia in urban South Indians

The intestinal fatty acid-binding protein gene is proposed as a candidate gene for diabetes because the protein it codes is involved in fatty acid absorption and metabolism. This study investigates the association of the Ala54Thr variant of the intestinal fatty acid-binding protein gene on type 2 diabetes mellitus and other related metabolic traits in Asian Indians. Ala54Thr polymorphism was genotyped by using polymerase chain reaction-restriction fragment length polymorphism in unrelated 773 type 2 diabetic and 899 normal glucose-tolerant (NGT) subjects, randomly chosen from the Chennai Urban Rural Epidemiology Study, an ongoing population-based study in South India. The Ala54Thr polymorphism was not associated with type 2 diabetes mellitus or obesity. However, genotype-phenotype study revealed that the NGT subjects carrying the Thr54 allele had significantly higher 2-hour plasma glucose (P = .007), glycated hemoglobin (P = .004), 2-hour insulin (P = .027), and fasting low-density lipoprotein cholesterol (P = .032) levels compared with those with the Ala54 allele. Normal glucose-tolerant subjects with Ala54Thr and Thr54Thr genotypes had significantly higher fasting serum triglyceride levels (P = .003) compared with those with Ala54Ala. The subjects were stratified into those with hypertriglyceridemia (serum triglyceride levels >or=150 mg/dL) and those without. The odds ratio for hypertriglyceridemia for the individuals carrying the Ala54Thr genotype was 1.491 (95% confidence interval [CI], 1.22-1.83, P < .0001), and for those carrying the Thr54Thr genotype, it was 1.888 (95% CI, 1.34-2.67; P < .0001). Subjects were also stratified into those with metabolic syndrome (MS) and those without, according to modified Adult Treatment Panel III guidelines. The odds ratio (adjusted for age and sex) for MS for the individuals carrying the Ala54Thr genotype was 1.240 (95% CI, 1.02-1.51; P = .03), whereas for those carrying the Thr54Thr genotype, it was 1.812 (95% CI, 1.28-2.57; P = .001). Carriers of the Thr54 allele have associations with MS and hypertriglyceridemia in this urban South Indian population.

Effect of polymorphisms in the PPARGC1A gene on body fat in Asian Indians

OBJECTIVE: To evaluate whether polymorphisms in the peroxisome proliferator-activated receptor-gamma coactivator-1 alpha (PPARGC1A) gene were related to body fat in Asian Indians. METHODS: Three polymorphisms of PPARGC1A gene, the Thr394Thr, Gly482Ser and +A2962G, were genotyped on 82 type 2 diabetic and 82 normal glucose tolerant (NGT) subjects randomly chosen from the Chennai Urban Rural Epidemiology Study using PCR-RFLP, and the nature of the variants were confirmed using direct sequencing. Linkage disequilibrium (LD) was estimated from the estimates of haplotypic frequencies using an expectation-maximization algorithm. Visceral, subcutaneous and total abdominal fat were measured using computed tomography, whereas dual X-ray absorptiometry was used to measure central abdominal and total body fat. RESULTS: None of the three polymorphisms studied were in LD. The genotype (0.59 vs 0.32, P=0.001) and allele (0.30 vs 0.17, P=0.007) frequencies of Thr394Thr polymorphism were significantly higher in type 2 diabetic subjects compared to those in NGT subjects. The odds ratio for diabetes (adjusted for age, sex and body mass index) for the susceptible genotype, XA (GA+AA) of Thr394Thr polymorphism, was 2.53 (95% confidence intervals: 1.30-5.04, P=0.009). Visceral and subcutaneous fat were significantly higher in NGT subjects with XA genotype of the Thr394Thr polymorphism compared to those with GG genotype (visceral fat: XA 148.2+/-46.9 vs GG 106.5+/-51.9 cm(2), P=0.001; subcutaneous fat: XA 271.8+/-167.1 vs GG 181.5+/-78.5 cm(2), P=0.001). Abdominal (XA 4521.9+/-1749.6 vs GG 3445.2+/-1443.4 g, P=0.004), central abdominal (XA 1689.0+/-524.0 vs GG 1228.5+/-438.7 g, P<0.0001) and non-abdominal fat (XA 18763.8+/-8789.4 vs GG 13160.4+/-4255.3 g, P<0.0001) were also significantly higher in the NGT subjects with XA genotype compared to those with GG genotype. The Gly482Ser and +A2962G polymorphisms were not associated with any of the body fat measures. CONCLUSION: Among Asian Indians, the Thr394Thr (G --> A) polymorphism is associated with increased total, visceral and subcutaneous body fat.

Genetic variant rs17225178 in the ARNT2 gene is associated with Asperger Syndrome

Background Autism Spectrum Conditions (ASC) are neurodevelopmental conditions characterized by difficulties in communication and social interaction, alongside unusually repetitive behaviours and narrow interests. Asperger Syndrome (AS) is one subgroup of ASC and differs from classic autism in that in AS there is no language or general cognitive delay. Genetic, epigenetic and environmental factors are implicated in ASC and genes involved in neural connectivity and neurodevelopment are good candidates for studying the susceptibility to ASC. The aryl-hydrocarbon receptor nuclear translocator 2 (ARNT2) gene encodes a transcription factor involved in neurodevelopmental processes, neuronal connectivity and cellular responses to hypoxia. A mutation in this gene has been identified in individuals with ASC and single nucleotide polymorphisms (SNPs) have been nominally associated with AS and autistic traits in previous studies. Methods In this study, we tested 34 SNPs in ARNT2 for association with AS in 118 cases and 412 controls of Caucasian origin. P values were adjusted for multiple comparisons, and linkage disequilibrium (LD) among the SNPs analysed was calculated in our sample. Finally, SNP annotation allowed functional and structural analyses of the genetic variants in ARNT2. We tested the replicability of our result using the genome-wide association studies (GWAS) database of the Psychiatric Genomics Consortium (PGC). Results We report statistically significant association of rs17225178 with AS. This SNP modifies transcription factor binding sites and regions that regulate the chromatin state in neural cell lines. It is also included in a LD block in our sample, alongside other genetic variants that alter chromatin regulatory regions in neural cells. Conclusions These findings demonstrate that rs17225178 in the ARNT2 gene is associated with AS and support previous studies that pointed out an involvement of this gene in the predisposition to ASC.

Ascl1 coordinately regulates gene expression and the chromatin landscape during neurogenesis

The proneural transcription factor Ascl1 coordinates gene expression in both proliferating and differentiating progenitors along the neuronal lineage. Here, we used a cellular model of neurogenesis to investigate how Ascl1 interacts with the chromatin landscape to regulate gene expression when promoting neuronal differentiation. We find that Ascl1 binding occurs mostly at distal enhancers and is associated with activation of gene transcription. Surprisingly, the accessibility of Ascl1 to its binding sites in neural stem/progenitor cells remains largely unchanged throughout their differentiation, as Ascl1 targets regions of both readily accessible and closed chromatin in proliferating cells. Moreover, binding of Ascl1 often precedes an increase in chromatin accessibility and the appearance of new regions of open chromatin, associated with de novo gene expression during differentiation. Our results reveal a function of Ascl1 in promoting chromatin accessibility during neurogenesis, linking the chromatin landscape at Ascl1 target regions with the temporal progression of its transcriptional program.

Crosstalk between oestrogen receptors and thyroid hormone receptor isoforms results in differential regulation of the preproenkephalin gene

Nuclear receptors are ligand-activated transcription factors, which have the potential to integrate internal metabolic events in an organism, with consequences for control of behaviour. Previous studies from this laboratory have shown that thyroid hormone receptor (TR) isoforms can inhibit oestrogen receptor (ER)alpha-mediated induction of preproenkephalin (PPE) gene expression in the hypothalamus. Also, thyroid hormone administration inhibits lordosis, a behaviour facilitated by PPE expression. We have examined the effect of multiple ligand-binding TR isoforms on the ER-mediated induction of the PPE gene in transient transfection assays in CV-1 cells. On a natural PPE gene promoter fragment containing two putative oestrogen response elements (EREs), both ER alpha and beta isoforms mediate a four to five-fold induction by oestrogen. Cotransfection of TR alpha 1 along with ER alpha inhibited the ER alpha transactivation of PPE by approximately 50%. However, cotransfection with either TR beta 1 or TR beta 2 expression plasmids produced no effect on the ER alpha or ER beta mediated induction of PPE. Therefore, under these experimental conditions, interactions with a single ER isoform are specific to an individual TR isoform. Transfection with a TR alpha 1 DNA-binding mutant could also inhibit ER alpha transactivation, suggesting that competition for binding on the ERE may not be the exclusive mechanism for inhibition. Data with the coactivator, SRC-1, suggested that coactivator squelching may participate in the inhibition. In dramatic contrast, when ER beta is cotransfected, TR alpha 1 stimulated ER beta-mediated transactivation of PPE by approximately eight-fold over control levels. This is the first study revealing specific interactions among nuclear receptor isoforms on a neuroendocrine promoter. These data also suggest that the combinatorics of ER and TR isoforms allow multiple forms of flexible gene regulations in the service of neuroendocrine integration.

Eag 1, Eag 2 and Kcnn3 gene brain expression of isolated reared rats

The Eag1 and Eag2, voltage-dependent potassium channels, and the small-conductance calcium-activated potassium channel (Kcnn3) are highly expressed in limbic regions of the brain, where their function is still unknown. Eag1 co-localizes with tyrosine hydroxilase enzyme in the substantia nigra and ventral tegmental area. Kcnn3 deficiency leads to enhanced serotonergic and dopaminergic neurotransmission accompanied by distinct alterations in emotional behaviors. As exposure to stress is able to change the expression and function of several ion channels, suggesting that they might be involved in the consequences of stress, we aimed at investigating Eag 1, Eag2 and Kcnn3 mRNA expression in the brains of rats submitted to isolation rearing. As the long-lasting alterations in emotional and behavioral regulation after stress have been related to changes in serotonergic neurotransmission, expressions of serotonin Htr1a and Htr2a receptors in male Wistar rats` brain were also investigated. Rats were reared in isolation or in groups of five for nine weeks after weaning. Isolated and socially reared rats were tested for exploratory activity in the open field test for 5 min and brains were processed for reverse-transcription coupled to quantitative polymerase chain reaction (qRT-PCR). Isolated reared rats showed decreased exploratory activity in the open field. Compared to socially reared rats, isolated rats showed reduced Htr2a mRNA expression in the striatum and brainstem and reduced Eag2 mRNA expression in all examined regions except cerebellum. To our knowledge, this is the first work to show that isolation rearing can change Eag2 gene expression in the brain. The involvement of this channel in stress-related behaviors is discussed.

Taxonomic and Phylogenetic Relationships Between Species of the Genus Culex (Diptera: Culicidae) From Brazil Inferred From the Cytochrome c Oxidase I Mitochondrial Gene

Species of the genus Culex Linnaeus have been incriminated as the main vectors of lymphatic filariases and are important vectors of arboviruses, including West Nile virus. Sequences corresponding to a fragment of 478 bp of the cytochrome c oxidase subunit I gene, which includes part of the barcode region, of 37 individuals of 17 species of genus Culex were generated to establish relationships among five subgenera, Culex, Phenacomyia, Melanoconion, Microculex, and Carrollia, and one species of the genus Lutzia that occurs in Brazil. Bayesian methods were employed for the phylogenetic analyses. Results of sequence comparisons showed that individuals identified as Culex dolosus, Culex mollis, and Culex imitator possess high intraspecific divergence (3.1, 2.3, and 3.5%, respectively) when using the Kimura two parameters model. These differences were associated either with distinct morphological characteristics of the male genitalia or larval and pupal stages, suggesting that these may represent species complexes. The Bayesian topology suggested that the genus and subgenus Culex are paraphyletic relative to Lutzia and Phenacomyia, respectively. The cytochrome c oxidase subunit I sequences may be a useful tool to both estimate phylogenetic relationships and identify morphologically similar species of the genus Culex.

Phylogenetic relationships among species of Anopheles (Nyssorhynchus) (Diptera, Culicidae) based on nuclear and mitochondrial gene sequences

Phylogenetic relationships among 21 species of mosquitoes in subgenus Nyssorhynchus were inferred from the nuclear white and mitochondrial NADH dehydrogenase subunit 6 (ND6) genes. Bayestan phylogenetic methods found that none of the three Sections within Nyssorhynchus (Albimanus, Argyritarsis, Myzorhynchella) were supported in all analyses, although Myzorhynchella was found to be monophyletic at the combined genes Within the Albimanus Section the monophyly of the Stroder Subgroup was strongly supported and within the Myzorhynchella Section Anopheles anrunesi and An lutzu formed a strongly supported monophyletic group The epidemiologically significant Albitarsis Complex showed evidence of paraphyly (relative to An lanet-Myzorhynchella) and discordance across gene trees, and the previously synonomized species of An. dunhami and An goeldii were recovered as sister species Finally, there was evidence of complexes in several species, including An antunesi, An deaneorum, and An. strodei (c) 2010 Elsevier B.V. All rights reserved

Association of genetic variants in the adiponectin encoding gene (ADIPOQ) with type 2 diabetes in Japanese Brazilians

Aim: The objective of this study is to assess the contribution of ADIPOQ variants to type 2 diabetes in Japanese Brazilians. Methods: We genotyped 200 patients with diabetes mellitus (100 male and 100 female, aged 55.0 years [47.5-64.0 years]) and 200 control subjects with normal glucose tolerant (NGT) (72 male and 128 female, aged 52.0 years [43.5-64.5 years]). Results: Whereas each polymorphism studied (T45G, G276T, and A349G) was not significantly associated with type 2 diabetes mellitus, the haplotype GGA was overrepresented in our diabetic population (9.3% against 3.1% in NGT individuals, P=.0003). Also, this haplotype was associated with decreased levels of adiponectin. We also identified three mutations in exon 3: I164T, R221S, and H241P, but, owing to the low frequencies of them, associations with type 2 diabetes could not be evaluated. The subjects carrying the R221S mutation had plasma adiponectin levels lower than those without the mutation (2.10 mu g/ml [1.35-2.55 mu g/ml] vs. 6.68 mu g/ml [3.90-11.23 mu g/ml], P=.015). Similarly, the I164T mutation carriers had mean plasma adiponectin levels lower than those noncarriers (3.73 mu g/ml [3.10-4.35 mu g/ml] vs. 6.68 mu g/ml [3.90-11.23 mu g/ml]), but this difference was not significant (P=.17). Conclusions: We identified in the ADIPOQ gene a risk haplotype for type 2 diabetes in the Japanese Brazilian population. (C) 2010 Elsevier Inc. All rights reserved.

Clinical evaluation and COL2A1 gene analysis in 21 Brazilian families with Stickler syndrome: Identification of novel mutations, further genotype/phenotype correlation, and its implications for the diagnosis

We present clinical and molecular evaluation from a large cohort of patients with Stickler syndrome: 78 individuals from 21 unrelated Brazilian families. The patients were selected in a Hospital with a craniofacial dysmorphology assistance service and clinical diagnosis was based on the presence of cleft palate associated to facial and ocular anomalies of Stickler syndrome. Analysis of COL2A1 gene revealed 9 novel and 4 previously described pathogenic mutations. Except for the mutation c.556G>T (p.Gly186X), all the others were located in the triple helical domain. We did not find genotype/phenotype correlation in relation to type and position of the mutation in the triple helical domain. However, a significantly higher proportion of myopia in patients with mutations located in this domain was observed in relation to those with the mutation in the non-tripe helical domain (c.556G>T; P < 0.04). A trend towards a higher prevalence of glaucoma, although not statistically significant, was observed in the presence of the mutation c.556G>T. It is possible. that this mutation alters the splicing of the mRNA instead of only creating a premature stop codon and therefore it can lead to protein products of different ocular effects. One novel DNA variation (c.1266+7G>C) occurs near a splice site and it was observed to co-segregate with the phenotype in one of the two families with this DNA variation. As in silico analysis predicted that the c.1266+7G>C DNA variation can affect the efficiency of the splicing, we still cannot rule it out as non-pathogenic. Our study also showed that ascertainment through cleft palate associated to other craniofacial signs can be very efficient for identification of Stickler syndrome patients. Still, high frequency of familial cases and high frequency of underdevelopment of distal lateral tibial epiphyses observed in our patients suggested that the inclusion of this information can improve the clinical diagnosis of Stickler syndrome. (C) 2008 Elsevier Masson SAS. All rights reserved.

Axillary bud development in pineapple nodal segments correlates with changes on cell cycle gene expression, hormone level, and sucrose and glutamate contents

During the process of lateral organ development after plant decapitation, cell division and differentiation occur in a balanced manner initiated by specific signaling, which triggers the reentrance into the cell cycle. Here, we investigated short-term variations in the content of some endogenous signals, such as auxin, cytokinins (Cks), and other mitogenic stimuli (sucrose and glutamate), which are likely correlated with the cell cycle reactivation in the axillary bud primordium of pineapple nodal segments. Transcript levels of cell cycle-associated genes, CycD2;1, and histone H2A were analyzed. Nodal segments containing the quiescent axillary meristem cells were cultivated in vitro during 24 h after the apex removal and de-rooting. From the moment of stem apex and root removal, decapitated nodal segment (DNS) explants showed a lower indol-3-acetic acid (IAA) concentration than control explants, and soon after, an increase of endogenous sucrose and iP-type Cks were detected. The decrease of IAA may be the primary signal for cell cycle control early in G1 phase, leading to the upregulation of CycD2;1 gene in the first h. Later, the iP-type Cks and sucrose could have triggered the progression to S-phase since there was an increase in H2A expression at the eighth h. DNS explants revealed substantial increase in Z-type Cks and glutamate from the 12th h, suggesting that these mitogens could also operate in promoting pineapple cell cycle progression. We emphasize that the use of non-synchronized tissue rather than synchronous cell suspension culture makes it more difficult to interpret the results of a dynamic cell division process. However, pineapple nodal segments cultivated in vitro may serve as an interesting model to shed light on apical dominance release and the reentrance of quiescent axillary meristem cells into the cell cycle.

Human Stem Cell Cultures from Cleft Lip/Palate Patients Show Enrichment of Transcripts Involved in Extracellular Matrix Modeling By Comparison to Controls

Nonsyndromic cleft lip and palate (NSCL/P) is a complex disease resulting from failure of fusion of facial primordia, a complex developmental process that includes the epithelial-mesenchymal transition (EMT). Detection of differential gene transcription between NSCL/P patients and control individuals offers an interesting alternative for investigating pathways involved in disease manifestation. Here we compared the transcriptome of 6 dental pulp stem cell (DPSC) cultures from NSCL/P patients and 6 controls. Eighty-seven differentially expressed genes (DEGs) were identified. The most significant putative gene network comprised 13 out of 87 DEGs of which 8 encode extracellular proteins: ACAN, COL4A1, COL4A2, GDF15, IGF2, MMP1, MMP3 and PDGFa. Through clustering analyses we also observed that MMP3, ACAN, COL4A1 and COL4A2 exhibit co-regulated expression. Interestingly, it is known that MMP3 cleavages a wide range of extracellular proteins, including the collagens IV, V, IX, X, proteoglycans, fibronectin and laminin. It is also capable of activating other MMPs. Moreover, MMP3 had previously been associated with NSCL/P. The same general pattern was observed in a further sample, confirming involvement of synchronized gene expression patterns which differed between NSCL/P patients and controls. These results show the robustness of our methodology for the detection of differentially expressed genes using the RankProd method. In conclusion, DPSCs from NSCL/P patients exhibit gene expression signatures involving genes associated with mechanisms of extracellular matrix modeling and palate EMT processes which differ from those observed in controls. This comparative approach should lead to a more rapid identification of gene networks predisposing to this complex malformation syndrome than conventional gene mapping technologies.