Small extrachromosomal circular DNA harboring targeted tumor suppressor gene mutations supports intratumor heterogeneity in mouse liver cancer induced by multiplexed CRISPR/Cas9.

BACKGROUND: Primary liver cancer has significant intratumor genetic heterogeneity (IGH), which drives cancer evolution and prevents effective cancer treatment. CRISPR/Cas9-induced mouse liver cancer models can be used to elucidate how IGH is developed. However, as CRISPR/Cas9 could induce chromothripsis and extrachromosomal DNA in cells in addition to targeted mutations, we wondered whether this effect contributes to the development of IGH in CRISPR/Cas9-induced mouse liver cancer. METHODS: CRISPR/Cas9-based targeted somatic multiplex-mutagenesis was used to target 34 tumor suppressor genes (TSGs) for induction of primary liver tumors in mice. Target site mutations in tumor cells were analyzed and compared between single-cell clones and their subclones, between different time points of cell proliferation, and between parental clones and single-cell clones derived from mouse subcutaneous allografts. Genomic instability and generation of extrachromosomal circular DNA (eccDNA) was explored as a potential mechanism underlying the oscillation of target site mutations in these liver tumor cells. RESULTS: After efficiently inducing autochthonous liver tumors in mice within 30-60 days, analyses of CRISPR/Cas9-induced tumors and single-cell clones derived from tumor nodules revealed multiplexed and heterogeneous mutations at target sites. Many target sites frequently displayed more than two types of allelic variations with varying frequencies in single-cell clones, indicating increased copy number of these target sites. The types and frequencies of targeted TSG mutations continued to change at some target sites between single-cell clones and their subclones. Even the proliferation of a subclone in cell culture and in mouse subcutaneous graft altered the types and frequencies of targeted TSG mutations in the absence of continuing CRISPR/Cas9 genome editing, indicating a new source outside primary chromosomes for the development of IGH in these liver tumors. Karyotyping of tumor cells revealed genomic instability in these cells manifested by high levels of micronuclei and chromosomal aberrations including chromosomal fragments and chromosomal breaks. Sequencing analysis further demonstrated the generation of eccDNA harboring targeted TSG mutations in these tumor cells. CONCLUSIONS: Small eccDNAs carrying TSG mutations may serve as an important source supporting intratumor heterogeneity and tumor evolution in mouse liver cancer induced by multiplexed CRISPR/Cas9.

[Analysis of flavonoids and antitumor activity of transgenic Saussurea involucrate].

The aim of this paper was to investigate the flavonoids of callus of transgenic and non-transgenic Saussurea involucrate and its antitumor activity on the esophageal cancer CaEs-17 cells. The species and content of mono-phenols were detected by high performance liquid chromatography. The growth of human esophageal cancer CaEs-17 cells was detected using CCK-8 assay, apoptosis morphology observation and flow cytometry. Expression of related apoptotic genes Bax and Bcl-2 were determined by qPCR. The results showed that the content of total flavonoids in the transgenic callus was 2.72 times that of the non-transgenic callus. The cyanidin-galactoside was detected in the transgenic callus, but not in the non-transgenic callus. The inhibitory effect of the extracts from the transgenic callus on CaEs-17 cells was more significant than that of the non-transgenic callus, and the IC50 value had a decreased of 26.4%. Flow cytometry analysis results showed that the apoptosis induction effect of the extracts from the transgenic callus on CaEs-17 cells was significantly better than that of the non-transgenic callus. Fluorescence quantitative PCR analysis results showed that the extracts from the transgenic callus could up-regulate the expression of proapoptotic gene Bax and down-regulate the expression of apoptotic gene Bcl-2, and the regulation effect of the transgenic callus was more significant. Therefore, compared with the non-transgenic callus, the antitumor activity of the extracts from the callus of transgenic S. involucrate on the esophageal cancer CaEs-17 cells was significantly increased, which was closely related to the accumulation of cyanidin-galactoside and its metabolism-related flavonoid compounds in the transgenic callus.

MicroRNA profiling and prediction of recurrence/relapse-free survival in stage I lung cancer.

