Search

Supplementary Figure S1. Transcriptomic analyses for molecular subtypes Supplementary Figure S2. Identification of the subtypes in LUSC TCGA cohort (n=431) Supplementary Figure S3. The genomic difference across subtypes in LUSC. Supplementary Figure S4. Estimation of Stromal and immune score and distribution of infiltrating immune cells for tumor and noncancer control sample in LUSC Supplementary Figure S5. The immune landscape and micro-environmental signature in LUSC TCGA cohort (n=431) Supplementary Figure S6. The dual roles of stroma are described by distinct gene expression patterns in tumor(subtype A and B) and noncancer control samples. Supplementary Figure S7. Distribution of infiltrating immune cells for tumor and noncancer control sample in LUSC Supplementary Figure S8. Distinct expression pattern of M1 and M2 signature genes in LUSC subtypes Supplementary Figure S9. Comparison of immune contents between subtypes in our LUSC cohort (n=101) and TCGA LUSC cohort (n=431) Supplementary Figure S10. Immune related pathways enriched for genes which have copy number variants. Supplementary Figure S11. Comparison of the LUSC expression subtypes with the previously defined subtypes Supplementary Figure S12. Clinical association analysis across subtypes Supplementary Figure S13. Survival analysis across subtypes

The clinical features of 101 LUSC patients enrolled in this study.

The immune microenvironment in lung squamous cell carcinoma (LUSC) is not well understood, with interactions between the host immune system and the tumor, as well as the molecular pathogenesis of LUSC, awaiting better characterization. To date, no molecularly targeted agents have been developed for LUSC treatment. Identification of predictive and prognostic biomarkers for LUSC could help optimize therapy decisions. We sequenced whole exomes and RNA from 101 tumors and matched noncancer control Korean samples. We used the information to predict subtype-specific interactions within the LUSC microenvironment and to connect genomic alterations with immune signatures. Hierarchical clustering based on gene expression and mutational profiling revealed subtypes that were either immune defective or immune competent. We analyzed infiltrating stromal and immune cells to further characterize the tumor microenvironment. Elevated expression of macrophage 2 signature genes in the immune competent subtype confirmed that tumor-associated macrophages (TAM) linked inflammation and mutation-driven cancer. A negative correlation was evident between the immune score and the amount of somatic copy-number variation (SCNV) of immune genes (r = −0.58). The SCNVs showed a potential detrimental effect on immunity in the immune-deficient subtype. Knowledge of the genomic alterations in the tumor microenvironment could be used to guide design of immunotherapy options that are appropriate for patients with certain cancer subtypes. Cancer Immunol Res; 6(7); 848–59. ©2018 AACR.

DHFR gene amplification is commonly present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. In this study, we proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real-time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). We identified an amplification unit spanning 11 genes, from the DHFR gene to the ATP6AP1L gene position, with high adjusted interaction frequencies on chromosome 5 (~2.2 Mbp) and a twenty-fold tandemly amplified region, and novel inversions at the start and end positions of the amplified region as well as frameshift insertions in most of the MSH and MLH genes were detected. These mutations might stimulate chromosomal breakage and cause the dysregulation of mismatch repair. Characterizing the tandem gene-amplified unit may be critical for identifying the mechanisms that trigger genomic rearrangements. These findings may provide new insight into the mechanisms underlying the amplification process and the evolution of drug resistance., Cancer chemotherapy: Sequencing surplus gene copies Sequencing a large region of DNA containing many surplus copies of genes linked to drug resistance in colon cancer cells may illuminate how these genomic rearrangements arise. Such regions of gene amplification are highly repetitive, making them impossible to sequence using ordinary methods, and little is known about how they are generated. Using advanced methods, Jeong-Sun Seo at Seoul National University Bundang Hospital in South Korea and co-workers sequenced a region of gene amplification in colon cancer cells. The amplified region was approximately 20 times the length of that in healthy cells and contained many copies of an eleven-gene segment, including a gene implicated in drug resistance. The region also contained mutations in chromosomal repair genes which would disrupt repair pathways. These results illuminate the genetic changes that lead to gene amplification and drug resistance in cancer cells.

