Haplotype tagging for the identification of common disease genes

Association between SLCO1B1 genetic polymorphisms and bleeding risk in patients treated with edoxaban

Since SLCO1B1 encodes the uptake transporter OATP1B1, which can influence the pharmacokinetic and pharmacodynamic profiles of edoxaban, polymorphisms in SLCO1B1 may affect the edoxaban response. This study aimed to investigate the association between SLCO1B1 gene polymorphisms and the bleeding risk in patients receiving edoxaban. We genotyped 10 single-nucleotide polymorphisms (SNPs) from the SLCO1B1 gene in patients receiving edoxaban. We also analyzed rs3842 of ABCB1 as a confounder. The odds ratio (OR) and adjusted OR (AOR) were calculated from univariate and multivariable analysis, respectively. The area under the receiver operating characteristic curve (AUROC) was constructed for the discrimination of the model. A total of 159 patients receiving edoxaban were analyzed. Overdose and rs4149056 showed significant association with bleeding complications by around 11- and 5.5-fold, respectively. Additionally, patients with the rs4149057 variant allele (C) had a 3.9-fold increased bleeding risk compared with wild-type homozygote carriers (TT), whereas rs2306283 variant homozygote (GG) carriers had a 0.27-fold reduced bleeding risk compared with wild-type allele (A) carriers. Patients with the variant-type homozygote (CC) of ABCB1 rs3842 had a higher bleeding risk than T allele carriers (AOR = 5.3 and 5.9). The final models for multivariable analyses were acceptable based on the AUROC values (> 0.70). These findings may help predict bleeding risk in patients taking edoxaban and help personalize treatment.

LmTag: functional-enrichment and imputation-aware tag SNP selection for population-specific genotyping arrays

Despite the rapid development of sequencing technology, single-nucleotide polymorphism (SNP) arrays are still the most cost-effective genotyping solutions for large-scale genomic research and applications. Recent years have witnessed the rapid development of numerous genotyping platforms of different sizes and designs, but population-specific platforms are still lacking, especially for those in developing countries. SNP arrays designed for these countries should be cost-effective (small size), yet incorporate key information needed to associate genotypes with traits. A key design principle for most current platforms is to improve genome-wide imputation so that more SNPs not included in the array (imputed SNPs) can be predicted. However, current tag SNP selection methods mostly focus on imputation accuracy and coverage, but not the functional content of the array. It is those functional SNPs that are most likely associated with traits. Here, we propose LmTag, a novel method for tag SNP selection that not only improves imputation performance but also prioritizes highly functional SNP markers. We apply LmTag on a wide range of populations using both public and in-house whole-genome sequencing databases. Our results show that LmTag improved both functional marker prioritization and genome-wide imputation accuracy compared to existing methods. This novel approach could contribute to the next generation genotyping arrays that provide excellent imputation capability as well as facilitate array-based functional genetic studies. Such arrays are particularly suitable for under-represented populations in developing countries or non-model species, where little genomics data are available while investment in genome sequencing or high-density SNP arrays is limited. $\textrm{LmTag}$ is available at: https://github.com/datngu/LmTag.

What have we learned from genome-wide association studies (GWAS) in Parkinson's disease?

After fifteen years of genome-wide association studies (GWAS) in Parkinson's disease (PD), what have we learned? Addressing this question will help catalogue the progress made towards elucidating disease mechanisms, improving the clinical utility of the identified loci, and envisioning how we can harness the strides to develop translational GWAS strategies. Here we review the advances of PD GWAS made to date while critically addressing the challenges and opportunities for next-generation GWAS. Thus, deciphering the missing heritability in underrepresented populations is currently at the reach of hand for a truly comprehensive understanding of the genetics of PD across the different ethnicities. Moreover, state-of-the-art GWAS designs hold a true potential for enhancing the clinical applicability of genetic findings, for instance, by improving disease prediction (PD risk and progression). Lastly, advanced PD GWAS findings, alone or in combination with clinical and environmental parameters, are expected to have the capacity for defining patient enriched cohorts stratified by genetic risk profiles and readily available for neuroprotective clinical trials. Overall, envisioning future strategies for advanced GWAS is currently timely and can be instrumental in providing novel genetic readouts essential for a true clinical translatability of PD genetic findings.

Multiple haplotype-based analyses provide genetic and evolutionary insights into tomato fruit weight and composition

Improving fruit quality traits such as metabolic composition remains a challenge for tomato breeders. To better understand the genetic architecture of these traits and decipher the demographic history of the loci controlling tomato quality traits, we applied an innovative approach using multiple haplotype-based analyses, aiming to test the potentials of haplotype based study in association and genomic prediction studies. We performed and compared haplotype vs SNP-based associations (hapQTL) with multi-locus mixed model (MLMM), focusing on tomato fruit weight and metabolite contents (i.e. sugars, organic acids and amino acids). Using a panel of 163 tomato accessions genotyped with 5995 SNPs, we detected a total of 784 haplotype blocks, with an average size of haplotype blocks ~58 kb. A total of 108 significant associations for 26 traits were detected thanks to Haplotype/SNP-based Bayes models. Haplotype-based Bayes model (97 associations) outperformed SNP-based Bayes model (50 associations) and MLMM (53 associations) in identifying marker-trait associations as well as in genomic prediction (especially for those traits with moderate to low heritability). To decipher the demographic history, we identified 24 positive selective sweeps using the integrated haplotype score (iHS). Most of the significant associations for tomato quality traits were located within selective sweeps (54.63% and 71.7% in hapQTL and MLMM models, respectively). Promising candidate genes were identified controlling tomato fruit weight and metabolite contents. We thus demonstrated the benefits of using haplotypes for evolutionary and genetic studies, providing novel insights into tomato quality improvement and breeding history.

Tools in Pharmacogenomics Biomarker Identification for Cancer Patients

The understanding of the biological differences which underlie the inter-individual variability in drug response improved the efficacy of cancer therapy in the era of precision medicine. In fact molecularly targeted drugs and immunotherapy represent a revolution in cancer treatment. The identification of genetic predictive and/or prognostic biomarkers linked to drug pharmacokinetics (PK) and pharmacodynamics (PD) is allowed by the development of high-throughput omics tools for detecting and understanding biological differences among individuals, in order to improve drug efficacy and minimize risk of toxicity. Personalized medicine in cancer treatment reduces costs of the healthcare system. Unfortunately, pharmacogenomics biomarkers discovery is influenced by complexity, need of high-quality evidence, and a validation process for regulatory purposes. This chapter is focused on the critic analysis of presently available pharmacogenomics tools for discovering or testing genetic polymorphic variants in drug metabolizing enzyme to be introduced in clinical practice for the prospective stratification of cancer patients.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

Pharmacogenomics Biomarker Discovery and Validation for Translation in Clinical Practice

Interindividual variability in drug efficacy and toxicity is a major challenge in clinical practice. Variations in drug pharmacokinetics (PKs) and pharmacodynamics (PDs) can be, in part, explained by polymorphic variants in genes encoding drug metabolizing enzymes and transporters (absorption, distribution, metabolism, and excretion) or in genes encoding drug receptors. Pharmacogenomics (PGx) has allowed the identification of predictive biomarkers of drug PKs and PDs and the current knowledge of genome-disease and genome-drug interactions offers the opportunity to optimize tailored drug therapy. High-throughput PGx genotyping, from targeted to more comprehensive strategies, allows the identification of PK/PD genotypes to be developed as clinical predictive biomarkers. However, a biomarker needs a robust process of validation followed by clinical-grade assay development and must comply to stringent regulatory guidelines. We here discuss the methodological challenges and the emerging technological tools in PGx biomarker discovery and validation, at the crossroad among molecular genetics, bioinformatics, and clinical medicine.

