Genome Wide Association Mapping of Grain Arsenic, Copper, Molybdenum and Zinc in Rice (Oryza sativa L.) Grown at Four International Field Sites

The mineral concentrations in cereals are important for human health, especially for individuals who consume a cereal subsistence diet. A number of elements, such as zinc, are required within the diet, while some elements are toxic to humans, for example arsenic. In this study we carry out genome-wide association (GWA) mapping of grain concentrations of arsenic, copper, molybdenum and zinc in brown rice using an established rice diversity panel of,300 accessions and 36.9 k single nucleotide polymorphisms (SNPs). The study was performed across five environments: one field site in Bangladesh, one in China and two in the US, with one of the US sites repeated over two years. GWA mapping on the whole dataset and on separate subpopulations of rice revealed a large number of loci significantly associated with variation in grain arsenic, copper, molybdenum and zinc. Seventeen of these loci were detected in data obtained from grain cultivated in more than one field location, and six co-localise with previously identified quantitative trait loci. Additionally, a number of candidate genes for the uptake or transport of these elements were located near significantly associated SNPs (within 200 kb, the estimated global linkage disequilibrium previously employed in this rice panel). This analysis highlights a number of genomic regions and candidate genes for further analysis as well as the challenges faced when mapping environmentally-variable traits in a highly genetically structured diversity panel.

Genetic mapping of the rice ionome in leaves and grain:Identification of QTLs for 17 elements including arsenic, cadmium, iron and selenium

Research into the composition of cereal grains is motivated by increased interest in food quality. Here multi-element analysis is conducted on leaves and grain of the Bala x Azucena rice mapping population grown in the field. Quantitative trait loci (QTLs) for the concentration of 17 elements were detected, revealing 36 QTLs for leaves and 41 for grains. Epistasis was detected for most elements. There was very little correlation between leaf and grain element concentrations. For selenium, lead, phosphorus and magnesium QTLs were detected in the same location for both tissues. In general, there were no major QTL clusters, suggesting separate regulation of each element. QTLs for grain iron, zinc, molybdenum and selenium are potential targets for marker assisted selection to improve seed nutritional quality. An epistatic interaction for grain arsenic also looks promising to decrease the concentration of this carcinogenic element. © Springer Science+Business Media B.V. 2009.

Integration of copy number and transcriptomics provides risk stratification in prostate cancer: A discovery and validation cohort study

BACKGROUND: Understanding the heterogeneous genotypes and phenotypes of prostate cancer is fundamental to improving the way we treat this disease. As yet, there are no validated descriptions of prostate cancer subgroups derived from integrated genomics linked with clinical outcome.

METHODS: In a study of 482 tumour, benign and germline samples from 259 men with primary prostate cancer, we used integrative analysis of copy number alterations (CNA) and array transcriptomics to identify genomic loci that affect expression levels of mRNA in an expression quantitative trait loci (eQTL) approach, to stratify patients into subgroups that we then associated with future clinical behaviour, and compared with either CNA or transcriptomics alone.

FINDINGS: We identified five separate patient subgroups with distinct genomic alterations and expression profiles based on 100 discriminating genes in our separate discovery and validation sets of 125 and 103 men. These subgroups were able to consistently predict biochemical relapse (p = 0.0017 and p = 0.016 respectively) and were further validated in a third cohort with long-term follow-up (p = 0.027). We show the relative contributions of gene expression and copy number data on phenotype, and demonstrate the improved power gained from integrative analyses. We confirm alterations in six genes previously associated with prostate cancer (MAP3K7, MELK, RCBTB2, ELAC2, TPD52, ZBTB4), and also identify 94 genes not previously linked to prostate cancer progression that would not have been detected using either transcript or copy number data alone. We confirm a number of previously published molecular changes associated with high risk disease, including MYC amplification, and NKX3-1, RB1 and PTEN deletions, as well as over-expression of PCA3 and AMACR, and loss of MSMB in tumour tissue. A subset of the 100 genes outperforms established clinical predictors of poor prognosis (PSA, Gleason score), as well as previously published gene signatures (p = 0.0001). We further show how our molecular profiles can be used for the early detection of aggressive cases in a clinical setting, and inform treatment decisions.