Background Low-pass sequencing (LPS) has been extensively investigated for applicability to various genetic studies due to its advantages over genotype array data including cost-effectiveness. Predicting the risk of complex diseases such as Parkinson’s disease (PD) using polygenic risk score (PRS) based on the genetic variations has shown decent prediction accuracy. Although ultra-LPS has been shown to be effective in PRS calculation, array data has been favored to the majority of PRS analysis, especially for PD. Results Using eight high-coverage WGS, we assessed imputation approaches for downsampled LPS data ranging from 0.5 × to 7.0 × . We demonstrated that uncertain genotype calls of LPS diminished imputation accuracy, and an imputation approach using genotype likelihoods was plausible for LPS. Additionally, comparing imputation accuracies between LPS and simulated array illustrated that LPS had higher accuracies particularly at rare frequencies. To evaluate ultra-low coverage data in PRS calculation for PD, we prepared low-coverage WGS and genotype array of 87 PD cases and 101 controls. Genotype imputation of array and downsampled LPS were conducted using a population-specific reference panel, and we calculated risk scores based on the PD-associated SNPs from an East Asian meta-GWAS. The PRS models discriminated cases and controls as previously reported when both LPS and genotype array were used. Also strong correlations in PRS models for PD between LPS and genotype array were discovered. Conclusions Overall, this study highlights the potentials of LPS under 1.0 × followed by genotype imputation in PRS calculation and suggests LPS as attractive alternatives to genotype array in the area of precision medicine for PD.

The underrepresentation of non-Europeans in human genetic studies so far has limited the diversity of individuals in genomic datasets and led to reduced medical relevance for a large proportion of the world’s population. Population-specific reference genome datasets as well as genome-wide association studies in diverse populations are needed to address this issue. Here we describe the pilot phase of the GenomeAsia 100K Project. This includes a whole-genome sequencing reference dataset from 1,739 individuals of 219 population groups and 64 countries across Asia. We catalogue genetic variation, population structure, disease associations and founder effects. We also explore the use of this dataset in imputation, to facilitate genetic studies in populations across Asia and worldwide., Using whole-genome sequencing data from 1,739 individuals, the GenomeAsia 100K Project catalogues genetic variation, population structure and disease associations to facilitate genetic studies in Asian populations and increase representation in genetics studies worldwide.

Background We characterized the microbiome composition from various samples in pregnant women and their neonates and evaluated the association between the distribution of microbiome and obstetric factors. Results A prospective study was performed in live births delivered at a single tertiary center between March 2020 and January 2021. Samples collected from pregnant women and their newborns included maternal cervicovaginal discharge (VD), amniotic fluid (AF), umbilical cord blood (CB), neonatal gastric liquid (GL), and meconium. Next-generation sequencing was used to analyze the composition of the microbiome from the samples. We identified 19,597,239 bacterial sequences from the 641 samples. The alpha diversity of the negative control group was slightly higher than that of the VD group. The samples were separated well based on the microbiome pattern unique to their body sites despite the significant batch effects present within our dataset. We observed a higher relative abundance of Lactobacillus at the genus level but lower abundances of Aerococcus, Fusobacterium, Gardnerella, Peptoniphilus, Porphyromonas, and Prevotella in VD. In evaluating the association with clinical data, Staphylococcus in M showed a strong association with preterm birth. The lists of bacterial taxa were connected to the weights of newborn babies in the CB group. The microbiome pattern of GL was similar to that of AF, but not technically identical. Twins did not demonstrate significant differences in the composition of microbiomes from neonate samples. Conclusions In this study, samples known to be traditionally sterile, such as AF and CB, showed specific bacterial taxa although they were extremely close to negative controls. The microbiome patterns were grouped according to the body compartments where the samples had been obtained, not according to the mother-fetus/infant pairs. Further studies lowering the possibility of contamination and redeeming the low-biomass specimens are necessary to determine the relationship between microbiomes and specific pathological complications during pregnancy.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in “Charrúas” and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

High-grade serous ovarian cancer (HGSOC) is one of the deadliest cancers that can occur in women. This study aimed to investigate the molecular characteristics of HGSOC through integrative analysis of multi-omics data. We used fresh-frozen, chemotherapy-naïve primary ovarian cancer tissues and matched blood samples of HGSOC patients and conducted next-generation whole-exome sequencing (WES) and RNA sequencing (RNA-seq). Genomic and transcriptomic profiles were comprehensively compared between patients with germline BRCA1/2 mutations and others with wild-type BRCA1/2. HGSOC samples initially divided into two groups by the presence of germline BRCA1/2 mutations showed mutually exclusive somatic mutation patterns, yet the implementation of high-dimensional analysis of RNA-seq and application of epithelial-to-mesenchymal (EMT) index onto the HGSOC samples revealed that they can be divided into two subtypes, homologous recombination repair (HRR)-activated type and mesenchymal type. Patients with mesenchymal HGSOC, characterized by the activation of the EMT transcriptional program, low genomic alteration and diverse cell-type compositions, exhibited significantly worse overall survival than did those with HRR-activated HGSOC (p = 0.002). In validation with The Cancer Genome Atlas (TCGA) HGSOC data, patients with a high EMT index (≥the median) showed significantly worse overall survival than did those with a low EMT index (&lt, the median) (p = 0.030). In conclusion, through a comprehensive multi-omics approach towards our HGSOC cohorts, two distinctive types of HGSOC (HRR-activated and mesenchymal) were identified. Our novel EMT index seems to be a potential prognostic biomarker for HGSOC.

