Interleukin-34 Limits the Therapeutic Effects of Immune Checkpoint Blockade

Interleukin-34 (IL-34) is an alternative ligand to colony-stimulating factor-1 (CSF-1) for the CSF-1 receptor that acts as a key regulator of monocyte/macrophage lineage. In this study, we show that tumor-derived IL-34 mediates resistance to immune checkpoint blockade regardless of CSF-1 existence in various murine cancer models. Consistent with its immunosuppressive characteristics, the expression of IL-34 in tumors correlates with decreased frequencies of cellular (such as CD8+ and CD4+ T cells and M1-biased macrophages) and molecular (including various cytokines and chemokines) effectors at the tumor microenvironment. Then, a neutralizing antibody against IL-34 improved the therapeutic effects of the immune checkpoint blockade in combinatorial therapeutic models, including a patient-derived xenograft model. Collectively, we revealed that tumor-derived IL-34 inhibits the efficacy of immune checkpoint blockade and proposed the utility of IL-34 blockade as a new strategy for cancer therapy.

PDX models of human lung squamous cell carcinoma: consideration of factors in preclinical and co-clinical applications

Background

Treatment of human lung squamous cell carcinoma (LUSC) using current targeted therapies is limited because of their diverse somatic mutations without any specific dominant driver mutations. These mutational diversities preventing the use of common targeted therapies or the combination of available therapeutic modalities would require a preclinical animal model of this tumor to acquire improved clinical responses. Patient-derived xenograft (PDX) models have been recognized as a potentially useful preclinical model for personalized precision medicine. However, whether the use of LUSC PDX models would be appropriate enough for clinical application is still controversial.

Methods

In the process of developing PDX models from Korean patients with LUSC, the authors investigated the factors influencing the successful initial engraftment of tumors in NOD scid gamma mice and the retainability of the pathological and genomic characteristics of the parental patient tumors in PDX tumors.

Conclusions

The authors have developed 62 LUSC PDX models that retained the pathological and genomic features of parental patient tumors, which could be used in preclinical and co-clinical studies.

Trial registration Tumor samples were obtained from 139 patients with LUSC between November 2014 and January 2019. All the patients provided signed informed consents. This study was approved by the institutional review board (IRB) of Samsung Medical Center (2018-03-110)

The clinical efficacy of PSMA PET/MRI in biochemically recurrent prostate cancer compared with standard of care imaging modalities and confirmatory histopathology: results of a single-centre, prospective clinical trial

Prospective evidence for the clinical role and efficacy of prostate specific membrane antigen (PSMA) positron emission tomography (PET)/magnetic resonance imaging (MRI) combining MRI characterization and localization of lesions with PET avidity in comparison to conventional imaging is limited. In a prospective clinical trial, we aimed to evaluate the diagnostic yield and therapeutic impact of PSMA PET/MRI in men with biochemical recurrence (BCR) following curative therapy. A single-centre, prospective clinical trial at the Princess Alexandra Hospital recruited 30 patients with BCR. Patients underwent PSMA PET/MRI and concurrent conventional CT chest, abdomen, pelvis and whole-body bone scan. Biopsy was performed when safety possible for histological correlation of identified lesions. Clinical efficacy and impact of PSMA PET findings were evaluated. 30 patients with BCR were recruited (median PSA 0.69 ng/ml). PSMA avid lesions were present in 21 patients (70%). 23 patients were previously treated with definitive surgery, 6 patients received external beam radiotherapy and 1 patient had low dose rate brachytherapy. A total of 8 of 9 lesions biopsied were positive (88.9% histological correlation). PSMA PET/MRI detected local recurrence (p = 0.005) and pelvic lesions (p = 0.06) more accurately than conventional imaging. PSMA PET/MRI may be useful in staging men with biochemical recurrence, especially when PSA is low. Our data demonstrates a high detection rate, especially for locally recurrent disease, and highlights the role of this modality when PSA is low. This modality has the potential to significantly improve prostate cancer detection and may have implications for earlier salvage treatment, avoidance of futile local therapy and change patient management to lead to improved outcomes.

Validity of patient-derived xenograft mouse models for lung cancer based on exome sequencing data

Patient-derived xenograft (PDX) mouse models are frequently used to test the drug efficacy in diverse types of cancer. They are known to recapitulate the patient characteristics faithfully, but a systematic survey with a large number of cases is yet missing in lung cancer. Here we report the comparison of genomic characters between mouse and patient tumor tissues in lung cancer based on exome sequencing data. We established PDX mouse models for 132 lung cancer patients and performed whole exome sequencing for trio samples of tumor-normal-xenograft tissues. Then we computed the somatic mutations and copy number variations, which were used to compare the PDX and patient tumor tissues. Genomic and histological conclusions for validity of PDX models agreed in most cases, but we observed eight (~7%) discordant cases. We further examined the changes in mutations and copy number alterations in PDX model production and passage processes, which highlighted the clonal evolution in PDX mouse models. Our study shows that the genomic characterization plays complementary roles to the histological examination in cancer studies utilizing PDX mouse models.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.

