BAC-Pool Sequencing and Analysis of Large Segments of A12 and D12 Homoeologous Chromosomes in Upland Cotton

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.     

Comprehensive Draft Genome Analyses of Three Rockfishes (Scorpaeniformes,Sebastiscus) via Genome Survey Sequencing

Sebastiscus species, marine rockfishes, are of essential economic value. However, the genomic data of this genus is lacking and incomplete. Here, whole genome sequencing of all species of Sebastiscus was conducted to provide fundamental genomic information. The genome sizes were estimated to be 802.49 Mb (S. albofasciatus), 786.79 Mb (S. tertius), and 776.00 Mb (S. marmoratus) by using k-mer analyses. The draft genome sequences were initially assembled, and genome-wide microsatellite motifs were identified. The heterozygosity, repeat ratios, and numbers of microsatellite motifs all suggested possibly that S. tertius is more closely related to S. albofasciatus than S. marmoratus at the genetic level. Moreover, the complete mitochondrial genome sequences were assembled from the whole genome data and the phylogenetic analyses genetically supported the validation of Sebastiscus species. This study provides an important genome resource for further studies of Sebastiscus species.     

Long Span DNA Paired-End-Tag (DNA-PET) Sequencing Strategy for the Interrogation of Genomic Structural Mutations and Fusion-Point-Guided Reconstruction of Amplicons

Structural variations (SVs) contribute significantly to the variability of the human genome and extensive genomic rearrangements are a hallmark of cancer. While genomic DNA paired-end-tag (DNA-PET) sequencing is an attractive approach to identify genomic SVs, the current application of PET sequencing with short insert size DNA can be insufficient for the comprehensive mapping of SVs in low complexity and repeat-rich genomic regions. We employed a recently developed procedure to generate PET sequencing data using large DNA inserts of 10–20 kb and compared their characteristics with short insert (1 kb) libraries for their ability to identify SVs. Our results suggest that although short insert libraries bear an advantage in identifying small deletions, they do not provide significantly better breakpoint resolution. In contrast, large inserts are superior to short inserts in providing higher physical genome coverage for the same sequencing cost and achieve greater sensitivity, in practice, for the identification of several classes of SVs, such as copy number neutral and complex events. Furthermore, our results confirm that large insert libraries allow for the identification of SVs within repetitive sequences, which cannot be spanned by short inserts. This provides a key advantage in studying rearrangements in cancer, and we show how it can be used in a fusion-point-guided-concatenation algorithm to study focally amplified regions in cancer.     

Molecular and genomic characterisation of a panel of human anal cancer cell lines

Anal cancer is a rare disease that has doubled in incidence over the last four decades. Current treatment and survival of patients with this disease has not changed substantially over this period of time, due, in part, to a paucity of preclinical models to assess new therapeutic options. To address this hiatus, we set-out to establish, validate and characterise a panel of human anal squamous cell carcinoma (ASCC) cell lines by employing an explant technique using fresh human ASCC tumour tissue. The panel of five human ASCC cell lines were validated to confirm their origin, squamous features and tumourigenicity, followed by molecular and genomic (whole-exome sequencing) characterisation. This panel recapitulates the genetic and molecular characteristics previously described in ASCC including phosphoinositide-3-kinase (PI3K) mutations in three of the human papillomavirus (HPV) positive lines and TP53 mutations in the HPV negative line. The cell lines demonstrate the ability to form tumouroids and retain their tumourigenic potential upon xenotransplantation, with varied inducible expression of major histocompatibility complex class I (MHC class I) and Programmed cell death ligand 1 (PD-L1). We observed differential responses to standard chemotherapy, radiotherapy and a PI3K specific molecular targeted agent in vitro, which correlated with the clinical response of the patient tumours from which they were derived. We anticipate this novel panel of human ASCC cell lines will form a valuable resource for future studies into the biology and therapeutics of this rare disease.Subject terms: Targeted therapies, Anal cancer, Immune evasion, Cancer models     

Analytical and Clinical Validation of a Digital Sequencing Panel for Quantitative,Highly Accurate Evaluation of Cell-Free Circulating Tumor DNA

Next-generation sequencing of cell-free circulating solid tumor DNA addresses two challenges in contemporary cancer care. First this method of massively parallel and deep sequencing enables assessment of a comprehensive panel of genomic targets from a single sample, and second, it obviates the need for repeat invasive tissue biopsies. Digital Sequencing^TM is a novel method for high-quality sequencing of circulating tumor DNA simultaneously across a comprehensive panel of over 50 cancer-related genes with a simple blood test. Here we report the analytic and clinical validation of the gene panel. Analytic sensitivity down to 0.1% mutant allele fraction is demonstrated via serial dilution studies of known samples. Near-perfect analytic specificity (> 99.9999%) enables complete coverage of many genes without the false positives typically seen with traditional sequencing assays at mutant allele frequencies or fractions below 5%. We compared digital sequencing of plasma-derived cell-free DNA to tissue-based sequencing on 165 consecutive matched samples from five outside centers in patients with stage III-IV solid tumor cancers. Clinical sensitivity of plasma-derived NGS was 85.0%, comparable to 80.7% sensitivity for tissue. The assay success rate on 1,000 consecutive samples in clinical practice was 99.8%. Digital sequencing of plasma-derived DNA is indicated in advanced cancer patients to prevent repeated invasive biopsies when the initial biopsy is inadequate, unobtainable for genomic testing, or uninformative, or when the patient’s cancer has progressed despite treatment. Its clinical utility is derived from reduction in the costs, complications and delays associated with invasive tissue biopsies for genomic testing.     

