Candidate genes versus genome-wide associations: which are better for detecting genetic susceptibility to infectious disease?

The potential role of RNA sequencing in diagnosing unexplained insensitivity to conventional chemotherapy in pediatric patients with B-cell acute lymphoblastic leukemia

Pediatric B-cell acute lymphoblastic leukemia (B-ALL) is a highly heterogeneous disease. According to large-scale RNA sequencing (RNA-seq) data, B-ALL patients can be divided into more than 10 subgroups. However, many genomic defects associated with resistance mechanisms have not yet been identified. As an individual clinical tool for molecular diagnostic risk classification, RNA-seq and gene expression pattern-based therapy could be potential upcoming strategies. In this study, we retrospectively analyzed the RNA-seq gene expression profiles of 45 children whose molecular diagnostic classifications were inconsistent with the response to chemotherapy. The relationship between the transcriptome and chemotherapy response was analyzed. Fusion gene identification was conducted for the included patients who did not have known high-risk associated fusion genes or gene mutations. The most frequently detected fusion gene pair in the high-risk group was the DHRSX duplication, which is a novel finding. Fusions involving ABL1, LMNB2, NFATC1, PAX5, and TTYH3 at onset were more frequently detected in the high-risk group, while fusions involving LFNG, TTYH3, and NFATC1 were frequently detected in the relapse group. According to the pathways involved, the underlying drug resistance mechanism is related to DNA methylation, autophagy, and protein metabolism. Overall, the implementation of an RNA-seq diagnostic system will identify activated markers associated with chemotherapy response, and guide future treatment adjustments.

Sex differences in the genetic and molecular mechanisms of coronary artery disease

Sex differences in coronary artery disease (CAD) presentation, risk factors and prognosis have been widely studied. Similarly, studies on atherosclerosis have shown prominent sex differences in plaque biology. Our understanding of the underlying genetic and molecular mechanisms that drive these differences remains fragmented and largely understudied. Through reviewing genetic and epigenetic studies, we identified more than 40 sex-differential candidate genes (13 within known CAD loci) that may explain, at least in part, sex differences in vascular remodeling, lipid metabolism and endothelial dysfunction. Studies with transcriptomic and single-cell RNA sequencing data from atherosclerotic plaques highlight potential sex differences in smooth muscle cell and endothelial cell biology. Especially, phenotypic switching of smooth muscle cells seems to play a crucial role in female atherosclerosis. This matches the known sex differences in atherosclerotic phenotypes, with men being more prone to lipid-rich plaques, while women are more likely to develop fibrous plaques with endothelial dysfunction. To unravel the complex mechanisms that drive sex differences in CAD, increased statistical power and adjustments to study designs and analysis strategies are required. This entails increasing inclusion rates of women, performing well-defined sex-stratified analyses and the integration of multi-omics data.

Multi-phasic gene profiling using candidate gene approach predict the capacity of specific antibody production and maintenance following COVID-19 vaccination in Japanese population

Background

Vaccination against severe acute respiratory syndrome coronavirus type 2 is highly effective in preventing infection and reducing the severity of coronavirus disease (COVID-19). However, acquired humoral immunity wanes within six months. Focusing on the different tempo of acquisition and attenuation of specific antibody titers in individuals, we investigated the impact of genetic polymorphisms on antibody production after COVID-19 vaccination.

Methods

In total 236 healthcare workers from a Japanese municipal hospital, who received two doses of the vaccine were recruited. We employed a candidate gene approach to identify the target genetic polymorphisms affecting antibody production after vaccination. DNA samples from the study populations were genotyped for 33 polymorphisms in 15 distinct candidate genes encoding proteins involved in antigen-presenting cell activation, T cell activation, T-B interaction, and B cell survival. We measured total anti-SARS-Cov2 spike IgG antibody titers and analyzed the association with genetic polymorphisms at several time points after vaccination using an unbiased statistical method, and stepwise logistic regression following multivariate regression.

Results

Significant associations were observed between seven SNPs in NLRP3, OAS1, IL12B, CTLA4, and IL4, and antibody titers at 3 weeks after the first vaccination as an initial response. Six SNPs in NLRP3, TNF, OAS1, IL12B, and CTLA4 were associated with high responders with serum antibody titer > 4000 BAU/ml as boosting effect at 3 weeks after the second vaccination. Analysis of long-term maintenance showed the significance of the three SNPs in IL12B, IL7R, and MIF for the maintenance of antibody titers and that in BAFF for attenuation of neutralizing antibodies. Finally, we proposed a predictive model composed of gene profiles to identify the individuals with rapid antibody attenuation by receiver operating characteristic (ROC) analysis (area under the curve (AUC)= 0.76, sensitivity = 82.5%, specificity=67.8%).

Conclusions

The candidate gene approach successfully showed shifting responsible gene profiles and initial and boosting effect mainly related to the priming phase into antibody maintenance including B cell survival, which traces the phase of immune reactions. These gene profiles provide valuable information for further investigation of humoral immunity against COVID-19 and for building a strategy for personalized vaccine schedules.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using F_ST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

Development of a robust and quantitative high-throughput screening method for assessing phenotypic variation in large Neisseria meningitidis isolate collections

Genome-wide association studies are a powerful approach for identifying determinants of disease. For infectious diseases, high throughput assays are required for measuring the variance in multiple virulence-related phenotypes of large bacterial isolate collections and for association of this phenotypic variance with genotype. The primary limiting factors are cost, effectiveness and a standardized inoculum. A method was developed to create an inoculum array of multiple isolates that could be used for a series of high-throughput multi-isolate phenotypic investigations in a laboratory setting. A key starting point was the standardisation of the inoculum by production of identical batches of each isolate from cells grown to mid-log phase. Cultures with pre-determined optical densities were aliquoted in set patterns into multiple multi-well plates containing 50% glycerol and stored at -80 °C. Prior to a specific assay, an inoculum plate was defrosted and subjected to a brief period of incubation. Control strains can be placed on each plate in order to control for intra-assay variability. A high throughput screen is described in detail for quantification of biofilm formation. This example utilised the crystal violet staining method and multi-assay stock plates containing 16 meningococcal isolates.•Multi-assay stock plate of exponentially growing isolates is cost-effective and simple to implement in a laboratory setting.•This method would predict realistic standard deviations for multiple isolates in phenotypic assays and generate data for performance of power calculations for genotyping.•This method has the potential to identify both known and unknown genetic determinants of phenotypic variability for each tested isolate when paired with genetic analysis of whole genome sequencing data.

A Systematic Review of Polygenic Models for Predicting Drug Outcomes

Polygenic models have emerged as promising prediction tools for the prediction of complex traits. Currently, the majority of polygenic models are developed in the context of predicting disease risk, but polygenic models may also prove useful in predicting drug outcomes. This study sought to understand how polygenic models incorporating pharmacogenetic variants are being used in the prediction of drug outcomes. A systematic review was conducted with the aim of gaining insights into the methods used to construct polygenic models, as well as their performance in drug outcome prediction. The search uncovered 89 papers that incorporated pharmacogenetic variants in the development of polygenic models. It was found that the most common polygenic models were constructed for drug dosing predictions in anticoagulant therapies (n = 27). While nearly all studies found a significant association with their polygenic model and the investigated drug outcome (93.3%), less than half (47.2%) compared the performance of the polygenic model against clinical predictors, and even fewer (40.4%) sought to validate model predictions in an independent cohort. Additionally, the heterogeneity of reported performance measures makes the comparison of models across studies challenging. These findings highlight key considerations for future work in developing polygenic models in pharmacogenomic research.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Genetic Predisposition to Persistent Human Papillomavirus-Infection and Virus-Induced Cancers

