Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Médecine génomique et Maladies infectieuses [Genomic medicine and infectious diseases].

Relentless progress in our knowledge of the nature and functional consequences of human genetic variation allows for a better understanding of the protracted battle between pathogens and their human hosts. Multiple polymorphisms have been identified that impact our response to infections or to anti-infective drugs, and some of them are already used in the clinic. However, to make personalized medicine a reality in infectious diseases, a sustained effort is needed not only in research but also in genomic education.

Natural variability of in vitro adherence to fibrinogen and fibronectin does not correlate with in vivo infectivity of Staphylococcus aureus.

Adherence to fibrinogen and fibronectin plays a crucial role in Staphylococcus aureus experimental endocarditis. Previous genetic studies have shown that infection and carriage isolates do not systematically differ in their virulence-related genes, including genes conferring adherence, such as clfA and fnbA. We set out to determine the range of adherence phenotypes in carriage isolates of S. aureus, to compare the adherence of these isolates to the adherence of infection isolates, and to determine the relationship between adherence and infectivity in a rat model of experimental endocarditis. A total of 133 healthy carriage isolates were screened for in vitro adherence to fibrinogen and fibronectin, and 30 isolates were randomly chosen for further investigation. These 30 isolates were compared to 30 infective endocarditis isolates and 30 blood culture isolates. The infectivities of the carriage isolates, which displayed either extremely low or high adherence to fibrinogen and fibronectin, were tested using a rat model of experimental endocarditis. The levels of adherence to both fibrinogen and fibronectin were very similar for isolates from healthy carriers and members of the two groups of infection isolates. All three groups of isolates showed a wide range of adherence to fibrinogen and fibronectin. Moreover, the carriage isolates that showed minimal adherence and the carriage isolates that showed strong adherence had the same infectivity in experimental endocarditis. Adherence was proven to be important for pathogenesis in experimental endocarditis, but even the least adherent carriage strains had the ability to induce infection. We discuss the roles of differential gene expression, human host factors, and gene redundancy in resolving this apparent paradox.

Genes mirror geography within Europe.

Understanding the genetic structure of human populations is of fundamental interest to medical, forensic and anthropological sciences. Advances in high-throughput genotyping technology have markedly improved our understanding of global patterns of human genetic variation and suggest the potential to use large samples to uncover variation among closely spaced populations. Here we characterize genetic variation in a sample of 3,000 European individuals genotyped at over half a million variable DNA sites in the human genome. Despite low average levels of genetic differentiation among Europeans, we find a close correspondence between genetic and geographic distances; indeed, a geographical map of Europe arises naturally as an efficient two-dimensional summary of genetic variation in Europeans. The results emphasize that when mapping the genetic basis of a disease phenotype, spurious associations can arise if genetic structure is not properly accounted for. In addition, the results are relevant to the prospects of genetic ancestry testing; an individual's DNA can be used to infer their geographic origin with surprising accuracy-often to within a few hundred kilometres.

Intraspecific competition reveals conditional fitness effects of single gene polymorphism at the Arabidopsis root growth regulator BRX.

Intraspecific genetic variation for morphological traits is observed in many organisms. In Arabidopsis thaliana, alleles responsible for intraspecific morphological variation are increasingly being identified. However, the fitness consequences remain unclear in most cases. Here, the fitness effects of alleles of the BRX gene are investigated. A brx loss-of-function allele, which was found in a natural accession, results in a highly branched but poorly elongated root system. Comparison between the control accession Sav-0 and an introgression of brx into this background (brxS) indicated that, surprisingly, brx loss of function did not negatively affect fitness in pure stands. However, in mixed, well-watered stands brxS performance and reproductive output decreased significantly, as the proportion of Sav-0 neighbors increased. Additional comparisons between brxS and a brxS line that was complemented by a BRX transgene confirmed a direct effect of the loss-of-function allele on plant performance, as indicated by restored competitive ability of the transgenic genotype. Further, because plant height was very similar across genotypes and because the experimental setup largely excluded shading effects, the impaired competitiveness of the brx loss-of-function genotype likely reflects below-ground competition. In summary, these data reveal conditional fitness effects of a single gene polymorphism in response to intraspecific competition in Arabidopsis.

