New computational biology tools for the systematic analysis of the structure and expression of human genes

From the late 1980s, the automation of sequencing techniques and the computer spread gave rise to a flourishing number of new molecular structures and sequences and to proliferation of new databases in which to store them. Here are presented three computational approaches able to analyse the massive amount of publicly avalilable data in order to answer to important biological questions. The first strategy studies the incorrect assignment of the first AUG codon in a messenger RNA (mRNA), due to the incomplete determination of its 5' end sequence. An extension of the mRNA 5' coding region was identified in 477 in human loci, out of all human known mRNAs analysed, using an automated expressed sequence tag (EST)-based approach. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1, QARS and TDP2 and the consequences for the functional studies are discussed. The second approach analyses the codon bias, the phenomenon in which distinct synonymous codons are used with different frequencies, and, following integration with a gene expression profile, estimates the total number of codons present across all the expressed mRNAs (named here "codonome value") in a given biological condition. Systematic analyses across different pathological and normal human tissues and multiple species shows a surprisingly tight correlation between the codon bias and the codonome bias. The third approach is useful to studies the expression of human autism spectrum disorder (ASD) implicated genes. ASD implicated genes sharing microRNA response elements (MREs) for the same microRNA are co-expressed in brain samples from healthy and ASD affected individuals. The different expression of a recently identified long non coding RNA which have four MREs for the same microRNA could disrupt the equilibrium in this network, but further analyses and experiments are needed.

Mechanistic study of the effects of 8-oxo-7,8-dihydroguanine (8-oxoG) on reporter gene expression in mammalian cells

DNA damage causes replication errors, leading to genetic instability or cell death. Besides that, many types of DNA base modifications have been shown to interfere with transcriptional elongation if they are located in the transcribed DNA strand of active genes, acting as roadblocks for RNA polymerases. It is widely assumed that transcription blockage by endogenous DNA damage is responsible for the early cell senescence in organs and accelerated ageing observed in individuals with compromised nucleotide excision repair.rnThe aims of this work were to design new experimental systems for testing transcription blocking potentials of DNA base modifications in an individual gene and to apply these test systems to the investigation of the effects of a frequent endogenously generated base modification, namely 8-oxo-7,8-hydroxyguanine (8-oxoG), on the gene transcription in cells. Several experimental strategies were employed for this purpose. First, I constructed an episomal vector encoding for a short-lived EGFP-ODC fusion protein and measured expression of the reporter gene in permanently transfected clonal cell lines exposed to DNA damaging agents. Second, the expression of plasmid-borne EGFP gene damaged with photosensitisers to obtain one or several oxidative purine modifications per plasmid molecule was determined in transiently transfected human and mouse host cells in an approach known as “host cell reactivation”. As a prerequisite for these experiments, a robust method of precise quantitative measurement of the EGFP gene expression in transiently transfected cells by flow cytometry was developed and validated. Third, I elaborated a very efficient procedure for insertion of synthetic oligonucleotides carrying 8-oxoG into plasmid DNA, avoiding any unwanted base damage and strand breaks. The consequences of 8-oxoG placed in defined positions in opposing DNA strands of the EGFP gene for transcription were measured by host cell reactivation in cells with functional 8-oxoguanine DNA glycosylase (OGG1) gene and in OGG1 null cells.rnThe results obtained in Ogg1-/- cells demonstrated that unrepaired 8-oxoG, even if situated in the transcribed DNA strand, does not have any negative effect on the reporter gene transcription. On the other hand, as few as one 8-oxoG was sufficient to cause a significant decrease of the gene expression in OGG1-proficient cell lines, i.e. in the presence of base excision repair. For two analysed positions of 8-oxoG in the plasmid DNA, the inhibition of gene transcription by the base modification correlated with the efficiency of its excision by purified OGG1 protein under cell-free conditions. Based on these findings, it has to be concluded that the observed decrease of transcription is mediated by excision of the base modification by OGG1 and probably caused by the repair-induced single-strand breaks. The mechanism of transcription inhibition by 8-oxoG is therefore clearly distinct from stalling of elongating RNA polymerase II complexes at the modified base.