Effect of functional variant rs11466313 on breast cancer susceptibility and TGFB1 promoter activity

PURPOSE

This study aimed to investigate whether genetic polymorphisms in TGFB1 contribute to breast cancer (BC) susceptibility, and explore the mechanism of action.

METHODS

A total of 7 tagging SNPs (tSNPs) were genotyped in 1161 BC cases and 1337 age-matched controls among Chinese Han population. Bioinformatics analysis was used to predict functional SNP closely linked to tSNPs. Luciferase gene reporter assay was performed to determine the effect of genetic variants on promoter activity. DNA pull-down assay and mass spectrometry were used to identify the differentially binding proteins to genetic variants.

RESULTS

Genotyping analysis showed that rs1800469 (C>T) in the 5' regulatory region of TGFB1 was associated with reduced BC risk. Bioinformatics analysis predicted that rs11466313 (-2389_-2391 Del/AGG) in the 5' regulatory region of TGFB1, was closely linked to tSNP rs1800469 and could be functional. The genotyping of rs11466313 by PCR-SSCP showed that rs11466313 also conferred decreased BC risk. Luciferase assays demonstrated that rs11466313 minor allele reduced over ninefold of promoter activity compared with its major allele (p < 0.001). DNA pull-down assay and mass spectrometry revealed that rs11466313 minor allele lost the binding ability with FAM98B and HSP90B. Knocking down FAM98B but not HSP90B, the enhanced promoter activity driven by TGFB1 rs11466313 major allele was attenuated.

CONCLUSIONS

This study elucidates the impact of functional polymorphism rs11466313 in the regulatory region of TGFB1 on breast cancer susceptibility and gene expression, and could be helpful for future research to determine the value of this TGFB1 variant in the clinical setting.

Aggressive prostate cancer phenotype and genome-wide association studies: where are we now?

The majority of prostate cancer (PCa) is indolent, however, a percentage of patients are initially diagnosed with metastatic disease, for which there is a worse prognosis. There is a lack of biomarkers to identify men at greater risk for developing aggressive PCa. Genome-wide association studies (GWAS) scan the genome to search associations of SNPs with specific traits, like cancer. To date, eight GWAS have resulted in the reporting of 16 SNPs associated with aggressive PCa (p < 5.00 × 10^-2). Still, validation studies need to be conducted to confirm the obtained results as GWAS can generate false-positive results. Furthermore, post-GWAS studies provide a better understanding of the functional consequences.

Study of the Relationship between Polymorphisms in the IL-8 Gene Promoter Region and Coccidiosis Resistance Index in Jinghai Yellow Chickens

Interleukin 8 (IL-8) participates in the immune response and has the function of inducing neutrophils to release lysosomal enzymes and eliminate pathogens. This study was to investigate the effect of single nucleotide mutations in the IL-8 gene promoter region on the coccidiosis resistance index. In this study, 180 infected Eimeria tenella (E. tenella) Jinghai yellow chickens were used as experimental samples. DNA sequencing technology was used to detect single nucleotide polymorphisms (SNPs) in the IL-8 gene promoter region. The association between these SNPs and coccidiosis resistance indexes (including superoxide dismutase (SOD), malondialdehyde (MDA), glutathione peroxidase (GSH-PX), catalase (CAT), nitric oxide (NO), interleukin-1β (IL-1β), interleukin-2 (IL-2), interleukin-6 (IL-6), IL-8, and interferon-γ (IFN-γ)) were analyzed. Three SNPs (T-550C, G-398T, and T-360C) were detected. Significant associations were found between each genotype at the T-550C site with NO (p-value = 0.006) and IL-8 (p-value = 0.034) indexes. Significant associations were found between each genotype at the G-398T site with SOD (p-value = 0.042), CAT (p-value = 0.049), NO (p-value = 0.008), and IL-2 (p-value = 0.044) indexes. Significant associations were found between each genotype at the T-360C site with SOD (p-value = 0.007), NO (p-value = 0.046), IL-2 (p-value = 0.041), IL-8 (p-value = 0.039), and IFN-γ (p-value = 0.042) indexes. Haplotype analysis showed that multiple indexes of the H1H3 haplotype combination were significantly higher than other haplotype combinations. Therefore, mutation of the IL-8 gene promoter region has a significant regulatory effect on the coccidiosis resistance index, with a change in transcription factor binding potentially altering IL-8 gene expression, thereby further affecting the IL-8 level in plasma. However, the specific mechanism needs further study.

Association Analysis of Single Nucleotide Polymorphisms in the 5' Regulatory Region of the IL-6 Gene with Eimeria tenella Resistance in Jinghai Yellow Chickens

Interleukin 6 (IL-6) is an immunoregulatory cytokine involved in various inflammatory and immune responses. To investigate the effects of single nucleotide polymorphisms (SNPs) and haplotypes of IL-6 on resistance to Eimeria tenella (E. tenella), SNPs in the 5' regulatory region of IL-6 were detected with direct sequencing, and the effects of SNPs and haplotypes on resistance to E. tenella were analyzed by the least square model in Jinghai yellow chickens. Nineteen SNPs were identified in the 5' regulation region of IL-6, among which three SNPs were newly discovered. The SNP association analysis results showed that nine of the SNPs were significantly associated with E. tenella resistance indexes; the A-483G locus was significantly associated with the GSH-Px, IL-2, and IL-17 indexes (p < 0.05); the C-447G locus was significantly associated with the SOD, GSH-Px, IL-17, and IL-2 indexes (p < 0.05); and the G-357A locus had significant effects on the CAT and IL-16 indexes (p < 0.05). Haplotype analysis showed that H2H3 and H2H5 were favorable haplotype combinations with good coccidium resistance. Furthermore, we used qRT-PCR and observed that the expression of IL-6 in the infection group was higher than that in the control group in the liver, proventriculus, small intestine, thymus, kidney, and bursa of Fabricius and extremely significantly different than that in the cecum especially (p < 0.01). In summary, SNPs and haplotypes in the 5' regulatory region of IL-6 have important effects on E. tenella resistance, and the results will provide a reference for molecular marker selection of E. tenella resistance in Jinghai yellow chickens.

The prospect of discovering new biomarkers for ovarian cancer based on current knowledge of susceptibility loci and genetic variation (Review)

Ovarian cancer is the most lethal gynaecological malignancy. The cancer initially presents with non‑specific symptoms; thus, it is typically not discovered until the patient has reached the late, considerably more lethal, stages of the disease. Research focus is currently on finding novel biomarkers, especially for early detection and stratification of the disease. One promising approach has been to focus on mutations or variations in the genetic code that are associated with the risk of developing ovarian cancer. A certain heritable component is already known regarding genes such as BRCA1/2, TP53, MSH6, BRIP1 and RAD51C, yet these are estimated to only account for ~3.1% of the total risk. Recent advances in sequencing technologies have enabled the investigation of hundreds of thousands of genetic variants in genome‑wide association studies in tens of thousands of patients, which has led to the discovery of 108 (39 loci with P<5.0x10‑8) novel susceptibility loci for ovarian cancer, presented in this review. Using the published variants in a patient cohort screening, together with variants identified in our ongoing whole exome sequencing project, future aims are to ascertain whether certain of the novel variants could be used as biomarkers for early diagnosis and/or treatment decisions.

Sparse Convolutional Denoising Autoencoders for Genotype Imputation

