A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

Correction: Copy number of pancreatic polypeptide receptor gene NPY4R correlates with body mass index and waist circumference.

[This corrects the article DOI: 10.1371/journal.pone.0194668.].

Identification and rescue of a tRNA wobble inosine deficiency causing intellectual disability disorder.

The deamination of adenosine to inosine at the wobble position of tRNA is an essential post-transcriptional RNA modification required for wobble decoding in bacteria and eukaryotes. In humans, the wobble inosine modification is catalyzed by the heterodimeric ADAT2/3 complex. Here, we describe novel pathogenic ADAT3 variants impairing adenosine deaminase activity through a distinct mechanism that can be corrected through expression of the heterodimeric ADAT2 subunit. The variants were identified in a family in which all three siblings exhibit intellectual disability linked to biallelic variants in the ADAT3 locus. The biallelic ADAT3 variants result in a missense variant converting alanine to valine at a conserved residue or the introduction of a premature stop codon in the deaminase domain. Fibroblast cells derived from two ID-affected individuals exhibit a reduction in tRNA wobble inosine levels and severely diminished adenosine tRNA deaminase activity. Notably, the ADAT3 variants exhibit impaired interaction with the ADAT2 subunit and alterations in ADAT2-dependent nuclear localization. Based upon these findings, we find that tRNA adenosine deaminase activity and wobble inosine modification can be rescued in patient cells by overexpression of the ADAT2 catalytic subunit. These results uncover a key role for the inactive ADAT3 deaminase domain in proper assembly with ADAT2 and demonstrate that ADAT2/3 nuclear import is required for maintaining proper levels of the wobble inosine modification in tRNA.

Linkage and exome analysis implicate multiple genes in non-syndromic intellectual disability in a large Swedish family.

BACKGROUND: Non-syndromic intellectual disability is genetically heterogeneous with dominant, recessive and complex forms of inheritance. We have performed detailed genetic studies in a large multi-generational Swedish family, including several members diagnosed with non-syndromic intellectual disability. Linkage analysis was performed on 22 family members, nine affected with mild to moderate intellectual disability and 13 unaffected family members. METHODS: Family members were analyzed with Affymetrix Genome-Wide Human SNP Array 6.0 and the genetic data was used to detect copy number variation and to perform genome wide linkage analysis with the SNP High Throughput Linkage analysis system and the Merlin software. For the exome sequencing, the samples were prepared using the Sure Select Human All Exon Kit (Agilent Technologies, Santa Clara, CA, USA) and sequenced using the Ion Proton™ System. Validation of identified variants was performed with Sanger sequencing. RESULTS: The linkage analysis results indicate that intellectual disability in this family is genetically heterogeneous, with suggestive linkage found on chromosomes 1q31-q41, 4q32-q35, 6p25 and 14q24-q31 (LOD scores of 2.4, simulated p-value of 0.000003 and a simulated genome-wide p-value of 0.06). Exome sequencing was then performed in 14 family members and 7 unrelated individuals from the same region. The analysis of coding variation revealed a pathogenic and candidate variants in different branches of the family. In three patients we find a known homozygous pathogenic mutation in the Homo sapiens solute carrier family 17 member 5 (SLC17A5), causing Salla disease. We also identify a deletion overlapping KDM3B and a duplication overlapping MAP3K4 and AGPAT4, both overlapping variants previously reported in developmental disorders. CONCLUSIONS: DNA samples from the large family analyzed in this study were initially collected based on a hypothesis that affected members shared a major genetic risk factor. Our results show that a complex phenotype such as mild intellectual disability in large families from genetically isolated populations may show considerable genetic heterogeneity.

Copy number determination of the gene for the human pancreatic polypeptide receptor NPY4R using read depth analysis and droplet digital PCR.

