Fundamental immune-oncogenicity trade-offs define driver mutation fitness.

Missense driver mutations in cancer are concentrated in a few hotspots1. Various mechanisms have been proposed to explain this skew, including biased mutational processes2, phenotypic differences3-6 and immunoediting of neoantigens7,8; however, to our knowledge, no existing model weighs the relative contribution of these features to tumour evolution. We propose a unified theoretical 'free fitness' framework that parsimoniously integrates multimodal genomic, epigenetic, transcriptomic and proteomic data into a biophysical model of the rate-limiting processes underlying the fitness advantage conferred on cancer cells by driver gene mutations. Focusing on TP53, the most mutated gene in cancer1, we present an inference of mutant p53 concentration and demonstrate that TP53 hotspot mutations optimally solve an evolutionary trade-off between oncogenic potential and neoantigen immunogenicity. Our model anticipates patient survival in The Cancer Genome Atlas and patients with lung cancer treated with immunotherapy as well as the age of tumour onset in germline carriers of TP53 variants. The predicted differential immunogenicity between hotspot mutations was validated experimentally in patients with cancer and in a unique large dataset of healthy individuals. Our data indicate that immune selective pressure on TP53 mutations has a smaller role in non-cancerous lesions than in tumours, suggesting that targeted immunotherapy may offer an early prophylactic opportunity for the former. Determining the relative contribution of immunogenicity and oncogenic function to the selective advantage of hotspot mutations thus has important implications for both precision immunotherapies and our understanding of tumour evolution.

Magnitude and Dynamics of the T-Cell Response to SARS-CoV-2 Infection at Both Individual and Population Levels.

T cells are involved in the early identification and clearance of viral infections and also support the development of antibodies by B cells. This central role for T cells makes them a desirable target for assessing the immune response to SARS-CoV-2 infection. Here, we combined two high-throughput immune profiling methods to create a quantitative picture of the T-cell response to SARS-CoV-2. First, at the individual level, we deeply characterized 3 acutely infected and 58 recovered COVID-19 subjects by experimentally mapping their CD8 T-cell response through antigen stimulation to 545 Human Leukocyte Antigen (HLA) class I presented viral peptides (class II data in a forthcoming study). Then, at the population level, we performed T-cell repertoire sequencing on 1,815 samples (from 1,521 COVID-19 subjects) as well as 3,500 controls to identify shared "public" T-cell receptors (TCRs) associated with SARS-CoV-2 infection from both CD8 and CD4 T cells. Collectively, our data reveal that CD8 T-cell responses are often driven by a few immunodominant, HLA-restricted epitopes. As expected, the T-cell response to SARS-CoV-2 peaks about one to two weeks after infection and is detectable for at least several months after recovery. As an application of these data, we trained a classifier to diagnose SARS-CoV-2 infection based solely on TCR sequencing from blood samples, and observed, at 99.8% specificity, high early sensitivity soon after diagnosis (Day 3-7 = 85.1% [95% CI = 79.9-89.7]; Day 8-14 = 94.8% [90.7-98.4]) as well as lasting sensitivity after recovery (Day 29+/convalescent = 95.4% [92.1-98.3]). These results demonstrate an approach to reliably assess the adaptive immune response both soon after viral antigenic exposure (before antibodies are typically detectable) as well as at later time points. This blood-based molecular approach to characterizing the cellular immune response has applications in clinical diagnostics as well as in vaccine development and monitoring.

A copy number variant is associated with a spectrum of pigmentation patterns in the rock pigeon (Columba livia).

Rock pigeons (Columba livia) display an extraordinary array of pigment pattern variation. One such pattern, Almond, is characterized by a variegated patchwork of plumage colors that are distributed in an apparently random manner. Almond is a sex-linked, semi-dominant trait controlled by the classical Stipper (St) locus. Heterozygous males (ZStZ+ sex chromosomes) and hemizygous Almond females (ZStW) are favored by breeders for their attractive plumage. In contrast, homozygous Almond males (ZStZSt) develop severe eye defects and often lack plumage pigmentation, suggesting that higher dosage of the mutant allele is deleterious. To determine the molecular basis of Almond, we compared the genomes of Almond pigeons to non-Almond pigeons and identified a candidate St locus on the Z chromosome. We found a copy number variant (CNV) within the differentiated region that captures complete or partial coding sequences of four genes, including the melanosome maturation gene Mlana. We did not find fixed coding changes in genes within the CNV, but all genes are misexpressed in regenerating feather bud collar cells of Almond birds. Notably, six other alleles at the St locus are associated with depigmentation phenotypes, and all exhibit expansion of the same CNV. Structural variation at St is linked to diversity in plumage pigmentation and gene expression, and thus provides a potential mode of rapid phenotypic evolution in pigeons.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon.

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Plant organ evolution revealed by phylotranscriptomics in Arabidopsis thaliana.

