Comprehensive Sequence Analysis Research Articles

Genome-wide characterization of the SPL gene family involved in the age development of Jatropha curcas

BackgroundSPL (SQUAMOSA-promoter binding protein-like) proteins form a large family of plant-specific transcription factors that play essential roles in various aspects of plant growth and development. They are potentially important candidates for genetic improvement of agronomic traits. However, there were limited information about the SPL genes in Jatropha curcas, an important biofuel plant.ResultsIn Jatropha, 15 JcSPL genes were identified. Phylogenetic analysis revealed that most of the JcSPLs were closely related to SPLs from woody plant rather than herbaceous plant and distantly related to monocotyledon SPLs. Gene structure, conserved motif and repetitive sequence analysis indicated diverse and specific functions of some JcSPL genes. By combination of target prediction and degradome sequencing analysis, 10 of the 15 JcSPLs were shown to be targets of JcmiR156. Quantitative PCR analysis showed diversified spatial-temporal expression patterns of JcSPLs. It is interesting that the expression levels of JcSPL3 were the highest in all tissues examined in 7- or 10-year-old plants and exhibited increasing trend with plant age, suggesting its important role in the regulation of age development in Jatropha. Overexpression of JcSPL3 in Arabidopsis resulted in earlier flowering time, shorter silique length and reduced biomass of roots.ConclusionsThrough comprehensive and systematic analysis of phylogenetic relationships, conserved motifs, gene structures, chromosomal locations, repetitive sequence and expression patterns, 15 JcSPL genes were identified in Jatropha and characterized in great detail. These results provide deep insight into the evolutionary origin and biological significance of plant SPLs and lay the foundation for further functional characterization of JcSPLs with the purpose of genetic improvement in Jatropha.

A Landscape of Murine Long Non-Coding RNAs Reveals the Leading Transcriptome Alterations in Adipose Tissue during Aging

SUMMARYAging is an inevitable process that involves profound physiological changes. Long non-coding RNAs (lncRNAs) are emerging as important regulators in various biological processes but are not systemically studied in aging. To provide an organism-wide lncRNA landscape during aging, we conduct comprehensive RNA sequencing (RNA-seq) analyses across the mouse lifespan. Of the 1,675 aging-regulated lncRNAs (AR-lncRNAs) identified, the majority are connected to inflammation-related biological pathways. AR-lncRNAs exhibit high tissue specificity; conversely, those with higher tissue specificity are preferentially regulated during aging. White adipose tissue (WAT) displays the highest number of AR-lncRNAs and develops the most dynamic crosstalk between AR-lncRNA and AR-mRNA during aging. An adipose-enriched AR-lncRNA, lnc-adipoAR1, is negatively correlated with aging, and knocking it down inhibits adipogenesis, phenocopying the compromised adipogenic capacity of aged fat. Our works together reveal AR-lncRNAs as essential components in aging and suggest that although each tissue ages in a distinct manner, WAT is a leading contributor to aging-related health decline.

Identification and Quantitation of Amylase Trypsin Inhibitors Across Cultivars Representing the Diversity of Bread Wheat.

α-Amylase/trypsin inhibitors (ATIs) may have a role in nonceliac wheat sensitivity (NCWS) and celiac disease (CD), but the ATI content and diversity across a range of wheat cultivars are not well characterized. Discovery proteomics was used to detect ATIs across two wheat cultivars: Chara and Magenta. Comprehensive mapping of detected ATIs with the ATIs from the recently published wheat genome RefSeq v1.0 shows the presence of three major subclasses: monomeric (9%), dimeric (61%), and chloroform-methanol (CM) type (30%). Subsequently, the level of 18 ATI isoforms (63 peptides) grouped into four subtypes was monitored across 15 commercial wheat cultivars and the eight parental lines from a multiparent advanced-generation intercross (MAGIC) population using liquid chromatography-multiple reaction monitoring-mass spectrometry (LC-MRM-MS). The ATI content of wheat cultivars Janz, Sunvale, Diamond Bird, and Longreach Scout was significantly lower than that of other wheat cultivars. The MAGIC parental cultivars Baxter and Xiaoyan 54 contain higher levels (∼115% relative to the average wheat ATI content), whereas cultivar Pastor contained the lowest levels (∼87%). Comprehensive sequence analysis, annotation, chromosomal locations, and epitope mapping enabled us to build an LC-MRM-MS method to monitor and quantify the immunostimulatory ATI proteins potentially related to NCWS, autoimmune diseases, and metabolic disorders. This provides an opportunity to select wheat cultivars with significantly lower levels of ATIs.