Metabolic Syndrome Prediction Using Machine Learning Models with Genetic and Clinical Information from a Nonobese Healthy Population

The prevalence of metabolic syndrome (MS) in the nonobese population is not low. However, the identification and risk mitigation of MS are not easy in this population. We aimed to develop an MS prediction model using genetic and clinical factors of nonobese Koreans through machine learning methods. A prediction model for MS was designed for a nonobese population using clinical and genetic polymorphism information with five machine learning algorithms, including naïve Bayes classification (NB). The analysis was performed in two stages (training and test sets). Model A was designed with only clinical information (age, sex, body mass index, smoking status, alcohol consumption status, and exercise status), and for model B, genetic information (for 10 polymorphisms) was added to model A. Of the 7,502 nonobese participants, 647 (8.6%) had MS. In the test set analysis, for the maximum sensitivity criterion, NB showed the highest sensitivity: 0.38 for model A and 0.42 for model B. The specificity of NB was 0.79 for model A and 0.80 for model B. In a comparison of the performances of models A and B by NB, model B (area under the receiver operating characteristic curve [AUC] = 0.69, clinical and genetic information input) showed better performance than model A (AUC = 0.65, clinical information only input). We designed a prediction model for MS in a nonobese population using clinical and genetic information. With this model, we might convince nonobese MS individuals to undergo health checks and adopt behaviors associated with a preventive lifestyle.

Initial Experience With Hand-Assisted Laparoscopic Living Donor Nephrectomy: Training and Clinical Practice as a General Surgeon

BACKGROUND

To analyze our initial results of hand-assisted laparoscopic living donor nephrectomy, executed by a skilled gastrointestinal surgeon.

METHODS

A total of 22 consecutive patients underwent the hand-assisted laparoscopic living donor nephrectomy between December 2014 and January 2017. We retrospectively analyze the patient's perioperative clinical data, which were collected prospectively.

RESULTS

The right kidney was harvested in 12 patients. The mean operative time and intraoperative blood loss was 241.0 ± 43.4 minutes (range, 140-310 min) and 293.2 ± 203.1 mL (range, 50-700 mL), respectively. The mean warm ischemic time was 288.4 ± 103.4 seconds (range, 179-610 s). Postoperative complications included chyle leakage in 2 patients who were left kidney donors and oliguria in 1 patient who was a right kidney donor. All patients recovered with conservative care, and the mean hospital stay was 7.5 ± 1.7 days. The mean creatinine level was 0.7 ± 0.2 mg/dL before surgery, 1.1 ± 0.3 mg/dL at postoperative day (POD) 1, and 1.0 ± 0.2 mg/dL after discharge. The mean glomerular filtration rate was 97.9 ± 18.2 mL/min/1.73 m² before surgery, 60.7 ± 10.4 at POD 1, and 67.3 ± 11.1 after discharge. Operation time was not associated with patient body mass index and case number. No significant differences, other than postoperative complications, were found in the perioperative data for the side of kidney donation.

CONCLUSION

A skilled surgeon with experience in laparoscopic abdominal surgery (such as gastrectomy or colectomy) might safely perform hand-assisted donor nephrectomy. However, we could not identify a clear case number to complete the learning curve.

The genome of common long-arm octopus Octopus minor

Background

The common long-arm octopus (Octopus minor) is found in mudflats of subtidal zones and faces numerous environmental challenges. The ability to adapt its morphology and behavioral repertoire to diverse environmental conditions makes the species a promising model for understanding genomic adaptation and evolution in cephalopods.

Findings

The final genome assembly of O. minor is 5.09 Gb, with a contig N50 size of 197 kb and longest size of 3.027 Mb, from a total of 419 Gb raw reads generated using the Pacific Biosciences RS II platform. We identified 30,010 genes; 44.43% of the genome is composed of repeat elements. The genome-wide phylogenetic tree indicated the divergence time between O. minor and Octopus bimaculoides was estimated to be 43 million years ago based on single-copy orthologous genes. In total, 178 gene families are expanded in O. minor in the 14 bilaterian species.

Conclusions

We found that the O. minor genome was larger than that of closely related O. bimaculoides, and this difference could be explained by enlarged introns and recently diversified transposable elements. The high-quality O. minor genome assembly provides a valuable resource for understanding octopus genome evolution and the molecular basis of adaptations to mudflats.

Importance of Timed and Detailed Evaluation of Kidney Transplantation Candidates

Kidney transplant recipients are at increased risk of cardiovascular morbidity and malignant neoplasm, and meticulous evaluation of potential recipients is needed to minimize risks of complications after transplantation. The purpose of this study was to analyze the results of preoperative assessments and document the importance of timed and detailed examinations.

METHODS

Medical records of patients evaluated as kidney transplant candidates from January 2015 to September 2017 were retrospectively collected and analyzed.