About 30% stage I non-small cell lung cancer (NSCLC) patients undergoing resection will recur. Robust prognostic markers are required to better manage therapy options. MicroRNAs (miRNAs) are a class of small non-coding RNAs of 19-25 nt and play important roles in gene regulation in human cancers. The purpose of this study is to identify miRNA expression profiles that would better predict prognosis of stage I NSCLC. MiRNAs extracted from 527 stage I NSCLC patients were profiled on the human miRNA expression profiling v2 panel (Illumina). The expression profiles were analyzed for their association with cancer subtypes, lung cancer brain metastasis and recurrence/relapse free survival (RFS). MiRNA expression patterns between lung adenocarcinoma and squamous cell carcinoma differed significantly with 171 miRNAs, including Let-7 family members and miR-205. Ten miRNAs associated with brain metastasis were identified including miR-145*, which inhibit cell invasion and metastasis. Two miRNA signatures that are highly predictive of RFS were identified. The first contained 34 miRNAs derived from 357 stage I NSCLC patients independent of cancer subtype, whereas the second containing 27 miRNAs was adenocarcinoma specific. Both signatures were validated using formalin-fixed paraffin embedded and/or fresh frozen tissues in independent data set with 170 stage I patients. Our findings have important prognostic or therapeutic implications for the management of stage I lung cancer patients. The identified miRNAs hold great potential as targets for histology-specific treatment or prevention and treatment of recurrent disease.

Identification of Las2, a major modifier gene affecting the Pas1 mouse lung tumor susceptibility locus.

Lung cancer is the leading cause of cancer death worldwide. Here, we describe a genome-wide association study of chemically induced lung tumorigenesis on 593 mice from 21 inbred strains using 115,904 genotyped and 1,952,918 imputed single nucleotide polymorphisms (SNPs). Using a genetic background-controlled genome search, we identified a novel lung tumor susceptibility gene Las2 (Lung adenoma susceptibility 2) on distal chromosome 18. Las2 showed strong association with resistance to tumor induction (rs30245983; P = 1.87 x 10(-9)) as well as epistatic interactions (P = 1.71 x 10(-3)) with the pulmonary adenoma susceptibility 1 locus, a major locus affecting mouse lung tumor development (rs13459098, P = 5.64 x 10(-27)). Sequencing analysis revealed four nonsynonymous SNPs and two insertions/deletions in the susceptible allele of Las2, resulting in the loss of tumor suppressor activities in both cell colony formation and nude mouse tumorigenicity assays. Deletion of LAS2 was observed in approximately 40% of human lung adenocarcinomas, implying that loss of function of LAS2 may be a key step for lung tumorigenesis.

A gene expression signature predicts survival of patients with stage I non-small cell lung cancer.

BACKGROUND: Lung cancer is the leading cause of cancer-related death in the United States. Nearly 50% of patients with stages I and II non-small cell lung cancer (NSCLC) will die from recurrent disease despite surgical resection. No reliable clinical or molecular predictors are currently available for identifying those at high risk for developing recurrent disease. As a consequence, it is not possible to select those high-risk patients for more aggressive therapies and assign less aggressive treatments to patients at low risk for recurrence. METHODS AND FINDINGS: In this study, we applied a meta-analysis of datasets from seven different microarray studies on NSCLC for differentially expressed genes related to survival time (under 2 y and over 5 y). A consensus set of 4,905 genes from these studies was selected, and systematic bias adjustment in the datasets was performed by distance-weighted discrimination (DWD). We identified a gene expression signature consisting of 64 genes that is highly predictive of which stage I lung cancer patients may benefit from more aggressive therapy. Kaplan-Meier analysis of the overall survival of stage I NSCLC patients with the 64-gene expression signature demonstrated that the high- and low-risk groups are significantly different in their overall survival. Of the 64 genes, 11 are related to cancer metastasis (APC, CDH8, IL8RB, LY6D, PCDHGA12, DSP, NID, ENPP2, CCR2, CASP8, and CASP10) and eight are involved in apoptosis (CASP8, CASP10, PIK3R1, BCL2, SON, INHA, PSEN1, and BIK). CONCLUSIONS: Our results indicate that gene expression signatures from several datasets can be reconciled. The resulting signature is useful in predicting survival of stage I NSCLC and might be useful in informing treatment decisions.