Sequence characterization and N-glycoproteomics of secretory immunoglobulin A from donkey milk

Secretory immunoglobulin A (sIgA) is the major antibody present in the human milk where it confers passive immunity to neonates. Other than human, non-ruminants such as equine, swine etc., also possess sIgA in milk but detailed characterization is limited. In the present study, we characterized sIgA from donkey milk for amino acid sequence and N-glycosylation through LC-MS/MS analysis. The complete amino acid sequence of alpha chain constant region (CH) was elucidated. The sequence analysis of variable regions (VH and VL) and light chain constant region (CL) showed several amino acid substitutions indicating the presence of diverse immunoglobulin repertoire. Glycoproteomic analysis of secretory component revealed bi-antennary complex and hybrid types with differential core fucosylation at site N83LT, only complex glycans at N135GT, N423GT and N530LT with mainly NeuAc whereas N291QT harbors high mannose glycans. Heavy chain possesses majorly bi-antennary complex with differential core fucosylation at sites N139AS and N338VS, in which N338VS shows partial occupancy. Joining chain harbors only complex type at N72IS, with core fucosylation and terminal NeuGc to some extent. N-glycan repertoire in part is similar to human sIgA. This comprehensive analysis of sequence and glycan pattern of donkey milk sIgA would be beneficial for its potential applications.

In Vitro Selection of New DNA Aptamers for Human Vascular Endothelial Growth Factor 165

Two DNA aptamers that bind the heparin-binding domain (HBD) of the human vascular endothelial growth factor 165 (VEGF-165) have been previously reported. Although VEGF-165 is a homodimeric protein and the two aptamers have different sequences and secondary structures, the aptamers appear to occupy the same binding site and cannot form a 2 : 1 aptamer/protein complex, thus making them unsuitable for creating a higher-affinity dimeric DNA aptamer. This has motivated us to conduct a new in vitro selection experiment to search for new VEGF-165-binding DNA aptamers with different properties. We undertook a multistream selection strategy in which the concentration of VEGF-165 was varied significantly. We carried out 11 rounds of selection, and next-generation sequencing was conducted for every round in each stream. From comprehensive sequence analysis, we identified four classes of DNA aptamers, of which two were reported before, but two are new DNA aptamers. One of the new aptamers exhibits a unique property that has never been observed before: it is capable of forming the 2 : 1 aptamer/protein complex with VEGF-165. This work has expanded the repertoire of VEGF-165-binding DNA aptamers and creates a possibility to engineer a higher affinity homodimeric aptamer for VEGF-165.

Indian hedgehog in craniofacial neural crest cells links to skeletal malocclusion by regulating associated cartilage formation and gene expression.