RESULTS

Of the 216 patients evaluated during the study period, 135 (62.5%) were male, 112 (51.9%) had diabetes mellitus, 163 (75.5%) had hypertension, 31 (14.4%) had a cardiovascular event history, and 7 (3.2%) had previous history of malignant neoplasms. Mean (SD) patient age was 50.7 (10.8) years. All 216 recipient candidates underwent echocardiography. Mean (SD) ejection fraction was 57.8% (5.9%), and 48 candidates (22.2%) showed regional wall motional abnormality. Coronary angiography was performed on 81 candidates, and in 57 (70.4%) of these, coronary artery disease was detected. Malignant neoplasms were detected in 10 (4.6%) candidates. Kidney transplantation was performed on 55 candidates. One recipient died of Pneumocystis jirovecii pneumonia at 15 months after kidney transplant, but there was no death-censored graft failure, newly detected malignant neoplasm, or cardiovascular event over a mean (SD) follow-up duration of 15.5 (8.6) months.

CONCLUSION

Evaluation of kidney transplant candidates resulted in diagnoses of malignant neoplasms in 4.6% of patients and coronary artery disease in 26.4% of patients. The results of this study demonstrate candidates for kidney transplant should undergo detailed preoperative evaluation.

Abstract 1036: Establishment and molecular profiling of patient-derived xenograft models of Korean patients with primary lung cancer

Non-small cell lung cancer is the most prevalent type of lung cancer, and one of its subtype, squamous cell carcinoma (SCC) has the genomic complexity and high overall mutation due to the chemical carcinogens so that limited therapeutic options are available in comparison with adenocarcinoma. We constructed the patient-derived xenograft models of Korean primary lung cancer and then focused on the SCC PDX models for this particular ethnic subgroup. PDX models were established using the tissues of patients who underwent surgery as primary lung cancer at Samsung Medical Center during the period between October, 2014 and September, 2017. Briefly, tumor tissues from patients were subcutaneously engrafted and passaged two more times in NOD-scid-IL2Rγnull mouse. Models were selected whose tumor tissues showed the histopathology and genomics similar to the patient tumor tissues. Statistics of established PDX models are summarized in the Table.

Successful engraftmentSuccess rate (%)Primary lung cancer (n=532)19937.41Histology--Non-small cell lung cancer (n=524)19437.02- Adenocarcinoma (n=355)8925.07- Squamous cell carcinoma (n=135)8361.48- Other NSCLC (n=34)2264.71Small cell carcinoma (n=8)562.5Type of resection--Curative lung resection (n=458)17437.99Biopsy (n=74)2533.78Pathologic staging of ADC (n=355)--I (n=162)4024.69II (n=58)1525.86III (n=68)1927.94IV (n=67)1522.39

In order to help facilitate the target drug development of SCC, we concentrated on the construction of comprehensive database of fourteen SCC PDX models integrating the information from whole exome sequencing, whole transcriptome sequencing and clinical profiles. We detected mutations commonly reported in SCC genes including EGFR, BRAF, CUL3, DDR2, TP53, MYD88, PIK3CA, PTEN, MED12, SMAD4, AKT1, CTNNB1, IDH1, HRAS, KRAS, NRAS, SMARCA4, GNAS, CDKN2A, APC, CDKN2A and RRAS2. Specifically potential therapeutic targets such as EFGR, BRAF, CUL3, DDR2, TP53 and KRAS which are currently under several clinical trials were found in our SCC PDX models, and intriguingly one model had five mutations in CUL3, KRAS, TP53, and EGFR simultaneously. We expect these SCC PDX models could be valuable resource for target drug development and interpretation of their therapeutic responses in the pre- or co-clinical trials. Detailed information on pathology and genomic information on all PDX models will be available at DNA Link PDX webpage.

Citation Format: Tae Ho Kim, Hyunjin Heo, Soo Jung Lee, Seungje Lee, Eunjoo Hwang, Jinseon Lee, Heekyoung Lee, Kevin Koo, Hanna Lee, Seungjae Lee, Hwanseok Rhee, Jong Eun Lee, Jhingook Kim. Establishment and molecular profiling of patient-derived xenograft models of Korean patients with primary lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1036.

Case Report: Spontaneous Hemorrhage of a Rare Renal Tumor in the Native Kidney of a Renal Transplant Recipient

Renal cancers are some of the most common solid organ malignancies found during follow-up of patients who have undergone renal transplantation (RT). In this case report, we describe a life-threatening spontaneous hemorrhage of a rare subtype of renal cell carcinoma in the native kidney of a 27-year-old man, 4 years after RT. After fluid resuscitation and stabilization, the patient underwent emergent open radical nephrectomy with the final histopathology reporting T1bN0Mx mucinous tubular and spindle cell (MTSC) carcinoma. This case report highlights the need to consider an underlying malignancy in patients who presents with spontaneous hemorrhage of native kidneys after RT.

