Evolutionary epigenetics of modern human populations

Epigenetic variability is a new mechanism for the study of human microevolution, because it creates both phenotypic diversity within an individual and within population. This mechanism constitutes an important reservoir for adaptation in response to new stimuli and recent studies have demonstrated that selective pressures shape not only the genetic code but also DNA methylation profiles. The aim of this thesis is the study of the role of DNA methylation changes in human adaptive processes, considering the Italian peninsula and macro-geographical areas. A whole-genome analysis of DNA methylation profile across the Italian penisula identified some genes whose methylation levels differ between individuals of different Italian districts (South, Centre and North of Italy). These genes are involved in nitrogen compound metabolism and genes involved in pathogens response. Considering individuals with different macro-geographical origins (individuals of Asians, European and African ancestry) more significant DMRs (differentially methylated regions) were identified and are located in genes involved in glucoronidation, in immune response as well as in cell comunication processes. A "profile" of each ancestry (African, Asian and European) was described. Moreover a deepen analysis of three candidate genes (KRTCAP3, MAD1L and BRSK2) in a cohort of individuals of different countries (Morocco, Nigeria, China and Philippines) living in Bologna, was performed in order to explore genetic and epigenetic diversity. Moreover this thesis have paved the way for the application of DNA methylation for the study of hystorical remains and in particular for the age-estimation of individuals starting from biological samples (such as teeth or blood). Noteworthy, a mathematical model that considered methylation values of DNA extracted from cementum and pulp of living individuals can estimate chronological age with high accuracy (median absolute difference between age estimated from DNA methylation and chronological age was 1.2 years).

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Studio del possibile ruolo del complesso polycomb nel determinare la schisi del labbro e del palato

La labioschisi con o senza palatoschisi non-sindromica (NSCL/P) è tra le più frequenti alterazioni dello sviluppo embrionale, causata dall’interazione di fattori genetici e ambientali, moti dei quali ancora ignoti. L'obiettivo del mio progetto di Dottorato consiste nell’identificazione di fattori di rischio genetico in un processo a due stadi che prevede la selezione di geni candidati e la verifica del loro coinvolgimento nella determinazione della malformazione mediante studi di associazione. Ho analizzato alcuni polimorfismi a singolo nucleotide (SNPs) dei geni RFC1 e DHFR, appartenenti alla via metabolica dell’acido folico, evidenziando una debole associazione tra alcuni degli SNPs indagati e la NSCL/P nella popolazione italiana. Presso il laboratorio della Dott.ssa Mangold dell’Università di Bonn, ho valutato il ruolo di 15 diverse regioni cromosomiche nel determinare la suscettibilità alla malattia, evidenziando una significativa associazione per i marcatori localizzati in 8q24 e 1p22. Ho quindi rivolto la mia attenzione al ruolo del complesso Polycomb nell’insorgenza della schisi. Nell’uomo i due complessi Polycomb, PRC1 e PRC2, rimodellano la cromatina agendo da regolatori dei meccanismi trascrizionali alla base della differenziazione cellulare e dello sviluppo embrionale. Ho ipotizzato che mutazioni a carico di geni appartenenti a PRC2 possano essere considerati potenziali fattori di rischio genetico nel determinare la NSCL/P. Il razionale consiste nel fatto che JARID2, una proteina che interagisce con PRC2, è associata all’insorgenza della NSCL/P ed espressa a livello delle cellule epiteliali delle lamine palatine che si approssimano alla fusione. L’indagine condotta analizzando i geni di elementi o partner dei due complessi Polycomb, ha evidenziato un’associazione significativa con alcuni polimorfismi dei geni indagati, associazione ulteriormente confermata dall’analisi degli aplotipi. Le analisi condotte sui geni candidati mi hanno permesso di raccogliere dati interessanti sull’eziologia della malformazione. Studi indipendenti saranno necessari per poter validare l'associazione tra le varianti genetiche di questi geni candidati e la NSCL/P.

Fine Mapping of qroot-yield-1.06, a QTL for Root, Plant Vigor and Yield in Maize

