Purpose: The intratumoral heterogeneity (ITH) and the evolution of genomic architectures associated with the development of distant metastases are not well understood in colorectal cancers.Experimental Design: We performed multiregion biopsies of primary and liver metastatic regions from five colorectal cancers with whole-exome sequencing and copy number profiling.Results: In addition to a substantial level of genetic ITH, multiregion genetic profiling identifies the subclonal mutational architecture, leading to the region-based or spatial categorization of somatic mutations and the inference of intratumoral evolutionary history of cancers. The universal mutations (those observed in all the regional biopsies) are enriched in known cancer genes such as APC and TP53 with distinct mutational spectra compared with biopsy-or region-specific mutations, suggesting that major operative mutational mechanisms and their selective pressures are not constant across the metastatic progression. The phylogenies inferred from genomic data show branching evolutionary patterns where some primary biopsies are often segregated with metastastic lesions. Our analyses also revealed that copy number changes such as the chromosomal gains of c-MYC and chromothripsis can be region specific and the potential source of genetic ITH.Conclusions: Our data show that the genetic ITH is prevalent in colorectal cancer serving as a potential driving force to generate metastasis-initiating clones and also as a means to infer the intratumoral evolutionary history of cancers. The paucity of recurrent metastasis-clonal events suggests that colorectal cancer distant metastases may not follow a uniform course of genomic evolution, which should be considered in the genetic diagnosis and the selection of therapeutic targets for the advanced colorectal cancer.

Objective. To evaluate associations of genetic polymorphisms in cytochrome P450 (CYP) isoforms 2D6, 3A5, and 3A4 with blood concentrations of hydroxychloroquine (HCQ) and its metabolite, N-desethyl HCQ (DHCQ), in patients with systemic lupus erythematosus (SLE).Methods. SLE patients taking HCQ for >3 months were recruited and were genotyped for 4 single-nucleotide polymorphisms in CYP2D6*10, CYP3A5*3, and CYP3A4* 18B. Blood HCQ and DHCQ concentrations ([HCQ] and [DHCQ]) were measured and their association with corresponding genotypes was investigated.Results. A total of 194 patients were included in the analysis. CYP2D6*10 polymorphisms (rs1065852 and rs1135840) were significantly associated with the [DHCQ]:[HCQ] ratio after adjustment for age, sex, dose per weight per day, and SLE Disease Activity Index score (P 5 0.03 and P < 0.01, respectively). In adjusted models, the [DHCQ]:[HCQ] ratio was highest in patients with the G/G genotype of the CYP2D6*10 (rs1065852) polymorphism and lowest in those with the A/A genotype (P 5 0.03). Similarly, the [DHCQ]:[HCQ] ratio was highest in patients with the C/C genotype of the CYP2D6*10 (rs1135840) polymorphism and lowest in those with the G/G genotype (P < 0.01). The CYP2D6*10 (rs1065852) polymorphism was significantly related to the [DHCQ] (P 5 0.01). However, the polymorphisms of CYP3A5*3 and CYP3A4*18B did not show any significant association with the Hydroxychloroquine (HCQ) is an antimalarial drug that is proven to be a safe and effective treatment for systemic lupus erythematosus (SLE) (1). Despite its wide application, only a few studies have previously measured blood HCQ levels in patients taking the drug in the long term. Interestingly, blood HCQ levels vary widely between patients, even those taking the same dose at the same frequency (2-7). This interpersonal variation is not well understood. However, the blood HCQ concentration is closely related to the treatment response in autoimmune diseases such as SLE (4,5,8,9). Therefore, identifying factors related to variations in blood levels is critical to maximize the benefit of HCQ in SLE patients.Cytochrome P450 (CYP) enzymes play major roles in drug metabolism. Certain single-nucleotide polymorphisms (SNPs) in CYP genes may have a large impact on CYP enzyme activity. HCQ is metabolized to N-desethyl HCQ (DHCQ) in the liver through the Ndesethylation pathway (10,11). This reaction is mediated by CYP 2D6, 3A4, 3A5, and 2C8 isoforms (11-13).

Although ductal carcinoma in situ (DCIS) precedes invasive ductal carcinoma (IDC), the related genomic alterations remain unknown. To identify the genomic landscape of DCIS and better understand the mechanisms behind progression to IDC, we performed whole-exome sequencing and copy number profiling for six cases of pure DCIS and five pairs of synchronous DCIS and IDC. Pure DCIS harbored well-known mutations (e.g., TP53, PIK3CA and AKT1), copy number alterations (CNAs) and chromothripses, but had significantly fewer driver genes and co-occurrence of mutation/CNAs than synchronous DCIS-IDC. We found neither recurrent nor significantly mutated genes with synchronous DCIS-IDC compared to pure DCIS, indicating that there may not be a single determinant for pure DCIS progression to IDC. Of note, synchronous DCIS genomes were closer to IDC than pure DCIS. Among the clinicopathologic parameters, progesterone receptor (PR)-negative status was associated with increased mutations, CNAs, co-occurrence of mutations/CNAs and driver mutations. Our results indicate that although pure DCIS has already acquired some drivers, more changes are needed to progress to IDC. In addition, IDC-associated DCIS is more aggressive than pure DCIS at genomic level and should really be considered IDC. Finally, the data suggest that PR-negativity could be used to predict aggressive breast cancer genotypes.

