The case for resequencing studies of Arabidopsis thaliana accessions: mining the dark matter of natural genetic variation.

Ultra-high-throughput sequencing (UHTS) techniques are evolving rapidly and may soon become an affordable and routine tool for sequencing plant DNA, even in smaller plant biology labs. Here we review recent insights into intraspecific genome variation gained from UHTS, which offers a glimpse of the rather unexpected levels of structural variability among Arabidopsis thaliana accessions. The challenges that will need to be addressed to efficiently assemble and exploit this information are also discussed.

Common genetic variation and the control of HIV-1 in humans.

To extend the understanding of host genetic determinants of HIV-1 control, we performed a genome-wide association study in a cohort of 2,554 infected Caucasian subjects. The study was powered to detect common genetic variants explaining down to 1.3% of the variability in viral load at set point. We provide overwhelming confirmation of three associations previously reported in a genome-wide study and show further independent effects of both common and rare variants in the Major Histocompatibility Complex region (MHC). We also examined the polymorphisms reported in previous candidate gene studies and fail to support a role for any variant outside of the MHC or the chemokine receptor cluster on chromosome 3. In addition, we evaluated functional variants, copy-number polymorphisms, epistatic interactions, and biological pathways. This study thus represents a comprehensive assessment of common human genetic variation in HIV-1 control in Caucasians.

Chromosomal rearrangement and genetic structure in the shrews of the Sorex araneus group

