Integrative analysis identifies potential DNA methylation biomarkers for pan-cancer diagnosis and prognosis

DNA methylation status is closely associated with diverse diseases, and is generally more stable than gene expression, thus abnormal DNA methylation could be important biomarkers for tumor diagnosis, treatment and prognosis. However, the signatures regarding DNA methylation changes for pan-cancer diagnosis and prognosis are less explored. Here we systematically analyzed the genome-wide DNA methylation patterns in diverse TCGA cancers with machine learning. We identified seven CpG sites that could effectively discriminate tumor samples from adjacent normal tissue samples for 12 main cancers of TCGA (1216 samples, AUC > 0.99). Those seven potential diagnostic biomarkers were further validated in the other 9 different TCGA cancers and 4 independent datasets (AUC > 0.92). Three out of the seven CpG sites were correlated with cell division, DNA replication and cell cycle. We also identified 12 CpG sites that can effectively distinguish 26 different cancers (7605 samples), and the result was repeatable in independent datasets as well as two disparate tumors with metastases (micro-average AUC > 0.89). Furthermore, a series of potential signatures that could significantly predict the prognosis of tumor patients for 7 different cancer were identified via survival analysis (p-value < 1e-4). Collectively, DNA methylation patterns vary greatly between tumor and adjacent normal tissues, as well as among different types of cancers. Our identified signatures may aid the decision of clinical diagnosis and prognosis for pan-cancer and the potential cancer-specific biomarkers could be used to predict the primary site of metastatic breast and prostate cancers.     

The Relationship of Ethnicity,Socioeconomic Factors,and Overweight in U.S. Adolescents

Objective: To examine the extent to which race/ethnic differences in income and education account for sex‐specific disparities in overweight prevalence in white, African American, Hispanic, and Asian U.S. teens. Research Methods and Procedures: We used nationally representative data collected from 13,113 U.S. adolescents enrolled in the National Longitudinal Study of Adolescent Health. Logistic regression models were used to examine the relationship of family income and parental education to overweight prevalence (body mass index ≥ 85th percentile of age and sex‐specific cutoff points from the 2000 Centers for Disease Control and Prevention/National Center for Health Statistics growth charts). In addition, we used coefficients from our logistic regression models to project the effects on overweight prevalence of equalizing the socioeconomic status (SES) differences between race/ethnic groups. Results: Keeping adolescents in their same environments and changing only family income and parental education had a limited effect on the disparities in overweight prevalence. Ethnicity‐SES‐overweight differences were greater among females than males. Given that overweight prevalence decreased with increasing SES among white females and remained elevated and even increased among higher SES African‐American females, African‐American/white disparity in overweight prevalence increased at the highest SES. Conversely, disparity was lessened at the highest SES for white, Hispanic, and Asian females. Among males, disparity was lowest at the average SES level. Discussion: One cannot automatically assume that the benefits of increased SES found among white adults will transfer to other gender‐age‐ethnic groups. Our findings suggest that efforts to reduce overweight disparities between ethnic groups must look beyond income and education and focus on other factors, such as environmental, contextual, biological, and sociocultural factors.     

The Aminoacyl-tRNA Synthetase and tRNA Expression Levels Are Deregulated in Cancer and Correlate Independently with Patient Survival

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes that load amino acids to their cognate tRNA molecules. The expression of certain ARSs and tRNAs has been shown to be deregulated in cancer, presumably to accommodate elevated protein synthesis requirements. In this work, the expression of cytoplasmic ARSs and tRNAs in ten TCGA cancers has been systematically examined. ARSs were found to be mostly upregulated in tumours and their upregulation often correlated with worse patient survival. tRNAs were found to be either upregulated or downregulated in tumours and their expression sometimes correlated to worse survival outcomes. However, although the expression of most ARSs and tRNAs was deregulated in tumours when compared to healthy adjacent tissues, only in a few cases, and independently, did it correlate to patient survival. These data point to the general uncoupling of concomitant ARS and tRNA expression deregulation and patient survival, highlighting the different ARS and tRNA requirements in cancers.     

