Functional impact of genetic variation on gene expression

Numerous links between genetic variants and phenotypes are known and genome-wide association studies dramatically increased the number of genetic variants associated with traits during the last decade. However, how changes in the DNA perturb the molecular mechanisms and impact on the phenotype of an organism remains elusive. Studies suggest that many traitassociated variants are in the non-coding region of the genome and probably act through regulation of gene expression. During my thesis I investigated how genetic variants affect gene expression through gene regulatory mechanisms. The first chapter was a collaborative project with a pharmaceutical company, where we investigated genome-wide copy number variation (CNVs) among Cynomolgus monkeys (Macaca fascicularis) used in pharmaceutical studies, and associated them to changes in gene expression. We found substantial copy number variation and identified CNVs linked to tissue-specific expression changes of proximal genes. The second and third chapters focus on genetic variation in humans and its effects on gene regulatory mechanisms and gene expression. The second chapter studies two human trios, where the allelic effects of genetic variation on genome-wide gene expression, protein-DNA binding and chromatin modifications were investigated. We found abundant allele specific activity across all measured molecular phenotypes and show extended coordinated behavior among them. In the third chapter, we investigated the impact of genetic variation on these phenotypes in 47 unrelated individuals. We found that chromatin phenotypes are organized into local variable modules, often linked to genetic variation and gene expression. Our results suggest that chromatin variation emerges as a result of perturbations of cis-regulatory elements by genetic variants, leading to gene expression changes. The work of this thesis provides novel insights into how genetic variation impacts gene expression by perturbing regulatory mechanisms. -- De nombreux liens entre variations génétiques et phénotypes sont connus. Les études d'association pangénomique ont considérablement permis d'augmenter le nombre de variations génétiques associées à des phénotypes au cours de la dernière décennie. Cependant, comprendre comment ces changements perturbent les mécanismes moléculaires et affectent le phénotype d'un organisme nous échappe encore. Des études suggèrent que de nombreuses variations, associées à des phénotypes, sont situées dans les régions non codantes du génome et sont susceptibles d'agir en modifiant la régulation d'expression des gènes. Au cours de ma thèse, j'ai étudié comment les variations génétiques affectent les niveaux d'expression des gènes en perturbant les mécanismes de régulation de leur expression. Le travail présenté dans le premier chapitre est un projet en collaboration avec une société pharmaceutique. Nous avons étudié les variations en nombre de copies (CNV) présentes chez le macaque crabier (Macaca fascicularis) qui est utilisé dans les études pharmaceutiques, et nous les avons associées avec des changements d'expression des gènes. Nous avons découvert qu'il existe une variabilité substantielle du nombre de copies et nous avons identifié des CNVs liées aux changements d'expression des gènes situés dans leur voisinage. Ces associations sont présentes ou absentes de manière spécifique dans certains tissus. Les deuxième et troisième chapitres se concentrent sur les variations génétiques dans les populations humaines et leurs effets sur les mécanismes de régulation des gènes et leur expression. Le premier se penche sur deux trios humains, père, mère, enfant, au sein duquel nous avons étudié les effets alléliques des variations génétiques sur l'expression des gènes, les liaisons protéine-ADN et les modifications de la chromatine. Nous avons découvert que l'activité spécifique des allèles est abondante abonde dans tous ces phénotypes moléculaires et nous avons démontré que ces derniers ont un comportement coordonné entre eux. Dans le second, nous avons examiné l'impact des variations génétiques de ces phénotypes moléculaires chez 47 individus, sans lien de parenté. Nous avons observé que les phénotypes de la chromatine sont organisés en modules locaux, qui sont liés aux variations génétiques et à l'expression des gènes. Nos résultats suggèrent que la variabilité de la chromatine est due à des variations génétiques qui perturbent des éléments cis-régulateurs, et peut conduire à des changements dans l'expression des gènes. Le travail présenté dans cette thèse fournit de nouvelles pistes pour comprendre l'impact des différentes variations génétiques sur l'expression des gènes à travers les mécanismes de régulation.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps.