Evidence for the evolution of multiple genomes in arbuscular mycorrhizal fungi.

Ancient asexuals directly contradict the evolutionary theories that explain why organisms should evolve a sexual life history. The mutualistic, arbuscular mycorrhizal fungi are thought to have been asexual for approximately 400 million years. In the absence of sex, highly divergent descendants of formerly allelic nucleotide sequences are thought to evolve in a genome. In mycorrhizal fungi, where individual offspring receive hundreds of nuclei from the parent, it has been hypothesized that a population of genetically different nuclei should evolve within one individual. Here we use DNA-DNA fluorescent in situ hybridization to show that genetically different nuclei co-exist in individual arbuscular mycorrhizal fungi. We also show that the population genetics techniques used in other organisms are unsuitable for detecting recombination because the assumptions and underlying processes do not fit the fungal genomic structure shown here. Instead we used a phylogenetic approach to show that the within-individual genetic variation that occurs in arbuscular mycorrhizal fungi probably evolved through accumulation of mutations in an essentially clonal genome, with some infrequent recombination events. We conclude that mycorrhizal fungi have evolved to be multi-genomic.

Testing demographic models of effective population size.

Due to practical difficulties in obtaining direct genetic estimates of effective sizes, conservation biologists have to rely on so-called 'demographic models' which combine life-history and mating-system parameters with F-statistics in order to produce indirect estimates of effective sizes. However, for the same practical reasons that prevent direct genetic estimates, the accuracy of demographic models is difficult to evaluate. Here we use individual-based, genetically explicit computer simulations in order to investigate the accuracy of two such demographic models aimed at investigating the hierarchical structure of populations. We show that, by and large, these models provide good estimates under a wide range of mating systems and dispersal patterns. However, one of the models should be avoided whenever the focal species' breeding system approaches monogamy with no sex bias in dispersal or when a substructure within social groups is suspected because effective sizes may then be strongly overestimated. The timing during the life cycle at which F-statistics are evaluated is also of crucial importance and attention should be paid to it when designing field sampling since different demographic models assume different timings. Our study shows that individual-based, genetically explicit models provide a promising way of evaluating the accuracy of demographic models of effective size and delineate their field of applicability.

The HSR on chromosome 1 of the house mouse, Mus domesticus: distribution and frequency in Switzerland.

A total of 357 house mice (Mus domesticus) from 83 localities uniformly distributed throughout Switzerland were screened for the presence of a homogenously staining region (HSR) on chromosome 1. Altogether 47 mice from 11 localities were HSR/+ or HSR/HSR. One sample of 11 individuals all had an HSR/HSR karyotype. Almost all mice with the variant were collected from the Rhone valley (HSR frequency: 61%) and Val Bregaglia (HSR frequency: 81%). For samples from most of the area of Switzerland, the HSR was absent. There was no strong association between the geographic distribution of the HSR and the areas of occurrence of metacentrics. However, at Chiggiogna the HSR was found on Rb (1.3). Possible explanations for the HSR polymorphism are discussed.

Fanconi-Bickel syndrome and autosomal recessive proximal tubulopathy with hypercalciuria (ARPTH) are allelic variants caused by GLUT2 mutations.