Background Rare diseases are pathologies that affect less than 1 in 2000 people. They are difficult to diagnose due to their low frequency and their often highly heterogeneous symptoms. Rare diseases have in general a high impact on the quality of life and life expectancy of patients, which are in general children or young people. The advent of high-throughput sequencing techniques has improved diagnosis in several different areas, from pediatrics, achieving a diagnostic rate of 41% with whole genome sequencing (WGS) and 36% with whole exome sequencing, to neurology, achieving a diagnostic rate between 47 and 48.5% with WGS. This evidence has encouraged our group to pursue a molecular diagnosis using WGS for this and several other patients with rare diseases. Results We used whole genome sequencing to achieve a molecular diagnosis of a 7-year-old girl with a severe panvascular artery disease that remained for several years undiagnosed. We found a frameshift variant in one copy and a large deletion involving two exons in the other copy of a gene called YY1AP1. This gene is related to Grange syndrome, a recessive rare disease, whose symptoms include stenosis or occlusion of multiple arteries, congenital heart defects, brachydactyly, syndactyly, bone fragility, and learning disabilities. Bioinformatic analyses propose these mutations as the most likely cause of the disease, according to its frequency, in silico predictors, conservation analyses, and effect on the protein product. Additionally, we confirmed one mutation in each parent, supporting a compound heterozygous status in the child. Conclusions In general, we think that this finding can contribute to the use of whole genome sequencing as a diagnosis tool of rare diseases, and in particular, it can enhance the set of known mutations associated with different diseases.

Anaplastic thyroid cancer (ATC) and advanced differentiated thyroid cancers (DTCs) show fatal outcomes, unlike DTCs. Here, we demonstrate mutational landscape of 27 ATCs and 86 advanced DTCs by massively-parallel DNA sequencing, and transcriptome of 13 ATCs and 12 advanced DTCs were profiled by RNA sequencing. TERT, AKT1, PIK3CA, and EIF1AX were frequently co-mutated with driver genes (BRAFV600E and RAS) in advanced DTCs as well as ATC, but tumor suppressors (e.g., TP53 and CDKN2A) were predominantly altered in ATC. CDKN2A loss was significantly associated with poor disease-specific survival in patients with ATC or advanced DTCs, and up-regulation of CD274 (PD-L1) and PDCD1LG2 (PD-L2). Transcriptome analysis revealed a fourth molecular subtype of thyroid cancer (TC), ATC-like, which hardly reflects the molecular signatures in DTC. Furthermore, the activation of JAK-STAT signaling pathway could be a potential druggable target in RAS-positive ATC. Our findings provide insights for precision medicine in patients with advanced TCs., Anaplastic thyroid cancer (ATC) and advanced differentiated thyroid cancers (DTCs) come with a dismal prognosis. Here, Yoo and colleagues reveal the genomic and transcriptomic landscape of ATC and DTC, highlighting potential therapeutic vulnerabilities.

Understanding the tumor microenvironment is important to efficiently identify appropriate patients for immunotherapies in a variety of cancers. Here, we presented the tumor microenvironmental analysis of 2,033 cancer samples across 7 cancer types: colon adenocarcinoma, skin cutaneous melanoma, kidney renal papillary cell carcinoma, sarcoma, pancreatic adenocarcinoma, glioblastoma multiforme, and pheochromocytoma / paraganglioma from The Cancer Genome Atlas cohort. Unsupervised hierarchical clustering based on the gene expression profiles separated the cancer samples into two distinct clusters, and characterized those into immune-competent and immune-deficient subtypes using the estimated abundances of infiltrated immune and stromal cells. We demonstrated differential tumor microenvironmental characteristics of immune-competent subtypes across 7 cancer types, particularly immunosuppressive tumor microenvironment features in kidney renal papillary cell carcinoma with significant poorer survival rates and immune-supportive features in sarcoma and skin cutaneous melanoma. Additionally, differential genomic instability patterns between the subtypes were found across the cancer types, and discovered that immune-competent subtypes in most of cancer types had significantly higher immune checkpoint gene expressions. Overall, this study suggests that our subtyping approach based on transcriptomic data could contribute to precise prediction of immune checkpoint inhibitor responses in a wide range of cancer types.