INTERPRETATION: For the first time in prostate cancer this study demonstrates the importance of integrated genomic analyses incorporating both benign and tumour tissue data in identifying molecular alterations leading to the generation of robust gene sets that are predictive of clinical outcome in independent patient cohorts.

Gene regulatory mechanisms underpinning prostate cancer susceptibility

Molecular characterization of genome-wide association study (GWAS) loci can uncover key genes and biological mechanisms underpinning complex traits and diseases. Here we present deep, high-throughput characterization of gene regulatory mechanisms underlying prostate cancer risk loci. Our methodology integrates data from 295 prostate cancer chromatin immunoprecipitation and sequencing experiments with genotype and gene expression data from 602 prostate tumor samples. The analysis identifies new gene regulatory mechanisms affected by risk locus SNPs, including widespread disruption of ternary androgen receptor (AR)-FOXA1 and AR-HOXB13 complexes and competitive binding mechanisms. We identify 57 expression quantitative trait loci at 35 risk loci, which we validate through analysis of allele-specific expression. We further validate predicted regulatory SNPs and target genes in prostate cancer cell line models. Finally, our integrated analysis can be accessed through an interactive visualization tool. This analysis elucidates how genome sequence variation affects disease predisposition via gene regulatory mechanisms and identifies relevant genes for downstream biomarker and drug development.

An integrated and comparative approach to genetic selection of the aquaculture species Dicentrarchus labrax

The European sea bass, Dicentrarchus labrax, is one of the most important marine species cultivated in Southern Europe and has not benefited from selective breeding. One of the major goals in the sea bass (D. labrax) aquaculture industry is to understand and control the complexity of growth associated traits. The aim of the methodology developed for the studies reported in the thesis was not only to establish genetic and genomic resources for sea bass, but to also develop a conceptual strategy to efficiently create knowledge in a research environment that can easily be transferred to the aquaculture industry. The strategy involved; i) establishing an annotated sea bass transcriptome and then using it to, ii) identify new genetic markers for target QTL regions so that, iii) new QTL analysis could be performed and marker based resolution of the DNA regions of interest increased, and then iv) to merge the linkage map and the physical map in order to map the QTL confidence intervals to the sea bass genome and identify genes underlying the targeted traits. Finally to test if genes in the QTL regions that are candidates for divergent growth phenotypes have modified patterns of transcription that reflects the modified whole organism physiology SuperSAGE-SOLiD4 gene expression was used with sea bass with high growth heterogeneity. The SuperSAGE contributed to significantly increase the transcriptome information for sea bass muscle, brain and liver and also led to the identification of putative candidate genes lying in the genomic region of growth related QTL. Lastly all differentially expressed transcripts in brain, liver and muscle of the European sea bass with divergent specific growth rates were mapped to gene pathways and networks and the regulatory pathways most affected identified and established the tissue specific changes underlying the divergent SGR. Owing to the importance of European sea bass to Mediterranean aquaculture and the developed genomics resources from the present thesis and from other studies it should be possible to implement genetic selection programs using marker assisted selection.

Variants in MTNR1B influence fasting glucose levels.

To identify previously unknown genetic loci associated with fasting glucose concentrations, we examined the leading association signals in ten genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95% CI = 0.06-0.08) mmol/l in fasting glucose levels (P = 3.2 x 10(-50)) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P = 1.1 x 10(-15)). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05-1.12), per G allele P = 3.3 x 10(-7)) in a meta-analysis of 13 case-control studies totaling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P = 1.1 x 10(-57)) and GCK (rs4607517, P = 1.0 x 10(-25)) loci.

Study of a Swiss dopa-responsive dystonia family with a deletion in GCH1: redefining DYT14 as DYT5.