Predictive factors for age at menopause in Caucasian females.

OBJECTIVE: Early onset of menopause results in the premature exposure to low estrogen levels and is associated with a number of postmenopausal health problems and higher risk of mortality. The aim of this study was to determine genetic and environmental factors associated with age at natural and surgical menopause. METHODS: Multiple regression analysis using a sample of Caucasians composed of 154 females with surgical and 248 with natural menopause. RESULTS: Breastfeeding is a significant predictor of earlier natural menopause (P<0.05). Use of oral contraceptives and smoking were not significantly associated with age at menopause. Females who did not have history of pregnancies are at significantly higher risk (P<0.001) of getting early surgical menopause than those who did. We also tested the association of seven single nucleotide polymorphisms (SNPs) of the estrogen receptor alpha (ER-alpha) gene with age at menopause. No association was observed with age at menopause but the PvuII p allele was overrepresented in women with surgical menopause and associated with menopause per se (P=0.029; OR=1.8, 95% CI=1.1-3.0). CONCLUSIONS: Breastfeeding and alcohol consumption are significantly associated with earlier natural menopause. No significant effects of the ER-alpha genotypes were observed on the age of menopause. Given the important role of the ER-alpha in estrogen signaling, which directly influences the menopausal process, further studies are required to better define the relationship between this gene and age at menopause.

Association analysis of estrogen receptor alpha gene polymorphisms with cross-sectional geometry of the femoral neck in Caucasian nuclear families.

Bone geometry is a key factor in bone strength, which is the ultimate intrinsic determinant of fracture risk. Though the heritability of bone geometry is high, little effort has been spent on searching for the underlying genes. In this study, employing a sample of 1,873 subjects from 405 Caucasian nuclear families, we studied seven single nucleotide polymorphisms (SNPs) and their haplotypes of the ER-alpha gene for association with six hip geometric variables, namely, cross-sectional area (CSA), cortical thickness (CT), endocortical diameter (ED), subperiosteal width (W), sectional modulus (Z) and buckling ratio (BR). The major method used was the quantitative transmission disequilibrium test (QTDT). Our major findings were summarized below. The within-family association between SNP4 (rs1801132) in exon 4 with endocortical diameter and subperiosteal width was detected in single locus analyses (P=0.008 and 0.021, respectively) and verified in haplotype analyses (P=0.034 and 0.058, respectively). The total association of SNP4 with these two diameters was also observed in both single locus and haplotype analyses (P=0.005 and 0.031 for ED, plus P=0.003 and 0.070 for W). In addition, the total association between SNP5 (rs932477) in intron 4 with cortical thickness and buckling ratio was detected (single locus analyses: P=0.035 and 0.041, respectively). Haplotype analyses further supported the above association (P=0.010 and 0.004, respectively). Similar patterns of associations with the studied SNPs and their haplotypes were present in subsamples stratified by sex, too. However, after permutation tests, the empirical significance level was set as P<0.011, which renders most associations insignificant. Therefore, we concluded that polymorphisms in the ER-alpha gene were nominally associated with femoral neck (FN) geometry variables estimated from DXA. Such genetic effects on hip geometry were not sex specific.

Contribution of genotype and ethnicity to bone mineral density variation in Caucasians and Chinese: a test for five candidate genes for bone mass.

BACKGROUND: Ethnicity is shown to be one of important factors affecting bone mineral density (BMD). The present study was performed to compare the association of six markers for five candidate genes with BMD variation in two populations of different ethnicity, Caucasian and Chinese, and the contribution of genotype and ethnicity to this variation in the populations. METHODS: The studied restriction fragment length polymorphisms were BsaH I of the calcium-sensing receptor gene, SacI of the alpha2HS-glycoprotein (AHSG) gene, PvuII and XbaI of the oestrogen receptor alpha gene, ApaI of the vitamin D receptor (VDR) gene and BstBI of the parathyroid hormone gene. The association of these markers with BMD was analysed by one-way and two-way ANOVA with adjustment for covariates. RESULTS: Two polymorphisms, AHSG-SacI and VDR-ApaI, showed no association with BMD, while the others were associated with BMD variation at some skeletal sites in either males or females. The polymorphisms indicated clear distinctions between the associations depending on ethnicity, gender and skeletal site. Similar patterns were observed in their contribution to the total population BMD variation. Ethnicity appears to have a larger effect on the total population BMD variation in females than in males. It may account, on the average, for about 2% total population BMD variation at the spine of females and about 1% at the hip of males and females. CONCLUSION: The results of the present study suggest that significant interethnic differentiation at some loci may contribute to the significant interethnic difference in BMD. However, this contribution apparently is not large.