BACKGROUND: Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation. RESULTS: We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages. CONCLUSIONS: Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Genetic variation among pathogens causing "Helminthosporium" diseases of rice, maize and wheat

Twenty isolates of four fungal species, agents of "Helminthosporium" diseases in cereals, were collected from different regions: nine Bipolarisoryzae isolated from rice (Oryza sativa), seven B.sorokiniana from wheat (Triticum aestivum), two B. maydis, and two Exserohilumturcicum from maize (Zea mays). The strains were compared by PCR-RFLP and RAPD analysis. Size polymorphism among the isolates in the ITS region comprising the 5.8 S rDNA indicated genetic differences among the isolates, while a UPGMA phenogram constructed after the digestion of this region with restriction enzymes showed inter- and intra-specific polymorphism. The RAPD profiles indicated an expressive level of polymorphism among different species, compared with a low level of polymorphism among isolates of the same species. A UPGMA phenogram grouped the isolates according to the species and their host plant. RAPD profiles did not reveal polymorphism that directly correlated climatic factors with geographic source of the isolates of B. sorokiniana, and B. oryzae. Teleomorphic species revealed high similarity with their correspondent anamorphs.

Genetic variation and evolutionary divergence within and among populations, species, and genera of the Cambarinae

Seven crayfish species from three genera of the subfamily Cambarinae were electrophoretically examined for genetic variation at a total of twenty-six loci. Polymorphism was detected primarily at three loci: Ao-2, Lap, and Pgi. The average heterozygosities over-all loci for each species were found to be very low when compared to most other invertebrate species that have been examined electrophoretically. With the exception of Cambarus bartoni, the interpopulation genetic identities are high within any given species. The average interspecific identities are somewhat lower and the average intergeneric identities are lower still. Populations, species and genera conform to the expected taxonomic progression. The two samples of ~ bartoni show high genetic similarity at only 50 percent of the loci compared. Locus by locus identity comparisons among species yield U-shaped distributions of genetic identities. Construction of a phylogenetic dendrogram using species mean genetic distances values shows that species grouping is in agreement with morphological taxonomy with the exception of the high similarity between Orconectespropinquus and Procambarus pictus. This high similarity suggests the possibility of a regulatory change between the two species. It appears that the low heterozygosities, high interpopulation genetic identities, and taxonomic mispositioning can all be explained on the basis of low mutation rates.

Genetic variation of and environmental effects on inducibility of resistance in tomatoes (Solanum lycopersicum L.) to Phytophthora infestans (Mont) de Bary

Many plant strengtheners are promoted for their supposed effects on nutrient uptake and/or resistance induction (IR). In addition, many organic fertilizers are supposed to enhance plant health and several studies have shown that tomatoes grown organically are more resistant to late blight, caused by Phytophthora infestans to tomatoes grown conventionally. Much is known about the mechanisms underlying IR. In contrast, there is no systematic knowledge about genetic variation for IR. Therefore, the following questions were addressed in the presented dissertation: (i) Is there genetic variation among tomato genotypes for inducibility of resistance to P. infestans? (ii) How do different PS compare with the chemical inducer BABA in their ability to IR? (iii) Does IR interact with the inducer used and different organic fertilizers? A varietal screening showed that contrary to the commonly held belief IR in tomatoes is genotype and isolate specific. These results indicate that it should be possible to select for inducibility of resistance in tomato breeding. However, isolate specificity also suggests that there could be pathogen adaptation. The three tested PS as well as two of the three tested organic fertilisers all induced resistance in the tomatoes. Depending on PS or BABA variety and isolate effects varied. In contrast, there were no variety and isolate specific effects of the fertilisers and no interactions with the PS and fertilisers. This suggests that the different PS should work independent of the soil substrate used. In contrast the results were markedly different when isolate mixtures were used for challenge inoculations. Plants were generally less susceptible to isolate mixtures than to single isolates. In addition, the effectiveness of the PS was greater and more similar to BABA when isolate mixtures were used. The fact that the different PS and BABA differed in their ability to induce resistance in different host genotype -pathogen isolate combinations puts the usefulness of IR as a breeding goal in question. This would result in varieties depending on specific inducers. The results with the isolate mixtures are highly relevant. On the one hand they increase the effectiveness of the resistance inducers. On the other hand, measures that increase the pathogen diversity such as the use of diversified host populations will also increase the overall resistance of the hosts. For organic tomato production the results indicate that it is possible to enhance the tomato growing system with respect to plant health management by using optimal fertilisers, plant strengtheners and any measures that increase system diversity.

