Genetic Variants Associated with Gestational Hypertriglyceridemia and Pancreatitis

Severe hypertriglyceridemia is a well-known cause of pancreatitis. Usually, there is a moderate increase in plasma triglyceride level during pregnancy. Additionally, certain pre-existing genetic traits may render a pregnant woman susceptible to development of severe hypertriglyceridemia and pancreatitis, especially in the third trimester. To elucidate the underlying mechanism of gestational hypertriglyceridemic pancreatitis, we undertook DNA mutation analysis of the lipoprotein lipase (LPL), apolipoprotein C2 (APOC2), apolipoprotein A5 (APOA5), lipase maturation factor 1 (LMF1), and glycosylphosphatidylinositol-anchored high-density lipoprotein-binding protein 1 (GPIHBP1) genes in five unrelated pregnant Chinese women with severe hypertriglyceridemia and pancreatitis. DNA sequencing showed that three out of five patients had the same homozygous variation, p.G185C, in APOA5 gene. One patient had a compound heterozygous mutation, p.A98T and p.L279V, in LPL gene. Another patient had a compound heterozygous mutation, p.A98T & p.C14F in LPL and GPIHBP1 gene, respectively. No mutations were seen in APOC2 or LMF1 genes. All patients were diagnosed with partial LPL deficiency in non-pregnant state. As revealed in our study, genetic variants appear to play an important role in the development of severe gestational hypertriglyceridemia, and, p.G185C mutation in APOA5 gene appears to be the most common variant implicated in the Chinese population. Antenatal screening for mutations in susceptible women, combined with subsequent interventions may be invaluable in the prevention of potentially life threatening gestational hypertriglyceridemia-induced pancreatitis.     

Integrative Analysis Reveals Relationships of Genetic and Epigenetic Alterations in Osteosarcoma

Background

Osteosarcomas are the most common non-haematological primary malignant tumours of bone, and all conventional osteosarcomas are high-grade tumours showing complex genomic aberrations. We have integrated genome-wide genetic and epigenetic profiles from the EuroBoNeT panel of 19 human osteosarcoma cell lines based on microarray technologies.

Principal Findings

The cell lines showed complex patterns of DNA copy number changes, where genomic copy number gains were significantly associated with gene-rich regions and losses with gene-poor regions. By integrating the datasets, 350 genes were identified as having two types of aberrations (gain/over-expression, hypo-methylation/over-expression, loss/under-expression or hyper-methylation/under-expression) using a recurrence threshold of 6/19 (>30%) cell lines. The genes showed in general alterations in either DNA copy number or DNA methylation, both within individual samples and across the sample panel. These 350 genes are involved in embryonic skeletal system development and morphogenesis, as well as remodelling of extracellular matrix. The aberrations of three selected genes, CXCL5, DLX5 and RUNX2, were validated in five cell lines and five tumour samples using PCR techniques. Several genes were hyper-methylated and under-expressed compared to normal osteoblasts, and expression could be reactivated by demethylation using 5-Aza-2′-deoxycytidine treatment for four genes tested; AKAP12, CXCL5, EFEMP1 and IL11RA. Globally, there was as expected a significant positive association between gain and over-expression, loss and under-expression as well as hyper-methylation and under-expression, but gain was also associated with hyper-methylation and under-expression, suggesting that hyper-methylation may oppose the effects of increased copy number for detrimental genes.

Conclusions

Integrative analysis of genome-wide genetic and epigenetic alterations identified dependencies and relationships between DNA copy number, DNA methylation and mRNA expression in osteosarcomas, contributing to better understanding of osteosarcoma biology.     

Genetic Variants Associated with Increased Risk of Malignant Pleural Mesothelioma: A Genome-Wide Association Study

Asbestos exposure is the main risk factor for malignant pleural mesothelioma (MPM), a rare aggressive tumor. Nevertheless, only 5–17% of those exposed to asbestos develop MPM, suggesting the involvement of other environmental and genetic risk factors.To identify the genetic risk factors that may contribute to the development of MPM, we conducted a genome-wide association study (GWAS; 370,000 genotyped SNPs, 5 million imputed SNPs) in Italy, among 407 MPM cases and 389 controls with a complete history of asbestos exposure. A replication study was also undertaken and included 428 MPM cases and 1269 controls from Australia.Although no single marker reached the genome-wide significance threshold, several associations were supported by haplotype-, chromosomal region-, gene- and gene-ontology process-based analyses. Most of these SNPs were located in regions reported to harbor aberrant alterations in mesothelioma (SLC7A14, THRB, CEBP350, ADAMTS2, ETV1, PVT1 and MMP14 genes), causing at most a 2–3-fold increase in MPM risk. The Australian replication study showed significant associations in five of these chromosomal regions (3q26.2, 4q32.1, 7p22.2, 14q11.2, 15q14).Multivariate analysis suggested an independent contribution of 10 genetic variants, with an Area Under the ROC Curve (AUC) of 0.76 when only exposure and covariates were included in the model, and of 0.86 when the genetic component was also included, with a substantial increase of asbestos exposure risk estimation (odds ratio, OR: 45.28, 95% confidence interval, CI: 21.52–95.28).These results showed that genetic risk factors may play an additional role in the development of MPM, and that these should be taken into account to better estimate individual MPM risk in individuals who have been exposed to asbestos.     