The pathogenesis of Alzheimer's disease (AD) and the commonest cause of dementia in the elderly remain incompletely understood. Recently, epigenetic modifications have been shown to play a potential role in neurodegeneration, but the specific involvement of epigenetic signatures landscaped by heterochromatin has not been studied in AD. Herein, we discovered that H3K9me3‐mediated heterochromatin condensation is elevated in the cortex of sporadic AD postmortem brains. In order to identify which epigenomes are modulated by heterochromatin, we performed H3K9me3‐chromatin immunoprecipitation (ChIP)‐sequencing and mRNA‐sequencing on postmortem brains from normal subjects and AD patients. The integrated analyses of genome‐wide ChIP‐ and mRNA‐sequencing data identified epigenomes that were highly occupied by H3K9me3 and inversely correlated with their mRNA expression levels in AD. Biological network analysis further revealed H3K9me3‐landscaped epigenomes to be mainly involved in synaptic transmission, neuronal differentiation, and cell motility. Together, our data show that the abnormal heterochromatin remodeling by H3K9me3 leads to down‐regulation of synaptic function‐related genes, suggesting that the epigenetic alteration by H3K9me3 is associated with the synaptic pathology of sporadic AD., Histone H3K9me3 is densely stained in nuclear foci where heterochromatins are spatially localized in neurons of AD temporal cortex. H3k9me3‐heterochromatin‐landscaped epigenomes are closely associated with the synaptic dysfunction in AD. Together, the epigenetically disrupted biological pathways could serve as novel biomarkers in AD.

DHFRgene amplification is present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. However, little is known aboutDHFRgene amplification due to difficulties in quantifying amplification size and recognizing the repetitive rearrangements involved in the process. In this study, we have proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). Amplification of theDHFRgene was optimized to generate homogenously amplified patterns. The amplification units of 11 genes, from theDHFRgene to theATP6AP1Lgene position on chromosome 5 (~2.2Mbp), and a twenty-fold tandemly amplified region were verified using long-range genome and RNA sequencing data. In doing so, a novel inversion at the start and end positions of the amplified region as well as frameshift insertions in most of theMSHandMLHgenes were detected. These might stimulate chromosomal breakage and cause the dysregulation of mismatch repair pathways. Using Hi-C technology, high adjusted interaction frequencies were detected on the amplified unit and unsuspected position on 5q, which could have a complex network of spatial contacts to harbor gene amplification. Characterizing the tandem gene-amplified unit and genomic variants as well as chromosomal interactions on intra-chromosome 5 can be critical in identifying the mechanisms behind genomic rearrangements. These findings may give new insight into the mechanisms underlying the amplification process and evolution of drug resistance.

Recent work has demonstrated that two archaic human groups (Neanderthals and Denisovans) interbred with modern humans and contributed to the contemporary human gene pool. These findings relied on the availability of high-coverage genomes from both Neanderthals and Denisovans. Here we search for evidence of archaic admixture from a worldwide panel of 1,667 individuals using an approach that does not require the presence of an archaic human reference genome. We find no evidence for archaic admixture in the Andaman Islands, as previously claimed, or on the island of Flores, where Homo floresiensis fossils have been found. However, we do find evidence for at least one archaic admixture event in sub-Saharan Africa, with the strongest signal in Khoesan and Pygmy individuals from Southern and Central Africa. The locations of these putative archaic admixture tracts are weighted against functional regions of the genome, consistent with the long-term effects of purifying selection against introgressed genetic material.

Objective The aim of this study was to develop a common targeted massively parallel sequencing platform for the noninvasive prenatal diagnosis (NIPD) of multiple X-linked diseases. Method The custom capture probe was designed to target 33 genes and recombination hotspots. We tested the carrier mother and male proband pair of 6 families. Plasma DNA of the pregnant carrier mother was collected at different gestational weeks and sequenced. The fetal genotype of each family was determined by estimating the imbalance between the 2 maternal haplotypes constructed using a common custom-designed platform. Results The targeted sequencing of the maternal, proband, and fetal genomic DNAs and maternal plasma DNAs resulted in uniform coverage across the target region. Three to 5 recombination points were observed in each sample. However, these recombination points did not affect the haplotype dosage analysis for fetal genotype prediction. Consequently, all fetal genotypes in the 6 families obtained from haplotype dosage analysis of maternal plasma sequencing data were predicted correctly. Conclusions Since a single platform that covers multiple diseases may prevent the need for disease-specific probes for the NIPD of individual disorders, this approach may provide a practical advantage for clinically implementing the NIPD of X-linked diseases.

We aimed to develop forest therapy programs (FTPs) to prevent dementia and related health problems in the elderly population, with the assumption that health benefits are FTP-type specific and depend on the participant&rsquo, s psychophysiological traits. For this purpose, we developed two distinct FTPs, namely, a guided-breathing meditation program (BP) and a walking program (WP), we adopted the approach of Sasang constitutional (SC) medicine, which categorizes individuals into one of three SC types (SC1, SC2, or SC3) for medical care. The FTPs ran 11 sessions over 11 weeks. We recruited 29/31/28 participants who were 65 years of age or older for the BP/WP/control groups, respectively, obtained electrophysiological measurements via electroencephalogram (EEG), heart rate variability (HRV), and bioimpedance, and analyzed the intervention effects with analysis of covariance. Compared with the control, the BP and WP resulted in benefits for neural activity and parasympathetic nervous activity (PNA), respectively, and both FTPs yielded distinct beneficial effects on bioimpedance. Constitution-specific effects were also present. The SC1- and SC2-type participants gained positive effects in neural activity from the BP and WP, respectively. The SC3-type participants showed improvements in PNA from the WP. In conclusion, for older individuals, both programs conferred health benefits that would help prevent dementia, and the benefits were program-specific and constitution-specific.