Root-yield-1.06 is a major QTL affecting root system architecture (RSA) and other agronomic traits in maize. The effect of this QTL has been evaluated with the development of near isogenic lines (NILs) differing at the QTL position. The objective of this study was to fine map qroot-yield-1.06 by marker-assisted searching for chromosome recombinants in the QTL interval and concurrent root phenotyping in both controlled and field conditions, through successive generations. Complementary approaches such as QTL meta-analysis and RNA-seq were deployed in order to help prioritizing candidate genes within the QTL target region. Using a selected group of genotypes, field based root analysis by ‘shovelomics’ enabled to accurately collect RSA information of adult maize plants. Shovelomics combined with software-assisted root imaging analysis proved to be an informative and relatively highly automated phenotyping protocol. A QTL interval mapping was conducted using a segregating population at the seedling stage grown in controlled environment. Results enabled to narrow down the QTL interval and to identify new polymorphic markers for MAS in field experiments. A collection of homozygous recombinant NILs was developed by screening segregating populations with markers flanking qroot-yield-1.06. A first set of lines from this collection was phenotyped based on the adapted shovelomics protocol. QTL analysis based on these data highlighted an interval of 1.3 Mb as completely linked with the target QTL but, a larger safer interval of 4.1 Mb was selected for further investigations. QTL meta-analysis allows to synthetize information on root QTLs and two mQTLs were identified in the qroot-yield-1.06 interval. Trascriptomics analysis based on RNA-seq data of the two contrasting QTL-NILs, confirmed alternative haplotypes at chromosome bin 1.06. qroot-yield-1.06 has now been delimited to a 4.1-Mb interval, and thanks to the availability of additional untested homozygous recombinant NILs, the potentially achievable mapping resolution at qroot-yield-1.06 is c. 50 kb.

Optimierung von Transformationsprotokollen und Entwicklung von SSR-Markern für das Kolumnarwachstum beim Apfel

Ziel war es, molekularbiologische Untersuchungen zum Kolumnarwachstum des Apfels durchzuführen. Anhand Sequenzdaten des ‘Golden Delicious’ Genoms (Velasco et al. 2010) wurden drei neue SSR Marker entwickelt. Sie konnten bei untersuchten Geisenheimer Nachkommenschaften zuverlässig den Kolumnarwuchs auf DNA-Ebene detektieren. Zusätzlich wurden von Bai et al. (2012) veröffentlichte Marker untersucht. Die von Bai et al. (2012) gefundenen Grenzen des co-Lokus konnten in dieser Arbeit anhand der Geisenheimer Nachkommenschaften nicht bestätigt werden. Die „linke“ Begrenzung der co-Region wird nach Untersuchungen dieser Arbeit am ehesten von dem Marker Mdo.chr10.11 (Moriya et al. 2012) bei 18,757 Mbp definiert. Die „rechte“ Begrenzung der co-Region wird vermutlich von den Markern Co04R13 (Baldi et al. 2012) und C1753-3520 (Bai et al. 2012) bei 18,905 Mbp definiert, wodurch die potentielle co-Region auf 148 kb auf Chromosom 10 eingegrenzt werden könnte. Für Funktionsanalysen möglicher Kandidatengene des co-Gens wurde ein Agrobakterien-vermitteltes Transformationssystem für die Geisenheimer Apfelselektionen ‘A 14’ und ‘Procats 28’ adaptiert. Zusätzlich wurde der bereits in der Literatur als transformierbar beschriebene Genotyp ‘Jonagold’ (Viss et al. 2003) transformiert. Bei Transformationen der Apfelselektion ‘A 14’ gelang es, transgene Zellen an den Explantaten, am Kallusgewebe und an den Regeneraten zu erzeugen. Bei Transformationen von ‘Jonagold’ wurde ein fast vollständig transgenes Regenerat erzeugt.

Regulatory T cell-mediated suppression of Th9 cell development and effector function

