Investigation of putative roles of smoking-associated salivary microbiome alterations on carcinogenesis by integrative in silico analysis

Growing evidence suggests that cigarette smoking alters the salivary microbiome composition and affects the risk of various complex diseases including cancer. However, the potential role of the smoking-associated microbiome in cancer development remains unexplained. Here, the putative roles of smoking-related microbiome alterations in carcinogenesis were investigated by in silico analysis and suggested evidence can be further explored by experimental methodologies. The Disbiome database was used to extract smoking-associated microbial taxa in saliva and taxon set enrichment analysis (TSEA) was conducted to identify the gene sets associated with extracted microbial taxa. We further analyzed the expression profiles of identified genes by using RNA-sequencing data from TCGA and GTEx projects. Associations of the genes with smoking-related phenotypes in cancer datasets were analyzed to prioritize genes for their interplay between smoking-related microbiome and carcinogenesis. Thirty-eight microbial taxa associated with smoking were included in the TSEA and this revealed sixteen genes that were significantly associated with smoking-associated microbial taxa. All genes were found to be differentially expressed in at least one cancer dataset, yet the ELF3 and CTSH were the most common differentially expressed genes giving significant results for several cancer types. Moreover, C2CD3, CTSH, DSC3, ELF3, RHOT2, and WSB2 showed statistically significant associations with smoking-related phenotypes in cancer datasets. This study provides in silico evidence for the potential roles of the salivary microbiome on carcinogenesis. The results shed light on the importance of smoking cessation strategies for cancer management and interventions to stratify smokers for their risk of smoking-induced carcinogenesis.

Integrative analyses of molecular pathways and key candidate biomarkers associated with colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the leading causes of cancer-related deaths and mining the molecular factors underlying CRC pathogenesis is imperative for alleviating the disease burden.

OBJECTIVE

To highlight key molecular pathways, prioritize hub genes and their regulators related to CRC.

METHODS

Data sets of TCGA-COAD and GTEx were used to identify differentially expressed genes (DEGs) and their functional enrichments in pathways and biological processes were analyzed using bioinformatics tools. Protein-protein interaction network was constructed and hub genes were identified using Cytoscape. Ingenuity Pathway Analysis was used to analyze the relations of the hub genes with diseases and canonical pathways. Key regulators targeting the hub genes such as TFs, miRNAs and their interactions were identified using in silico tools.

RESULTS

AURKA, CDK1, MYC, CDH1, CCNB1, CDC20, UBE2C, PLK1, KIF11, and CCNA2 were prioritized as hub genes based on their topological properties. Enrichment analyses emphasized the roles of DEGs and hub genes in the cell cycle process. Interactions of the hub genes with TFs and miRNAs suggested TP53, EZH2 and KLF4 as being promising candidate biomarkers for CRC.

CONCLUSIONS

Our results provide in silico evidence for candidate biomolecules that may have strong biomarker potential for CRC-related translational strategies.

Characterization and in silico analyses of the BRCA1/2 variants identified in individuals with personal and/or family history of BRCA-related cancers

Pathogenic variants in the coding regions of the BRCA1/2 lead dysfunctional or nonfunctional BRCA proteins however the contribution of non-coding BRCA1/2 variants to BRCA-related disease risk has not been fully elucidated. Thus, we characterized the functional impact of both coding and non-coding BRCA1/2 variants identified in individuals with personal and/or family history of BRCA-related cancers. The data were produced by resequencing the exons and exon-intron junctions of the BRCA1/2 in 125 individuals and were comprehensively analyzed by using bioinformatics tools and databases. A total of 96 variants (59 coding and 37 non-coding) including 7 novel variants were identified and analyzed for their functional importance. We identified 11 missense variants that potentially affect protein function; 22 variants were likely to alter different types of posttranslational modifications. Also, multiple non-coding BRCA1/2 variants were found to reside in the critical regulatory regions that have the potential to act as eQTLs and affect alternative splicing. The results of our study shed light on the possible contributions of not only coding variants but also non-coding BRCA1/2 variants in BRCRA-related cancers. Further investigation is required to fully understand their potential associations with phenotypes which may ultimately lead their utilization on cancer management as a biomarker.

Hsa-miR-584-5p as a novel candidate biomarker in Turkish men with severe coronary artery disease