The frontal craniofacial skeleton derived from neural crest cells is vital for facial structure and masticatory functions. The exact role of Indian hedgehog (Ihh) in facial and masticatory development has not been fully explored. In this study, we generated craniofacial neural crest cells-specific Ihh deletion mice (Wnt1-Cre;Ihhfl/fl ;Tomatofl/+ ) and found the gradual dwarfism without perinatal lethality. Morphological and histological analyses revealed unambiguous craniofacial phenotypes in mutants, where we observed skeletal malocclusion accompanied by markedly hypoplastic nasomaxillary complex and reversed incisor occlusion. Both the replacement of nasal concha cartilage by turbinate bones and the endochondral ossification of nasal septum ethmoid bone were substantially delayed. We also observed hypoplastic mandibles in mutants where the mandibular ramus was unexpectedly the most affected. Both the condylar process and mandibular angle cartilages were distorted. However, dental examination showed no significant changes in teeth and dentition. Finally, a comprehensive RNA sequence analysis utilizing condylar cartilage identified Ihh-associated gene network including several cell cycle genes and 16 genes related to the extracellular matrix, sulfate transporters, transcription factors, receptors, a ciliogenesis factor, and an adhesion molecule. Our data provide direct in vivo evidence that Ihh plays crucial roles in midface and masticatory system formation, likely by activating key genes.

Lipidomics and RNA sequencing reveal a novel subpopulation of nanovesicle within extracellular matrix biomaterials.

Biomaterials composed of extracellular matrix (ECM) provide both mechanical support and a reservoir of constructive signaling molecules that promote functional tissue repair. Recently, matrix-bound nanovesicles (MBVs) have been reported as an integral component of ECM bioscaffolds. Although liquid-phase extracellular vesicles (EVs) have been the subject of intense investigation, their similarity to MBV is limited to size and shape. Liquid chromatography-mass spectrometry (LC-MS)-based lipidomics and redox lipidomics were used to conduct a detailed comparison of liquid-phase EV and MBV phospholipids. Combined with comprehensive RNA sequencing and bioinformatic analysis of the intravesicular cargo, we show that MBVs are a distinct and unique subpopulation of EV and a distinguishing feature of ECM-based biomaterials. The results begin to identify the differential biologic activities mediated by EV that are secreted by tissue-resident cells and deposited within the ECM.

Comprehensive Sequence Analysis of IQD Gene Family and their Expression Profiling in Grapevine (Vitis vinifera).

The plant-specific IQ67-domain (IQD) protein family members are downstream targets of calcium sensors, known to regulate plant growth and lateral organ polarity, and basal defense response against environmental cues. No systematic study of IQD gene family has been performed on grapevine. The public availability of grapevine genome enables us to perform identification, phylogeny, chromosomal orientation, and gene structure analysis of the IQD genes in grapevine. We identified 49 VvIQD genes (VvIQD1–VvIQD49) and further classified them into eight subgroups based on phylogenetic relationships. The 49 VvIQD genes were assigned to 19 different chromosomal positions. The collinear relationship between grapevine and Arabidopsis IQDs (VvIQD and AtIQD), and within grapevine VvIQDs, was highly conserved. In addition, most of duplicated gene pairs showed Ka/Ks ratio less than 1.00, indicating purifying selection within these gene pairs, implying functional discrepancy after duplication. Transcription profiling of VvIQD genes shed light on their specific role in grapevine tissue and organ development. The qRT-PCR validation of the 49 VvIQD genes in grape berry tissue from cultivars with distinct berry shape during developmental phases suggested candidate genes involved in the shape of grape berries. The subcellular prediction of VvIQD22, VvIQD23, VvIQD38, and VvIQD49 genes validated their localization in the nucleus and plasma membrane. The VvIQD49 protein interaction with VvCaM2 was also verified by bimolecular fluorescence complementation (BiFC) analysis in the plasma membrane. Our findings will be valuable for the functional genomic studies for desirable shape development of grape berries.

Cross-species transmission of the newly identified coronavirus 2019-nCoV.

