Accumulation of potential driver genes with genomic alterations predicts survival of high-risk neuroblastoma patients

Background

Neuroblastoma is the most common pediatric malignancy with heterogeneous clinical behaviors, ranging from spontaneous regression to aggressive progression. Many studies have identified aberrations related to the pathogenesis and prognosis, broadly classifying neuroblastoma patients into high- and low-risk groups, but predicting tumor progression and clinical management of high-risk patients remains a big challenge.

Results

We integrate gene-level expression, array-based comparative genomic hybridization and functional gene-interaction network of 145 neuroblastoma patients to detect potential driver genes. The drivers are summarized into a driver-gene score (DGscore) for each patient, and we then validate its clinical relevance in terms of association with patient survival. Focusing on a subset of 48 clinically defined high-risk patients, we identify 193 recurrent regions of copy number alterations (CNAs), resulting in 274 altered genes whose copy-number gain or loss have parallel impact on the gene expression. Using a network enrichment analysis, we detect four common driver genes, ERCC6, HECTD2, KIAA1279, EMX2, and 66 patient-specific driver genes. Patients with high DGscore, thus carrying more copy-number-altered genes with correspondingly up- or down-regulated expression and functional implications, have worse survival than those with low DGscore (P?=?0.006). Furthermore, Cox proportional-hazards regression analysis shows that, adjusted for age, tumor stage and MYCN amplification, DGscore is the only significant prognostic factor for high-risk neuroblastoma patients (P?=?0.008).

Conclusions

Integration of genomic copy number alteration, expression and functional interaction-network data reveals clinically relevant and prognostic putative driver genes in high-risk neuroblastoma patients. The identified putative drivers are potential drug targets for individualized therapy.

Reviewers

This article was reviewed by Armand Valsesia, Susmita Datta and Aleksandra Gruca.

     

Analysis of the conservation of synteny between Fugu and human chromosome 12

Background

The pufferfish Fugu rubripes (Fugu) with its compact genome is increasingly recognized as an important vertebrate model for comparative genomic studies. In particular, large regions of conserved synteny between human and Fugu genomes indicate its utility to identify disease-causing genes. The human chromosome 12p12 is frequently deleted in various hematological malignancies and solid tumors, but the actual tumor suppressor gene remains unidentified.

Results

We investigated approximately 200 kb of the genomic region surrounding the ETV6 locus in Fugu (fETV6) in order to find conserved functional features, such as genes or regulatory regions, that could give insight into the nature of the genes targeted by deletions in human cancer cells. Seven genes were identified near the fETV6 locus. We found that the synteny with human chromosome 12 was conserved, but extensive genomic rearrangements occurred between the Fugu and human ETV6 loci.

Conclusion

This comparative analysis led to the identification of previously uncharacterized genes in the human genome and some potentially important regulatory sequences as well. This is a good indication that the analysis of the compact Fugu genome will be valuable to identify functional features that have been conserved throughout the evolution of vertebrates.

     

Genome rearrangements in <Emphasis Type="Italic">Escherichia coli</Emphasis> during de novo acquisition of resistance to a single antibiotic or two antibiotics successively

Background

The ability of bacteria to acquire resistance to antibiotics relies to a large extent on their capacity for genome modification. Prokaryotic genomes are highly plastic and can utilize horizontal gene transfer, point mutations, and gene deletions or amplifications to realize genome expansion and rearrangements. The contribution of point mutations to de novo acquisition of antibiotic resistance is well-established. In this study, the internal genome rearrangement of Escherichia coli during to de novo acquisition of antibiotic resistance was investigated using whole-genome sequencing.

Results

Cells were made resistant to one of the four antibiotics and subsequently to one of the three remaining. This way the initial genetic rearrangements could be documented together with the effects of an altered genetic background on subsequent development of resistance. A DNA fragment including ampC was amplified by a factor sometimes exceeding 100 as a result of exposure to amoxicillin. Excision of prophage e14 was observed in many samples with a double exposure history, but not in cells exposed to a single antibiotic, indicating that the activation of the SOS stress response alone, normally the trigger for excision, was not sufficient to cause excision of prophage e14. Partial deletion of clpS and clpA occurred in strains exposed to enrofloxacin and tetracycline. Other deletions were observed in some strains, but not in replicates with the exact same exposure history. Various insertion sequence transpositions correlated with exposure to specific antibiotics.