T helper (Th) 9 cells are an important subpopulation of the CD4+ T helper cells. Due to their ability to secrete Interleukin-(IL-)9, Th9 cells essentially contribute to the expulsion of parasitic helminths from the intestinal tract but they play also an immunopathological role in the course of asthma. Recently, a beneficial function of Th9 cells in anti-tumor immune responses was published. In a murine melanoma tumor model Th9 cells were shown to enhance the anti-melanoma immune response via the recruitment of CD8+ T cells, dendritic cells and mast cells. In contrast to Th9 effector cells regulatory T cells (Tregs) are able to control an immune response with the aid of different suppressive mechanisms. Based on their ability to suppress an immune response Tregs are believed to be beneficial in asthma by diminishing excessive allergic reactions. However, concerning cancer they can have a detrimental function because Tregs inhibit an effective anti-tumor immune reaction. Thus, the analysis of Th9 suppression by Tregs is of central importance concerning the development of therapeutic strategies for the treatment of cancer and allergic diseases and was therefore the main objective of this PhD thesis.rnIn general it could be demonstrated that the development of Th9 cells can be inhibited by Tregs in vitro. The production of the lineage-specific cytokine IL-9 by developing Th9 cells was completely suppressed at a Treg/Th9 ratio of 1:1 on the transcriptional (qRT-PCR) as well as on the translational level (ELISA). In contrast, the expression of IRF4 that was found to strongly promote Th9 development was not reduced in the presence of Tregs, suggesting that IRF4 requires additional transcription factors to induce the differentiation of Th9 cells. In order to identify such factors, which regulate Th9 development and therefore represent potential targets for Treg-mediated suppressive mechanisms, a transcriptome analysis using “next-generation sequencing” was performed. The expression of some genes which were found to be up- or downregulated in Th9 cells in the presence of Tregs was validated with qRT-PCR. Time limitations prevented a detailed functional analysis of these candidate genes. Nevertheless, the analysis of the suppressive mechanisms revealed that Tregs probably suppress Th9 cells via the increase of the intracellular cAMP concentration. In contrast, IL-9 production by differentiated Th9 cells was only marginally affected by Tregs in vitro and in vivo analysis (asthma, melanoma model). Hence, Tregs represent very effective inhibitors of Th9 development whereas they have only a minimal suppressive influence on differentiated Th9 cells.rn

Pharmacogenetics in palliative care

Response to analgesics, anticancer pharmacotherapy and pharmacotherapy of other cancer related symptoms vary broadly between individuals. Age, disease, comorbidities, concomitant medication, organ function and patients' compliance may partly explain the differences. However, the focus of ongoing research has shifted towards genomic variants of phase I and II drug metabolizing enzymes with one important goal being an individual dose adjustment according to a patient's genotype. Polymorphisms of the cytochrome P 450 2D6 influence the metabolism of many drugs including the analgesics codeine, tramadol, hydrocodone and oxycodone, as well as the metabolism of tricyclic antidepressants and the anticancer drug tamoxifen. Other candidate genes such as (opioid)-receptors, transporters and other molecules important for pharmacotherapy in pain management are discussed. Although pharmacogenetics as a diagnostic tool has the potential to improve patient therapy, study results are often equivocal and limited by small sample sizes and often by their retrospective design. Well designed studies are needed to demonstrate superiority of pharmoacogenetics to conventional dosing regimes.

Personalized therapy in pain management: where do we stand?

Genomic variations influencing response to pharmacotherapy of pain are currently under investigation. Drug-metabolizing enzymes represent a major target of ongoing research in order to identify associations between an individual's drug response and genetic profile. Polymorphisms of the cytochrome P450 enzymes (CYP2D6) influence metabolism of codeine, tramadol, hydrocodone, oxycodone and tricyclic antidepressants. Blood concentrations of some NSAIDs depend on CYP2C9 and/or CYP2C8 activity. Genomic variants of these genes associate well with NSAIDs' side effect profile. Other candidate genes, such as those encoding (opioid) receptors, transporters and other molecules important for pharmacotherapy in pain management, are discussed; however, study results are often equivocal. Besides genetic variants, further variables, for example, age, disease, comorbidity, concomitant medication, organ function as well as patients' compliance, may have an impact on pharmacotherapy and need to be addressed when pain therapists prescribe medication. Although pharmacogenetics as a diagnostic tool has the potential to improve patient therapy, well-designed studies are needed to demonstrate superiority to conventional dosing regimes.

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohn's disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10 ? ? ). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohn's disease.

Genome-wide association study in German patients with attention deficit/hyperactivity disorder