Genotype imputation, where missing genotypes can be computationally imputed, is an essential tool in genomic analysis ranging from genome wide associations to phenotype prediction. Traditional genotype imputation methods are typically based on haplotype-clustering algorithms, hidden Markov models (HMMs), and statistical inference. Deep learning-based methods have been recently reported to suitably address the missing data problems in various fields. To explore the performance of deep learning for genotype imputation, in this study, we propose a deep model called a sparse convolutional denoising autoencoder (SCDA) to impute missing genotypes. We constructed the SCDA model using a convolutional layer that can extract various correlation or linkage patterns in the genotype data and applying a sparse weight matrix resulted from the L₁ regularization to handle high dimensional data. We comprehensively evaluated the performance of the SCDA model in different scenarios for genotype imputation on the yeast and human genotype data, respectively. Our results showed that SCDA has strong robustness and significantly outperforms popular reference-free imputation methods. This study thus points to another novel application of deep learning models for missing data imputation in genomic studies.

Natural variations of FT family genes in soybean varieties covering a wide range of maturity groups

BACKGROUND

Flowering time and maturity are among the most important adaptive traits in soybean (Glycine max (L.) Merill). Flowering Locus T (FT) family genes function as key flowering integrators, with flowering-promoting members GmFT2a/GmFT5a and flowering-inhibiting members GmFT4/GmFT1a antagonistically regulating vegetative and reproductive growth. However, to date, the relations between natural variations of FT family genes and the diversity of flowering time and maturity in soybean are not clear. Therefore, we conducted this study to discover natural variations in FT family genes in association with flowering time and maturity.

RESULTS

Ten FT family genes, GmFT1a, GmFT1b, GmFT2a, GmFT2b, GmFT3a, GmFT3b, GmFT4, GmFT5a, GmFT5b and GmFT6, were cloned and sequenced in the 127 varieties evenly covering all 14 known maturity groups (MG0000-MGX). They were diversified at the genome sequence polymorphism level. GmFT3b and GmFT5b might have experienced breeding selection in the process of soybean domestication and breeding. Haplotype analysis showed that a total of 17 haplotypes had correlative relationships with flowering time and maturity among the 10 FT genes, namely, 1a-H3, 1b-H1, 1b-H6, 1b-H7, 2a-H1, 2a-H3, 2a-H4, 2a-H9, 2b-H3, 2b-H4, 2b-H6, 2b-H7, 3b-H4, 5a-H1, 5a-H2, 5a-H4 and 5b-H1. Based on the association analysis, 38 polymorphic sites had a significant association with flowering time at the level of p < 0.01.

CONCLUSIONS

Some natural variations exist within the 10 FT family genes, which might be involved in soybean adaptation to different environments and have an influence on diverse flowering time and maturity. This study will facilitate the understanding of the roles of FTs in flowering and maturity.

Association between Fbxl5 gene polymorphisms and partial economic traits in Jinghai Yellow chickens

The Fbxl5 gene is a member of the F-BOX family and plays an important role in maintaining iron homeostasis in cells. In order to reveal the genetic effects of Fbxl5 gene polymorphisms on body weight (BW) traits and reproductive performance in chickens, Fbxl5 gene polymorphisms were detected in 363 Jinghai Yellow chickens by PCR-single-strand conformation polymorphism (SSCP) and DNA sequencing methods using three primers. With primer 1, three genotypes (BB, bb, Bb) were detected in the Jinghai Yellow chicken population and two mutations (g. 14257 T > C and g. 14262 T > C) were revealed by gene sequencing. With primer 2, two genotypes (EE, Ee) were detected in the same population and one mutation (g. 19018 G > A), and for primer 3, three genotypes (FF, ff, Ff) and one mutation (g. 19018 G > A) were detected. Four single nucleotide polymorphisms (SNPs) were used to estimate the frequency distributions of the eight haplotypes with PHASE 2.1 software. CTCG was the major haplotype with a frequency of 37.93 %, while the least frequent was TCTA with a frequency of 2.98 %. The BW of haplotype combination H1H8 was higher than that of the other haplotypes and was a dominant combination. In terms of reproductive performance, the age at the first egg of the haplotype combination H9H1 was later than in the other haplotypes, but the mean egg weight at 300 days was relatively optimal. The H1H2 haplotype produced the highest mean egg weight in 300 days, although the total number of eggs in 300 days was smaller in the H2H4 haplotype with the highest at first egg. Therefore, we can consider using the haplotype combination H1H2 for selection. The findings of this study expand the theoretical basis of the use of the Fbxl5 gene in the molecular breeding of poultry.

QTL Mapping and Validation of Adult Plant Resistance to Stripe Rust in Chinese Wheat Landrace Humai 15

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a devastating foliar disease that affects common wheat and barley throughout the world. The reasonable deployment of adult plant resistance (APR) wheat varieties is one of the best methods for controlling this disease. Wheat landraces are valuable resources for identifying the genes/QTLs responsible for disease resistance. Humai 15 is a Chinese spring wheat landrace and it has exhibited adequate levels of APR to the prevalent Pst races in field environments for many years. In this study, a population of 177 recombinant inbred lines (RILs) was derived from Humai 15 × Mingxian 169. After screening based on a 90K chip array using 45 RILs and Kompetitive Allelic Specific PCR marker genotyping for the population of RILs, a major effect QTL in Humai 15 was located on the centromere of chromosome 2B, where it accounted for up to 47.2% of the phenotypic variation. Two other minor QTL genes from Humai 15 were located on chromosome arms 3BS and 4BL. The Yr18 gene was identified on chromosome arm 7DS in Mingxian 169.

Tagging SNP-set selection with maximum information based on linkage disequilibrium structure in genome-wide association studies

Motivation

Effective tagging single-nucleotide polymorphism (SNP)-set selection is crucial to SNP-set analysis in genome-wide association studies (GWAS). Most of the existing tagging SNP-set selection methods cannot make full use of the information hidden in common or rare variants associated diseases. It is noticed that some SNPs have overlapping genetic information owing to linkage disequilibrium (LD) structure between SNPs. Therefore, when testing the association between SNPs and disease susceptibility, it is sufficient to elect the representative SNPs (called tag SNP-set or tagSNP-set) with maximum information.

Results

It is proposed a new tagSNP-set selection method based on LD information between SNPs, namely TagSNP-Set with Maximum Information. Compared with classical SNP-set analytical method, our method not only has higher power, but also can minimize the number of selected tagSNPs and maximize the information provided by selected tagSNPs with less genotyping cost and lower time complexity.

Contact

hesicheng12@163.com.

Supplementary information

Supplementary data are available at Bioinformatics online.

A rapidly photo-activatable light-up fluorescent nucleoside and its application in DNA base variation sensing

A new DNA building block (d(Tet)U) bearing a tetrazole and allyloxy group at N-phenyl ring linked through an aminopropynyl linker to the 5-position of 2'-deoxyuridine was synthesized. The modified DNA can be lit up via a photoinduced intramolecular tetrazole-alkene cycloaddition reaction, but quenched when the fully-matched double strand is formed. This conspicuous difference in fluorescence could open a door for DNA single nucleotide polymorphism (SNP) typing.

Identification and validation of regulatory SNPs that modulate transcription factor chromatin binding and gene expression in prostate cancer