Genetic Variants and Increased Expression of Parascaris equorum P-glycoprotein-11 in Populations with Decreased Ivermectin Susceptibility

Macrocyclic lactones (MLs) represent the major drug class for control of parasitic infections in humans and animals. However, recently reports of treatment failures became more frequent. In addition to human and ruminant parasitic nematodes this also is the case for the horse-nematode Parascaris equorum. Nevertheless, to date the molecular basis of ML resistance is still not understood. Unspecific resistance mechanisms involving transporters such as P-glycoproteins (Pgps) are expected to contribute to ML resistance in nematodes. Here, complete sequences of two P. equorum Pgps were cloned and identified as orthologs of Caenorhabditis elegans Ppg-11 and an unnamed Caenorhabditis briggsae Pgp designated as Pgp-16 using phylogenetic analysis. Quantitative real-time PCR was used to compare expression between tissues. Significantly higher PeqPgp-11 expression was found in the gut for both genders, whereas for PeqPgp-16 the body wall was identified as predominant expression site. Furthermore, Pgps were analyzed regarding their participation in resistance development. Using SeqDoC analyses, Pgp-sequences of P. equorum populations with different ML susceptibility were compared. This approach revealed three single nucleotide polymorphisms (SNPs) causing missense mutations in the PeqPgp-11 sequence which correlated with decreased ML susceptibility. However, no resistance associated differences in mRNA expression levels were detected between embryonated eggs of these populations. In contrast, comparison of two pre-adult groups with different ivermectin (IVM) susceptibility revealed the presence of the three SNPs and in addition statistically significant PeqPgp-11 overexpression in the group of worms with reduced susceptibility. These results indicate that Pgp-11 might be involved in IVM resistance in P. equorum as it shows increased expression in an IVM exposed life-cycle stage of an IVM resistant population as well as occurrence of putatively resistance associated SNPs in populations with reduced IVM susceptibility. These SNPs are promising diagnostic candidates for detection of ML resistance with potential also for other parasitic nematode species.     

Meta-analysis of Dense Genecentric Association Studies Reveals Common and Uncommon Variants Associated with Height

Height is a classic complex trait with common variants in a growing list of genes known to contribute to the phenotype. Using a genecentric genotyping array targeted toward cardiovascular-related loci, comprising 49,320 SNPs across approximately 2000 loci, we evaluated the association of common and uncommon SNPs with adult height in 114,223 individuals from 47 studies and six ethnicities. A total of 64 loci contained a SNP associated with height at array-wide significance (p < 2.4 × 10⁻⁶), with 42 loci surpassing the conventional genome-wide significance threshold (p < 5 × 10⁻⁸). Common variants with minor allele frequencies greater than 5% were observed to be associated with height in 37 previously reported loci. In individuals of European ancestry, uncommon SNPs in IL11 and SMAD3, which would not be genotyped with the use of standard genome-wide genotyping arrays, were strongly associated with height (p < 3 × 10⁻¹¹). Conditional analysis within associated regions revealed five additional variants associated with height independent of lead SNPs within the locus, suggesting allelic heterogeneity. Although underpowered to replicate findings from individuals of European ancestry, the direction of effect of associated variants was largely consistent in African American, South Asian, and Hispanic populations. Overall, we show that dense coverage of genes for uncommon SNPs, coupled with large-scale meta-analysis, can successfully identify additional variants associated with a common complex trait.     

Genetic Variants of CD209 Associated with Kawasaki Disease Susceptibility

Background

Kawasaki disease (KD) is a systemic vasculitis with unknown etiology mainly affecting children in Asian countries. Dendritic cell-specific intercellular adhesion molecule-3 grabbing non-integrin (DC-SIGN, CD209) in humans was showed to trigger an anti-inflammatory cascade and associated with KD susceptibility. This study was conducted to investigate the association between genetic polymorphisms of CD209 and the risk KD.