CONTEXT: Many inherited disorders of calcium and phosphate homeostasis are unexplained at the molecular level. OBJECTIVE: The objective of the study was to identify the molecular basis of phosphate and calcium abnormalities in two unrelated, consanguineous families. PATIENTS: The affected members in family 1 presented with rickets due to profound urinary phosphate-wasting and hypophosphatemic rickets. In the previously reported family 2, patients presented with proximal renal tubulopathy and hypercalciuria yet normal or only mildly increased urinary phosphate excretion. METHODS: Genome-wide linkage scans and direct nucleotide sequence analyses of candidate genes were performed. Transport of glucose and phosphate by glucose transporter 2 (GLUT2) was assessed using Xenopus oocytes. Renal sodium-phosphate cotransporter 2a and 2c (Npt2a and Npt2c) expressions were evaluated in transgenically rescued Glut2-null mice (tgGlut2-/-). RESULTS: In both families, genetic mapping and sequence analysis of candidate genes led to the identification of two novel homozygous mutations (IVS4-2A>G and R124S, respectively) in GLUT2, the gene mutated in Fanconi-Bickel syndrome, a rare disease usually characterized by renal tubulopathy, impaired glucose homeostasis, and hepatomegaly. Xenopus oocytes expressing the [R124S]GLUT2 mutant showed a significant reduction in glucose transport, but neither wild-type nor mutant GLUT2 facilitated phosphate import or export; tgGlut2-/- mice demonstrated a profound reduction of Npt2c expression in the proximal renal tubules. CONCLUSIONS: Homozygous mutations in the facilitative glucose transporter GLUT2, which cause Fanconi-Bickel syndrome, can lead to very different clinical and biochemical findings that are not limited to mild proximal renal tubulopathy but can include significant hypercalciuria and highly variable degrees of urinary phosphate-wasting and hypophosphatemia, possibly because of the impaired proximal tubular expression of Npt2c.

Flexible social organization and high incidence of drifting in the sweat bee, Halictus scabiosae.

The very diverse social systems of sweat bees make them interesting models to study social evolution. Here we focus on the dispersal behaviour and social organization of Halictus scabiosae, a common yet poorly known species of Europe. By combining field observations and genetic data, we show that females have multiple reproductive strategies, which generates a large diversity in the social structure of nests. A detailed microsatellite analysis of 60 nests revealed that 55% of the nests contained the offspring of a single female, whereas the rest had more complex social structures, with three clear cases of multiple females reproducing in the same nest and frequent occurrence of unrelated individuals. Drifting among nests was surprisingly common, as 16% of the 122 nests in the overall sample and 44% of the nests with complex social structure contained females that had genotypes consistent with being full-sisters of females sampled in other nests of the population. Drifters originated from nests with an above-average productivity and were unrelated to their nestmates, suggesting that drifting might be a strategy to avoid competition among related females. The sex-specific comparison of genetic differentiation indicated that dispersal was male-biased, which would reinforce local resource competition among females. The pattern of genetic differentiation among populations was consistent with a dynamic process of patch colonization and extinction, as expected from the unstable, anthropogenic habitat of this species. Overall, our data show that H. scabiosae varies greatly in dispersal behaviour and social organization. The surprisingly high frequency of drifters echoes recent findings in wasps and bees, calling for further investigation of the adaptive basis of drifting in the social insects.

Genomewide association study using a high-density single nucleotide polymorphism array and case-control design identifies a novel essential hypertension susceptibility locus in the promoter region of endothelial NO synthase.

Essential hypertension is a multifactorial disorder and is the main risk factor for renal and cardiovascular complications. The research on the genetics of hypertension has been frustrated by the small predictive value of the discovered genetic variants. The HYPERGENES Project investigated associations between genetic variants and essential hypertension pursuing a 2-stage study by recruiting cases and controls from extensively characterized cohorts recruited over many years in different European regions. The discovery phase consisted of 1865 cases and 1750 controls genotyped with 1M Illumina array. Best hits were followed up in a validation panel of 1385 cases and 1246 controls that were genotyped with a custom array of 14 055 markers. We identified a new hypertension susceptibility locus (rs3918226) in the promoter region of the endothelial NO synthase gene (odds ratio: 1.54 [95% CI: 1.37-1.73]; combined P=2.58 · 10(-13)). A meta-analysis, using other in silico/de novo genotyping data for a total of 21 714 subjects, resulted in an overall odds ratio of 1.34 (95% CI: 1.25-1.44; P=1.032 · 10(-14)). The quantitative analysis on a population-based sample revealed an effect size of 1.91 (95% CI: 0.16-3.66) for systolic and 1.40 (95% CI: 0.25-2.55) for diastolic blood pressure. We identified in silico a potential binding site for ETS transcription factors directly next to rs3918226, suggesting a potential modulation of endothelial NO synthase expression. Biological evidence links endothelial NO synthase with hypertension, because it is a critical mediator of cardiovascular homeostasis and blood pressure control via vascular tone regulation. This finding supports the hypothesis that there may be a causal genetic variation at this locus.