Conclusions

Many of the genome rearrangements have not been reported before to occur during resistance development. The observed correlation between genome rearrangements and specific antibiotic pressure, as well as their presence in independent replicates indicates that these events do not occur randomly. Taken together, the observed genome rearrangements illustrate the plasticity of the E. coli genome when exposed to antibiotic stress.

     

Intraosseous inflammatory myofibroblastic tumor of the mandible with a novel <Emphasis Type="Italic">ATIC-ALK</Emphasis> fusion mutation: a case report

Background

Inflammatory myofibroblastic tumor (IMT) is a rare low-grade malignant neoplasm with a predilection for children and young adults, and typically arises in the lung, abdominopelvic region, and retroperitoneum. IMTs in the maxillofacial region are extreme rare. Approximately 50% of IMT harbor rearrangements of the anaplastic lymphoma kinase (ALK) gene at 2p23 with various fusion partners.

Case presentation

We herein report a case of intraosseous IMT of the mandible with a novel ATIC-ALK fusion. Tooth 43 did not erupt after the loss of tooth 83 in an 11-year-old girl with no previous history of trauma. Panoramic tomography showed a unilocular radiolucent lesion in the right anterior mandible resorbing the root of tooth 42 and the medial side of the root of tooth 44. Computed tomography revealed a well- circumscribed 3-cm osteolytic lesion of the right anterior mandible eroding the buccal cortical plate. The entire lesion was curetted out. A histopathological examination revealed the proliferation of plump spindle cells with a storiform architecture and lymphocytes scattered around spindle cells. The spindle cells showed diffuse cytoplasmic staining for ALK by immunohistochemistry. A fluorescence in situ hybridization analysis revealed the translocation of a part of the ALK gene locus at chromosome 2p23. A rapid amplification of cDNA ends analysis confirmed the rearrangement of ALK and identified ATIC as a partner of this ALK fusion mutant.

Conclusion

To the best of our knowledge, this is the first case of intraosseous IMT of the mandible with a novel ATIC-ALK fusion. We also herein reviewed similar tumors reported in the literature.

     

Whole-exon sequencing of human myeloma cell lines shows mutations related to myeloma patients at relapse with major hits in the DNA regulation and repair pathways

Background

Human myeloma cell lines (HMCLs) are widely used for their representation of primary myeloma cells because they cover patient diversity, although not fully. Their genetic background is mostly undiscovered, and no comprehensive study has ever been conducted in order to reveal those details.

Methods

We performed whole-exon sequencing of 33 HMCLs, which were established over the last 50?years in 12 laboratories. Gene expression profiling and drug testing for the 33 HMCLs are also provided and correlated to exon-sequencing findings.

Results

Missense mutations were the most frequent hits in genes (92%). HMCLs harbored between 307 and 916 mutations per sample, with TP53 being the most mutated gene (67%). Recurrent bi-allelic losses were found in genes involved in cell cycle regulation (RB1, CDKN2C), the NFκB pathway (TRAF3, BIRC2), and the p53 pathway (TP53, CDKN2A). Frequency of mutations/deletions in HMCLs were either similar to that of patients (e.g., DIS3, PRDM1, KRAS) or highly increased (e.g., TP53, CDKN2C, NRAS, PRKD2). MAPK was the most altered pathway (82% of HMCLs), mainly by RAS mutants. Surprisingly, HMCLs displayed alterations in epigenetic (73%) and Fanconi anemia (54%) and few alterations in apoptotic machinery. We further identified mutually exclusive and associated mutations/deletions in genes involved in the MAPK and p53 pathways as well as in chromatin regulator/modifier genes. Finally, by combining the gene expression profile, gene mutation, gene deletion, and drug response, we demonstrated that several targeted drugs overcome or bypass some mutations.

Conclusions

With this work, we retrieved genomic alterations of HMCLs, highlighting that they display numerous and unprecedented abnormalities, especially in DNA regulation and repair pathways. Furthermore, we demonstrate that HMCLs are a reliable model for drug screening for refractory patients at diagnosis or at relapse.

     