Prostate cancer (PCa) is the second most common solid tumor for cancer related deaths in American men. Genome wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with the increased risk of PCa. Because most of the susceptibility SNPs are located in noncoding regions, little is known about their functional mechanisms. We hypothesize that functional SNPs reside in cell type-specific regulatory elements that mediate the binding of critical transcription factors (TFs), which in turn result in changes in target gene expression. Using PCa-specific functional genomics data, here we identify 38 regulatory candidate SNPs and their target genes in PCa. Through risk analysis by incorporating gene expression and clinical data, we identify 6 target genes (ZG16B, ANKRD5, RERE, FAM96B, NAALADL2 and GTPBP10) as significant predictors of PCa biochemical recurrence. In addition, 5 SNPs (rs2659051, rs10936845, rs9925556, rs6057110 and rs2742624) are selected for experimental validation using Chromatin immunoprecipitation (ChIP), dual-luciferase reporter assay in LNCaP cells, showing allele-specific enhancer activity. Furthermore, we delete the rs2742624-containing region using CRISPR/Cas9 genome editing and observe the drastic downregulation of its target gene UPK3A. Taken together, our results illustrate that this new methodology can be applied to identify regulatory SNPs and their target genes that likely impact PCa risk. We suggest that similar studies can be performed to characterize regulatory variants in other diseases.

Genetic predisposition to lung adenocarcinoma among never-smoking Chinese with different epidermal growth factor receptor mutation status

OBJECTIVES

The inconsistent findings from genetic association studies may be related to the heterogeneity in different molecular subtypes of lung cancer. This study evaluated the predisposing single-nucleotide polymorphisms (SNPs) in epidermal growth factor receptor (EGFR) mutant and EGFR wild-type lung adenocarcinoma separately among never-smokers.

MATERIALS AND METHODS

This was a two-stage case-control study. Never-smokers with pathologically confirmed lung adenocarcinoma and healthy controls were recruited in Hong Kong and Macau. Genomic DNA was extracted and genotyped by MassARRAY. In the discovery stage, 51 SNPs were investigated at the SNP, gene and pathway level among 103 EGFR mutant and 78 EGFR wild-type lung adenocarcinoma cases compared with matched controls. In the validation stage, SNPs that were identified with significant lung cancer risk were replicated in a separate cohort of 84 lung adenocarcinoma cases and compared with 103 Chinese Han, Beijing and 105 Chinese Han, Southern public controls from the 1000 genome database.

RESULTS AND CONCLUSION

The genetic association of IL-6 rs2069840 with EGFR mutant lung adenocarcinoma was ascertained. In the discovery stage, haplotype GGG in three SNPs (rs2069840, rs2069852, rs2066992) of IL-6, synergetic effects of IL-6 rs2069840 and environmental tobacco smoke in the workplace were found to be related to EGFR mutant lung adenocarcinoma. ERCC2 rs238406 showed a marginally significant association with EGFR mutant lung adenocarcinoma in the validation stage (P=0.096). ERCC2 rs50871 and ATM rs611646 showed significant association with EGFR wild-type lung adenocarcinoma in the discovery stage. In conclusion, IL-6 rs2069840 conferred susceptibility to EGFR mutant lung adenocarcinoma in a Hong Kong and Macau never-smoking Chinese population.

Association between SLC19A1 Gene Polymorphism and High Dose Methotrexate Toxicity in Childhood Acute Lymphoblastic Leukaemia and Non Hodgkin Malignant Lymphoma: Introducing a Haplotype based Approach

Background

We investigated the clinical relevance of SLC 19A1 genetic variability for high dose methotrexate (HD-MTX) related toxicities in children and adolescents with acute lymphoblastic leukaemia (ALL) and non Hodgkin malignant lymphoma (NHML).

Patients and methods

Eighty-eight children and adolescents with ALL/NHML were investigated for the influence of SLC 19A1 single nucleotide polymorphisms (SNPs) and haplotypes on HD-MTX induced toxicities.

Results

Patients with rs2838958 TT genotype had higher probability for mucositis development as compared to carriers of at least one rs2838958 C allele (OR 0.226 (0.071–0.725), p < 0.009). Haplotype TGTTCCG (H4) statistically significantly reduced the risk for the occurrence of adverse events during treatment with HD-MTX (OR 0.143 (0.023–0.852), p = 0.030).

Conclusions

SLC 19A1 SNP and haplotype analysis could provide additional information in a personalized HD-MTX therapy for children with ALL/NHML in order to achieve better treatment outcome. However further studies are needed to validate the results.

Haplotypes for Type, Degree, and Rate of Marbling in Cattle Are Syntenic with Human Muscular Dystrophy

Traditional analyses of a QTL on Bota 19 implicate a surfeit of candidates, but each is of marginal significance in explaining the deposition of healthy, low melting temperature fat within marbled muscle of Wagyu cattle. As an alternative approach, we have used genomic, multigenerational segregation to identify 14 conserved, ancestral 20 Mb haplotypes. These determine the degree and rate of marbling in Wagyu and other breeds of cattle. The melting temperature of intramuscular fat is highly heritable and traceable by haplotyping. Fortunately, for the production of healthy beef, some of these haplotypes are sufficiently penetrant to be expressed in heterozygous crossbreds, thereby allowing selection of sires which will improve the healthiness of beef produced under even harsh climatic conditions. The region of Bota 19 is syntenic to a region of Hosa 17 known to be important in muscle metabolism and in determining susceptibility to a form of human muscular dystrophy.

The Molecular Revolution in Cutaneous Biology: The Era of Genome-Wide Association Studies and Statistical, Big Data, and Computational Topics

The investigation of biological systems involving all organs of the body including the skin is in era of big data. This requires heavy-duty computational tools, and novel statistical methods. Microarrays have allowed the interrogation of thousands of common genetic markers in thousands of individuals from the same population (termed genome wide association studies or GWAS) to reveal common variation associated with disease or phenotype. These markers are usually single nucleotide polymorphisms (SNPs) that are relatively common in the population. In the case of dermatological diseases such as alopecia areata, vitiligo, psoriasis and atopic dermatitis, common variants have been identified that are associated with disease, and these provide insights into biological pathways and reveal possible novel drug targets. Other skin phenotypes such as acne, color and skin cancers are also being investigated with GWAS. Analyses of such large GWAS datasets require a consideration of a number of statistical issues including the testing of multiple markers, population substructure, and ultimately a requirement for replication. There are also issues regarding the missing heritability of disease that cannot be entirely explained with current GWAS approaches. Next generation sequencing technologies such as exome and genome sequencing of similar patient cohorts will reveal additional variants contributing to disease susceptibility. However, the data generated with these approaches will be orders of magnitude greater than that those generated with arrays, with concomitant challenges in the identification of disease causing variants.

New multilocus linkage disequilibrium measure for tag SNP selection

Numerous approaches have been proposed for selecting an optimal tag single-nucleotide polymorphism (SNP) set. Most of these approaches are based on linkage disequilibrium (LD). Classical LD measures, such as D′ and r2, are frequently used to quantify the relationship between two marker (pairwise) linkage disequilibria. Despite of their successful use in many applications, these measures cannot be used to measure the LD between multiple-marker. These LD measures need information about the frequencies of alleles collected from haplotype dataset. In this study, a cluster algorithm is proposed to cluster SNPs according to multilocus LD measure which is based on information theory. After that, tag SNPs are selected in each cluster optimized by the number of tag SNPs, prediction accuracy and so on. The experimental results show that this new LD measure can be directly applied to genotype dataset collected from the HapMap project, so that it saves the cost of haplotyping. More importantly, the proposed method significantly improves the efficiency and prediction accuracy of tag SNP selection.

CTLA-4 polymorphisms associate with breast cancer susceptibility in Asians: a meta-analysis