Coronary artery disease (CAD) is still the preliminary cause of mortality and morbidity in the developed world. Identification of novel predictive and therapeutic biomarkers is crucial for accurate diagnosis, prognosis and treatment of the CAD. The aim of this study was to detect novel candidate miRNA biomarker that may be used in the management of CAD. We performed miRNA profiling in whole blood samples of angiographically confirmed Turkish men with CAD and non-CAD controls with insignificant coronary stenosis. Validation of microarray results was performed by qRT-PCR in a larger cohort of 62 samples. We subsequently assessed the diagnostic value of the miRNA and correlations of miRNA with clinical parameters. miRNA-target identification and network analyses were conducted by Ingenuity Pathway Analysis (IPA) software. Hsa-miR-584-5p was one of the top significantly dysregulated miRNA observed in miRNA microarray. Men-specific down-regulation (p = 0.040) of hsa-miR-584-5p was confirmed by qRT-PCR. ROC curve analysis highlighted the potential diagnostic value of hsa-miR-584-5p with a power area under the curve (AUC) of 0.714 and 0.643 in men and in total sample, respectively. The expression levels of hsa-miR-584-5p showed inverse correlation with stenosis and Gensini scores. IPA revealed CDH13 as the only CAD related predicted target for the miRNA with biological evidence of its involvement in CAD. This study suggests that hsa-miR-584-5p, known to be tumor suppressor miRNA, as a candidate biomarker for CAD and highlighted its putative role in the CAD pathogenesis. The validation of results in larger samples incorporating functional studies warrant further research.

Apolipoprotein E-C1-C4-C2 gene cluster region and inter-individual variation in plasma lipoprotein levels: a comprehensive genetic association study in two ethnic groups

The apolipoprotein E-C1-C4-C2 gene cluster at 19q13.32 encodes four amphipathic apolipoproteins. The influence of APOE common polymorphisms on plasma lipid/lipoprotein profile, especially on LDL-related traits, is well recognized; however, little is known about the role of other genes/variants in this gene cluster. In this study, we evaluated the role of common and uncommon/rare genetic variation in this gene region on inter-individual variation in plasma lipoprotein levels in non-Hispanic Whites (NHWs) and African blacks (ABs). In the variant discovery step, the APOE, APOC1, APOC4, APOC2 genes were sequenced along with their flanking and hepatic control regions (HCR1 and HCR2) in 190 subjects with extreme HDL-C/TG levels. The next step involved the genotyping of 623 NHWs and 788 ABs for the identified uncommon/rare variants and common tagSNPs along with additional relevant SNPs selected from public resources, followed by association analyses with lipid traits. A total of 230 sequence variants, including 15 indels were identified, of which 65 were novel. A total of 70 QC-passed variants in NHWs and 108 QC-passed variants in ABs were included in the final association analyses. Single-site association analysis of SNPs with MAF>1% revealed 20 variants in NHWs and 24 variants in ABs showing evidence of association with at least one lipid trait, including several variants exhibiting independent associations from the established APOE polymorphism even after multiple-testing correction. Overall, our study has confirmed known associations and also identified novel associations in this genomic region with various lipid traits. Our data also support the contribution of both common and uncommon/rare variation in this gene region in affecting plasma lipid profile in the general population.

Resequencing of the CETP gene in American whites and African blacks: Association of rare and common variants with HDL-cholesterol levels

BACKGROUND

Cholesteryl ester transfer protein (CETP) plays a crucial role in lipid metabolism. Associations of common CETP variants with variation in plasma lipid levels, and/or CETP mass/activity have been extensively studied and well-documented; however, the effects of uncommon/rare CETP variants on plasma lipid profile remain undefined. Hence, resequencing of the gene in extreme phenotypes and follow-up rare-variant association analyses are essential to fill this gap.

OBJECTIVE

To identify common and uncommon/rare variants in the CETP gene by resequencing the entire gene and test the effects of both common and uncommon/rare CETP variants on plasma lipid traits in two genetically distinct populations.

METHODS AND RESULTS

The entire CETP gene plus flanking regions were resequenced in 190 individuals comprising 95 non-Hispanic whites (NHWs) and 95 African blacks with extreme HDL-C levels. A total of 279 sequence variants were identified, of which 25 were novel. Selected variants were genotyped in the entire samples of 623 NHWs and 788 African blacks and 184 QC-passed variants were tested in relation to plasma lipid traits by using gene-based, single-site, haplotype and rare variant association analyses (SKAT-O). Two novel and independent associations of rs1968905 and rs289740 with HDL-C were identified in African blacks. Using SKAT-O analysis, we also identified rare variants with minor allele frequency <0.01 to be associated with HDL-C in both NHWs (P=0.024) and African blacks (P=0.009).

CONCLUSIONS

Our results point out that in addition to the common CETP variants, rare genetic variants in the CETP gene also contribute to the phenotypic variation of HDL-C in the general population.

Genetic contribution of SCARB1 variants to lipid traits in African Blacks: a candidate gene association study

Background

High-density lipoprotein cholesterol (HDL-C) exerts many anti-atherogenic properties including its role in reverse cholesterol transport (RCT). Scavenger receptor class B member 1 (SCARB1) plays a key role in RCT by selective uptake of HDL cholesteryl esters. We aimed to explore the genetic contribution of SCARB1 to affecting lipid levels in African Blacks from Nigeria.