Comparison of human chromosome 21 conserved nongenic sequences (CNGs) with the mouse and dog genomes shows that their selective constraint is independent of their genic environment.

The analysis of conservation between the human and mouse genomes resulted in the identification of a large number of conserved nongenic sequences (CNGs). The functional significance of this nongenic conservation remains unknown, however. The availability of the sequence of a third mammalian genome, the dog, allows for a large-scale analysis of evolutionary attributes of CNGs in mammals. We have aligned 1638 previously identified CNGs and 976 conserved exons (CODs) from human chromosome 21 (Hsa21) with their orthologous sequences in mouse and dog. Attributes of selective constraint, such as sequence conservation, clustering, and direction of substitutions were compared between CNGs and CODs, showing a clear distinction between the two classes. We subsequently performed a chromosome-wide analysis of CNGs by correlating selective constraint metrics with their position on the chromosome and relative to their distance from genes. We found that CNGs appear to be randomly arranged in intergenic regions, with no bias to be closer or farther from genes. Moreover, conservation and clustering of substitutions of CNGs appear to be completely independent of their distance from genes. These results suggest that the majority of CNGs are not typical of previously described regulatory elements in terms of their location. We propose models for a global role of CNGs in genome function and regulation, through long-distance cis or trans chromosomal interactions.

Regional variations in ABC transporter expression along the mouse intestinal tract.

The ATP-binding cassette (ABC) family of proteins comprise a group of membrane transporters involved in the transport of a wide variety of compounds, such as xenobiotics, vitamins, lipids, amino acids, and carbohydrates. Determining their regional expression patterns along the intestinal tract will further characterize their transport functions in the gut. The mRNA expression levels of murine ABC transporters in the duodenum, jejunum, ileum, and colon were examined using the Affymetrix MuU74v2 GeneChip set. Eight ABC transporters (Abcb2, Abcb3, Abcb9, Abcc3, Abcc6, Abcd1, Abcg5, and Abcg8) displayed significant differential gene expression along the intestinal tract, as determined by two statistical models (a global error assessment model and a classic ANOVA, both with a P < 0.01). Concordance with semiquantitative real-time PCR was high. Analyzing the promoters of the differentially expressed ABC transporters did not identify common transcriptional motifs between family members or with other genes; however, the expression profile for Abcb9 was highly correlated with fibulin-1, and both genes share a common complex promoter model involving the NFkappaB, zinc binding protein factor (ZBPF), GC-box factors SP1/GC (SP1F), and early growth response factor (EGRF) transcription binding motifs. The cellular location of another of the differentially expressed ABC transporters, Abcc3, was examined by immunohistochemistry. Staining revealed that the protein is consistently expressed in the basolateral compartment of enterocytes along the anterior-posterior axis of the intestine. Furthermore, the intensity of the staining pattern is concordant with the expression profile. This agrees with previous findings in which the mRNA, protein, and transport function of Abcc3 were increased in the rat distal intestine. These data reveal regional differences in gene expression profiles along the intestinal tract and demonstrate that a complete understanding of intestinal ABC transporter function can only be achieved by examining the physiologically distinct regions of the gut.

Influence of learning on range expansion and adaptation to novel habitats.

Learning has been postulated to 'drive' evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual-based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically-based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.

Evolutionary forces affect multigenomic arbuscular mycorhizal fungi : an experimental approach using Glomus intraradices