Methods

A total of 948 subjects (381 KD and 567 controls) were recruited. Nine tagging SNPs (rs8112310, rs4804800, rs11465421, rs1544766, rs4804801, rs2287886, rs735239, rs735240, rs4804804) were selected for TaqMan allelic discrimination assay. Clinical phenotypes, coronary artery lesions (CAL) and intravenous immunoglobulin (IVIG) treatment outcomes were collected for analysis.

Results

Significant associations were found between CD209 polymorphisms (rs4804800, rs2287886, rs735240) and the risk of KD. Haplotype analysis for CD209 polymorphisms showed that A/A/G haplotype (P = 0.0002, OR = 1.61) and G/A/G haplotype (P = 0.0365, OR = 1.52) had higher risk of KD as compared with G/G/A haplotype in rs2287886/rs735239/rs735240 pairwise allele analysis. There were no significant association in KD with regards to CAL formation and IVIG treatment responses.

Conclusion

CD209 polymorphisms were responsible for the susceptibility of KD, but not CAL formation and IVIG treatment responsiveness.     

Host Genetic Variants and Gene Expression Patterns Associated with Epstein-Barr Virus Copy Number in Lymphoblastoid Cell Lines

Lymphoblastoid cell lines (LCLs) are commonly used in molecular genetics, supplying DNA for the HapMap and 1000 Genomes Projects, used to test chemotherapeutic agents, and informing the basis of a number of population genetics studies of gene expression. The process of transforming human B cells into LCLs requires the presence of Epstein-Barr virus (EBV), a double-stranded DNA virus which through B-cell immortalisation maintains an episomal virus genome in every cell of an LCL at variable copy numbers. Previous studies have reported that EBV alters host-gene expression and EBV copy number may be under host genetic control. We performed a genome-wide association study of EBV genome copy number in LCLs and found the phenotype to be highly heritable, although no individual SNPs achieved a significant association with EBV copy number. The expression of two host genes (CXCL16 and AGL) was positively correlated and expression of ADARB2 was negatively correlated with EBV copy number in a genotype-independent manner. This study shows an association between EBV copy number and the gene expression profile of LCLs, and suggests that EBV copy number should be considered as a covariate in future studies of host gene expression in LCLs.     

Analysis of Pharmacogenomic Variants Associated with Population Differentiation

In the present study, we systematically investigated population differentiation of drug-related (DR) genes in order to identify common genetic features underlying population-specific responses to drugs. To do so, we used the International HapMap project release 27 Data and Pharmacogenomics Knowledge Base (PharmGKB) database. First, we compared four measures for assessing population differentiation: the chi-square test, the analysis of variance (ANOVA) F-test, F_st, and Nearest Shrunken Centroid Method (NSCM). F_st showed high sensitivity with stable specificity among varying sample sizes; thus, we selected F_st for determining population differentiation. Second, we divided DR genes from PharmGKB into two groups based on the degree of population differentiation as assessed by F_st: genes with a high level of differentiation (HD gene group) and genes with a low level of differentiation (LD gene group). Last, we conducted a gene ontology (GO) analysis and pathway analysis. Using all genes in the human genome as the background, the GO analysis and pathway analysis of the HD genes identified terms related to cell communication. “Cell communication” and “cell-cell signaling” had the lowest Benjamini-Hochberg’s q-values (0.0002 and 0.0006, respectively), and “drug binding” was highly enriched (16.51) despite its relatively high q-value (0.0142). Among the 17 genes related to cell communication identified in the HD gene group, five genes (STX4, PPARD, DCK, GRIK4, and DRD3) contained single nucleotide polymorphisms with F_st values greater than 0.5. Specifically, the F_st values for rs10871454, rs6922548, rs3775289, rs1954787, and rs167771 were 0.682, 0.620, 0.573, 0.531, and 0.510, respectively. In the analysis using DR genes as the background, the HD gene group contained six significant terms. Five were related to reproduction, and one was “Wnt signaling pathway,” which has been implicated in cancer. Our analysis suggests that the HD gene group from PharmGKB is associated with cell communication and drug binding.     

Evaluation of Replication of Variants Associated with Genetic Risk of Otitis Media