The evolution of phenotypes occurs through changes both in protein sequence and gene expression levels. Though much of plant morphological evolution can be explained by changes in gene expression, examining its evolution has challenges. To gain a new perspective on organ evolution in plants, we applied a phylotranscriptomics approach. We combined a phylostratigraphic approach with gene expression based on the strand-specific RNA-seq data from seedling, floral bud, and root of 19 Arabidopsis thaliana accessions to examine the age and sequence divergence of transcriptomes from these organs and how they adapted over time. Our results indicate that, among the sense and antisense transcriptomes of these organs, the sense transcriptomes of seedlings are the evolutionarily oldest across all accessions and are the most conserved in amino acid sequence for most accessions. In contrast, among the sense transcriptomes from these same organs, those from floral bud are evolutionarily youngest and least conserved in sequence for most accessions. Different organs have adaptive peaks at different stages in their evolutionary history; however, all three show a common adaptive signal from the Magnoliophyta to Brassicale stage. Our research highlights how phylotranscriptomic analyses can be used to trace organ evolution in the deep history of plant species.

Immunosequencing identifies signatures of cytomegalovirus exposure history and HLA-mediated effects on the T cell repertoire.

An individual's T cell repertoire dynamically encodes their pathogen exposure history. To determine whether pathogen exposure signatures can be identified by documenting public T cell receptors (TCRs), we profiled the T cell repertoire of 666 subjects with known cytomegalovirus (CMV) serostatus by immunosequencing. We developed a statistical classification framework that could diagnose CMV status from the resulting catalog of TCRß sequences with high specificity and sensitivity in both the original cohort and a validation cohort of 120 different subjects. We also confirmed that three of the identified CMV-associated TCRß molecules bind CMV in vitro, and, moreover, we used this approach to accurately predict the HLA-A and HLA-B alleles of most subjects in the first cohort. As all memory T cell responses are encoded in the common format of somatic TCR recombination, our approach could potentially be generalized to a wide variety of disease states, as well as other immunological phenotypes, as a highly parallelizable diagnostic strategy.

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits.

To understand the population genetics of structural variants and their effects on phenotypes, we developed an approach to mapping structural variants that segregate in a population sequenced at low coverage. We avoid calling structural variants directly. Instead, the evidence for a potential structural variant at a locus is indicated by variation in the counts of short-reads that map anomalously to that locus. These structural variant traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between a structural variant trait at one locus, and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3×) population sequence data from 488 recombinant inbred Arabidopsis thaliana genomes, we identified 6502 segregating structural variants. Remarkably, 25% of these were transpositions. While many structural variants cannot be delineated precisely, we validated 83% of 44 predicted transposition breakpoints by polymerase chain reaction. We show that specific structural variants may be causative for quantitative trait loci for germination and resistance to infection by the fungus Albugo laibachii, isolate Nc14. Further we show that the phenotypic heritability attributable to read-mapping anomalies differs from, and, in the case of time to germination and bolting, exceeds that due to standard genetic variation. Genes within structural variants are also more likely to be silenced or dysregulated. This approach complements the prevalent strategy of structural variant discovery in fewer individuals sequenced at high coverage. It is generally applicable to large populations sequenced at low-coverage, and is particularly suited to mapping transpositions.

Limited Contribution of DNA Methylation Variation to Expression Regulation in Arabidopsis thaliana.

The extent to which epigenetic variation affects complex traits in natural populations is not known. We addressed this question using transcriptome and DNA methylation data from a sample of 135 sequenced A. thaliana accessions. Across individuals, expression was significantly associated with cis-methylation for hundreds of genes, and many of these associations remained significant after taking SNP effects into account. The pattern of correlations differed markedly between gene body methylation and transposable element methylation. The former was usually positively correlated with expression, and the latter usually negatively correlated, although exceptions were found in both cases. Finally, we developed graphical models of causality that adapt to a sample with heavy population structure, and used them to show that while methylation appears to affect gene expression more often than expression affects methylation, there is also strong support for both being independently controlled. In conclusion, although we find clear evidence for epigenetic regulation, both the number of loci affected and the magnitude of the effects appear to be small compared to the effect of SNPs.

Wham: Identifying Structural Variants of Biological Consequence.