Functional characterisation of in vitro synthesised G-protein coupled receptors in polymersomes

Membrane proteins play an indispensable role in physiological processes. It is, therefore, not surprising that many diseases are based on the malfunction of membrane proteins. Hence membrane proteins and especially G-protein coupled receptors(GPCRs)- the largest subfamily- have become an important drug target. Due to their high selectivity and sensitivity membrane proteins are also feasible for the detection of small quantities of substances with biosensors. Despite this widespread interest in GPCRs due to their importance as drug targets and biosensors there is still a lack of knowledge of structure, function and endogenous ligands for quiet a few of the previously identified receptors.rnBottlenecks in over-expression, purification, reconstitution and handling of membrane proteins arise due to their hydrophobic nature. Therefore the production of reasonable amounts of functional membrane proteins for structural and functional studies is still challenging. Also the limited stability of lipid based membrane systems hampers their application as platforms forrnscreening applications and biosensors.rnIn recent years the in vitro protein synthesis became a promising alternative to gain better yields for expression of membrane proteins in bio-mimetic membrane systems. These expression systems are based on cell extracts. Therefore cellular effects on protein expression are reduced. The open nature of the cell-free expression systems easily allows for the adjustment of reactionrnconditions for the protein of interest. The cell-free expression in the presence of bio-mimetic membrane systems allows the direct incorporation of the membrane proteins and therefore skips the time-consuming purification and reconstitution processes. Amphiphilic block-copolymers emerged as promising alternative for the less stable lipid-based membrane systems. They, likernlipids, form membraneous structures in aqueous solutions but exhibit increased mechanical and chemical stability.rnThe aim of this work was the generation of a GPCR-functionalised membrane system by combining both promising alternatives: in vitro synthesis and polymeric membrane systems. This novel platform should be feasible for the characterisation of the incorporated GPCR. Immunodetection of Dopamine receptor 1 and 2 expressed in diblock- and triblock-polymersomes demonstrated the successful in vitro expression of GPCRs in polymeric membranes. Antibodyrnbinding studies suggested a favoured orientation of dopamine receptors in triblockpolymersomes.rnA dopamine-replacement assay on DRD2-functionalised immobilised triblockpolymersomes confirmed functionality of the receptor in the polymersomes. The altered binding curve suggests an effect of the altered hydrophobic environment presented by the polymer membrane on protein activity.

Functional characterization of the placenta specific protein PLAC1 and its use for cancer immunotherapy