ABSTRACT The role of chromosomal rearrangements in the speciation process is much debated and many theoretical models have been developed. The shrews of the Sorex araneus group offer extraordinary opportunities to study the relationship between chromosomal variation and speciation. Indeed, this group of morphologically very similar species received a great deal of attention due to its karyotypic variability, which is mainly attributed to Robertsonian fusions. To explore the impact of karyotypic changes on genetic differentiation, we first studied the relationship between genetic and karyotypic structure among Alpine species and among chromosome races of the S. araneus group using Bayesian admixture analyses. The results of these analyses confirmed the taxonomic status of the studied species even though introgression can still be detected between species. Moreover, the strong spatial sub-structure highlighted the role of historical factors (e.g. geographical isolation) on genetic structure. Next, we studied gene flow at the chromosome level to address the question of the impact of chromosomal rearrangements on genetic differentiation. We used flow sorted chromosomes from three different karyotypic taxa of the S. araneus group to map microsatellite markers at the chromosóme arm level. We have been able to map 24 markers and to show that the karyotypic organisation of these taxa is well conserved, which suggests that these markers can be used for further inter-taxa studies. A general prediction of chromosomal speciation models is that genetic differentiation between two taxa should be larger across rearranged chromosomes than across chromosomes common to both taxa. We combined two approaches using mapped microsatellites to test this prediction. First, we studied the genetic differentiation among five shrew taxa placed at different evolutionary levels (i.e. within and among species). In this large scale study, we detected an overall significant difference in genetic structure between rearranged vs. common chromosomes. Moreover, this effect varied among pairwise comparisons, which allowed us to differentiate the role of the karyotypic complexity of hybrids and of the evolutionary divergence between taxa. Secondly, we compared the levels of gene flow measured across common vs. rearranged chromosomes in two karyotypically different hybrid zones (strong vs. low complexity of hybrids), which show similar levels of genetic structure. We detected a significantly stronger genetic structure across rearranged chromosomes in the hybrid zone showing the highest level of hybrid complexity. The large variance observed among loci suggested that other factors, such as the position of markers within the chromosome, also certainly affects genetic structure. In conclusion, our results strongly support the role of chromosomal rearrangements in the reproductive barrier and suggest their importance in speciation process of the S. araneus group. RESUME Le rôle des réarrangements chromosomiques dans les processus de spéciation est fortement débattu et de nombreux modèles théoriques ont été développés sur le sujet. Les musaraignes du groupe Sorex araneus présentent de nombreuses opportunités pour étudier les relations entre les variations chromosomiques et la spéciation. En effet, ce groupe d'espèces morphologiquement très proches a attiré l'attention des chercheurs en raison de sa variabilité caryotypique principalement attribuée à des fusions Robertsoniennes. Pour explorer l'impact des changements caryotypiques sur la différenciation génétique, nous avons tout d'abord étudié les relations entre la structure génétique et caryotypique de races chromosomiques et d'espèces alpine du groupe S. araneus en utilisant des analyses Bayesiennes d' « admixture ». Les résultats de ces analyses ont confirmé le statut taxonomique des espèces étudiées bien que nous ayons détecté de l'introgression entre espèces. L'observation d'une sous structure spatiale relativement forte souligne l'importance des facteurs historiques (telle que l'isolation géographique) sur la structure génétique de ce groupe. Ensuite, nous avons étudié le flux de gène au niveau des chromosomes pour aborder de manière directe la question de l'impact des réarrangements chromosomiques sur la différenciation génétique. En conséquence, nous avons utilisé des tris de chromosomes de trois taxons du groupe S. araneus pour localiser des marqueurs microsatellites au niveau du bras chromosomique. Au cours de cette étude, nous avons pu localiser 24 marqueurs et montrer une forte conservation dans l'organisation du caryotype de ces taxa. Ce résultat suggère que leur utilisation est appropriée pour des études entre taxa. Une prédiction générale à tous les modèles de spéciation chromosomique correspond à la plus grande différenciation génétique des chromosomes réarrangés que des chromosomes communs. Nous avons combiné deux approches utilisant des microsatellites localisés au niveau du bras chromosomique pour tester cette prédiction. Premièrement, nous avons étudié la différenciation génétique entre cinq taxa du groupe S. araneus se trouvant à des niveaux évolutifs différents (i.e. à l'intérieur et entre espèce). Au cours de cette étude, nous avons détecté une différenciation globale significativement plus élevée sur les chromosomes réarrangés. Cet effet varie entre les comparaisons, ce qui nous a permis de souligner le rôle de la complexité caryotypique des hybrides et du niveau de divergence évolutive entre taxa. Deuxièmement, nous avons comparé le flux de gènes des chromosomes communs et réarrangés dans deux zones d'hybridation caryotypiquement différentes (forte vs. Faible complexité des hybrides) mais présentant un niveau de différenciation génétique similaire. Ceci nous a permis de détecter une structure génétique significativement plus élevée sur les chromosomes réarrangés au centre de la zone d'hybridation présentant la plus grande complexité caryotypic. La forte variance observée entre loci souligne en outre le fait que d'autres facteurs, tel que la position du marqueur sur le chromosome, affectent probablement aussi la structure génétique mesurée. En conclusion, nos résultats supportent fortement le rôle des réarrangements chromosomiques dans la barrière reproductive entre espèces ainsi que leur importance dans les processus de spéciation des musaraignes du groupe S. araneus.

Microsatellite variation reveals low genetic subdivision in a chromosome race of Sorex araneus (Mammalia, Insectivora)

Two hundred and forty-five individuals of the common shrew (Sorex araneus, Insectivora, Mammalia) from 24 sampling localities situated in four different valleys of the western European Alps were genotyped for six microsatellite loci. Allelic variability ranged from 3 to 32 different alleles at a single locus and the average gene diversity over all loci was 0.69. An analysis for F and R statistics revealed weak genetic population subdivision (Fst = 0.032; Rst = 0.016). This suggests considerable gene flow and little phylogeographic structure within and between valleys. We tested whether a stepwise mutation model (SMM) better explained variation at the microsatellite loci than an infinite allele model (IAM). No trend in favor of either model was detected.

Consequences of genetic erosion on fitness and phenotypic plasticity in European tree frog populations (Hyla arborea).