The immune microenvironment in lung squamous cell carcinoma (LUSC) is not well understood, with interactions between the host immune system and the tumor, as well as the molecular pathogenesis of LUSC, awaiting better characterization. To date, no molecularly targeted agents have been developed for LUSC treatment. Identification of predictive and prognostic biomarkers for LUSC could help optimize therapy decisions. We sequenced whole exomes and RNA from 101 tumors and matched noncancer control Korean samples. We used the information to predict subtype-specific interactions within the LUSC microenvironment and to connect genomic alterations with immune signatures. Hierarchical clustering based on gene expression and mutational profiling revealed subtypes that were either immune defective or immune competent. We analyzed infiltrating stromal and immune cells to further characterize the tumor microenvironment. Elevated expression of macrophage 2 signature genes in the immune competent subtype confirmed that tumor-associated macrophages (TAM) linked inflammation and mutation-driven cancer. A negative correlation was evident between the immune score and the amount of somatic copy-number variation (SCNV) of immune genes (r = −0.58). The SCNVs showed a potential detrimental effect on immunity in the immune-deficient subtype. Knowledge of the genomic alterations in the tumor microenvironment could be used to guide design of immunotherapy options that are appropriate for patients with certain cancer subtypes. Cancer Immunol Res; 6(7); 848–59. ©2018 AACR.

Background Fusion genes are good candidates of molecular targets for cancer therapy. However, there is insufficient research on the clinical implications and functional characteristics of fusion genes in colorectal cancer (CRC). Methods In this study, we analysed RNA sequencing data of CRC patients (147 tumour and 47 matched normal tissues) to identify oncogenic fusion genes and evaluated their role in CRC. Results We validated 24 fusion genes, including novel fusions, by three algorithms and Sanger sequencing. Fusions from most patients were mutually exclusive CRC oncogenes and included tumour suppressor gene mutations. Eleven fusion genes from 13 patients (8.8%) were determined as oncogenic fusion genes by analysing their gene expression and function. To investigate their oncogenic impact, we performed proliferation and migration assays of CRC cell lines expressing fusion genes of GTF3A-CDK8, NAGLU- IKZF3, RNF121- FOLR2, and STRN-ALK. Overexpression of these fusion genes increased cell proliferation except GTF3A-CDK8. In addition, overexpression of NAGLU-IKZF3 enhanced migration of CRC cells. We demonstrated that NAGLU-IKZF3, RNF121-FOLR2, and STRN-ALK had tumourigenic effects in CRC. Conclusion In summary, we identified and characterised oncogenic fusion genes and their function in CRC, and implicated NAGLU-IKZF3 and RNF121-FOLR2 as novel molecular targets for personalised medicine development.

For the noninvasive prenatal diagnosis (NIPD) of X-linked recessive diseases such as Duchenne muscular dystrophy (DMD), maternal haplotype phasing is a critical step for dosage analysis of the inherited allele. Until recently, the proband-based indirect haplotyping method has been preferred despite its limitations for use in clinical practice. Here, we describe a method for directly determining the maternal haplotype without requiring the proband’s DNA in DMD families. We used targeted linked-read deep sequencing (mean coverage of 692×) of gDNA from 5 mothers to resolve their haplotypes and predict the mutation status of the fetus. The haplotype of DMD alleles in the carrier mother was successfully phased through a targeted linked-read sequencing platform. Compared with the proband-based phasing method, linked-read sequencing was more accurate in differentiating whether the recombination events occurred in the proband or in the fetus. The predicted inheritance of the DMD mutation was diagnosed correctly in all 5 families in which the mutation had been confirmed using amniocentesis or chorionic villus sampling. Direct haplotyping by this targeted linked-read sequencing method could be used as a phasing method for the NIPD of DMD, especially when the genomic DNA of the proband is unavailable.

In the original article, part of Table 1 headings and entries were missing. The correct Table is as shown below. The original article has been corrected.