Comparing genetic ancestry and self-reported race/ethnicity in a multiethnic population in New York City

Self-reported race/ethnicity is frequently used in epidemiological studies to assess an individual’s background origin. However, in admixed populations such as Hispanic, self-reported race/ethnicity may not accurately represent them genetically because they are admixed with European, African and Native American ancestry. We estimated the proportions of genetic admixture in an ethnically diverse population of 396 mothers and 188 of their children with 35 ancestry informative markers (AIMs) using the STRUCTURE version 2.2 program. The majority of the markers showed significant deviation from Hardy-Weinberg equilibrium in our study population. In mothers self-identified as Black and White, the imputed ancestry proportions were 77.6% African and 75.1% European respectively, while the racial composition among self-identified Hispanics was 29.2% European, 26.0% African, and 44.8% Native American. We also investigated the utility of AIMs by showing the improved fitness of models in paraoxanase-1 genotype-phenotype associations after incorporating AIMs; however, the improvement was moderate at best. In summary, a minimal set of 35 AIMs is sufficient to detect population stratification and estimate the proportion of individual genetic admixture; however, the utility of these markers remains questionable.     

Two distinct phenotypes of immunologically hot gastric cancer subtypes

An in-depth understanding of the tumor microenvironment (TME) is required for the development of improved combination immunotherapies for gastric cancer. Recently, we classified these cancers into four main types defined by their immunological attributes, namely Hot 1, Hot 2, Intermediate and Cold. Of these, the T cell-inflamed “Hot” tumors were further divided into Hot 1 and Hot 2 with different clinical outcomes. Thus, overall survival and progression-free survival of patients with Hot 1 tumors were shorter than with Hot 2. In the present study, we re-evaluated RNA-Seq data of 6 Hot 1 and 6 Hot 2 gastric cancers to elucidate the underlying reason for the poor prognosis and T cell dysfunction in the former. In addition, 56 Hot 1 and 27 Hot 2 tumors in The Cancer Genome Atlas (TCGA) were analyzed. We report that single sample Gene Set Enrichment Analysis (ssGSEA) and differential gene expression analysis identified differences between Hot 1 and Hot 2 tumors involved in metabolism and cell adhesion pathways. Therefore, it is suggested that strategies to modulate active metabolism in Hot 1 tumors should be integrated into the treatment of these gastric cancers.     

JOINED AT THE HIP: LINKED CHARACTERS AND THE PROBLEM OF MISSING DATA IN STUDIES OF DISPARITY

Paleontological investigations into morphological diversity, or disparity, are often confronted with large amounts of missing data. We illustrate how missing discrete data affect disparity using a novel simulation for removing data based on parameters from published datasets that contain both extinct and extant taxa. We develop an algorithm that assesses the distribution of missing characters in extinct taxa, and simulates data loss by applying that distribution to extant taxa. We term this technique “linkage.” We compare differences in disparity metrics and ordination spaces produced by linkage and random character removal. When we incorporated linkage among characters, disparity metrics declined and ordination spaces shrank at a slower rate with increasing missing data, indicating that correlations among characters govern the sensitivity of disparity analysis. We also present and test a new disparity method that uses the linkage algorithm to correct for the bias caused by missing data. We equalized proportions of missing data among time bins before calculating disparity, and found that estimates of disparity changed when missing data were taken into account. By removing the bias of missing data, we can gain new insights into the morphological evolution of organisms and highlight the detrimental effects of missing data on disparity analysis.     

Ethnic differences in allelic distribution of IFN-g in South African women but no link with cervical cancer