Human papillomaviruses (HPVs) are the most common sexually transmitted pathogens worldwide and among the more than 200 identified HPV types, approximately 15 high risk (HR-HPV) types are oncogenic, being strongly associated with the development of cervical cancer, anogenital cancers and an increasing fraction of head and neck squamous cell carcinomas (HNSCC). HPV-associated cervix cancer accounts for 83% of HPV-attributable cancers, and more than two-thirds of those cases occur in developing countries. Despite the high frequency of HPV infections, in most cases, the virus is cleared by the host immune response and only a small proportion of infected individuals develop persistent infections that can result in malignant transformation, indicating that other elements, including biological, genetic and environmental factors may influence the individual susceptibility to HPV-associated cancers. Previous studies have quantified that heritability, in the form of genetic variants, common in the general population, is implicated in nearly 30% of cervical cancers and a large number of studies conducted across various populations have identified genetic variants that appear to be associated with genes that predispose or protect the host to HPV infections thereby affecting individual susceptibility to HPV-associated cancers. In this article, we provide an overview of gene association studies on HPV-associated cancers with emphasis on genome-wide association study (GWAS) that have identified novel genetic factors linked to HPV infection or HPV-associated cancers.

Genome Wide Prediction, Mapping and Development of Genomic Resources of Mastitis Associated Genes in Water Buffalo

Water buffalo (Bubalus bubalis) are an important animal resource that contributes milk, meat, leather, dairy products, and power for plowing and transport. However, mastitis, a bacterial disease affecting milk production and reproduction efficiency, is most prevalent in populations having intensive selection for higher milk yield, especially where the inbreeding level is also high. Climate change and poor hygiene management practices further complicate the issue. The management of this disease faces major challenges, like antibiotic resistance, maximum residue level, horizontal gene transfer, and limited success in resistance breeding. Bovine mastitis genome wide association studies have had limited success due to breed differences, sample sizes, and minor allele frequency, lowering the power to detect the diseases associated with SNPs. In this work, we focused on the application of targeted gene panels (TGPs) in screening for candidate gene association analysis, and how this approach overcomes the limitation of genome wide association studies. This work will facilitate the targeted sequencing of buffalo genomic regions with high depth coverage required to mine the extremely rare variants potentially associated with buffalo mastitis. Although the whole genome assembly of water buffalo is available, neither mastitis genes are predicted nor TGP in the form of web-genomic resources are available for future variant mining and association studies. Out of the 129 mastitis associated genes of cattle, 101 were completely mapped on the buffalo genome to make TGP. This further helped in identifying rare variants in water buffalo. Eighty-five genes were validated in the buffalo gene expression atlas, with the RNA-Seq data of 50 tissues. The functions of 97 genes were predicted, revealing 225 pathways. The mastitis proteins were used for protein-protein interaction network analysis to obtain additional cross-talking proteins. A total of 1,306 SNPs and 152 indels were identified from 101 genes. Water Buffalo-MSTdb was developed with 3-tier architecture to retrieve mastitis associated genes having genomic coordinates with chromosomal details for TGP sequencing for mining of minor alleles for further association studies. Lastly, a web-genomic resource was made available to mine variants of targeted gene panels in buffalo for mastitis resistance breeding in an endeavor to ensure improved productivity and the reproductive efficiency of water buffalo.

Somatic cell score: gene polymorphisms and other effects in Holstein and Simmental cows

Objective

The aim of the study was to evaluate the influence of gene polymorphisms and nongenetic factors on the somatic cell score (SCS) in the milk of Holstein (n = 148) and Simmental (n = 73) cows and their crosses (n = 6).

Methods

The SCS was calculated by the formula SCS = log₂(SCC/100,000)+3, where SCC is the somatic cell count. Polymorphisms in the casein alpha S1 (CSN1S1), beta-casein (CSN2), kappa-casein (CSN3), beta-lactoglobulin (LGB), acyl-CoA diacylglycerol transferase 1 (DGAT1), leptin (LEP), fatty acid synthase (FASN), stearoyl CoA desaturase 1 (SCD1), and 1-acylglycerol-3-phosphate O-acyltransferase 6 (AGPAT6) genes were genotyped, and association analysis to the SCS in the cow’s milk was performed. Further, the impact of breed, farm, year, month of the year, lactation stage and parity on the SCS were analysed. Phenotype correlations among SCS and milk constituents were computed by Pearson correlation coefficients.

Results

Only CSN2 genotypes A¹/A² were found to have significant association with the SCS (p<0.05), and alleles of CSN1S1 and DGAT1 genes (p<0.05). Other polymorphisms were not found to be significant. SCS had significant association with the combined effect of farm and year, lactation stage and month of the year. Lactation parity and breed had not significant association with SCS. The phenotypic correlation of SCS to lactose content was negative and significant, while the correlation to protein content was positive and significant. The correlations of SCS to fat, casein, nonfat solids, urea, citric acid, acetone and ketones contents were very low and not significant.

Conclusion

Only CSN2 genotypes, CSN1S1 and DGAT1 alleles did show an obvious association to the SCS. The results confirmed the importance of general quality management of farms on the microbial milk quality, and effects of lactation stage and month of the year. The lactose content in milk reflects the health status of the udder.

Role of ACE2 genetic polymorphisms in susceptibility to SARS-CoV-2 among highly exposed but non infected healthcare workers

We aim to evaluate the role of single-nucleotide polymorphisms of the angiotensin-converting enzyme 2 in susceptibility to SARS-CoV-2 infection. We included 28 uninfected but highly exposed healthcare workers and 39 hospitalized patients with COVID-19. Thirty-five SNPs were rationally selected. Two variants were associated with increased risk of being susceptible to SARS-CoV-2: the minor A allele in the rs2106806 variant (OR 3.75 [95% CI 1.23-11.43]) and the minor T allele in the rs6629110 variant (OR 3.39 [95% CI 1.09-10.56]). Evaluating the role of genetic variants in susceptibility to SARS-CoV-2 infection could help identify more vulnerable individuals and suggest potential drug targets for COVID-19 patients.

Impact of PSCA Polymorphisms on the Risk of Duodenal Ulcer

BACKGROUND

While duodenal ulcer (DU) and gastric cancer (GC) are both H. pylori infection-related diseases, individuals with DU are known to have lower risk for GC. Many epidemiological studies have identified the PSCA rs2294008 T-allele as a risk factor of GC, while others have found an association between the rs2294008 C-allele and risk of DU and gastric ulcer (GU). Following these initial reports, however, few studies have since validated these associations. Here, we aimed to validate the association between variations in PSCA and the risk of DU/GU and evaluate its interaction with environmental factors in a Japanese population.

METHODS

Six PSCA SNPs were genotyped in 584 DU cases, 925 GU cases, and 8,105 controls from the Japan Multi-Institutional Collaborative Cohort (J-MICC). Unconditional logistic regression models were applied to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between the SNPs and risk of DU/GU.

RESULTS

PSCA rs2294008 C-allele was associated with per allele OR of 1.34 (95% CI, 1.18-1.51; P = 2.28 × 10-6) for the risk of DU. This association was independent of age, sex, study site, smoking habit, drinking habit, and H. pylori status. On the other hand, we did not observe an association between the risk of GU and PSCA SNPs.

CONCLUSIONS

Our study confirms an association between the PSCA rs2294008 C-allele and the risk of DU in a Japanese population.

SNPs from BCHE and DRD3 genes associated to cocaine abuse amongst violent individuals from Sao Paulo, Brazil