Methods

We resequenced 13 exons and exon-intron boundaries of SCARB1 in 95 individuals with extreme HDL-C levels using Sanger method. Then, we genotyped 147 selected variants (78 sequence variants, 69 HapMap tagSNPs, and 2 previously reported relevant variants) in the entire sample of 788 African Blacks using either the iPLEX Gold or TaqMan methods. A total of 137 successfully genotyped variants were further evaluated for association with major lipid traits.

Results

The initial gene-based analysis demonstrated evidence of association with HDL-C and apolipoprotein A-I (ApoA-I). The follow-up single-site analysis revealed nominal evidence of novel associations of nine common variants with HDL-C and/or ApoA-I (P < 0.05). The strongest association was between rs11057851 and HDL-C (P = 0.0043), which remained significant after controlling for multiple testing using false discovery rate. Rare variant association testing revealed a group of 23 rare variants (frequencies ≤1 %) associated with HDL-C (P = 0.0478). Haplotype analysis identified four SCARB1 regions associated with HDL-C (global P < 0.05).

Conclusions

To our knowledge, this is the first report of a comprehensive association study of SCARB1 variations with lipid traits in an African Black population. Our results showed the consistent association of SCARB1 variants with HDL-C across various association analyses, supporting the role of SCARB1 in lipoprotein-lipid regulatory mechanism.

Electronic supplementary material

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

Resequencing of LPL in African Blacks and associations with lipoprotein-lipid levels

Genome-wide association studies have identified several loci associated with plasma lipid levels but those common variants together account only for a small proportion of the genetic variance of lipid traits. It has been hypothesized that the remaining heritability may partly be explained by rare variants with strong effect sizes. Here, we have comprehensively investigated the associations of both common and uncommon/rare variants in the lipoprotein lipase (LPL) gene in relation to plasma lipoprotein-lipid levels in African Blacks (ABs). For variant discovery purposes, the entire LPL gene and flanking regions were resequenced in 95 ABs with extreme high-density lipoprotein cholesterol (HDL-C) levels. A total of 308 variants were identified, of which 64 were novel. Selected common tagSNPs and uncommon/rare variants were genotyped in the entire sample (n=788), and 126 QC-passed variants were evaluated for their associations with lipoprotein-lipid levels by using single-site, haplotype and rare variant (SKAT-O) association analyses. We found eight not highly correlated (r(2)<0.40) signals (rs1801177:G>A, rs8176337:G>C, rs74304285:G>A, rs252:delA, rs316:C>A, rs329:A>G, rs12679834:T>C, and rs4921684:C>T) nominally (P<0.05) associated with lipid traits (HDL-C, LDL-C, ApoA1 or ApoB levels) in our sample. The most significant SNP, rs252:delA, represented a novel association observed with LDL-C (P=0.002) and ApoB (P=0.012). For TG and LDL-C, the haplotype analysis was more informative than the single-site analysis. The SKAT-O analysis revealed that the bin (group) containing 22 rare variants with MAF≤0.01 exhibited nominal association with TG (P=0.039) and LDL-C (P=0.027). Our study indicates that both common and uncommon/rare LPL variants/haplotypes may affect plasma lipoprotein-lipid levels in general African population.

Comprehensive evaluation of the association of APOE genetic variation with plasma lipoprotein traits in U.S. whites and African blacks

Although common APOE genetic variation has a major influence on plasma LDL-cholesterol, its role in affecting HDL-cholesterol and triglycerides is not well established. Recent genome-wide association studies suggest that APOE also affects plasma variation in HDL-cholesterol and triglycerides. It is thus important to resequence the APOE gene to identify both common and uncommon variants that affect plasma lipid profile. Here, we have sequenced the APOE gene in 190 subjects with extreme HDL-cholesterol levels selected from two well-defined epidemiological samples of U.S. non-Hispanic Whites (NHWs) and African Blacks followed by genotyping of identified variants in the entire datasets (623 NHWs, 788 African Blacks) and association analyses with major lipid traits. We identified a total of 40 sequence variants, of which 10 are novel. A total of 32 variants, including common tagSNPs (≥5% frequency) and all uncommon variants (<5% frequency) were successfully genotyped and considered for genotype-phenotype associations. Other than the established associations of APOE*2 and APOE*4 with LDL-cholesterol, we have identified additional independent associations with LDL-cholesterol. We have also identified multiple associations of uncommon and common APOE variants with HDL-cholesterol and triglycerides. Our comprehensive sequencing and genotype-phenotype analyses indicate that APOE genetic variation impacts HDL-cholesterol and triglycerides in addition to affecting LDL-cholesterol.

Impact of genetic variants in human scavenger receptor class B type I (SCARB1) on plasma lipid traits

BACKGROUND

