Comparative analysis of immune infiltrates in head and neck cancers across anatomical sites.

BACKGROUND: The response rate to immune checkpoint inhibitors targeting programmed cell death 1 (PD-1) receptor is 13%-18% for patients with recurrent or metastatic head and neck squamous cell carcinoma (HNSCC). Detailed understanding of the tumor immune microenvironment (TIME) is crucial in order to explain and improve this response rate. HNSCCs arise at various anatomical locations including the oral cavity, hypopharynx, larynx and oropharynx. Studies directly comparing immune infiltration between anatomical sites are scarce. Since the distinct locations could drive deviating microenvironments, we questioned whether the immune composition varies across these HNSCC sites. METHODS: Here, we characterized the TIME of 76 fresh tumor specimens using flow cytometry and performed single-cell RNA-sequencing on nine head and neck tumor samples. RESULTS: We found major differences in the composition of the TIME between patients. When comparing anatomical sites: tumors originating from the oral cavity had higher T cell infiltrates than tumors from other anatomical sites. The percentage of tumor-infiltrating T-lymphocytes positive for the immune checkpoint PD-1 varied considerably between patients, with the highest fraction of PD-1+ T cells found in larynx squamous cell carcinomas (SCCs). While we had hypothesized that the anatomical sites of tumor origin would drive sample clustering, our data showed that the type of TIME was more dominant and was particularly driven by the fraction of T cells positive for PD-1. Moreover, a high proportion of PD-1+ CD8+ T cells associated with an improved overall survival. Using single-cell RNA-sequencing, we observed that PD-1 expression was highest in the CD8-ENTPD1 tissue resident memory T cell/exhausted T cell and CD4-CXCL13 type 1 T helper cell clusters. CONCLUSIONS: We found that oral cavity SCCs had the highest frequencies of T cells. We also observed considerable interpatient heterogeneity for PD-1 on T cells, with noticeably higher frequencies of PD-1+ CD4+ T helper cells in larynx SCCs. Within the entire cohort, a higher fraction of CD8+ T cells positive for PD-1 was linked to improved overall survival. Whether the fraction of PD-1+ T cells within the TIME enables immune checkpoint inhibitor response prediction for patients with head and neck cancer remains to be determined.

T Cells Engineered to Express Immunoreceptors Targeting the Frequently Expressed Medullary Thyroid Cancer Antigens Calcitonin, CEA, and RET M918T.

Background: Medullary thyroid cancer (MTC) is a rare malignancy originating from the calcitonin-producing C cells of the thyroid. Despite recent therapeutic advances, metastatic MTC remains incurable. Adoptive cell therapy (ACT) using genetically engineered T cells targeting either tissue-restricted tumor-associated antigens or mutated neoantigens has led to durable remissions in other metastatic solid tumors. The majority of MTC express the tumor-associated antigens calcitonin and carcinoembryonic antigen (CEA), and ∼40% of MTC harbor the RET M918T oncogenic driver mutation. Methods: We developed and characterized three immunoreceptors that recognize extracellular CEA, a calcitonin epitope presented by HLA-A*24:02, or an RET M918T neoepitope restricted by HLA-DPB1*04:01/02. The chimeric antigen receptor (CAR) targeting CEA was synthetically designed, while the T cell receptors (TCRs) targeting calcitonin and RET M918T were isolated from a transgenic mouse and patient with MTC, respectively. These immunoreceptors were genetically engineered into peripheral blood T cells and tested for antigen specificity and antitumor activity. Results: T cells expressing the anti-CEA CAR or the calcitonin-reactive TCR produced effector cytokines and displayed cytotoxicity against cell lines expressing their cognate antigen in vitro. In immunodeficient mice harboring a human MTC cell line, the adoptive transfer of T cells engineered to express the anti-CEA CAR or calcitonin-reactive TCR led to complete tumor regression. T cells expressing the HLA-DPB1*04:01/02-restricted TCR targeting RET M918T, which was cloned from peripheral blood CD4+ T cells of a patient with MTC, demonstrated specific reactivity against cells pulsed with the mutated peptide and MTC tumor cells that expressed HLA-DPB1*04:01 and RET M918T. Conclusion: The preclinical data presented herein demonstrate the potential of using genetically engineered T cells targeting CEA, calcitonin, and/or RET M918T to treat metastatic MTC.