BACKGROUND: The failure of specific types of human papillomaviruses (HPV) to raise effective immune responses may be important in the pathogenesis of cervical cancer, the second most common cancer in South African women. Polymorphisms of a number of cytokine genes have been implicated in inducing susceptibility or resistance to cancers caused by infectious agents owing to their role in determining host immune response. Polymorphisms of IL-10 and IFN-gamma genes are believed to influence the expression and/or secretion levels of their respective cytokines. METHODS AND RESULTS: In this study, women with histologically proven cancer of the cervix (n = 458) and hospital-based controls (n = 587) were investigated for bi-allelic -1082 (A/G) polymorphisms of IL-10 and the bi-allelic +874(A/T) polymorphisms of IFN-gamma. In addition, the distributions of the allelic frequencies were stratified in both the African and mixed race population groups of South Africa. We found striking differences in the allele distribution of IFN-gamma (X2 = 0.02) among the two ethnic groups. A significant increase in the allele distribution of the IFN-gamma AA genotype was found in the African group compared to the mixed population group (OR, 0.5; 95% CI, 0.2-1.0). For IL-10 there were no significant allelic differences between the two South African ethnic groups. Furthermore, when the ethnic groups were combined the IL-10 allelic frequencies in the combined South African data were similar to those observed in an Oriental population from Southern China and in an Italian population. However, the allele frequencies of the IFN-gamma genotype among the two South African ethnic groups were different when compared to an Italian Caucasoid group. While crude analysis of these data showed both statistically significantly increased and diminished risks of cervical cancer among high producers of INF-gamma and low producers of IL-10 respectively, these associations were no longer significant when the data were adjusted for confounding factors. CONCLUSION: These findings demonstrate a clear correlation between ethnicity and IFN-gamma polymorphism across different population groups. However, these differences in ethnicity and gene polymorphisms in the aforementioned cytokines are suggested not to influence the development of invasive cervical cancer but may represent an important susceptibility biomarker for other diseases and should be explored further.     

Genome-wide local ancestry approach identifies genes and variants associated with chemotherapeutic susceptibility in African Americans

Chemotherapeutic agents are used in the treatment of many cancers, yet variable resistance and toxicities among individuals limit successful outcomes. Several studies have indicated outcome differences associated with ancestry among patients with various cancer types. Using both traditional SNP-based and newly developed gene-based genome-wide approaches, we investigated the genetics of chemotherapeutic susceptibility in lymphoblastoid cell lines derived from 83 African Americans, a population for which there is a disparity in the number of genome-wide studies performed. To account for population structure in this admixed population, we incorporated local ancestry information into our association model. We tested over 2 million SNPs and identified 325, 176, 240, and 190 SNPs that were suggestively associated with cytarabine-, 5'-deoxyfluorouridine (5'-DFUR)-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-4)). Importantly, some of these variants are found only in populations of African descent. We also show that cisplatin-susceptibility SNPs are enriched for carboplatin-susceptibility SNPs. Using a gene-based genome-wide association approach, we identified 26, 11, 20, and 41 suggestive candidate genes for association with cytarabine-, 5'-DFUR-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-3)). Fourteen of these genes showed evidence of association with their respective chemotherapeutic phenotypes in the Yoruba from Ibadan, Nigeria (p<0.05), including TP53I11, COPS5 and GAS8, which are known to be involved in tumorigenesis. Although our results require further study, we have identified variants and genes associated with chemotherapeutic susceptibility in African Americans by using an approach that incorporates local ancestry information.     

MEXPRESS: visualizing expression,DNA methylation and clinical TCGA data

Background

