Coding and noncoding drivers of mantle cell lymphoma identified through exome and genome sequencing.

Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma (NHL) that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established, and the features that contribute to differences in clinical course remain limited. To extend our understanding of the biological pathways involved in this malignancy, we performed a large-scale genomic analysis of MCL using data from 51 exomes and 34 genomes alongside previously published exome cohorts. To confirm our findings, we resequenced the genes identified in the exome cohort in 191 MCL tumors, each having clinical follow-up data. We confirmed the prognostic association of TP53 and NOTCH1 mutations. Our sequencing revealed novel recurrent noncoding mutations surrounding a single exon of the HNRNPH1gene. In RNA-seq data from 103 of these cases, MCL tumors with these mutations had a distinct imbalance of HNRNPH1 isoforms. This altered splicing of HNRNPH1 was associated with inferior outcomes in MCL and showed a significant increase in protein expression by immunohistochemistry. We describe a functional role for these recurrent noncoding mutations in disrupting an autoregulatory feedback mechanism, thereby deregulating HNRNPH1 protein expression. Taken together, these data strongly imply a role for aberrant regulation of messenger RNA processing in MCL pathobiology.

The double-hit signature identifies double-hit diffuse large B-cell lymphoma with genetic events cryptic to FISH.

High-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/THs) include a group of diffuse large B-cell lymphomas (DLBCLs) with inferior outcomes after standard chemoimmunotherapy. We recently described a gene expression signature that identifies 27% of germinal center B-cell DLBCLs (GCB-DLBCLs) as having a double-hit-like expression pattern (DHITsig) and inferior outcomes; however, only half of these cases have both MYC and BCL2 translocations identifiable using standard breakapart fluorescence in situ hybridization (FISH). Here, 20 DHITsig+ GCB-DLBCLs apparently lacking MYC and/or BCL2 rearrangements underwent whole-genome sequencing. This revealed 6 tumors with MYC or BCL2 rearrangements that were cryptic to breakapart FISH. Copy-number analysis identified 3 tumors with MYC and 6 tumors with MIR17HG gains or amplifications, both of which may contribute to dysregulation of MYC and its downstream pathways. Focal deletions of the PVT1 promoter were observed exclusively among DHITsig+ tumors lacking MYC translocations; this may also contribute to MYC overexpression. These results highlight that FISH fails to identify all HGBL-DH/THs, while revealing a range of other genetic mechanisms potentially underlying MYC dysregulation in DHITsig+ DLBCL, suggesting that gene expression profiling is more sensitive for identifying the biology underlying poor outcomes in GCB-DLBCL.

Genomic divergence in a ring species complex.

Ring species provide particularly clear demonstrations of how one species can gradually evolve into two, but are rare in nature. In the greenish warbler (Phylloscopus trochiloides) species complex, a ring of populations wraps around Tibet. Two reproductively isolated forms co-exist in central Siberia, with a gradient of genetic and phenotypic characteristics through the southern chain of populations connecting them. Previous genetic evidence has proven inconclusive, however, regarding whether species divergence took place in the face of continuous gene flow and whether hybridization between the terminal forms of the ring ever occurred. Here we use genome-wide analyses to show that, although spatial patterns of genetic variation are currently mostly as expected of a ring species, historical breaks in gene flow have existed at more than one location around the ring, and the two Siberian forms have occasionally interbred. Substantial periods of geographical isolation occurred not only in the north but also in the western Himalayas, where there is now an extensive hybrid zone between genetically divergent forms. Limited asymmetric introgression has occurred directly between the Siberian forms, although it has not caused a blending of those forms, suggesting selection against introgressed genes in the novel genetic background. Levels of reproductive isolation and genetic introgression are consistent with levels of phenotypic divergence around the ring, with phenotypic similarity and extensive interbreeding across the southwestern contact zone and strong phenotypic divergence and nearly complete reproductive isolation across the northern contact zone. These results cast doubt on the hypothesis that the greenish warbler should be viewed as a rare example of speciation by distance, but demonstrate that the greenish warbler displays a continuum from slightly divergent neighbouring populations to almost fully reproductively isolated species.

Ecology can inform genetics: Disassortative mating contributes to MHC polymorphism in Leach's storm-petrels (Oceanodroma leucorhoa).