Transcriptomic profiles of neoantigen-reactive T cells in human gastrointestinal cancers.

Tumor-infiltrating neoantigen-reactive T cells can mediate regression of metastatic gastrointestinal cancers yet remain poorly characterized. We performed immunological screening against personalized neoantigens in combination with single-cell RNA sequencing on tumor-infiltrating lymphocytes from bile duct and pancreatic cancer patients to characterize the transcriptomic landscape of neoantigen-reactive T cells. We found that most neoantigen-reactive CD8+ T cells displayed an exhausted state with significant CXCL13 and GZMA co-expression compared with non-neoantigen-reactive bystander cells. Most neoantigen-reactive CD4+ T cells from a patient with bile duct cancer also exhibited an exhausted phenotype but with overexpression of HOPX or ADGRG1 while lacking IL7R expression. Thus, neoantigen-reactive T cells infiltrating gastrointestinal cancers harbor distinct transcriptomic signatures, which may provide new opportunities for harnessing these cells for therapy.

Disruption of global hypothalamic microRNA (miR) profiles and associated behavioral changes in California mice (Peromyscus californicus) developmentally exposed to endocrine disrupting chemicals.

Developmental exposure to endocrine disrupting chemicals (EDCs), e.g., bisphenol A (BPA) or genistein (GEN), causes longstanding epigenome effects. MicroRNAs (miRs) regulate which mRNAs will be translated to proteins and thereby serve as the final checkpoint in epigenetic control. Scant amount is known, however, whether EDCs affect neural miRNA (miR) patterns. We aimed to test the hypothesis that developmental exposure of California mice (Peromyscus californicus) to GEN, BPA, or both chemicals influences hypothalamic miR/small RNA profiles and ascertain the extent such biomolecular alterations correlate with behavioral and metabolic changes. California mice were developmentally exposed to GEN (250 mg/kg feed weight, FW), GEN (250 mg/kg FW)+BPA (5 mg/kg FW), low dose (LD) BPA (5 mg/kg FW), or upper dose (UD) BPA (50 mg/kg FW). Adult offspring were tested in a battery of behavioral and metabolic tests; whereupon, mice were euthanized, brains were collected and frozen, small RNAs were isolated from hypothalamic punches, and subsequently sequenced. California mice exposed to one or both EDCs engaged in one or more repetitive behaviors. GEN, LD BPA, and UD BPA altered aspects of ultrasonic and audible vocalizations. Each EDC exposure led to sex-dependent differences in differentially expressed miR/small RNAs with miR7-2, miR146, and miR148a being increased in all female and male EDC exposed groups. Current findings reveal that developmental exposure to GEN and/or BPA affects hypothalamic miR/small RNA expression patterns, and such changes correlate with EDC-induced behavioral and metabolic alterations. miR146 is likely an important mediator and biomarker of EDC exposure in mammals, including humans.

MYBPC3 Haplotype Linked to Hypertrophic Cardiomyopathy in Rhesus Macaques (Macaca mulatta).