Genomic rearrangements at the FRA2H common fragile site frequently involve non-homologous recombination events across LTR and L1(LINE) repeats

Common fragile sites (cFSs) are non-random chromosomal regions that are prone to breakage under conditions of replication stress. DNA damage and chromosomal alterations at cFSs appear to be critical events in the development of various human diseases, especially carcinogenesis. Despite the growing interest in understanding the nature of cFS instability, only a few cFSs have been molecularly characterised. In this study, we fine-mapped the location of FRA2H using six-colour fluorescence in situ hybridisation and showed that it is one of the most active cFSs in the human genome. FRA2H encompasses approximately 530 kb of a gene-poor region containing a novel large intergenic non-coding RNA gene (AC097500.2). Using custom-designed array comparative genomic hybridisation, we detected gross and submicroscopic chromosomal rearrangements involving FRA2H in a panel of 54 neuroblastoma, colon and breast cancer cell lines. The genomic alterations frequently involved different classes of long terminal repeats and long interspersed nuclear elements. An analysis of breakpoint junction sequence motifs predominantly revealed signatures of microhomology-mediated non-homologous recombination events. Our data provide insight into the molecular structure of cFSs and sequence motifs affected by their activation in cancer. Identifying cFS sequences will accelerate the search for DNA biomarkers and targets for individualised therapies.     

Methylation-sensitive high resolution melting (MS-HRM): a new approach for sensitive and high-throughput assessment of methylation

In this article, we show that high resolution melting analysis (HRM) is a sensitive and specific method for the detection of methylation. Methylated DNA and unmethylated DNA acquire different sequences after bisulphite treatment resulting in PCR products with markedly different melting profiles. We used PCR to amplify both methylated and unmethylated sequences and assessed HRM for the determination of the methylation status of the MGMT promoter region. Reconstruction experiments showed that MGMT methylation could be detected at levels as low as 0.1%. Moreover, MS-HRM allows for estimation of the methylation level by comparing the melting profiles of unknown PCR products to the melting profiles of PCR products derived from standards with a known unmethylated to methylated template ratio. We used MS-HRM for the analysis of eight cell lines of known methylation status and a panel of colorectal cancer specimens. The simplicity and high reproducibility of the MS-HRM protocol makes MS-HRM the method of choice for methylation assessment in many diagnostic and research applications.     

Draft Sequences of the Radish (Raphanus sativus L.) Genome

Radish (Raphanus sativus L., n = 9) is one of the major vegetables in Asia. Since the genomes of Brassica and related species including radish underwent genome rearrangement, it is quite difficult to perform functional analysis based on the reported genomic sequence of Brassica rapa. Therefore, we performed genome sequencing of radish. Short reads of genomic sequences of 191.1 Gb were obtained by next-generation sequencing (NGS) for a radish inbred line, and 76,592 scaffolds of ≥300 bp were constructed along with the bacterial artificial chromosome-end sequences. Finally, the whole draft genomic sequence of 402 Mb spanning 75.9% of the estimated genomic size and containing 61,572 predicted genes was obtained. Subsequently, 221 single nucleotide polymorphism markers and 768 PCR-RFLP markers were used together with the 746 markers produced in our previous study for the construction of a linkage map. The map was combined further with another radish linkage map constructed mainly with expressed sequence tag-simple sequence repeat markers into a high-density integrated map of 1,166 cM with 2,553 DNA markers. A total of 1,345 scaffolds were assigned to the linkage map, spanning 116.0 Mb. Bulked PCR products amplified by 2,880 primer pairs were sequenced by NGS, and SNPs in eight inbred lines were identified.     

Comprehensive whole genome survey analyses of male and female brown-spotted flathead fish Platycephalus sp.1

The flathead fish Platycephalus sp.1 is an ecologically and commercially important marine fish in the northwestern Pacific with notable sexual differences in growth and development. Yet the genomic data of this species is lacking. In the present study, whole genome sequencing of two individuals (one male and one female) of Platycephalus sp.1 were conducted to provide fundamental genomic information. The genome sizes were estimated to be 674.96 Mb (male) and 684.15 Mb (female) by using k-mer analyses. The heterozygosity and repeat ratios suggested possible male heterogamety of this species. The draft genome sequences were initially assembled and genome-wide microsatellite motifs were identified. Besides, the complete mitochondrial genome sequences were assembled and the phylogenetic analyses genetically supported the validation of Platycephalus sp.1. The reported genomic data and genetic markers in this study could be useful in future comparative genomics and evolutionary biology studies.