The respiratory syncytial virus small hydrophobic protein is phosphorylated via a mitogen-activated protein kinase p38-dependent tyrosine kinase activity during virus infection

The phosphorylation status of the small hydrophobic (SH) protein of respiratory syncytial virus (RSV) was examined in virus-infected Vero cells. The SH protein v.,as isolated from [S-35]methionine- and [P-33]orthophosphate-labelled IRSV-infected cells and analysed by SDS-PAGE. In each case, a protein product of the expected size for the SH protein was observed. Phosphoamino acid analysis and reactivity with the phosphotyrosine specific antibody PY20 showed that the SH protein was modified by tyrosine phosphorylation. The role or tyrosine kinase activity in SH protein phosphorylation was confirmed by the use of genistein, a broad-spectrum tyrosine kinase inhibitor, to inhibit SH protein phosphorylation. Further analysis showed that the different glycosylated forms of the SH protein were phosphorylated, as was the oligomeric form of the protein. Phosphorylation of the SH protein was specifically inhibited by the mitogen-activated protein kinase (MAPK) p38 inhibitor SB203580, suggesting that SH protein phosphorylation occurs via a MAPK p38-dependent pathway. Analysis of virus-infected cells using fluorescence microscopy showed that, although the SH protein was distributed throughout the cytoplasm, it appeared to accumulate, at low levels, in the endoplasmic reticulum/Golgi complex, confirming recent observations. However, in the presence of SB203580. an increased accumulation of the SH protein in the Golgi complex was observed, although other virus structures, such as virus filaments and inclusion bodies, remained largely unaffected. These results showed that during RSV infection, the SH protein is modified by an MAPK p38-dependant tyrosine kinase activity and that this modification influences its cellular distribution.

The Colorectal cancer disease-specific transcriptome may facilitate the discovery of more biologically and clinically relevant information

BACKGROUND: To date, there are no clinically reliable predictive markers of response to the current treatment regimens for advanced colorectal cancer. The aim of the current study was to compare and assess the power of transcriptional profiling using a generic microarray and a disease-specific transcriptome-based microarray. We also examined the biological and clinical relevance of the disease-specific transcriptome.

METHODS: DNA microarray profiling was carried out on isogenic sensitive and 5-FU-resistant HCT116 colorectal cancer cell lines using the Affymetrix HG-U133 Plus2.0 array and the Almac Diagnostics Colorectal cancer disease specific Research tool. In addition, DNA microarray profiling was also carried out on pre-treatment metastatic colorectal cancer biopsies using the colorectal cancer disease specific Research tool. The two microarray platforms were compared based on detection of probesets and biological information.

RESULTS: The results demonstrated that the disease-specific transcriptome-based microarray was able to out-perform the generic genomic-based microarray on a number of levels including detection of transcripts and pathway analysis. In addition, the disease-specific microarray contains a high percentage of antisense transcripts and further analysis demonstrated that a number of these exist in sense:antisense pairs. Comparison between cell line models and metastatic CRC patient biopsies further demonstrated that a number of the identified sense:antisense pairs were also detected in CRC patient biopsies, suggesting potential clinical relevance.

CONCLUSIONS: Analysis from our in vitro and clinical experiments has demonstrated that many transcripts exist in sense:antisense pairs including IGF2BP2, which may have a direct regulatory function in the context of colorectal cancer. While the functional relevance of the antisense transcripts has been established by many studies, their functional role is currently unclear; however, the numbers that have been detected by the disease-specific microarray would suggest that they may be important regulatory transcripts. This study has demonstrated the power of a disease-specific transcriptome-based approach and highlighted the potential novel biologically and clinically relevant information that is gained when using such a methodology.

Relative quantification of the proteomeic changes associated with the mycotoxin zearalenone in the H295R steroidogenesis model