OBJECTIVE: To report the study of a multigenerational Swiss family with dopa-responsive dystonia (DRD). METHODS: Clinical investigation was made of available family members, including historical and chart reviews. Subject examinations were video recorded. Genetic analysis included a genome-wide linkage study with microsatellite markers (STR), GTP cyclohydrolase I (GCH1) gene sequencing, and dosage analysis. RESULTS: We evaluated 32 individuals, of whom 6 were clinically diagnosed with DRD, with childhood-onset progressive foot dystonia, later generalizing, followed by parkinsonism in the two older patients. The response to levodopa was very good. Two additional patients had late onset dopa-responsive parkinsonism. Three other subjects had DRD symptoms on historical grounds. We found suggestive linkage to the previously reported DYT14 locus, which excluded GCH1. However, further study with more stringent criteria for disease status attribution showed linkage to a larger region, which included GCH1. No mutation was found in GCH1 by gene sequencing but dosage methods identified a novel heterozygous deletion of exons 3 to 6 of GCH1. The mutation was found in seven subjects. One of the patients with dystonia represented a phenocopy. CONCLUSIONS: This study rules out the previously reported DYT14 locus as a cause of disease, as a novel multiexonic deletion was identified in GCH1. This work highlights the necessity of an accurate clinical diagnosis in linkage studies as well as the need for appropriate allele frequencies, penetrance, and phenocopy estimates. Comprehensive sequencing and dosage analysis of known genes is recommended prior to genome-wide linkage analysis.

Common variants near MC4R are associated with fat mass, weight and risk of obesity.

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Évaluation neurobiologique des souris spontanément hypertendues : Du vieillissement à la génomique comparative

Le but de cette thèse est premièrement d’évaluer l’effet du vieillissement sur les fonctions psychomotrices des souches de souris sélectionnées génétiquement en fonction de leur tension artérielle (TA); deuxièmement, de localiser les déterminants génétiques des phénotypes psychophysiologiques à partir de souches recombinantes congéniques (RCS). Ces travaux ont mené à la publication de 4 articles. Le premier article décrit l’évaluation des fonctions psychomotrices des souches avec une tension artérielle élevée (HBP), basse (LBP) et normale (NBP). La performance aux épreuves d’exploration, d’habiletés motrices et d’apprentissage spatial, a été mesurée sur deux cohortes âgées respectivement de 12 mois et de trois mois. Indépendamment de l’âge, les HBPs sont hyperactives dans l’open-field (OF), mais pas dans le test d’exploration de trous. Inversement, les LBP explorent moins d’espaces que les NBP et, à trois mois seulement, sont hypoactives dans l’OF. Par ailleurs, les HBPs et les LBP présentent des déficits précoces de coordination motrice et des fonctions visuo-motrices. Le second article concerne l’évaluation longitudinale de la coordination motrice, de l’anxiété et de l’apprentissage spatial des souches HBP, LBP et NBP, à l’âge de deux mois et de 12 mois. Le vieillissement accentue l’hyperactivité des HBPs dans l’OF. Par contre, l’hypoactivité des souris LBP est détectable seulement à l’âge de deux mois. Indépendamment de l’âge, les souris HBP et LBP montrent une perception réduite du danger dans l’épreuve d’anxiété et des dysfonctions visuo-motrices au labyrinthe aquatique. Enfin, des déficits précoces de coordination motrice se manifestent seulement chez les HBPs. Il reste à déterminer si les déficits observés sont liés à des déterminants génétiques indépendants ou secondaires aux altérations de la tension artérielle. Le troisième article présente la comparaison entre les souches consanguines A/J et C57Bl/6J (B6) aux épreuves de l’OF, de la planche à trous, du labyrinthe aquatique et du cintre (coordination motrice). Les B6 explore d’avantage l’OF et la planche à trous. Les B6 sont moins rapides sur le cintre, mais supérieurs aux A/J dans le labyrinthe aquatique, avec une plate-forme invisible ou visible. Ces résultats démontrent l’implication de déterminants génétiques. Cette thèse se termine par un quatrième article sur la localisation des déterminants génétiques de la susceptibilité au stress dans les RCS, dérivées de A/J et B6, et présentant un agencement spécifique de 12.5% du génome. La réactivité émotionnelle est évaluée dans l’OF et le plus-maze; la réponse de stress est mesurée par radio télémétrie de la température interne pendant le stress d’immobilisation (SI) sous diète régulière et riche en sel; l’excrétion des électrolytes urinaires est dosée après 24 heures de diète salée. Les loci les plus significatifs sont situés dans les régions suivantes: de l’émotionalité dans l’OF (Emo1) sur le chr. 1 (LOD=4.6) correspondant à la région homologue impliquée dans la cohorte d’hypertension familiale du Saguenay; de la dopa décarboxylase (ddc) sur le chr. 11 pour l’émergence du plus-maze (LOD=4.7); de la protéine liant l’endotoxine (lbp) sur le chr. 2 pour l’hypothermie initiale en réponse au SI (LOD=4); et de HSP90 sur le chr. 12 pour l’excrétion de Ca++ (LOD=4.6). Des banques de données sont ensuite interrogées pour recenser les polymorphismes des régions régulatrices ou codantes des gènes candidats chez les souches ancestrales A/J et B6, dont les séquences sont disponibles pour le génome entier. Des utilitaires web permettent de dévoiler les changements dans la structure secondaire de l’ARNm, l’interférence avec des microARN ou avec d’autres motifs de liaison. Plusieurs SNPs fonctionnels ont été identifiés pour le QTL du chr. 1, particulièrement dans les éléments de régulation; ceux-ci impliquant des gènes reliés avec les réponses inflammatoire/immunitaire ou avec le système cardiovasculaire. La quantification par la PCR confirme une régulation à la baisse d’atp1a2 dans le cœur et le cerveau des souches susceptibles à l’anxiété. Ces résultats confirment l’intrication des altérations de la susceptibilité au stress et de la régulation de la TA.