Violence and drug abuse are highly destructive phenomena found world-wide, especially in Brazil. They seem to rise proportionally to one another and possibly related. Additionally, genetics may also play a role in drug abuse. This study has focused on identifying the use of cocaine within postmortem cases arriving at the Institute of Legal Medicine of Sao Paulo as well as the presence of certain single nucleotide polymorphisms (SNPs) to better understand one's susceptibility to abuse the drug. Both hair and blood samples have been extracted through a simple methanol overnight incubation or a rapid dilute-and-shoot method, respectively. The samples were then analyzed using an UPLC-ESI-MS/MS and genotyped through RT-PCR. Statistical analyses were performed via SPSS software. From 105 postmortem cases, 53% and 51% of the cases shown to be positive for cocaine in hair and blood, respectively. Genetic wise, a significant difference has been observed for SNP rs4263329 from the BCHE gene with higher frequencies of the genotypes A/G and G/G seen in cocaine users (OR=8.91; 95%CI=1.58-50.21; p=0.01). Likewise, also SNP rs6280 from the DRD3 gene presented a significant association, with both genotypes T/C and C/C being more frequent in users (OR=4.96; 95% CI=1.07-23.02; p=0.04). To conclude, a rather high proportion of cocaine has been found, which may suggest a connotation between the use of the drug and risky/violent behaviors. Additionally, significant associations were also found within two SNPs related to cocaine use, however, due to several inherent limitations, these must be confirmed.

Identification of single-nucleotide variants associated with susceptibility to Salmonella in pigs using a genome-wide association approach

BACKGROUND

Salmonella enterica serovars are a major cause of foodborne illness and have a substantial impact on global human health. In Canada, Salmonella is commonly found on swine farms and the increasing concern about drug use and antimicrobial resistance associated with Salmonella has promoted research into alternative control methods, including selecting for pig genotypes associated with resistance to Salmonella. The objective of this study was to identify single-nucleotide variants in the pig genome associated with Salmonella susceptibility using a genome-wide association approach. Repeated blood and fecal samples were collected from 809 pigs in 14 groups on farms and tonsils and lymph nodes were collected at slaughter. Sera were analyzed for Salmonella IgG antibodies by ELISA and feces and tissues were cultured for Salmonella. Pig DNA was genotyped using a custom 54 K single-nucleotide variant oligo array and logistic mixed-models used to identify SNVs associated with IgG seropositivity, shedding, and tissue colonization.

RESULTS

Variants in/near PTPRJ (p = 0.0000066), ST6GALNAC3 (p = 0.0000099), and DCDC2C (n = 3, p < 0.0000086) were associated with susceptibility to Salmonella, while variants near AKAP12 (n = 3, p < 0.0000358) and in RALGAPA2 (p = 0.0000760) may be associated with susceptibility.

CONCLUSIONS

Further study of the variants and genes identified may improve our understanding of neutrophil recruitment, intracellular killing of bacteria, and/or susceptibility to Salmonella and may help future efforts to reduce Salmonella on-farm through genetic approaches.

Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment

Methotrexate (MTX) is one of the staples of pediatric acute lymphoblastic leukemia (ALL) treatment. MTX targets the folate metabolic pathway (FMP). Abnormal function of the enzymes in FMP, due to genetic aberrations, leads to adverse drug reactions. The aim of this study was to investigate variants in pharmacogenes involved in FMP and their association with MTX pharmacokinetics (MTX elimination profile) and toxicity in the consolidation therapy phase of pediatric ALL patients. Eleven variants in the thymidylate synthetase (TYMS), methylenetetrahydrofolate reductase (MTHFR), dihydrofolate reductase (DHFR), SLC19A1 and SLCO1B genes were analyzed in 148 patients, using PCR- and sequencing-based methodology. For the Serbian and European control groups, data on allele frequency distribution were extracted from in-house and public databases. Our results show that the A allele of SLC19A1 c.80 variant contributes to slow MTX elimination. Additionally, the AA genotype of the same variant is a predictor of MTX-related hepatotoxicity. Patients homozygous for TYMS 6bp deletion were more likely to experience gastrointestinal toxicity. No allele frequency dissimilarity was found for the analyzed variants between Serbian and European populations. Statistical modelling did not show a joint effect of analyzed variants. Our results indicate that SLC19A1 c.80 variant and TYMS 6bp deletion are the most promising pharmacogenomic markers of MTX response in pediatric ALL patients.

Allelic variants in the zinc transporter-3 gene, SLC30A3, a candidate gene identified from gene expression studies, show gender-specific association with schizophrenia

Previous microarray analysis of gene expression in frontal cortex showed differential expression of genes associated with synaptic function in schizophrenia compared to matched-controls in two independent cohorts. One of these genes validated in both cohorts, SLC30A3, which encodes the Zinc Transporter 3 (ZNT3), is localised to synaptic vesicles in glutamate synapses and known to be involved in cognitive function. In view of the robust depletion of SLC30A3 mRNA in two independent studies and the importance of this gene in cognitive function, we investigated whether single nucleotide polymorphism (SNP) associations with schizophrenia could be detected in a UK case controlled schizophrenia cohort. Four SNPs were selected across this gene and genotyped in a cohort of cases and controls from East UK. We found significant associations with schizophrenia at the allelic (ORs: 1.51 to 1.57), genotype (ORs: 1.46 to 1.53) and haplotype level (P = 2.15 × 10−4). These associations proved to be gender-specific with significant effects of allele (ORs: 1.74 to 2.11), genotype (ORs: 1.78 to 2.14) and haplotype (P = 3.51 × 10−5) observed in female schizophrenia cases but not males, when split by gender. In conclusion, SNPs in SLC30A3 showed a gender-specific association with schizophrenia in this East UK cohort, which merits further investigation in other population samples.

Polymorphisms in Fatty Acid Desaturase 2 Gene are Associated with Milk Production Traits in Chinese Holstein Cows

Simple Summary

Searching for causative polymorphisms underlying the variability of milk production traits and then incorporating them into breeding programs are very effective ways to improve the efficiency and reliability of conventional dairy cattle breeding. Fatty acid desaturase 2 (FADS2) plays a pivotal role in the biosynthesis of polyunsaturated fatty acids. Previous studies provided evidence that FADS2 was one of the most downregulated genes during negative energy balance in the liver of postpartum dairy cattle. Genes involved in the energetic pathways may influence other production traits, such as protein, fat and milk yields. Therefore, in the present study, we investigated the common genetic variants of the FADS2 gene in Chinese Holstein cows. Our results provided direct evidence that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Abstract

This study investigated the single nucleotide polymorphisms (SNPs) of Fatty acid desaturase 2 (FADS2) gene and further explored their genetic effects on conventionally collected milk production traits in Chinese Holstein cows using 18,264 test-day records of 841 cows. One missense mutation c. 908 C > T (SNP site in the complementary DNA sequence), which caused an amino acid change from alanine to valine (294Ala > Val), and two 3’ untranslated region (UTR) SNPs, c.1571 G > A and c.2776 A > G were finally identified. The SNP c.908 C > T was significantly associated with test-day milk yield, fat percentage and 305-day milk, fat and protein yield. In particular, the T allele of the SNP c.908 C > T showed a significant association with decreased somatic cell score (SCS) in the investigated population. Significant relationship between the SNP c.1571 G > A and 305-day milk yield showed that genotype GG was linked to the highest milk yield. Substituting the allele G for A at the c.2776 A > G locus resulted in a decrease of protein percentage. Our results demonstrated that FADS2 was an interesting candidate for selection to increase milk production and improve resistance against mastitis.