Zearalenone (ZEN) is a mycotoxin with endocrine disrupting effects having vast economic implications in e.g. pig farming. Structurally, ZEN resembles 17b-estradiol, and thus is able to bind to estrogen receptors (ER) in target cells. Because of this, it is also classified as a non-steroidal estrogen, a phytoestrogen, a mycoestrogen, and a growth promoter. Quantitative proteomic analysis was undertaken using stable-isotope labeling by amino acids in cell culture (SILAC) upon exposure of the steroidogenesis cell model H295R with ZEN to elucidate its effect on protein regulation. ZEN significantly regulated 21 proteins, including proteins with known endocrine disrupting effects and several oncogenes. In addition, network analysis using Ingenuity Pathway Analysis showed that ZEN affected the oxidative phosphorylation pathway and the mitochondrial dysfunction pathway, both previously reported to be involved in endocrine dysfunction.

Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types

Background: MicroRNAs (miRNAs) are a class of small RNA molecules that regulate expression of specific mRNA targets. They can be released from cells, often encapsulated within extracellular vesicles (EVs), and therefore have the potential to mediate intercellular communication. It has been suggested that certain miRNAs may be selectively exported, although the mechanism has yet to be identified. Manipulation of the miRNA content of EVs will be important for future therapeutic applications. We therefore wished to assess which endogenous miRNAs are enriched in EVs and how effectively an overexpressed miRNA would be exported.

Results: Small RNA libraries from HEK293T cells and vesicles before or after transfection with a vector for miR-146a overexpression were analysed by deep sequencing. A subset of miRNAs was found to be enriched in EVs; pathway analysis of their predicted target genes suggests a potential role in regulation of endocytosis. RT-qPCR in additional cell types and analysis of publicly available data revealed that many of these miRNAs tend to be widely preferentially exported. Whilst overexpressed miR-146a was highly enriched both in transfected cells and their EVs, the cellular:EV ratios of endogenous miRNAs were not grossly altered. MiR-451 was consistently the most highly exported miRNA in many different cell types. Intriguingly, Argonaute2 (Ago2) is required for miR-451 maturation and knock out of Ago2 has been shown to decrease expression of other preferentially exported miRNAs (eg miR-150 and miR-142-3p).

Conclusion: The global expression data provided by deep sequencing confirms that specific miRNAs are enriched in EVs released by HEK293T cells. Observation of similar patterns in a range of cell types suggests that a common mechanism for selective miRNA export may exist.

Progressive lung cancer determined by expression profiling and transcriptional regulation

Clinically, our ability to predict disease outcome for patients with early stage lung cancer is currently poor. To address this issue, tumour specimens were collected at surgery from non-small cell lung cancer (NSCLC) patients as part of the European Early Lung Cancer (EUELC) consortium. The patients were followed-up for three years post-surgery and patients who suffered progressive disease (PD, tumour recurrence, metastasis or a second primary) or remained disease-free (DF) during follow-up were identified. RNA from both tumour and adjacent-normal lung tissue was extracted from patients and subjected to microarray expression profiling. These samples included 36 adenocarcinomas and 23 squamous cell carcinomas from both PD and DF patients. The microarray data was subject to a series of systematic bioinformatics analyses at gene, network and transcription factor levels. The focus of these analyses was 2-fold: firstly to determine whether there were specific biomarkers capable of differentiating between PD and DF patients, and secondly, to identify molecular networks which may contribute to the progressive tumour phenotype. The experimental design and analyses performed permitted the clear differentiation between PD and DF patients using a set of biomarkers implicated in neuroendocrine signalling and allowed the inference of a set of transcription factors whose activity may differ according to disease outcome. Potential links between the biomarkers, the transcription factors and the genes p21/CDKN1A and Myc, which have previously been implicated in NSCLC development, were revealed by a combination of pathway analysis and microarray meta-analysis. These findings suggest that neuroendocrine-related genes, potentially driven through p21/CDKN1A and Myc, are closely linked to whether or not a NSCLC patient will have poor clinical outcome.

Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes

Aims/hypothesis: An abnormal urinary albumin excretion rate (AER) is often the first clinically detectable manifestation of diabetic nephropathy. Our aim was to estimate the heritability and to detect genetic variation associated with elevated AER in patients with type 1 diabetes. 
 Methods: The discovery phase genome-wide association study (GWAS) included 1,925 patients with type 1 diabetes and with data on 24 h AER. AER was analysed as a continuous trait and the analysis was stratified by the use of antihypertensive medication. Signals with a p value <10^{−4 were followed up in 3,750 additional patients with type 1 diabetes from seven studies.} 
 Results: The narrow-sense heritability, captured with our genotyping platform, was estimated to explain 27.3% of the total AER variability, and 37.6% after adjustment for covariates. In the discovery stage, five single nucleotide polymorphisms in the GLRA3 gene were strongly associated with albuminuria (p < 5 × 10^{−8). In the replication group, a nominally significant association (p = 0.035) was observed between albuminuria and rs1564939 in GLRA3, but this was in the opposite direction. Sequencing of the surrounding genetic region in 48 Finnish and 48 UK individuals supported the possibility that population-specific rare variants contribute to the synthetic association observed at the common variants in GLRA3. The strongest replication (p = 0.026) was obtained for rs2410601 between the PSD3 and SH2D4A genes. Pathway analysis highlighted natural killer cell mediated immunity processes.} 
 Conclusions/interpretation: This study suggests novel pathways and molecular mechanisms for the pathogenesis of albuminuria in type 1 diabetes.

Using SILAC proteomics to investigate the effect of the mycotoxin, alternariol, in the human H295R steroidogenesis model

The mycotoxin alternariol (AOH) is an important contaminant of fruits and cereal products. The current study sought to address the effect of a non-toxic AOH concentration on the proteome of the steroidogenic H295R cell model. Quantitative proteomics based on stable isotope labeling by amino acids in cell culture (SILAC) coupled to 1D-SDS-PAGE-LC-MS/MS was applied to subcellular-enriched protein samples. Gene ontology (GO) and ingenuity pathway analysis (IPA) were further carried out for functional annotation and identification of protein interaction networks. Furthermore, the effect of AOH on apoptosis and cell cycle distribution was also determined by the use of flow cytometry analysis. This work identified 22 proteins that were regulated significantly. The regulated proteins are those involved in early stages of steroid biosynthesis (SOAT1, NPC1, and ACBD5) and C21-steroid hormone metabolism (CYP21A2 and HSD3B1). In addition, several proteins known to play a role in cellular assembly, organization, protein synthesis, and cell cycle were regulated. These findings provide a new framework for studying the mechanisms by which AOH modulates steroidogenesis in H295R cell model. 

Convergent genetic and expression data implicate immunity in Alzheimer's disease

Background: Late-onset Alzheimer's disease (AD) is heritable with 20 genes showing genome-wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease, we extended these genetic data in a pathway analysis.

Methods: The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain.

Results: ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response (P = 3.27 X 10(-12) after multiple testing correction for pathways), regulation of endocytosis (P = 1.31 X 10(-11)), cholesterol transport (P = 2.96 X 10(-9)), and proteasome-ubiquitin activity (P = 1.34 X 10(-6)). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected P = .002-.05).

Conclusions: The immime response, regulation of endocytosis, cholesterol transport, and protein ubiquitination represent prime targets for AD therapeutics. (C) 2015 Published by Elsevier Inc. on behalf of The Alzheimer's Association.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

Genetically Determined Height and Coronary Artery Disease