In humans, abnormal thickening of the left ventricle of the heart clinically defines hypertrophic cardiomyopathy (HCM), a common inherited cardiovascular disorder that can precede a sudden cardiac death event. The wide range of clinical presentations in HCM obscures genetic variants that may influence an individual's susceptibility to sudden cardiac death. Although exon sequencing of major sarcomere genes can be used to detect high-impact causal mutations, this strategy is successful in only half of patient cases. The incidence of left ventricular hypertrophy (LVH) in a managed research colony of rhesus macaques provides an excellent comparative model in which to explore the genomic etiology of severe HCM and sudden cardiac death. Because no rhesus HCM-associated mutations have been reported, we used a next-generation genotyping assay that targets 7 sarcomeric rhesus genes within 63 genomic sites that are orthologous to human genomic regions known to harbor HCM disease variants. Amplicon sequencing was performed on 52 macaques with confirmed LVH and 42 unrelated, unaffected animals representing both the Indian and Chinese rhesus macaque subspecies. Bias-reduced logistic regression uncovered a risk haplotype in the rhesus MYBPC3 gene, which is frequently disrupted in both human and feline HCM; this haplotype implicates an intronic variant strongly associated with disease in either homozygous or carrier form. Our results highlight that leveraging evolutionary genomic data provides a unique, practical strategy for minimizing population bias in complex disease studies.

Viral genome integration of canine papillomavirus 16.

Papillomaviruses infect humans and animals, most often causing benign proliferations on skin or mucosal surfaces. Rarely, these infections persist and progress to cancer. In humans, this transformation most often occurs with high-risk papillomaviruses, where viral integration is a critical event in carcinogenesis. The first aim of this study was to sequence the viral genome of canine papillomavirus (CPV) 16 from a pigmented viral plaque that progressed to metastatic squamous cell carcinoma in a dog. The second aim was to characterize multiple viral genomic deletions and translocations as well as host integration sites. The full viral genome was identified using a combination of PCR and high throughput sequencing. CPV16 is most closely related to chipapillomaviruses CPV4, CPV9, and CPV12 and we propose CPV16 be classified as a chipapillomavirus. Assembly of the full viral genome enabled identification of deletion of portions of the E1 and E2/E4 genes and two viral translocations within the squamous cell carcinoma. Genome walking was performed which identified four sites of viral integration into the host genome. This is the first description of integration of a canine papillomavirus into the host genome, raising the possibility that CPV16 may be a potential canine high-risk papillomavirus type.

Refinement of the canine CD1 locus topology and investigation of antibody binding to recombinant canine CD1 isoforms.

CD1 molecules are antigen-presenting glycoproteins primarily found on dendritic cells (DCs) responsible for lipid antigen presentation to CD1-restricted T cells. Despite their pivotal role in immunity, little is known about CD1 protein expression in dogs, notably due to lack of isoform-specific antibodies. The canine (Canis familiaris) CD1 locus was previously found to contain three functional CD1A genes: canCD1A2, canCD1A6, and canCD1A8, where two variants of canCD1A8, canCD1A8.1 and canCD1A8.2, were assumed to be allelic variants. However, we hypothesized that these rather represented two separate genes. Sequencing of three overlapping bacterial artificial chromosomes (BACs) spanning the entire canine CD1 locus revealed canCD1A8.2 and canCD1A8.1 to be located in tandem between canCD1A7 and canCD1C, and canCD1A8.1 was consequently renamed canCD1A9. Green fluorescent protein (GFP)-fused canine CD1 transcripts were recombinantly expressed in 293T cells. All proteins showed a highly positive GFP expression except for canine CD1d and a splice variant of canine CD1a8 lacking exon 3. Probing with a panel of anti-CD1 monoclonal antibodies (mAbs) showed that Ca13.9H11 and Ca9.AG5 only recognized canine CD1a8 and CD1a9 isoforms, and Fe1.5F4 mAb solely recognized canine CD1a6. Anti-CD1b mAbs recognized the canine CD1b protein, but also bound CD1a2, CD1a8, and CD1a9. Interestingly, Ca9.AG5 showed allele specificity based on a single nucleotide polymorphism (SNP) located at position 321. Our findings have refined the structure of the canine CD1 locus and available antibody specificity against canine CD1 proteins. These are important fundamentals for future investigation of the role of canine CD1 in lipid immunity.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes.

We used the Affymetrix(®) Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.