Inbreeding and loss of genetic variation in a reintroduced population of Mauritius Kestrel

Many populations have recovered from severe bottlenecks either naturally or through intensive conservation management. In the past, however, few conservation programs have monitored the genetic health of recovering populations. We conducted a conservation genetic assessment of a small, reintroduced population of Mauritius Kestrel (Falco punctatus) to determine whether genetic deterioration has occurred since its reintroduction. We used pedigree analysis that partially accounted for individuals of unknown origin to document that (1) inbreeding occurred frequently (2.6% increase per generation; N-el = 18.9), (2) 25% of breeding pairs were composed of either closely or moderately related individuals, (3) genetic diversity has been lost from the population (1,6% loss per generation; N-ev = 32.1) less rapidly than the corresponding increase in inbreeding, and (4) ignoring the contribution of unknown individuals to a pedigree will bias the metrics derived from that pedigree, ultimately obscuring the prevailing genetic dynamics. The rates of inbreeding and loss of genetic variation in the subpopulation of Mauritius Kestrel we examined were extreme and among the highest yet documented in a wild vertebrate population. Thus, genetic deterioration may affect this population's long-term viability. Remedial conservation strategies are needed to reduce the impact of inbreeding and loss of genetic variation in this species, We suggest that schemes to monitor genetic variation after reintroduction should be an integral component of endangered species recovery programs

Genetic variation in a chilean endangered endemic: Gomortega keule (Molina) Baillon

Gomortega keule (Molina) Baillon is an endangered, rare species, the only representative of its genus, and endemic to Central Chile. Populations of this tree are now fragmented and few individuals can be found in any of them. Genetic diversity was studied in 33 individuals from three populations in Cauquenes, a coastal mountain area (35°58'S-72°41'W). Fifteen InterSimple Sequence Repeat primers were used to determine the degree of similarity between and within populations. This revealed that 30% of the variation exhibited was between populations while 70% was within; nevertheless individuals were clearly clustered in a pattern which reflected a narrow base of diversity. Three other species from the Laurales order were used in order to provide an external reference as to the degree of diversity. In addition, an external wild population from the native species, Peumus boldus, was used to verify the utility of the markers. We show that the primers are effective in quickly giving an estimate of the degree of diversity of a population, thus giving important topical information relevant to preserving endangered species. Aspects of the conservation and management policy for the species in order to maintain the remaining populations and to preserve the genetic resources are discussed.

Gene-nutrient interactions with dietary fat modulate the association between genetic variation of the ACSL1 gene and metabolic syndrome