Background: The nature and underlying mechanisms of an inverse association between adult height and the risk of coronary artery disease (CAD) are unclear. 
Methods: We used a genetic approach to investigate the association between height and CAD, using 180 height-associated genetic variants. We tested the association between a change in genetically determined height of 1 SD (6.5 cm) with the risk of CAD in 65,066 cases and 128,383 controls. Using individual-level genotype data from 18,249 persons, we also examined the risk of CAD associated with the presence of various numbers of height-associated alleles. To identify putative mechanisms, we analyzed whether genetically determined height was associated with known cardiovascular risk factors and performed a pathway analysis of the height-associated genes. 
Results: We observed a relative increase of 13.5% (95% confidence interval [CI], 5.4 to 22.1; P<0.001) in the risk of CAD per 1-SD decrease in genetically determined height. There was a graded relationship between the presence of an increased number of height-raising variants and a reduced risk of CAD (odds ratio for height quartile 4 versus quartile 1, 0.74; 95% CI, 0.68 to 0.84; P<0.001). Of the 12 risk factors that we studied, we observed significant associations only with levels of low-density lipoprotein cholesterol and triglycerides (accounting for approximately 30% of the association). We identified several overlapping pathways involving genes associated with both development and atherosclerosis. 
Conclusions: There is a primary association between a genetically determined shorter height and an increased risk of CAD, a link that is partly explained by the association between shorter height and an adverse lipid profile. Shared biologic processes that determine achieved height and the development of atherosclerosis may explain some of the association. (Funded by the British Heart Foundation and others.)

Alzheimer's disease-like pathology has transient effects on the brain and blood metabolome

The pathogenesis of Alzheimer's disease (AD) is complex involving multiple contributing factors. The extent to which AD pathology impacts upon the metabolome is still not understood, nor is it known how disturbances change as the disease progresses. For the first time we have profiled longitudinally (6, 8, 10, 12 and 18 months) both the brain and plasma metabolome of APP/PS1 double transgenic and wild type (WT) mice. A total of 187 metabolites were quantified using a targeted metabolomics methodology. Multivariate statistical analysis produced models that distinguished APP/PS1 from WT mice at 8, 10 and 12 months.Metabolic pathway analysis found perturbed polyamine metabolism in both brain and blood plasma. There were other disturbances in essential amino acids,branched chain amino acids and also in the neurotransmitter serotonin.Pronounced imbalances in phospholipid and acylcarnitine homeostasis was evident in two age groups. AD-like pathology therefore impacts greatly on both the brain and blood metabolomes, although there appears to be a clear temporal sequence whereby changes to brain metabolites precede those in blood.

Pine wilt disease: a worldwide threat to forest ecosystems

The pinewood nematode (PWN), Bursaphelenchus xylophilus, the causal agent of pine wilt disease (PWD), is a serious pest and pathogen of forest tree species, in particular among the genus Pinus. It was first reported from Japan in the beginning of the XXth century, where it became the major ecological catastrophe of pine forests, with losses reaching over 2 million m3/ year in the 1980s. It has since then spread to other Asian countries such as China, Taiwan and Korea, causing serious losses and economic damage. In 1999, the PWN was first detected in the European Union (EU), in Portugal, and immmediately prompted several government (national and EU) actions to assess the extent of the nematode’s presence, and to contain B. xylophilus and its insect vector (Monochamus galloprovincialis) to an area with a 30km radius in the Setúbal Peninsula, 20 km south of Lisbon. International wood trade, with its political as well as economic ramifications, has been seriously jeopardized. The origin of the population of PWN found in Portugal remains elusive. Several hypotheses may be considered regarding pathway analysis, basically from two general origins: North America or the Far East (Japan or China). World trade of wood products such as timber, wooden crates, palettes, etc… play an important role in the potential dissemination of the pinewood nematode. In fact, human activities involving the movement of wood products may be considered the single most important factor in spreading of the PWN. Despite the dedicated and concerted actions of government agencies, this disease continues to spread. Very recently (2006), in Portugal, forestry and phytosanitary authorities (DGRF and DGPC) have announced a new strategy for the control and ultimately the erradication of the nematode, under the coordination of the national program for the control of the pinewood nematode (PROLUNP). Research regarding the bioecology of the nematode and insect as well as new detection methods, e.g., involving real-time PCR, has progressed since 1999. International agreements (GATT, WTO) and sharing of scientific information is of paramount importance to effectively control the nematode and its vector, and thus protect our forest ecosystems and forest economy.

Modeling and optimization of extracellular polysaccharides production by Enterobacter A47

