Graphical tests for Hardy-Weinberg equilibrium based on the ternary plot

A reduced SNP panel optimised for non-invasive genetic assessment of a genetically impoverished conservation icon, the European bison

The European bison was saved from the brink of extinction due to considerable conservation efforts since the early twentieth century. The current global population of > 9500 individuals is the result of successful ex situ breeding based on a stock of only 12 founders, resulting in an extremely low level of genetic variability. Due to the low allelic diversity, traditional molecular tools, such as microsatellites, fail to provide sufficient resolution for accurate genetic assessments in European bison, let alone from non-invasive samples. Here, we present a SNP panel for accurate high-resolution genotyping of European bison, which is suitable for a wide variety of sample types. The panel accommodates 96 markers allowing for individual and parental assignment, sex determination, breeding line discrimination, and cross-species detection. Two applications were shown to be utilisable in further Bos species with potential conservation significance. The new SNP panel will allow to tackle crucial tasks in European bison conservation, including the genetic monitoring of reintroduced populations, and a molecular assessment of pedigree data documented in the world's first studbook of a threatened species.

Accurate phenotypic classification and exome sequencing allow identification of novel genes and variants associated with adult-onset hearing loss

Adult-onset progressive hearing loss is a common, complex disease with a strong genetic component. Although to date over 150 genes have been identified as contributing to human hearing loss, many more remain to be discovered, as does most of the underlying genetic diversity. Many different variants have been found to underlie adult-onset hearing loss, but they tend to be rare variants with a high impact upon the gene product. It is likely that combinations of more common, lower impact variants also play a role in the prevalence of the disease. Here we present our exome study of hearing loss in a cohort of 532 older adult volunteers with extensive phenotypic data, including 99 older adults with normal hearing, an important control set. Firstly, we carried out an outlier analysis to identify genes with a high variant load in older adults with hearing loss compared to those with normal hearing. Secondly, we used audiometric threshold data to identify individual variants which appear to contribute to different threshold values. We followed up these analyses in a second cohort. Using these approaches, we identified genes and variants linked to better hearing as well as those linked to worse hearing. These analyses identified some known deafness genes, demonstrating proof of principle of our approach. However, most of the candidate genes are novel associations with hearing loss. While the results support the suggestion that genes responsible for severe deafness may also be involved in milder hearing loss, they also suggest that there are many more genes involved in hearing which remain to be identified. Our candidate gene lists may provide useful starting points for improved diagnosis and drug development.

Insecticide resistance levels and associated mechanisms in three Aedes aegypti populations from Venezuela

BACKGROUND

The massive use of insecticides in public health has exerted selective pressure resulting in the development of resistance in Aedes aegypti to different insecticides in Venezuela. Between 2010 and 2020, the only insecticides available for vector control were the organophosphates (Ops) fenitrothion and temephos which were focally applied.

OBJECTIVES

To determine the state of insecticide resistance and to identify the possible biochemical and molecular mechanisms involved in three populations of Ae. aegypti from Venezuela.

METHODS

CDC bottle bioassays were conducted on Ae. aegypti collected between October 2019 and February 2020 in two hyperendemic localities for dengue in Aragua State and in a malaria endemic area in Bolívar State. Insecticide resistance mechanisms were studied using biochemical assays and polymerase chain reaction (PCR) to detect kdr mutations.

FINDINGS

Bioassays showed contrasting results among populations; Las Brisas was resistant to malathion, permethrin and deltamethrin, Urbanización 19 de Abril was resistant to permethrin and Nacupay to malathion. All populations showed significantly higher activity of mixed function oxidases and glutathione-S-transferases (GSTs) in comparison with the susceptible strain. The kdr mutations V410L, F1534C, and V1016I were detected in all populations, with F1534C at higher frequencies.

MAIN CONCLUSION

Insecticide resistance persists in three Ae. aegypti populations from Venezuela even in the relative absence of insecticide application.

Tracing the introduction of the invasive common myna using population genomics

The common myna (Acridotheres tristis) is one of the most invasive bird species in the world, yet its colonisation history is only partly understood. We identified the introduction history and population structure, and quantified the genetic diversity of myna populations from the native range in India and introduced populations in New Zealand, Australia, Fiji, Hawaii, and South Africa, based on thousands of single nucleotide polymorphism markers in 814 individuals. We were able to identify the source population of mynas in several invasive locations: mynas from Fiji and Melbourne, Australia, were likely founded by individuals from a subpopulation in Maharashtra, India, while mynas in Hawaii and South Africa were likely independently founded by individuals from other localities in India. Our findings suggest that New Zealand mynas were founded by individuals from Melbourne, which, in turn, were founded by individuals from Maharashtra. We identified two genetic clusters among New Zealand mynas, divided by New Zealand's North Island's axial mountain ranges, confirming previous observations that mountains and thick forests may form barriers to myna dispersal. Our study provides a foundation for other population and invasion genomic studies and provides useful information for the management of this invasive species.

Accurate phenotypic classification and exome sequencing allow identification of novel genes and variants associated with adult-onset hearing loss

Adult-onset progressive hearing loss is a common, complex disease with a strong genetic component. Although to date over 150 genes have been identified as contributing to human hearing loss, many more remain to be discovered, as does most of the underlying genetic diversity. Many different variants have been found to underlie adult-onset hearing loss, but they tend to be rare variants with a high impact upon the gene product. It is likely that combinations of more common, lower impact variants also play a role in the prevalence of the disease. Here we present our exome study of hearing loss in a cohort of 532 older adult volunteers with extensive phenotypic data, including 99 older adults with normal hearing, an important control set. Firstly, we carried out an outlier analysis to identify genes with a high variant load in older adults with hearing loss compared to those with normal hearing. Secondly, we used audiometric threshold data to identify individual variants which appear to contribute to different threshold values. We followed up these analyses in a second cohort. Using these approaches, we identified genes and variants linked to better hearing as well as those linked to worse hearing. These analyses identified some known deafness genes, demonstrating proof of principle of our approach. However, most of the candidate genes are novel associations with hearing loss. While the results support the suggestion that genes responsible for severe deafness may also be involved in milder hearing loss, they also suggest that there are many more genes involved in hearing which remain to be identified. Our candidate gene lists may provide useful starting points for improved diagnosis and drug development.

Temporal and biting dynamics of the chromosomal inversion 2La in the malaria vectors Anopheles gambiae and Anopheles coluzzii in Bangui, Central African Republic

The chromosomal rearrangement 2La has been directly involved in the ecological and deadly epidemiological success of the malaria mosquitoes Anopheles gambiae and Anopheles coluzzii in sub-Saharan Africa. However, little is known about the biological and ecological factors that drive the local and temporal dynamics of this inversion in both species. Here, we performed a year-round longitudinal survey in Bangui, Central African Republic. We monthly sampled A. gambiae and A. coluzzii mosquitoes indoor and outdoor using human landing catches (HLC) for 48 h non-stop. We molecularly karyotyped all specimens to study the 2La inversion frequency variations, and monitored the mosquito spatial and temporal biting behavior throughout the year. In total, we successfully karyotyped 5121 A. gambiae and 986 A. coluzzii specimens. The 2La inversion frequency was higher in A. coluzzii than in A. gambiae across the year. In A. gambiae and A. coluzzii, the inversion frequency or karyotypes did not influence the biting behavior, either location or time. Moreover, the inversion frequency variation in both species was also independent of local climatic changes. Overall, our results revealed that in Bangui, the 2La inversion segregates at different frequency in each species, but this is not influenced by their trophic behavior. Studying the impact of urban settings and the population genetic structure of these two A. gambiae complex members could bring insights into the intrinsic relationship between 2La inversion and local conditions. More studies are needed to understand the polymorphic equilibrium of this inversion in Bangui.