RNA-seq analyses of changes in the Anopheles gambiae transcriptome associated with resistance to pyrethroids in Kenya: identification of candidate-resistance genes and candidate-resistance SNPs.

BACKGROUND: The extensive use of pyrethroids for control of malaria vectors, driven by their cost, efficacy and safety, has led to widespread resistance. To favor their sustainable use, the World Health Organization (WHO) formulated an insecticide resistance management plan, which includes the identification of the mechanisms of resistance and resistance surveillance. Recognized physiological mechanisms of resistance include target site mutations in the para voltage-gated sodium channel, metabolic detoxification and penetration resistance. Such understanding of resistance mechanisms has allowed the development of resistance monitoring tools, including genotyping of the kdr mutation L1014F/S in the para gene. METHODS: The sequence-based technique RNA-seq was applied to study changes in the transcriptome of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from the Western Province of Kenya. The resulting gene expression profiles were compared to data in the most recent literature to derive a list of candidate resistance genes. RNA-seq data were analyzed also to identify sequence polymorphisms linked to resistance. RESULTS: A total of five candidate-resistance genes (AGAP04177, AGAP004572, AGAP008840, AGAP007530 and AGAP013036) were identified with altered expression between resistant and susceptible mosquitoes from West and East Africa. A change from G to C at position 36043997 of chromosome 3R resulting in A101G of the sulfotransferase gene AGAP009551 was significantly associated with the resistance phenotype (odds ratio: 5.10). The kdr L1014S mutation was detected at similar frequencies in both phenotypically resistant and susceptible mosquitoes, suggesting it is no longer fully predictive of the resistant phenotype. CONCLUSIONS: Overall, these results support the conclusion that resistance to pyrethroids is a complex and evolving phenotype, dependent on multiple gene functions including, but not limited to, metabolic detoxification. Functional convergence among metabolic detoxification genes may exist, with the role of each gene being modulated by the life history and selection pressure on mosquito populations. As a consequence, biochemical assays that quantify overall enzyme activity may be a more suitable method for predicting metabolic resistance than gene-based assays.

Complete genome sequence of canine papillomavirus type 16.

Papillomaviruses are epitheliotropic, circular, double-stranded DNA viruses within the family Papillomaviridae that are associated with benign and malignant tumors in humans and animals. We report the complete genome sequence of canine papillomavirus type 16 identified within multiple pigmented cutaneous plaques and squamous cell carcinoma from an intact female Basenji dog.

Probing functional polymorphisms in the dengue vector, Aedes aegypti.

BACKGROUND: Dengue is the most prevalent arboviral disease world-wide and its primary vector is the mosquito Aedes aegypti. The current lack of commercially-available vaccines makes control of vector populations the only effective strategy to prevent dengue transmission. Aedes aegypti geographic populations exhibit great variability in insecticide resistance and susceptibility to dengue infection. The characterization of single nucleotide polymorphisms (SNPs) as molecular markers to study quantitatively this variation is needed greatly because this species has a low abundance of microsatellite markers and limited known restriction fragments length polymorphisms (RFLPs) and single-strand conformation polymorphism (SSCP) markers. RESULTS: We used RNA-seq to characterize SNPs in three Ae. aegypti strains, including the Liverpool (LVP) strain, from which the current genome annotation is derived. We identified 131,764 unique genome locations with at least one alternative nucleotide to what is reported in the reference annotation. These comprised changes in both open-reading frames (ORFs) and untranslated regions (UTRs) of transcripts. An in depth-look at sequence variation in immunity genes revealed that those associated with autophagy, MD2-like receptors and Peptidoglycan Recognition Proteins had more sequence variation in their 3'UTRs than mutations associated with non-synonymous changes. This supports the conclusion that these genes had maintained their functional specificity while being adapted to different regulatory domains. In contrast, a number of peroxidases, serpins and Clip-domain serine proteases exhibited conservation of putative UTR regulatory sequences while displaying diversification of the ORFs. Transcriptome evidence also was found for ~2500 novel transcriptional units (NTUs) not annotated in the reference genome. CONCLUSIONS: The transcriptome-wide assessment of within and inter-strain polymorphisms in Ae. aegypti adds considerably to the number of molecular markers available for genetic studies in this mosquito. Additionally, data supporting NTU discovery emphasizes the need for continuous amendments of the reference genome annotation.

Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species.