Polysaccharides are gaining increasing attention as potential environmental friendly and sustainable building blocks in many fields of the (bio)chemical industry. The microbial production of polysaccharides is envisioned as a promising path, since higher biomass growth rates are possible and therefore higher productivities may be achieved compared to vegetable or animal polysaccharides sources. This Ph.D. thesis focuses on the modeling and optimization of a particular microbial polysaccharide, namely the production of extracellular polysaccharides (EPS) by the bacterial strain Enterobacter A47. Enterobacter A47 was found to be a metabolically versatile organism in terms of its adaptability to complex media, notably capable of achieving high growth rates in media containing glycerol byproduct from the biodiesel industry. However, the industrial implementation of this production process is still hampered due to a largely unoptimized process. Kinetic rates from the bioreactor operation are heavily dependent on operational parameters such as temperature, pH, stirring and aeration rate. The increase of culture broth viscosity is a common feature of this culture and has a major impact on the overall performance. This fact complicates the mathematical modeling of the process, limiting the possibility to understand, control and optimize productivity. In order to tackle this difficulty, data-driven mathematical methodologies such as Artificial Neural Networks can be employed to incorporate additional process data to complement the known mathematical description of the fermentation kinetics. In this Ph.D. thesis, we have adopted such an hybrid modeling framework that enabled the incorporation of temperature, pH and viscosity effects on the fermentation kinetics in order to improve the dynamical modeling and optimization of the process. A model-based optimization method was implemented that enabled to design bioreactor optimal control strategies in the sense of EPS productivity maximization. It is also critical to understand EPS synthesis at the level of the bacterial metabolism, since the production of EPS is a tightly regulated process. Methods of pathway analysis provide a means to unravel the fundamental pathways and their controls in bioprocesses. In the present Ph.D. thesis, a novel methodology called Principal Elementary Mode Analysis (PEMA) was developed and implemented that enabled to identify which cellular fluxes are activated under different conditions of temperature and pH. It is shown that differences in these two parameters affect the chemical composition of EPS, hence they are critical for the regulation of the product synthesis. In future studies, the knowledge provided by PEMA could foster the development of metabolically meaningful control strategies that target the EPS sugar content and oder product quality parameters.

Performances de la puce exon et son application dans l’analyse de l’épissage alternatif associé à la métastase du cancer de sein

Nous montrons l’utilisation de la puce exon d’Affymetrix pour l’analyse simultanée de l’expression des gènes et de la variation d’isoformes. Nous avons utilisé les échantillons d’ARN du cerveau et des tissus de référence qui ont été antérieurement utilisés dans l’étude du consortium MicroArray Quality Control (MAQC). Nous démontrons une forte concordance de la quantification de l’expression des gènes entre trois plateformes d’expression populaires à savoir la puce exon d’Affymetrix, la puce Illumina et la puce U133A d’Affymetrix. Plus intéressant nous montrons que la majorité des discordances entre les trois plateformes résulterait des positions différentes des sondes à travers les plateformes et que les variations d’isoforme exactes ne peuvent être identifiées que par la puce exon. Nous avons détecté avec succès, entre les tissus de référence et ceux du cerveau, une centaine de cas d’évènements d’épissage alternatif. La puce exon est requise dans l’analyse de l’épissage alternatif associé aux pathologies telles que les cancers et les troubles neurologiques. Comme application de cette technologie, nous avons analysé les variations d’épissage dans la métastase du cancer de sein développé dans le model de la souris. Nous avons utilisé une gamme bien définie de trois lignées de tumeur mammaire ayant différents potentiels métastatiques. Par des analyses statistiques, nous avons répertorié 2623 transcripts présentant des variations d’expression et d’isoformes entre les types de tumeur. Une analyse du réseau de gènes montre qu’environ la moitié d’entre eux est impliquée dans plusieurs activités cellulaires, ainsi que dans nombreux cancers et désordres génétiques.

Genome wide search for genetic determinants of habitual alcohol, tobacco and coffee use, obesity-related traits, response to mental and physical stress and hemodynamic traits