Combining genome-wide association study and FST -based approaches to identify targets of Borrelia-mediated selection in natural rodent hosts

Recent advances in high-throughput sequencing technologies provide opportunities to gain novel insights into the genetic basis of phenotypic trait variation. Yet to date, progress in our understanding of genotype-phenotype associations in nonmodel organisms in general and natural vertebrate populations in particular has been hampered by small sample sizes typically available for wildlife populations and a resulting lack of statistical power, as well as a limited ability to control for false-positive signals. Here we propose to combine a genome-wide association study (GWAS) and F_ST -based approach with population-level replication to partly overcome these limitations. We present a case study in which we used this approach in combination with genotyping-by-sequencing (GBS) single nucleotide polymorphism (SNP) data to identify genomic regions associated with Borrelia afzelii resistance or susceptibility in the natural rodent host of this Lyme disease-causing spirochete, the bank vole (Myodes glareolus). Using this combined approach we identified four consensus SNPs located in exonic regions of the genes Slc26a4, Tns3, Wscd1 and Espnl, which were significantly associated with the voles' Borrelia infectious status within and across populations. Functional links between host responses to bacterial infections and most of these genes have previously been demonstrated in other rodent systems, making them promising new candidates for the study of evolutionary host responses to Borrelia emergence. Our approach is applicable to other systems and may facilitate the identification of genetic variants underlying disease resistance or susceptibility, as well as other ecologically relevant traits, in wildlife populations.

Genetic variability of serotonin pathway associated with schizophrenia onset, progression, and treatment

Schizophrenia (SZ) onset and treatment outcome have important genetic components, however individual genes do not have strong effects on SZ phenotype. Therefore, it is important to use the pathway-based approach and study metabolic and signaling pathways, such as dopaminergic and serotonergic. Serotonin pathway has an important role in brain signaling, nevertheless, its role in SZ is not as thoroughly examined as that of dopamine pathway. In this study, we reviewed serotonin pathway genes and genetic variations associated with SZ, including variations at DNA, RNA, and epigenetic level. We obtained 30 serotonin pathway genes from Kyoto encyclopedia of genes and genomes and used these genes for the literature review. We extracted 20 protein coding serotonin pathway genes with genetic variations associated with SZ onset, development, and treatment from 31 research papers. Genes associated with SZ are present on all levels of serotonin pathway: serotonin synthesis, transport, receptor binding, intracellular signaling, and reuptake; however, regulatory genes are poorly researched. We summarized common challenges of genetic association studies and presented some solutions. The analysis of reported serotonin pathway-SZ associations revealed lack of information about certain serotonin pathway genes potentially associated with SZ. Furthermore, it is becoming clear that interactions among serotonin pathway genes and their regulators may bring further knowledge about their involvement in SZ.

Association mapping of wheat Fusarium head blight resistance-related regions using a candidate-gene approach and their verification in a biparental population

Markers, located in Dicer1 and Ara6 genes, which are likely involved in cross-kingdom RNA trafficking, are associated with FHB resistance in GABI wheat population and were validated in biparental population. Association studies are a common approach to detect marker-trait associations for Fusarium head blight (FHB) resistance in wheat (Triticum aestivum), although verification of detected associations is exceptional. In the present study, candidate-gene association mapping (CG) of genes from silencing and secretory pathways, which may be involved in wheat resistance against FHB and cross-kingdom RNA trafficking, was performed. Fourteen markers, located in nine genes, were tested for association with FHB resistance in 356 lines from the GABI (genome analysis of the biological system of plants) wheat population. Three markers located in the genes Dicer1 and Ara6 were shown to be significantly associated with the studied trait. Verification of this finding was performed using the recombinant inbred lines (RILs) population 'Apache × Biscay', segregating for four of our 14 selected markers. We could show association of the Ara6 marker with plant height as well as association with FHB resistance for three markers located in Rab5-like GTPase gene Ara6 and Dicer1. These results confirmed the trait-marker associations detected also in the CG approach. Gene products of the associated genes are involved in response of the plant to pathogens, plant metabolism and may be involved in cross-kingdom RNA trafficking efficiency. The markers detected in the GABI wheat population, which were also validated in the biparental population, can potentially be used in wheat breeding.

Systematic Review and Meta-Analysis of Genetic Risk of Developing Chronic Postsurgical Pain

Chronic postsurgical pain (CPSP) is a significant detriment to postsurgical recovery and a risk factor for prolonged opioid use. Emerging evidence suggests the estimated heritability for chronic pain is 45% and that genetic factors partially explain individual susceptibility to CPSP. The aim of this study was to systematically review, assess quality, and summarize the studies in humans that have investigated genetic factors associated with CPSP. We also conducted a meta-analysis to derive a single effect size for evaluable genetic associations with CPSP. Our comprehensive literature search included review of 21 full-text articles evaluating variants of 69 genes for association with CPSP. We found significant gene variant associations reported for variants/haplotypes of 26 genes involved in neurotransmission, pain signaling, immune responses and neuroactive ligand-receptor interaction, with CPSP. Six variants of 5 genes (COMT: rs4680 and rs6269, OPRM1: rs1799971, GCH1: rs3783641, KCNS1: rs734784 and TNFA: rs1800629), were evaluated by more than one study and were included in the meta-analysis. At rs734784 (A>G) of KCNS1, presence of G allele marginally increased risk of CPSP (Additive genetic model; Odds ratio: 1.511; 95% CI 1-2.284; P value: .050), while the other variants did not withstand meta-analyses criteria. Our findings demonstrate the role of genetic factors with different functions in CPSP, and also emphasize that single genetic factors have small effect sizes in explaining complex conditions like CPSP. Heterogeneity in surgical cohorts, population structure, and outcome definitions, as well as small number of available studies evaluating same variants, limit the meta-analysis. There is a need for large-scale, homogenous, replication studies to validate candidate genes, and understand the underlying biological networks underpinning CPSP. PERSPECTIVE: Our systematic review comprehensively describes 21 studies evaluating genetic association with CPSP, and limitations thereof. A meta-analysis of 6 variants (5 genes) found marginally increased risk for CPSP associated with rs734784 A>G of the potassium voltage-gated channel gene (KCNS1). Understanding genetic predisposition for CPSP will enable prediction and personalized management.

Pharmacogenomic and Pharmacotranscriptomic Profiling of Childhood Acute Lymphoblastic Leukemia: Paving the Way to Personalized Treatment

Personalized medicine is focused on research disciplines which contribute to the individualization of therapy, like pharmacogenomics and pharmacotranscriptomics. Acute lymphoblastic leukemia (ALL) is the most common malignancy of childhood. It is one of the pediatric malignancies with the highest cure rate, but still a lethal outcome due to therapy accounts for 1%⁻3% of deaths. Further improvement of treatment protocols is needed through the implementation of pharmacogenomics and pharmacotranscriptomics. Emerging high-throughput technologies, including microarrays and next-generation sequencing, have provided an enormous amount of molecular data with the potential to be implemented in childhood ALL treatment protocols. In the current review, we summarized the contribution of these novel technologies to the pharmacogenomics and pharmacotranscriptomics of childhood ALL. We have presented data on molecular markers responsible for the efficacy, side effects, and toxicity of the drugs commonly used for childhood ALL treatment, i.e., glucocorticoids, vincristine, asparaginase, anthracyclines, thiopurines, and methotrexate. Big data was generated using high-throughput technologies, but their implementation in clinical practice is poor. Research efforts should be focused on data analysis and designing prediction models using machine learning algorithms. Bioinformatics tools and the implementation of artificial i Lack of association of the CEP72 rs924607 TT genotype with intelligence are expected to open the door wide for personalized medicine in the clinical practice of childhood ALL.

