Using genome-wide associations to identify metabolic pathways involved in maize aflatoxin accumulation resistance

Background

Aflatoxin is a potent carcinogen that can contaminate grain infected with the fungus Aspergillus flavus. However, resistance to aflatoxin accumulation in maize is a complex trait with low heritability. Here, two complementary analyses were performed to better understand the mechanisms involved. The first coupled results of a genome-wide association study (GWAS) that accounted for linkage disequilibrium among single nucleotide polymorphisms (SNPs) with gene-set enrichment for a pathway-based approach. The rationale was that the cumulative effects of genes in a pathway would give insight into genetic differences that distinguish resistant from susceptible lines of maize. The second involved finding non-pathway genes close to the most significant SNP-trait associations with the greatest effect on reducing aflatoxin in multiple environments. Unlike conventional GWAS, the latter analysis emphasized multiple aspects of SNP-trait associations rather than just significance and was performed because of the high genotype x environment variability exhibited by this trait.

Results

The most significant metabolic pathway identified was jasmonic acid (JA) biosynthesis. Specifically, there was at least one allelic variant for each step in the JA biosynthesis pathway that conferred an incremental decrease to the level of aflatoxin observed among the inbred lines in the GWAS panel. Several non-pathway genes were also consistently associated with lowered aflatoxin levels. Those with predicted functions related to defense were: leucine-rich repeat protein kinase, expansin B3, reversion-to-ethylene sensitivity1, adaptor protein complex2, and a multidrug and toxic compound extrusion protein.

Conclusions

Our genetic analysis provided strong evidence for several genes that were associated with aflatoxin resistance. Inbred lines that exhibited lower levels of aflatoxin accumulation tended to share similar haplotypes for genes specifically in the pathway of JA biosynthesis, along with several non-pathway genes with putative defense-related functions. Knowledge gained from these two complementary analyses has improved our understanding of population differences in aflatoxin resistance.

Electronic supplementary material

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

SCtool:识别复杂性状和疾病间遗传关联的工具

识别复杂性状和疾病间遗传关联可以提供有用的病因学见解,并有助于确定可能的因果关系的优先级。尽管已有很多工具可以实现复杂性状和疾病间遗传关联,但是某些工具代码可读性差、并且不同工具基于不同的计算机语言、工具间的串联性较差。因此,本研究基于全基因组关联研究(GWAS)数据,提出了SCtool,一个开源、跨平台和用户友好的软件工具。SCtool整合了ldsc, TwosampleMR和MR-BMA三种软件,其主要功能是基于GWAS汇总水平的数据,识别复杂性状和疾病、复杂性状和复杂性状以及疾病与疾病间的遗传相关性并探究其间潜在的因果关联。最后,使用SCtool揭示了全身性铁状态(铁蛋白,血清铁,转铁蛋白,转铁蛋白饱和度)与表观遗传时钟GrimAge之间的遗传关联。     

The use of genetic correlations to evaluate associations between SNP markers and quantitative traits

Open-pollinated progeny of Corymbia citriodora established in replicated field trials were assessed for stem diameter, wood density, and pulp yield prior to genotyping single nucleotide polymorphisms (SNP) and testing the significance of associations between markers and assessment traits. Multiple individuals within each family were genotyped and phenotyped, which facilitated a comparison of standard association testing methods and an alternative method developed to relate markers to additive genetic effects. Narrow-sense heritability estimates indicated there was significant additive genetic variance within this population for assessment traits ( $ {\widehat{h}^{{2}}} = 0.{28}\;{\text{to}}\;0.{44} $ ) and genetic correlations between the three traits were negligible to moderate (r _G?=?0.08 to 0.50). The significance of association tests (p values) were compared for four different analyses based on two different approaches: (1) two software packages were used to fit standard univariate mixed models that include SNP-fixed effects, (2) bivariate and multivariate mixed models including each SNP as an additional selection trait were used. Within either the univariate or multivariate approach, correlations between the tests of significance approached +1; however, correspondence between the two approaches was less strong, although between-approach correlations remained significantly positive. Similar SNP markers would be selected using multivariate analyses and standard marker-trait association methods, where the former facilitates integration into the existing genetic analysis systems of applied breeding programs and may be used with either single markers or indices of markers created with genomic selection processes.     

DNA repair in lymphoblastoid cell lines established from human genetic disorders