Long-chain acyl CoA synthetase 1 (ACSL1) plays an important role in fatty acid metabolism and triacylglycerol (TAG) synthesis. Disturbance of these pathways may result in dyslipidemia and insulin resistance, hallmarks of the metabolic syndrome (MetS). Dietary fat is a key environmental factor that may interact with genetic determinants of lipid metabolism to affect MetS risk. We investigated the relationship between ACSL1 polymorphisms (rs4862417, rs6552828, rs13120078, rs9997745, and rs12503643) and MetS risk and determined potential interactions with dietary fat in the LIPGENE-SU.VI.MAX study of MetS cases and matched controls (n = 1,754). GG homozygotes for rs9997745 had increased MetS risk {odds ratio (OR) 1.90 [confidence interval (CI) 1.15, 3.13]; P = 0.01}, displayed elevated fasting glucose (P = 0.001) and insulin concentrations (P = 0.002) and increased insulin resistance (P = 0.03) relative to the A allele carriers. MetS risk was modulated by dietary fat, whereby the risk conferred by GG homozygosity was abolished among individuals consuming either a low-fat (<35% energy) or a high-PUFA diet (>5.5% energy). In conclusion, ACSL1 rs9997745 influences MetS risk, most likely via disturbances in fatty acid metabolism, which was modulated by dietary fat consumption, particularly PUFA intake, suggesting novel gene-nutrient interactions.

Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles

Quaternary climatic fluctuations have had profound effects on the phylogeographic structure of many species. Classically, species were thought to have become isolated in peninsular refugia, but there is limited evidence that large, non-polar species survived outside traditional refugial areas. We examined the phylogeographic structure of the red fox (Vulpes vulpes), a species that shows high ecological adaptability in the western Palaearctic region. We compared mitochondrial DNA sequences (cytochrome b and control region) from 399 modern and 31 ancient individuals from across Europe. Our objective was to test whether red foxes colonised the British Isles from mainland Europe in the late Pleistocene, or whether there is evidence that they persisted in the region through the Last Glacial Maximum. We found red foxes to show a high degree of phylogeographic structuring across Europe and, consistent with palaeontological and ancient DNA evidence, confirmed via phylogenetic indicators that red foxes were persistent in areas outside peninsular refugia during the last ice age. Bayesian analyses and tests of neutrality indicated population expansion. We conclude that there is evidence that red foxes from the British Isles derived from central European populations that became isolated after the closure of the landbridge with Europe.

Genetic variation at flowering time loci in wild and cultivated barley

The worldwide spread of barley cultivation required adaptation to agricultural environments far distant from those found in its centre of domestication. An important component of this adaptation is the timing of flowering, achieved predominantly in response to day length and temperature. Here, we use a collection of cultivars, landraces and wild barley accessions to investigate the origins and distribution of allelic diversity at four major flowering time loci, mutations at which have been under selection during the spread of barley cultivation into Europe. Our findings suggest that while mutant alleles at the PPD-H1 and PPD-H2 photoperiod loci occurred pre-domestication, the mutant vernalization non-responsive alleles utilized in landraces and cultivars at the VRN-H1 and VRN-H2 loci occurred post-domestication. The transition from wild to cultivated barley is associated with a doubling in the number of observed multi-locus flowering-time haplotypes, suggesting that the resulting phenotypic variation has aided adaptation to cultivation in the diverse ecogeographic locations encountered. Despite the importance of early-flowering alleles during the domestication of barley in Europe, we show that novel VRN alleles associated with early flowering in wild barley have been lost in domesticates, highlighting the potential of wild germplasm as a source of novel allelic variation for agronomic traits.

Natural genetic variation in caesium (Cs) accumulation by Arabidopsis thaliana

Ingestion of caesium (Cs) radioisotopes poses a health risk to humans. Crop varieties that accumulate less Cs in their edible tissues may provide a useful countermeasure. This study was performed to determine whether quantitative genetics on a model plant (Arabidopsis thaliana) might inform such 'safe'-crop strategies. Arabidopsis accessions and recombinant inbred lines (RILs), from Landsberg erecta (Ler) x Cape Verdi Island (Cvi), Ler x Columbia (Col), and Niederzenz (Nd) x Col mapping populations, were grown on agar supplemented with subtoxic levels of Cs. Shoot Cs concentration varied up to three-fold, and shoot f. wt varied up to 25-fold within populations. The heritability of growth and Cs accumulation traits ranged from 0.06 to 0.28. Four quantitative trait loci (QTL) accounted for > 80 of the genetic contribution to the total phenotypic variation in shoot Cs concentration in the Ler x Col population. QTL identified in this study, in particular, QTL co-localizing to the top and bottom regions of Chromosomes I and V in two different mapping populations, are amenable to positional cloning and, through collinearity, may inform selection or breeding strategies for the development of 'safe' crops.