BackgroundHan Chinese, Japanese and Korean, the three major ethnic groups of East Asia, share many similarities in appearance, language and culture etc., but their genetic relationships, divergence times and subsequent genetic exchanges have not been well studied.ResultsWe conducted a genome-wide study and evaluated the population structure of 182 Han Chinese, 90 Japanese and 100 Korean individuals, together with the data of 630 individuals representing 8 populations wordwide. Our analyses revealed that Han Chinese, Japanese and Korean populations have distinct genetic makeup and can be well distinguished based on either the genome wide data or a panel of ancestry informative markers (AIMs). Their genetic structure corresponds well to their geographical distributions, indicating geographical isolation played a critical role in driving population differentiation in East Asia. The most recent common ancestor of the three populations was dated back to 3000 ~ 3600 years ago. Our analyses also revealed substantial admixture within the three populations which occurred subsequent to initial splits, and distinct gene introgression from surrounding populations, of which northern ancestral component is dominant.ConclusionsThese estimations and findings facilitate to understanding population history and mechanism of human genetic diversity in East Asia, and have implications for both evolutionary and medical studies.Electronic supplementary materialThe online version of this article (10.1186/s41065-018-0057-5) contains supplementary material, which is available to authorized users.

To elucidate the pathogenesis of hepatocellular carcinoma (HCC) and develop useful prognosis predictors, it is necessary to identify biologically relevant genomic alterations in HCC. In our study, we defined recurrently altered regions (RARs) common to many cases of HCCs, which may contain tumor-related genes, using whole-genome array-CGH and explored their associations with the clinicopathologic features. Gene set enrichment analysis was performed to investigate functional implication of RARs. On an average, 23.1% of the total probes were altered per case. Mean numbers of altered probes are significantly higher in high-grade, bigger and microvascular invasion (MVI) positive tumors. In total, 32 RARs (14 gains and 18 losses) were defined and 4 most frequent RARs are gains in 1q21.1-q32.1 (64.5%), 1q32.1-q44 (59.2%), 8q11.21-q24.3 (48.7%) and a loss in 17p13.3-p12 (51.3%). Through focusing on RARs, we identified genes and functional pathways likely to be involved in hepatocarcinogenesis. Among genes in the recurrently gained regions on 1q, expression of KIF14 and TPM3 was significantly increased, suggesting their oncogenic potential in HCC. Some RARs showed the significant associations with the clinical features. Especially, the recurrent loss in 9p24.2-p21.1 and gain in 8q11.21-q24.3 are associated with the high tumor grade and MVI, respectively. Functional analysis showed that cytokine receptor binding and defense response to virus pathways are significantly enriched in high grade-related RARs. Taken together, our results and the strategy of analysis will help to elucidate pathogenesis of HCC and to develop biomarkers for predicting behaviors of HCC.

Recent discovery of the copy number variation (CNV) in normal individuals has widened our understanding of genomic variation. However, most of the reported CNVs have been identified in Caucasians, which may not be directly applicable to people of different ethnicities. To profile CNV in East-Asian population, we screened CNVs in 3578 healthy, unrelated Korean individuals, using the Affymetrix Genome-Wide Human SNP array 5.0. We identified 144 207 CNVs using a pooled data set of 100 randomly chosen Korean females as a reference. The average number of CNVs per genome was 40.3, which is higher than that of CNVs previously reported using lower resolution platforms. The median size of CNVs was 18.9 kb (range 0.2–5406 kb). Copy number losses were 4.7 times more frequent than copy number gains. CNV regions (CNVRs) were defined by merging overlapping CNVs identified in two or more samples. In total, 4003 CNVRs were defined encompassing 241.9 Mb accounting for ∼8% of the human genome. A total of 2077 CNVRs (51.9%) were potentially novel. Known CNVRs were larger and more frequent than novel CNVRs. Sixteen percent of the CNVRs were observed in ≥1% of study subjects and 24% overlapped with the OMIM genes. A total of 476 (11.9%) CNVRs were associated with segmental duplications. CNVS/CNVRs identified in this study will be valuable resources for studying human genome diversity and its association with disease.

Mutator phenotype tumors provide unique opportunities to unravel malignant progression because of various gene alterations acquired during clonal tumor evolution. Gastric carcinomas, which have been known to show frequent genetic instability, would be composed of initial gene alterations shared by most tumor areas and subsequent alterations restricted to particular tumor sites. To analyse the timing of genetic events, we examined separate sites of tumor tissue obtained from a given gastric carcinoma patient with microsatellite instability (MSI). Our study included 95 normal/tumor area pairs from 25 patients. Six of the 25 patients (24%) demonstrated various levels of MSI ranging from 7% (two of 30) to 97% (28 of 29) of markers tested in multiple tumor sites. Of the six patients, ®ve manifested frameshift mutations in a tract of ten deoxyadenosines within transforming growth factor b receptor type II and four demonstrated frameshift mutations in a tract of eight deoxyguanosines within BAX. These mutations were common to all tumor sites regardless of the various level of MSI phenotype, indicating initial events. Two of the six patients exhibited frameshift mutations in mononucleotide repeats of mismatch repair genes, hMSH3 and hMSH6, and the insulin-like growth factor II receptor in restricted tumor areas, indicating additional alterations. Insulin-like growth factor II receptor mutations appear to be caused by hMSH3 and hMSH6 mutations because the former mutations were con®ned to tumor portions with the latter two mismatch repair lesions. These results provide genetic progression evidence for gastric carcinomas of the mutator pathway. In this pathway, mismatch repair insuciency initially targets mononucleotide tracts of transforming growth factor b receptor type II and BAX. During tumorigenesis, primary mismatch repair failure may give rise to the secondary mismatch repair lesions, frameshift mutations of hMSH3 and hMSH6, which result in another tumorigenic mutation in the insulin-like growth factor II receptor.

CNVRuler software is available with an online manual at the website, www.ircgp.com/CNVRuler/index.html.