Structure, property, and function of sheepshead (Archosargus probatocephalus) teeth

OBJECTIVES

This paper studies A. probatocephalus teeth and investigates the mechanical properties and chemical composition of the enameloid and dentin.

DESIGN

Nanoindentation tests with a max load of 1000 μN and X-ray Energy Dispersive Spectroscopy (EDS) were performed along the diameter of the polished sample. Microstructural analysis of the dentin tubules was performed from SEM images.

RESULTS

From nanoindentation testing, the dentin of the sheepshead teeth has a nanoindentation hardness of 0.89 ± 0.21 (mean ± S.D.) GPa and a reduced Young's modulus of 23.29 ± 5.30 GPa. The enameloid of A. probatocephalus has a hardness of 4.36 ± 0.44 GPa and a mean reduced Young's modulus of 98.14 ± 6.91 GPa. Additionally, nanoindentation tests showed that the enameloid's hardness and modulus increased closer to the surface of the tooth. X-ray Energy Dispersive Spectroscopy (EDS) data further suggests that the gradient may be a result of the wt% fluoride within the enameloid, where an increase in fluoride results in an increase in reduced Young's modulus and hardness.

CONCLUSION

The microstructural characterization of the number density and area of the dentin tubules were used to address the porosity effect in the dentin to achieve the experimentally validated microhardness. The mechanical properties of the sheepshead teeth were also compared with previous nanoindentation tests from other aquatic species. The sheepshead teeth exhibit a greater reduced Young's modulus and hardness compared to shark and piranha teeth.

miR-204 downregulates EphB2 in aging mouse hippocampal neurons

Hippocampal synaptic function and plasticity deteriorate with age, often resulting in learning and memory deficits. As MicroRNAs (miRNAs) are important regulators of neuronal protein expression, we examined whether miRNAs may contribute to this age‐associated decline in hippocampal function. We first compared the small RNA transcriptome of hippocampal tissues from young and old mice. Among 269 hippocampal miRNAs, 80 were differentially expressed (≥ twofold) among the age groups. We focused on 36 miRNAs upregulated in the old mice compared with those in the young mice. The potential targets of these 36 miRNAs included 11 critical Eph/Ephrin synaptic signaling components. The expression levels of several genes in the Eph/Ephrin pathway, including EphB2, were significantly downregulated in the aged hippocampus. EphB2 is a known regulator of synaptic plasticity in hippocampal neurons, in part by regulating the surface expression of the NMDA receptor NR1 subunit. We found that EphB2 is a direct target of miR‐204 among miRNAs that were upregulated with age. The transfection of primary hippocampal neurons with a miR‐204 mimic suppressed both EphB2 mRNA and protein expression and reduced the surface expression of NR1. Transfection of miR‐204 also decreased the total expression of NR1. miR‐204 induces senescence‐like phenotype in fully matured neurons as evidenced by an increase in p16‐positive cells. We suggest that aging is accompanied by the upregulation of miR‐204 in the hippocampus, which downregulates EphB2 and results in reduced surface and total NR1 expression. This mechanism may contribute to age‐associated decline in hippocampal synaptic plasticity and the related cognitive functions.

The genome of Diuraphis noxia, a global aphid pest of small grains

BACKGROUND

The Russian wheat aphid, Diuraphis noxia Kurdjumov, is one of the most important pests of small grains throughout the temperate regions of the world. This phytotoxic aphid causes severe systemic damage symptoms in wheat, barley, and other small grains as a direct result of the salivary proteins it injects into the plant while feeding.

RESULTS

We sequenced and de novo assembled the genome of D. noxia Biotype 2, the strain most virulent to resistance genes in wheat. The assembled genomic scaffolds span 393 MB, equivalent to 93% of its 421 MB genome, and contains 19,097 genes. D. noxia has the most AT-rich insect genome sequenced to date (70.9%), with a bimodal CpG(O/E) distribution and a complete set of methylation related genes. The D. noxia genome displays a widespread, extensive reduction in the number of genes per ortholog group, including defensive, detoxification, chemosensory, and sugar transporter groups in comparison to the Acyrthosiphon pisum genome, including a 65% reduction in chemoreceptor genes. Thirty of 34 known D. noxia salivary genes were found in this assembly. These genes exhibited less homology with those salivary genes commonly expressed in insect saliva, such as glucose dehydrogenase and trehalase, yet greater conservation among genes that are expressed in D. noxia saliva but not detected in the saliva of other insects. Genes involved in insecticide activity and endosymbiont-derived genes were also found, as well as genes involved in virus transmission, although D. noxia is not a viral vector.

CONCLUSIONS

This genome is the second sequenced aphid genome, and the first of a phytotoxic insect. D. noxia's reduced gene content of may reflect the influence of phytotoxic feeding in shaping the D. noxia genome, and in turn in broadening its host range. The presence of methylation-related genes, including cytosine methylation, is consistent with other parthenogenetic and polyphenic insects. The D. noxia genome will provide an important contrast to the A. pisum genome and advance functional and comparative genomics of insects and other organisms.

