Decades of population genetic research reveal the need for harmonization of molecular markers: the grey wolf Canis lupus as a case study

Following protection measures implemented since the 1970s, large carnivores are currently increasing in number and returning to areas from which they were absent for decades or even centuries. Monitoring programmes for these species rely extensively on non-invasive sampling and genotyping. However, attempts to connect results of such studies at larger spatial or temporal scales often suffer from the incompatibility of genetic markers implemented by researchers in different laboratories. This is particularly critical for long-distance dispersers, revealing the need for harmonized monitoring schemes that would enable the understanding of gene flow and dispersal dynamics. Based on a review of genetic studies on grey wolves Canis lupus from Europe, we provide an overview of the genetic markers currently in use, and identify opportunities and hurdles for studies based on continent-scale datasets. Our results highlight an urgent need for harmonization of methods to enable transnational research based on data that have already been collected, and to allow these data to be linked to material collected in the future. We suggest timely standardization of newly developed genotyping approaches, and propose that action is directed towards the establishment of shared single nucleotide polymorphism panels, next-generation sequencing of microsatellites, a common reference sample collection and an online database for data exchange. Enhanced cooperation among genetic researchers dealing with large carnivores in consortia would facilitate streamlining of methods, their faster and wider adoption, and production of results at the large spatial scales that ultimately matter for the conservation of these charismatic species.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree.

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics.

Integrating single nucleotide polymorphism (SNP) p-values from genome-wide association studies (GWAS) across genes and pathways is a strategy to improve statistical power and gain biological insight. Here, we present Pascal (Pathway scoring algorithm), a powerful tool for computing gene and pathway scores from SNP-phenotype association summary statistics. For gene score computation, we implemented analytic and efficient numerical solutions to calculate test statistics. We examined in particular the sum and the maximum of chi-squared statistics, which measure the strongest and the average association signals per gene, respectively. For pathway scoring, we use a modified Fisher method, which offers not only significant power improvement over more traditional enrichment strategies, but also eliminates the problem of arbitrary threshold selection inherent in any binary membership based pathway enrichment approach. We demonstrate the marked increase in power by analyzing summary statistics from dozens of large meta-studies for various traits. Our extensive testing indicates that our method not only excels in rigorous type I error control, but also results in more biologically meaningful discoveries.

Host Genetics of Response to Hepatitis B Vaccine: A Systematic Review

Background. Hepatitis B virus (HBV) is an important cause of chronic viral disease worldwide and can be life threatening. While a safe and effective vaccine is widely available, 5 to 10% of healthy vaccinees fail to achieve a protective anti-hepatitis B surface antigen antibody (anti-HBs) titer (>10mIU/ml). A limited number of studies investigated host genetics of the response to HBV vaccine. To our knowledge, no comprehensive overview of genetic polymorphisms both within and outside the HLA system has been done so far. Aim. The aim of this study was to perform a systematic review of the literature of human genetics influencing immune response after hepatitis B vaccination. Methods. Literature searches using keywords were conducted in the electronic databases Medline, Embase and ISI Web of Science the cut-off date being March 2014. After selection of papers according to stringent inclusion criteria, relevant information was systematically collected from the remaining articles, including demographic data, number of patients, schedule and type of vaccine, phenotypes, genes and single nucleotide polymorphisms (SNPs) genotyping results and their association with immune response to hepatitis B vaccine. Results. The literature search produced a total of 1968 articles from which 46 studies were kept for further analyses. From these studies, data was extracted for 19 alleles from the human leukocyte antigen (HLA) region that were reported as significant at least twice. Among those alleles, 9 were firmly associated with vaccine response outcome (DQ2 [DQB1*02 and DQB1*0201], DR3 [DRB1*03 and DRB1*0301], DR7 [DRB1*07 and DRB1*0701], C4AQ0, DPB1*0401, DQ3, DQB1*06, DRB1*01 and DRB1*13 [DRB1*1301]). In addition, data was extracted for 55 different genes from which 13 extra-HLA genes had polymorphisms that were studied by different group of investigators or by the same group with a replication study. Among the 13 genes allowing comparison, 4 genes (IL-1B, IL-2, IL-4R and IL- 6) revealed no significant data, 6 genes (IL-4, IL-10, IL-12B, IL-13, TNFA, IFNG and TLR2) were explored with inconsistent results and 2 genes (CD3Z and ITGAL) yielded promising results as their association with vaccine response was confirmed by a replication approach. Furthermore, this review produced a list of 46 SNPs from 26 genes that were associated with immune response to vaccine only once, providing novel candidates to be tested in datasets from existing genome-wide association studies (GWAS). Conclusion. To the best of our knowledge, this is the first systematic review of immunogenetic studies of response to hepatitis B vaccine. While this work reassesses the role of several HLA alleles on vaccine response outcome, the associations with polymorphisms in genes outside the HLA region were rather inconsistent. Moreover, this work produced a list of 46 significant SNPs that were reported by a single group of investigators, opening up some interesting possibilities for further research.