Genetic associations of hemoglobin in children with chronic kidney disease in the PediGFR Consortium

BACKGROUND

Genome-wide association studies (GWAS) in healthy populations have identified variants associated with erythrocyte traits, but genetic causes of hemoglobin variation in children with CKD are incompletely understood.

METHODS

The Pediatric Investigation of Genetic Factors Linked with Renal Progression (PediGFR) Consortium comprises three pediatric CKD cohorts: Chronic Kidney Disease in Children (CKiD), Effect of Strict Blood Pressure Control and ACE Inhibition on the Progression of CRF in Pediatric Patients (ESCAPE), and Cardiovascular Comorbidity in Children with CKD (4C). We performed cross-sectional and longitudinal association studies of single-nucleotide polymorphisms (SNPs) in 1125 patients.

RESULTS

Children of European (n = 725) or Turkish (n = 400) ancestry (EA or TA) were included. In cross-sectional analysis, two SNPs (rs10758658 and rs12718597) previously associated with RBC traits were significantly associated with hemoglobin levels in children of EA and TA. In longitudinal analysis, SNP rs2540917 was nominally associated with hemoglobin in EA and TA children.

CONCLUSIONS

SNPs associated with erythrocyte traits in healthy populations were marginally significant for an association with hemoglobin. Further analyses/replication studies are needed in larger CKD cohorts to investigate SNPs of unknown significance associated with hemoglobin. Functional studies will be required to confirm that the observed associations between SNPs and clinical phenotype are causal.

Gamete compatibility genes in mammals: candidates, applications and a potential path forward

Fertilization represents a critical stage in biology, where successful alleles of a previous generation are shuffled into new arrangements and subjected to the forces of selection in the next generation. Although much research has been conducted on how variation in morphological and behavioural traits lead to variation in fertilization patterns, surprisingly little is known about fertilization at a molecular level, and specifically about how genes expressed on the sperm and egg themselves influence fertilization patterns. In mammals, several genes have been identified whose products are expressed on either the sperm or the egg, and which influence the fertilization process, but the specific mechanisms are not yet known. Additionally, in 2014 an interacting pair of proteins was identified: 'Izumo' on the sperm, and 'Juno' on the egg. With the identification of these genes comes the first opportunity to understand the molecular aspects of fertilization in mammals, and to identify how the genetic characteristics of these genes influence fertilization patterns. Here, we review recent progress in our understanding of fertilization and gamete compatibility in mammals, which should provide a helpful guide to researchers interested in untangling the molecular mechanisms of fertilization and the resulting impacts on population biology and evolutionary processes.

Genetic factors associated with exercise performance in atmospheric hypoxia

Background and Objective

‘Natural selection’ has been shown to have enriched the genomes of high-altitude native populations with genetic variants of advantage in this hostile hypoxic environment. In lowlanders who ascend to altitude, genetic factors may also contribute to the substantial interindividual variation in exercise performance noted at altitude. We performed a systematic literature review to identify genetic variants of possible influence on human hypoxic exercise performance, commenting on the strength of any identified associations.

Criteria for considering studies for this review

All studies of the association of genetic factors with human hypoxic exercise performance, whether at sea level using ‘nitrogen dilution of oxygen’ (normobaric hypoxia), or at altitude or in low-pressure chambers (field or chamber hypobaric hypoxia, respectively) were sought for review.

Search strategy for identification of studies

Two electronic databases were searched (Ovid MEDLINE, Embase) up to 31 January 2014. We also searched the reference lists of relevant articles for eligible studies. All studies published in English were included, as were studies in any language for which the abstract was available in English.

Data collection and analysis

Studies were selected and data extracted independently by two reviewers. Differences regarding study inclusion were resolved through discussion. The quality of each study was assessed using a scoring system based on published guidelines for conducting and reporting genetic association studies.

Results

A total of 11 studies met all inclusion criteria and were included in the review. Subject numbers ranged from 20 to 1,931 and consisted of healthy individuals in all cases. The maximum altitude of exposure ranged from 2,690 to 8,848 m. The exercise performance phenotypes assessed were mountaineering performance (n = 5), running performance (n = 2), and maximum oxygen consumption (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \dot{V} $$\end{document}V˙O₂max) (n = 4). In total, 13 genetic polymorphisms were studied, four of which were associated with hypoxic exercise performance. The adenosine monophosphate deaminase (AMPD1) C34T (rs17602729), beta2-adrenergic receptor (ADRB2) Gly16Arg single nucleotide polymorphism (SNP) (rs1042713), and androgen receptor CAG repeat polymorphisms were associated with altitude performance in one study, and the angiotensin I-converting enzyme (ACE) insertion/deletion (I/D) (rs4646994) polymorphism was associated with performance in three studies. The median score achieved in the study quality analysis was 6 out of 10 for case–control studies, 8 out of 10 for cohort studies with a discrete outcome, 6 out of 9 for cohort studies with a continuous outcome, and 4.5 out of 8 for genetic admixture studies.

Conclusion

The small number of articles identified in the current review and the limited number of polymorphisms studied in total highlights that the influence of genetic factors on exercise performance in hypoxia has not been studied in depth, which precludes firm conclusions being drawn. Support for the association between the ACE-I allele and improved high-altitude performance was the strongest, with three studies identifying a relationship. Analysis of study quality highlights the need for future studies in this field to improve the conduct and reporting of genetic association studies.

Electronic supplementary material

The online version of this article (doi:10.1007/s40279-015-0309-8) contains supplementary material, which is available to authorized users.

Immunogenetics of IgG4-Related AIP

Autoimmune pancreatitis (AIP) is a unique form of chronic pancreatitis characterized by high serum IgG4 concentration and a variety of complicating extra-pancreatic lesions. AIP has the features of a complex disease that is caused by multifactorial genes. However, the genetic factors underlying AIP have not been elucidated conclusively. Association studies by the candidate-gene approach and genome-wide association studies (GWAS) have revealed several susceptibility genes for AIP, including HLA DRB1*04:05-DQB1*04:01, FCRL3, CTLA4, and KCNA3, albeit in small-scale analyses. Thus, GWAS of large sample sizes and multinational collaborative meta-analyses are needed to identify the precise genetic variants that are associated with AIP onset. Systems genetics approaches that integrate DNA sequencing, expression quantitative trait locus (eQTL) mapping, proteomics, and metabolomics will also be useful in clarifying the pathogenesis of AIP.

Renin-Angiotensin System Gene Variants and Type 2 Diabetes Mellitus: Influence of Angiotensinogen

Genome-wide association studies (GWAS) have been successfully used to call for variants associated with diseases including type 2 diabetes mellitus (T2DM). However, some variants are not included in the GWAS to avoid penalty in multiple hypothetic testing. Thus, candidate gene approach is still useful even at GWAS era. This study attempted to assess whether genetic variations in the renin-angiotensin system (RAS) and their gene interactions are associated with T2DM risk. We genotyped 290 T2DM patients and 267 controls using three genes of the RAS, namely, angiotensin converting enzyme (ACE), angiotensinogen (AGT), and angiotensin II type 1 receptor (AGTR1). There were significant differences in allele frequencies between cases and controls for AGT variants (P = 0.05) but not for ACE and AGTR1. Haplotype TCG of the AGT was associated with increased risk of T2DM (OR 1.92, 95% CI 1.15-3.20, permuted P = 0.012); however, no evidence of significant gene-gene interactions was seen. Nonetheless, our analysis revealed that the associations of the AGT variants with T2DM were independently associated. Thus, this study suggests that genetic variants of the RAS can modestly influence the T2DM risk.