Les habitudes de consommation de substances psychoactives, le stress, l’obésité et les traits cardiovasculaires associés seraient en partie reliés aux mêmes facteurs génétiques. Afin d’explorer cette hypothèse, nous avons effectué, chez 119 familles multi-générationnelles québécoises de la région du Saguenay-Lac-St-Jean, des études d’association et de liaison pangénomiques pour les composantes génétiques : de la consommation usuelle d’alcool, de tabac et de café, de la réponse au stress physique et psychologique, des traits anthropométriques reliés à l’obésité, ainsi que des mesures du rythme cardiaque (RC) et de la pression artérielle (PA). 58000 SNPs et 437 marqueurs microsatellites ont été utilisés et l’annotation fonctionnelle des gènes candidats identifiés a ensuite été réalisée. Nous avons détecté des corrélations phénotypiques significatives entre les substances psychoactives, le stress, l’obésité et les traits hémodynamiques. Par exemple, les consommateurs d’alcool et de tabac ont montré un RC significativement diminué en réponse au stress psychologique. De plus, les consommateurs de tabac avaient des PA plus basses que les non-consommateurs. Aussi, les hypertendus présentaient des RC et PA systoliques accrus en réponse au stress psychologique et un indice de masse corporelle (IMC) élevé, comparativement aux normotendus. D’autre part, l’utilisation de tabac augmenterait les taux corporels d’épinéphrine, et des niveaux élevés d’épinéphrine ont été associés à des IMC diminués. Ainsi, en accord avec les corrélations inter-phénotypiques, nous avons identifié plusieurs gènes associés/liés à la consommation de substances psychoactives, à la réponse au stress physique et psychologique, aux traits reliés à l’obésité et aux traits hémodynamiques incluant CAMK4, CNTN4, DLG2, DAG1, FHIT, GRID2, ITPR2, NOVA1, NRG3 et PRKCE. Ces gènes codent pour des protéines constituant un réseau d’interactions, impliquées dans la plasticité synaptique, et hautement exprimées dans le cerveau et ses tissus associés. De plus, l’analyse des sentiers de signalisation pour les gènes identifiés (P = 0,03) a révélé une induction de mécanismes de Potentialisation à Long Terme. Les variations des traits étudiés seraient en grande partie liées au sexe et au statut d’hypertension. Pour la consommation de tabac, nous avons noté que le degré et le sens des corrélations avec l’obésité, les traits hémodynamiques et le stress sont spécifiques au sexe et à la pression artérielle. Par exemple, si des variations ont été détectées entre les hommes fumeurs et non-fumeurs (anciens et jamais), aucune différence n’a été observée chez les femmes. Nous avons aussi identifié de nombreux traits reliés à l’obésité dont la corrélation avec la consommation de tabac apparaît essentiellement plus liée à des facteurs génétiques qu’au fait de fumer en lui-même. Pour le sexe et l’hypertension, des différences dans l’héritabilité de nombreux traits ont également été observées. En effet, des analyses génétiques sur des sous-groupes spécifiques ont révélé des gènes additionnels partageant des fonctions synaptiques : CAMK4, CNTN5, DNM3, KCNAB1 (spécifique à l’hypertension), CNTN4, DNM3, FHIT, ITPR1 and NRXN3 (spécifique au sexe). Ces gènes codent pour des protéines interagissant avec les protéines de gènes détectés dans l’analyse générale. De plus, pour les gènes des sous-groupes, les résultats des analyses des sentiers de signalisation et des profils d’expression des gènes ont montré des caractéristiques similaires à celles de l’analyse générale. La convergence substantielle entre les déterminants génétiques des substances psychoactives, du stress, de l’obésité et des traits hémodynamiques soutiennent la notion selon laquelle les variations génétiques des voies de plasticité synaptique constitueraient une interface commune avec les différences génétiques liées au sexe et à l’hypertension. Nous pensons, également, que la plasticité synaptique interviendrait dans de nombreux phénotypes complexes influencés par le mode de vie. En définitive, ces résultats indiquent que des approches basées sur des sous-groupes et des réseaux amélioreraient la compréhension de la nature polygénique des phénotypes complexes, et des processus moléculaires communs qui les définissent.