The current outbreak of viral pneumonia in the city of Wuhan, China, was caused by a novel coronavirus designated 2019‐nCoV by the World Health Organization, as determined by sequencing the viral RNA genome. Many initial patients were exposed to wildlife animals at the Huanan seafood wholesale market, where poultry, snake, bats, and other farm animals were also sold. To investigate possible virus reservoir, we have carried out comprehensive sequence analysis and comparison in conjunction with relative synonymous codon usage (RSCU) bias among different animal species based on the 2019‐nCoV sequence. Results obtained from our analyses suggest that the 2019‐nCoV may appear to be a recombinant virus between the bat coronavirus and an origin‐unknown coronavirus. The recombination may occurred within the viral spike glycoprotein, which recognizes a cell surface receptor. Additionally, our findings suggest that 2019‐nCoV has most similar genetic information with bat coronovirus and most similar codon usage bias with snake. Taken together, our results suggest that homologous recombination may occur and contribute to the 2019‐nCoV cross‐species transmission.

Repertoire analysis of memory T-cell receptors in Japanese patients with inflammatory bowel disease.

Background and AimThe T‐cell receptor (TCR) repertoire was assessed in response to various antigens and was considered to be associated with the pathogenesis of inflammatory bowel disease (IBD). Thus, we performed TCR repertoire analysis to examine the pathology of IBD from changes in the TCR repertoire of memory T cells in the intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs) of patients with IBD.MethodsLPMCs in the surgical specimens and PBMCs were isolated from 12 patients with IBD (5 patients with ulcerative colitis [UC] and 7 patients with Crohn's disease [CD]). PBMCs were collected from 10 healthy individuals as controls. Comprehensive TCR sequence analyses of adaptor‐ligation polymerase chain reaction (PCR) products were performed using MiSeq.ResultsThe diversity of TCR‐α and TCR‐β in PBMCs was significantly lower in patients with IBD than that in controls (P = 0.00084 and 0.0013, respectively). Comparisons of TCR diversity in LPMCs and PBMCs between CD and UC showed that the diversity in LPMC was not affected by diseases, whereas that in PBMCs was significantly lower in CD than in UC (P = 0.045 and 0.049, respectively). Some TCR clones may have shown a specific increase or decrease in CD and UC, and many clones were common to both LPMCs and PBMCs in the same patients.ConclusionThe diversity of TCR clones in LPMCs and PBMCs in patients with IBD was significantly lower than that of PBMCs in controls. TCR diversity in PBMCs was particularly low in patients with CD.

Genetic and evolutionary analysis of enterovirus 71 base dinucleotide.

In this paper, a comprehensive analysis of genome sequences of 99 EV71 virus isolates base dinucleotide odds ratio was conducted to explore the potential mechanisms of its evolution. The general correlation of dinucleotides at different sites of the codons reveals that mutational pressure is the main factor determining EV 71 evolution. Furthermore, the apparent geographic segmentation of isolates before 2008 and low-frequency appearance of TpA and CpG also suggest that the geographic and host factors play important roles in the evolution of the virus. Finally, obvious fluctuation and divergence of the dinucleotide usage pattern was observed in the EV71 isolated in China after 2008, which might reveal that the virus entered a rapid evolution period, especially after 2010. These results indicate that EV71 has more initiative and adaptability under external pressures.

P027 Repertoire analysis of memory T-cell receptors in Japanese patients with inflammatory bowel disease

Abstract Background The T-cell receptor (TCR) repertoire was assessed in response to various antigens and was considered to be associated with the pathogenesis of inflammatory bowel disease (IBD). Thus, we performed TCR repertoire analysis to examine the pathology of IBD from changes in the TCR repertoire of memory T cells in the intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs) of patients with IBD. Methods LPMCs in the surgical specimens and PBMCs were isolated from 12 patients with IBD [five patients with ulcerative colitis (UC) and seven patients with Crohn’s disease (CD)]. PBMCs were collected from 10 healthy individuals as controls. Comprehensive TCR sequence analyses of adaptor-ligation PCR products were performed using MiSeq. Results The diversity of TCR-α and -β in PBMCs was significantly lower in patients with IBD than that in controls (p = 0.00084 and 0.0013, respectively). Comparisons of TCR diversity in LPMCs and PBMCs between CD and UC showed that the diversity in LPMC was not affected by diseases, whereas that in PBMCs was significantly lower in CD than in UC (p = 0.045 and 0.049, respectively). Some TCR clones may have shown a specific increase or decrease in CD and UC and many clones were common to both LPMCs and PBMCs in the same patients. Conclusion The diversity of TCR clones in LPMCs and PBMCs in patients with IBD was significantly lower than that of PBMCs in controls. TCR diversity in PBMCs was particularly low in patients with CD.