Genomic evidence for adaptive inversion clines in Drosophila melanogaster.

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Transgenic mice overexpressing neuropeptide Y: An experimental model of metabolic and cardiovascular diseases

Neuropeptide Y (NPY) is an abundant neurotransmitter in the brain and sympathetic nervous system (SNS). Hypothalamic NPY is known to be a key player in food intake and energy expenditure. NPY’s role in cardiovascular regulation has also been shown. In humans, a Leucine 7 to Proline 7 single nucleotide polymorphism (p.L7P) in the signal peptide of the NPY gene has been associated with traits of metabolic syndrome. The p.L7P subjects also show increased stress-related release of NPY, which suggests that more NPY is produced and released from SNS. The main objective of this study was to create a novel mouse model with noradrenergic cell-targeted overexpression of NPY, and to characterize the metabolic and vascular phenotype of this model. The mouse model was named OE-NPY^DBH mouse. Overexpression of NPY in SNS and brain noradrenergic neurons led to increased adiposity without significant weight gain or increased food intake. The mice showed lipid accumulation in the liver at young age, which together with adiposity led to impaired glucose tolerance and hyperinsulinemia with age. The mice displayed stress-related increased mean arterial blood pressure, increased plasma levels of catecholamines and enhanced SNS activity measured by GDP binding activity to brown adipose tissue mitochondria. Sexual dimorphism in NPY secretion pattern in response to stress was also seen. In an experimental model of vascular injury, the OE-NPY^DBH mice developed more pronounced neointima formation compared with wildtype controls. These results together with the clinical data indicate that NPY in noradrenergic cells plays an important role in the pathogenesis of metabolic syndrome and related diseases. Furthermore, new insights on the role of the extrahypothalamic NPY in the process have been obtained. The OE-NPY^DBH model provides an important tool for further stress and metabolic syndrome-related studies.

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8. Results Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion The 'isolated population' label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

Polygenic determinants of white matter volume derived from GWAS lack reproducibility in a replicate sample

A recent publication reported an exciting polygenic effect of schizophrenia (SCZ) risk variants, identified by a large genome-wide association study (GWAS), on total brain and white matter volumes in schizophrenic patients and, even more prominently, in healthy subjects. The aim of the present work was to replicate and then potentially extend these findings. According to the original publication, polygenic risk scores using single nucleotide polymorphism (SNP) information of SCZ GWAS (polygenic SCZ risk scores; PSS) were calculated in 122 healthy subjects, enrolled in a structural magnetic resonance imaging (MRI) study. These scores were computed based on P-values and odds ratios available through the Psychiatric GWAS Consortium. In addition, polygenic white matter scores (PWM) were calculated, using the respective SNP subset in the original publication. None of the polygenic scores, either PSS or PWM, were found to be associated with total brain, white matter or gray matter volume in our replicate sample. Minor differences between the original and the present study that might have contributed to lack of reproducibility (but unlikely explain it fully), are number of subjects, ethnicity, age distribution, array technology, SNP imputation quality and MRI scanner type. In contrast to the original publication, our results do not reveal the slightest signal of association of the described sets of GWAS-identified SCZ risk variants with brain volumes in adults. Caution is indicated in interpreting studies building on polygenic risk scores without replication sample.

Metodologia baseada em técnicas de mineração de dados para suporte à certificação de raças de ovinos

RESUMO O objetivo deste trabalho foi desenvolver uma metodologia baseada em técnicas de mineração de dados para selecionar os principais marcadores SNP (Single Nucleotide Polymorphism) para as raças de ovinos: Crioula, Morada Nova e Santa Inês. Os dados utilizados foram obtidos do Consórcio Internacional de Ovinos e são compostos por 72 animais das raças citadas, e cada animal possui 49.034 marcadores SNP. Considerando que o número de atributos (marcadores) é muito maior que o de observações (animais), foram aplicadas as técnicas de predição LASSO (Least Absolute Shrinkage and Selection Operator), Random Forest e Boosting para a geração de modelos preditivos que incorporam métodos de seleção de atributos. Os resultados revelaram que os modelos preditivos selecionaram os principais marcadores SNP para identificação das raças estudadas. O modelo LASSO selecionou um total de 29 marcadores relevantes. A partir dos modelos Random Forest e Boosting, foram obtidos 27 e 20 marcadores importantes, respectivamente. Por meio da intersecção dos modelos gerados, identificou-se um subconjunto de 18 marcadores com maior potencial de identificação das raças.