BACKGROUND: Copy number variation (CNV) plays an important role in human genetic diversity and has been associated with multiple complex disorders. Here we investigate a CNV on chromosome 10q11.22 that spans NPY4R, the gene for the appetite-regulating pancreatic polypeptide receptor Y4. This genomic region has been challenging to map due to multiple repeated elements and its precise organization has not yet been resolved. Previous studies using microarrays were interpreted to show that the most common copy number was 2 per genome. RESULTS: We have investigated 18 individuals from the 1000 Genomes project using the well-established method of read depth analysis and the new droplet digital PCR (ddPCR) method. We find that the most common copy number for NPY4R is 4. The estimated number of copies ranged from three to seven based on read depth analyses with Control-FREEC and CNVnator, and from four to seven based on ddPCR. We suggest that the difference between our results and those published previously can be explained by methodological differences such as reference gene choice, data normalization and method reliability. Three high-quality archaic human genomes (two Neanderthal and one Denisova) display four copies of the NPY4R gene indicating that a duplication occurred prior to the human-Neanderthal/Denisova split. CONCLUSIONS: We conclude that ddPCR is a sensitive and reliable method for CNV determination, that it can be used for read depth calibration in CNV studies based on already available whole-genome sequencing data, and that further investigation of NPY4R copy number variation and its consequences are necessary due to the role of Y4 receptor in food intake regulation.

Expression profiling and in situ screening of circular RNAs in human tissues.

Circular RNAs (circRNAs) were recently discovered as a class of widely expressed noncoding RNA and have been implicated in regulation of gene expression. However, the function of the majority of circRNAs remains unknown. Studies of circRNAs have been hampered by a lack of essential approaches for detection, quantification and visualization. We therefore developed a target-enrichment sequencing method suitable for screening of circRNAs and their linear counterparts in large number of samples. We also applied padlock probes and in situ sequencing to visualize and determine circRNA localization in human brain tissue at subcellular levels. We measured circRNA abundance across different human samples and tissues. Our results highlight the potential of this RNA class to act as a specific diagnostic marker in blood and serum, by detection of circRNAs from genes exclusively expressed in the brain. The powerful and scalable tools we present will enable studies of circRNA function and facilitate screening of circRNA as diagnostic biomarkers.

Copy number of pancreatic polypeptide receptor gene NPY4R correlates with body mass index and waist circumference.