Genomic basis of insularity and ecological divergence in barn owls (Tyto alba) of the Canary Islands

Islands, and the particular organisms that populate them, have long fascinated biologists. Due to their isolation, islands offer unique opportunities to study the effect of neutral and adaptive mechanisms in determining genomic and phenotypical divergence. In the Canary Islands, an archipelago rich in endemics, the barn owl (Tyto alba), present in all the islands, is thought to have diverged into a subspecies (T. a. gracilirostris) on the eastern ones, Fuerteventura and Lanzarote. Taking advantage of 40 whole-genomes and modern population genomics tools, we provide the first look at the origin and genetic makeup of barn owls of this archipelago. We show that the Canaries hold diverse, long-standing and monophyletic populations with a neat distinction of gene pools from the different islands. Using a new method, less sensitive to structure than classical F_ST, to detect regions involved in local adaptation to insular environments, we identified a haplotype-like region likely under selection in all Canaries individuals and genes in this region suggest morphological adaptations to insularity. In the eastern islands, where the subspecies is present, genomic traces of selection pinpoint signs of adapted body proportions and blood pressure, consistent with the smaller size of this population living in a hot arid climate. In turn, genomic regions under selection in the western barn owls from Tenerife showed an enrichment in genes linked to hypoxia, a potential response to inhabiting a small island with a marked altitudinal gradient. Our results illustrate the interplay of neutral and adaptive forces in shaping divergence and early onset speciation.

Interaction Effect between Physical Activity and the BDNF Val66Met Polymorphism on Depression in Women from the PISMA-ep Study

The relationship between depression and the Val66Met polymorphism at the brain-derived neurotrophic factor gene (BDNF), has been largely studied. It has also been related to physical activity, although the results remain inconclusive. The aim of this study is to investigate the relationship between this polymorphism, depression and physical activity in a thoroughly characterised sample of community-based individuals from the PISMA-ep study. A total of 3123 participants from the PISMA-ep study were genotyped for the BDNF Val66Met polymorphism, of which 209 had depression. Our results are in line with previous studies reporting a protective effect of physical activity on depression, specifically in light intensity. Interestingly, we report a gene-environment interaction effect in which Met allele carriers of the BDNF Val66Met polymorphism who reported more hours of physical activity showed a decreased prevalence of depression. This effect was observed in the total sample (OR = 0.95, 95%CI = 0.90–0.99, p = 0.027) and was strengthened in women (OR = 0.93, 95%CI = 0.87–0.98, p = 0.019). These results highlight the potential role of physical activity as a promising therapeutic strategy for preventing and adjuvant treatment of depression and suggest molecular and genetic particularities of depression between sexes.

Landscape and climatic variations shaped secondary contacts amid barn owls of the Western Palearctic

The combined actions of climatic variations and landscape barriers shape the history of natural populations. When organisms follow their shifting niches, obstacles in the landscape can lead to the splitting of populations, on which evolution will then act independently. When two such populations are reunited, secondary contact occurs in a broad range of admixture patterns, from narrow hybrid zones to the complete dissolution of lineages. A previous study suggested that barn owls colonized the Western Palearctic after the last glaciation in a ring-like fashion around the Mediterranean Sea, and conjectured an admixture zone in the Balkans. Here, we take advantage of whole-genome sequences of 94 individuals across the Western Palearctic to reveal the complex history of the species in the region using observational and modeling approaches. Even though our results confirm that two distinct lineages colonized the region, one in Europe and one in the Levant, they suggest that it predates the last glaciation and identify a secondary contact zone between the two in Anatolia. We also show that barn owls recolonized Europe after the glaciation from two distinct glacial refugia: a previously identified western one in Iberia and a new eastern one in Italy. Both glacial lineages now communicate via eastern Europe, in a wide and permeable contact zone. This complex history of populations enlightens the taxonomy of Tyto alba in the region, highlights the key role played by mountain ranges and large water bodies as barriers and illustrates the power of population genomics in uncovering intricate demographic patterns.

TIGER: The gene expression regulatory variation landscape of human pancreatic islets

Genome-wide association studies (GWASs) identified hundreds of signals associated with type 2 diabetes (T2D). To gain insight into their underlying molecular mechanisms, we have created the translational human pancreatic islet genotype tissue-expression resource (TIGER), aggregating >500 human islet genomic datasets from five cohorts in the Horizon 2020 consortium T2DSystems. We impute genotypes using four reference panels and meta-analyze cohorts to improve the coverage of expression quantitative trait loci (eQTL) and develop a method to combine allele-specific expression across samples (cASE). We identify >1 million islet eQTLs, 53 of which colocalize with T2D signals. Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identify eight cASE colocalizations, among which we found a T2D-associated SLC30A8 variant. We make all data available through the TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translates into therapeutic insight and precision medicine for T2D.

Genetic Differentiation and Demographic Trajectory of the Insular Formosan and Orii's Flying Foxes

Insular flying foxes are keystone species in island ecosystems due to their critical roles in plant pollination and seed dispersal. These species are vulnerable to population decline because of their small populations and low reproductive rates. The Formosan flying fox (Pteropus dasymallus formosus) is one of the 5 subspecies of the Ryukyu flying fox. Pteropus dasymallus formosus has suffered from a severe decline and is currently recognized as a critically endangered population in Taiwan. On the contrary, the Orii's flying fox (Pteropus dasymallus inopinatus) is a relatively stable population inhabiting Okinawa Island. Here, we applied a genomic approach called double digest restriction-site associated DNA sequencing to study these 2 subspecies for a total of 7 individuals. We detected significant genetic structure between the 2 populations. Despite their contrasting contemporary population sizes, both populations harbor very low degrees of genetic diversity. We further inferred their demographic history based on the joint folded site frequency spectrum and revealed that both P. d. formosus and P. d. inopinatus had maintained small population sizes for a long period of time after their divergence. Recently, these populations experienced distinct trajectories of demographic changes. While P. d. formosus suffered from a drastic ~10-fold population decline not long ago, P. d. inopinatus underwent a ~4.5-fold population expansion. Our results suggest separate conservation management for the 2 populations-population recovery is urgently needed for P. d. formosus while long-term monitoring for adverse genetic effects should be considered for P. d. inopinatus.

Validation of a prostate cancer polygenic risk score

BACKGROUND

Genome-wide association studies have identified over 100 single-nucleotide polymorphisms (SNPs) associated with prostate cancer (PrCa), and polygenic risk scores (PRS) based on their combined genotypes have been developed for risk stratification. We aimed to assess the contribution of PRS to PrCa risk in a large multisite study.

METHODS

The sample included 1972 PrCa cases and 1919 unaffected controls. Next-generation sequencing was used to assess pathogenic variants in 14 PrCa-susceptibility genes and 72 validated PrCa-associated SNPs. We constructed a population-standardized PRS and tested its association with PrCa using logistic regression adjusted for age and family history of PrCa.

RESULTS

