Prior anti-CTLA-4 therapy impacts molecular characteristics associated with anti-PD-1 response in advanced melanoma.

Immune checkpoint inhibitors (ICIs), including CTLA-4- and PD-1-blocking antibodies, can have profound effects on tumor immune cell infiltration that have not been consistent in biopsy series reported to date. Here, we analyze seven molecular datasets of samples from patients with advanced melanoma (N = 514) treated with ICI agents to investigate clinical, genomic, and transcriptomic features of anti-PD-1 response in cutaneous melanoma. We find that prior anti-CTLA-4 therapy is associated with differences in genomic, individual gene, and gene signatures in anti-PD-1 responders. Anti-CTLA-4-experienced melanoma tumors that respond to PD-1 blockade exhibit increased tumor mutational burden, inflammatory signatures, and altered cell cycle processes compared with anti-CTLA-4-naive tumors or anti-CTLA-4-experienced, anti-PD-1-nonresponsive melanoma tumors. We report a harmonized, aggregate resource and suggest that prior CTLA-4 blockade therapy is associated with marked differences in the tumor microenvironment that impact the predictive features of PD-1 blockade therapy response.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB.

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Impact of checkpoint blockade on cancer vaccine-activated CD8+ T cell responses.

Immune and molecular profiling of CD8 T cells of patients receiving DC vaccines expressing three full-length melanoma antigens (MAs) was performed. Antigen expression levels in DCs had no significant impact on T cell or clinical responses. Patients who received checkpoint blockade before DC vaccination had higher baseline MA-specific CD8 T cell responses but no evidence for improved functional responses to the vaccine. Patients who showed the best clinical responses had low PD-1 expression on MA-specific T cells before and after DC vaccination; however, blockade of PD-1 during antigen presentation by DC had minimal functional impact on PD-1high MA-specific T cells. Gene and protein expression analyses in lymphocytes and tumor samples identified critical immunoregulatory pathways, including CTLA-4 and PD-1. High immune checkpoint gene expression networks correlated with inferior clinical outcomes. Soluble serum PD-L2 showed suggestive positive association with improved outcome. These findings show that checkpoint molecular pathways are critical for vaccine outcomes and suggest specific sequencing of vaccine combinations.

Detecting differential copy number variation between groups of samples.

We present a method to detect copy number variants (CNVs) that are differentially present between two groups of sequenced samples. We use a finite-state transducer where the emitted read depth is conditioned on the mappability and GC-content of all reads that occur at a given base position. In this model, the read depth within a region is a mixture of binomials, which in simulations matches the read depth more closely than the often-used negative binomial distribution. The method analyzes all samples simultaneously, preserving uncertainty as to the breakpoints and magnitude of CNVs present in an individual when it identifies CNVs differentially present between the two groups. We apply this method to identify CNVs that are recurrently associated with postglacial adaptation of marine threespine stickleback (Gasterosteus aculeatus) to freshwater. We identify 6664 regions of the stickleback genome, totaling 1.7 Mbp, which show consistent copy number differences between marine and freshwater populations. These deletions and duplications affect both protein-coding genes and cis-regulatory elements, including a noncoding intronic telencephalon enhancer of DCHS1 The functions of the genes near or included within the 6664 CNVs are enriched for immunity and muscle development, as well as head and limb morphology. Although freshwater stickleback have repeatedly evolved from marine populations, we show that freshwater stickleback also act as reservoirs for ancient ancestral sequences that are highly conserved among distantly related teleosts, but largely missing from marine stickleback due to recent selective sweeps in marine populations.

Dorsal spine evolution in threespine sticklebacks via a splicing change in MSX2A.

BACKGROUND: Dorsal spine reduction in threespine sticklebacks (Gasterosteus aculeatus) is a classic example of recurrent skeletal evolution in nature. Sticklebacks in marine environments typically have long spines that form part of their skeletal armor. Many derived freshwater populations have evolved shorter spines. Changes in spine length are controlled in part by a quantitative trait locus (QTL) previously mapped to chromosome 4, but the causative gene and mutations underlying the repeated evolution of this interesting skeletal trait have not been identified. RESULTS: Refined mapping of the spine length QTL shows that it lies near the MSX2A transcription factor gene. MSX2A is expressed in developing spines. In F1 marine × freshwater fish, the marine allele is preferentially expressed. Differences in expression can be attributed to splicing regulation. Due to the use of an alternative 5 ' splice site within the first exon, the freshwater allele produces greater amounts of a shortened, non-functional transcript and makes less of the full-length transcript. Sequence changes in the MSX2A region are shared by many freshwater fish, suggesting that repeated evolution occurs by reuse of a spine-reduction variant. To demonstrate the effect of full-length MSX2A on spine length, we produced transgenic freshwater fish expressing a copy of marine MSX2A. The spines of the transgenic fish were significantly longer on average than those of their non-transgenic siblings, partially reversing the reduced spine lengths that have evolved in freshwater populations. CONCLUSIONS: MSX2A is a major gene underlying dorsal spine reduction in freshwater sticklebacks. The gene is linked to a separate gene controlling bony plate loss, helping explain the concerted effects of chromosome 4 on multiple armor-reduction traits. The nature of the molecular changes provides an interesting example of morphological evolution occurring not through a simple amino acid change, nor through a change only in gene expression levels, but through a change in the ratio of splice products encoding both normal and truncated proteins.