Abstract Arbuscular Mycorhizal Fungi (AMF) are important plant symbionts that can improve floristic diversity and ecosystem productivity. These important fungi are obligate biotrophs and form symbioses with roots of the majority of plant species, improving plant nutrient acquisition in exchange of photosynthates. AM fungi are successful both ecologically as they occupy a very large spectrum of environments as well as host range and evolutionarily, as this symbiosis is over 400 million years old. These fungi grow and reproduce clonally by hyphae and multinucleate spores. AMF are coenocytic and recent work has shown that they harbor genetically different nuclei and that AMF populations are genetically diverse. How AMF species diversity is maintained has been addressed theoretically and experimentally at the community level. Much less attention has been drawn to understand how genetic diversity is maintained within populations although closely related individuals are more likely to compete for the same resources and occupy similar niches. How infra-individual genetic diversity is shaped and maintained has received even less attention. In Chapter 2, we show that individuals from a field population may differ in their symbiotic efficiency under reduced phosphate availability: We show there is genetic variation in an AMF field population for fitness-related growth traits in response to different phosphate availability acid host species. Furthermore, AFLP fingerprints of the same individuals growing in contrasting environments diverged suggesting that the composition in nuclei of AMF is dynamical and affected by environmental factors. Thus environmental heterogeneity is likely to play an important role for the maintenance of genetic diversity at the population level. In Chapter 3 we show that single spores do not inherit necessarily the same genetic material. We have found genetic divergences using two different types of molecular marker, as well as phenotypic divergences among single spore lines. Our results stress the importance of considering these organisms as a multilevel hierarchical system and of better knowing their life cycle. They have important consequences for the understanding of AMF genetics, ecology and the development of commercial AMF inocculum. Résumé Les champignons endomycorhiziens arbusculaires (CEA) sont d'importants symbiontes pour les plantes, car ils augmentent la diversité et la productivité des écosystèmes. Ces importants symbiontes sont des biotrophes obligatoires et forment une symbiose avec la plupart des plantes terrestres. Ils améliorent l'acquisition de substances nutritives de leurs hôtes en échange de sucres obtenus par photosynthèse. Ces champignons ont un grand succès écologique, ils colonisent une grande rangée d'environnements ainsi que d'hôtes. Ils ont aussi un succès évolutif certain de part le fait que cette symbiose existe depuis plus de 400 millions d'années. Les CEA sont asexués et croissent clonalement en formant des hyphes et des spores multinuclées. Les CEA sont des coenocytes et des travaux de recherche récents ont montré qu'ils possèdent des noyaux génétiquement différents. D'autres travaux ont aussi révélé que les populations de CEA sont génétiquement diversifiées. Comment la diversité des CEA est maintenue a seulement été adressée par des études théoriques et expérimentalement au niveau des communautés. Très peu d'attention a été portée sur le maintien de la diversité génétique infra et inter populationnelle, or ce sont les individus les plus proches génétiquement qui vont entrer en compétition pour des ressources et niches similaires. La formation et le maintien de la diversité intra-individu des CEA a reçu très peu d'attention. Dans le chapitre 2, nous montrons que des individus CEA d'un même champ différent dans leur efficacité symbiotique lorsque la concentration en phosphoré est réduite. Nous montrons qu'il existe de la variance génétique dans une population de CEA provenant d'un même champ en réponse à différentes concentrations de phosphore, ainsi qu'en réponse à différentes espèces d'hôtes, et ceci pour des traits de croissance vraisemblablement liés au succès reproducteur. De plus grâce à des AFLP nous avons pu montrer que le génome de ces individus subissent des changements lorsqu'ils croissent dans des environnements contrastés. Ceci suggère que les noyaux génétiquement différents des CEA sont des entités dynamiques. Il est fort probable que l'hétérogénéité environnementale joue un rôle dans le maintien de la diversité génétique des populations de CEA. Dans le chapitre 3, nous montrons que toutes les spores d'un même mycélium parental de CEA ne reçoivent pas exactement le même contenu génétique. Nous avons mis en évidence des divergences entre des Lignées monosporales en utilisant deux types de marqueur moléculaires, ainsi que des différences phénotypiques. Nos résutats soulignent l'importance de considézer ces organismes comme dés systëmes hiérarchiques mufti-niveaux, ainsi que de mieux connaître leur cycle de vie. Nos résultats ont d'importantes conséquences pour la compréhension du système génétique des CEA, ainsi que de leur évolution, leur écológie, mais également des conséquences pour la production d' inoccultim commercial.