Identification de QTL à pression artérielle dans le chromosome 18 du rat et analyse des gènes candidats Adrb2 et Nedd4l associés à l’hypertension essentielle

L’hypertension constitue un facteur majeur de risque de maladies cardiovasculaires et touche à un pourcentage important de la population humaine. Il s’agit d’une maladie complexe étant donné son caractère multifactoriel. La régulation de la pression artérielle (PA) est sous contrôle de plusieurs gènes, appelés loci pour traits quantitatifs ou QTL, et elle est le résultat de leur interaction. Étant donné que la PA est un trait quantitatif sous contrôle de plusieurs composantes comme les facteurs génétiques et environnementaux, l’étude de l’hypertension est limitée chez les populations humaines. Ainsi la stratégie principale pour l’étude de cette maladie est l’identification de QTL à PA chez des souches congéniques de rat construites à partir des lignées hyper- et normotendues; à savoir les souches Dahl salt-sensitive [1] et Lewis, respectivement. Des études précédentes dans notre laboratoire ont localisé trois QTL à PA au niveau du chromosome 18 chez le rat. Au cours de ce projet, de nouvelles sous-souches ont été construites afin de raffiner la cartographie de ces QTL. Ainsi, les C18QTL1, C18QTL3 et C18QTL4 ont été définis. Des analyses moléculaires ont été effectuées sur deux gènes candidats pour le C18QTL3; à savoir, Adrb2 et Nedd4l associés précédemment à l’hypertension. La comparaison des résultats de séquençage des régions régulatrices et codantes de ces deux gènes, ainsi que leur analyse d’expression par qRT-PCR chez les souches contrastantes DSS et Lewis, n’ont pas montré de différence significative pouvant expliquer la variation du phénotype observé. Des études plus poussées devront être effectuées sur ces deux gènes et, le cas échéant, l’analyse d’autres gènes contenus dans le C18QTL3 devra être entamée afin d’identifier le gène responsable de ce QTL.

Études de réseaux d’expression génique : utilité pour l’élucidation des déterminants génétiques des traits complexes