In recent years, increasing amounts of genomic and clinical cancer data have become publically available through large-scale collaborative projects such as The Cancer Genome Atlas (TCGA). However, as long as these datasets are difficult to access and interpret, they are essentially useless for a major part of the research community and their scientific potential will not be fully realized. To address these issues we developed MEXPRESS, a straightforward and easy-to-use web tool for the integration and visualization of the expression, DNA methylation and clinical TCGA data on a single-gene level (http://mexpress.be).

Results

In comparison to existing tools, MEXPRESS allows researchers to quickly visualize and interpret the different TCGA datasets and their relationships for a single gene, as demonstrated for GSTP1 in prostate adenocarcinoma. We also used MEXPRESS to reveal the differences in the DNA methylation status of the PAM50 marker gene MLPH between the breast cancer subtypes and how these differences were linked to the expression of MPLH.

Conclusions

We have created a user-friendly tool for the visualization and interpretation of TCGA data, offering clinical researchers a simple way to evaluate the TCGA data for their genes or candidate biomarkers of interest.

Electronic supplementary material

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

Genetic admixture, self-reported ethnicity, self-estimated admixture, and skin pigmentation among Hispanics and Native Americans

The relationship between ethnicity and biology is of interest to anthropologists, biomedical scientists, and historians in understanding how human groups are constructed. Ethnic self-identification in recently admixed groups such as Hispanics, African Americans, and Native Americans (NA) is likely to be complex due to the heterogeneity in individual admixture proportions and social environments within these groups. This study examines the relationships between self-identified ethnicity, self-estimated admixture proportions, skin pigmentation, and genetic marker estimated admixture proportions. These measures were assessed using questionnaires, skin color measurements, and genotyping of a panel of 76 ancestry informative markers, among 170 Hispanics and NAs from New Mexico, a state known for its complex history of interactions between people of NA and European (EU) ancestry. Results reveal that NAs underestimate their degree of EU admixture, and that Hispanics underestimate their degree of NA admixture. Within Hispanics, genetic-marker estimated admixture is better predicted by forehead skin pigmentation than by self-estimated admixture. We also find that Hispanic individuals self-identified as "half-White, half Hispanic" and "Spanish" have lower levels of NA admixture than those self-identified as "Mexican" and "Mexican American." Such results highlight the interplay between culture and biology in how individuals identify and view themselves, and have implications for how ethnicity and disease risk are assessed in a medical setting.     

WWOX at the crossroads of cancer,metabolic syndrome related traits and CNS pathologies

WWOX was cloned as a putative tumor suppressor gene mapping to chromosomal fragile site FRA16D. Deletions affecting WWOX accompanied by loss of expression are frequent in various epithelial cancers. Translocations and deletions affecting WWOX are also common in multiple myeloma and are associated with worse prognosis. Metanalysis of gene expression datasets demonstrates that low WWOX expression is significantly associated with shorter relapse-free survival in ovarian and breast cancer patients. Although somatic mutations affecting WWOX are not frequent, analysis of TCGA tumor datasets led to identifying 44 novel mutations in various tumor types. The highest frequencies of mutations were found in head and neck cancers and uterine and gastric adenocarcinomas.     

A yeast replication origin consists of multiple copies of a small conserved sequence

Autonomously replicating sequences (ARSs) of the yeast S. cerevisiae function as replication origins on plasmids and probably also on chromosomes. ARS function requires a copy of the ARS core consensus (5'-[A/T]TTTAT[A/G]TTT[A/T]-3') and additional sequences 3' to the T-rich strand of the consensus. Our analysis of an ARS from chromosome III, the C2G1 ARS, suggests that ARS function depends on the presence of an exact match to the core consensus and the presence of additional near matches in the 3' flanking region. We have demonstrated that ARS function can be mediated by multiple matches to the core consensus by constructing synthetic ARS elements from oligonucleotides containing copies of the consensus sequence. We find that two copies of the core consensus are sufficient for ARS activity and that an artificial ARS as efficient as a natural chromosomal ARS can be constructed from multiple core consensus elements in a specific orientation.     

Racial admixture and its impact on BMI and blood pressure in African and Mexican Americans

Admixed populations such as African Americans and Hispanic Americans present both challenges and opportunities in genetic epidemiologic research. Because of variation in admixture levels among individuals, case-control association studies may be subject to stratification bias. On the other hand, admixed populations also present special opportunities both for examining the role of genetic and environmental factors for observed racial/ethnic differences, and for possibly mapping alleles that contribute to such differences. Here we examined the distribution and relationship of individual admixture (IA) estimates with BMI and three measures of blood pressure in two admixed populations in the NHLBI Family Blood Pressure Program (FBPP): African Americans and Mexican Americans. For the African Americans, we observed modest but significant differences in average African IA among four recruitment sites. We observed a slight excess of African IA among hypertensives compared to normotensives, and a positive (non-significant) regression of African IA on blood pressure in untreated participants. Within Mexican Americans, we found no difference in average IA between hypertensives and normotensives, but a positive (marginally significant) regression of African IA on diastolic blood pressure. We also observed a significant positive regression of Caucasian IA (and negative regression of Native American IA) on BMI. Our results are suggestive of genetic differences between Africans and non-Africans that influence blood pressure, but such effects are likely to be modest compared to environmental ones. Excess obesity among Native Americans compared to whites is not consistent with a simple genetic explanation.     

APOA5 polymorphisms influence plasma triglycerides in young, healthy African Americans and whites of the CARDIA Study

Genetic variation in the apolipoprotein A-V gene (APOA5) has been associated with variation in plasma triglyceride (TG) levels in African American and white females and males older than 40 years and/or at increased risk of coronary artery disease. We have examined whether plasma TG levels are associated with 16 APOA5 polymorphisms in young (18-30 years) African American (1,075 females and 783 males) and white (1,041 females and 932 males) individuals of the Coronary Artery Risk Development in Young Adults (CARDIA) Study selected without regard to health. Plasma TG was significantly (P < 0.01) associated with markers 27376 and 28837 (-3A/G) in both white females and males, with 27709 (-1131T/C) and 29085 in white males, with 29009 (S19W) in African American females and white males, and with 30966 in African American females. No statistically significant associations were observed in African American males. These six single-nucleotide polymorphisms individually accounted for 0-0.78% of lnTG variation among white females, 0-2.46% among white males, and 0-0.69% among African American females. The results of our study suggest a small but replicable context-dependent influence of the APOA5 gene region on plasma TG levels in young, healthy individuals.     

A comprehensive analysis of IDO1 expression with tumour‐infiltrating immune cells and mutation burden in gynaecologic and breast cancers

Gynaecologic and breast cancers share some similarities at the molecular level. The aims of our study are to highlight the similarities and differences about IDO1, an important immune‐related gene in female cancers. The NGS data from TCGA of cervical squamous cell carcinoma (CESC), ovarian serous cystadenocarcinoma (OV), uterine corpus endometrial carcinoma (UCEC), uterine carcinosarcoma (UCS) and breast invasive carcinoma (BRCA) were analysed to identify molecular features, and clinically significant and potential therapeutic targets of IDO1. We found IDO1 was significantly up‐regulated in four gynaecologic cancers and breast cancer. According to breast cancer PAM50 classification scheme, IDO1 expression was higher in tumours of basal than other subtypes and showed better survival prognosis in BRCA and OV. Through immune infiltration analysis, we found a strong correlation between IDO1 and immune cell populations especially for dendritic cells and T cells. In addition, we investigated the association between IDO1 and tumour mutation burden (TMB) and found that IDO1 was significantly correlated with TMB in BRCA and CESC. GSVA revealed that hallmarks significantly correlated with IDO1 were involved in interferon gamma response, allograft rejection and inflammatory response. We also found PD‐L1 and LAG3 were highly positive related to IDO1 in gynaecologic cancers when comparing with their corresponding normal tissues. Our results indicated that IDO1 participated in anti‐tumour immune process and is correlated with mutation burden. These findings may expand our outlook of potential anti‐IDO1 treatments.     

Examination of ancestry and ethnic affiliation using highly informative diallelic DNA markers: application to diverse and admixed populations and implications for clinical epidemiology and forensic medicine

We and others have identified several hundred ancestry informative markers (AIMs) with large allele frequency differences between different major ancestral groups. For this study, a panel of 199 widely distributed AIMs was used to examine a diverse set of 796 DNA samples including self-identified European Americans, West Africans, East Asians, Amerindians, African Americans, Mexicans, Mexican Americans, Puerto Ricans and South Asians. Analysis using a Bayesian clustering algorithm (STRUCTURE) showed grouping of individuals with similar ethnic identity without any identifier other than the AIMs genotyping and showed admixture proportions that clearly distinguished different individuals of mixed ancestry. Additional analyses showed that, for the majority of samples, the predicted ethnic identity corresponded with the self-identified ethnicity at high probability (P > 0.99). Overall, the study demonstrates that AIMs can provide a useful adjunct to forensic medicine, pharmacogenomics and disease studies in which major ancestry or ethnic affiliation might be linked to specific outcomes.Electronic Supplementary Material Supplementary material is available for this article at