Physical activity attenuates the body mass index-increasing influence of genetic variation in the FTO gene 1-3.

BACKGROUND: Intronic variation in the FTO (fat mass and obesity-associated) gene has been unequivocally associated with increased body mass index (BMI; in kg/m(2)) and the risk of obesity in populations of different ethnicity. OBJECTIVE: We examined whether this robust genetic predisposition to obesity can be attenuated by being more physically active. DESIGN: The FTO variant rs1121980 was genotyped in 20,374 participants (39-79 y of age) from the European Prospective Investigation into Cancer and Nutrition-Norfolk Study, an ethnically homogeneous population-based cohort. Physical activity (PA) was assessed with a validated self-reported questionnaire. The interaction between rs1121980 and PA on BMI and waist circumference (WC) was examined by including the interaction term in mixed-effect models. RESULTS: We confirmed that the risk (T) allele of rs1121980 was significantly associated with BMI (0.31-unit increase per allele; P < 0.001) and WC (0.77-cm increase per allele; P < 0.001). The PA level attenuated the effect of rs1121980 on BMI and WC; ie, whereas in active individuals the risk allele increased BMI by 0.25 per allele, the increase in BMI was significantly (P for interaction = 0.004) more pronounced (76%) in inactive individuals (0.44 per risk allele). We observed similar effects for WC (P for interaction = 0.02): the risk allele increased WC by 1.04 cm per allele in inactive individuals but by only 0.64 cm in active individuals. CONCLUSIONS: Our results showed that PA attenuates the effect of the FTO rs1121980 genotype on BMI and WC. This observation has important public health implications because we showed that a genetic susceptibility to obesity induced by FTO variation can be overcome, at least in part, by adopting a physically active lifestyle.

Genetic variation in GABRB3 is associated with Asperger syndrome and multiple endophenotypes relevant to autism.

BACKGROUND: Autism spectrum conditions (ASC) are associated with deficits in social interaction and communication, alongside repetitive, restricted, and stereotyped behavior. ASC is highly heritable. The gamma-aminobutyric acid (GABA)-ergic system has been associated consistently with atypicalities in autism, in both genetic association and expression studies. A key component of the GABA-ergic system is encoded by the GABRB3 gene, which has been previously implicated both in ASC and in individual differences in empathy. METHODS: In this study, 45 genotyped single nucleotide polymorphisms (SNPs) within GABRB3 were tested for association with Asperger syndrome (AS), and related quantitative traits measured through the following tests: the Empathy Quotient (EQ), the Autism Spectrum Quotient (AQ), the Systemizing Quotient-Revised (SQ-R), the Embedded Figures Test (EFT), the Reading the Mind in the Eyes Test (RMET), and the Mental Rotation Test (MRT). Two-loci, three-loci, four-loci haplotype analyses, and one seven-loci haplotype analysis were also performed in the AS case--control sample. RESULTS: Three SNPs (rs7180158, rs7165604, rs12593579) were significantly associated with AS, and two SNPs (rs9806546, rs11636966) were significantly associated with EQ. Two SNP-SNP pairs, rs12438141-rs1035751 and rs12438141-rs7179514, showed significant association with variation in the EFT scores. One SNP-SNP pair, rs7174437-rs1863455, was significantly associated with variation in the MRT scores. Additionally, a few haplotypes, including a 19 kb genomic region that formed a linkage disequilibrium (LD) block in our sample and contained several nominally significant SNPs, were found to be significantly associated with AS. CONCLUSION: The current study confirms the role of GABRB3 as an important candidate gene in both ASC and normative variation in related endophenotypes.