Summary Antibody-based cancer therapies have been successfully introduced into the clinic and have emerged as the most promising therapeutics in oncology. The limiting factor regarding the development of therapeutical antibody vaccines is the identification of tumor-associated antigens. PLAC1, the placenta-specific protein 1, was categorized for the first time by the group of Prof. Sahin as such a tumor-specific antigen. Within this work PLAC1 was characterized using a variety of biochemical methods. The protein expression profile, the cellular localization, the conformational state and especially the interacting partners of PLAC1 and its functionality in cancer were analyzed. Analysis of the protein expression profile of PLAC1 in normal human tissue confirms the published RT-PCR data. Except for placenta no PLAC1 expression was detectable in any other normal human tissue. Beyond, an increased PLAC1 expression was detected in several cancer cell lines derived of trophoblastic, breast and pancreatic lineage emphasizing its properties as tumor-specific antigen. rnThe cellular localization of PLAC1 revealed that PLAC1 contains a functional signal peptide which conducts the propeptide to the endoplasmic reticulum (ER) and results in the secretion of PLAC1 by the secretory pathway. Although PLAC1 did not exhibit a distinct transmembrane domain, no unbound protein was detectable in the cell culture supernatant of overexpressing cells. But by selective isolation of different cellular compartments PLAC1 was clearly enriched within the membrane fraction. Using size exclusion chromatography PLAC1 was characterized as a highly aggregating protein that forms a network of high molecular multimers, consisting of a mixture of non-covalent as well as covalent interactions. Those interactions were formed by PLAC1 with itself and probably other cellular components and proteins. Consequently, PLAC1 localize outside the cell, where it is associated to the membrane forming a stable extracellular coat-like structure.rnThe first mechanistic hint how PLAC1 promote cancer cell proliferation was achieved identifying the fibroblast growth factor FGF7 as a specific interacting partner of PLAC1. Moreover, it was clearly shown that PLAC1 as well as FGF7 bind to heparin, a glycosaminoglycan of the ECM that is also involved in FGF-signaling. The participation of PLAC1 within this pathway was approved after co-localizing PLAC1, FGF7 and the FGF7 specific receptor (FGFR2IIIb) and identifying the formation of a trimeric complex (PLAC1, FGF7 and the specific receptor FGFR2IIIb). Especially this trimeric complex revealed the role of PLAC1. Binding of PLAC1 together with FGF7 leads to the activation of the intracellular tyrosine kinase of the FGFR2IIIb-receptor and mediate the direct phosphorylation of the AKT-kinase. In the absence of PLAC1, no FGF7 mediated phosphorylation of AKT was observed. Consequently the function of PLAC1 was clarified: PLAC1 acts as a co-factor by stimulating proliferation by of the FGF7-FGFR2 signaling pathway.rnAll together, these novel biochemical findings underline that the placenta specific protein PLAC1 could be a new target for cancer immunotherapy, especially considering its potential applicability for antibody therapy in tumor patients.

Preparation and functional characterisation of recombinant silicatein-alpha from the sponge Suberites domuncula

The marine world is an immense source of biodiversity that provides substances with striking potentials in medicinal chemistry and biotechnology. Sponges (Porifera) are marine animals that represent the most impressive example of organisms possessing the ability to metabolise silica through a family of enzymes known as silicateins. Complex skeletal structures (spicules) made of pure biogenic silica (biosilica) are produced under physiological conditions. Biosilica is a natural material comprising inorganic and organic components with unique mechanical, optical, and physico-chemical properties, including promising potential to be used for development of therapeutic agents in regenerative medicine. Unravelling the intimate physiological mechanisms occurring in sponges during the construction of their siliceous spicules is an on-going project, and several questions have been addressed by the studies proposed by our working group. In this doctoral work, the recombinant DNA technology is exploited for functional and structural characterisation of silicatein. Its precursors are produced as fusion proteins with a chaperone tag (named TF-Ps), and a robust method for the overexpression of native soluble proteins in high concentrations has been developed. In addition, it is observed and proven experimentally that the maturation of silicatein is an autocatalytic event that: (i) can be modulated by rational use of protease inhibitors; (ii) is influenced by the temperature of the environment; (iii) only slightly depends on the pH. In the same experimental framework, observations on the dynamics in the maturation of silicateins allow a better understanding of how the axial filaments form during the early stages of spicule construction. In addition, the definition of new distinct properties of silicatein (termed “structure-guiding” and “structure-forming”) is introduced. By homology models and through comparisons with similar proteins (the cathepsins), domains with significant surface hydrophobicity are identified as potential self-assembly mediators. Moreover, a high-throughput screening showed that TF-Ps could generate crystals under certain conditions, becoming promising for further structural studies. With the goal of optimise the properties of the recombinant silicatein, implementation of new production systems are tried for the first time. Success in the expression of silicatein-type proteins in insect and yeast cells, constitute a promising basis for further development, towards the establishment of an efficient method for the production of a high-value pure and soluble protein.

Comparative functional analysis of factors controlling glial differentiation in Drosophila and mouse

