Establishment of an eHAP1 human haploid cell line hybrid reference genome assembled from short and long reads

Abstract

Background: Haploid cell lines are a valuable research tool with broad applicability for genetic assays. As such the fully haploid human cell line, eHAP1, has been used in a wide array of studies. However, the absence of a corresponding reference genome sequence for this cell line has limited the potential for more widespread applications to experiments dependent on available sequence, like capture-clone methodologies.Results: We generated ~15x coverage Nanopore long reads from ten GridION flowcells. We utilized this data to assemble a de novo draft genome using minimap and miniasm and subsequently polished using Racon. This assembly was further polished using previously generated, low-coverage, Illumina short reads with Pilon and ntEdit. This resulted in a hybrid eHAP1 assembly with >90% complete BUSCO scores. We further assessed the eHAP1 long read data for structural variants using Sniffles and identify a variety of rearrangements, including a previously established Philadelphia translocation. Finally, we demonstrate how some of these variants overlap open chromatin regions, potentially impacting regulatory regions.Conclusions: By integrating both long and short reads, we generated a high-quality reference assembly for eHAP1 cells. We identify structural variants using long reads, including some that may impact putative regulatory elements. The union of long and short reads demonstrates the utility in combining sequencing platforms to generate a high-quality reference genome de novo solely from low coverage data. We expect the resulting eHAP1 genome assembly to provide a useful resource to enable novel experimental applications in this important model cell line. Introduction:The vast majority of eukaryotic cells are diploid and many cellular models used experimentally are either diploid or polyploid. The presence of additional alleles, while evolutionarily beneficial,

Identification of a Stable, Non-Canonically Regulated Nrf2 Form in Lung Cancer Cells

“…To assess the expression of specific NFE2L2 transcripts in A549 cell line, we have used the RNA sequencing data from the project available at NCBI Gene Expression Omnibus (GEO) public database (accession numbers GSM2308412, where A549 cell line was sequenced with paired Illumina protocol. Primary analysis of RNA-seq data included the quality control of sequenced reads with the use of FastQC (Andrews, 2010), reads trimming with the usage of Trimmomatic [ 16 ] and mapping to the reference genome based on NCBI reference human genome (assembly GRCh38.p13) and annotation (release 109) [ 17 ] with the Hisat2 aligner [ 18 ]. Further data preparations was performed with SAMtools software [ 19 ] and R software [ 20 ] together with Bioconductor platform.…”

G-SAIP: Graphical Sequence Alignment Through Parallel Programming in the Post-Genomic Era

“…The second dataset used was a WGS Illumina paired-end sequence reads of Drosophila melanogaster with 4 Gb of size per file 65 and accessible in SRA repository with SRR10735526 accession number. Homo sapiens genome 62 3.1 Gb GCA_000001405.28…”