Association of a single nucleotide polymorphism in growth differentiate factor 5 with congenital dysplasia of the hip: a case-control study

Introduction

Congenital dysplasia of the hip is an abnormal seating of the femoral head in the acetabulum, mainly caused by shallow acetabulum and lax joint capsule. Genetic factors play a considerable role in the pathogenesis of congenital dysplasia of the hip. The gene growth differentiate factor 5 (GDF5) has been implicated in skeletal development and joint morphogenesis in humans and mice. A functional single nucleotide polymorphism (SNP) in the 5'-untranslated region of GDF5 (rs143383) was reported to be associated with osteoarthritis susceptibility. As a key regulator in morphogenesis of skeletal components and soft tissues in and around the joints, GDF5 may be involved in the aetiology and pathogenesis of congenital dysplasia of the hip. Our objective is to evaluate if the GDF5 SNP is associated with congenital dysplasia of the hip in people of Han Chinese origin.

Methods

The GDF5 SNP was genotyped in 338 children with congenital dysplasia of the hip and 622 control subjects.

Results

The SNP was significantly associated with congenital dysplasia of the hip (p = 0.0037; odds ration (OR) = 1.40; 95% confidence interval (CI) = 1.11 to 1.75). A significant difference was detected in female samples when stratified by gender (p = 0.0053; OR = 1.46; 95% CI = 1.21 to 1.91), and in hip dislocation when stratified by severity (p = 0.0078; OR = 1.43; 95% CI = 1.11 to 1.85).

Conclusions

Our results indicate that GDF5 is important in the aetiology of congenital dysplasia of the hip. To the authors' knowledge this is the first time that a definite association with the congenital dysplasia of the hip susceptibility has been detected.

     

Rapid identification of an adult plant stripe rust resistance gene in hexaploid wheat by high-throughput SNP array genotyping of pooled extremes

Key message

High-throughput SNP array analysis of pooled extreme phenotypes in a segregating population by KASP marker genotyping permitted rapid, cost-effective location of a stripe rust resistance QTL in wheat.

Abstract

German wheat cultivar “Friedrichswerther” has exhibited high levels of adult plant resistance (APR) to stripe rust in field environments for many years. F_2:3 lines and F₆ recombinant inbred line (RILs) populations derived from a cross between Friedrichswerther and susceptible landrace Mingxian 169 were evaluated in the field in 2013, 2016 and 2017. Illumina 90K iSelect SNP arrays were used to genotype bulked extreme pools and parents; 286 of 1135 polymorphic SNPs were identified on chromosome 6B. Kompetitive Allele-Specific PCR (KASP) markers were used to verify the chromosome region associated with the resistance locus. A linkage map was constructed with 18 KASP-SNP markers, and a major effect QTL was identified within a 1.4 cM interval flanked by KASP markers IWB71602 and IWB55937 in the region 6BL3-0-0.36. The QTL, named QYr.nwafu-6BL, was stable across environments, and explained average 54.4 and 47.8% of the total phenotypic variation in F_2:3 lines and F₆ RILs, respectively. On the basis of marker genotypes, pedigree analysis and relative genetic distance QYr.nwafu-6BL is likely to be a new APR QTL. Combined high-throughput SNP array genotyping of pooled extremes and validation by KASP assays lowers sequencing costs compared to genome-wide association studies with SNP arrays, and more importantly, permits rapid isolation of major effect QTL in hexaploid wheat as well as improving accuracy of mapping in the QTL region. QYr.nwafu-6BL with flanking KASP markers developed and verified in a subset of 236 diverse lines can be used in marker-assisted selection to improve stripe rust resistance in breeding programs.

     

ClpXP affects the cell metabolism of <Emphasis Type="Italic">Salmonella typhimurium</Emphasis> partially in an RpoS-dependent manner

Introduction

ClpXP protease is an important proteolytic system in Salmonella enterica serovar typhimurium (S. typhimurium). Inactivation of ClpXP by deletion of clpP resulted in overproduction of RpoS and a growth defect phenotype. Only one report has indicated that deleting rpoS can restore the growth of a S. typhimurium clpP mutant to the wild-type level. Whether overproduction of RpoS is responsible for the growth deficiency resulting from clpP disruption and how ClpXP affects the cell metabolism of S. typhimurium remain to be elucidated.