The mean age of PrCa cases at diagnosis and age of controls at testing/last clinic visit was 59.5 ± 7.2 and 57.2 ± 13.0 years, respectively. Among 1740 cases with pathology data, 57.4% had Gleason score ≤ 6, while 42.6% had Gleason score ≥ 8. In addition, 39.6% cases and 20.1% controls had a family history of PrCa. The PRS was significantly higher in cases than controls (mean ± SD: 1.42 ± 1.11 vs 1.02 ± 0.76; P < .0001). Compared with men in the 1st quartile of age-adjusted PRS, those in the 2nd, 3rd, and 4th quartile were 1.58 (95% confidence interval [CI]: 1.31-1.90), 2.36 (95% CI: 1.96-2.84), and 3.98 (95% CI: 3.29-4.82) times as likely to have PrCa (all P < .0001). Adjustment for family history yielded similar results. PRS predictive performance was consistent with prior literature (area under the receiver operating curve = 0.64; 95% CI: 0.62-0.66).

CONCLUSIONS

These data suggest that a 72-SNP PRS is predictive of PrCa, supporting its potential use in clinical risk assessment.

Contemporary genetic structure affects genetic stock identification of steelhead trout in the Snake River basin

Genetic stock identification is a widely applied tool for the mixed‐stock management of salmonid species throughout the North Pacific Rim. The effectiveness of genetic stock identification is dependent on the level of differentiation among stocks which is often high due to the life history of these species that involves high homing fidelity to their natal streams. However, the utility of this tool can be reduced when natural genetic structuring has been altered by hatchery translocation and/or supplementation. We examined the genetic population structure of ESA‐listed steelhead in the Snake River basin of the United States. We analyzed 9,613 natural‐origin adult steelhead returning to Passive Integrated Transponder detection sites throughout the basin from 2010 through 2017. Individuals were genotyped at 180 single nucleotide polymorphic genetic markers and grouped into 20 populations based on their return location. While we expected to observe a common pattern of hierarchical genetic structuring due to isolation by distance, we observed low genetic differentiation between populations in the upper Salmon River basin compared to geographically distant populations in the lower Snake River basin. These results were consistent with lower genetic stock assignment probabilities observed for populations in this upper basin. We attribute these patterns of reduced genetic structure to the translocation of lower basin steelhead stocks and ongoing hatchery programs in the upper Salmon River basin. We discuss the implications of these findings on the utility of genetic stock identification in the basin and discuss opportunities for increasing assignment probabilities in the face of low genetic structure.

Comparative Analysis of Soil Microbiome Profiles in the Companion Planting of White Clover and Orchard Grass Using 16S rRNA Gene Sequencing Data

Companion planting is one of the most common and effective planting methods in modern agriculture. White clover (Trifolium repens L.) and orchard grass (Dactylis glomerata L.) are two typical pastures planted together to promote each other's growth. However, the detailed biological foundations of companion planting remain unclear. In this study, we screened typical microbiome profiles under separate and combination planting conditions using 16s RNA gene sequencing techniques. We identified the typical distinctive microorganism subtypes based on the microbiome profiles and recognized the enriched functions of top abundant microorganisms in soil using different planting strategies with the help of Kyoto Encyclopedia of Genes and Genomes and Clusters of Orthologous Groups annotation. This analysis confirmed that the optimal microorganisms and screened functional annotations are correlated with nitrogen fixation; thus, companion planting may improve the yield and efficacy of plants by improving the efficiency of nitrogen fixation.

Effect of interaction between a specific subtype of child abuse and the FKBP5 rs1360780 SNP on DNA methylation among patients with bipolar disorder

BACKGROUND

Child abuse is a risk factor for mood disorders, and linked to decreased DNA methylation (DNAm) of FKBP5 intron 7 through interactions with the single nucleotide polymorphism (SNP) rs1360780. However, no study has investigated which specific subtypes of child abuse are related to decreased DNAm of FKBP5 intron 7 in mood disorders. We therefore aimed to examine the relationship among various subtypes of child abuse, rs1360780, and the DNAm level of FKBP5 intron 7.

METHODS

A total of 190 subjects (87 patients with major depressive disorder [MDD], 61 patients with bipolar disorder [BD], and 42 healthy controls) participated. The Child Abuse and Trauma Scale (CATS) was used to evaluate child abuse. Whole blood was processed for genotyping, and pyrosequencing was conducted to assess the DNAm level of FKBP5 intron 7. A multiple regression analysis was used to analyze the DNAm level as a dependent variable, and the CATS subtypes and rs1360780 were used as independent variables.

RESULTS

Emotional abuse/neglect, one of the specific subtypes of child abuse, was related to lower DNAm of FKBP5 intron 7 interacting with rs1360780 in the BD patients. There were no significant results in the MDD patients or the controls.

LIMITATIONS

Since the study was limited to Japanese individuals, particularly those with MDD and BD, the findings are not generalizable. Furthermore, as child abuse was measured retrospectively, there may be recall bias.

CONCLUSIONS

This finding indicates that a specific subtype of child abuse may play an important role in the development of BD.

The geographic mosaic of arms race coevolution is closely matched to prey population structure

Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of "trait remixing" in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

HW_TEST, a program for comprehensive HARDY-WEINBERG equilibrium testing

This article deals with a Windows (© Microsoft Inc.) executable, user-friendly program that tests the hypothesis of Hardy-Weinberg (HW) proportions from autosomal multiallelic data using different methods that include parametric, nonparametric and exact bootstrap tests, the latter obtained through computer simulations. The program can be obtained free of charge directly from the internet repository https://github.com/Lemes-RenanB/HardyWeinbergTesting.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis)

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Na_v proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Na_v1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

Machine learning approach to single nucleotide polymorphism-based asthma prediction

Machine learning (ML) is poised as a transformational approach uniquely positioned to discover the hidden biological interactions for better prediction and diagnosis of complex diseases. In this work, we integrated ML-based models for feature selection and classification to quantify the risk of individual susceptibility to asthma using single nucleotide polymorphism (SNP). Random forest (RF) and recursive feature elimination (RFE) algorithm were implemented to identify the SNPs with high implication to asthma. K-nearest neighbor (kNN) and support vector machine (SVM) algorithms were trained to classify the identified SNPs whether associated with non-asthmatic or asthmatic samples. Feature selection step showed that RF outperformed RFE and the feature importance score derived from RF was consistently high for a subset of SNPs, indicating the robustness of RF in selecting relevant features associated with asthma. Model comparison showed that the integration of RF-SVM obtained the highest model performance with an accuracy, precision, and sensitivity of 62.5%, 65.3%, and 69%, respectively, when compared to the baseline, RF-kNN, and an external MeanDiff-kNN models. Furthermore, results show that the occurrence of asthma can be predicted with an Area under the Curve (AUC) of 0.62 and 0.64 for RF-SVM and RF-kNN models, respectively. This study demonstrates the integration of ML models to augment traditional methods in predicting genetic predisposition to multifactorial diseases such as asthma.

The Impact of CYP3A4*22 on Tacrolimus Pharmacokinetics and Outcome in Clinical Practice at a Single Kidney Transplant Center