Existing methods for identifying structural variants (SVs) from short read datasets are inaccurate. This complicates disease-gene identification and efforts to understand the consequences of genetic variation. In response, we have created Wham (Whole-genome Alignment Metrics) to provide a single, integrated framework for both structural variant calling and association testing, thereby bypassing many of the difficulties that currently frustrate attempts to employ SVs in association testing. Here we describe Wham, benchmark it against three other widely used SV identification tools-Lumpy, Delly and SoftSearch-and demonstrate Wham's ability to identify and associate SVs with phenotypes using data from humans, domestic pigeons, and vaccinia virus. Wham and all associated software are covered under the MIT License and can be freely downloaded from github (https://github.com/zeeev/wham), with documentation on a wiki (http://zeeev.github.io/wham/). For community support please post questions to https://www.biostars.org/.

DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation.

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation.

High resolution genetic mapping uncovers chitin synthase-1 as the target-site of the structurally diverse mite growth inhibitors clofentezine, hexythiazox and etoxazole in Tetranychus urticae.

The acaricides clofentezine, hexythiazox and etoxazole are commonly referred to as 'mite growth inhibitors', and clofentezine and hexythiazox have been used successfully for the integrated control of plant mite pests for decades. Although they are still important today, their mode of action has remained elusive. Recently, a mutation in chitin synthase 1 (CHS1) was linked to etoxazole resistance. In this study, we identified and investigated a Tetranychus urticae strain (HexR) harboring recessive, monogenic resistance to each of hexythiazox, clofentezine, and etoxazole. To elucidate if there is a common genetic basis for the observed cross-resistance, we adapted a previously developed bulk segregant analysis method to map with high resolution a single, shared resistance locus for all three compounds. This finding indicates that the underlying molecular basis for resistance to all three compounds is identical. This locus is centered on the CHS1 gene, and as supported by additional genetic and biochemical studies, a non-synonymous variant (I1017F) in CHS1 associates with resistance to each of the tested acaricides in HexR. Our findings thus demonstrate a shared molecular mode of action for the chemically diverse mite growth inhibitors clofentezine, hexythiazox and etoxazole as inhibitors of an essential, non-catalytic activity of CHS1. Given the previously documented cross-resistance between clofentezine, hexythiazox and the benzyolphenylurea (BPU) compounds flufenoxuron and cycloxuron, CHS1 should be also considered as a potential target-site of insecticidal BPUs.

Reciprocal responses in the interaction between Arabidopsis and the cell-content-feeding chelicerate herbivore spider mite.

Most molecular-genetic studies of plant defense responses to arthropod herbivores have focused on insects. However, plant-feeding mites are also pests of diverse plants, and mites induce different patterns of damage to plant tissues than do well-studied insects (e.g. lepidopteran larvae or aphids). The two-spotted spider mite (Tetranychus urticae) is among the most significant mite pests in agriculture, feeding on a staggering number of plant hosts. To understand the interactions between spider mite and a plant at the molecular level, we examined reciprocal genome-wide responses of mites and its host Arabidopsis (Arabidopsis thaliana). Despite differences in feeding guilds, we found that transcriptional responses of Arabidopsis to mite herbivory resembled those observed for lepidopteran herbivores. Mutant analysis of induced plant defense pathways showed functionally that only a subset of induced programs, including jasmonic acid signaling and biosynthesis of indole glucosinolates, are central to Arabidopsis's defense to mite herbivory. On the herbivore side, indole glucosinolates dramatically increased mite mortality and development times. We identified an indole glucosinolate dose-dependent increase in the number of differentially expressed mite genes belonging to pathways associated with detoxification of xenobiotics. This demonstrates that spider mite is sensitive to Arabidopsis defenses that have also been associated with the deterrence of insect herbivores that are very distantly related to chelicerates. Our findings provide molecular insights into the nature of, and response to, herbivory for a representative of a major class of arthropod herbivores.

Population bulk segregant mapping uncovers resistance mutations and the mode of action of a chitin synthesis inhibitor in arthropods.

Because of its importance to the arthropod exoskeleton, chitin biogenesis is an attractive target for pest control. This point is demonstrated by the economically important benzoylurea compounds that are in wide use as highly specific agents to control insect populations. Nevertheless, the target sites of compounds that inhibit chitin biogenesis have remained elusive, likely preventing the full exploitation of the underlying mode of action in pest management. Here, we show that the acaricide etoxazole inhibits chitin biogenesis in Tetranychus urticae (the two-spotted spider mite), an economically important pest. We then developed a population-level bulk segregant mapping method, based on high-throughput genome sequencing, to identify a locus for monogenic, recessive resistance to etoxazole in a field-collected population. As supported by additional genetic studies, including sequencing across multiple resistant strains and genetic complementation tests, we associated a nonsynonymous mutation in the major T. urticae chitin synthase (CHS1) with resistance. The change is in a C-terminal transmembrane domain of CHS1 in a highly conserved region that may serve a noncatalytic but essential function. Our finding of a target-site resistance mutation in CHS1 shows that at least one highly specific chitin biosynthesis inhibitor acts directly to inhibit chitin synthase. Our work also raises the possibility that other chitin biogenesis inhibitors, such as the benzoylurea compounds, may also act by inhibition of chitin synthases. More generally, our genetic mapping approach should be powerful for high-resolution mapping of simple traits (resistance or otherwise) in arthropods.

