quantiNemo: an individual-based program to simulate quantitative traits with explicit genetic architecture in a dynamic metapopulation.

quantiNemo is an individual-based, genetically explicit stochastic simulation program. It was developed to investigate the effects of selection, mutation, recombination and drift on quantitative traits with varying architectures in structured populations connected by migration and located in a heterogeneous habitat. quantiNemo is highly flexible at various levels: population, selection, trait(s) architecture, genetic map for QTL and/or markers, environment, demography, mating system, etc. quantiNemo is coded in C++ using an object-oriented approach and runs on any computer platform. Availability: Executables for several platforms, user's manual, and source code are freely available under the GNU General Public License at http://www2.unil.ch/popgen/softwares/quantinemo.

Mutations in the cdc10 start gene of Schizosaccharomyces pombe implicate the region of homology between cdc10 and SWI6 as important for p85cdc10 function.

The cdc10 gene of the fission yeast Schizosaccharomyces pombe is required for traverse of start and commitment to the mitotic cell division cycle rather than other fates. The product of the gene, p85cdc10, is a component of a factor that is thought to be involved in regulating the transcription of genes that are required for DNA synthesis. In order to define regions of the p85cdc10 protein that are important for its function a fine structure genetic map of the cdc10 gene was derived and the sequences of 13 cdc10ts mutants determined. The 13 mutants tested define eight alleles. Eleven of the mutants are located in the region that contains the two copies of the cdc10/SWI6 repeat motif, implicating it as important for p85cdc10 function.

Comparative genetic analyses point to HCP5 as susceptibility locus for HCV-associated hepatocellular carcinoma.

BACKGROUND & AIMS: Recently, genetic variations in MICA (lead single nucleotide polymorphism [SNP] rs2596542) were identified by a genome-wide association study (GWAS) to be associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) in Japanese patients. In the present study, we sought to determine whether this SNP is predictive of HCC development in the Caucasian population as well. METHODS: An extended region around rs2596542 was genotyped in 1924 HCV-infected patients from the Swiss Hepatitis C Cohort Study (SCCS). Pair-wise correlation between key SNPs was calculated both in the Japanese and European populations (HapMap3: CEU and JPT). RESULTS: To our surprise, the minor allele A of rs2596542 in proximity of MICA appeared to have a protective impact on HCC development in Caucasians, which represents an inverse association as compared to the one observed in the Japanese population. Detailed fine-mapping analyses revealed a new SNP in HCP5 (rs2244546) upstream of MICA as strong predictor of HCV-related HCC in the SCCS (univariable p=0.027; multivariable p=0.0002, odds ratio=3.96, 95% confidence interval=1.90-8.27). This newly identified SNP had a similarly directed effect on HCC in both Caucasian and Japanese populations, suggesting that rs2244546 may better tag a putative true variant than the originally identified SNPs. CONCLUSIONS: Our data confirms the MICA/HCP5 region as susceptibility locus for HCV-related HCC and identifies rs2244546 in HCP5 as a novel tagging SNP. In addition, our data exemplify the need for conducting meta-analyses of cohorts of different ethnicities in order to fine map GWAS signals.

Fas ligand-induced proinflammatory transcriptional responses in reconstructed human epidermis. Recruitment of the epidermal growth factor receptor and activation of MAP kinases.

Fas ligand (FasL) exerts potent proapoptotic and proinflammatory actions on epidermal keratinocytes and has been implicated in the pathogenesis of eczema, toxic epidermal necrolysis, and drug-induced skin eruptions. We used reconstructed human epidermis to investigate the mechanisms of FasL-induced inflammatory responses and their relationships with FasL-triggered caspase activity. Caspase activity was a potent antagonist of the pro-inflammatory gene expression triggered by FasL prior to the onset of cell death. Furthermore, we found that FasL-stimulated autocrine production of epidermal growth factor receptor (EGFR) ligands, and the subsequent activation of EGFR and ERK1 and ERK2 mitogen-activated protein kinases, were obligatory extracellular steps for the FasL-induced expression of a subset of inflammatory mediators, including CXCL8/interleukin (IL)-8, ICAM-1, IL-1alpha, IL-1beta, CCL20/MIP-3alpha, and thymic stromal lymphopoietin. These results expand the known physiological role of EGFR and its ligands from promoting keratinocyte mitogenesis and survival to mediating FasL-induced epidermal inflammation.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations.

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.

A novel vasopressin-induced transcript promotes MAP kinase activation and ENaC downregulation.

In the principal cell of the renal collecting duct, vasopressin regulates the expression of a gene network responsible for sodium and water reabsorption through the regulation of the water channel and the epithelial sodium channel (ENaC). We have recently identified a novel vasopressin-induced transcript (VIT32) that encodes for a 142 amino acid vasopressin-induced protein (VIP32), which has no homology with any protein of known function. The Xenopus oocyte expression system revealed two functions: (i) when injected alone, VIT32 cRNA rapidly induces oocyte meiotic maturation through the activation of the maturation promoting factor, the amphibian homolog of the universal M phase trigger Cdc2/cyclin; and (ii) when co-injected with the ENaC, VIT32 cRNA selectively downregulates channel activity, but not channel cell surface expression. In the kidney principal cell, VIP32 may be involved in the downregulation of transepithelial sodium transport observed within a few hours after vasopressin treatment. VIP32 belongs to a novel gene family ubiquitously expressed in oocyte and somatic cells that may be involved in G to M transition and cell cycling.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.

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.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species.