The detrimental effects of genetic erosion on small isolated populations are widely recognized contrary to their interactions with environmental changes. The ability of genotypes to plastically respond to variability is probably essential for the persistence of these populations. Genetic erosion impact may be exacerbated if inbreeding affects plastic responses or if their maintenance were at higher phenotypic costs. To understand the interplay 'genetic erosion-fitness-phenotypic plasticity', we experimentally compared, in different environments, the larval performances and plastic responses to predation of European tree frogs (Hyla arborea) from isolated and connected populations. Tadpoles from isolated populations were less performant, but the traits affected were environmental dependant. Heterosis observed in crosses between isolated populations allowed attributing their low fitness to inbreeding. Phenotypic plasticity can be maintained in the face of genetic erosion as inducible defences in response to predator were identical in all populations. However, the higher survival and developmental costs for isolated populations in harsh conditions may lead to an additional fitness loss for isolated populations.

Microsatellite markers reveal shallow genetic differentiation between cohorts of the common sea urchin Paracentrotus lividus (Lamarck) in northwest Mediterranean

Temporal variability was studied in the common sea urchin Paracentrotus lividus through the analysis of the genetic composition of three yearly cohorts sampled over two consecutive springs in a locality in northwestern Mediterranean. Individuals were aged using growth ring patterns observed in tests and samples were genotyped for five microsatellite loci. No reduction of genetic diversity was observed relative to a sample of the adult population from the same location or within cohorts across years. FST and amova results indicated that the differentiation between cohorts is rather shallow and not significant, as most variability is found within cohorts and within individuals. This mild differentiation translated into estimates of effective population size of 90100 individuals. When the observed excess of homozygotes was taken into account, the estimate of the average number of breeders increased to c. 300 individuals. Given our restricted sampling area and the known small-scale heterogeneity in recruitment in this species, our results suggest that at stretches of a few kilometres of shoreline, large numbers of progenitors are likely to contribute to the larval pool at each reproduction event. Intercohort variation in our samples is six times smaller than spatial variation between adults of four localities in the western Mediterranean. Our results indicate that, notwithstanding the stochastic events that take place during the long planktonic phase and during the settlement and recruitment processes, reproductive success in this species is high enough to produce cohorts genetically diverse and with little differentiation between them. Further research is needed before the link between genetic structure and underlying physical and biological processes can be well established.

Genetic polymorphisms and drug interactions modulating CYP2D6 and CYP3A activities have a major effect on oxycodone analgesic efficacy and safety

Background and purpose: The major drug-metabolizing enzymes for the oxidation of oxycodone are CYP2D6 and CYP3A. A high interindividual variability in the activity of these enzymes because of genetic polymorphisms and/or drug-drug interactions is well established. The possible role of an active metabolite in the pharmacodynamics of oxycodone has been questioned and the importance of CYP3A-mediated effects on the pharmacokinetics and pharmacodynamics of oxycodone has been poorly explored. Experimental approach: We conducted a randomized crossover (five arms) double-blind placebo-controlled study in 10 healthy volunteers genotyped for CYP2D6. Oral oxycodone (0.2 mg·kg−1) was given alone or after inhibition of CYP2D6 (with quinidine) and/or of CYP3A (with ketoconazole). Experimental pain (cold pressor test, electrical stimulation, thermode), pupil size, psychomotor effects and toxicity were assessed. Key results: CYP2D6 activity was correlated with oxycodone experimental pain assessment. CYP2D6 ultra-rapid metabolizers experienced increased pharmacodynamic effects, whereas cold pressor test and pupil size were unchanged in CYP2D6 poor metabolizers, relative to extensive metabolizers. CYP2D6 blockade reduced subjective pain threshold (SPT) for oxycodone by 30% and the response was similar to placebo. CYP3A4 blockade had a major effect on all pharmacodynamic assessments and SPT increased by 15%. Oxymorphone Cmax was correlated with SPT assessment (ρS= 0.7) and the only independent positive predictor of SPT. Side-effects were observed after CYP3A4 blockade and/or in CYP2D6 ultra-rapid metabolizers. Conclusions and implications: The modulation of CYP2D6 and CYP3A activities had clear effects on oxycodone pharmacodynamics and these effects were dependent on CYP2D6 genetic polymorphism.

Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity.