Previous studies have investigated the association between cytotoxic T-lymphocyte antigen-4 (CTLA-4) polymorphisms and breast cancer susceptibility, but the results remained inconsistent. Therefore, we evaluated the relationship between four common CTLA-4 polymorphisms and breast cancer risk by a meta-analysis, aiming to derive a comprehensive and precise conclusion. We searched EMBASE, Pubmed, Web of Science, CNKI, and Wanfang databases until July 18th, 2016. Finally, ten eligible studies involving 4,544 breast cancer patients and 4,515 cancer-free controls were included; all these studies were from Asia. Odds ratio (OR) and 95% confidence interval (CI) were used to evaluate the breast cancer risk in five genetic models. The results indicated that the CTLA-4 +49A>G (rs231775) polymorphism had a significant association with decreased breast cancer risk in allelic, homozygous, dominant and recessive models. Also, the +6230G>A (rs3087243) polymorphism reduced breast cancer risk especially in the Chinese population under homozygous and recessive models. In contrast, the -1661A>G (rs4553808) polymorphism increased breast cancer risk in allelic, heterozygous and dominant models, whereas -1722 T>C (rs733618) did not relate to breast cancer risk. In conclusion, CTLA-4 polymorphisms significantly associate with breast cancer susceptibility in Asian populations, and different gene loci may have different effects on breast cancer development. Further large-scale studies including multi-racial populations are required to confirm our findings.

Discovering Genome-Wide Tag SNPs Based on the Mutual Information of the Variants

Exploring linkage disequilibrium (LD) patterns among the single nucleotide polymorphism (SNP) sites can improve the accuracy and cost-effectiveness of genomic association studies, whereby representative (tag) SNPs are identified to sufficiently represent the genomic diversity in populations. There has been considerable amount of effort in developing efficient algorithms to select tag SNPs from the growing large-scale data sets. Methods using the classical pairwise-LD and multi-locus LD measures have been proposed that aim to reduce the computational complexity and to increase the accuracy, respectively. The present work solves the tag SNP selection problem by efficiently balancing the computational complexity and accuracy, and improves the coverage in genomic diversity in a cost-effective manner. The employed algorithm makes use of mutual information to explore the multi-locus association between SNPs and can handle different data types and conditions. Experiments with benchmark HapMap data sets show comparable or better performance against the state-of-the-art algorithms. In particular, as a novel application, the genome-wide SNP tagging is performed in the 1000 Genomes Project data sets, and produced a well-annotated database of tagging variants that capture the common genotype diversity in 2,504 samples from 26 human populations. Compared to conventional methods, the algorithm requires as input only the genotype (or haplotype) sequences, can scale up to genome-wide analyses, and produces accurate solutions with more information-rich output, providing an improved platform for researchers towards the subsequent association studies.

Association of SNP Haplotypes of HKT Family Genes with Salt Tolerance in Indian Wild Rice Germplasm

BACKGROUND

Rice is one of the most important crops for global food security but its productivity is adversely affected by salt stress prevalent in about 30 % of the cultivated land. For developing salt-tolerant rice varieties through conventional breeding or biotechnological interventions, there is an urgent need to identify natural allelic variants that may confer salt tolerance. Here, 299 wild rice accessions collected from different agro-climatic regions of India were evaluated during growth under salt stress. Of these 95 representative accessions were sequenced for members of HKT ion transporter family genes by employing Ion Torrent PGM sequencing platform.

RESULTS

Haplotype analysis revealed haplotypes H5 and H1 of HKT1;5 and HKT2;3, respectively associated with high salinity tolerance. This is the first study of allele mining of eight members of HKT gene family from Indian wild rice reporting a salt tolerant allele of HKT2;3. HKT1;5 also showed a salt tolerant allele from wild rice. Phylogenetic analysis based on the nucleotide sequences showed different grouping of the HKT family genes as compared to the prevailing protein sequence based classification.

CONCLUSIONS

The salt tolerant alleles of the HKT genes from wild rice may be introgressed into modern high yielding cultivars to widen the existing gene pool and enhance rice production in the salt affected areas.

Association of SNP Haplotypes of HKT Family Genes with Salt Tolerance in Indian Wild Rice Germplasm

BACKGROUND

Rice is one of the most important crops for global food security but its productivity is adversely affected by salt stress prevalent in about 30 % of the cultivated land. For developing salt-tolerant rice varieties through conventional breeding or biotechnological interventions, there is an urgent need to identify natural allelic variants that may confer salt tolerance. Here, 299 wild rice accessions collected from different agro-climatic regions of India were evaluated during growth under salt stress. Of these 95 representative accessions were sequenced for members of HKT ion transporter family genes by employing Ion Torrent PGM sequencing platform.

RESULTS

Haplotype analysis revealed haplotypes H5 and H1 of HKT1;5 and HKT2;3, respectively associated with high salinity tolerance. This is the first study of allele mining of eight members of HKT gene family from Indian wild rice reporting a salt tolerant allele of HKT2;3. HKT1;5 also showed a salt tolerant allele from wild rice. Phylogenetic analysis based on the nucleotide sequences showed different grouping of the HKT family genes as compared to the prevailing protein sequence based classification.

CONCLUSIONS

The salt tolerant alleles of the HKT genes from wild rice may be introgressed into modern high yielding cultivars to widen the existing gene pool and enhance rice production in the salt affected areas.

Associations between CD36 gene polymorphisms, fat tolerance and oral fat preference in a young-adult population

BACKGROUND/OBJECTIVES

CD36 is known to be an orosensory receptor for dietary long-chain fatty acids, as well as being involved in the chemosensory mechanisms within the human gut. Recent data have demonstrated an association between CD36 single-nucleotide polymorphisms (SNPs) and lipid consumption behaviours in humans. This study aimed to test for associations between CD36 SNPs and response to a high-fat meal in a young healthy Australian cohort. Secondary associations were tested between CD36 gene variants and fasting lipid parameters, body composition, cardiovascular disease (CVD) risk factors and measures of oral fat preference.

SUBJECTS/METHODS

Two SNPs (rs1527479 and rs1984112) were assessed for associations with response to a 75 g saturated fat oral fat tolerance test (OFTT), whole-body substrate oxidation, fasting plasma lipids, CVD risk factors and self-reported habitual diet questionnaires. Genotyping was performed using real-time polymerase chain reaction.

RESULTS

Cross-sectional data were collected on 56 individuals (28 m, 28 f; 24.9±3.3 years), with 42 completing participation in a high-fat OFTT. No genotypic associations were evident in anthropometric data or self-reported fat preference measures. AA SNP carriers at rs1984112 exhibited significantly elevated fasting triglyceride when compared with non-carriers (P=0.024). This group also tended to have an elevated response to a high-fat meal (P=0.078).

CONCLUSIONS

Although these data show the potential pleiotropic influence of CD36 SNP rs1984112 on lipoprotein accumulation in a young healthy cohort, further assessment in a larger cohort is warranted.

Association between MMP-3 and MMP-9 polymorphisms and coronary artery disease

Matrix metalloproteinase (MMP)-3 and MMP-9 polymorphisms are characterized by plaque stability in coronary arteries. The aim of the current study was to investigate the 5A/6A polymorphism in the MMP-3 gene and C/T polymorphism in the MMP-9 gene in patients with coronary artery disease (CAD). The study population consisted of 400 patients who underwent coronary angiography. There were two groups consisting of 200 consecutive patients with CAD, presenting with stable angina pectoris, and 200 consecutive patients exhibiting normal coronary arteries. Two single nucleotide polymorphisms in the MMP gene, MMP-3 and MMP-9, were detected using a polymerase chain reaction-restriction fragment length polymorphism assay. Mean age, gender distribution, smoking status, presence of diabetes mellitus and hypercholesterolemia were identified to be similar between the groups. One hundred and twenty seven (63.5%) patients had hypertension in the CAD group, whereas only 55 (27.5%) patients had hypertension in the control group (P<0.001). No significant difference in frequency of alleles and genotypes of MMP-9 C→T between the CAD and control groups was identified. The 5A allele frequency of MMP-3 in the CAD group was significantly higher when compared with the control group (P<0.001; odds ratio=2.18). The genotype frequency of MMP-3 5A/5A in the CAD group was significantly higher when compared with the controls (P=0.005). When compared with the homozygous wild-type (6A/6A) genotype of the MMP-3 gene, the cumulative frequency of heterozygote and homozygote genotypes of the MMP-3 gene was significantly higher in the CAD compared with the control group (P<0.001). Thus, the present study demonstrated that the 5A/5A and 6A/5A+5A/5A genotypes of the MMP-3 gene were associated with an increased risk of CAD.