The heritability of attention deficit hyperactivity disorder (ADHD) is approximately 0.8. Despite several larger scale attempts, genome-wide association studies (GWAS) have not led to the identification of significant results. We performed a GWAS based on 495 German young patients with ADHD (according to DSM-IV criteria; Human660W-Quadv1; Illumina, San Diego, CA) and on 1,300 population-based adult controls (HumanHap550v3; Illumina). Some genes neighboring the single nucleotide polymorphisms (SNPs) with the lowest P-values (best P-value: 8.38 × 10(-7)) have potential relevance for ADHD (e.g., glutamate receptor, metabotropic 5 gene, GRM5). After quality control, the 30 independent SNPs with the lowest P-values (P-values ≤ 7.57 × 10(-5) ) were chosen for confirmation. Genotyping of these SNPs in up to 320 independent German families comprising at least one child with ADHD revealed directionally consistent effect-size point estimates for 19 (10 not consistent) of the SNPs. In silico analyses of the 30 SNPs in the largest meta-analysis so far (2,064 trios, 896 cases, and 2,455 controls) revealed directionally consistent effect-size point estimates for 16 SNPs (11 not consistent). None of the combined analyses revealed a genome-wide significant result. SNPs in previously described autosomal candidate genes did not show significantly lower P-values compared to SNPs within random sets of genes of the same size. We did not find genome-wide significant results in a GWAS of German children with ADHD compared to controls. The second best SNP is located in an intron of GRM5, a gene located within a recently described region with an infrequent copy number variation in patients with ADHD.

Determinants of sustained viral suppression in HIV-infected patients with self-reported poor adherence to antiretroviral therapy

Background Good adherence to antiretroviral therapy (ART) is critical for successful HIV treatment. However, some patients remain virologically suppressed despite suboptimal adherence. We hypothesized that this could result from host genetic factors influencing drug levels. Methods Eligible individuals were Caucasians treated with efavirenz (EFV) and/or boosted lopinavir (LPV/r) with self-reported poor adherence, defined as missing doses of ART at least weekly for more than 6 months. Participants were genotyped for single nucleotide polymorphisms (SNPs) in candidate genes previously reported to decrease EFV (rs3745274, rs35303484, rs35979566 in CYP2B6) and LPV/r clearance (rs4149056 in SLCO1B1, rs6945984 in CYP3A, rs717620 in ABCC2). Viral suppression was defined as having HIV-1 RNA <400 copies/ml throughout the study period. Results From January 2003 until May 2009, 37 individuals on EFV (28 suppressed and 9 not suppressed) and 69 on LPV/r (38 suppressed and 31 not suppressed) were eligible. The poor adherence period was a median of 32 weeks with 18.9% of EFV and 20.3% of LPV/r patients reporting missed doses on a daily basis. The tested SNPs were not determinant for viral suppression. Reporting missing >1 dose/week was associated with a lower probability of viral suppression compared to missing 1 dose/week (EFV: odds ratio (OR) 0.11, 95% confidence interval (CI): 0.01–0.99; LPV/r: OR 0.29, 95% CI: 0.09–0.94). In both groups, the probability of remaining suppressed increased with the duration of continuous suppression prior to the poor adherence period (EFV: OR 3.40, 95% CI: 0.62–18.75; LPV/r: OR 5.65, 95% CI: 1.82–17.56). Conclusions The investigated genetic variants did not play a significant role in the sustained viral suppression of individuals with suboptimal adherence. Risk of failure decreased with longer duration of viral suppression in this population.

Discovering imaging endophenotypes for major depression

Psychiatry research lacks an in-depth understanding of mood disorders phenotypes, leading to limited success of genetics studies of major depressive disorder (MDD). The dramatic progress in safe and affordable magnetic resonance-based imaging methods has the potential to identify subtle abnormalities of neural structures, connectivity and function in mood disordered subjects. This review paper presents strategies to improve the phenotypic definition of MDD by proposing imaging endophenotypes derived from magnetic resonance spectroscopy measures, such as cortical gamma-amino butyric acid (GABA) and glutamate/glutamine concentrations, and from measures of resting-state activity and functional connectivity. The proposed endophenotypes are discussed regarding specificity, mood state-independence, heritability, familiarity, clinical relevance and possible associations with candidate genes. By improving phenotypic definitions, the discovery of new imaging endophenotypes will increase the power of candidate gene and genome-wide associations studies. It will also help to develop and evaluate novel therapeutic treatments and enable clinicians to apply individually tailored therapeutic approaches. Finally, improvements of the phenotypic definition of MDD based on neuroimaging measures will contribute to a new classification system of mood disorders based on etiology and pathophysiology.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.

Influence of pregnancy on the adipocytokine and peroxisome proliferator-activated receptor pathways in peripheral blood mononuclear cells from healthy donors and rheumatoid arthritis patients

To identify candidate genes that are regulated by human pregnancy and have the potential to modulate rheumatoid arthritis (RA) disease activity.

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.