The present study is a comparative functional analysis of three factors controlling glial differentiation in mouse (Fyn Src kinase, hnRNPF/H and NG2) and their homologues in Drosophila (Src42A and 64B, Glorund and Kon-tiki (Kon)). In Drosophila, mutations in any of these genes were not associated with major embryonic neurodevelopmental phenotypes. Src kinases and Glorund were shown to be ubiquitously expressed, whereas kon mRNA showed selective expression in muscles as well as in central and peripheral glia. Kon was also shown to be expressed in L3 larvae with high levels of protein accumulation at the neuromuscular junction (NMJ) and in muscles in the form of speckles. Knockdown of kon in glia resulted in NMJ phenotypes, mainly characterized by a significant increase in bouton number and a reduction in α-Konecto staining intensity at the NMJ. From the three glial layers ensheathing the peripheral nervous system, subperineurial glial showed to be the one contributing the most to kon knockdown dependent NMJ phenotypes, while perineurial glia only had a minor role. The knockdown of kon in glia also showed to affect Glutamate receptor subunit (α-GluRIIA) clustering in the postsynapse, same as microtubule arrangement in the presynapse, as seen by α-Futsch pattern interruptions and alterations. kon knockdown in glia also resulted in impaired axonal transport, as seen by the accumulation of Bruchpilot-positive vesicles along the nerves, abnormal formation of neuronal derived protrusions and swellings, filled with vacuole-like structures. Glia number along the peripheral nerves is also reduced as consequence of kon knockdown. Muscle derived Kon was shown to accumulate at the NMJ and play a role in bouton consolidation and to interfere with phagocytosis of ghost boutons. NMJ bouton and branch number was also significantly increased in Kon overexpression in glia. The overexpression of Kon in glia also resulted in a massive elongation of the ventral nerve cord, which served in a suppressor screen to identify intracellular interaction partners of Kon in glia. It was shown that Kon is processed in glia and preliminary results indicate that the metalloendopeptidase Kuzbanian (the fly homologue of ADAM10) may play a role in the shedding of Konecto. In the present work, Kon is shown as a multifunctional gene with various roles in glia-neuron and glia-neuron-muscle interaction.

Pheochromocytoma in rats with multiple endocrine neoplasia (MENX) shares gene expression patterns with human pheochromocytoma

Pheochromocytomas are rare neoplasias of neural crest origin arising from chromaffin cells of the adrenal medulla and sympathetic ganglia (extra-adrenal pheochromocytoma). Pheochromocytoma that develop in rats homozygous for a loss-of-function mutation in p27Kip1 (MENX syndrome) show a clear progression from hyperplasia to tumor, offering the possibility to gain insight into tumor pathobiology. We compared the gene-expression signatures of both adrenomedullary hyperplasia and pheochromocytoma with normal rat adrenal medulla. Hyperplasia and tumor show very similar transcriptome profiles, indicating early determination of the tumorigenic signature. Overrepresentation of developmentally regulated neural genes was a feature of the rat lesions. Quantitative RT-PCR validated the up-regulation of 11 genes, including some involved in neural development: Cdkn2a, Cdkn2c, Neurod1, Gal, Bmp7, and Phox2a. Overexpression of these genes precedes histological changes in affected adrenal glands. Their presence at early stages of tumorigenesis indicates they are not acquired during progression and may be a result of the lack of functional p27Kip1. Adrenal and extra-adrenal pheochromocytoma development clearly follows diverged molecular pathways in MENX rats. To correlate these findings to human pheochromocytoma, we studied nine genes overexpressed in the rat lesions in 46 sporadic and familial human pheochromocytomas. The expression of GAL, DGKH, BMP7, PHOX2A, L1CAM, TCTE1, EBF3, SOX4, and HASH1 was up-regulated, although with different frequencies. Immunohistochemical staining detected high L1CAM expression selectively in 27 human pheochromocytomas but not in 140 nonchromaffin neuroendocrine tumors. These studies reveal clues to the molecular pathways involved in rat and human pheochromocytoma and identify previously unexplored biomarkers for clinical use.

Expression of human cathepsin L or human cathepsin V in mouse thymus mediates positive selection of T helper cells in cathepsin L knock-out mice