Hotelling's T2 multivariate profiling for detecting differential expression in microarrays.

The most widely used statistical methods for finding differentially expressed genes (DEGs) are essentially univariate. In this study, we present a new T(2) statistic for analyzing microarray data. We implemented our method using a multiple forward search (MFS) algorithm that is designed for selecting a subset of feature vectors in high-dimensional microarray datasets. The proposed T2 statistic is a corollary to that originally developed for multivariate analyses and possesses two prominent statistical properties. First, our method takes into account multidimensional structure of microarray data. The utilization of the information hidden in gene interactions allows for finding genes whose differential expressions are not marginally detectable in univariate testing methods. Second, the statistic has a close relationship to discriminant analyses for classification of gene expression patterns. Our search algorithm sequentially maximizes gene expression difference/distance between two groups of genes. Including such a set of DEGs into initial feature variables may increase the power of classification rules. We validated our method by using a spike-in HGU95 dataset from Affymetrix. The utility of the new method was demonstrated by application to the analyses of gene expression patterns in human liver cancers and breast cancers. Extensive bioinformatics analyses and cross-validation of DEGs identified in the application datasets showed the significant advantages of our new algorithm.

Association tests of interleukin-6 (IL-6) and type II tumor necrosis factor receptor (TNFR2) genes with bone mineral density in Caucasians using a re-sampling approach.

Interleukin 6 (IL-6) and tumor necrosis factor (TNF) are important cytokines for bone turnover. In this study, a promoter C-174G single-nucleotide polymorphism (SNP) within the IL-6 gene affecting the transcription rate of IL-6 and an exon 6 T676G SNP of the TNF receptor 2 (TNFR2) gene causing an M196R amino-acid change were examined for their relationship with bone mineral density (BMD). Four hundred and five multi-offspring Caucasian families, including 389 male children and 744 female children, were used. One thousand re-samplings were conducted and in each data set, one child was randomly chosen from each family. For each data set, one-way analysis of variance (ANOVA) test was independently implemented using age, age2, sex, height and weight as covariates. There were 523, 288, 204 and 369 significant results out of 1,000-replicate re-samplings of the data of the IL-6 SNP (P<0.05) for one-third, mid-distal, ultradistal radius BMD, and the first principal component (PC1) extracted from the three radial BMDs, respectively, which means that the confidences for associations of the C-174G SNP in the IL-6 gene with one-third, mid-distal, ultradistal radius (totally called distal forearm) BMDs, and PC1, were 52.3, 28.8, 20.4 and 36.9%, respectively. For this SNP with BMD at other skeletal sites and the TNFR2 T676G SNP with BMD at any site, significant results were far less than 200 times out of 1,000 re-sampling replicates. The exceedingly consistent permutation results further improved the confidence of the associations. It may imply that the IL-6 C-174G SNP is associated with distal forearm BMD, but there is no evidence that the TNFR2 T676G SNP is related with BMD in US Caucasians. This is the first attempt to conduct association test utilizing a re-sampling approach. Our results may be more informative than other association analyses that were only based on one sampling result. The results also suggest that different samplings could produce significantly diverse results even for the same population and the results from one sampling are unlikely to be conclusive. Our results have significant implications for association studies and interpretation of non-reproducible association findings.

Test of linkage and/or association between the estrogen receptor alpha gene with bone mineral density in Caucasian nuclear families.