Multiple genetic studies have linked copy number variation (CNV) in different genes to body mass index (BMI) and obesity. A CNV on chromosome 10q11.22 has been associated with body weight. This CNV region spans NPY4R, the gene encoding the pancreatic polypeptide receptor Y4, which has been described as a satiety-stimulating receptor. We have investigated CNV of the NPY4R gene and analysed its relationship to BMI, waist circumference and self-reported dietary intake from 558 individuals (216 men and 342 women) representing a wide BMI range. The copy number for NPY4R ranged from 2 to 8 copies (average 4.6±0.8). Rather than the expected negative correlation, we observed a positive correlation between NPY4R copy number and BMI as well as waist circumference in women (Pearson's r = 0.267, p = 2.65×10−7 and r = 0.256, p = 8×10−7, respectively). Each additional copy of NPY4R correlated with 2.6 kg/m2 increase in BMI and 5.67 cm increase in waist circumference (p = 2.8×10−5 and p = 6.2×10−5, respectively) for women. For men, there was no statistically significant correlation between CNV and BMI. Our results suggest that NPY4R genetic variation influences body weight in women, but the exact role of this receptor appears to be more complex than previously proposed.

Genome-wide association study of susceptibility loci for cervical cancer.

BACKGROUND: Cervical carcinoma has a heritable genetic component, but the genetic basis of cervical cancer is still not well understood. METHODS: We performed a genome-wide association study of 731 422 single nucleotide polymorphisms (SNPs) in 1075 cervical cancer case subjects and 4014 control subjects and replicated it in 1140 case subjects and 1058 control subjects. The association between top SNPs and cervical cancer was estimated by odds ratios (ORs) and 95% confidence intervals (CIs) with unconditional logistic regression. All statistical tests were two-sided. RESULTS: Three independent loci in the major histocompatibility complex (MHC) region at 6p21.3 were associated with cervical cancer: the first is adjacent to the MHC class I polypeptide-related sequence A gene (MICA) (rs2516448; OR = 1.42, 95% CI = 1.31 to 1.54; P = 1.6×10(-18)); the second is between HLA-DRB1 and HLA-DQA1 (rs9272143; OR = 0.67, 95% CI = 0.62 to 0.72; P = 9.3×10(-24)); and the third is at HLA-DPB2 (rs3117027; OR=1.25, 95% CI = 1.15 to 1.35; P = 4.9×10(-8)). We also confirmed previously reported associations of B*0702 and DRB1*1501-DQB1*0602 with susceptibility to and DRB1*1301-DQA1*0103-DQB1*0603 with protection against cervical cancer. The three new loci are statistically independent of these specific human leukocyte antigen alleles/haplotypes. MICA encodes a membrane-bound protein that acts as a ligand for NKG2D to activate antitumor effects. The risk allele of rs2516448 is in perfect linkage disequilibrium with a frameshift mutation (A5.1) of MICA, which results in a truncated protein. Functional analysis shows that women carrying this mutation have lower levels of membrane-bound MICA. CONCLUSIONS: Three novel loci in the MHC may affect susceptibility to cervical cancer in situ, including the MICA-A5.1 allele that may cause impaired immune activation and increased risk of tumor development.

Infantile cerebellar-retinal degeneration associated with a mutation in mitochondrial aconitase, ACO2.

Degeneration of the cerebrum, cerebellum, and retina in infancy is part of the clinical spectrum of lysosomal storage disorders, mitochondrial respiratory chain defects, carbohydrate glycosylation defects, and infantile neuroaxonal dystrophy. We studied eight individuals from two unrelated families who presented at 2-6 months of age with truncal hypotonia and athetosis, seizure disorder, and ophthalmologic abnormalities. Their course was characterized by failure to acquire developmental milestones and culminated in profound psychomotor retardation and progressive visual loss, including optic nerve and retinal atrophy. Despite their debilitating state, the disease was compatible with survival of up to 18 years. Laboratory investigations were normal, but the oxidation of glutamate by muscle mitochondria was slightly reduced. Serial brain MRI displayed progressive, prominent cerebellar atrophy accompanied by thinning of the corpus callosum, dysmyelination, and frontal and temporal cortical atrophy. Homozygosity mapping followed by whole-exome sequencing disclosed a Ser112Arg mutation in ACO2, encoding mitochondrial aconitase, a component of the Krebs cycle. Specific aconitase activity in the individuals' lymphoblasts was severely reduced. Under restrictive conditions, the mutant human ACO2 failed to complement a yeast ACO1 deletion strain, whereas the wild-type human ACO2 succeeded, indicating that this mutation is pathogenic. Thus, a defect in mitochondrial aconitase is associated with an infantile neurodegenerative disorder affecting mainly the cerebellum and retina. In the absence of noninvasive biomarkers, determination of the ACO2 sequence or of aconitase activity in lymphoblasts are warranted in similarly affected individuals, based on clinical and neuroradiologic grounds.

Characterizing and interpreting genetic variation from personal genome sequencing.

Since the completion of the human genome project, there has been enormous progress in the development of novel technologies for DNA sequencing. The advent of next-generation sequencing technologies now makes it possible to sequence an entire human genome in one or a few experiments. As a consequence, several individual human genomes have now been fully sequenced, using different experimental strategies. Although the protocols differ between the various sequencing technologies, the challenges of analyzing the data, calling variation, and interpreting the results are similar for all platforms. Here, we give an overview of the human genome sequencing projects completed to date. The strategies for aligning sequence reads and extracting information about different types of genetic variation from the sequence data are discussed. Identification of structural variation, such as copy number variation and insertion-deletion variants, can be complex, and there are a plethora of algorithms and analysis tools available. We also give an overview of the challenge of interpreting the whole-genome sequence data both from a technical and clinical perspective.

Multiple recurrent genetic events converge on control of histone lysine methylation in medulloblastoma.

We used high-resolution SNP genotyping to identify regions of genomic gain and loss in the genomes of 212 medulloblastomas, malignant pediatric brain tumors. We found focal amplifications of 15 known oncogenes and focal deletions of 20 known tumor suppressor genes (TSG), most not previously implicated in medulloblastoma. Notably, we identified previously unknown amplifications and homozygous deletions, including recurrent, mutually exclusive, highly focal genetic events in genes targeting histone lysine methylation, particularly that of histone 3, lysine 9 (H3K9). Post-translational modification of histone proteins is critical for regulation of gene expression, can participate in determination of stem cell fates and has been implicated in carcinogenesis. Consistent with our genetic data, restoration of expression of genes controlling H3K9 methylation greatly diminishes proliferation of medulloblastoma in vitro. Copy number aberrations of genes with critical roles in writing, reading, removing and blocking the state of histone lysine methylation, particularly at H3K9, suggest that defective control of the histone code contributes to the pathogenesis of medulloblastoma.

Excessive genomic DNA copy number variation in the Li-Fraumeni cancer predisposition syndrome.

DNA copy number variations (CNVs) are a significant and ubiquitous source of inherited human genetic variation. However, the importance of CNVs to cancer susceptibility and tumor progression has not yet been explored. Li-Fraumeni syndrome (LFS) is an autosomal dominantly inherited disorder characterized by a strikingly increased risk of early-onset breast cancer, sarcomas, brain tumors and other neoplasms in individuals harboring germline TP53 mutations. Known genetic determinants of LFS do not fully explain the variable clinical phenotype in affected family members. As part of a wider study of CNVs and cancer, we conducted a genome-wide profile of germline CNVs in LFS families. Here, by examining DNA from a large healthy population and an LFS cohort using high-density oligonucleotide arrays, we show that the number of CNVs per genome is well conserved in the healthy population, but strikingly enriched in these cancer-prone individuals. We found a highly significant increase in CNVs among carriers of germline TP53 mutations with a familial cancer history. Furthermore, we identified a remarkable number of genomic regions in which known cancer-related genes coincide with CNVs, in both LFS families and healthy individuals. Germline CNVs may provide a foundation that enables the more dramatic chromosomal changes characteristic of TP53-related tumors to be established. Our results suggest that screening families predisposed to cancer for CNVs may identify individuals with an abnormally high number of these events.

Further evidence for role of a promoter variant in the TNFRSF6 gene in Alzheimer disease.

The tumor necrosis factor receptor superfamily, member 6 gene (TNFRSF6) is situated on chromosome 10q, near the region indicated in several AD linkage studies. In a previous study a significant association was found between a promoter variant within the TNFRSF6 gene and Alzheimer disease (AD). To further investigate the TNFRSF6 region, 34 SNPs within 180 kb were first genotyped in an exploratory set of 121 early-onset dementia cases and 152 controls in order to establish the extent of linkage disequilibrium (LD) and the major haplotypes in the region. This analysis showed that the LD was clustered in two large blocks within each of which a limited number of haplotypes represented most of the diversity. Haplotype tagging markers were chosen and genotyped in an additional 204 late onset AD cases and 177 controls. Tests of association were performed for single markers for case/control status and for a quantitative measure of cognitive ability [Mini-Mental State Examination (MMSE) scores] within cases only. The previously associated marker, located in the promoter of TNFRSF6, now gave significant association with cognitive status in the Scottish early-onset dementia samples (p = 0.005) with the strongest signals being evident in the APOE-e4 carrier subgroup, thus providing a second independent positive result for this marker. The same marker was also significantly associated with cognitive performance in the Swedish APOE-e4 carriers (p = 0.014), providing a third independent signal. Association analysis was also performed using the major haplotypes in the region, employing both case/control status and MMSE scores. The haplotype analysis did not give further significant findings. These results together with previous data suggest that a promoter marker in TNFRSF6 plays a moderate but demonstrable role in AD etiology.