What is the role of the thioredoxin antioxidant complex in relation to HAM/TSP?

We have little information about the definite role of the thioredoxin antioxidant complex system during viral infection, particularly during human T-cell lymphotropic virus type 1 (HTLV-1) infection and the HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) state. Therefore, we conducted comprehensive next-generation sequencing (NGS) analysis to determine Trx system expression changes in three categories of subjects: sero-negative HTLV-1 individuals, asymptomatic HTLV-1 people and HAM/TSP patients. We found that Trx capacity is reduced in the HAM/TSP state compared to healthy individuals, which indicates increasing inflammation and reduction of apoptosis, which might contribute to the progression of inflammation in the spinal cord, which in turn may develop into the HAM/TSP state.

The proinflammatory cytokines TNF-α and IL-6 in lumpfish (Cyclopterus lumpus L.) -identification, molecular characterization, phylogeny and gene expression analyses

The proinflammatory cytokines TNF-α and IL-6 are important mediators of inflammatory reactions and orchestrators of the immune system in vertebrate. In this study, we have identified TNF-α and IL-6 in lumpfish, molecular characterized them at mRNA and gene level, performed homology modelling and measured their gene expression in different tissues and upon in vitro stimulation. A comprehensive phylogenetic analysis of TNF-α teleost sequences give novel insight into the TNF –α biology. Interestingly, we identified two isoforms of luIL-6. In normal tissue and leukocyte, the level of luTNF-α transcripts was higher than luIL-6. The expression pattern were parallel, except for brain, eye and gonad, and they displayed a similar induction pattern upon exposure to PAMPs, being most highly upregulated by flagellin. This is the first in-depth characterization of TNF and IL-6 in lumpfish. In recent years, lumpfish has become an important species for the aquaculture industry and establishment of qPCR-assays of luTNF-α and luIL-6 provide a valuable tool to measure effect of immune modulation, such as vaccination, microbiological disease and physiological trials. Lumpfish is also interesting for comparative studies as it represent a phylogenetic group that is poorly described immunologically.

Fatty Acid Amide Hydrolases: An Expanded Capacity for Chemical Communication?

Fatty acid amide hydrolase (FAAH) is an enzyme that belongs to the amidase signature (AS) superfamily and is widely distributed in multicellular eukaryotes. FAAH hydrolyzes lipid signaling molecules - namely, N-acylethanolamines (NAEs) - which terminates their actions. Recently, the crystal structure of Arabidopsis thaliana FAAH was solved and key residues were identified for substrate-specific interactions. Here, focusing on residues surrounding the substrate-binding pocket, a comprehensive analysis of FAAH sequences from angiosperms reveals a distinctly different family of FAAH-like enzymes. We hypothesize that FAAH, in addition to its role in seedling development, also acts in an N-acyl amide communication axis to facilitate plant-microbe interactions and that structural diversity provides for the flexible use of a wide range of small lipophilic signaling molecules.

Unlocking COVID therapeutic targets: A structure-based rationale against SARS-CoV-2, SARS-CoV and MERS-CoV Spike

