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.

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.

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.

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.

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.

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.

Genome-Scale Analysis of Programmed DNA Elimination Sites in Tetrahymena thermophila.

Genetically programmed DNA rearrangements can regulate mRNA expression at an individual locus or, for some organisms, on a genome-wide scale. Ciliates rely on a remarkable process of whole-genome remodeling by DNA elimination to differentiate an expressed macronucleus (MAC) from a copy of the germline micronucleus (MIC) in each cycle of sexual reproduction. Here we describe results from the first high-throughput sequencing effort to investigate ciliate genome restructuring, comparing Sanger long-read sequences from a Tetrahymena thermophila MIC genome library to the MAC genome assembly. With almost 25% coverage of the unique-sequence MAC genome by MIC genome sequence reads, we created a resource for positional analysis of MIC-specific DNA removal that pinpoints MAC genome sites of DNA elimination at nucleotide resolution. The widespread distribution of internal eliminated sequences (IES) in promoter regions and introns suggests that MAC genome restructuring is essential not only for what it removes (for example, active transposons) but also for what it creates (for example, splicing-competent introns). Consistent with the heterogeneous boundaries and epigenetically modulated efficiency of individual IES deletions studied to date, we find that IES sites are dramatically under-represented in the ∼25% of the MAC genome encoding exons. As an exception to this general rule, we discovered a previously unknown class of small (<500 bp) IES with precise elimination boundaries that can contribute the 3' exon of an mRNA expressed during genome restructuring, providing a new mechanism for expanding mRNA complexity in a developmentally regulated manner.

A label-free differential quantitative mass spectrometry method for the characterization and identification of protein changes during citrus fruit development.

BACKGROUND: Citrus is one of the most important and widely grown commodity fruit crops. In this study a label-free LC-MS/MS based shot-gun proteomics approach was taken to explore three main stages of citrus fruit development. These approaches were used to identify and evaluate changes occurring in juice sac cells in various metabolic pathways affecting citrus fruit development and quality. RESULTS: Protein changes in citrus juice sac cells were identified and quantified using label-free shotgun methodologies. Two alternative methods, differential mass-spectrometry (dMS) and spectral counting (SC) were used to analyze protein changes occurring during earlier and late stages of fruit development. Both methods were compared in order to develop a proteomics workflow that could be used in a non-model plant lacking a sequenced genome. In order to resolve the bioinformatics limitations of EST databases from species that lack a full sequenced genome, we established iCitrus. iCitrus is a comprehensive sequence database created by merging three major sources of sequences (HarvEST:citrus, NCBI/citrus/unigenes, NCBI/citrus/proteins) and improving the annotation of existing unigenes. iCitrus provided a useful bioinformatics tool for the high-throughput identification of citrus proteins. We have identified approximately 1500 citrus proteins expressed in fruit juice sac cells and quantified the changes of their expression during fruit development. Our results showed that both dMS and SC provided significant information on protein changes, with dMS providing a higher accuracy. CONCLUSION: Our data supports the notion of the complementary use of dMS and SC for label-free comparative proteomics, broadening the identification spectrum and strengthening the identification of trends in protein expression changes during the particular processes being compared.

Tensile force-dependent neurite elicitation via anti-beta1 integrin antibody-coated magnetic beads.

Previous work using glass microneedles to apply calibrated, localized force to neurons showed that tensile force is a sufficient signal for neurite initiation and elongation. However, previous studies did not examine the kinetics or probability of neurite initiation as a function of force or the rate of force application. Here we report the use of a new technique-magnetic bead force application-to systematically investigate the role of force in these phenomena with better ease of use and control over force than glass microneedles. Force-induced neurite initiation from embryonic chick forebrain neurons appeared to be a first-order random process whose rate increased with increasing force, and required the presence of peripheral microtubules. In addition, the probability of initiation was more than twofold lower for neurons exposed to rapid initial force ramps (450 pN/s) than for neurons exposed to slower ramps (1.5 and 11 pN/s). We observed a low force threshold for elongation (15-100 pN), which was not previously detected in chick forebrain neurites elongated by glass microneedles. Finally, neurites subjected to constant force elongated at variable instantaneous rates, and switched abruptly between elongation and retraction, similar to spontaneous, growth-cone-mediated outgrowth and microtubule dynamic instability.