Lymphoblastoid cell lines (LCLs) established from chromosomal breakage syndromes or related genetic disorders have been used to study the effects of mutagens on human lymphoid cells. The disorders studied include xeroderma pigmentosum, ataxia telangiectasia, Fanconi's anemia, Bloom's syndrome and Cockayne's syndrome. Three approaches were used to assess the cells' ability to cope with a particular mutagen: (1) assaying recovery of DNA snythetic capabilities as measured by [³H]thymidine (dT) incorporation; (2) measurements of classical excision DNA repair by isopyknic sedimentation of DNA density labeled with 5-bromo-2-deoxyuridine (BrdU); (3) determining cell survival by colony formation in microtiter plates. LCLs established from xeroderma pigmentosum showed increased sensitivities to ultraviolet (354 nm) light and N-acetoxy-2-acetylaminofluorene (AAAF) as determined by DNA synthesis or colony formation and had diminished levels of excision-repair. Cockayne's syndrome LCLs, on the other hand, had increased sensitivities to ultraviolet (UV) light, AAAF and N-methyl-N′-nitro-N-nitrosoguanidine (MNNG) while showing near normal levels of DNA-repair after treatment with each agent. An LCL established from ataxia telangiectasia had decreased DNA repair synthesis and defective colony-forming ability following treatment with MNNG. LCLs, in addition to ease of establishment, appear likely to provide useful material for the study of DNA repair replication and its relationship to carcinogenesis.     

The use of LC-MS to identify differentially expressed proteins in docetaxel-resistant prostate cancer cell lines

Docetaxel is a taxane-derived chemotherapy drug that has been approved for treatment of prostate cancer. While docetaxel is frequently used as a treatment for hormone-refractory prostate cancer, a subset of patients either do not respond to this treatment or those that do respond eventually become resistant to the drug over time. Resistance to docetaxel is complex and multi-factoral and further understanding of the cellular biochemistry underlying resistance is vital to improve treatment efficacy. To identify proteins altered in the resistant phenotype, three parental cell lines DU145, 22RV1 and PC-3, as well as their docetaxel resistant sub-lines, were subjected to quantitative label-free LC-MS proteomic profiling. A total of 189 significant (p < 0.05) protein abundance changes were identified in the DU145 resistant sub-lines, 254 in the 22RV1 sub-lines, and 51 and 72 in the 8 and 12 nM resistant PC-3 sub-lines, respectively. From these, 29 proteins demonstrated a significant (p < 0.05) fold change across two or more resistant variants. These included proteins indicative of an epithelial-to-mesenchemyl transition as well as altered heat shock response elements.     

MicroRNA signatures associated with immortalization of EBV-transformed lymphoblastoid cell lines and their clinical traits

Objective: MicroRNAs (miRNAs) are negative regulators of gene expression that play important roles in cell processes such as proliferation, development and differentiation. Recently, it has been reported that miRNAs are related to development of carcinogenesis. The aim of this study was to identify miRNAs associated with terminal immortalization of Epstein–Barr virus (EBV)‐transformed lymphoblastoid cell line (LCL) and associated clinical traits. Material and Methods: Hence, we performed miRNA microarray approach with early‐ (p6) and late‐passage (p161) LCLs. Results and Conclusion: Microarray data showed that nine miRNAs (miR‐20b*, miR‐28‐5p, miR‐99a, miR‐125b, miR‐151‐3p, miR‐151:9.1, miR‐216a, miR‐223* and miR‐1296) were differentially expressed in most LCLs during long‐term culture. In particular, miR‐125b was up‐regulated in all the tested late‐passage LCLs. miR‐99a, miR‐125b, miR‐216a and miR‐1296 were putative negative regulators of RASGRP3, GPR160, PRKCH and XAF1, respectively, which were found to be differentially expressed in LCLs during long‐term culture in a previous study. Linear regression analysis showed that miR‐200a and miR‐296‐3p correlated with triglyceride and HbA1C levels, respectively, suggesting that miRNA signatures of LCLs could provide information on the donor’s health. In conclusion, our study suggests that expression changes of specific miRNAs may be required for terminal immortalization of LCLs. Thus, differentially expressed miRNAs would be a potential marker for completion of cell immortalization during EBV‐mediated tumorigenesis.     

Molecular signatures of cardiac defects in down syndrome lymphoblastoid cell lines suggest altered ciliome and hedgehog pathways