Background: Although there is evidence that the CYP3A4*22 variant should be considered in tacrolimus dosing in renal transplantation, its impact beyond tacrolimus dose requirements remains controversial. Methods: In a cohort of 121 kidney transplant recipients, we analyzed the CYP3A4*1B, CYP3A4*22, and CYP3A5*3 alleles and the ABCB1 variants 1236C>T, 2677G>T/A, and 3435C>T for their impact on exposure and dose requirement. Relevant clinical outcome measures such as acute rejection within the first year after transplantation, delayed graft function, and renal function at discharge (estimated glomerular filtration rate) were evaluated. Results: Extensive metabolizer (n = 17, CYP3A4*1/*1 carriers with at least one CYP3A5*1 allele) showed significantly higher tacrolimus dose requirement (P = 0.004) compared with both intermediate metabolizer (IM, n = 93, CYP3A5*3/*3 plus CYP3A4*1/*1 or CYP3A4*22 carriers plus one CYP3A5*1 allele), and poor metabolizer (n = 11, CYP3A4*22 allele in combination with CYP3A5*3/*3) after onset of therapy. Significantly higher dose requirement was observed in CYP3A5 expressers (P = 0.046) compared with non-expressers again at onset of therapy. Using the log additive genetic model, the area under the curve for the total observation period up to 16 days was significantly associated with the CYP3A5*3 genotype (P = 3.34 × 10-4) as well as with the IM or extensive metabolizer phenotype (P = 1.54 × 10-4), even after adjustment for multiple testing. Heterozygous carriers for CYP3A4*22 showed significantly higher areas under the curve than the CYP3A4*1/*1 genotype in the second week post-transplantation (adjusted P = 0.016). Regarding clinical outcomes, acute rejection was significantly associated with human leukocyte antigen mismatch (≥3 alleles; OR = 12.14, 95% CI 1.76, 525.21, P = 0.019 after correction for multiple testing). Graft recipients from deceased donors showed higher incidende of delayed graft function (OR 7.15, 95% CI 2.23, 30.46, adjusted P = 0.0008) and a lower estimated glomerular filtration rate at discharge (P = 0.0001). Tested CYP3A4 or CYP3A5 variants did not show any effects on clinical outcome parameters. ABCB1 variants did neither impact on pharmacokinetics nor on clinical endpoints. Conclusion: At our transplantation center, both CYP3A5*3 and, to a lesser extent, CYP3A4*22 affect tacrolimus pharmacokinetics early after onset of therapy with consequences for steady-state treatment in routine clinical practice.

[A comparative analysis of allele frequencies of rs1801133 and rs1801131 of MTHFR in patients with stroke and healthy people from the Moscow region]

AIM

To study genetic characteristics of the population of the Moscow region and analyze the association of rs1801133 and rs1801131 of MTHFR with the risk of ischemic stroke (IS).

MATERIAL AND METHODS

A sample of 170 and 115 patients with atherothrombotic and cardioembolic subtypes of IS and 360 residents of the Moscow region without IS were examined. MTHFR alleles were determined by a multiplex real-time polymerase chain reaction.

RESULTS AND CONCLUSION

No association between the frequencies of MTHFR alleles and the risk of ischemic stroke was found. The comparison of allele frequencies with those in Caucasian populations published in the dbSNP (NCBI) and 1000 Genomes Project databases revealed significant differences for rs1801133 from the EUR 1000 Genomes Project. The allele frequency data for MTHFR could increase the accuracy and reliability of the individual risk calculation for multifactorial diseases in the Russian population.

Bayesian model selection for the study of Hardy-Weinberg proportions and homogeneity of gender allele frequencies

Standard statistical tests for Hardy-Weinberg equilibrium assume the equality of allele frequencies in the sexes, whereas tests for the equality of allele frequencies in the sexes assume Hardy-Weinberg equilibrium. This produces a circularity in the testing of genetic variants, which has recently been resolved with new frequentist likelihood and exact procedures. In this paper, we tackle the same problem by posing it as a Bayesian model comparison problem. We formulate an exhaustive set of ten alternative scenarios for biallelic genetic variants. Using Dirichlet and Beta priors for genotype and allele frequencies, we derive marginal likelihoods for all scenarios, and select the most likely scenario using the posterior probabilities that each of these scenarios is the one in place. Different from the usual frequentist testing approach, the Bayesian approach allows one to compare any number of models, and not just two at a time, and the models compared do not have to be nested. We illustrate our Bayesian approach with genetic data from the 1,000 genomes project and through a simulation study.

Genomic analysis of MHC-based mate choice in the monogamous California mouse

Variation at Major Histocompatibility Complex (MHC) genes is thought to be an important mechanism underlying mate choice in vertebrates, with individuals typically predicted to prefer MHC-dissimilar reproductive partners. However, analyses based on individual MHC loci have generated contradictory results regarding the role of these genes in mate-choice decisions. To provide a more comprehensive assessment of relationships between MHC variation and mating behavior, we used an exome capture strategy to characterize variability at 13 MHC loci, 312 innate immune system genes, and 1044 nonimmune genes in 25 obligate monogamous pairs of California mice (Peromyscus californicus) from 2 free-living populations of this species in Monterey County, California. Pairwise genotypic comparisons and analyses of SNP-based allelic differences failed to detect disassortative mating based on MHC variability; reproductive partners were not more dissimilar than randomly generated male-female pairs at MHC, innate or nonimmune loci. Within populations, individuals tended to be more closely related at MHC genes than at innate or nonimmune genes. Consistent with the functional role of immunogenes, the 2 study populations were highly differentiated at MHC and innate genes but not at nonimmune loci. Collectively, our results suggest that MHC genetic variation in California mice reflects local differences in pathogen exposure rather than disassortative mating based on variability at MHC Class I and II genes.

Genetic variants in microRNA biogenesis genes as novel indicators for secondary growth in Populus

MicroRNAs (miRNAs) function as key regulators of complex traits, but how genetic alterations in miRNA biogenesis genes (miRBGs) affect quantitative variation has not been elucidated. We conducted transcript analyses and association genetics to investigate how miRBGs, miRNA genes (MIRNAs) and their respective targets contribute to secondary growth in a natural population of 435 Populus tomentosa individuals. This analysis identified 29 843 common single-nucleotide polymorphisms (SNPs; frequency > 0.10) within 682 genes (80 miRBGs, 152 MIRNAs, and 457 miRNA targets). Single-SNP association analysis found SNPs in 234 candidate genes exhibited significant additive/dominant effects on phenotypes. Among these, specific candidates that associated with the same traits produced 791 miRBG-MIRNA-target combinations, suggesting possible genetic miRBG-MIRNA and MIRNA-target interactions, providing an important clue for the regulatory mechanisms of miRBGs. Multi-SNP association found 4672 epistatic pairs involving 578 genes that showed significant associations with traits and identified 106 miRBG-MIRNA-target combinations. Two multi-hierarchical networks were constructed based on correlations of miRBG-miRNA and miRNA-target expression to further probe the mechanisms of trait diversity underlying changes in miRBGs. Our study opens avenues for the investigation of miRNA function in perennial plants and underscored miRBGs as potentially modulating quantitative variation in traits.

Common marmoset (Callithrix jacchus) personality, subjective well-being, hair cortisol level and AVPR1a, OPRM1, and DAT genotypes

We studied personality, subjective well-being, and hair cortisol level, in common marmosets Callithrix jacchus, a small, cooperatively breeding New World monkey, by examining their associations with one another and genotypes. Subjects were 68 males and 9 females that lived in the RIKEN Center for Life Science Technologies. Personality and subjective well-being were assessed by keeper ratings on two questionnaires, hair samples were obtained to assay cortisol level and buccal swabs were used to assess AVPR1a, OPRM1 and DAT genotypes. Three personality domains-Dominance, Sociability, and Neuroticism-were identified. Consistent with findings in other species, Sociability and Neuroticism were related to higher and lower subjective well-being, respectively. Sociability was also associated with higher hair cortisol levels. The personality domains and hair cortisol levels were heritable and associated with genotypes: the short form of AVPR1a was associated with lower Neuroticism and the AA genotype of the A111T SNP of OPRM1 was related to lower Dominance, lower Neuroticism, and higher hair cortisol level. Some genetic associations were not in directions that one would expect given findings in other species. These findings provide insights into the proximate and ultimate bases of personality in common marmosets, other primates and humans.

Association of Polymorphisms in CYBA, SOD1, and CAT Genes with Type 1 Diabetes and Diabetic Peripheral Neuropathy in Children and Adolescents

AIMS

The aim of our study was to investigate possible associations between three SNPs: rs4673 in the CYBA gene; rs1041740 in the SOD1 gene; and rs1001179 in the CAT gene, and type 1 diabetes (T1D) or diabetic peripheral neuropathy (DPN) in T1D patients.

MATERIALS AND METHODS

Allelic variants of the selected SNPs were determined by allelic discrimination assays in 114 T1D patients enrolled in the study group and in 90 healthy individuals from a control group. Associations between each of the three SNPs were tested in subgroups of T1D patients divided according to the presence of DPN.

RESULTS

The TT genotype of rs4673 in the CYBA gene was associated with DPN in T1D patients (OR 4.997, 95% CI 1.403-19.083, p = 0.016). Weak significance was observed for a protective effect of the TT genotype of rs1041740 in the SOD1 gene relative to T1D development (OR 0.318, 95% CI 0.092-0.959, p = 0.056). There was no significant association between the CAT gene SNP rs1001179 and T1D or DPN.

CONCLUSION

We showed a strong association of the CYBA polymorphism rs4673 with DPN in Slovak children and adolescents with T1D. Further studies are necessary to assess the relationship between rs1041740 and T1D or DPN.

Two complement receptor one alleles have opposing associations with cerebral malaria and interact with α+thalassaemia

Malaria has been a major driving force in the evolution of the human genome. In sub-Saharan African populations, two neighbouring polymorphisms in the Complement Receptor One (CR1) gene, named Sl2 and McCb, occur at high frequencies, consistent with selection by malaria. Previous studies have been inconclusive. Using a large case-control study of severe malaria in Kenyan children and statistical models adjusted for confounders, we estimate the relationship between Sl2 and McCb and malaria phenotypes, and find they have opposing associations. The Sl2 polymorphism is associated with markedly reduced odds of cerebral malaria and death, while the McCb polymorphism is associated with increased odds of cerebral malaria. We also identify an apparent interaction between Sl2 and α+thalassaemia, with the protective association of Sl2 greatest in children with normal α-globin. The complex relationship between these three mutations may explain previous conflicting findings, highlighting the importance of considering genetic interactions in disease-association studies.

MetaGenyo: a web tool for meta-analysis of genetic association studies

Background

Genetic association studies (GAS) aims to evaluate the association between genetic variants and phenotypes. In the last few years, the number of this type of study has increased exponentially, but the results are not always reproducible due to experimental designs, low sample sizes and other methodological errors. In this field, meta-analysis techniques are becoming very popular tools to combine results across studies to increase statistical power and to resolve discrepancies in genetic association studies. A meta-analysis summarizes research findings, increases statistical power and enables the identification of genuine associations between genotypes and phenotypes. Meta-analysis techniques are increasingly used in GAS, but it is also increasing the amount of published meta-analysis containing different errors. Although there are several software packages that implement meta-analysis, none of them are specifically designed for genetic association studies and in most cases their use requires advanced programming or scripting expertise.

Results

We have developed MetaGenyo, a web tool for meta-analysis in GAS. MetaGenyo implements a complete and comprehensive workflow that can be executed in an easy-to-use environment without programming knowledge. MetaGenyo has been developed to guide users through the main steps of a GAS meta-analysis, covering Hardy-Weinberg test, statistical association for different genetic models, analysis of heterogeneity, testing for publication bias, subgroup analysis and robustness testing of the results.

Conclusions

MetaGenyo is a useful tool to conduct comprehensive genetic association meta-analysis. The application is freely available at http://bioinfo.genyo.es/metagenyo/.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1990-4) contains supplementary material, which is available to authorized users.

On the testing of Hardy-Weinberg proportions and equality of allele frequencies in males and females at biallelic genetic markers

Standard statistical tests for equality of allele frequencies in males and females and tests for Hardy‐Weinberg equilibrium are tightly linked by their assumptions. Tests for equality of allele frequencies assume Hardy‐Weinberg equilibrium, whereas the usual chi‐square or exact test for Hardy‐Weinberg equilibrium assume equality of allele frequencies in the sexes. In this paper, we propose ways to break this interdependence in assumptions of the two tests by proposing an omnibus exact test that can test both hypotheses jointly, as well as a likelihood ratio approach that permits these phenomena to be tested both jointly and separately. The tests are illustrated with data from the 1000 Genomes project.

Genome-Wide SNP Discovery and Analysis of Genetic Diversity in Farmed Sika Deer (Cervus nippon) in Northeast China Using Double-Digest Restriction Site-Associated DNA Sequencing

Sika deer are an economically valuable species owing to their use in traditional Chinese medicine, particularly their velvet antlers. Sika deer in northeast China are mostly farmed in enclosure. Therefore, genetic management of farmed sika deer would benefit from detailed knowledge of their genetic diversity. In this study, we generated over 1.45 billion high-quality paired-end reads (288 Gbp) across 42 unrelated individuals using double-digest restriction site-associated DNA sequencing (ddRAD-seq). A total of 96,188 (29.63%) putative biallelic SNP loci were identified with an average sequencing depth of 23×. Based on the analysis, we found that the majority of the loci had a deficit of heterozygotes (F_IS >0) and low values of H_obs, which could be due to inbreeding and Wahlund effects. We also developed a collection of high-quality SNP probes that will likely be useful in a variety of applications in genotyping for cervid species in the future.

Genetically-Driven Enhancement of Dopaminergic Transmission Affects Moral Acceptability in Females but Not in Males: A Pilot Study

Moral behavior has been a key topic of debate for philosophy and psychology for a long time. In recent years, thanks to the development of novel methodologies in cognitive sciences, the question of how we make moral choices has expanded to the study of neurobiological correlates that subtend the mental processes involved in moral behavior. For instance, in vivo brain imaging studies have shown that distinct patterns of brain neural activity, associated with emotional response and cognitive processes, are involved in moral judgment. Moreover, while it is well-known that responses to the same moral dilemmas differ across individuals, to what extent this variability may be rooted in genetics still remains to be understood. As dopamine is a key modulator of neural processes underlying executive functions, we questioned whether genetic polymorphisms associated with decision-making and dopaminergic neurotransmission modulation would contribute to the observed variability in moral judgment. To this aim, we genotyped five genetic variants of the dopaminergic pathway [rs1800955 in the dopamine receptor D4 (DRD4) gene, DRD4 48 bp variable number of tandem repeat (VNTR), solute carrier family 6 member 3 (SLC6A3) 40 bp VNTR, rs4680 in the catechol-O-methyl transferase (COMT) gene, and rs1800497 in the ankyrin repeat and kinase domain containing 1 (ANKK1) gene] in 200 subjects, who were requested to answer 56 moral dilemmas. As these variants are all located in genes belonging to the dopaminergic pathway, they were combined in multilocus genetic profiles for the association analysis. While no individual variant showed any significant effects on moral dilemma responses, the multilocus genetic profile analysis revealed a significant gender-specific influence on human moral acceptability. Specifically, those genotype combinations that improve dopaminergic signaling selectively increased moral acceptability in females, by making their responses to moral dilemmas more similar to those provided by males. As females usually give more emotionally-based answers and engage the "emotional brain" more than males, our results, though preliminary and therefore in need of replication in independent samples, suggest that this increase in dopamine availability enhances the cognitive and reduces the emotional components of moral decision-making in females, thus favoring a more rationally-driven decision process.

Testing Departure from Hardy-Weinberg Proportions

The Hardy-Weinberg principle, one of the most important principles in population genetics, was originally developed for the study of allele frequency changes in a population over generations. It is now, however, widely used in studies of human diseases to detect inbreeding, population stratification, and genotyping errors. For assessment of deviation from Hardy-Weinberg proportions in data, the most popular approaches include the asymptotic Pearson's chi-squared goodness-of-fit test and the exact test. Pearson's chi-squared goodness-of-fit test is simple and straightforward, but is very sensitive to a small sample size or rare allele frequency. The exact test of Hardy-Weinberg proportions is preferable in these situations. The exact test can be performed through complete enumeration of heterozygote genotypes or on the basis of the Markov chain Monte Carlo procedure. In this chapter, we describe the Hardy-Weinberg principle and the commonly used Hardy-Weinberg proportion tests and their applications, and we demonstrate how the chi-squared test and exact test of Hardy-Weinberg proportions can be performed step-by-step using the popular software programs SAS, R, and PLINK, which have been widely used in genetic association studies, along with numerical examples. We also discuss approaches for testing Hardy-Weinberg proportions in case-control study designs that are better than traditional approaches for testing Hardy-Weinberg proportions in controls only. Finally, we note that deviation from the Hardy-Weinberg proportions in affected individuals can provide evidence for an association between genetic variants and diseases.

Chimpanzee Personality and the Arginine Vasopressin Receptor 1A Genotype

Polymorphisms of the arginine vasopressin receptor 1a (AVPR1a) gene have been linked to various measures related to human social behavior, including sibling conflict and agreeableness. In chimpanzees, AVPR1a polymorphisms have been associated with traits important for social interactions, including sociability, joint attention, dominance, conscientiousness, and hierarchical personality dimensions named low alpha/stability, disinhibition, and negative emotionality/low dominance. We examined associations between AVPR1a and six personality domains and hierarchical personality dimensions in 129 chimpanzees (Pan troglodytes) living in Japan or in a sanctuary in Guinea. We fit three linear and three animal models. The first model included genotype, the second included sex and genotype, and the third included genotype, sex, and sex × genotype. All personality phenotypes were heritable. Chimpanzees possessing the long form of the allele were higher in conscientiousness, but only in models that did not include the other predictors; however, additional analyses suggested that this may have been a consequence of study design. In animal models that included sex and sex × genotype, chimpanzees homozygous for the short form of the allele were higher in extraversion. Taken with the findings of previous studies of chimpanzees and humans, the findings related to conscientiousness suggest that AVPR1a may be related to lower levels of impulsive aggression. The direction of the association between AVPR1a genotype and extraversion ran counter to what one would expect if AVPR1a was related to social behaviors. These results help us further understand the genetic basis of personality in chimpanzees.

Rare sex or out of reach equilibrium? The dynamics of F IS in partially clonal organisms

BACKGROUND

Partially clonal organisms are very common in nature, yet the influence of partial asexuality on the temporal dynamics of genetic diversity remains poorly understood. Mathematical models accounting for clonality predict deviations only for extremely rare sex and only towards mean inbreeding coefficient [Formula: see text]. Yet in partially clonal species, both F IS  < 0 and F IS  > 0 are frequently observed also in populations where there is evidence for a significant amount of sexual reproduction. Here, we studied the joint effects of partial clonality, mutation and genetic drift with a state-and-time discrete Markov chain model to describe the dynamics of F IS over time under increasing rates of clonality.

RESULTS

Results of the mathematical model and simulations show that partial clonality slows down the asymptotic convergence to F IS  = 0. Thus, although clonality alone does not lead to departures from Hardy-Weinberg expectations once reached the final equilibrium state, both negative and positive F IS values can arise transiently even at intermediate rates of clonality. More importantly, such "transient" departures from Hardy Weinberg proportions may last long as clonality tunes up the temporal variation of F IS and reduces its rate of change over time, leading to a hyperbolic increase of the maximal time needed to reach the final mean [Formula: see text] value expected at equilibrium.

CONCLUSION

Our results argue for a dynamical interpretation of F IS in clonal populations. Negative values cannot be interpreted as unequivocal evidence for extremely scarce sex but also as intermediate rates of clonality in finite populations. Complementary observations (e.g. frequency distribution of multiloci genotypes, population history) or time series data may help to discriminate between different possible conclusions on the extent of clonality when mean [Formula: see text] values deviating from zero and/or a large variation of F IS over loci are observed.

Testing for Hardy-Weinberg equilibrium at biallelic genetic markers on the X chromosome

Testing genetic markers for Hardy–Weinberg equilibrium (HWE) is an important tool for detecting genotyping errors in large-scale genotyping studies. For markers at the X chromosome, typically the χ² or exact test is applied to the females only, and the hemizygous males are considered to be uninformative. In this paper we show that the males are relevant, because a difference in allele frequency between males and females may indicate HWE not to hold. The testing of markers on the X chromosome has received little attention, and in this paper we lay down the foundation for testing biallelic X-chromosomal markers for HWE. We develop four frequentist statistical test procedures for X-linked markers that take both males and females into account: the χ² test, likelihood ratio test, exact test and permutation test. Exact tests that include males are shown to have a better Type I error rate. Empirical data from the GENEVA project on venous thromboembolism is used to illustrate the proposed tests. Results obtained with the new tests differ substantially from tests that are based on female genotype counts only. The new tests detect differences in allele frequencies and seem able to uncover additional genotyping error that would have gone unnoticed in HWE tests based on females only.

Biomass traits and candidate genes for bioenergy revealed through association genetics in coppiced European Populus nigra (L.)

BACKGROUND

Second generation (2G) bioenergy from lignocellulosic feedstocks has the potential to develop as a sustainable source of renewable energy; however, significant hurdles still remain for large-scale commercialisation. Populus is considered as a promising 2G feedstock and understanding the genetic basis of biomass yield and feedstock quality are a research priority in this model tree species.

RESULTS

We report the first coppiced biomass study for 714 members of a wide population of European black poplar (Populus nigra L.), a native European tree, selected from 20 river populations ranging in latitude and longitude between 40.5 and 52.1°N and 1.0 and 16.4°E, respectively. When grown at a single site in southern UK, significant Site of Origin (SO) effects were seen for 14 of the 15 directly measured or derived traits including biomass yield, leaf area and stomatal index. There was significant correlation (p < 0.001) between biomass yield traits over 3 years of harvest which identified leaf size and cell production as strong predictors of biomass yield. A 12 K Illumina genotyping array (constructed from 10,331 SNPs in 14 QTL regions and 4648 genes) highlighted significant population genetic structure with pairwise FST showing strong differentiation (p < 0.001) between the Spanish and Italian subpopulations. Robust associations reaching genome-wide significance are reported for main stem height and cell number per leaf; two traits tightly linked to biomass yield. These genotyping and phenotypic data were also used to show the presence of significant isolation by distance (IBD) and isolation by adaption (IBA) within this population.

CONCLUSIONS

The three associations identified reaching genome-wide significance at p < 0.05 include a transcription factor; a putative stress response gene and a gene of unknown function. None of them have been previously linked to bioenergy yield; were shown to be differentially expressed in a panel of three selected genotypes from the collection and represent exciting, novel candidates for further study in a bioenergy tree native to Europe and Euro-Asia. A further 26 markers (22 genes) were found to reach putative significance and are also of interest for biomass yield, leaf area, epidermal cell expansion and stomatal patterning. This research on European P. nigra provides an important foundation for the development of commercial native trees for bioenergy and for advanced, molecular breeding in these species.

No Association between Personality and Candidate Gene Polymorphisms in a Wild Bird Population

Consistency of between-individual differences in behaviour or personality is a phenomenon in populations that can have ecological consequences and evolutionary potential. One way that behaviour can evolve is to have a genetic basis. Identifying the molecular genetic basis of personality could therefore provide insight into how and why such variation is maintained, particularly in natural populations. Previously identified candidate genes for personality in birds include the dopamine receptor D4 (DRD4), and serotonin transporter (SERT). Studies of wild bird populations have shown that exploratory and bold behaviours are associated with polymorphisms in both DRD4 and SERT. Here we tested for polymorphisms in DRD4 and SERT in the Seychelles warbler (Acrocephalus sechellensis) population on Cousin Island, Seychelles, and then investigated correlations between personality and polymorphisms in these genes. We found no genetic variation in DRD4, but identified four polymorphisms in SERT that clustered into five haplotypes. There was no correlation between bold or exploratory behaviours and SERT polymorphisms/haplotypes. The null result was not due to lack of power, and indicates that there was no association between these behaviours and variation in the candidate genes tested in this population. These null findings provide important data to facilitate representative future meta-analyses on candidate personality genes.

FTO Is Associated with Aortic Valve Stenosis in a Gender Specific Manner of Heterozygote Advantage: A Population-Based Case-Control Study

Background

Single nucleotide polymorphisms (SNPs) within the Fat mass and obesity associated (FTO) gene have been linked with increased body weight. However, the data on an association of FTO with cardiovascular diseases remains conflicting. Therefore, we ascertained whether FTO is associated with aortic valve stenosis (AVS), one of the most frequent cardiovascular diseases in the Western world.

Methods and Findings

In this population-based case-control study the FTO SNP rs9939609 was analyzed in 300 German patients with AVS and 429 German controls of the KORA survey S4, representing a random population. Blood samples were collected prior to aortic valve replacement in AVS cases and FTO rs9939609 was genotyped via ARMS-PCR. Genotype frequencies differed significantly between AVS cases and KORA controls (p = 0.004). Separate gender-analyses uncovered an association of FTO with AVS exclusively in males; homozygote carriers for the risk-allele (A) had a higher risk to develop AVS (p = 0.017, odds ratio (OR) 1.727; 95% confidence interval (CI) 1.087–2.747, recessive model), whereas heterozygote carriers for the risk-allele showed a lower risk (p = 0.002, OR 0.565, 95% CI 0.384–0.828, overdominant model). After adjustment for multiple co-variables, the odds ratios of heterozygotes remained significant for an association with AVS (p = 0.008, OR 0.565, 95% CI 0.369–0.861).

Conclusions

This study revealed an association of FTO rs9939609 with AVS. Furthermore, this association was restricted to men, with heterozygotes having a significantly lower chance to develop AVS. Lastly, the association between FTO and AVS was independent of BMI and other variables such as diabetes mellitus.

No association between genetic or epigenetic variation in insulin growth factors and antipsychotic-induced metabolic disturbances in a cross-sectional sample

AIM

Second-generation antipsychotics (SGA) are known to induce metabolic disturbances. Genetic pathways, such as the IGF pathway could be associated with increased metabolic syndrome (MetS). Additionally, IGF2 methylation varies as a function of environmental influences and is associated with schizophrenia and MetS. The current study aims to evaluate whether genetic and epigenetic variation in genes of the IGF pathway are associated with metabolic disturbances in patients under treatment with SGAs.

METHODS

Cross-sectional metabolic data from 438 patients with schizophrenia spectrum disorder was analyzed. Using the Sequenom MassARRAY iPLEX(TM) platform, 27 SNPs of the IGF1 and IGF2 genes and the IGF receptors IGF1R and IGF2R were genotyped. Methylation status of seven IGF2 CpG dinucleotides was evaluated using a Sequenom MALDI-TOF spectrometer.

RESULTS & CONCLUSION

There was a significant association between IGF2 methylation and genotype, but no significant association between genetic or epigenetic variability and metabolic parameters in the present study.

Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P<1.7×10(-6)). Markers were found across 19 chromosomes, explained 1-13% of trait variation, and implicated 275 unique genes in trait associations. Phenology had the largest number of associated genes (240 genes), followed by biomass (53 genes) and ecophysiology traits (25 genes). The GWAS results propose numerous loci for further investigation. Many traits had significant associations with multiple genes, underscoring their genetic complexity. Genes were also identified with multiple trait associations within and/or across trait categories. In some cases, traits were genetically correlated while in others they were not.

Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa

• Populus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities. • We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29,354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H(2)) and overall population genetic structure from principal component analysis. • Populus trichocarpa had high phenotypic variation and moderate/high H(2) for many traits. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. • Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

Statistical inference for Hardy-Weinberg proportions in the presence of missing genotype information

In genetic association studies, tests for Hardy-Weinberg proportions are often employed as a quality control checking procedure. Missing genotypes are typically discarded prior to testing. In this paper we show that inference for Hardy-Weinberg proportions can be biased when missing values are discarded. We propose to use multiple imputation of missing values in order to improve inference for Hardy-Weinberg proportions. For imputation we employ a multinomial logit model that uses information from allele intensities and/or neighbouring markers. Analysis of an empirical data set of single nucleotide polymorphisms possibly related to colon cancer reveals that missing genotypes are not missing completely at random. Deviation from Hardy-Weinberg proportions is mostly due to a lack of heterozygotes. Inbreeding coefficients estimated by multiple imputation of the missings are typically lowered with respect to inbreeding coefficients estimated by discarding the missings. Accounting for missings by multiple imputation qualitatively changed the results of 10 to 17% of the statistical tests performed. Estimates of inbreeding coefficients obtained by multiple imputation showed high correlation with estimates obtained by single imputation using an external reference panel. Our conclusion is that imputation of missing data leads to improved statistical inference for Hardy-Weinberg proportions.

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

Association analysis identifies Melampsora ×columbiana poplar leaf rust resistance SNPs

Populus species are currently being domesticated through intensive time- and resource-dependent programs for utilization in phytoremediation, wood and paper products, and conversion to biofuels. Poplar leaf rust disease can greatly reduce wood volume. Genetic resistance is effective in reducing economic losses but major resistance loci have been race-specific and can be readily defeated by the pathogen. Developing durable disease resistance requires the identification of non-race-specific loci. In the presented study, area under the disease progress curve was calculated from natural infection of Melampsora ×columbiana in three consecutive years. Association analysis was performed using 412 P. trichocarpa clones genotyped with 29,355 SNPs covering 3,543 genes. We found 40 SNPs within 26 unique genes significantly associated (permutated P<0.05) with poplar rust severity. Moreover, two SNPs were repeated in all three years suggesting non-race-specificity and three additional SNPs were differentially expressed in other poplar rust interactions. These five SNPs were found in genes that have orthologs in Arabidopsis with functionality in pathogen induced transcriptome reprogramming, Ca²⁺/calmodulin and salicylic acid signaling, and tolerance to reactive oxygen species. The additive effect of non-R gene functional variants may constitute high levels of durable poplar leaf rust resistance. Therefore, these findings are of significance for speeding the genetic improvement of this long-lived, economically important organism.

Nicotinic receptor gene variants interact with attention deficient hyperactive disorder symptoms to predict smoking trajectories from early adolescence to adulthood

OBJECTIVE

To examine the association of single nucleotide polymorphisms (SNPs) of the CHRNB3 (rs13280604) and CHRNA6 (rs892413) nicotinic acetylcholine receptor (nAChR) genes and symptoms of attention deficit hyperactivity disorder (ADHD) in predicting smoking patterns from early adolescence to adulthood.

METHOD

A longitudinal cohort of 1137 unrelated youths from the National Longitudinal Study of Adolescent Health provided responses to four surveys from Waves I to IV, and a genetic sample in Wave III. Growth mixture modeling was used to identify smoking patterns and to assess the effects of the two SNPs and ADHD symptoms on cigarette use over time.

RESULTS

There were significant main effects of ADHD symptoms and CHRNA6 variants in predicting the number of cigarettes smoked and the pattern of use over time, respectively. There were no main effects of the CHRNB3 variants. However, a significant CHRNB3 variant×ADHD symptom interaction was observed, such that individuals with elevated ADHD symptoms and a particular CHRNB3 variant were at increased risk of cigarette use over time.

CONCLUSIONS

These findings demonstrate that a SNP in a nicotinic receptor gene may interact with ADHD symptoms to link with increased cigarette use across adolescence and young adulthood. Unique associations between specific variants and patterns of ADHD symptoms were identified which may be useful for targeting prevention efforts to individuals at greatest risk for cigarette smoking.

Plasma Dimethylglycine and Risk of Incident Acute Myocardial Infarction in Patients With Stable Angina Pectoris

Objective—

Dimethylglycine is linked to lipid metabolism, and increased plasma levels may be associated with adverse prognosis in patients with coronary artery disease. We evaluated the relationship between plasma dimethylglycine and risk of incident acute myocardial infarction in a large prospective cohort of patients with stable angina pectoris, of whom approximately two thirds were participants in a B-vitamin intervention trial. Model discrimination and reclassification when adding plasma dimethylglycine to established risk factors were obtained. We also explored temporal changes and the test–retest reliability of plasma dimethylglycine.

Approach and Results—

Four thousand one hundred fifty patients (72% men; median age 62 years) were included. Plasma dimethylglycine was associated with several traditional coronary artery disease risk factors. During a median follow-up of 4.6 years, 343 (8.3%) patients experienced an acute myocardial infarction. The hazard ratio (95% confidence interval) for acute myocardial infarction was 1.95 (1.42–2.68; P <0.001) when comparing plasma dimethylglycine quartile 4 to 1 in a Cox regression model adjusted for age, sex, and fasting status. Adjusting for traditional coronary artery disease risk factors only slightly modified the estimates, which were particularly strong among nonsmokers and among patients with serum triglyceride or apolipoprotein B100 levels ≤median ( P for interaction=0.004, 0.004, and 0.03, respectively). Plasma dimethylglycine improved discrimination and reclassification and had high test–retest reliability.

Conclusions—

Plasma dimethylglycine is independently related to incident acute myocardial infarction and enhances risk prediction in patients with stable angina pectoris. Our results motivate further studies on the relationship between 1-carbon metabolism and atherothrombosis. A potential interplay with lipid and energy metabolism merits particular attention.

Reducing bias of allele frequency estimates by modeling SNP genotype data with informative missingness

The presence of missing single-nucleotide polymorphism (SNP) genotypes is common in genetic studies. For studies with low-density SNPs, the most commonly used approach to dealing with genotype missingness is to simply remove the observations with missing genotypes from the analyses. This naïve method is straightforward but is valid only when the missingness is random. However, a given assay often has a different capability in genotyping heterozygotes and homozygotes, causing the phenomenon of "differential dropout" in the sense that the missing rates of heterozygotes and homozygotes are different. In practice, differential dropout among genotypes exists in even carefully designed studies, such as the data from the HapMap project and the Wellcome Trust Case Control Consortium. Under the assumption of Hardy-Weinberg equilibrium and no genotyping error, we here propose a statistical method to model the differential dropout among different genotypes. Compared with the naïve method, our method provides more accurate allele frequency estimates when the differential dropout is present. To demonstrate its practical use, we further apply our method to the HapMap data and a scleroderma data set.

Testing departure from Hardy-Weinberg proportions

The Hardy-Weinberg principle, one of the most important principles in population genetics, was originally developed for the study of allele frequency changes in a population over generations. It is now, however, widely used in studies of human diseases to detect inbreeding, populations stratification, and genotyping errors. For assessment of deviation from the Hardy-Weinberg proportions in data, the most popular approaches include the asymptotic Pearson's chi-square goodness-of-fit test and the exact test. The Pearson's chi-square goodness-of-fit test is simple and straightforward, but it is very sensitive to small sample size or rare allele frequency. The exact test of Hardy-Weinberg proportions is preferable in these situations. The exact test can be performed through complete enumeration of heterozygote genotypes or on the basis of the Markov chain Monte Carlo procedure. In this chapter, we describe the Hardy-Weinberg principle and the commonly used Hardy-Weinberg proportions tests and their applications, and we demonstrate how the chi-square test and exact test of Hardy-Weinberg proportions can be performed step-by-step using the popular software programs SAS, R, and PLINK, which have been widely used in genetic association studies, along with numerical examples. We also discuss recent approaches for testing Hardy-Weinberg proportions in case-control study designs that are better than traditional approaches for testing Hardy-Weinberg proportions in controls only. Finally, we note that deviation from the Hardy-Weinberg proportions in affected individuals can provide evidence for an association between genetic variants and diseases.

Possible genetic sign of heat stress adaptation in Hungarian Grey Bos taurus breed

The synthesis of Heat Shock Protein 70.2 mRNA is also regulated by the Upper Promoter elements of the gene. This promoter region is polymorphic in cattle. These polymorphisms have a major effect on the activity of the mRNA transcription. In a comparison of quantity of transcribed mRNA from the wild type and AP2 mutant allele the wild type can produce 2-3-fold more transcripts.The Hungarian Grey Cattle (HG) and Norwegian Red (NFR) as control breed were genotyped with PCR-RFLP method. Our results showed that the frequencies of alleles in breeds (p(wt)HG = 0.859419, p(wt)NFR = 0.5) are different. The effective response to heat stress in the Norwegian Red seems to be less important than in the Hungarian Grey breed. The extensive keeping in hot and arid region during centuries could have been proved as selection pressure for the heat tolerance.Our results combined with the global climate forecasts emphasize the role of autochthonous, well adopted, heat tolerant breeds in the near future.