BACKGROUND: The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly. RESULTS: In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies. CONCLUSIONS: Many current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

Identifying human-rhesus macaque gene orthologs using heterospecific SNP probes.

We genotyped a Chinese and an Indian-origin rhesus macaque using the Affymetrix Genome-Wide Human SNP Array 6.0 and cataloged 85,473 uniquely mapping heterospecific SNPs. These SNPs were assigned to rhesus chromosomes according to their probe sequence alignments as displayed in the human and rhesus reference sequences. The conserved gene order (synteny) revealed by heterospecific SNP maps is in concordance with that of the published human and rhesus macaque genomes. Using these SNPs' original human rs numbers, we identified 12,328 genes annotated in humans that are associated with these SNPs, 3674 of which were found in at least one of the two rhesus macaques studied. Due to their density, the heterospecific SNPs allow fine-grained comparisons, including approximate boundaries of intra- and extra-chromosomal rearrangements involving gene orthologs, which can be used to distinguish rhesus macaque chromosomes from human chromosomes.

Multiple clonal MLL fusions in a patient receiving CHOP-based chemotherapy.

MLL rearrangements were analysed in the blood of a patient receiving chemotherapy for diffuse large B-cell lymphoma using inverse polymerase chain reaction targeting exon 12, parallel sequencing and a custom algorithm design. Of thirteen MLL rearrangements detected, five were capable of generating MLL fusion genes, including MLL-MLLT3, the most common fusion in acute myeloid leukaemia (AML). Other fusions, all previously clinically unobserved, included MLL-NKD1, a fusion to the negative regulator of Wnt/ß-catenin signaling, a pathway linked to leukaemic cell proliferation. The majority of the fusions exhibited clonal persistence from before treatment until 6 months post-chemotherapy, suggesting the fusions may confer a survival advantage to the mutant clone. MLL breakpoints were partly clustered at a specific location, indicating commonality in the process of their formation. Further, the same MLL breakpoint location exhibited a 50-100-fold increase in C to T transitions, consistent with attack by activation-induced cytidine deaminase (AICDA). As is also observed in AML and acute lymphoblastic leukaemia, in this single patient setting, MLL is capable of interacting with multiple fusion partners. This finding defines a discrete site of MLL susceptibility to fragmentation, linked to possible deregulation of AICDA function.

Genome-wide SNP discovery in walnut with an AGSNP pipeline updated for SNP discovery in allogamous organisms.