Objectives

The aim of this study is to investigate the effect of ClpXP on cell metabolism of S. typhimurium and explore the possible co-effect of RpoS associated with ClpXP in cell metabolism.

Method

We constructed a clpP rpoS double deletion mutant TT-19 (ΔclpP ΔrpoS TT-1) using a two-step phage transduction technique. We then compared the metabolite fingerprints of Salmonella rpoS deletion mutant TT-14 (ΔrpoS TT-1), clpP deletion mutant TT-16 (ΔclpP TT-1), and clpP rpoS double deletion mutant TT-19 (ΔclpP ΔrpoS TT-1) with those of the wild-type strain TT-1 by using gas chromatography coupled with mass spectrometry (GC–MS).

Results

Deletion of rpoS recovered only a part of the growth of Salmonella clpP mutant. Further metabolome analysis indicated that clpP disruption changed the levels of 16 extra- and 19 intracellular substances, while the extracellular concentrations of 4 compounds (serine, l-5-oxoproline, l-glutamic acid, and l-tryptophan) and intracellular concentrations of 10 compounds (l-isoleucine, glycine, serine, l-methionine, l-phenylalanine, malic acid, citric acid, urea, putrescine, and 6-hydroxypurine) returned to their wild-type levels when rpoS was also deleted.

Conclusion

ClpXP affects the cell metabolism of S. typhimurium partially in an RpoS-dependent manner.

     

A SNP in the 5′ flanking region of the myostatin-1b gene is associated with harvest traits in Atlantic salmon (Salmo salar)

Background

Myostatin (MSTN) belongs to the transforming growth factor-β superfamily and is a potent negative regulator of skeletal muscle development and growth in mammals. Most teleost fish possess two MSTN paralogues. However, as a consequence of a recent whole genome-duplication event, salmonids have four: MSTN-1 (?1a and -1b) and MSTN-2 (?2a and -2b). Evidence suggests that teleost MSTN plays a role in the regulation of muscle growth. In the current study, the MSTN-1b gene was re-sequenced and screened for SNP markers in a commercial population of Atlantic salmon. After genotyping 4,800 progeny for the discovered SNPs, we investigated their association with eight harvest traits - four body-weight traits, two ratios of weight traits, flesh colour and fat percentage - using a mixed model association analysis.

Results

Three novel SNPs were discovered in the MSTN-1b gene of Atlantic salmon. One of the SNPs, located within the 5′ flanking region (g.1086C?>?T), had a significant association with harvest traits (p?<?0.05), specifically for: Harvest Weight (kg), Gutted Weight (kg), Deheaded Weight (kg) and Fillet Weight (kg). The haplotype-based association analysis was consistent with this result because the two haplotypes that showed a significant association with body-weight traits, hap4 and hap5 (p?<?0.05 and p?<?0.01, respectively), differ by a single substitution at the g.1086C?>?T locus. The alleles at g.1086C?>?T act in an additive manner and explain a small percentage of the genetic variation of these phenotypes.

Conclusions

The association analysis revealed that g.1086C?>?T had a significant association with all body-weight traits under study. Although the SNP explains a small percentage of the variance, our results indicate that a variation in the 5′ flanking region of the myostatin gene is associated with the genetic regulation of growth in Atlantic salmon.

     

Nonrecurrent <Emphasis Type="Italic">PMP22</Emphasis>-<Emphasis Type="Italic">RAI1</Emphasis> contiguous gene deletions arise from replication-based mechanisms and result in Smith–Magenis syndrome with evident peripheral neuropathy