Genome-wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco-evolutionary processes in natural populations. Host-parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade-offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome-wide SNPs. We first use genome-wide association to identify individual SNPs associated with nematode burden. We then partition genome-wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome-wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large-effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).

Inherited variation in immune response genes in follicular lymphoma and diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) both depend on immune-mediated survival and proliferation signals from the tumor microenvironment. Inherited genetic variation influences this complex interaction. A total of 89 studies investigating immune-response genes in DLBCL and FL were critically reviewed. Relatively consistent association exists for variation in the tumor necrosis factor alpha (TNFA) and interleukin-10 loci and DLBCL risk; for DLBCL outcome association with the TNFA locus exists. Variations at chromosome 6p31-32 were associated with FL risk. Importantly, individual risk alleles have been shown to interact with each other. We suggest that the pathogenetic impact of polymorphic genes should include gene-gene interaction analysis and should be validated in preclinical model systems of normal B lymphopoiesis and B-cell malignancies. In the future, large cohort studies of interactions and genome-wide association studies are needed to extend the present findings and explore new risk alleles to be studied in preclinical models.

Developmental programming of hypothalamic neuronal circuits: impact on energy balance control

The prevalence of obesity in adults and children has increased globally at an alarming rate. Mounting evidence from both epidemiological studies and animal models indicates that adult obesity and associated metabolic disorders can be programmed by intrauterine and early postnatal environment- a phenomenon known as "fetal programming of adult disease." Data from nutritional intervention studies in animals including maternal under- and over-nutrition support the developmental origins of obesity and metabolic syndrome. The hypothalamic neuronal circuits located in the arcuate nucleus controlling appetite and energy expenditure are set early in life and are perturbed by maternal nutritional insults. In this review, we focus on the effects of maternal nutrition in programming permanent changes in these hypothalamic circuits, with experimental evidence from animal models of maternal under- and over-nutrition. We discuss the epigenetic modifications which regulate hypothalamic gene expression as potential molecular mechanisms linking maternal diet during pregnancy to the offspring's risk of obesity at a later age. Understanding these mechanisms in key metabolic genes may provide insights into the development of preventative intervention strategies.

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD)

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined P_ADJ = 1.15 × 10⁻⁶, OR 1.332 (1.187–1.496)]. This association was characterized by a highly significant sex difference (P_diff = 0.0032) in that it was clearly confined to females with genome-wide significance [P_ADJ = 2.62 × 10⁻⁸, OR 1.541 (1.324–1.796); males: P_ADJ = 0.382, OR 1.084 (0.905–1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMD pathogenesis.

COMT met allele differentially predicts risk versus severity of aberrant eating in a large community sample

Prefrontal dopamine (DA) transmission participates in the reinforcement of reward-driven behaviors like eating. Because catechol-O-methyltransferase (COMT) degrades DA and is expressed in the prefrontal cortex, variation in the COMT gene may modulate eating behavior. Previous studies have shown that the met allele of the COMT val158met single nucleotide polymorphism (SNP) is associated with Bulimia Nervosa (BN). The specific aim of this study was to test whether the met allele increased risk for, and severity of, eating disorder symptomatology in community volunteers. Caucasian adults (N=1003; 51.2% female) from the University of Pittsburgh Adult Health and Behavior Project (AHAB) were genotyped and completed the Eating Disorders Inventory (EDI). Logistic and Poisson regression analyses assessed genotype-dependent presence and severity of eating disorder symptomatology. The met allele was significantly associated with the presence of symptoms on the Bulimia subscale and the severity of Body Dissatisfaction scores. All EDI subscales significantly predicted BMI. To our knowledge, these are the first data indicating that the COMT met allele increases risk for some symptoms of disordered eating, while increasing severity of others, in a community sample. These novel findings may have important implications for understanding the etiology of heterogeneous disordered eating phenotypes.

Functional-mixed effects models for candidate genetic mapping in imaging genetic studies

The aim of this paper is to develop a functional-mixed effects modeling (FMEM) framework for the joint analysis of high-dimensional imaging data in a large number of locations (called voxels) of a three-dimensional volume with a set of genetic markers and clinical covariates. Our FMEM is extremely useful for efficiently carrying out the candidate gene approaches in imaging genetic studies. FMEM consists of two novel components including a mixed effects model for modeling nonlinear genetic effects on imaging phenotypes by introducing the genetic random effects at each voxel and a jumping surface model for modeling the variance components of the genetic random effects and fixed effects as piecewise smooth functions of the voxels. Moreover, FMEM naturally accommodates the correlation structure of the genetic markers at each voxel, while the jumping surface model explicitly incorporates the intrinsically spatial smoothness of the imaging data. We propose a novel two-stage adaptive smoothing procedure to spatially estimate the piecewise smooth functions, particularly the irregular functional genetic variance components, while preserving their edges among different piecewise-smooth regions. We develop weighted likelihood ratio tests and derive their exact approximations to test the effect of the genetic markers across voxels. Simulation studies show that FMEM significantly outperforms voxel-wise approaches in terms of higher sensitivity and specificity to identify regions of interest for carrying out candidate genetic mapping in imaging genetic studies. Finally, FMEM is used to identify brain regions affected by three candidate genes including CR1, CD2AP, and PICALM, thereby hoping to shed light on the pathological interactions between these candidate genes and brain structure and function.

Genetic variation in DNA repair pathways and risk of non-Hodgkin's lymphoma

Molecular and genetic evidence suggests that DNA repair pathways may contribute to lymphoma susceptibility. Several studies have examined the association of DNA repair genes with lymphoma risk, but the findings from these reports have been inconsistent. Here we provide the results of a focused analysis of genetic variation in DNA repair genes and their association with the risk of non-Hodgkin's lymphoma (NHL). With a population of 1,297 NHL cases and 1,946 controls, we have performed a two-stage case/control association analysis of 446 single nucleotide polymorphisms (SNPs) tagging the genetic variation in 81 DNA repair genes. We found the most significant association with NHL risk in the ATM locus for rs227060 (OR = 1.27, 95% CI: 1.13-1.43, p = 6.77×10(-5)), which remained significant after adjustment for multiple testing. In a subtype-specific analysis, associations were also observed for the ATM locus among both diffuse large B-cell lymphomas (DLBCL) and small lymphocytic lymphomas (SLL), however there was no association observed among follicular lymphomas (FL). In addition, our study provides suggestive evidence of an interaction between SNPs in MRE11A and NBS1 associated with NHL risk (OR = 0.51, 95% CI: 0.34-0.77, p = 0.0002). Finally, an imputation analysis using the 1,000 Genomes Project data combined with a functional prediction analysis revealed the presence of biologically relevant variants that correlate with the observed association signals. While the findings generated here warrant independent validation, the results of our large study suggest that ATM may be a novel locus associated with the risk of multiple subtypes of NHL.

Cancer stem cell marker Musashi-1 rs2522137 genotype is associated with an increased risk of lung cancer

Gene single nucleotide polymorphisms (SNPs) have been extensively studied in association with development and prognosis of various malignancies. However, the potential role of genetic polymorphisms of cancer stem cell (CSC) marker genes with respect to cancer risk has not been examined. We conducted a case-control study involving a total of 1000 subjects (500 lung cancer patients and 500 age-matched cancer-free controls) from northeastern China. Lung cancer risk was analyzed in a logistic regression model in association with genotypes of four lung CSC marker genes (CD133, ALDH1, Musashi-1, and EpCAM). Using univariate analysis, the Musashi-1 rs2522137 GG genotype was found to be associated with a higher incidence of lung cancer compared with the TT genotype. No significant associations were observed for gene variants of CD133, ALDH1, or EpCAM. In multivariate analysis, Musashi-1 rs2522137 was still significantly associated with lung cancer when environmental and lifestyle factors were incorporated in the model, including lower BMI; family history of cancer; prior diagnosis of chronic obstructive pulmonary disease, pneumonia, or pulmonary tuberculosis; occupational exposure to pesticide; occupational exposure to gasoline or diesel fuel; heavier smoking; and exposure to heavy cooking emissions. The value of the area under the receiver-operating characteristic (ROC) curve (AUC) was 0.7686. To our knowledge, this is the first report to show an association between a Musashi-1 genotype and lung cancer risk. Further, the prediction model in this study may be useful in determining individuals with high risk of lung cancer.

Challenges and prospects in genome-wide quantitative trait loci mapping of standing genetic variation in natural populations

A considerable challenge in evolutionary genetics is to understand the genetic mechanisms that facilitate or impede evolutionary adaptation in natural populations. For this, we must understand the genetic loci contributing to trait variation and the selective forces acting on them. The decreased costs and increased feasibility of obtaining genotypic data on a large number of individuals have greatly facilitated gene mapping in natural populations, particularly because organisms whose genetics have been historically difficult to study are now within reach. Here we review the methods available to evolutionary ecologists interested in dissecting the genetic basis of traits in natural populations. Our focus lies on standing genetic variation in outbred populations. We present an overview of the current state of research in the field, covering studies on both plants and animals. We also draw attention to particular challenges associated with the discovery of quantitative trait loci and discuss parallels to studies on crops, livestock, and humans. Finally, we point to some likely future developments in genetic mapping studies.

Lymphoblastoid cell lines models of drug response: successes and lessons from this pharmacogenomic model

A new standard for medicine is emerging that aims to improve individual drug responses through studying associations with genetic variations. This field, pharmacogenomics, is undergoing a rapid expansion due to a variety of technological advancements that are enabling higher throughput with reductions in cost. Here we review the advantages, limitations, and opportunities for using lymphoblastoid cell lines (LCL) as a model system for human pharmacogenomic studies. There are a wide range of publicly available resources with genome-wide data available for LCLs from both related and unrelated populations, removing the cost of genotyping the data for drug response studies. Furthermore, in contrast to human clinical trials or in vivo model systems, with high-throughput in vitro screening technologies, pharmacogenomics studies can easily be scaled to accommodate large sample sizes. An important component to leveraging genome-wide data in LCL models is association mapping. Several methods are discussed herein, and include multivariate concentration response modeling, issues with multiple testing, and successful examples of the 'triangle model' to identify candidate variants. Once candidate gene variants have been determined, their biological roles can be elucidated using pathway analyses and functionally confirmed using siRNA knockdown experiments. The wealth of genomics data being produced using related and unrelated populations is creating many exciting opportunities leading to new insights into the genetic contribution and heritability of drug response.

Imaging genetics in adult attention-deficit/hyperactivity disorder (ADHD): a way towards pathophysiological understanding?

Attention-deficit/hyperactivity disorder (ADHD) is a common, early-onset and enduring developmental disorder whose underlying etiological and neurobiological processes are the current focus of major research. Research strategies have made considerable effort in elucidating the complex genetic architecture of ADHD and indicate various pathways from genotype to phenotype. Understanding ADHD as a neuropsychiatric disorder enabled to investigate markers of neural activity as endophenotypes to better explain the link from gene to symptomatology (the so-called imaging genetics approach). Overcoming the originally rather restrictive requirements for an endophenotype, imaging genetics studies are supposed to offer a much more flexible and hypothesis-driven approach towards the etiology of ADHD. Although 1) ADHD often persists into adulthood, thus remaining a prevalent disorder, and 2) imaging genetics provides a promising research approach, a review on imaging genetics in adult ADHD - as available for childhood ADHD (Durston 2010) - is lacking. In this review, therefore, findings from the few available imaging genetics studies in adult ADHD will be summarized and complemented by relevant findings from healthy controls and children with ADHD that are considered important for the adult ADHD imaging genetics approach. The studies will be reviewed regarding implications for basic research and possible practical applications. Imaging genetics studies in adult ADHD have the potential to further clarify pathophysiological pathways and mechanisms, to derive new testable hypotheses, to investigate genetic interaction effects and to partly influence practical applications. In combination with other research strategies, they can incrementally foster the understanding of relevant processes in a more comprehensive way. Current limitations comprise the incapability to discover new genes, a high genetic load in patients potentially obscuring the effect of single candidate genes, the mostly unknown heritability of the endophenotype and the heterogeneous manifestation of ADHD.