Gene expression defines natural changes in mammalian lifespan

Mammals differ more than 100-fold in maximum lifespan, which can be altered in either direction during evolution, but the molecular basis for natural changes in longevity is not understood. Divergent evolution of mammals also led to extensive changes in gene expression within and between lineages. To understand the relationship between lifespan and variation in gene expression, we carried out RNA-seq-based gene expression analyses of liver, kidney, and brain of 33 diverse species of mammals. Our analysis uncovered parallel evolution of gene expression and lifespan, as well as the associated life-history traits, and identified the processes and pathways involved. These findings provide direct insights into how nature reversibly adjusts lifespan and other traits during adaptive radiation of lineages.

Brain somatic mutations in MTOR cause focal cortical dysplasia type II leading to intractable epilepsy

Focal cortical dysplasia type II (FCDII) is a sporadic developmental malformation of the cerebral cortex characterized by dysmorphic neurons, dyslamination and medically refractory epilepsy. It has been hypothesized that FCD is caused by somatic mutations in affected regions. Here, we used deep whole-exome sequencing (read depth, 412-668×) validated by site-specific amplicon sequencing (100-347,499×) in paired brain-blood DNA from four subjects with FCDII and uncovered a de novo brain somatic mutation, mechanistic target of rapamycin (MTOR) c.7280T>C (p.Leu2427Pro) in two subjects. Deep sequencing of the MTOR gene in an additional 73 subjects with FCDII using hybrid capture and PCR amplicon sequencing identified eight different somatic missense mutations found in multiple brain tissue samples of ten subjects. The identified mutations accounted for 15.6% of all subjects with FCDII studied (12 of 77). The identified mutations induced the hyperactivation of mTOR kinase. Focal cortical expression of mutant MTOR by in utero electroporation in mice was sufficient to disrupt neuronal migration and cause spontaneous seizures and cytomegalic neurons. Inhibition of mTOR with rapamycin suppressed cytomegalic neurons and epileptic seizures. This study provides, to our knowledge, the first evidence that brain somatic activating mutations in MTOR cause FCD and identifies mTOR as a treatment target for intractable epilepsy in FCD.

Geometric effects on stress wave propagation

The present study, through finite element simulations, shows the geometric effects of a bioinspired solid on pressure and impulse mitigation for an elastic, plastic, and viscoelastic material. Because of the bioinspired geometries, stress wave mitigation became apparent in a nonintuitive manner such that potential real-world applications in human protective gear designs are realizable. In nature, there are several toroidal designs that are employed for mitigating stress waves; examples include the hyoid bone on the back of a woodpecker's jaw that extends around the skull to its nose and a ram's horn. This study evaluates four different geometries with the same length and same initial cross-sectional diameter at the impact location in three-dimensional finite element analyses. The geometries in increasing complexity were the following: (1) a round cylinder, (2) a round cylinder that was tapered to a point, (3) a round cylinder that was spiraled in a two dimensional plane, and (4) a round cylinder that was tapered and spiraled in a two-dimensional plane. The results show that the tapered spiral geometry mitigated the greatest amount of pressure and impulse (approximately 98% mitigation) when compared to the cylinder regardless of material type (elastic, plastic, and viscoelastic) and regardless of input pressure signature. The specimen taper effectively mitigated the stress wave as a result of uniaxial deformational processes and an induced shear that arose from its geometry. Due to the decreasing cross-sectional area arising from the taper, the local uniaxial and shear stresses increased along the specimen length. The spiral induced even greater shear stresses that help mitigate the stress wave and also induced transverse displacements at the tip such that minimal wave reflections occurred. This phenomenon arose although only longitudinal waves were introduced as the initial boundary condition (BC). In nature, when shearing occurs within or between materials (friction), dissipation usually results helping the mitigation of the stress wave and is illustrated in this study with the taper and spiral geometries. The combined taper and spiral optimized stress wave mitigation in terms of the pressure and impulse; thus providing insight into the ram's horn design and woodpecker hyoid designs found in nature.

An epigenomic roadmap to induced pluripotency reveals DNA methylation as a reprogramming modulator

Reprogramming of somatic cells to induced pluripotent stem cells involves a dynamic rearrangement of the epigenetic landscape. To characterize this epigenomic roadmap, we have performed MethylC-seq, ChIP-seq (H3K4/K27/K36me3) and RNA-Seq on samples taken at several time points during murine secondary reprogramming as part of Project Grandiose. We find that DNA methylation gain during reprogramming occurs gradually, while loss is achieved only at the ESC-like state. Binding sites of activated factors exhibit focal demethylation during reprogramming, while ESC-like pluripotent cells are distinguished by extension of demethylation to the wider neighbourhood. We observed that genes with CpG-rich promoters demonstrate stable low methylation and strong engagement of histone marks, whereas genes with CpG-poor promoters are safeguarded by methylation. Such DNA methylation-driven control is the key to the regulation of ESC-pluripotency genes, including Dppa4, Dppa5a and Esrrb. These results reveal the crucial role that DNA methylation plays as an epigenetic switch driving somatic cells to pluripotency.