There are no approved target therapeutics against SARS-CoV-2 or other beta-CoVs. The beta-CoV Spike protein is a promising target considering the critical role in viral infection and pathogenesis and its surface exposed features. We performed a structure-based strategy targeting highly conserved druggable regions resulting from a comprehensive large-scale sequence analysis and structural characterization of Spike domains across SARSr- and MERSr-CoVs. We have disclosed 28 main consensus druggable pockets within the Spike. The RBD and SD1 (S1 subunit); and the CR, HR1 and CH (S2 subunit) represent the most promising conserved druggable regions. Additionally, we have identified 181 new potential hot spot residues for the hSARSr-CoVs and 72 new hot spot residues for the SARSr- and MERSr-CoVs, which have not been described before in the literature. These sites/residues exhibit advantageous structural features for targeted molecular and pharmacological modulation. This study establishes the Spike as a promising anti-CoV target using an approach with a potential higher resilience to resistance development and directed to a broad spectrum of Beta-CoVs, including the new SARS-CoV-2 responsible for COVID-19. This research also provides a structure-based rationale for the design and discovery of chemical inhibitors, antibodies or other therapeutic modalities successfully targeting the Beta-CoV Spike protein.

Comprehensive Analysis of Pre-Charge Sequence in Automotive Battery Systems

<p>Electric vehicles (EV) have brought promising technologies for future mobility solutions. As one of the key components of EVs, battery systems have fundamental functions which disconnect the battery during parking and in case of failure. To provide a safe system, specialized high voltage (HV) electromechanical switches are used to perform these major functions such as switch on, switch off or pre-charging. Due to these components can be easily damaged, expensive, heavy and bulky, a solution based on pure semiconductors may be desired to accomplish these operations. Many studies were exhibited on EV battery systems regarding developing solid-state systems for HV switchgear. Developing technology on semiconductor devices allows to make a safety concept based on only solidstate components. This study presents a comprehensive analysis off pre-charge sequences between conventional and semiconductor switchgear to be used in electric vehicle battery systems. Spice simulations are presented to investigate advantages and drawbacks of these systems.</p>

Single-molecule real-time sequencing identifies massive full-length cDNAs and alternative-splicing events that facilitate comparative and functional genomics study in the hexaploid crop sweet potato.

BackgroundSweet potato (Ipomoea batatas (L.) Lam.) is one of the most important crops in many developing countries and provides a candidate source of bioenergy. However, neither a complete reference genome nor large-scale full-length cDNA sequences for this outcrossing hexaploid crop are available, which in turn impedes progress in research studies in I. batatas functional genomics and molecular breeding.MethodsIn this study, we sequenced full-length transcriptomes in I. batatas and its diploid ancestor I. trifida by single-molecule real-time sequencing and Illumina second-generation sequencing technologies. With the generated datasets, we conducted comprehensive intraspecific and interspecific sequence analyses and experimental characterization.ResultsA total of 53,861/51,184 high-quality long-read transcripts were obtained, which covered about 10,439/10,452 loci in the I. batatas/I. trifida genome. These datasets enabled us to predict open reading frames successfully in 96.83%/96.82% of transcripts and identify 34,963/33,637 full-length cDNA sequences, 1,401/1,457 transcription factors, 25,315/27,090 simple sequence repeats, 1,656/1,389 long non-coding RNAs, and 5,251/8,901 alternative splicing events. Approximately, 32.34%/38.54% of transcripts and 46.22%/51.18% multi-exon transcripts underwent alternative splicing in I. batatas/I. trifida. Moreover, we validated one alternative splicing event in each of 10 genes and identified tuberous-root-specific expressed isoforms from a starch-branching enzyme, an alpha-glucan phosphorylase, a neutral invertase, and several ABC transporters. Overall, the collection and analysis of large-scale long-read transcripts generated in this study will serve as a valuable resource for the I. batatas research community, which may accelerate the progress in its structural, functional, and comparative genomics studies.

TMOD-25. GLIOBLASTOMA ORGANOIDS: A MODEL SYSTEM FOR PATIENT-SPECIFIC THERAPEUTIC TESTING