Structure-activity relationships among DNA ligase inhibitors: Characterization of a selective uncompetitive DNA ligase I inhibitor.

In human cells, there are three genes that encode DNA ligase polypeptides with distinct but overlapping functions. Previously small molecule inhibitors of human DNA ligases were identified using a structure-based approach. Three of these inhibitors, L82, a DNA ligase I (LigI)-selective inhibitor, and L67, an inhibitor of LigI and DNA ligases III (LigIII), and L189, an inhibitor of all three human DNA ligases, have related structures that are composed of two 6-member aromatic rings separated by different linkers. Here we have performed a structure-activity analysis to identify determinants of activity and selectivity. The majority of the LigI-selective inhibitors had a pyridazine ring whereas the LigI/III- and LigIII-selective inhibitors did not. In addition, the aromatic rings in LigI-selective inhibitors had either arylhydrazone or acylhydrazone, but not vinyl linkers. Among the LigI-selective inhibitors, L82-G17 exhibited increased activity against and selectivity for LigI compared with L82. Notably. L82-G17 is an uncompetitive inhibitor of the third step of the ligation reaction, phosphodiester bond formation. Cells expressing LigI were more sensitive to L82-G17 than isogenic LIG1 null cells. Furthermore, cells lacking nuclear LigIIIα, which can substitute for LigI in DNA replication, were also more sensitive to L82-G17 than isogenic parental cells. Together, our results demonstrate that L82-G17 is a LigI-selective inhibitor with utility as a probe of the catalytic activity and cellular functions of LigI and provide a framework for the future design of DNA ligase inhibitors.

DNA ligases as therapeutic targets.

During DNA replication, DNA joining events link Okazaki fragments on the lagging strand. In addition, they are required to repair DNA single- and double-strand breaks and to complete repair events initiated by the excision of mismatched and damaged bases. In human cells, there are three genes encoding DNA ligases. These enzymes are ATP-dependent and contain a conserved catalytic region. Biophysical studies have shown that the catalytic region contains three domains that, in the absence of DNA, are in an extended conformation. When the catalytic region engages a DNA nick, it adopts a compact, ring structure around the DNA nick with each of the three domains contacting the DNA. Protein-protein interactions involving the regions flanking the conserved catalytic regions of human DNA ligases play a major role in directing these enzymes to participate in specific DNA transactions. Among the human LIG genes, the LIG3 gene is unique in that it encodes multiple DNA ligase polypeptides with different N- and C-termini. One of these polypeptides is targeted to mitochondria where it plays an essential role in the maintenance of the mitochondrial genome. In the nucleus, DNA ligases I, III and IV have distinct but overlapping functions in DNA replication and repair. Small molecule inhibitors of human DNA ligases have been identified using structure-based approaches. As expected, these inhibitors are cytotoxic and also potentiate the cytotoxicity of DNA damaging agents. The results of preclinical studies with human cancer cell lines and mouse models of human cancer suggest that DNA ligase inhibitors may have utility as anti-cancer agents.

DNA ligase I, the replicative DNA ligase.

Multiple DNA ligation events are required to join the Okazaki fragments generated during lagging strand DNA synthesis. In eukaryotes, this is primarily carried out by members of the DNA ligase I family. The C-terminal catalytic region of these enzymes is composed of three domains: a DNA binding domain, an adenylation domain and an OB-fold domain. In the absence of DNA, these domains adopt an extended structure but transition into a compact ring structure when they engage a DNA nick, with each of the domains contacting the DNA. The non-catalytic N-terminal region of eukaryotic DNA ligase I is responsible for the specific participation of these enzymes in DNA replication. This proline-rich unstructured region contains the nuclear localization signal and a PCNA interaction motif that is critical for localization to replication foci and efficient joining of Okazaki fragments. DNA ligase I initially engages the PCNA trimer via this interaction motif which is located at the extreme N-terminus of this flexible region. It is likely that this facilitates an additional interaction between the DNA binding domain and the PCNA ring. The similar size and shape of the rings formed by the PCNA trimer and the DNA ligase I catalytic region when it engages a DNA nick suggest that these proteins interact to form a double-ring structure during the joining of Okazaki fragments. DNA ligase I also interacts with replication factor C, the factor that loads the PCNA trimeric ring onto DNA. This interaction, which is regulated by phosphorylation of the non-catalytic N-terminus of DNA ligase I, also appears to be critical for DNA replication.