Differential expression of microRNAs and their target genes in non-small-cell lung cancer

MicroRNAs (miRNAs) are single‑stranded RNA species that constitute a class of non‑coding RNAs, and are emerging as key regulators of gene expression. Since each miRNA is capable of regulating multiple genes, miRNAs are attractive markers for studies of coordinated gene expression. In this study, we investigated miRNA expression profiling using a massively parallel sequencing technique to compare non‑small‑cell lung cancer (NSCLC) tissue and normal lung tissue. Lung cancer tissue and normal lung tissue were obtained from nine NSCLC patients. RNA isolated from these samples was processed using RNA sequencing (RNA Seq) and the HiSeq 2000 system. Differentially expressed miRNAs and mRNAs were analyzed using a t‑test. We selected target pairs that showed a negative correlation among significantly differentially expressed miRNAs and their putative target mRNAs using miRBase Targets. The differences in the expression levels of 222 miRNAs and 1,597 genes were statistically significant, as indicated by an absolute fold change ≥1.5 and P<0.05. miR‑577, miR‑301b, miR‑944, miR‑891a and miR‑615‑3p were generally upregulated, and miR‑338‑3p was generally downregulated. miRNA‑mRNA target pair analysis revealed that 49 miRNAs had 696 target mRNAs. There were significantly differentially expressed miRNAs and mRNAs between lung cancer and normal tissue. Further investigation of miRNAs and their target genes is warranted to better understand NSCLC.

Targeted exon sequencing fails to identify rare coding variants with large effect in rheumatoid arthritis

Introduction

Although it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea.

Methods

We analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722).

Results

After stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10⁻⁴). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA.

Conclusions

We were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-014-0447-7) contains supplementary material, which is available to authorized users.

RNA sequencing identifies novel markers of non-small cell lung cancer

INTRODUCTION

The development of reliable gene expression profiling technology increasingly impacts our understanding of lung cancer biology. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomes of non-small cell lung cancer (NSCLC) and normal lung tissues and to investigate expression in lung cancer tissues.

METHODS

We enrolled 88 male patients (mean age, 61.2 years) with NSCLC. RNA-Seq was performed on 88 pairs of NSCLC tumor tissue and non-tumor tissue from 54 patients with adenocarcinoma and 34 patients with squamous cell carcinoma. Immunohistochemistry was performed to validate differential candidate gene expression in a different NSCLC group.

RESULTS

RNA-Seq produced 25.41 × 10(6) (± 8.90 × 10(6)) reads in NSCLC tissues and 24.70×10(6) (± 4.70 × 10(6)) reads in normal lung tissues [mean (± standard deviation)]. Among the genes expressed in both tissues, 335 were upregulated and 728 were downregulated ≥ 2-fold (p < 0.001). Four upregulated genes - CBX3, GJB2, CRABP2, and DSP - not previously reported in lung cancer were studied further. Their altered expression was verified by immunohistochemistry in a different set of NSCLC tissues (n = 154). CBX3 was positive in 90.3% (139 cases) of the samples; GJB2, in 22.7% (35 cases); CRABP2, in 72.1% (111 cases); and DSP, in 17.5% (27 cases). The positive rate of CRABP2 was higher in adenocarcinoma than squamous cell carcinoma (p < 0.01).

CONCLUSIONS

CBX3 and CRABP2 expression was markedly increased in lung cancer tissues and especially CRABP2 may be promising candidate genes in lung adenocarcinoma.

Comprehensive analysis to improve the validation rate for single nucleotide variants detected by next-generation sequencing

Next-generation sequencing (NGS) has enabled the high-throughput discovery of germline and somatic mutations. However, NGS-based variant detection is still prone to errors, resulting in inaccurate variant calls. Here, we categorized the variants detected by NGS according to total read depth (TD) and SNP quality (SNPQ), and performed Sanger sequencing with 348 selected non-synonymous single nucleotide variants (SNVs) for validation. Using the SAMtools and GATK algorithms, the validation rate was positively correlated with SNPQ but showed no correlation with TD. In addition, common variants called by both programs had a higher validation rate than caller-specific variants. We further examined several parameters to improve the validation rate, and found that strand bias (SB) was a key parameter. SB in NGS data showed a strong difference between the variants passing validation and those that failed validation, showing a validation rate of more than 92% (filtering cutoff value: alternate allele forward [AF]≥20 and AF<80 in SAMtools, SB<–10 in GATK). Moreover, the validation rate increased significantly (up to 97–99%) when the variant was filtered together with the suggested values of mapping quality (MQ), SNPQ and SB. This detailed and systematic study provides comprehensive recommendations for improving validation rates, saving time and lowering cost in NGS analyses.