Les traits quantitatifs complexes sont des caractéristiques mesurables d’organismes vivants qui résultent de l’interaction entre plusieurs gènes et facteurs environnementaux. Les locus génétiques liés à un caractère complexe sont appelés «locus de traits quantitatifs » (QTL). Récemment, en considérant les niveaux d’expression tissulaire de milliers de gènes comme des traits quantitatifs, il est devenu possible de détecter des «QTLs d’expression» (eQTL). Alors que ces derniers ont été considérés comme des phénotypes intermédiaires permettant de mieux comprendre l’architecture biologique des traits complexes, la majorité des études visent encore à identifier une mutation causale dans un seul gène. Cette approche ne peut remporter du succès que dans les situations où le gène incriminé a un effet majeur sur le trait complexe, et ne permet donc pas d’élucider les situations où les traits complexes résultent d’interactions entre divers gènes. Cette thèse propose une approche plus globale pour : 1) tenir compte des multiples interactions possibles entre gènes pour la détection de eQTLs et 2) considérer comment des polymorphismes affectant l’expression de plusieurs gènes au sein de groupes de co-expression pourraient contribuer à des caractères quantitatifs complexes. Nos contributions sont les suivantes : Nous avons développé un outil informatique utilisant des méthodes d’analyse multivariées pour détecter des eQTLs et avons montré que cet outil augmente la sensibilité de détection d’une classe particulière de eQTLs. Sur la base d’analyses de données d’expression de gènes dans des tissus de souris recombinantes consanguines, nous avons montré que certains polymorphismes peuvent affecter l’expression de plusieurs gènes au sein de domaines géniques de co-expression. En combinant des études de détection de eQTLs avec des techniques d’analyse de réseaux de co-expression de gènes dans des souches de souris recombinantes consanguines, nous avons montré qu’un locus génétique pouvait être lié à la fois à l’expression de plusieurs gènes au niveau d’un domaine génique de co-expression et à un trait complexe particulier (c.-à-d. la masse du ventricule cardiaque gauche). Au total, nos études nous ont permis de détecter plusieurs mécanismes par lesquels des polymorphismes génétiques peuvent être liés à l’expression de plusieurs gènes, ces derniers pouvant eux-mêmes être liés à des traits quantitatifs complexes.

Domestication of plants in the americas: Insights from mendelian and molecular genetics

Background Plant domestication occurred independently in four different regions of the Americas. In general, different species were domesticated in each area, though a few species were domesticated independently in more than one area. The changes resulting from human selection conform to the familiar domestication syndrome, though different traits making up this syndrome, for example loss of dispersal, are achieved by different routes in crops belonging to different families. Genetic and Molecular Analyses of Domestication Understanding of the genetic control of elements of the domestication syndrome is improving as a result of the development of saturated linkage maps for major crops, identification and mapping of quantitative trait loci, cloning and sequencing of genes or parts of genes, and discoveries of widespread orthologies in genes and linkage groups within and between families. As the modes of action of the genes involved in domestication and the metabolic pathways leading to particular phenotypes become better understood, it should be possible to determine whether similar phenotypes have similar underlying genetic controls, or whether human selection in genetically related but independently domesticated taxa has fixed different mutants with similar phenotypic effects. Conclusions Such studies will permit more critical analysis of possible examples of multiple domestications and of the origin(s) and spread of distinctive variants within crops. They also offer the possibility of improving existing crops, not only major food staples but also minor crops that are potential export crops for developing countries or alternative crops for marginal areas.

QTLs for shelf life in lettuce co-locate with those for leaf biophysical properties but not with those for leaf developmental traits

Developmental and biophysical leaf characteristics that influence post-harvest shelf life in lettuce, an important leafy crop, have been examined. The traits were studied using 60 informative F-9 recombinant inbed lines (RILs) derived from a cross between cultivated lettuce (Lactuca sativa cv. Salinas) and wild lettuce (L. serriola acc. UC96US23). Quantitative trait loci (QTLs) for shelf life co-located most closely with those for leaf biophysical properties such as plasticity, elasticity, and breakstrength, suggesting that these are appropriate targets for molecular breeding for improved shelf life. Significant correlations were found between shelf life and leaf size, leaf weight, leaf chlorophyll content, leaf stomatal index, and epidermal cell number per leaf, indicating that these pre-harvest leaf development traits confer post-harvest properties. By studying the population in two contrasting environments in northern and southern Europe, the genotype by environment interaction effects of the QTLs relevant to leaf development and shelf life were assessed. In total, 107 QTLs, distributed on all nine linkage groups, were detected from the 29 traits. Only five QTLs were common in both environments. Several areas where many QTLs co-located (hotspots) on the genome were identified, with relatively little overlap between developmental hotspots and those relating to shelf life. However, QTLs for leaf biophysical properties (breakstrength, plasticity, and elasticity) and cell area correlated well with shelf life, confirming that the ideal ideotype lettuce should have small cells with strong cell walls. The identification of QTLs for leaf development, strength, and longevity will lead to a better understanding of processability at a genetic and cellular level, and allow the improvement of salad leaf quality through marker-assisted breeding.