Divergence, convergence, and the ancestry of feral populations in the domestic rock pigeon.

Domestic pigeons are spectacularly diverse and exhibit variation in more traits than any other bird species [1]. In The Origin of Species, Charles Darwin repeatedly calls attention to the striking variation among domestic pigeon breeds-generated by thousands of years of artificial selection on a single species by human breeders-as a model for the process of natural divergence among wild populations and species [2]. Darwin proposed a morphology-based classification of domestic pigeon breeds [3], but the relationships among major groups of breeds and their geographic origins remain poorly understood [4, 5]. We used a large, geographically diverse sample of 361 individuals from 70 domestic pigeon breeds and two free-living populations to determine genetic relationships within this species. We found unexpected relationships among phenotypically divergent breeds as well as convergent evolution of derived traits among several breed groups. Our findings also illuminate the geographic origins of breed groups in India and the Middle East and suggest that racing breeds have made substantial contributions to feral pigeon populations.

The genome of Tetranychus urticae reveals herbivorous pest adaptations.

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant-herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana.

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

The onset temperature of the heat-shock response and whole-organism thermal tolerance are tightly correlated in both laboratory-acclimated and field-acclimatized tidepool sculpins (Oligocottus maculosus).

We examined the relationship between thermal tolerance, measured as critical thermal maximum (CT(max)), and aspects of the heat-shock response in tidepool sculpins (Oligocottus maculosus) acclimated to constant laboratory temperatures or acclimatized to field conditions. The CT(max) of fish laboratory acclimated to 6°, 13°, and 20°C were 27.6° ± 0.1°C, 29.5° ± 0.1°C, and 30.8° ± 0.1°C, respectively, increasing linearly by 0.2°C for each 1°C increase in acclimation temperature. The CT(max) of field-acclimatized fish from the low intertidal (29.9° ± 0.1°C) was significantly lower than that of fish from the mid- (30.5° ± 0.1°C) and high (30.4° ± 0.1°C) intertidal. CT(max) and the onset temperature of hsp70 induction in gill (T(on)) were highly correlated in both laboratory-acclimated and field-acclimatized sculpins, with T(on) occurring at 2°C below CT(max) in all cases. However, there was no consistent relationship between CT(max) and the maximum levels of gill hsp70 mRNA. Predicted "acclimation" temperature (15.9° ± 0.3°C) and mean habitat temperature (15.9° ± 1.6°C) were similar for sculpins from low intertidal pools, but this relationship was not apparent in mid- and high intertidal fish. Mark-recapture experiments indicated that approximately 80% of fish from low intertidal pools were residents of that pool, but residency rates were less than 50% in mid- and high intertidal pools, which may explain the lack of correlation between CT(max) and habitat variables in these groups. These data indicate that gill hsp70 T(on) and CT(max) are highly correlated indicators of the thermal performance of tidepool sculpins in both laboratory and field settings.

Heat shock response of killifish (Fundulus heteroclitus): candidate gene and heterologous microarray approaches.

Northern and southern subspecies of the Atlantic killifish, Fundulus heteroclitus, differ in maximal thermal tolerance. To determine whether these subspecies also differ in their heat shock response (HSR), we exposed 20°C acclimated killifish to a 2 h heat shock at 34°C and examined gene expression in fish from both subspecies during heat shock and recovery using real-time quantitative PCR and a heterologous cDNA microarray designed for salmonid fishes. The heat shock proteins Hsp70-1, hsp27, and hsp30 were upregulated to a greater extent in the high temperature-tolerant southern subspecies than in the less tolerant northern subspecies, whereas hsp70-2 (which showed the largest upregulation of all the heat shock proteins) in both gill and muscle and hsp90α in muscle was upregulated to a greater extent in northern than in southern fish. These data demonstrate that differences in the HSR between subspecies cannot be due to changes in a single global regulator but must occur via gene-specific mechanisms. They also suggest that the role, if any, of hsps in establishing thermal tolerance is complex and varies from gene to gene. Heterologous microarray hybridization provided interpretable gene expression signatures, detecting differential regulation of genes known to be involved in the heat shock response in other species. Under control conditions, a variety of genes were differentially expressed in muscle between subspecies that suggest differences in muscle fiber type and could relate to previously observed differences between subspecies in the thermal sensitivity of swimming performance and metabolism.