BACKGROUND: A genome-wide set of single nucleotide polymorphisms (SNPs) is a valuable resource in genetic research and breeding and is usually developed by re-sequencing a genome. If a genome sequence is not available, an alternative strategy must be used. We previously reported the development of a pipeline (AGSNP) for genome-wide SNP discovery in coding sequences and other single-copy DNA without a complete genome sequence in self-pollinating (autogamous) plants. Here we updated this pipeline for SNP discovery in outcrossing (allogamous) species and demonstrated its efficacy in SNP discovery in walnut (Juglans regia L.). RESULTS: The first step in the original implementation of the AGSNP pipeline was the construction of a reference sequence and the identification of single-copy sequences in it. To identify single-copy sequences, multiple genome equivalents of short SOLiD reads of another individual were mapped to shallow genome coverage of long Sanger or Roche 454 reads making up the reference sequence. The relative depth of SOLiD reads was used to filter out repeated sequences from single-copy sequences in the reference sequence. The second step was a search for SNPs between SOLiD reads and the reference sequence. Polymorphism within the mapped SOLiD reads would have precluded SNP discovery; hence both individuals had to be homozygous. The AGSNP pipeline was updated here for using SOLiD or other type of short reads of a heterozygous individual for these two principal steps. A total of 32.6X walnut genome equivalents of SOLiD reads of vegetatively propagated walnut scion cultivar 'Chandler' were mapped to 48,661 'Chandler' bacterial artificial chromosome (BAC) end sequences (BESs) produced by Sanger sequencing during the construction of a walnut physical map. A total of 22,799 putative SNPs were initially identified. A total of 6,000 Infinium II type SNPs evenly distributed along the walnut physical map were selected for the construction of an Infinium BeadChip, which was used to genotype a walnut mapping population having 'Chandler' as one of the parents. Genotyping results were used to adjust the filtering parameters of the updated AGSNP pipeline. With the adjusted filtering criteria, 69.6% of SNPs discovered with the updated pipeline were real and could be mapped on the walnut genetic map. A total of 13,439 SNPs were discovered by BES re-sequencing. BESs harboring SNPs were in 677 FPC contigs covering 98% of the physical map of the walnut genome. CONCLUSION: The updated AGSNP pipeline is a versatile SNP discovery tool for a high-throughput, genome-wide SNP discovery in both autogamous and allogamous species. With this pipeline, a large set of SNPs were identified in a single walnut cultivar.

Transcriptome profiling of citrus fruit response to huanglongbing disease.

Huanglongbing (HLB) or "citrus greening" is the most destructive citrus disease worldwide. In this work, we studied host responses of citrus to infection with Candidatus Liberibacter asiaticus (CaLas) using next-generation sequencing technologies. A deep mRNA profile was obtained from peel of healthy and HLB-affected fruit. It was followed by pathway and protein-protein network analysis and quantitative real time PCR analysis of highly regulated genes. We identified differentially regulated pathways and constructed networks that provide a deep insight into the metabolism of affected fruit. Data mining revealed that HLB enhanced transcription of genes involved in the light reactions of photosynthesis and in ATP synthesis. Activation of protein degradation and misfolding processes were observed at the transcriptomic level. Transcripts for heat shock proteins were down-regulated at all disease stages, resulting in further protein misfolding. HLB strongly affected pathways involved in source-sink communication, including sucrose and starch metabolism and hormone synthesis and signaling. Transcription of several genes involved in the synthesis and signal transduction of cytokinins and gibberellins was repressed while that of genes involved in ethylene pathways was induced. CaLas infection triggered a response via both the salicylic acid and jasmonic acid pathways and increased the transcript abundance of several members of the WRKY family of transcription factors. Findings focused on the fruit provide valuable insight to understanding the mechanisms of the HLB-induced fruit disorder and eventually developing methods based on small molecule applications to mitigate its devastating effects on fruit production.

Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing.

BACKGROUND: The availability of low cost sequencing has spurred its application to discovery and typing of variation, including variation induced by mutagenesis. Mutation discovery is challenging as it requires a substantial amount of sequencing and analysis to detect very rare changes and distinguish them from noise. Also challenging are the cases when the organism of interest has not been sequenced or is highly divergent from the reference. RESULTS: We describe the development of a simple method for reduced representation sequencing. Input DNA was digested with a single restriction enzyme and ligated to Y adapters modified to contain a sequence barcode and to provide a compatible overhang for ligation. We demonstrated the efficiency of this method at SNP discovery using rice and arabidopsis. To test its suitability for the discovery of very rare SNP, one control and three mutagenized rice individuals (1, 5 and 10 mM sodium azide) were used to prepare genomic libraries for Illumina sequencers by ligating barcoded adapters to NlaIII restriction sites. For genome-dependent discovery 15-30 million of 80 base reads per individual were aligned to the reference sequence achieving individual sequencing coverage from 7 to 15×. We identified high-confidence base changes by comparing sequences across individuals and identified instances consistent with mutations, i.e. changes that were found in a single treated individual and were solely GC to AT transitions. For genome-independent discovery 70-mers were extracted from the sequence of the control individual and single-copy sequence was identified by comparing the 70-mers across samples to evaluate copy number and variation. This de novo "genome" was used to align the reads and identify mutations as above. Covering approximately 1/5 of the 380 Mb genome of rice we detected mutation densities ranging from 0.6 to 4 per Mb of diploid DNA depending on the mutagenic treatment. CONCLUSIONS: The combination of a simple and cost-effective library construction method, with Illumina sequencing, and the use of a bioinformatic pipeline allows practical SNP discovery regardless of whether a genomic reference is available.

Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection.