The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10(-8)) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits.

Variability of root traits in common bean genotypes at different levels of phosphorus supply and ontogenetic stages

Selection of common bean (Phaseolus vulgaris L.) cultivars with enhanced root growth would be a strategy for increasing P uptake and grain yield in tropical soils, but the strong plasticity of root traits may compromise their inclusion in breeding programs. The aim of this study was to evaluate the magnitude of the genotypic variability of root traits in common bean plants at two ontogenetic stages and two soil P levels. Twenty-four common bean genotypes, comprising the four growth habits that exist in the species and two wild genotypes, were grown in 4 kg pots at two levels of applied P (20 and 80 mg kg-1) and harvested at the stages of pod setting and early pod filling. Root area and root length were measured by digital image analysis. Significant genotype × P level and genotype × harvest interactions in analysis of variance indicate that the genotypic variation of root traits depended on soil nutrient availability and the stage at which evaluation was made. Genotypes differed for taproot mass, basal and lateral root mass, root area and root length at both P levels and growth stages; differences in specific root area and length were small. Genotypes with growth habits II (upright indeterminate) and III (prostrate indeterminate) showed better adaptation to limited P supply than genotypes of groups I (determinate) and IV (indeterminate climbing). Between the two harvests, genotypes of groups II and III increased the mass of basal and lateral roots by 40 and 50 %, respectively, whereas genotypes of groups I and IV by only 7 and 19 %. Values of the genotypic coefficient of determination, which estimates the proportion of phenotypic variance resulting from genetic effects, were higher at early pod filling than at pod setting. Correlations between shoot mass and root mass, which could indicate indirect selection of root systems via aboveground biomass, were higher at early pod filling than at pod setting. The results indicate that selection for root traits in common bean genotypes should preferentially be performed at the early pod-filling stage.

Genetic polymorphisms and drug interactions modulating CYP2D6 and CYP3A activities have a major effect on oxycodone analgesic efficacy and safety.

BACKGROUND AND PURPOSE: The major drug-metabolizing enzymes for the oxidation of oxycodone are CYP2D6 and CYP3A. A high interindividual variability in the activity of these enzymes because of genetic polymorphisms and/or drug-drug interactions is well established. The possible role of an active metabolite in the pharmacodynamics of oxycodone has been questioned and the importance of CYP3A-mediated effects on the pharmacokinetics and pharmacodynamics of oxycodone has been poorly explored. EXPERIMENTAL APPROACH: We conducted a randomized crossover (five arms) double-blind placebo-controlled study in 10 healthy volunteers genotyped for CYP2D6. Oral oxycodone (0.2 mg x kg(-1)) was given alone or after inhibition of CYP2D6 (with quinidine) and/or of CYP3A (with ketoconazole). Experimental pain (cold pressor test, electrical stimulation, thermode), pupil size, psychomotor effects and toxicity were assessed. KEY RESULTS: CYP2D6 activity was correlated with oxycodone experimental pain assessment. CYP2D6 ultra-rapid metabolizers experienced increased pharmacodynamic effects, whereas cold pressor test and pupil size were unchanged in CYP2D6 poor metabolizers, relative to extensive metabolizers. CYP2D6 blockade reduced subjective pain threshold (SPT) for oxycodone by 30% and the response was similar to placebo. CYP3A4 blockade had a major effect on all pharmacodynamic assessments and SPT increased by 15%. Oxymorphone C(max) was correlated with SPT assessment (rho(S)= 0.7) and the only independent positive predictor of SPT. Side-effects were observed after CYP3A4 blockade and/or in CYP2D6 ultra-rapid metabolizers. CONCLUSIONS AND IMPLICATIONS: The modulation of CYP2D6 and CYP3A activities had clear effects on oxycodone pharmacodynamics and these effects were dependent on CYP2D6 genetic polymorphism.