Quantitative haplotyping of PCR products by nonsynchronous pyrosequencing with di-base addition

Molecular haplotyping is becoming increasingly important for studying the disease association of a specific allele because of its ability of providing more information than any single nucleotide polymorphism (SNP). Computational analysis and experimental techniques are usually performed for haplotypic determination. However, established methods are not suitable for analyzing haplotypes of massive natural DNA samples. Here we present a simple molecular approach to analyze haplotypes of conventional polymerase chain reaction (PCR) products quantitatively in a single sequencing run. In this approach, specific types and proportions of haplotypes in both individual and pooled samples could be determined by solving equations constructed from nonsynchronous pyrosequencing with di-base addition. Two SNPs (rs11176013 and rs11564148) in the gene for leucine-rich repeat kinase 2 (LRRK2) related to Parkinson's disease were selected as experimental sites. A series of DNA samples, including these two heterozygous loci, were investigated. This approach could accurately identify multiple DNA samples indicating that the approach is likely to be applied for haplotyping of unrestricted conventional PCR products from natural samples, and be especially applicable for analyzing short sequences in clinical diagnosis. Graphical Abstract One DNA sample consisting of 4 different DNA templates with different proportion are sequenced by nonsynchronous pyrosequencing with di-base addition. The number of incorporated nucleotides produced by a single sequencing reaction equals to the total of incorporated nucleotides. Four independent equations are constructed from the pyrograms of nonsynchronous pyrosequencing data. Molecular haplotypes of two adjacent SNPs can be quantitatively identified by solving these equations.

Lack of association between cytotoxic T-lymphocyte antigen-4 gene polymorphisms and lymphoid malignancy risk: evidence from a meta-analysis

Cytotoxic T-lymphocyte antigen-4 (CTLA-4) polymorphisms have been associated with susceptibility to lymphoid malignancies. However, results from the published single studies are inconsistent. Therefore, the present meta-analysis was conducted to get a more accurate estimation of the relationship between CTLA-4 gene polymorphisms and the lymphoid malignancy risk. We identified nine independent studies accounting for 3090 subjects up to January 30, 2016. Summary odds ratios (OR) and 95 % confidence intervals (CI) were used to evaluate the risk of lymphoid malignancies. Overall, no significant association was found between +49A/G (rs231775), -318C/T (rs5742909), and +6230A/G (rs3087243) CTLA-4 gene polymorphisms and lymphoid malignancies. Furthermore, ethnicity (Asian and Caucasian) and histopathology subgroup analyses (non-Hodgkin's lymphoma) also failed to detect an association between the studied polymorphisms and lymphoid malignancy risk. Our study shows that common CTLA-4 gene polymorphisms may not contribute to lymphoid malignancy susceptibility based on the current evidence.

Variations at regulatory regions of the milk protein genes are associated with milk traits and coagulation properties in the Sarda sheep

Regulatory variation at the ovine casein genes could have important effects on the composition and coagulation properties of milk. Herewith, we have partially resequenced the promoters and the 3'-UTR of the four casein genes in 25 Sarda sheep. Alignment of these sequences allowed us to identify a total of 29 SNPs. This level of polymorphism (one SNP every 250 bp) is remarkably high if compared with SNP densities estimated in human genic regions (approximately one SNP per bp). The 29 SNPs identified in our resequencing experiment, plus three previously reported SNPs mapping to the lactalbumin, alpha (LALBA) and β-lactoglobulin (BLG, also known as progestagen-associated endometrial protein, PAEP) genes, were genotyped with a multiplex TaqMan Open Array Real-Time PCR assay in 760 Sarda sheep with records for milk composition and coagulation properties. Association analysis revealed the existence of significant associations of CSN1S2 and CSN3 genotypes with milk protein and casein contents. Moreover, genotypes at CSN1S1 were significantly associated with rennet coagulation time, curd firming time and curd firmness, whereas CSN2 was associated with curd firming time. These results suggest that SNPs mapping to the promoters and 3'-UTRs of ovine casein genes may exert regulatory effects on gene expression and that they could be used for improving sheep milk quality and technological traits at the population level through marker assisted selection.

A bioinformatics strategy for detecting the complexity of Chronic Obstructive Pulmonary Disease in Northern Chinese Han Population

Chronic Obstructive Pulmonary Disease (COPD) is a complex human disease which is driven not only by genetic factors, but also by various environmental variables, such as gender, age and smoking. Therefore, there is a demand for investigating the complexity among various risk factors involved in COPD. In this study, 44 tagging SNPs from EPHX1, GSTP1, SERPINE2 and TGFB1 were selected and genotyped in 310 COPD cases and 203 controls, all of which belong to the Han from North China. We integrated functional prediction algorithms of nonsynonymous SNPs (nsSNPs) into Bayesian network to explore the complex regulatory relationships among disease traits and various risk factors. The results showed that three basic variables (age, sex and smoking) were risk factors of COPD-related trait and phenotype. Besides these environmental risk factors, deleterious nsSNPs were found to perform better than those of significant synonymous SNPs when used as variables to make risk prediction of disease outcome. This study provides further evidences for detecting the complexity of COPD in Northern Chinese Han Population.

CTLA-4+49 A/G polymorphism and antiglutamic acid decarboxylase antibody-associated encephalopathy in Taiwanese children

BACKGROUND

Anti-glutamic acid decarboxylase antibodies are associated with encephalopathy, an autoimmune central nervous system inflammatory disease. The cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4)+49 A/G polymorphism has been shown to confer genetic susceptibility to positive anti-glutamic acid decarboxylase antibodies in patients with type 1 diabetes mellitus in Japan. We aimed to investigate the association of the CTLA-4+49 A/G (rs231775) polymorphism in Taiwanese children with anti-glutamic acid decarboxylase antibody-associated encephalopathy.

METHODS

This was a case-control study from July 2011 to June 2012 performed at Chang Gung Children's Hospital in Taiwan. Genotyping of the CTLA-4+49 A/G polymorphism was performed by polymerase chain reaction-restriction fragment length polymorphism.

RESULTS

Seventeen patients with anti-glutamic acid decarboxylase antibody-associated encephalopathy and 97 controls were enrolled. The genotype, allele and carrier frequencies of the CTLA-4+49 A/G polymorphism were equally distributed in the patients and controls, with no significant differences between the two groups. In addition, we found a positive trend between the level of anti-glutamic acid decarboxylase antibodies and the G allele of the CTLA-4+49 A/G polymorphism, although this trend was not statistically significant.

CONCLUSIONS

Our results suggest that the CTLA-4+49 A/G (rs231775) polymorphism does not confer an increased susceptibility to anti-glutamic acid decarboxylase antibody-associated encephalopathy in Taiwanese children.

MyD88 Polymorphisms and Association with Susceptibility to Salmonella Pullorum