Alpha2-Adrenoceptor-mediated vascular responses induced by dexmedetomidine

Alpha2-Adrenoceptors are cell-surface G protein coupled receptors that mediate many of the effects of the catecholamines noradrenaline and adrenaline. The three human α2-adrenoceptor subtypes are widely expressed in different tissues and organs, and they mediate many different physiological and pharmacological effects in the central and peripheral nervous system and as postsynaptic receptors in target organs. Previous studies have demonstrated that α2-adrenoceptors mediate both vascular constriction and dilatation in humans. Large inter-individual variation has been observed in the vascular responses to α2-adrenoceptor activation in clinical studies. All three receptor subtypes are potential drug targets. It was therefore considered important to further elucidate the details of adrenergic vascular regulation and its genetic variation, since such knowledge may help to improve the development of future cardiovascular drugs and intensive care therapies. Dexmedetomidine is the most selective and potent α2-adrenoceptor agonist currently available for clinical use. When given systemically, dexmedetomidine induces nearly complete sympatholysis already at low concentrations, and postsynaptic effects, such vasoconstriction, can be observed with increasing concentrations. Thus, local infusions of small doses of dexmedetomidine into dorsal hand veins and the application of pharmacological sympathectomy with brachial plexus block provide a means to assess drug-induced peripheral vascular responses without interference from systemic pharmacological effects and autonomic nervous system regulation. Dexmedetomidine was observed to have biphasic effects on haemodynamics, with an initial decrease in blood pressure at low concentrations followed by substantial increases in blood pressure and coronary vascular resistance at high concentrations. Plasma concentrations of dexmedetomidine that significantly exceeded the recommended therapeutic level did not reduce myocardial blood flow below the level that is observed with the usual therapeutic concentrations and did not induce any evident myocardial ischaemia in healthy subjects. Further, it was demonstrated that dexmedetomidine also had significant vasodilatory effects through activation of endothelial nitric oxide synthesis, and thus when the endothelial component of the blood vessel response to dexmedetomidine was inhibited, peripheral vasoconstriction was augmented. Hand vein constriction responses to α2-adrenoceptor activation by dexmedetomidine were only weakly associated with the constriction responses to α1-adrenoceptor activation, pointing to independent cellular regulation by these two adrenoceptor classes. Substantial inter-individual variation was noted in the venous constriction elicited by activation of α2-adrenoceptors by dexmedetomidine. In two study populations from two different continents, a single nucleotide polymorphism in the PRKCB gene was found to be associated with the dorsal hand vein constriction response to dexmedetomidine, suggesting that protein kinase C beta may have an important role in the vascular α2-adrenoceptor signalling pathways activated by dexmedetomidine.

Prevalence study of SNP c.421G>T in the ADAMTS2 gene responsible for dermatosparaxis in White Dorper sheep in Brazil

Abstract: Dermatosparaxis is an autosomal recessive disorder of connective tissue; the disorder is clinically characterized by skin fragility and hyperextensibility. Dermatosparaxis in White Dorper sheep is caused by a single nucleotide polymorphism (SNP) (c.421G>T) in the ADAM metalloproteinase with thrombospondin type 1 motif, 2 (ADAMTS2) gene. The aim of this study was to investigate the prevalence of this SNP in a White Dorper herd in São Paulo state, Brazil. In this study, we collected blood DNA samples from 303 White Dorper sheep and performed polymerase chain reaction to amplify the SNP region. The samples were sequenced to determine the presence of the SNP in the ADAMTS2 gene. The SNP prevalence in the studied population was 15.5%; this finding indicates that more effective control measures should be used to prevent the inheritance of SNP c.421G>T in the ADAMTS2 gene in Brazilian White Dorper herds.

High baseline serum total and LDL cholesterol levels are associated with MDR1 haplotypes in Brazilian hypercholesterolemic individuals of European descent