We aimed to directly convert adult human dermal fibroblasts (aHDFs) into functional endothelial cells (ECs). Lentiviral vectors encoding endothelial transcription factors (TFs) were constructed. We examined whether five TFs (FOXO1, ER71, KLF2, TAL1, and LMO2) used for the generation of mouse induced ECs (iECs) could convert the aHDFs into human iECs. Twenty-eight days after transduction with lentiviral constructs, 32.1 ± 5.1% cells expressed vascular endothelial (VE)-cadherin. Factor screening revealed that only three factors (3F: ER71, KLF2, and TAL1) were necessary to induce VE-cadherin (+) cells (49.4 ± 3.5%). However, whole transcriptome sequencing showed that VE-cadherin (+) cells were not completely reprogrammed. Mature iECs double-positive for VE-cadherin/Pecam1 (DP cells) with a cobblestone appearance were obtained at a frequency of only 5.1 ± 0.6%. Using whole transcriptome analysis, the potential factors that could block the conversion were screened. Among candidates TWIST1-knockdown enhanced efficiency of conversion. Rosiglitazone, an inhibitor of epithelial-mesenchymal transition (EMT), also improved the conversion efficiency. Moreover, a 2nd second-stage conversion process, in which VE-cadherin (+) cells were incubated for additional two weeks, further enhanced the efficiency. The final protocol for 6 weeks yielded a conversion rate of 19.6 ± 3.0% iECs, defined by DP cells depicting the nature of mature ECs in various analyses. Further analyses revealed that the genetic and epigenetic profiles of iECs resembled those of functional ECs. Collectively, aHDFs can be converted into functional ECs through the transduction of ER71, KLF2, and TAL1, combined with two EMT inhibitors (siTWIST1 and rosiglitazone), followed by 2nd stage conversion.

Chronic traumatic encephalopathy (CTE) is a neurodegenerative disease that is associated with repetitive traumatic brain injury (TBI). CTE is known to share similar neuropathological features with Alzheimer’s disease (AD), but little is known about the molecular properties in CTE. To better understand the neuropathological mechanism of TBI-related disorders, we conducted transcriptome sequencing analysis of CTE including AD and CTE with AD (CTE/AD) post-mortem human brain samples. Through weighted gene co-expression network analysis (WGCNA) and principal component analysis (PCA), we characterized common and unique transcriptome signatures among CTE, CTE/AD, and AD. Interestingly, synapse signaling-associated gene signatures (such as synaptotagmins) were commonly down-regulated in CTE, CTE/AD, and AD. Quantitative real-time PCR (qPCR) and Western blot analyses confirmed that the levels of synaptotagmin 1 (SYT1) were markedly decreased in CTE and AD compared to normal. In addition, calcium/calmodulin-dependent protein kinase II (CaMKII), protein kinase A (PKA), protein kinase C (PKC), and AMPA receptor genes that play a pivotal role in memory function, were down-regulated in head trauma-related disorders. On the other hand, up-regulation of cell adhesion molecules (CAMs) associated genes was only found in CTE. Our results indicate that dysregulation of synaptic transmission- and memory function-related genes are closely linked to the pathology of head injury-related disorder and AD. Alteration of CAMs-related genes may be specific pathological markers for the CTE pathology.

Genotype imputation using the reference panel is a cost-effective strategy to fill millions of missing genotypes for the purpose of various genetic analyses. Here, we present the Northeast Asian Reference Database (NARD), including whole-genome sequencing data of 1,781 individuals from Korea, Mongolia, Japan, China, and Hong Kong. NARD provides the genetic diversities of Korean (n=850) and Mongolian (n=386) ancestries that were not present in the 1000 Genomes Project Phase 3 (1KGP3). We combined and re-phased the genotypes from NARD and 1KGP3 to construct a union set of haplotypes. This approach established a robust imputation reference panel for the Northeast Asian populations, which yields the greatest imputation accuracy of rare and low-frequency variants compared with the existing panels. Also, we illustrate that NARD can potentially improve disease variant discovery by reducing pathogenic candidates. Overall, this study provides a decent reference panel for the genetic studies in Northeast Asia.

Additional file 1: Figure S1. Geographic map of the study area in the NARD. Figure S2. Correlation between the sequencing depth and number of variants. Figure S3. Transition to transversion ratio of the populations in the NARD. Figure S4. Heterozygous to homozygous ratio of the global populations. Figure S5. Number of loss-of-function variants. Figure S6. Hardy-Weinberg Equilibrium test of variants in the NARD. Figure S7. Novel variant statistics. Figure S8. Differential genetic composition of the two MNG groups. Figure S9. Imputation performance evaluation of FRA individuals. Figure S10. Imputation performance evaluation of CHN and JPN individuals. Figure S11. Length distribution of shared IBD tracts between the two individuals in each population. Figure S12. The flow chart of the pipeline consisting of four major steps for NARD imputation server. Figure S13. The cross-validation error inferred by ADMIXTURE algorithm.