BACKGROUND: The application of next generation sequencing technologies and bioinformatic scripts to identify high frequency SNPs distributed throughout the peach genome is described. Three peach genomes were sequenced using Roche 454 and Illumina/Solexa technologies to obtain long contigs for alignment to the draft 'Lovell' peach sequence as well as sufficient depth of coverage for 'in silico' SNP discovery. DESCRIPTION: The sequences were aligned to the 'Lovell' peach genome released April 01, 2010 by the International Peach Genome Initiative (IPGI). 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle' were sequenced to add SNPs segregating in two breeding populations, Pop DF ('Dr. Davis' × 'F8, 1-42') and Pop DG ('Dr. Davis' × 'Georgia Belle'). Roche 454 sequencing produced 980,000 total reads with 236 Mb sequence for 'Dr. Davis' and 735,000 total reads with 172 Mb sequence for 'F8, 1-42'. 84 bp × 84 bp paired end Illumina/Solexa sequences yielded 25.5, 21.4, 25.5 million sequences for 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle', respectively. BWA/SAMtools were used for alignment of raw reads and SNP detection, with custom PERL scripts for SNP filtering. Velvet's Columbus module was used for sequence assembly. Comparison of aligned and overlapping sequences from both Roche 454 and Illumina/Solexa resulted in the selection of 6654 high quality SNPs for 'Dr. Davis' vs. 'F8, 1-42' and 'Georgia Belle', distributed on eight major peach genome scaffolds as defined from the 'Lovell' assembly. CONCLUSION: The eight scaffolds contained about 215-225 Mb of peach genomic sequences with one SNP/~ 40,000 bases. All sequences from Roche 454 and Illumina/Solexa have been submitted to NCBI for public use in the Short Read Archive database. SNPs have been deposited in the NCBI SNP database.

Assemblathon 1: a competitive assessment of de novo short read assembly methods.

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/.

Exploring the Transcriptome Landscape of Pomegranate Fruit Peel for Natural Product Biosynthetic Gene and SSR Marker Discovery(F).

Pomegranate fruit peel is rich in bioactive plant natural products, such as hydrolyzable tannins and anthocyanins. Despite their documented roles in human nutrition and fruit quality, genes involved in natural product biosynthesis have not been cloned from pomegranate and very little sequence information is available on pomegranate in the public domain. Shotgun transcriptome sequencing of pomegranate fruit peel cDNA was performed using RNA-Seq on the Illumina Genome Analyzer platform. Over 100 million raw sequence reads were obtained and assembled into 9,839 transcriptome assemblies (TAs) (>200 bp). Candidate genes for hydrolyzable tannin, anthocyanin, flavonoid, terpenoid and fatty acid biosynthesis and/or regulation were identified. Three lipid transfer proteins were obtained that may contribute to the previously reported IgE reactivity of pomegranate fruit extracts. In addition, 115 SSR markers were identified from the pomegranate fruit peel transcriptome and primers were designed for 77 SSR markers. The pomegranate fruit peel transcriptome set provides a valuable platform for natural product biosynthetic gene and SSR marker discovery in pomegranate. This work also demonstrates that next-generation transcriptome sequencing is an economical and effective approach for investigating natural product biosynthesis, identifying genes controlling important agronomic traits, and discovering molecular markers in non-model specialty crop species.