The first Genome Wide Association Study (GWAS) of otitis media (OM) found evidence of association in the Western Australian Pregnancy Cohort (Raine) study, but lacked replication in an independent OM population. The aim of this study was to investigate association at these loci in our family-based sample of chronic otitis media with effusion and recurrent otitis media (COME/ROM). Autosomal SNPs were selected from the Raine OM GWAS results. SNPs from the Raine cohort GWAS genotyped in our GWAS of COME/ROM had P-values ranging from P = 0.06–0.80. After removal of SNPs previously genotyped in our GWAS of COME/ROM (N = 21) and those that failed Fluidigm assay design (N = 1), 26 SNPs were successfully genotyped in 716 individuals from our COME/ROM family population. None of the SNP associations replicated in our family-based population (unadjusted P = 0.03–0.93). Replication in an independent sample would confirm that these represent novel OM loci, and that further investigation is warranted.     

Proteomic Analysis of MG132-Treated Germinating Pollen Reveals Expression Signatures Associated with Proteasome Inhibition

Chemical inhibition of the proteasome has been previously found to effectively impair pollen germination and tube growth in vitro. However, the mediators of these effects at the molecular level are unknown. By performing 2DE proteomic analysis, 24 differentially expressed protein spots, representing 14 unique candidate proteins, were identified in the pollen of kiwifruit (Actinidia deliciosa) germinated in the presence of the MG132 proteasome inhibitor. qPCR analysis revealed that 11 of these proteins are not up-regulated at the mRNA level, but are most likely stabilized by proteasome inhibition. These differentially expressed proteins are predicted to function in various pathways including energy and lipid metabolism, cell wall synthesis, protein synthesis/degradation and stress responses. In line with this evidence, the MG132-induced changes in the proteome were accompanied by an increase in ATP and ROS content and by an alteration in fatty acid composition.     

Genetic Variants Associated with Myocardial Infarction and the Risk Factors in Chinese Population

Background

Recent genome-wide association (GWA) studies in Caucasians identified multiple single nucleotide polymorphisms (SNPs) associated with coronary artery disease (CAD). The associations of those SNPs with myocardial infarction (MI) have not been replicated in Asian populations. Among those previously identified SNPs, we selected nine (rs10953541, rs1122608, rs12190287, rs12413409, rs1412444, rs1746048, rs3798220, rs4977574, rs579459, in or near genes 7q22, LDLR, TCF21, CYP17A1, LIPA, CXCL12, LPA, CDKN2A, ABO, respectively) because of the relatively high minor allele frequencies in Chinese individuals and tested the associations of the SNPs with MI and MI related risk factors in Chinese population.

Methods and Results

We conducted a case–control association study on a cohort of 2365 MI patients and 2678 unrelated controls from the Chinese population. Genotyping of 9 SNPs were performed by the TaqMan Real Time PCR method. After age, sex, and BMI adjustment, we observed the SNPs rs12190287, rs12413409, rs1412444, rs1746048 and rs4977574, were significantly associated with MI in additive models and rs12190287, rs12413409, rs4977574 were significantly associated with phenotypes of MI at the same time. We also found three SNPs rs1122608, rs3798220 and rs579459 were significantly associated with risk factors of MI, although they had no association with MI in Chinese population.

Conclusion

Results of this study indicate that 5 SNPs were associated with MI and 3 SNPs were associated with associated with lipoprotein levels but not with MI in a Chinese population. The present study supports some CAD-related genes in Caucasian as important genes for MI in a Chinese population.     

Mitochondrial GWA Analysis of Lipid Profile Identifies Genetic Variants to Be Associated with HDL Cholesterol and Triglyceride Levels

It has been suggested that mitochondrial dysfunction has an influence on lipid metabolism. The fact that mitochondrial defects can be accumulated over time as a normal part of aging may explain why cholesterol levels often are altered with age. To test the hypothesis whether mitochondrial variants are associated with lipid profile (total cholesterol, LDL, HDL, and triglycerides) we analyzed a total number of 978 mitochondrial single nucleotide polymorphisms (mtSNPs) in a sample of 2,815 individuals participating in the population-based KORA F4 study. To assess mtSNP association while taking the presence of heteroplasmy into account we used the raw signal intensity values measured on the microarray and applied linear regression. Ten mtSNPs (mt3285, mt3336, mt5285, mt6591, mt6671, mt9163, mt13855, mt13958, mt14000, and mt14580) were significantly associated with HDL cholesterol and one mtSNP (mt15074) with triglycerides levels. These results highlight the importance of the mitochondrial genome among the factors that contribute to the regulation of lipid levels. Focusing on mitochondrial variants may lead to further insights regarding the underlying physiological mechanisms, or even to the development of innovative treatments. Since this is the first mitochondrial genome-wide association analysis (mtGWAS) for lipid profile, further analyses are needed to follow up on the present findings.