The identification and clinical validation of cancer driver genes are essential to accelerate the translational transition of cancer genomics, as well as to find clinically confident targets for the therapeutic intervention of cancers. Here we identified recurrent LMNA-NTRK1 and TPM3-NTRK1 fusions in Korean patients with colon cancer (3 out of 147, 2%) through next-generation RNA sequencing (RNA-seq). NTRK1 fusions were mutually exclusive oncogenic drivers of colon cancer that were accompanied with in vitro potential of colony formation and in vivo tumorigenicity comparable to KM12, a human colon cancer cell line harboring TPM3-NTRK1 fusion. NTRK1-encoded TrkA protein was prevalent in 11 out of 216 Korean (5.1%) and 28 out of 472 Chinese patients (5.9%) from independent cohorts, respectively. The expression level of TrkA was significantly correlated with NTRK1 fusion (p = 0.0192), which was verified by a fluorescence in situ hybridization (FISH). Korean patients with TrkA-positive colon cancer had a marginal but significant shorter overall survival time than TrkA-negative colon cancer [hazard ratio (HR) = 0.5346, 95% confidential interval (CI) = 0.2548-0.9722, p = 0.0411]. In addition, KM12 cell line was sensitive to selective TrkA inhibitors. These results demonstrate that NTRK1 fusion is granted as a clinically relevant target for therapeutic intervention of colon cancer.