BACKGROUND: Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications. FINDINGS: We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement. CONCLUSIONS: The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

MAP/ERK kinase kinase 1 (MEKK1) mediates transcriptional repression by interacting with polycystic kidney disease-1 (PKD1) promoter-bound p53 tumor suppressor protein.

Mitogen-activated protein kinase (MAPK) cascades regulate a wide variety of cellular processes that ultimately depend on changes in gene expression. We have found a novel mechanism whereby one of the key MAP3 kinases, Mekk1, regulates transcriptional activity through an interaction with p53. The tumor suppressor protein p53 down-regulates a number of genes, including the gene most frequently mutated in autosomal dominant polycystic kidney disease (PKD1). We have discovered that Mekk1 translocates to the nucleus and acts as a co-repressor with p53 to down-regulate PKD1 transcriptional activity. This repression does not require Mekk1 kinase activity, excluding the need for an Mekk1 phosphorylation cascade. However, this PKD1 repression can also be induced by the stress-pathway stimuli, including TNFα, suggesting that Mekk1 activation induces both JNK-dependent and JNK-independent pathways that target the PKD1 gene. An Mekk1-p53 interaction at the PKD1 promoter suggests a new mechanism by which abnormally elevated stress-pathway stimuli might directly down-regulate the PKD1 gene, possibly causing haploinsufficiency and cyst formation.

High-density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family.

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well-differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG2 ), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin 21 serum level, hemoglobin drop and therapeutic response in patients with HCV infection

Background: T reatment o f chronic hepatitis C i s evolving, a nd direct acting antivirals ( DAAs) are now a dded to p egylated interferon-α ( Peg- INF-α) and ribavirin (RBV) for the treatment o f hepatitis C v irus ( HCV) genotype 1 infection. DAAs c ause d ifferent side effects and can even worsen RBV induced hemolytic anemia. T herefore, identifying host genetic d eterminants of R BV bioavailability and therapeutic e fficacy will remain crucial for individualized treatment. Recent d ata showed associations between R BV induced h emolytic anemia and genetic polymorphisms o f concentrative nucleoside transporters s uch as C NT3 (SLC28A3) and i nosine t riphosphatase (ITPA). T o analyze t he association of genetic variants of SLC28 transporters and ITPA with RBV induced hemolytic anemia and treatment o utcome. Methods: I n our study, 173 patients f rom t he S wiss Hepatitis C C ohort Study and 2 2 patients from Swiss Association for the Study of the Liver study 24 (61% HCV g enotype 1, 3 9% genotypes 2 o r 3) were analyzed for SLC28A2 single nucleotide p olymorphism (SNP) rs11854484, SLC28A3 rs56350726 and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101. RBV serum levels during treatment were measured in 49 patients. Results: SLC28A2 r s11854484 genotype TT was associated with significantly higher dosage- and body weight-adjusted RBV levels as compared to genotypes TC and CC (p=0.04 and p=0.02 at weeks 4 and 8, respectively). ITPA SNPs rs1127354 and rs7270101 were associated with h emolytic a nemia both in genotype as w ell as i n allelic a nalyses. SLC28A3 rs56350726 genotype TT (vs. AT/AA, RR=2.1; 95% CI 1.1-4.1) as well as the T allele (vs. A; RR=1.8, 95% CI 1.1-3.2) were associated with increased SVR rates. The combined analysis of overall ITPA activity and SLC28 v ariants together revealed n o significant a dditive effects on either treatment-related anemia or SVR. Conclusions: T he newly identified association between RBV serum levels a nd SLC28A2 rs11854484 genotype as well as the replicated association of ITPA and SLC28A3 g enetic p olymorphisms w ith RBV induced hemolytic anemia and treatment r esponse underpin the need for further studies on host genetic d eterminants of R BV bioavailability and therapeutic e fficacy f or individualized treatment of chronic hepatitis C.

Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study.

AIMS/HYPOTHESIS: Several susceptibility genes for type 2 diabetes have been discovered recently. Individually, these genes increase the disease risk only minimally. The goals of the present study were to determine, at the population level, the risk of diabetes in individuals who carry risk alleles within several susceptibility genes for the disease and the added value of this genetic information over the clinical predictors. METHODS: We constructed an additive genetic score using the most replicated single-nucleotide polymorphisms (SNPs) within 15 type 2 diabetes-susceptibility genes, weighting each SNP with its reported effect. We tested this score in the extensively phenotyped population-based cross-sectional CoLaus Study in Lausanne, Switzerland (n = 5,360), involving 356 diabetic individuals. RESULTS: The clinical predictors of prevalent diabetes were age, BMI, family history of diabetes, WHR, and triacylglycerol/HDL-cholesterol ratio. After adjustment for these variables, the risk of diabetes was 2.7 (95% CI 1.8-4.0, p = 0.000006) for individuals with a genetic score within the top quintile, compared with the bottom quintile. Adding the genetic score to the clinical covariates improved the area under the receiver operating characteristic curve slightly (from 0.86 to 0.87), yet significantly (p = 0.002). BMI was similar in these two extreme quintiles. CONCLUSIONS/INTERPRETATION: In this population, a simple weighted 15 SNP-based genetic score provides additional information over clinical predictors of prevalent diabetes. At this stage, however, the clinical benefit of this genetic information is limited.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.