Studies of MHC-based mate choice in wild populations often test hypotheses on species exhibiting female choice and male-male competition, which reflects the general prevalence of females as the choosy sex in natural systems. Here, we examined mutual mate-choice patterns in a small burrow-nesting seabird, the Leach's storm-petrel (Oceanodroma leucorhoa), using the major histocompatibility complex (MHC). The life history and ecology of this species are extreme: both partners work together to fledge a single chick during the breeding season, a task that requires regularly travelling hundreds of kilometres to and from foraging grounds over a 6- to 8-week provisioning period. Using a 5-year data set unprecedented for this species (n = 1078 adults and 925 chicks), we found a positive relationship between variation in the likelihood of female reproductive success and heterozygosity at Ocle-DAB2, a MHC class IIB locus. Contrary to previous reports rejecting disassortative mating as a mechanism for maintaining genetic polymorphism in this species, here we show that males make significant disassortative mate-choice decisions. Variability in female reproductive success suggests that the most common homozygous females (Ocle-DAB2*01/Ocle-DAB2*01) may be physiologically disadvantaged and, therefore, less preferred as lifelong partners for choosy males. The results from this study support the role of mate choice in maintaining high levels of MHC variability in a wild seabird species and highlight the need to incorporate a broader ecological framework and sufficient sample sizes into studies of MHC-based mating patterns in wild populations in general.

A comparison of genomic islands of differentiation across three young avian species pairs.