Abstract 1113: The comprehensive mutational profiling of lung adenocarcinoma through large-scale RNA sequencing analysis

All cancers harbor molecular alterations in their genomes. The transcriptional consequences of these somatic mutations have not yet been comprehensively explored in lung cancer. Here we present the first large scale RNA sequencing study of lung adenocarcinoma, demonstrating its power to identify somatic point mutations as well as transcriptional variants such as gene fusions, alternative splicing events and expression outliers. Our results reveal the genetic basis of 200 lung adenocarcinomas in Koreans including deep characterization of 87 surgical specimens by transcriptome sequencing. We identified driver somatic mutations in cancer genes including EGFR, KRAS, NRAS, BRAF, PIK3CA, MET and CTNNB1. Candidates for novel driver mutations were also identified in genes newly implicated in lung adenocarcinoma such as LMTK2, ARID1A, NOTCH2 and SMARCA4. We found 45 fusion genes, 8 of which were chimeric tyrosine kinases involving ALK, RET, ROS1, FGFR2, AXL and PDGFRA. Among 17 recurrent alternative splicing events, we identified exon 14 skipping in the proto-oncogene MET as highly likely to be a cancer driver. The number of somatic mutations and expression outliers varied markedly between individual cancers and was strongly correlated with smoking history of patients. We identified genomic blocks within which gene expression levels were consistently increased or decreased that could be explained by copy number alterations in samples. We also found an association between lymph node metastasis and somatic mutations in TP53. These findings broaden our understanding of lung adenocarcinoma and may also lead to new diagnostic and therapeutic approaches.

Citation Format: Won-Chul Lee, Young Seok Ju, Jong-Yeon Shin, June Koo Lee, Thomas Bleazard, Jun-Ho Lee, Jung-Oh Kim, Saet-Byeol Yu, Jihye Kim, Hwanseok Rhee, Se-Hoon Lee, Jong-Il Kim, Jin-Hyoung Kang, Young Tae Kim, Jeong-Sun Seo. The comprehensive mutational profiling of lung adenocarcinoma through large-scale RNA sequencing analysis. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1113. doi:10.1158/1538-7445.AM2013-1113

The transcriptional landscape and mutational profile of lung adenocarcinoma

All cancers harbor molecular alterations in their genomes. The transcriptional consequences of these somatic mutations have not yet been comprehensively explored in lung cancer. Here we present the first large scale RNA sequencing study of lung adenocarcinoma, demonstrating its power to identify somatic point mutations as well as transcriptional variants such as gene fusions, alternative splicing events, and expression outliers. Our results reveal the genetic basis of 200 lung adenocarcinomas in Koreans including deep characterization of 87 surgical specimens by transcriptome sequencing. We identified driver somatic mutations in cancer genes including EGFR, KRAS, NRAS, BRAF, PIK3CA, MET, and CTNNB1. Candidates for novel driver mutations were also identified in genes newly implicated in lung adenocarcinoma such as LMTK2, ARID1A, NOTCH2, and SMARCA4. We found 45 fusion genes, eight of which were chimeric tyrosine kinases involving ALK, RET, ROS1, FGFR2, AXL, and PDGFRA. Among 17 recurrent alternative splicing events, we identified exon 14 skipping in the proto-oncogene MET as highly likely to be a cancer driver. The number of somatic mutations and expression outliers varied markedly between individual cancers and was strongly correlated with smoking history of patients. We identified genomic blocks within which gene expression levels were consistently increased or decreased that could be explained by copy number alterations in samples. We also found an association between lymph node metastasis and somatic mutations in TP53. These findings broaden our understanding of lung adenocarcinoma and may also lead to new diagnostic and therapeutic approaches.

Exome sequencing is an efficient tool for genetic screening of Charcot-Marie-Tooth disease

Charcot-Marie-Tooth disease (CMT) is one of the most common inherited neuropathies and is a genetically and clinically heterogeneous disorder with variable inheritance modes. As several molecules have been reported to have therapeutic effects on CMT, depending on the underlying genetic causes, exact genetic diagnostics have become very important for executing personalized therapy. Whole-exome sequencing has recently been introduced as an available method to identify rare or novel genetic defects from genetic disorders. Particularly, CMT is a model disease to apply exome sequencing because more than 50 genes (loci) are involved in its development with weak genotype-phenotype correlation. This study performed the exome sequencing in 25 unrelated CMT patients who revealed neither 17p12 duplication/deletion nor several major CMT genes. This study identified eight causative heterozygous mutations (32%). This detection rate seems rather high because each sample was tested before the study for major genetic causes. Therefore, this study suggests that the exome sequencing can be a highly exact, rapid, and economical molecular diagnostic tool for CMT patients who are tested for major genetic causes.

Differentiation of early hepatocellular carcinoma from benign hepatocellular nodules on gadoxetic acid-enhanced MRI