A genetic deficiency of the cysteine protease cathepsin L (Ctsl) in mice results in impaired positive selection of conventional CD4+ T helper cells as a result of an incomplete processing of the MHC class II associated invariant chain or incomplete proteolytic generation of positively selecting peptide ligands. The human genome encodes, in contrast to the mouse genome, for two cathepsin L proteases, namely cathepsin L (CTSL) and cathepsin V (CTSV; alternatively cathepsin L2). In the human thymic cortex, CTSV is the predominately expressed protease as compared to CTSL or other cysteine cathepsins. In order to analyze the functions of CTSL and CTSV in the positive selection of CD4+ T cells we employed Ctsl knock-out mice crossed either with transgenic mice expressing CTSL under the control of its genuine human promoter or with transgenic mice expressing CTSV under the control of the keratin 14 (K14) promoter, which drives expression to the cortical epithelium. Both human proteases are expressed in the thymus of the transgenic mice, and independent expression of both CTSL and CTSV rescues the reduced frequency of CD4+ T cells in Ctsl-deficient mice. Moreover, the expression of the human cathepsins does not change the number of CD4+CD25+Foxp3+ regulatory T cells, but the normalization of the frequency of conventional CD4+ T cell in the transgenic mice results in a rebalancing of conventional T cells and regulatory T cells. We conclude that the functional differences of CTSL and CTSV in vivo are not mainly determined by their inherent biochemical properties, but rather by their tissue specific expression pattern.

Drug antigenicity, immunogenicity, and costimulatory signaling: evidence for formation of a functional antigen through immune cell metabolism

Recognition of drugs by immune cells is usually explained by the hapten model, which states that endogenous metabolites bind irreversibly to protein to stimulate immune cells. Synthetic metabolites interact directly with protein-generating antigenic determinants for T cells; however, experimental evidence relating intracellular metabolism in immune cells and the generation of physiologically relevant Ags to functional immune responses is lacking. The aim of this study was to develop an integrated approach using animal and human experimental systems to characterize sulfamethoxazole (SMX) metabolism-derived antigenic protein adduct formation in immune cells and define the relationship among adduct formation, cell death, costimulatory signaling, and stimulation of a T cell response. Formation of SMX-derived adducts in APCs was dose and time dependent, detectable at nontoxic concentrations, and dependent on drug-metabolizing enzyme activity. Adduct formation above a threshold induced necrotic cell death, dendritic cell costimulatory molecule expression, and cytokine secretion. APCs cultured with SMX for 16 h, the time needed for drug metabolism, stimulated T cells from sensitized mice and lymphocytes and T cell clones from allergic patients. Enzyme inhibition decreased SMX-derived protein adduct formation and the T cell response. Dendritic cells cultured with SMX and adoptively transferred to recipient mice initiated an immune response; however, T cells were stimulated with adducts derived from SMX metabolism in APCs, not the parent drug. This study shows that APCs metabolize SMX; subsequent protein binding generates a functional T cell Ag. Adduct formation above a threshold stimulates cell death, which provides a maturation signal for dendritic cells.

Vaccination-induced functional competence of circulating human tumor-specific CD8 T-cells

T-cells specific for foreign (e.g., viral) antigens can give rise to strong protective immune responses, whereas self/tumor antigen-specific T-cells are thought to be less powerful. However, synthetic T-cell vaccines composed of Melan-A/MART-1 peptide, CpG and IFA can induce high frequencies of tumor-specific CD8 T-cells in PBMC of melanoma patients. Here we analyzed the functionality of these T-cells directly ex vivo, by multiparameter flow cytometry. The production of multiple cytokines (IFNγ, TNFα, IL-2) and upregulation of LAMP-1 (CD107a) by tumor (Melan-A/MART-1) specific T-cells was comparable to virus (EBV-BMLF1) specific CD8 T-cells. Furthermore, phosphorylation of STAT1, STAT5 and ERK1/2, and expression of CD3 zeta chain were similar in tumor- and virus-specific T-cells, demonstrating functional signaling pathways. Interestingly, high frequencies of functionally competent T-cells were induced irrespective of patient's age or gender. Finally, CD8 T-cell function correlated with disease-free survival. However, this result is preliminary since the study was a Phase I clinical trial. We conclude that human tumor-specific CD8 T-cells can reach functional competence in vivo, encouraging further development and Phase III trials assessing the clinical efficacy of robust vaccination strategies.