Myeloid differentiation primary response gene 88 (MYD88), a universal adapter protein, plays an important role in activating the nuclear factor-κB (NF-κB) and regulating the expression of proinflammatory genes like tumor necrosis factor (TNF) and interleukin-1 (IL-1), which were highly involved in Salmonella Pullorum infection. To detect the relationship between polymorphisms of the MyD88 gene and Salmonella Pullorum disease, we screened the coding region (CDS) of the MYD88 gene by DNA pool construction and sequencing based on case-control study. Eight single nucleotide polymorphisms (SNPs) in the sequenced fragment (5 exons), 7 known loci and one novel mutation named G4810372T (SNP8), were found in the fifth exon. In addition, we found 7 nonsynonymous substitutions. The allele frequency of only one SNP, g.4810191C > T (SNP1), was significantly different (P < 0.05) between case and control groups. The genotype frequencies of SNP1 (g.4810191C > T) and SNP3 (g.4810257G > T) were of significant difference between the case and the control groups (P < 0.05). Collectively, SNPs of the MyD88 gene were significantly associated with susceptibility to Salmonella Pullorum infection, which can be used as a disease-resistant marker in chicken. These results provided a theoretical basis for future research on chicken breeding by marker-assisted selection.

Whole genome sequences are required to fully resolve the linkage disequilibrium structure of human populations

BACKGROUND

An understanding of linkage disequilibrium (LD) structures in the human genome underpins much of medical genetics and provides a basis for disease gene mapping and investigating biological mechanisms such as recombination and selection. Whole genome sequencing (WGS) provides the opportunity to determine LD structures at maximal resolution.

RESULTS

We compare LD maps constructed from WGS data with LD maps produced from the array-based HapMap dataset, for representative European and African populations. WGS provides up to 5.7-fold greater SNP density than array-based data and achieves much greater resolution of LD structure, allowing for identification of up to 2.8-fold more regions of intense recombination. The absence of ascertainment bias in variant genotyping improves the population representativeness of the WGS maps, and highlights the extent of uncaptured variation using array genotyping methodologies. The complete capture of LD patterns using WGS allows for higher genome-wide association study (GWAS) power compared to array-based GWAS, with WGS also allowing for the analysis of rare variation. The impact of marker ascertainment issues in arrays has been greatest for Sub-Saharan African populations where larger sample sizes and substantially higher marker densities are required to fully resolve the LD structure.

CONCLUSIONS

WGS provides the best possible resource for LD mapping due to the maximal marker density and lack of ascertainment bias. WGS LD maps provide a rich resource for medical and population genetics studies. The increasing availability of WGS data for large populations will allow for improved research utilising LD, such as GWAS and recombination biology studies.

Fine mapping of Msv1, a major QTL for resistance to Maize Streak Virus leads to development of production markers for breeding pipelines

Msv1 , the major QTL for MSV resistance was delimited to an interval of 0.87 cM on chromosome 1 at 87 Mb and production markers with high prediction accuracy were developed. Maize streak virus (MSV) disease is a devastating disease in the Sub-Saharan Africa (SSA), which causes significant yield loss in maize. Resistance to MSV has previously been mapped to a major QTL (Msv1) on chromosome 1 that is germplasm and environment independent and to several minor loci elsewhere in the genome. In this study, Msv1 was fine-mapped through QTL isogenic recombinant strategy using a large F 2 population of CML206 × CML312 to an interval of 0.87 cM on chromosome 1. Genome-wide association study was conducted in the DTMA (Drought Tolerant Maize for Africa)-Association mapping panel with 278 tropical/sub-tropical breeding lines from CIMMYT using the high-density genotyping-by-sequencing (GBS) markers. This study identified 19 SNPs in the region between 82 and 93 Mb on chromosome 1(B73 RefGen_V2) at a P < 1.00E-04, which coincided with the fine-mapped region of Msv1. Haplotype trend regression identified a haplotype block significantly associated with response to MSV. Three SNPs in this haplotype block at 87 Mb on chromosome 1 had an accuracy of 0.94 in predicting the disease reaction in a collection of breeding lines with known responses to MSV infection. In two biparental populations, selection for resistant Msv1 haplotype demonstrated a reduction of 1.03-1.39 units on a rating scale of 1-5, compared to the susceptible haplotype. High-throughput KASP assays have been developed for these three SNPs to enable routine marker screening in the breeding pipeline for MSV resistance.

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.

CTLA4 polymorphisms and de novo malignancy risk after renal transplantation in Chinese recipients

Genetic polymorphisms in cytotoxic T lymphocyte-associated antigen 4 (CTLA4) play an influential role in graft rejection and the long-term clinical outcome of organ transplantation. We investigated the association of five CTLA4 single-nucleotide polymorphisms (SNPs) (rs733618 C/T, rs4553808 A/G, rs5742909 C/T, rs231775 A/G, and rs3087243 G/A) with de novo malignancy in 1463 Chinese renal transplantation (RT) recipients who underwent a 192-month follow-up. Multivariate analyses revealed that recipient rs231775 genotype is significantly associated with tumorigenesis (P = 0.012). Multiplicative interaction between rs231775 AA and possible risk factors of malignancy revealed two significant results: rs231775 AA × primary diseases and rs231775 AA × number of HLA-mismatch. The frequency of haplotype TACAG was significantly higher in the tumor group (17.07%) than that in the nontumor group (1.53%). In addition, aristolochic acid nephropathy (P = 0.003) and the time of discovery of tumor (P = 0.000) also were independently associated with tumorigenesis. Our data show that the CTLA4 genotype rs231775 AA may be one of risk factors for the development of malignancy and haplotype TACAG was susceptible haplotype in Chinese kidney transplant recipients.

Cloud computing-based TagSNP selection algorithm for human genome data

Single nucleotide polymorphisms (SNPs) play a fundamental role in human genetic variation and are used in medical diagnostics, phylogeny construction, and drug design. They provide the highest-resolution genetic fingerprint for identifying disease associations and human features. Haplotypes are regions of linked genetic variants that are closely spaced on the genome and tend to be inherited together. Genetics research has revealed SNPs within certain haplotype blocks that introduce few distinct common haplotypes into most of the population. Haplotype block structures are used in association-based methods to map disease genes. In this paper, we propose an efficient algorithm for identifying haplotype blocks in the genome. In chromosomal haplotype data retrieved from the HapMap project website, the proposed algorithm identified longer haplotype blocks than an existing algorithm. To enhance its performance, we extended the proposed algorithm into a parallel algorithm that copies data in parallel via the Hadoop MapReduce framework. The proposed MapReduce-paralleled combinatorial algorithm performed well on real-world data obtained from the HapMap dataset; the improvement in computational efficiency was proportional to the number of processors used.

Phenotypic effects of concomitant insensitive acetylcholinesterase (ace-1(R)) and knockdown resistance (kdr(R)) in Anopheles gambiae: a hindrance for insecticide resistance management for malaria vector control

Background

Malaria is endemic in sub-Saharan Africa with considerable burden for human health. Major insecticide resistance mechanisms such as kdr^R and ace-1^Ralleles constitute a hindrance to malaria vector control programs. Anopheles gambiae bearing both kdr and ace-1 resistant alleles are increasingly recorded in wild populations. In order to maintain the efficacy of vector control strategies, the characterization of concomitant kdr and ace-1 resistance, and their pleiotropic effects on malaria vector phenotype on insecticide efficacy are important.

Methods

Larval and adult bioassays were performed with different insecticide classes used in public health following WHO standard guidelines on four laboratory Anopheles gambiae strains, sharing the same genetic background but harboring distinct resistance status: KISUMU with no resistance allele; ACERKIS with ace-1^R allele; KISKDR with kdr^R allele and ACERKDRKIS with both resistance alleles’ ace-1^R and kdr^R.