The MDR1 gene encodes the P-glycoprotein, an efflux transporter with broad substrate specificity. P-glycoprotein has raised great interest in pharmacogenetics because it transports a variety of structurally divergent drugs, including lipid-lowering drugs. The synonymous single-nucleotide polymorphism C3435T and the nonsynonymous single-nucleotide polymorphism G2677T/A in MDR1 have been indicated as potential determinants of variability in drug disposition and efficacy. In order to evaluate the effect of G2677T/A and C3435T MDR1 polymorphisms on serum levels of lipids before and after atorvastatin administration, 69 unrelated hypercholesterolemic individuals from São Paulo city, Brazil, were selected and treated with 10 mg atorvastatin orally once daily for four weeks. MDR1 polymorphisms were analyzed by PCR-RFLP. C3435T and G2677T polymorphisms were found to be linked. The allelic frequencies for C3435T polymorphism were 0.536 and 0.464 for the 3435C and 3435T alleles, respectively, while for G2677T/A polymorphism allele frequencies were 0.580 for the 2677G allele, 0.384 for the 2677T allele and 0.036 for the 2677A allele. There was no significant relation between atorvastatin response and MDR1 polymorphisms (repeated measures ANOVA; P > 0.05). However, haplotype analysis revealed an association between T/T carriers and higher basal serum total (TC) and LDL cholesterol levels (TC: 303 ± 56, LDL-C: 216 ± 57 mg/dl, respectively) compared with non-T/T carriers (TC: 278 ± 28, LDL-C: 189 ± 24 mg/dl; repeated measures ANOVA/Tukey test; P < 0.05). These data indicate that MDR1 polymorphism may have an important contribution to the control of basal serum cholesterol levels in Brazilian hypercholesterolemic individuals of European descent.

The search for circadian clock components in humans: new perspectives for association studies

Individual circadian clocks entrain differently to environmental cycles (zeitgebers, e.g., light and darkness), earlier or later within the day, leading to different chronotypes. In human populations, the distribution of chronotypes forms a bell-shaped curve, with the extreme early and late types _ larks and owls, respectively _ at its ends. Human chronotype, which can be assessed by the timing of an individual's sleep-wake cycle, is partly influenced by genetic factors - known from animal experimentation. Here, we review population genetic studies which have used a questionnaire probing individual daily timing preference for associations with polymorphisms in clock genes. We discuss their inherent limitations and suggest an alternative approach combining a short questionnaire (Munich ChronoType Questionnaire, MCTQ), which assesses chronotype in a quantitative manner, with a genome-wide analysis (GWA). The advantages of these methods in comparison to assessing time-of-day preferences and single nucleotide polymorphism genotyping are discussed. In the future, global studies of chronotype using the MCTQ and GWA may also contribute to understanding the influence of seasons, latitude (e.g., different photoperiods), and climate on allele frequencies and chronotype distribution in different populations.

Simultaneous detection of the C282Y, H63D and S65C mutations in the hemochromatosis gene using quenched-FRET real-time PCR

Hereditary hemochromatosis (HH) is a common autosomal disorder of iron metabolism mainly affecting Caucasian populations. Three recurrent disease-associated mutations have been detected in the hemochromatosis gene (HFE): C282Y, H63D, and S65C. Although HH phenotype has been associated with all three mutations, C282Y is considered the most relevant mutation responsible for hemochromatosis. Clinical complications of HH include cirrhosis of the liver, congestive cardiac failure and cardiac arrhythmias, endocrine pancreatic disease, which can be prevented by early diagnosis and treatment. Therefore, a reliable genotyping method is required for presymptomatic diagnosis. We describe the simultaneous detection of the C282Y, H63D and S65C mutations in the hemochromatosis gene by real-time PCR followed by melting curve analysis using fluorescence resonance energy transfer (FRET) probes. The acceptor fluorophore may be replaced by a quencher, increasing multiplex possibilities. Real-time PCR results were compared to the results of sequencing and conventional PCR followed by restriction digestion and detection by agarose gel electrophoresis (PCR-RFLP). Genotypes from 80 individuals obtained both by the conventional PCR-RFLP method and quenched-FRET real-time PCR were in full agreement. Sequencing also confirmed the results obtained by the new method, which proved to be an accurate, rapid and cost-effective diagnostic assay. Our findings demonstrate the usefulness of real-time PCR for the simultaneous detection of mutations in the HFE gene, which allows a reduction of a significant amount of time in sample processing compared to the PCR-RFLP method, eliminates the use of toxic reagents, reduces the risk of contamination in the laboratory, and enables full process automation.

Genetics of homocysteine metabolism and associated disorders

Homocysteine is a sulfur-containing amino acid derived from the metabolism of methionine, an essential amino acid, and is metabolized by one of two pathways: remethylation or transsulfuration. Abnormalities of these pathways lead to hyperhomocysteinemia. Hyperhomocysteinemia is observed in approximately 5% of the general population and is associated with an increased risk for many disorders, including vascular and neurodegenerative diseases, autoimmune disorders, birth defects, diabetes, renal disease, osteoporosis, neuropsychiatric disorders, and cancer. We review here the correlation between homocysteine metabolism and the disorders described above with genetic variants on genes coding for enzymes of homocysteine metabolism relevant to clinical practice, especially common variants of the MTHFR gene, 677C>T and 1298A>C. We also discuss the management of hyperhomocysteinemia with folic acid supplementation and fortification of folic acid and the impact of a decrease in the prevalence of congenital anomalies and a decline in the incidence of stroke mortality.