A comparison of three methods of Mendelian randomization when the genetic instrument, the risk factor and the outcome are all binary.

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Enzyme polymorphism and clinical variability of diseases: a study of diabetes mellitus.

We investigated possible relations among four common neonatal manifestations of diabetic pregnancy (macrosomia, hypoglycemia, hypocalcemia, jaundice) and four enzyme polymorphisms (PGM1, ADA, AK1, ACP1 in a sample of infants born of diabetic mothers. The pattern of associations observed between the two sets of variables is consistent with known differences in enzymatic activity within phenotypes of each system, suggesting that low enzymatic activity may have unfavorable effects on fetal development and on adaptability of the neonate to the extrauterine environment, Some of the polymorphic enzymes studied influence fetal growth in normal pregnancy as well. Analysis of relations between genetic polymorphisms and the clinical pattern of common diseases may provide a better understanding of the genetic basis of the clinical variability of diseases within and between human populations.

A genetic basis for functional hypothalamic amenorrhea.

Background: Functional hypothalamic amenorrhea is a reversible form of gonadotropin-releasing hormone (GnRH) deficiency commonly triggered by stressors such as excessive exercise, nutritional deficits, or psychological distress. Women vary in their susceptibility to inhibition of the reproductive axis by such stressors, but it is unknown whether this variability reflects a genetic predisposition to hypothalamic amenorrhea. We hypothesized that mutations in genes involved in idiopathic hypogonadotropic hypogonadism, a congenital form of GnRH deficiency, are associated with hypothalamic amenorrhea. Methods: We analyzed the coding sequence of genes associated with idiopathic hypogonadotropic hypogonadism in 55 women with hypothalamic amenorrhea and performed in vitro studies of the identified mutations. Results: Six heterozygous mutations were identified in 7 of the 55 patients with hypothalamic amenorrhea: two variants in the fibroblast growth factor receptor 1 gene FGFR1 (G260E and R756H), two in the prokineticin receptor 2 gene PROKR2 (R85H and L173R), one in the GnRH receptor gene GNRHR (R262Q), and one in the Kallmann syndrome 1 sequence gene KAL1 (V371I). No mutations were found in a cohort of 422 controls with normal menstrual cycles. In vitro studies showed that FGFR1 G260E, FGFR1 R756H, and PROKR2 R85H are loss-of-function mutations, as has been previously shown for PROKR2 L173R and GNRHR R262Q. Conclusions: Rare variants in genes associated with idiopathic hypogonadotropic hypogonadism are found in women with hypothalamic amenorrhea, suggesting that these mutations may contribute to the variable susceptibility of women to the functional changes in GnRH secretion that characterize hypothalamic amenorrhea. Our observations provide evidence for the role of rare variants in common multifactorial disease. (Funded by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and others; ClinicalTrials.gov number, NCT00494169.)

Etude clinique et analyse de liaison au locus 3q28 de deux familles suisses avec atrophie optique dominante de Kjer (OPA1) [Clinical study and genetic 3q28 locus linkage in 2 Swiss families with Kjer dominant optic atrophy (OPA1)]

METHODS: We examined 20 patients from 2 unrelated Swiss families to describe their clinical phenotype. In addition, a linkage analysis was performed in an attempt to confirm the reported genetic homogeneity of this condition as well as to refine its genomic localization. RESULTS: Two point analysis provided a cumulative LOD-score of 3.03 with marker D3S 2305. The absence of recombination precluded further refinement of the disease interval. CONCLUSIONS: Our data confirm the genetic homogeneity and the extreme variability of expression, occasionally mimicking low tension glaucoma.