Scavenger receptor class B type 1 (SCARB1) plays an important role in high-density lipoprotein cholesterol (HDL-C) metabolism in selective cholesteryl ester uptake and in free cholesterol cellular efflux.

METHODS AND RESULTS

This study aims to identify common (minor allele frequency ≥5%) and low-frequency/rare (minor allele frequency <5%) variants, using resequencing all 13 exons and exon-intron boundaries of SCARB1 in 95 individuals with extreme HDL-C levels selected from a population-based sample of 623 US non-Hispanic whites. The sequencing step identified 44 variants, of which 11 were novel with minor allele frequency <1%. Seventy-six variants (40 sequence variants, 32 common HapMap tag single nucleotide polymorphisms, and 4 relevant variants) were selected for genotyping in the total sample of 623 subjects followed by association analyses with lipid traits. Seven variants were nominally associated with apolipoprotein B (apoB; n=4) or HDL-C (n=3; P<0.05). Three variants associated with apoB remained significant after controlling false discovery rate. The most significant association was observed between rs4765615 and apoB (P=0.0059), while rs11057844 showed the strongest association with HDL-C (P=0.0035). A set of 17 rare variants (minor allele frequency ≤1%) showed significant association with apoB (P=0.0284). Haplotype analysis revealed 4 regions significantly associated with either apoB or HDL-C.

CONCLUSIONS

Our findings provide new information about the genetic role of SCARB1 in affecting plasma apoB levels in addition to its established role in HDL-C metabolism.

Lipoprotein lipase gene sequencing and plasma lipid profile

Lipoprotein lipase (LPL) plays a crucial role in lipid metabolism by hydrolyzing triglyceride (TG)-rich particles and affecting HDL cholesterol (HDL-C) levels. In this study, the entire LPL gene plus flanking regions were resequenced in individuals with extreme HDL-C/TG levels (n = 95), selected from a population-based sample of 623 US non-Hispanic White (NHW) individuals. A total of 176 sequencing variants were identified, including 28 novel variants. A subset of 64 variants [common tag single nucleotide polymorphisms (tagSNP) and selected rare variants] were genotyped in the total sample, followed by association analyses with major lipid traits. A gene-based association test including all genotyped variants revealed significant association with HDL-C (P = 0.024) and TG (P = 0.006). Our single-site analysis revealed seven independent signals (P < 0.05; r² < 0.40) with either HDL-C or TG. The most significant association was for the SNP rs295 exerting opposite effects on TG and HDL-C levels with P values of 7.5.10⁻⁴ and 0.002, respectively. Our work highlights some common variants and haplotypes in LPL with significant associations with lipid traits; however, the analysis of rare variants using burden tests and SKAT-O method revealed negligible effects on lipid traits. Comprehensive resequencing of LPL in larger samples is warranted to further test the role of rare variants in affecting plasma lipid levels.

Abstract 324: Resequencing of the Lipoprotein Lipase Gene in Individuals with Extreme HDL-Cholesterol Levels

In the United States, coronary heart disease (CHD) is the most common cause of death and number one killer of American males and females. Several risk factors for CHD have been identified, including low high-density lipoprotein cholesterol (HDL-C), elevated total cholesterol and low-density lipoprotein cholesterol (LDL-C), and high triglycerides (TG), but the genetic variations underlying inter-individual differences in plasma lipid profiles still need to be further explored. Thus, our interest has been focused on the lipoprotein lipase ( LPL ) which is one of the major genes involved in lipid metabolism. Previously, its sequence variation has been associated with the risk of CHD and other complex traits such as dyslipidemia, type 2 diabetes, essential hypertension, and Alzheimer’s disease. In order to further investigate the common and rare variation in LPL , we resequenced the 95 individuals having lower (n=48) and upper (n=48) 5 th percentiles of HDL-C levels selected from a well defined population-based non-Hispanic White (NHW) sample of 623 individuals. As a result of resequencing the entire LPL gene and additional 1kb at the 5’ and 3’ end, a total of 179 variants (substitutions or indels) plus one microsatellite were identified, including 91 relatively uncommon or rare variants [minor allele frequency (MAF) <0.05)] and 88 common variants (MAF ≥0.05). We have observed differences in the distribution of the 91 relatively uncommon or rare variants in two groups; 21 were present only in the low HDL-C group and 25 were present only in the high HDL-C group. The prevalence of uncommon or rare variants was found to be higher in the high HDL-C group. Thirty two common variants showed significantly different (p-value <0.05) allele frequency between the high and low HDL-C groups. To date, we have analyzed 39 common variants (MAF ≥0.05) in the entire sample and our results confirm the previous association findings as well as reveal novel variants that are associated with inter-individual variation in HDL and TG levels. The genotyping effort of relevant rare variants in the entire sample has also been completed and their analysis will help to unravel the role of rare variants in the LPL gene in regulation of plasma HDL-C levels and other lipids.