We report a novel strategy to directly convert adult human dermal fibroblasts (aHDFs) into authentic endothelial cells (ECs). Lentiviral vectors encoding endothelial transcription factors (TFs) were constructed, used to transduce aHDFs. We examined whether 5 TFs (FOXO1, ER71, KLF2, TAL1, and LMO2

Peroxiredoxin5 (Prdx5), a thioredoxin peroxidase, is an antioxidant enzyme that is widely studied for its antioxidant properties and protective roles in neurological and cardiovascular disorders. This study is aimed at investigating the functional significance of Prdx5 in mitochondria and at analyzing its roles in ciliogenesis during the process of vertebrate development.We found that several Prdx genes were strongly expressed in multiciliated cells in developing Xenopus embryos, and their peroxidatic functions were crucial for normal cilia development. Depletion of Prdx5 increased levels of cellular reactive oxygen species (ROS), consequently leading to mitochondrial dysfunction and abnormal cilia formation. Proteomic and transcriptomic approaches revealed that excessive ROS accumulation on Prdx5 depletion subsequently reduced the expression level of pyruvate kinase (PK), a key metabolic enzyme in energy production. We further confirmed that the promotor activity of PK was significantly reduced on Prdx5 depletion and that the reduction in PK expression and its promoter activity led to ciliary defects observed in Prdx5-depleted cells.Our data revealed the novel relationship between ROS and Prdx5 and the consequent effects of this interaction on vertebrate ciliogenesis. The normal process of ciliogenesis is interrupted by the Prdx5 depletion, resulting in excessive ROS levels and suggesting cilia as vulnerable targets of ROS.Prdx5 plays protective roles in mitochondria and is critical for normal cilia development by regulating the levels of ROS. The loss of Prdx5 is associated with excessive production of ROS, resulting in mitochondrial dysfunction and aberrant ciliogenesis.

Characterizing the molecular immune subtype and micro-environment of lung cancer is necessary to understand immunogenic interactions between infiltrating immune and stromal cells, and how tumor cells overcome immune checkpoint blockades. This study seeks to identify computational methodologies for subtyping gene expression-based tumor-immune micro-environment interactions, which differentiate non-small cell lung cancer (NSCLC) into immune-defective and immune-competent subtypes. Here, 101 lung squamous cell carcinomas (LUSCs) and 87 lung adenocarcinomas (LUADs) tumor samples have been analyzed. Several micro-environmental factors differentially induce LUAD or LUSC immune subtypes, as well as immune checkpoint expression. In particular, tumor-associated macrophages (TAMs) are key immune cells play a vital role in inflammation and cancer micro-environments of LUSCs; whereas, regulatory B cells are immunosuppressive and tumorigenic in LUADs. Additionally, cytolytic activity upon CD8+ T cell activation is decreased by the abundance of B cells and macrophages in immune-competent subtypes. Therefore, identifying immune subtypes in lung cancer and their impact on tumor micro-environment will lead to clinical tools for assessing LUADs and LUSCs in patients, as well as maximize the efficacy of immune checkpoint inhibitors.

Pluripotency is defined by the ability of a cell to differentiate to the derivatives of all the three embryonic germ layers: ectoderm, mesoderm and endoderm. Pluripotent cells can be captured via the archetypal derivation of embryonic stem cells or via somatic cell reprogramming. Somatic cells are induced to acquire a pluripotent stem cell (iPSC) state through the forced expression of key transcription factors, and in the mouse these cells can fulfil the strictest of all developmental assays for pluripotent cells by generating completely iPSC-derived embryos and mice. However, it is not known whether there are additional classes of pluripotent cells, or what the spectrum of reprogrammed phenotypes encompasses. Here we explore alternative outcomes of somatic reprogramming by fully characterizing reprogrammed cells independent of preconceived definitions of iPSC states. We demonstrate that by maintaining elevated reprogramming factor expression levels, mouse embryonic fibroblasts go through unique epigenetic modifications to arrive at a stable, Nanog-positive, alternative pluripotent state. In doing so, we prove that the pluripotent spectrum can encompass multiple, unique cell states.

Fluorescence-guided surgery using 5-aminolevulinic acid (5-ALA) is now a widely-used modality for glioblastoma (GBM) treatment. However, intratumoral heterogeneity of fluorescence intensity may reflect different onco-metabolic programs. Here, we investigated the metabolic mechanism underlying the heterogeneity of 5-ALA fluorescence in GBM. Using an in-house developed fluorescence quantification system for tumor tissues, we collected 3 types of GBM tissues on the basis of their fluorescence intensity, which was characterized as strong, weak, and none. Expression profiling by RNA-sequencing revealed 77 genes with a proportional relationship and 509 genes with an inverse relationship between gene expression and fluorescence intensity. Functional analysis and in vitro experiments confirmed glutaminase 2 (GLS2) as a key gene associated with the fluorescence heterogeneity. Subsequent metabolite profiling discovered that insufficient NADPH due to GLS2 underexpression was responsible for the delayed metabolism of 5-ALA and accumulation of protoporphyrin IX (PpIX) in the high fluorescence area. The expression level of GLS2 and related NADPH production capacity is associated with the regional heterogeneity of 5-ALA fluorescence in GBM.

Although p21(WAF1/CIP1) is known to be elevated during replicative senescence of human embryonic fibroblasts (HEFs), the mechanism for p21 up-regulation has not been elucidated clearly. In order to explore the mechanism, we analyzed expression of p21 mRNA and protein and luciferase activity of full-length p21 promoter. The result demonstrated that p21 up-regulation was accomplished largely at transcription level. The promoter assay using serially-deleted p21 promoter constructs revealed that p53 binding site was the most important site and Sp1 binding sites were necessary but not sufficient for transcriptional activation of p21. In addition, p53 protein was shown to interact with Sp1 protein. The interaction was increased in aged fibroblasts and was regulated by phosphorylation of p53 and Sp1. DNA binding activity of p53 was significantly elevated in aged fibroblasts but that of Sp1 was not. DNA binding activities of p53 and Sp1 were also regulated by phosphorylation. Phosphorylation of p53 at serine-15 and of Sp1 at serines appears to be involved. Taken together, the result demonstrated that p21 transcription during replicative senescence of HEFs is up-regulated by increase in DNA binding activity and interaction between p53 and Sp1 via phosphorylation.

Anaplastic lymphoma receptor tyrosine kinase ( ALK ) is located on chromosome 2p23; the chromosomal rearrangements of this gene are common genetic alterations, resulting in the creation of multiple fusion genes involved in tumorigenesis. However, the presence of an ALK fusion in myeloid malignancies is extremely rare. We report a case of acute myelomonocytic leukemia in a 31-year-old woman with an unusual rearrangement between RAN-binding protein 2 ( RANBP2 ) and ALK and a karyotype of 45,XX,inv(2)(p23q21),–7[20]. We detected an ALK rearrangement using fluorescence in situ hybridization, identified the ALK fusion partner by using RNA transcriptome sequencing, and demonstrated the RANBP2 – ALK fusion transcript by reverse transcriptase–PCR and Sanger sequencing. Immunohistochemistry for ALK showed strong staining of the nuclear membrane in leukemic cells. The patient had an unfavorable clinical course. Our results, together with a literature review, suggest the RANBP2 - ALK fusion combined with monosomy 7 may be related to a unique clonal hematologic disorder of childhood and adolescence, characterized by myelomonocytic leukemia and a poor prognosis.

Targeted resequencing using next-generation sequencing technology is being rapidly applied to the molecular diagnosis of human genetic diseases. The group of muscular dystrophies may be an appropriate candidate for this approach because these diseases exhibit genotype–phenotype heterogeneity. To perform a proof-of-concept study, we selected four patients with congenital muscular dystrophies with defective glycosylation of alpha-dystroglycan. A custom-solution-based target enrichment kit was designed to capture whole-genic regions of the 26 muscular-dystrophy-related genes, including six genes implicated in alpha-dystroglycanopathies. Although approximately 95% of both coding and noncoding regions were covered with at least 15-read depth, parts of the coding exons of FKRP and POMT2 were insufficiently covered. Homozygous and compound heterozygous POMGnT1 mutations were found in two patients. Two novel noncoding variants of FKTN were identified in one patient who had a retrotransposon insertion mutation of FKTN in only one allele. The current targeted resequencing strategy yielded promising results for the extension of this method to other muscular dystrophies. As suboptimal coverage in a small subset of coding regions may affect the sensitivity of the method, complementary Sanger sequencing may be required.