Forty percent of people with Down syndrome exhibit heart defects, most often an atrioventricular septal defect (AVSD) and less frequently a ventricular septal defect (VSD) or atrial septal defect (ASD). Lymphoblastoid cell lines (LCLs) were established from lymphocytes of individuals with trisomy 21, the chromosomal abnormality causing Down syndrome. Gene expression profiles generated from DNA microarrays of LCLs from individuals without heart defects (CHD(-); n?=?22) were compared with those of LCLs from patients with cardiac malformations (CHD(+); n?=?21). After quantile normalization, principal component analysis revealed that AVSD carriers could be distinguished from a combined group of ASD or VSD (ASD+VSD) carriers. From 9,758 expressed genes, we identified 889 and 1,016 genes differentially expressed between CHD(-) and AVSD and CHD(-) and ASD+VSD, respectively, with only 119 genes in common. A specific chromosomal enrichment was found in each group of affected genes. Among the differentially expressed genes, more than 65% are expressed in human or mouse fetal heart tissues (GEO dataset). Additional LCLs from new groups of AVSD and ASD+VSD patients were analyzed by quantitative PCR; observed expression ratios were similar to microarray results. Analysis of GO categories revealed enrichment of genes from pathways regulating clathrin-mediated endocytosis in patients with AVSD and of genes involved in semaphorin-plexin-driven cardiogenesis and the formation of cytoplasmic microtubules in patients with ASD-VSD. A pathway-oriented search revealed enrichment in the ciliome for both groups and a specific enrichment in Hedgehog and Jak-stat pathways among ASD+VSD patients. These genes or related pathways are therefore potentially involved in normal cardiogenesis as well as in cardiac malformations observed in individuals with trisomy 21.     

PEGylation-based strategy to identify pathways involved in the activation of apoptotic BAX protein

BackgroundBAX activation is a crucial step for commitment to apoptosis. Several activators, such as BimBH3-based therapeutic peptides and cleaved Bid (cBid) protein, can trigger BAX-mediated apoptosis, but it is unclear whether they proceed through the same pathway.MethodsHere we utilize PEGylation-based approach, which is shown to efficiently shield individual binding grooves in BAX from activators, to investigate and reveal that the activators take different routes to induce BAX-mediated apoptosis. Various spectroscopic/biochemical tools, including electron spin resonance, circular dichroism, fluorescence recovery after photobleaching, and label-transfer assay, were employed to reveal details in the processes.ResultsWe observe a key mutant BAX 164-PEG that acts differently in response to cBid and BimBH3 stimuli. While BimBH3 directly interacts with the trigger groove (TG) to induce the conformational changes in BAX that includes the release of α9 from the canonical groove (CG) and oligomerization, cBid engages with CG and works with mitochondrial lipids to fully activate BAX.ConclusionPEGylation-based approach is proven useful to shield individual binding grooves of BAX from apoptotic stimuli. Groove engagement in CG of BAX is required for a full cBid-induced BAX activation. This study has identified differences in the pathways involved during the initiation of BAX activation by full-length cBid protein versus synthetic BimBH3-based peptides.General significanceOur finding is potentially valuable for therapeutic application as the pore-forming activity of 164-PEG is independent from the cBid-mediated apoptotic pathways, but can be administrated by the synthetic short peptides.     

Instability of Mex- phenotype in human lymphoblastoid cell lines

Three lymphoblastoid cell lines (LCLs) had extremely low activities of O6-alkylguanine-DNA alkyltransferase (O6-AGT), and were classified as Mex-. They were highly sensitive to cell killing by 1-(4-amino-2-methyl-5-pyrimidinyl)-methyl-3-(2-chloroethyl)-3-nitrosoure a hydrochloride (ACNU), whereas NMO2, a Mex+ LCL with a high O6-AGT activity, was resistant to the agent. Small fractions of these Mex- LCLs survived the treatment with 10 micrograms/ml of ACNU for 24 h, and the surviving cells were found to be resistant to subsequent treatments with the agent. In addition, they contained O6-AGT activities comparable to that of NMO2 and were therefore regarded as Mex+. These results suggest that the Mex- phenotype in LCLs is unstable.     

Binding of gliadin to lymphoblastoid,myeloid and epithelial cell lines

The aim of our work was to investigate thein vitro reactivity of gliadin peptides of natural and synthetic origin with various cell lines. We have found that all tested cell lines of human, mouse and rat origin were agglutinated by enzymically digested gliadin (peptic-tryptic-and peptic-tryptic pancreatic digest of α-gliadin) in a concentration dependent manner. In order to test the specificity of binding, inhibition studies were performed using a panel of sugars as well as natural and synthetic peptides derived from gliadin. We have found that among twelve tested sugars only fetuin and phosphomannan were able to inhibition by gliadin peptides and most of the saccharides suggests that agglutinating activity of gliadin is the result of a nonspecific binding of gliadin to the cell membrane. Dedicated to Professor J. Šterzl on the occasion of his 70th birthday     

The use of enzyme polymorphisms to identify genetic differences in the genus Trichinella

Seven isolates of the genus Trichinella were evaluated for their genetic identity and variability using starch gel electrophoresis. The various strains exhibited polymorphism in 8 of 11 enzyme systems tested and the strains could be segregated into at least 3 basic patterns for geographical isolates of Trichinella: pigs, wild carnivores and Trichinella spiralis var. pseudospiralis. There was also evidence for further differences between isolates from wild carnivores. The isozyme technique appears to have potential for zoogeographic studies in the genus Trichinella.