OBJECTIVE

To test new diagnostic criteria for the discrimination of early hepatocellular carcinoma (HCC) from benign hepatocellular nodules on gadoxetic acid-enhanced MRI (Gd-EOB-MRI).

METHODS

We retrospectively analysed 34 patients with 29 surgically diagnosed early HCCs and 31 surgically diagnosed benign hepatocellular nodules. Two radiologists reviewed Gd-EOB-MRI, including diffusion-weighted imaging (DWI), and the signal intensity at each sequence, presence of arterial enhancement and washout were recorded. We composed new diagnostic criteria based on the lesion size and MRI findings, and then the diagnostic performance was compared with that of conventional imaging criteria with logistic regression and a generalised estimating equation method.

RESULTS

A size cut-off value (≥1.5 cm diameter) and MRI findings of T(1) hypointensity, T(2) hyperintensity, DWI hyperintensity on both low and high b-value images (b=50 and 800 s mm(-2), respectively), arterial enhancement, late washout and hepatobiliary hypointensity were selected as the diagnostic criteria. When lesions were considered malignant if they satisfied three or more of the above criteria, the sensitivity was significantly higher than when making a diagnosis based on arterial enhancement and washout alone (58.6% vs 13.8%, respectively; p=0.0002), while the specificity was 100.0% for both criteria.

CONCLUSION

Our new diagnostic criteria on Gd-EOB-MRI may help to improve the discrimination of early HCC from benign hepatocellular nodules.

FX: an RNA-Seq analysis tool on the cloud

FX is an RNA-Seq analysis tool, which runs in parallel on cloud computing infrastructure, for the estimation of gene expression levels and genomic variant calling. In the mapping of short RNA-Seq reads, FX uses a transcriptome-based reference primarily, generated from ~160 000 mRNA sequences from RefSeq, UCSC and Ensembl databases. This approach reduces the misalignment of reads originating from splicing junctions. Unmapped reads not aligned on known transcripts are then mapped on the human genome reference. FX allows analysis of RNA-Seq data on cloud computing infrastructures, supporting access through a user-friendly web interface.

AVAILABILITY

FX is freely available on the web at (http://fx.gmi.ac.kr), and can be installed on local Hadoop clusters. Guidance for the installation and operation of FX can be found under the 'Documentation' menu on the website.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exome sequencing and subsequent association studies identify five amino acid-altering variants influencing human height

Height is a highly heritable trait that involves multiple genetic loci. To identify causal variants that influence stature, we sequenced whole exomes of four children with idiopathic short stature. Ninety-five nonsynonymous single-nucleotide polymorphisms (nsSNPs) were selected as potential candidate variants. We performed association analysis in 740 cohort individuals and identified 11 nsSNPs in 10 loci (DIS3L2, ZBTB38, FAM154A, PTCH1, TSSC4, KIF18A, GPR133, ACAN, FAM59A, and NINL) associated with adult height (P < 0.05), including five novel loci. Of these, two nsSNPs (TSSC4 and KIF18A loci) were significant at P < 0.05 in the replication study (n = 1,000) and five (ZBTB38, FAM154A, TSSC4, KIF18A, and FAM59A loci) were significant at P < 0.01 in the combined analysis (n = 1,740). Together, the five nsSNPs accounted for approximately 2.5% of the height variation. This study demonstrated the utility of next-generation sequencing in identifying genetic variants and loci associated with complex traits.

Gene expression changes in patient-matched gastric normal mucosa, adenomas, and carcinomas

A subset of gastric carcinomas shows histologic evidence of a multistep process, progressing from gastric adenoma to gastric carcinoma. We examined gene expression changes during the gastric adenoma-carcinoma sequence in 26 snap-frozen samples (normal mucosa, adenoma, and carcinoma samples from eight patients and two additional carcinomas) by oligonucleotide microarray. Unsupervised hierarchical clustering analysis demonstrated differential gene expression between gastric normal mucosa, adenomas and carcinomas. We identified 319 and 422 genes differentially regulated in adenoma and carcinoma, respectively, relative to normal mucosa, using a combination of Welch's t-test and fold-change analysis. Applying a combination of robust multi-category support vector machines to the data, reveal that 39 and 21 genes were gradually up- and down-regulated, respectively, in succession in normal mucosa, adenoma, and carcinoma samples. We validated gene expression levels of four genes: hydroxyprostaglandin dehydrogenase 15 (HPGD), follistatin-like 1, trefoil factor 1 (TTF1) and trefoil factor 2 (TFF2) by RT-PCR and found direct correlation with microarray results. The expressions of the TFF2 and HPGD genes were further evaluated by immunohistochemistry in 103 adenomas and 70 carcinomas; expression of both proteins was decreased in these tissues. The progressive alteration in gene expression in the transition from normal mucosa to carcinoma suggests that these changes may play critical roles in gastric carcinogenesis.