Hereditary neuropathy with liability to pressure palsies (HNPP) and Smith–Magenis syndrome (SMS) are genomic disorders associated with deletion copy number variants involving chromosome 17p12 and 17p11.2, respectively. Nonallelic homologous recombination (NAHR)-mediated recurrent deletions are responsible for the majority of HNPP and SMS cases; the rearrangement products encompass the key dosage-sensitive genes PMP22 and RAI1, respectively, and result in haploinsufficiency for these genes. Less frequently, nonrecurrent genomic rearrangements occur at this locus. Contiguous gene duplications encompassing both PMP22 and RAI1, i.e., PMP22-RAI1 duplications, have been investigated, and replication-based mechanisms rather than NAHR have been proposed for these rearrangements. In the current study, we report molecular and clinical characterizations of six subjects with the reciprocal phenomenon of deletions spanning both genes, i.e., PMP22-RAI1 deletions. Molecular studies utilizing high-resolution array comparative genomic hybridization and breakpoint junction sequencing identified mutational signatures that were suggestive of replication-based mechanisms. Systematic clinical studies revealed features consistent with SMS, including features of intellectual disability, speech and gross motor delays, behavioral problems and ocular abnormalities. Five out of six subjects presented clinical signs and/or objective electrophysiologic studies of peripheral neuropathy. Clinical profiling may improve the clinical management of this unique group of subjects, as the peripheral neuropathy can be more severe or of earlier onset as compared to SMS patients having the common recurrent deletion. Moreover, the current study, in combination with the previous report of PMP22-RAI1 duplications, contributes to the understanding of rare complex phenotypes involving multiple dosage-sensitive genes from a genetic mechanistic standpoint.     

Production of β-carotene by expressing a heterologous multifunctional carotene synthase in <Emphasis Type="Italic">Yarrowia lipolytica</Emphasis>

Objectives

To obtain functional expression of a heterologous multifunctional carotene synthase containing phytoene synthase, phytoene dehydrogenase, and lycopene β-cyclase activities encoded by carS from Schizochytrium sp. in order to allow Yarrowia lipolytica to produce β-carotene.

Results

To increase the integration efficiency of a 3.8 kb carS under the control of P _GPD promoter with a 2 kb selection marker, ura3, along with a geranylgeranyl diphosphate synthase (GGS1) expression cassette (~10 kb in total), was inserted into the Y. lipolytica chromosome, and the DNA assembler method was combined with double chromosomal deletions of ku70 and ku80. This method resulted in a 13.4-fold increase in integration efficiency compared with the original method, reaching 63% (10/16). The resulting recombinant Y. lipolytica produced 0.41 mg β-carotene per g dry cell weight, while the wild type did not produce any indicating the functionality of the multifunctional carotene synthase in Y. lipolytica.

Conclusion

Expression of GGS1 and a multifunctional carotene synthase from Schizochytrium sp. in Y. lipolytica led to β-carotene production. DNA assembler efficiency was greatly increased by the deletion of ku70 and ku80, which resulted in decreased in vivo nonhomologous end-joining (NHEJ) in Y. lipolytica.

     

Genetic diagnosis of neurofibromatosis type 1: targeted next- generation sequencing with Multiple Ligation-Dependent Probe Amplification analysis

Background

Neurofibromatosis type 1 (NF1) is a dominantly inherited tumor predisposition syndrome that targets the peripheral nervous system. It is caused by mutations of the NF1 gene which serve as a negative regulator of the cellular Ras/MAPK (mitogen-activated protein kinases) signaling pathway. Owing to the complexity in some parts of clinical diagnoses and the need for better understanding of its molecular relationships, a genetic characterization of this disorder will be helpful in the clinical setting.

Methods

In this study, we present a customized targeted gene panel of NF1/KRAS/BRAF/p53 and SPRED1 genes combined with Multiple Ligation-Dependent Probe Amplification analysis for the NF1 mutation screening in a cohort of patients clinically suspected as NF1.

Results

In this study, we identified 73 NF1 mutations and two BRAF novel variants from 100 NF1 patients who were suspected as having NF1. These genetic alterations are heterogeneous and distribute in a complicated way without clustering in either cysteine–serine-rich domain or within the GAP-related domain. We also detected fifteen multi-exon deletions within the NF1 gene by MLPA Analysis.

Conclusions

Our results suggested that a genetic screening using a NGS panel with high coverage of Ras–signaling components combined with Multiple Ligation-Dependent Probe Amplification analysis will enable differential diagnosis of patients with overlapping clinical features.

     

Genetic sex determination of mice by simplex PCR

Background