Maternal undernutrition programs tissue-specific epigenetic changes in the glucocorticoid receptor in adult offspring

Epidemiological data indicate that an adverse maternal environment during pregnancy predisposes offspring to metabolic syndrome with increased obesity, and type 2 diabetes. The mechanisms are still unclear although epigenetic modifications are implicated and the hypothalamus is a likely target. We hypothesized that maternal undernutrition (UN) around conception in sheep would lead to epigenetic changes in hypothalamic neurons regulating energy balance in the offspring, up to 5 years after the maternal insult. We found striking evidence of decreased glucocorticoid receptor (GR) promoter methylation, decreased histone lysine 27 trimethylation, and increased histone H3 lysine 9 acetylation in hypothalami from male and female adult offspring of UN mothers. These findings are entirely compatible with the increased GR mRNA and protein observed in the hypothalami. The increased GR predicted the decreased hypothalamic proopiomelanocortin expression and increased obesity that we observed in the 5-year-old adult males. The epigenetic and expression changes in GR were specific to the hypothalamus. Hippocampal GR mRNA and protein were decreased in UN offspring, whereas pituitary GR was altered in a sex-specific manner. In peripheral polymorphonuclear leukocytes there were no changes in GR methylation or protein, indicating that this epigenetic analysis did not predict changes in the brain. Overall, these results suggest that moderate changes in maternal nutrition, around the time of conception, signal life-long and tissue-specific epigenetic alterations in a key gene regulating energy balance in the hypothalamus.

Functionally relevant microsatellite markers from chickpea transcription factor genes for efficient genotyping applications and trait association mapping

We developed 1108 transcription factor gene-derived microsatellite (TFGMS) and 161 transcription factor functional domain-associated microsatellite (TFFDMS) markers from 707 TFs of chickpea. The robust amplification efficiency (96.5%) and high intra-specific polymorphic potential (34%) detected by markers suggest their immense utilities in efficient large-scale genotyping applications, including construction of both physical and functional transcript maps and understanding population structure. Candidate gene-based association analysis revealed strong genetic association of TFFDMS markers with three major seed and pod traits. Further, TFGMS markers in the 5' untranslated regions of TF genes showing differential expression during seed development had higher trait association potential. The significance of TFFDMS markers was demonstrated by correlating their allelic variation with amino acid sequence expansion/contraction in the functional domain and alteration of secondary protein structure encoded by genes. The seed weight-associated markers were validated through traditional bi-parental genetic mapping. The determination of gene-specific linkage disequilibrium (LD) patterns in desi and kabuli based on single nucleotide polymorphism-microsatellite marker haplotypes revealed extended LD decay, enhanced LD resolution and trait association potential of genes. The evolutionary history of a strong seed-size/weight-associated TF based on natural variation and haplotype sharing among desi, kabuli and wild unravelled useful information having implication for seed-size trait evolution during chickpea domestication.

Bovine TB in livestock and wildlife: what's in the genes?