Extensive studies have been performed on the association between the estrogen receptor alpha (ER-alpha) gene and bone mineral density (BMD). Despite considerable efforts, the studies using limited markers and relatively small sample size have yielded largely inconsistent results. In this study, 1873 Caucasian subjects from 405 nuclear families containing 1512 sib pairs were recruited. BMD at the lumbar spine (LS) and femoral neck (FN) was measured by dual-energy X-ray absorptiometry (DXA). Seven single-nucleotide polymorphisms (SNPs) spanning from exon 1 to 8 in the ER-alpha gene were genotyped. The program QTDT (quantitative transmission disequilibrium test) was applied to test linkage and/or association of the ER-alpha gene and BMD variation using individual SNP markers and reconstructed haplotypes. Linkage disequilibrium (LD) was generally detected for SNPs in the ER-a gene (P < 0.05). Associations were observed between SNP rs932477 and FN BMD (P = 0.028), and between the most predominant three-marker haplotype (GCG) containing SNP rs932477 and FN BMD (P = 0.010). Within-family association (present only with both linkage and association) between SNP rs2228480 (G2014A) and FN BMD (P = 0.015) was observed. The most predominant seven-SNP haplotype (TCGCGGG) was associated with higher LS BMD (P = 0.015). However, after correction for multiple testing, these associations did not reach statistical significance. Denser markers may be necessary to better define the relationship between the ER-alpha gene and BMD variation in our sample.

Patterns of linkage disequilibrium and haplotype distribution in disease candidate genes.

BACKGROUND: The adequacy of association studies for complex diseases depends critically on the existence of linkage disequilibrium (LD) between functional alleles and surrounding SNP markers. RESULTS: We examined the patterns of LD and haplotype distribution in eight candidate genes for osteoporosis and/or obesity using 31 SNPs in 1,873 subjects. These eight genes are apolipoprotein E (APOE), type I collagen alpha1 (COL1A1), estrogen receptor-alpha (ER-alpha), leptin receptor (LEPR), parathyroid hormone (PTH)/PTH-related peptide receptor type 1 (PTHR1), transforming growth factor-beta1 (TGF-beta1), uncoupling protein 3 (UCP3), and vitamin D (1,25-dihydroxyvitamin D3) receptor (VDR). Yin yang haplotypes, two high-frequency haplotypes composed of completely mismatching SNP alleles, were examined. To quantify LD patterns, two common measures of LD, D' and r2, were calculated for the SNPs within the genes. The haplotype distribution varied in the different genes. Yin yang haplotypes were observed only in PTHR1 and UCP3. D' ranged from 0.020 to 1.000 with the average of 0.475, whereas the average r2 was 0.158 (ranging from 0.000 to 0.883). A decay of LD was observed as the intermarker distance increased, however, there was a great difference in LD characteristics of different genes or even in different regions within gene. CONCLUSION: The differences in haplotype distributions and LD patterns among the genes underscore the importance of characterizing genomic regions of interest prior to association studies.

Estrogen receptor alpha and vitamin D receptor gene polymorphisms and bone mineral density: association study of healthy pre- and postmenopausal Chinese women.

In the present study, we tested the association between the estrogen receptor alpha (ER-alpha) and vitamin D receptor (VDR) genes with bone mineral density (BMD). A total of 649 healthy Chinese women, classified as pre-menopausal (N=388) and post-menopausal (N=261) groups, were genotyped at the ER-alpha PvuII, XbaI, and VDR ApaI sites. BMDs at the lumbar spine (L(1)-L(4)) and total hip were measured by dual-energy X-ray absorptiometry. For the VDR ApaI locus, AA carriers had lower spine BMD than Aa (p=0.02) and aa carriers (p<0.01) in the pre-menopausal group. For the ER-alpha gene, carriers of haplotype px had lower spine BMD than the non-carriers (p=0.03) in the pre-menopausal group. Furthermore, we observed significant interaction between the ER-alpha and VDR genes in the post-menopausal group: with AA genotype (or A allele) at the VDR ApaI locus, pX carriers had higher spine BMD than the non-carriers (p=0.02), and PX carriers had lower hip BMD than the non-carriers (p=0.04). Our data suggest that the ER-alpha and VDR genes may be associated with the BMD variation in Chinese women.