Investigating fetal development in mice necessitates the determination of fetal sex. However, whilst the sex of adult and juvenile mice can be readily distinguished from anogenital distance, the sex of fetal and neonatal mice cannot be identified visually. Instead, genetic sex must be determined by PCR amplification of X chromosome genes with divergent Y chromosome gametologs. Existing simplex PCR methods are confounded by small size differences between amplicons, amplification of unexpected products, and biased amplification of the shorter amplicon.

Results

Primers were designed flanking an 84 bp deletion of the X-linked Rbm31x gene relative to its Y-linked gametolog Rbm31y. A single product was amplified from XX samples, with two products amplified from XY samples. Amplicons were resolved by gel electrophoresis for 20 min, with unbiased amplification of both products observed in XY samples.

Conclusion

This method achieves rapid and unequivocal genetic sex determination of mice in low volume PCR reactions, reducing reagent usage and simultaneously eliminating shortcomings of previous methods.

     

CRISPR-Cas9 mediated gene deletions in lager yeast <Emphasis Type="Italic">Saccharomyces pastorianus</Emphasis>

Background

The ease of use of CRISPR-Cas9 reprogramming, its high efficacy, and its multiplexing capabilities have brought this technology at the forefront of genome editing techniques. Saccharomyces pastorianus is an aneuploid interspecific hybrid of Saccharomyces cerevisiae and Saccharomyces eubayanus that has been domesticated for centuries and is used for the industrial fermentation of lager beer. For yet uncharacterised reasons, this hybrid yeast is far more resilient to genetic alteration than its ancestor S. cerevisiae.

Results

This study reports a new CRISPR-Cas9 method for accurate gene deletion in S. pastorianus. This method combined the Streptococcus pyogenes cas9 gene expressed from either a chromosomal locus or from a mobile genetic element in combination with a plasmid-borne gRNA expression cassette. While the well-established gRNA expression system using the RNA polymerase III dependent SNR52 promoter failed, expression of a gRNA flanked with Hammerhead and Hepatitis Delta Virus ribozymes using the RNA polymerase II dependent TDH3 promoter successfully led to accurate deletion of all four alleles of the SeILV6 gene in strain CBS1483. Furthermore the expression of two ribozyme-flanked gRNAs separated by a 10-bp linker in a polycistronic array successfully led to the simultaneous deletion of SeATF1 and SeATF2, genes located on two separate chromosomes. The expression of this array resulted in the precise deletion of all five and four alleles mediated by homologous recombination in the strains CBS1483 and Weihenstephan 34/70 respectively, demonstrating the multiplexing abilities of this gRNA expression design.

Conclusions

These results firmly established that CRISPR-Cas9 significantly facilitates and accelerates genome editing in S. pastorianus.

     

A 55 K SNP array-based genetic map and its utilization in QTL mapping for productive tiller number in common wheat

Key message

A high-density genetic map constructed with a wheat 55 K SNP array was highly consistent with the physical map of this species and it facilitated the identification of a novel major QTL for productive tiller number.

Abstract

Productive tiller number (PTN) plays a key role in wheat grain yield. In this study, a recombinant inbred line population with 199 lines derived from a cross between ‘20828’ and ‘Chuannong16’ was used to construct a high-density genetic map using wheat 55 K single nucleotide polymorphism (SNP) array. The constructed genetic map contains 12,109 SNP markers spanning 3021.04 cM across the 21 wheat chromosomes. The orders of the genetic and physical positions of these markers are generally in agreement, and they also match well with those based on the 660 K SNP array from which the one used in this study was derived. The ratios of SNPs located in each of the wheat deletion bins were similar among the wheat 9 K, 55 K, 90 K, 660 K and 820 K SNP arrays. Based on the constructed maps, a novel major quantitative trait locus QPtn.sau-4B for PTN was detected across multi-environments in a 0.55 cM interval on 4B and it explained 17.23–45.46% of the phenotypic variance. Twenty common genes in the physical interval between the flanking markers were identified on chromosome 4B of ‘Chinese Spring’ and wild emmer. These results indicate that wheat 55 K SNP array could be an ideal tool in primary mapping of target genes and the identification of QPtn.sau-4B laid a foundation for the following fine mapping and cloning work.