Abstract Glioblastoma treatment options remain limited due to its aggressive and invasive nature. It is increasingly appreciated that molecular heterogeneity between tumors and within tumors likely contributes to the lack of therapeutic advances. To maintain the inherent heterogeneity of glioblastoma, we employed a novel method to rapidly culture glioblastoma organoids (GBOs) directly from neurosurgical resection. GBOs are routinely generated around two weeks following initial resection. Comprehensive histologic and sequencing analyses demonstrated similarity to primary tumors. Leveraging clinical molecular and sequencing data, selected GBOs were treated with radiation/temozolamide and targeted inhibitor therapies. The effect on proliferation was measured by the percentage of KI67+ cells and gene set enrichment (GSEA) analysis was performed to compare the pre-treated expression signature amongst responsive and non-responsive tumors. Treatment of organoids with radiation/temozolamide led to a decrease in the percentage of KI67+ cells in four of eight patient-derived organoid lines with some evidence of correlative radiographic response Gene sets associated with radiation response and TNF signaling were enriched in radiation/temozolamide sensitive GBOs. GBO response to EGFR inhibition via gefitinib treatment was specific to EGFR altered tumors, whose expression also enriched for EGF signaling pathway expression. Two GBOs had downstream NF1 mutations that responded to the MEK inhibitor trametinib. On GSEA, gene expression of NF1 mutated GBOs enriched for RAS signaling. One GBO line was found to have a PI3K mutation and responded dramatically to mTOR inhibition via everolimus. Dichotomous efficacy of MEK and mTOR inhibition was also noted by tumor-specific changes in GBO diameter following treatment. This novel culturing method of GBOs maintains intertumoral and intratumoral heterogeneity and allows for therapeutic testing within two weeks of neurosurgical resection. As clinical sequencing because increasingly prevalent, GBOs may become a valuable tool to functionally test mutation-specific treatment strategies in a patient-specific manner within a clinically relevant time frame.

Cervical cancer patients that respond to chemoradiation therapy display an intense tumor infiltrating immune profile before treatment

Cervical cancer (CC) is a major cause of death and suffering to women globally with 570,000 new cases in 2017. It disproportionately affects those living in resource-constrained countries such as Brazil, with 90% of the deaths from CC happening in low and middle-income countries. Early detection is still the best strategy for improving response to therapy and survival and cases detected in advanced stages show variable response rates to the standard chemoradiation therapy protocols. Both the genetic landscape and the immune status of patients can dramatically affect cancer progression and response to therapy, as well as disease recurrence. Here we performed a comprehensive sequencing analysis using the cancer gene panel - Ion AmpliSeq™ Cancer hotspot Panel V2 CHPv2, as well as determined the immune infiltrate composition of a group of locally advanced CC patients with the goal of identifying genetic and immune characteristics associated with a clinical response to therapy. The expression levels of CD68+ tumor-associated macrophages (TAMs) and CD8+ tumor-infiltrating lymphocytes (TILs), as well as the immune checkpoint molecules PD-1, PD-L1 and PD-L2 in stroma and in tumor regions were analyzed by immunohistochemistry (IHC). The HPV infection status with high-risk strains was also determined. Twenty-one samples from patients with squamous cell carcinoma segregated into responder (11) and non-responder (10) groups according to standard chemoradiation therapy response were studied. Our findings indicate that responder patients showed an increase of an inflammatory tumor microenvironment as indicated by higher numbers of CD8+ and PD-L2+ TILs, as well as higher expression of PD-L1 immunoreactive area, as compared to the non-responder group. Additionally, our results demonstrate a correlation between the number of gene mutations and PD-L2+ TILs in the responder group. The genes PIK3CA and KDR/VEGFR were the most mutated genes, corroborating past findings. Together, these findings indicate an inflammatory tumor microenvironment present in patients that will respond to future chemoradiation treatment as compared to those that will not. This points to possible future predictors of response to therapy in CC patients.