Bovine tuberculosis (BTB) is a chronic, infectious disease found in domestic livestock and wildlife. It is caused predominantly by Mycobacterium bovis, which forms part of the Mycobacterium tuberculosis complex. BTB has serious implications for the movement of animals and animal products, biodiversity, and public health and is of significant economic concern. The existence of wildlife maintenance hosts makes it extremely difficult to eradicate BTB, even when established control strategies are in place, creating the need for alternative methods for controlling this disease. There are multiple factors that influence the outcome of infection by a pathogen, one of which is the host's genome. The identification of genetic variants involved in the susceptibility to BTB would supply a new selection of potential drug targets as well as the possibility for the breeding of animals with greater disease resistance. In this review, we collate the results of the BTB heritability and association studies performed in cattle and wildlife, discuss considerations and other methodologies (such as gene expression work) to be taken into account when performing genetic studies, and make some recommendations for future work in this area.

Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle

Background

Identification of single nucleotide polymorphisms (SNPs) for specific genes involved in reproduction might improve reliability of genomic estimates for these low-heritability traits. Semen from 550 Holstein bulls of high (≥ 1.7; n = 288) or low (≤ −2; n = 262) daughter pregnancy rate (DPR) was genotyped for 434 candidate SNPs using the Sequenom MassARRAY® system. Three types of SNPs were evaluated: SNPs previously reported to be associated with reproductive traits or physically close to genetic markers for reproduction, SNPs in genes that are well known to be involved in reproductive processes, and SNPs in genes that are differentially expressed between physiological conditions in a variety of tissues associated in reproductive function. Eleven reproduction and production traits were analyzed.

Results

A total of 40 SNPs were associated (P < 0.05) with DPR. Among these were genes involved in the endocrine system, cell signaling, immune function and inhibition of apoptosis. A total of 10 genes were regulated by estradiol. In addition, 22 SNPs were associated with heifer conception rate, 33 with cow conception rate, 36 with productive life, 34 with net merit, 23 with milk yield, 19 with fat yield, 13 with fat percent, 19 with protein yield, 22 with protein percent, and 13 with somatic cell score. The allele substitution effect for SNPs associated with heifer conception rate, cow conception rate, productive life and net merit were in the same direction as for DPR. Allele substitution effects for several SNPs associated with production traits were in the opposite direction as DPR. Nonetheless, there were 29 SNPs associated with DPR that were not negatively associated with production traits.

Conclusion

SNPs in a total of 40 genes associated with DPR were identified as well as SNPs for other traits. It might be feasible to include these SNPs into genomic tests of reproduction and other traits. The genes associated with DPR are likely to be important for understanding the physiology of reproduction. Given the large number of SNPs associated with DPR that were not negatively associated with production traits, it should be possible to select for DPR without compromising production.

Association studies indicate that protein disulfide isomerase is a risk factor in amyotrophic lateral sclerosis

Protein disulfide isomerase (PDI) plays an important role in the endoplasmic reticulum (ER) by facilitating the exchange of disulfide bonds and, together with other ER stress proteins, is induced in amyotrophic lateral sclerosis (ALS). However, genetic polymorphisms in the P4HB gene, which encodes PDI, have not been thoroughly investigated in ALS cases. In this study, we determined whether single-nucleotide polymorphisms (SNPs) in the P4HB gene were associated with familial ALS (FALS) and sporadic ALS (SALS). We report significant genotypic associations for two SNPs in P4HB with FALS, rs876016 (P=0.0198) and rs2070872 (P=0.0046), all values being FDR corrected. Significant allelic associations were also obtained for rs876016 with FALS (P=0.0155) and ALS (FALS and SALS) (P=0.0148). Four SNP haplotypes, which included two additional flanking SNPs, rs876017 and rs8324, were examined and rare haplotypes were found to be more common in ALS cases compared to controls. Seven haplotypes were significantly associated with FALS and one haplotype was significantly associated with SALS. One rare haplotype, which was present in controls, was overrepresented in a group of SOD1-positive FALS cases. Reduced survival was observed in FALS cases possessing at least one copy of the minor allele of rs2070872 (P=0.0059) and rs8324 (P=0.0167) and in individuals lacking the homozygous AAAC/AAAC diplotype (P=0.011). The results suggest that P4HB is a modifier gene in ALS susceptibility and may represent a potential therapeutic target for ALS.

Childhood adversity and DNA methylation of genes involved in the hypothalamus-pituitary-adrenal axis and immune system: whole-genome and candidate-gene associations

In recent years, translational research involving humans and animals has uncovered biological and physiological pathways that explain associations between early adverse circumstances and long-term mental and physical health outcomes. In this article, we summarize the human and animal literature demonstrating that epigenetic alterations in key biological systems, the hypothalamus-pituitary-adrenal axis and immune system, may underlie such disparities. We review evidence suggesting that changes in DNA methylation profiles of the genome may be responsible for the alterations in hypothalamus-pituitary-adrenal axis and immune system trajectories. Using some preliminary data, we demonstrate how explorations of genome-wide and candidate-gene DNA methylation profiles may inform hypotheses and guide future research efforts in these areas. We conclude our article by discussing the many important future directions, merging perspectives from developmental psychology, molecular genetics, neuroendocrinology, and immunology, that are essential for furthering our understanding of how early adverse circumstances may shape developmental trajectories, particularly in the areas of stress reactivity and physical or mental health.

Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population

Major histocompatibility complex (Mhc) genes are believed to play a key role in the genetic basis of disease control. Although numerous studies have sought links between Mhc and disease prevalence, many have ignored the ecological and epidemiological aspects of the host-parasite interaction. Consequently, interpreting associations between prevalence and Mhc has been difficult, whereas discriminating alleles for qualitative resistance, quantitative resistance and susceptibility remains challenging. Moreover, most studies to date have quantified associations between genotypes and disease status, overlooking the complex relationship between genotype and the properties of the Mhc molecule that interacts with parasites. Here, we address these problems and demonstrate avian malaria (Plasmodium) parasite species-specific associations with functional properties of Mhc molecules (Mhc supertypes) in a wild great tit (Parus major) population. We further show that correctly interpreting these associations depends crucially on understanding the spatial variation in risk of infection and the fitness effects of infection. We report that a single Mhc supertype confers qualitative resistance to Plasmodium relictum, whereas a different Mhc supertype confers quantitative resistance to Plasmodium circumflexum infections. Furthermore, we demonstrate common functional properties of Plasmodium-resistance alleles in passerine birds, suggesting this is a model system for parasite-Mhc associations in the wild.

Genetic predisposition to pass the standard SICCT test for bovine tuberculosis in British cattle

Bovine tuberculosis (bTB) imposes an important financial burden on the British cattle industry, yet despite intense efforts to control its spread, incidence is currently rising. Surveillance for bTB is based on a skin test that measures an immunological response to tuberculin. Cattle that fail the test are classified as "reactors" and slaughtered. Recent studies have identified genetic markers associated with the reaction of cattle to the tuberculin test. At marker INRA111 a relatively common '22' genotype occurs significantly more frequently in non-reactor cattle. Here we test the possibility that the putative protective '22' genotype does not confer resistance but instead causes cattle that carry it to react less strongly to the prescribed test, and hence avoid slaughter, potentially even though they are infected. We show that, after controlling for age and breed, '22' cattle react less strongly to the immunological challenge and may therefore be less likely to be classified as a reactor. These results highlight the potential discrepancy between infection and test status and imply that the effectiveness of the test-and-slaughter policy may be being compromised by selection for cattle that are genetically predisposed to react less strongly to tuberculin.