High expression of NPY receptors in the human testis

NPY receptors represent novel molecular therapeutic targets in cancer and obesity. However, the extent of NPY receptor expression in normal human tissues is poorly investigated. Based on the role of NPY in reproductive functions, the NPY receptor expression was studied in 25 normal human testes and, additionally, 24 testicular tumors using NPY receptor autoradiography. In the normal testis, Leydig cells strongly expressed NPY receptor subtype Y2, and small arterial blood vessels Y1. Y2 receptors were found to be functional with agonist-stimulated [(35)S]GTPγS binding autoradiography. Full functional integrity of the NPY system was further suggested by the immunohistochemical detection of NPY peptide in nerve fibers directly adjacent to Leydig cells and arteries. Germ cell tumors expressed Y1 and Y2 on tumor cells in 33% and Y1 on intratumoral blood vessels in 50%. Based on its strong NPY receptor expression in Leydig cells and blood vessels, the normal human testis represents a potentially important physiological and pharmalogical NPY target.

Functional Imaging of Pheochromocytoma with Ga-DOTATOC and C-HED in a Genetically Defined Rat Model of Multiple Endocrine Neoplasia

Rats affected by the MENX multitumor syndrome develop pheochromocytoma (100%). Pheochromocytomas are uncommon tumors and animal models are scarce, hence the interest in MENX rats to identify and preclinically evaluate novel targeted therapies. A prerequisite for such studies is a sensitive and noninvasive detection of MENXassociated pheochromocytoma. We performed positron emission tomography (PET) to determine whether rat pheochromocytomas are detected by tracers used in clinical practice, such as 68Ga-DOTATOC (somatostatin analogue) or (11)C-Hydroxyephedrine (HED), a norepinephrine analogue. We analyzed four affected and three unaffected rats. The PET scan findings were correlated to histopathology and immunophenotype of the tumors, their proliferative index, and the expression of genes coding for somatostatin receptors or the norepinephrine transporter. We observed that mean 68Ga-DOTATOC standard uptake value (SUV) in adrenals of affected animals was 23.3 ± 3.9, significantly higher than in control rats (15.4 ± 7.9; P = .03). The increase in mean tumor-to-liver ratio of (11)C-HED in the MENX-affected animals (1.6 ± 0.5) compared to controls (0.7 ± 0.1) was even more significant (P = .0016). In a unique animal model, functional imaging depicting two pathways important in pheochromocytoma biology discriminated affected animals from controls, thus providing the basis for future preclinical work with MENX rats.

Hormonally defined pancreatic and duodenal neuroendocrine tumors differ in their transcription factor signatures: expression of ISL1, PDX1, NGN3, and CDX2