Results

Larval bioassays indicate that the homozygote resistant strain harboring both alleles (ACERKDRKIS) displayed slightly but significantly higher resistance level to various insecticides like carbamates (bendiocarb, p < 0.001; propoxur, p = 0.02) and organophosphates (chlorpyriphos-methyl, p = 0.002; fenitrothion, p < 0.001) when compared to ACERKIS strain. However, no differences were recorded between ACERKDRKIS and KISKDR resistance level against permethrin (Pyrethroid, p = 0.7) and DDT (Organochlorine, p = 0.24). For adult bioassays, the percentages of mosquitoes knocked down were significantly lower for ACERKDRKIS than for KISKDR with permethrin (p = 0.003) but not with deltamethrin. The percentage of mortality from adult bioassays was similar between ACERKDRKIS and ACERKIS with carbamates and organophosphates, or between ACERKDRKIS and KISKDR with pyrethroid and DDT. Concerning acetylcholinesterase enzyme, ACERKDRKIS strain showed similarAChE1 activity than that of ACERKIS.

Conclusion

The presence of both kdr^R and ace-1^R alleles seems to increase the resistance levels to both carbamate and organophosphate insecticides and at operational level, may represent an important threat to malaria vector control programs in West Africa.

Influence of kynurenine 3-monooxygenase (KMO) gene polymorphism on cognitive function in schizophrenia

BACKGROUND

Cognitive deficits compromise quality of life and productivity for individuals with schizophrenia and have no effective treatments. Preclinical data point to the kynurenine pathway of tryptophan metabolism as a potential target for pro-cognitive drug development. We have previously demonstrated association of a kynurenine 3-monooxygenase (KMO) gene variant with reduced KMO gene expression in postmortem schizophrenia cortex, and neurocognitive endophenotypic deficits in a clinical sample. KMO encodes kynurenine 3-monooxygenase (KMO), the rate-limiting microglial enzyme of cortical kynurenine metabolism. Aberration of the KMO gene might be the proximal cause of impaired cortical kynurenine metabolism observed in schizophrenia. However, the relationship between KMO variation and cognitive function in schizophrenia is unknown. This study examined the effects of the KMO rs2275163C>T C (risk) allele on cognitive function in schizophrenia.

METHODS

We examined the association of KMO polymorphisms with general neuropsychological performance and P50 gating in a sample of 150 schizophrenia and 95 healthy controls.

RESULTS

Consistent with our original report, the KMO rs2275163C>T C (risk) allele was associated with deficits in general neuropsychological performance, and this effect was more marked in schizophrenia compared with controls. Additionally, the C (Arg452) allele of the missense rs1053230C>T variant (KMO Arg452Cys) showed a trend effect on cognitive function. Neither variant affected P50 gating.

CONCLUSIONS

These data suggest that KMO variation influences a range of cognitive domains known to predict functional outcome. Extensive molecular characterization of this gene would elucidate its role in cognitive function with implications for vertical integration with basic discovery.

Association analysis of the SLC22A11 (organic anion transporter 4) and SLC22A12 (urate transporter 1) urate transporter locus with gout in New Zealand case-control sample sets reveals multiple ancestral-specific effects

INTRODUCTION

There is inconsistent association between urate transporters SLC22A11 (organic anion transporter 4 (OAT4)) and SLC22A12 (urate transporter 1 (URAT1)) and risk of gout. New Zealand (NZ) Māori and Pacific Island people have higher serum urate and more severe gout than European people. The aim of this study was to test genetic variation across the SLC22A11/SLC22A12 locus for association with risk of gout in NZ sample sets.

METHODS

A total of 12 single nucleotide polymorphism (SNP) variants in four haplotype blocks were genotyped using TaqMan® and Sequenom MassArray in 1003 gout cases and 1156 controls. All cases had gout according to the 1977 American Rheumatism Association criteria. Association analysis of single markers and haplotypes was performed using PLINK and Stata.

RESULTS

A haplotype block 1 SNP (rs17299124) (upstream of SLC22A11) was associated with gout in less admixed Polynesian sample sets, but not European Caucasian (odds ratio; OR = 3.38, P = 6.1 × 10-4; OR = 0.91, P = 0.40, respectively) sample sets. A protective block 1 haplotype caused the rs17299124 association (OR = 0.28, P = 6.0 × 10-4). Within haplotype block 2 (SLC22A11) we could not replicate previous reports of association of rs2078267 with gout in European Caucasian (OR = 0.98, P = 0.82) sample sets, however this SNP was associated with gout in Polynesian (OR = 1.51, P = 0.022) sample sets. Within haplotype block 3 (including SLC22A12) analysis of haplotypes revealed a haplotype with trans-ancestral protective effects (OR = 0.80, P = 0.004), and a second haplotype conferring protection in less admixed Polynesian sample sets (OR = 0.63, P = 0.028) but risk in European Caucasian samples (OR = 1.33, P = 0.039).

CONCLUSIONS

Our analysis provides evidence for multiple ancestral-specific effects across the SLC22A11/SLC22A12 locus that presumably influence the activity of OAT4 and URAT1 and risk of gout. Further fine mapping of the association signal is needed using trans-ancestral re-sequence data.

Linkage disequilibrium and haplotype block structure in a composite beef cattle breed

Background

The development of linkage disequilibrium (LD) maps and the characterization of haplotype block structure at the population level are useful parameters for guiding genome wide association (GWA) studies, and for understanding the nature of non-linear association between phenotypes and genes. The elucidation of haplotype block structure can reduce the information of several single nucleotide polymorphisms (SNP) into the information of a haplotype block, reducing the number of SNPs in a coherent way for consideration in GWA and genomic selection studies.

Results

The maximum average LD, measured by r² varied between 0.33 to 0.40 at a distance of < 2.5 kb, and the minimum average values of r² varied between 0.05 to 0.07 at distances ranging from 400 to 500 kb, clearly showing that the average r² reduced with the increase in SNP pair distances. The persistence of LD phase showed higher values at shorter genomic distances, decreasing with the increase in physical distance, varying from 0.96 at a distance of < 2.5 kb to 0.66 at a distance from 400 to 500 kb. A total of 78% of all SNPs were clustered into haplotype blocks, covering 1,57 Mb of the total autosomal genome size.

Conclusions

This study presented the first high density linkage disequilibrium map and haplotype block structure for a composite beef cattle population, and indicates that the high density SNP panel over 700 k can be used for genomic selection implementation and GWA studies for Canchim beef cattle.

Antihypertensive therapy in preeclampsia is not modulated by VEGF polymorphisms

INTRODUCTION

Vascular endothelial growth factor (VEGF) is relevant for healthy pregnancy, and abnormalities in VEGF functions have been associated with hypertensive disorders of pregnancy. Our group recently demonstrated that VEGF genetic polymorphisms affect the susceptibility to preeclampsia (PE).

OBJECTIVE

Therefore, in this study our aim is to examine whether VEGF polymorphisms affect the antihypertensive responses in women with PE.

METHODS

We studied 113 white PE women who were stratified according to blood pressure levels after antihypertensive treatment (46 responsive, R group and 67 non-responsive, NR group). We then compared the frequencies of two VEGF genetic polymorphisms (C-2578A and G-634C) between R and NR groups.

RESULTS

We found no significant differences in genotype or allele distributions between R and NR groups (P > 0.05). In addition, no difference was observed in overall distribution of haplotypes (P > 0.05).

CONCLUSION

Our data suggest that VEGF polymorphisms do not affect responsiveness to the antihypertensive therapy in PE.