Detailed evaluations of genomic variation between sister species often reveal distinct chromosomal regions of high relative differentiation (i.e., "islands of differentiation" in FST ), but there is much debate regarding the causes of this pattern. We briefly review the prominent models of genomic islands of differentiation and compare patterns of genomic differentiation in three closely related pairs of New World warblers with the goal of evaluating support for the four models. Each pair (MacGillivray's/mourning warblers; Townsend's/black-throated green warblers; and Audubon's/myrtle warblers) consists of forms that were likely separated in western and eastern North American refugia during cycles of Pleistocene glaciations and have now come into contact in western Canada, where each forms a narrow hybrid zone. We show strong differences between pairs in their patterns of genomic heterogeneity in FST , suggesting differing selective forces and/or differing genomic responses to similar selective forces among the three pairs. Across most of the genome, levels of within-group nucleotide diversity (πWithin ) are almost as large as levels of between-group nucleotide distance (πBetween ) within each pair, suggesting recent common ancestry and/or gene flow. In two pairs, a pattern of the FST peaks having low πBetween suggests that selective sweeps spread between geographically differentiated groups, followed by local differentiation. This "sweep-before-differentiation" model is consistent with signatures of gene flow within the yellow-rumped warbler species complex. These findings add to our growing understanding of speciation as a complex process that can involve phases of adaptive introgression among partially differentiated populations.

Phase 2 study of panobinostat with or without rituximab in relapsed diffuse large B-cell lymphoma.

The majority of diffuse large B-cell lymphoma (DLBCL) tumors contain mutations in histone-modifying enzymes (HMEs), indicating a potential therapeutic benefit of histone deacetylase inhibitors (HDIs), and preclinical data suggest that HDIs augment the effect of rituximab. In this randomized phase 2 study, we evaluated the response rate and toxicity of panobinostat, a pan-HDI administered 30 mg orally 3 times weekly, with or without rituximab, in 40 patients with relapsed or refractory de novo (n = 27) or transformed (n = 13) DLBCL. Candidate genes and whole exomes were sequenced in relapse tumor biopsies to search for molecular correlates, and these data were used to quantify circulating tumor DNA (ctDNA) in serial plasma samples. Eleven of 40 patients (28%) responded to panobinostat (95% confidence interval [CI] 14.6-43.9) and rituximab did not increase responses. The median duration of response was 14.5 months (95% CI 9.4 to "not reached"). At time of data censoring, 6 of 11 patients had not progressed. Of the genes tested for mutations, only those in MEF2B were significantly associated with response. We detected ctDNA in at least 1 plasma sample from 96% of tested patients. A significant increase in ctDNA at day 15 relative to baseline was strongly associated with lack of response (sensitivity 71.4%, specificity 100%). We conclude that panobinostat induces very durable responses in some patients with relapsed DLBCL, and early responses can be predicted by mutations in MEF2B or a significant change in ctDNA level at 15 days after treatment initiation. This clinical trial was registered at www.ClinicalTrials.gov (#NCT01238692).

Genetic inactivation of TRAF3 in canine and human B-cell lymphoma.

Non-Hodgkin lymphomas (NHLs) are the most common cancer to affect pet dogs. In contrast to the many genes whose mutation contributes to lymphomagenesis in humans, relatively little is known about the acquired genetic alterations that lead to canine B-cell lymphomas (cBCLs). We performed a survey of 84 canine NHL tumors to identify genes affected by somatic point mutations. We found mutations affecting TRAF3, which encodes a negative regulator of nuclear factor (NF)-κB, to be a common feature of cBCLs, with mutations observed in 44% of tumors including a combination of somatic and rare germ-line variants. Overall, 30% of the tumors contained ≥1 somatic TRAF3 mutation. The majority of mutations are predicted to cause loss of TRAF3 protein including those impacting reading frame and splicing. To determine whether TRAF3 loss might be relevant to human NHL, we also analyzed 148 human diffuse large B-cell lymphoma (DLBCL) tumors and identified loss of TRAF3 as a common event, affecting ∼9% of DLBCLs, and reduced expression of TRAF3 among deleted cases. This study implicates mutations affecting NF-κB activity as a novel genetic commonality between human and canine NHLs and supports the potential utility of cBCLs with mutated TRAF3 as a model of the more aggressive activated B-cell subgroup of DLBCL.

Multiplex Droplet Digital PCR Quantification of Recurrent Somatic Mutations in Diffuse Large B-Cell and Follicular Lymphoma.

BACKGROUND: A plethora of options to detect mutations in tumor-derived DNA currently exist but each suffers limitations in analytical sensitivity, cost, or scalability. Droplet digital PCR (ddPCR) is an appealing technology for detecting the presence of specific mutations based on a priori knowledge and can be applied to tumor biopsies, including formalin-fixed paraffin embedded (FFPE) tissues. More recently, ddPCR has gained popularity in its utility in quantifying circulating tumor DNA. METHODS: We have developed a suite of novel ddPCR assays for detecting recurrent mutations that are prevalent in common B-cell non-Hodgkin lymphomas (NHLs), including diffuse large B-cell lymphoma, follicular lymphoma, and lymphoplasmacytic lymphoma. These assays allowed the differentiation and counting of mutant and wild-type molecules using one single hydrolysis probe. We also implemented multiplexing that allowed the simultaneous detection of distinct mutations and an "inverted" ddPCR assay design, based on employing probes matching wild-type alleles, capable of detecting the presence of multiple single nucleotide polymorphisms. RESULTS: The assays successfully detected and quantified somatic mutations commonly affecting enhancer of zeste 2 polycomb repressive complex 2 subunit (EZH2) (Y641) and signal transducer and activator of transcription 6 (STAT6) (D419) hotspots in fresh tumor, FFPE, and liquid biopsies. The "inverted" ddPCR approach effectively reported any single nucleotide variant affecting either of these 2 hotspots as well. Finally, we could effectively multiplex hydrolysis probes targeting 2 additional lymphoma-related hotspots: myeloid differentiation primary response 88 (MYD88; L265P) and cyclin D3 (CCND3; I290R). CONCLUSIONS: Our suite of ddPCR assays provides sufficient analytical sensitivity and specificity for either the invasive or noninvasive detection of multiple recurrent somatic mutations in B-cell NHLs.

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species.

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within-group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within-population selection-induced reductions in variation. An important implication of this 'sweep-before-differentiation' model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.

Digital PCR quantification of ultrahigh ERBB2 copy number identifies poor breast cancer survival after trastuzumab.

HER2/ERBB2 evaluation is necessary for treatment decision-making in breast cancer (BC), however current methods have limitations and considerable variability exists. DNA copy number (CN) evaluation by droplet digital PCR (ddPCR) has complementary advantages for HER2/ERBB2 diagnostics. In this study, we developed a single-reaction multiplex ddPCR assay for determination of ERBB2 CN in reference to two control regions, CEP17 and a copy-number-stable region of chr. 2p13.1, validated CN estimations to clinical in situ hybridization (ISH) HER2 status, and investigated the association of ERBB2 CN with clinical outcomes. 909 primary BC tissues were evaluated and the area under the curve for concordance to HER2 status was 0.93 and 0.96 for ERBB2 CN using either CEP17 or 2p13.1 as reference, respectively. The accuracy of ddPCR ERBB2 CN was 93.7% and 94.1% in the training and validation groups, respectively. Positive and negative predictive value for the classic HER2 amplification and non-amplification groups was 97.2% and 94.8%, respectively. An identified biological "ultrahigh" ERBB2 ddPCR CN group had significantly worse survival within patients treated with adjuvant trastuzumab for both recurrence-free survival (hazard ratio, HR: 3.3; 95% CI 1.1-9.6; p = 0.031, multivariable Cox regression) and overall survival (HR: 3.6; 95% CI 1.1-12.6; p = 0.041). For validation using RNA-seq data as a surrogate, in a population-based SCAN-B cohort (NCT02306096) of 682 consecutive patients receiving adjuvant trastuzumab, the ultrahigh-ERBB2 mRNA group had significantly worse survival. Multiplex ddPCR is useful for ERBB2 CN estimation and ultrahigh ERBB2 may be a predictive factor for decreased long-term survival after trastuzumab treatment.

Molecular evolution of the toll-like receptor multigene family in birds.

Toll-like receptors (TLR) are membrane-bound sensors of the innate immune system that recognize invariant and distinctive molecular features of invading microbes and are also essential for initiating adaptive immunity in vertebrates. The genetic variation at TLR genes has been directly related to differential pathogen outcomes in humans and livestock. Nonetheless, new insights about the impact of TLRs polymorphism on the evolutionary ecology of infectious diseases can be gained through the investigation of additional vertebrate groups not yet investigated in detail. In this study, we have conducted the first characterization of the entire TLR multigene family (N = 10 genes) in non-model avian species. Using primers targeting conserved coding regions, we aimed at amplifying large segments of the extracellular domains (275-435 aa) involved in pathogen recognition across seven phylogenetically diverse bird species. Our analyses suggest avian TLRs are dominated by stabilizing selection, suggesting that slow rates of nonsynonymous substitution help preserve biological function. Overall, mean values of ω (= d(n)/d(s)) at each TLR locus ranged from 0.196 to 0.517. However, we also found patterns of positive selection acting on specific amino acid sites that could be linked to species-specific differences in pathogen-associated molecular pattern recognition. Only 39 of 2,875 (∼1.35%) of the codons analyzed exhibited significant patterns of positive selection. At least one half of these positively selected codons can be mapped to putative ligand-binding regions, as suggested by crystallographic structures of TLRs and their ligands and mutagenic analyses. We also surveyed TLR polymorphism in wild populations of two bird species, the Lesser Kestrel Falco naumanni and the House Finch Carpodacus mexicanus. In general, avian TLRs displayed low to moderate single nucleotide polymorphism levels and an excess of silent nucleotide substitutions, but also conspicuous instances of positive directional selection. In particular, TLR5 and TLR15 exhibited the highest degree of genetic polymorphism and the highest occurrence of nonconservative amino acid substitutions. This study provides critical primers and a first look at the evolutionary patterns and implications of TLR polymorphism in non-model avian species and extends the list of candidate loci for avian eco-immunogenetics beyond the widely employed genes of the Major Histocompatibility Complex (MHC).

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird.

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.

Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection.

Insular populations have attracted the attention of evolutionary biologists because of their morphological and ecological peculiarities with respect to their mainland counterparts. Founder effects and genetic drift are known to distribute neutral genetic variability in these demes. However, elucidating whether these evolutionary forces have also shaped adaptive variation is crucial to evaluate the real impact of reduced genetic variation in small populations. Genes of the major histocompatibility complex (MHC) are classical examples of evolutionarily relevant loci because of their well-known role in pathogen confrontation and clearance. In this study, we aim to disentangle the partial roles of genetic drift and natural selection in the spatial distribution of MHC variation in insular populations. To this end, we integrate the study of neutral (22 microsatellites and one mtDNA locus) and MHC class II variation in one mainland (Iberia) and two insular populations (Fuerteventura and Menorca) of the endangered Egyptian vulture (Neophron percnopterus). Overall, the distribution of the frequencies of individual MHC alleles (n=17 alleles from two class II B loci) does not significantly depart from neutral expectations, which indicates a prominent role for genetic drift over selection. However, our results point towards an interesting co-evolution of gene duplicates that maintains different pairs of divergent alleles in strong linkage disequilibrium on islands. We hypothesize that the co-evolution of genes may counteract the loss of genetic diversity in insular demes, maximize antigen recognition capabilities when gene diversity is reduced, and promote the co-segregation of the most efficient allele combinations to cope with local pathogen communities.

Host-feeding patterns of native Culex pipiens and invasive Aedes albopictus mosquitoes (Diptera: Culicidae) in urban zones from Barcelona, Spain.

The feeding patterns of haematophagous arthropods are of major importance in the amplification and transmission of infectious disease agents to vertebrate hosts, including humans. The establishment of new vector populations in nonnative range might alter transmission networks. The Asian tiger mosquito Aedes albopictus (Skuse) represents an example of how an invasive species can alter the risk of viral transmission to humans. Blood meal molecular identification from two sympatric mosquito species (the invasive Ae. albopictus and the native Culex pipiens) was carried out by polymerase chain reaction-based methods. Samples were collected in Barcelona metropolitan area, Spain, from June to October 2009 as part of a monitoring-control program. Blood meals were identified to the species level in 30 Ae. albopictus and 43 Cx. pipiens. Ae. albopictus acquired blood exclusively from human hosts (100%), whereas Cx. pipiens fed on a diversity of avian and mammalian hosts, including 35.7% of blood meals from humans. Based on mosquito diet, our results suggest that the Ae. albopictus invasion in Spain might increase the risk of virus transmission to humans and could support local outbreaks of imported tropical viruses such as dengue and chikungunya. However, in the studied area, the presence of this invasive species would have a negligible effect on the transmission of zoonotic agents such as West Nile virus. However, Cx. pipiens could amplify and transmit West Nile virus, but avian contribution to its diet was lower than that reported in North America. Feeding patterns of these mosquito species may help to understand the flavivirus outbreaks recently reported in southwestern Europe.

Indigenous sex-selective salmon harvesting demonstrates pre-contact marine resource management in Burrard Inlet, British Columbia, Canada.

To gain insight into pre-contact Coast Salish fishing practices, we used new palaeogenetic analytical techniques to assign sex identifications to salmonid bones from four archaeological sites in Burrard Inlet (Tsleil-Waut), British Columbia, Canada, dating between about 2300-1000 BP (ca. 400 BCE-CE 1200). Our results indicate that male chum salmon (Oncorhynchus keta) were preferentially targeted at two of the four sampled archaeological sites. Because a single male salmon can mate with several females, selectively harvesting male salmon can increase a fishery's maximum sustainable harvest. We suggest such selective harvesting of visually distinctive male spawning chum salmon was a common practice, most effectively undertaken at wooden weirs spanning small salmon rivers and streams. We argue that this selective harvesting of males is indicative of an ancient and probably geographically widespread practice for ensuring sustainable salmon populations. The archaeological data presented here confirms earlier ethnographic accounts describing the selective harvest of male salmon.

On the relative roles of selection and genetic drift in shaping MHC variation.

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen-mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co-workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.

MHC diversity and differential exposure to pathogens in kestrels (Aves: Falconidae).

Pathogen diversity is thought to drive major histocompatibility complex (MHC) polymorphism given that host's immune repertories are dependent on antigen recognition capabilities. Here, we surveyed an extensive community of pathogens (n = 35 taxa) and MHC diversity in mainland versus island subspecies of the Eurasian kestrel Falco tinnunculus and in a sympatric mainland population of the phylogenetically related lesser kestrel Falco naumanni. Insular subspecies are commonly exposed to impoverished pathogen communities whilst different species' ecologies and contrasting life-history traits may lead to different levels of pathogen exposure. Although specific host traits may explain differential particular infections, overall pathogen diversity, richness and prevalence were higher in the truly cosmopolitan, euriphagous and long-distance disperser Eurasian kestrel than in the estenophagous, steppe-specialist, philopatric but long-distance migratory lesser kestrel. Accordingly, the continental population of Eurasian kestrels displayed a higher number (64 vs. 49) as well as more divergent alleles at both MHC class I and class II loci. Detailed analyses of amino acid diversity revealed that significant differences between both species were exclusive to those functionally important codons comprising the antigen binding sites. The lowest pathogen burdens and the smallest but still quite divergent set of MHC alleles (n = 16) were found in island Eurasian kestrels, where the rates of allele fixation at MHC loci seem to have occurred faster than at neutral markers. The results presented in this study would therefore support the role of pathogen diversity and abundance in shaping patterns of genetic variation at evolutionary relevant MHC genes.

Simultaneous analysis of multiple PCR amplicons enhances capillary SSCP discrimination of MHC alleles.

Major histocompatibility complex (MHC) genotyping still remains one of the most challenging issues for evolutionary ecologists. To date, none of the proposed methods have proven to be perfect, and all provide both important pros and cons. Although denaturing capillary electrophoresis has become a popular alternative, allele identification commonly relies upon conformational polymorphisms of two single-stranded DNA molecules at the most. Using the MHC class II (beta chain, exon 2) of the black kite (Aves: Accipitridae) as our model system, we show that the simultaneous analysis of overlapping PCR amplicons from the same target region substantially enhances allele discrimination. To cover this aim, we designed a multiplex PCR capable to generate four differentially sized and labeled amplicons from the same allele. Informative peaks to assist allele calling then fourfold those generated by the analysis of single PCR amplicons. Our approach proved successful to differentiate all the alleles (N=13) isolated from eight unrelated birds at a single optimal run temperature and electrophoretic conditions. In particular, we emphasize that this approach may constitute a straightforward and cost-effective alternative for the genotyping of single or duplicated MHC genes displaying low to moderate sets of divergent alleles.

Evaluating the quantity, quality and size distribution of cell-free DNA by multiplex droplet digital PCR.

Cell-free DNA (cfDNA) has become a comprehensive biomarker in the fields of non-invasive cancer detection and monitoring, organ transplantation, prenatal genetic testing and pathogen detection. While cfDNA samples can be obtained using a broad variety of approaches, there is an urgent need to standardize analytical tools aimed at assessing its basic properties. Typical methods to determine the yield and fragment size distribution of cfDNA samples are usually either blind to genomic DNA contamination or the presence of enzymatic inhibitors, which can confound and undermine downstream analyses. Here, we present a novel droplet digital PCR assay to identify suboptimal samples and aberrant cfDNA size distributions, the latter typically associated with high levels of circulating tumour DNA (ctDNA). Our assay was designed to promiscuously cross-amplify members of the human olfactory receptor (OR) gene family and includes a customizable diploid locus for the determination of absolute cfDNA concentrations. We demonstrate here the utility of our assay to estimate the yield and quality of cfDNA extracts and deduce fragment size distributions that correlate well with those inferred by capillary electrophoresis and high throughput sequencing. The assay described herein is a powerful tool to establish quality controls and stratify cfDNA samples based on presumed ctDNA levels, then facilitating the implementation of robust, cost-effective and standardized analytical workflows into clinical practice.

Genetic and evolutionary patterns of treatment resistance in relapsed B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) patients are typically treated with immunochemotherapy containing rituximab (rituximab, cyclophosphamide, hydroxydaunorubicin-vincristine (Oncovin), and prednisone [R-CHOP]); however, prognosis is extremely poor if R-CHOP fails. To identify genetic mechanisms contributing to primary or acquired R-CHOP resistance, we performed target-panel sequencing of 135 relapsed/refractory DLBCLs (rrDLBCLs), primarily comprising circulating tumor DNA from patients on clinical trials. Comparison with a metacohort of 1670 diagnostic DLBCLs identified 6 genes significantly enriched for mutations upon relapse. TP53 and KMT2D were mutated in the majority of rrDLBCLs, and these mutations remained clonally persistent throughout treatment in paired diagnostic-relapse samples, suggesting a role in primary treatment resistance. Nonsense and missense mutations affecting MS4A1, which encodes CD20, are exceedingly rare in diagnostic samples but show recurrent patterns of clonal expansion following rituximab-based therapy. MS4A1 missense mutations within the transmembrane domains lead to loss of CD20 in vitro, and patient tumors harboring these mutations lacked CD20 protein expression. In a time series from a patient treated with multiple rounds of therapy, tumor heterogeneity and minor MS4A1-harboring subclones contributed to rapid disease recurrence, with MS4A1 mutations as founding events for these subclones. TP53 and KMT2D mutation status, in combination with other prognostic factors, may be used to identify high-risk patients prior to R-CHOP for posttreatment monitoring. Using liquid biopsies, we show the potential to identify tumors with loss of CD20 surface expression stemming from MS4A1 mutations. Implementation of noninvasive assays to detect such features of acquired treatment resistance may allow timely transition to more effective treatment regimens.