We recently identified the transcription factor (TF) islet 1 gene product (ISL1) as a marker for well-differentiated pancreatic neuroendocrine tumors (P-NETs). In order to better understand the expression of the four TFs, ISL1, pancreatico-duodenal homeobox 1 gene product (PDX1), neurogenin 3 gene product (NGN3), and CDX-2 homeobox gene product (CDX2), that mainly govern the development and differentiation of the pancreas and duodenum, we studied their expression in hormonally defined P-NETs and duodenal (D-) NETs. Thirty-six P-NETs and 14 D-NETs were immunostained with antibodies against the four pancreatic hormones, gastrin, serotonin, calcitonin, ISL1, PDX1, NGN3, and CDX2. The TF expression pattern of each case was correlated with the tumor's hormonal profile. Insulin-positive NETs expressed only ISL1 (10/10) and PDX1 (9/10). Glucagon-positive tumors expressed ISL1 (7/7) and were almost negative for the other TFs. Gastrin-positive NETs, whether of duodenal or pancreatic origin, frequently expressed PDX1 (17/18), ISL1 (14/18), and NGN3 (14/18). CDX2 was mainly found in the gastrin-positive P-NETs (5/8) and rarely in the D-NETs (1/10). Somatostatin-positive NETs, whether duodenal or pancreatic in origin, expressed ISL1 (9/9), PDX1 (3/9), and NGN3 (3/9). The remaining tumors showed labeling for ISL1 in addition to NGN3. There was no association between a particular TF pattern and NET features such as grade, size, location, presence of metastases, and functional activity. We conclude from our data that there is a correlation between TF expression patterns and certain hormonally defined P-NET and D-NET types, suggesting that most of the tumor types originate from embryologically determined precursor cells. The observed TF signatures do not allow us to distinguish P-NETs from D-NETs.

miR-21 functionally interacts with the 3'UTR of chemokine CCL20 and down-regulates CCL20 expression in miR-21 transfected colorectal cancer cells

As deregulation of miRNAs and chemokine CCL20 was shown to play a role in colorectal cancer (CRC) pathogenesis, we analyzed the functional interactions of candidate miRNAs with CCL20 mRNA. After target prediction software programs indicated a role for miR-21 in CCL20 regulation, we applied the luciferase reporter assay system to demonstrate that miR-21 functionally interacts with the 3'UTR of CCL20 mRNA and down-regulates CCL20 in miR-21 mimic transfected CRC cell lines (Caco-2, SW480 and SW620). Thus, regulation of CCL20 expression by miR-21 might be a regulatory mechanism involved in progression of CRC.

Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score

BACKGROUND: During the past ten years many quantitative trait loci (QTL) affecting mastitis incidence and mastitis related traits like somatic cell score (SCS) were identified in cattle. However, little is known about the molecular architecture of QTL affecting mastitis susceptibility and the underlying physiological mechanisms and genes causing mastitis susceptibility. Here, a genome-wide expression analysis was conducted to analyze molecular mechanisms of mastitis susceptibility that are affected by a specific QTL for SCS on Bos taurus autosome 18 (BTA18). Thereby, some first insights were sought into the genetically determined mechanisms of mammary gland epithelial cells influencing the course of infection. METHODS: Primary bovine mammary gland epithelial cells (pbMEC) were sampled from the udder parenchyma of cows selected for high and low mastitis susceptibility by applying a marker-assisted selection strategy considering QTL and molecular marker information of a confirmed QTL for SCS in the telomeric region of BTA18. The cells were cultured and subsequently inoculated with heat-inactivated mastitis pathogens Escherichia coli and Staphylococcus aureus, respectively. After 1, 6 and 24 h, the cells were harvested and analyzed using the microarray expression chip technology to identify differences in mRNA expression profiles attributed to genetic predisposition, inoculation and cell culture. RESULTS: Comparative analysis of co-expression profiles clearly showed a faster and stronger response after pathogen challenge in pbMEC from less susceptible animals that inherited the favorable QTL allele 'Q' than in pbMEC from more susceptible animals that inherited the unfavorable QTL allele 'q'. Furthermore, the results highlighted RELB as a functional and positional candidate gene and related non-canonical Nf-kappaB signaling as a functional mechanism affected by the QTL. However, in both groups, inoculation resulted in up-regulation of genes associated with the Ingenuity pathways 'dendritic cell maturation' and 'acute phase response signaling', whereas cell culture affected biological processes involved in 'cellular development'. CONCLUSIONS: The results indicate that the complex expression profiling of pathogen challenged pbMEC sampled from cows inheriting alternative QTL alleles is suitable to study genetically determined molecular mechanisms of mastitis susceptibility in mammary epithelial cells in vitro and to highlight the most likely functional pathways and candidate genes underlying the QTL effect.