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.

Abstract 534: An ultrasensitive high-plex assay detecting 24 PIK3CA mutations using SAGAsafe® technology and 6-color Crystal Digital PCR™

Starting from a minimally invasive liquid biopsy sample such as a blood sample, one can determine genetic and biological properties of cancer by measuring circulating tumor DNA (ctDNA). ctDNA measurements require a highly sensitive and reliable detection technology to quantify often low-level genetic aberrations within a high background of wild-type sequences.

The focus of this study was to demonstrate the performance of the patented SAGAsafe digital PCR chemistry using 6-color Crystal Digital PCR, by adapting a SAGAsafe assay for a clinically relevant oncogene, PIK3CA. PIK3CA encodes the p110-alpha kinase and is the target of a recently approved drug, alpelisib, for metastatic breast cancer. By combining the power of the SAGAsafe technology with that of 6-color Crystal Digital PCR on the naica system®, we present a highly multiplexed, ultrasensitive proof-of-concept PIK3CA detection assay capable of detecting and quantifying 24 PIK3CA mutations in a single digital PCR reaction.

By evaluating the sensitivity, specificity and precision of the high-plex 6-color PIK3CA digital PCR assay, we show that despite the complexity of the assay, the data were extremely clean resulting in highly accurate quantifications down to mutant allele frequencies (MAFs) in the 0.01% range. In addition, superb PIK3CA mutation detection specificity and sensitivity were observed with Limits of Blank ranging from 0.003% to <0.001% MAF, and Limits of Detection ranging from 0.009% to 0.002% MAF in a high background of wild-type DNA. Moreover, a high concordance was observed between Crystal Digital PCR PIK3CA results and NGS data on breast cancer biopsies.

By joining 6-color Crystal Digital PCR and ultrasensitive SAGAsafe technology, high sensitivity and specificity can be maintained even in a complex multiplex digital PCR assay. This proof-of-concept assay sets the stage for the future development of ultrasensitive, highly multiplexed cancer detection panels to achieve better patient stratification and monitoring during therapy.

Citation Format: Miguel Alcaide, Anthony M. George, Yilun Chen, Cécile Jovelet, Sahbi Ben Kilani, Lao H. Saal, Allison C. Mallory. An ultrasensitive high-plex assay detecting 24 PIK3CA mutations using SAGAsafe® technology and 6-color Crystal Digital PCR™ [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 534.

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.

Abstract PO-32: NFKBIZ 3′ UTR mutations confer selective growth advantage and affect drug response in diffuse large B-cell lymphoma

Introduction: The activated B cell-like (ABC) molecular subgroup of diffuse large B-cell lymphoma (DLBCL) is characterized by activation of NF-κB signaling. Recurrent mutations affecting genes such as MYD88, CD79A/B, and TNFAIP3 contribute to this in some cases, but there remain tumors with no known genetic basis for this pathway activation. This suggests that our understanding of ABC DLBCL drivers remains incomplete. Previously, NFKBIZ was shown to be amplified in 10% of ABC DLBCLs and to contribute to activation of NF-κB signaling. We recently described a novel pattern of mutations affecting the 3′ UTR of NFKBIZ resulting in an overall mutation rate of 30% (UTR or AMP) in ABC DLBCL. These NFKBIZ UTR mutations are mutually exclusive with MYD88 mutations, thus suggesting they may also lead to activation of NF-κB signaling. The NFKBIZ protein interacts with NF-κB transcription factors and is thought to regulate canonical NF-κB signaling. We hypothesized that NFKBIZ UTR mutations affect the normally rapid degradation of this mRNA by disrupting secondary structures recognized by RNA-binding proteins such as ribonucleases. The resulting elevated NFKBIZ mRNA levels would lead to accumulation of protein and may be a novel mechanism to promote cell growth and survival in ABC DLBCL.

Methods: NFKBIZ 3′ UTR mutations were introduced into a DLBCL cell line using the CRISPR-Cas9 system. A competitive growth assay with wild-type (WT) and CRISPR-mutant lines was performed to assess whether UTR mutations provide a growth advantage in culture (in vitro) and in mouse xenografts (in vivo). RNA-sequencing was then performed on WT and a subset of CRISPR-mutant lines and analyses were performed to identify genes upregulated by IκB-ζ in mutant lines. The IC50 of relevant drugs was determined by WST-1 assays after drug treatment on WT and mutant lines.

Results: Introduction of NFKBIZ mutations into a DLBCL cell line confirmed that UTR deletions lead to increased mRNA and protein levels. NFKBIZ UTR deletions give DLBCL cells a selective growth advantage over WT both in vitro and in vivo. RNA-sequencing of mutant and WT lines revealed possible transcriptional targets of NFKBIZ, including NF-κB targets and genes commonly overexpressed in ABC DLBCL. Novel candidate NFKBIZ targets were also discovered through this analysis, including CD274, the gene encoding PD-L1. Mutant cell lines had significantly higher IC50 compared to WT for the drugs ibrutinib, idelalisib, and masitinib, but not bortezomib, suggesting that NKFBIZ UTR mutations make cell lines more resistant to specific NF-κB pathway-targeted drugs.

Conclusions: This work directly establishes a role for NFKBIZ amplifications and 3′ UTR mutations in driving ABC DLBCL through NF-κB signaling. We demonstrate that these mutations can lead to overexpression of NFKBIZ and provide a selective growth advantage to cells both in vitro and in vivo. In addition, we found that these mutant lines were more resistant to some targeted lymphoma drugs but not others.

Citation Format: Sarah E. Arthur, Nicole Thomas, Jeffrey Tang, Christopher K. Rushton, Miguel Alcaide, Adèle Telenius, Shannon Healy, Anja Mottok, David W. Scott, Christian Steidl, Ryan D. Morin. NFKBIZ 3′ UTR mutations confer selective growth advantage and affect drug response in diffuse large B-cell lymphoma [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr PO-32.

Abstract IA42: Detecting and quantifying mutations associated with treatment resistance in aggressive lymphomas using ctDNA

A significant proportion of diffuse large B-cell lymphoma (DLBCL) patients treated with immunochemotherapy containing rituximab (R-CHOP) exhibit either primary or acquired treatment resistance. The advancement of therapeutics in the relapse setting has likely been encumbered by our limited understanding of the molecular features that underlie resistance to R-CHOP. Unfortunately, our knowledge of DLBCL genetics is mostly limited to analyses conducted on diagnostic tissue biopsies, which have not been exposed to the selective pressures imposed by therapy. Identifying genetic alterations that contribute to treatment resistance may reveal additional treatment options and lead to biomarkers allowing patients to be paired with appropriate treatments. Genetic subgroups are gaining popularity as a new strategy to implement precision medicine in DLBCL (1). The relevance of these and other biomarkers in the relapse setting remains unclear due to limited genetic exploration of relapsed and refractory DLBCL (rrDLBCL). Progress has been limited, in part, by the requirement of tissue biopsies collected after relapse. It is well established that quantitative genomic techniques such as digital PCR and targeted sequencing can be used to determine the proportion of tumor DNA in plasma from lymphoma patients (2). With a sufficiently broad panel, sequencing affords additional opportunities including the ability to identify subclonal structure and population dynamics over time. This presentation will discuss our recent analysis of a large collection of ctDNA primarily comprising DLBCL patients on various clinical trials (3). Targeted sequencing of these samples and comparison to exome data from a meta-cohort of previously characterized untreated DLBCL biopsies revealed six genes significantly enriched for mutations upon relapse. We found both TP53 and KMT2D were mutated in the majority of rrDLBCLs, and these mutations persisted in the dominant clone following relapse, suggesting a role in primary treatment resistance. By inferring subclonal dynamics, we observed recurrent patterns of clonal expansion and contraction following rituximab-based therapy, with MS4A1 mutations representing the only example of consistent clonal expansion. MS4A1 missense mutations within the transmembrane domains led to loss of CD20 expression in vitro, and patient tumors harboring these mutations lacked CD20 protein expression. Our analysis nominates TP53 and KMT2D mutation status as novel prognostic factors that may facilitate the identification of high-risk patients prior to therapy. Moreover, we have demonstrated 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. In certain scenarios whole-exome sequencing (WES) or whole-genome sequencing (WGS) can be successfully applied to ctDNA, thereby allowing the identification of mutations, structural variation, and copy number changes. Low-pass sequencing of shotgun libraries can also be used to ascertain course estimates of ctDNA levels as well as the copy number landscape (4). Given the importance of copy number and structural alterations in the inference of genetic subgroups, these methods may allow the exploration of these groups and their stability over time. Through a series of illustrative examples, this presentation will explore the benefits of each of these techniques in the study of tumor evolution and acquired treatment resistance in DLBCL.

References: 1. Morin RD, Scott DW. DLBCL subclassification: Divide and conquer? Blood 2020;135:1722–4. 2. Rossi D et al. The development of liquid biopsy for research and clinical practice in lymphomas: Report of the 15-ICML workshop on ctDNA. Hematol Oncol 2020;38:34–7. 3. Rushton CK et al. Genetic and evolutionary patterns of treatment resistance in relapsed B-cell lymphoma. Blood Adv 2020;4:2886–98. 4. Adalsteinsson VA et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat Commun 2017;8:1324.

Citation Format: Nicole Thomas, Laura K. Hilton, Neil Michaud, Kevin Bushell, Ryan Rys, Michael Jain, Lois Shepherd, Marco A. Marra, John Kuruvilla, Michael Crump, Koren Mann, Sarit Assouline, Christian Steidl, Mark S. Cragg, David W. Scott, Nathalie Johnson, Ryan D. Morin, Christopher K. Rushton, Sarah E. Arthur, Miguel Alcaide, Matthew Cheung, Aixiang Jiang, Krysta M. Coyle, Kirstie L. S. Cleary. Detecting and quantifying mutations associated with treatment resistance in aggressive lymphomas using ctDNA [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr IA42.

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.

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.

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 F_ST ), 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 F_ST , 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 F_ST 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.

Abstract B22: Capturing L1 retrotransposon-mediated DNA transductions in endometriosis associated ovarian cancers as a way to track tumor development

Endometrioid ovarian cancer (ENOC) and clear cell ovarian cancer (CCOC) share a common precursor lesion, endometriosis (ectopic growth of uterine lining), hence the designation endometriosis-associated ovarian cancer (EAOC). Women with endometriosis have up to three-fold increased risk of developing ENOC and CCOC. Efforts have been made to look for biomarkers that can help identifying women at risk of developing cancer; however, there are currently no biomarkers that stratify risk of cancer development. We performed whole-genome sequencing (WGS) on 29 ENOC and 35 CCOC cases and observed a frequent transduction event originating from an active LINE-1 (L1) retrotransposable element in the TTC28 gene. Such event occurred in 34% (10/29) of ENOC, and 31% (11/35) of CCOC cases. L1 retrotransposons are repetitive, mobile genetic elements capable of taking downstream DNA fragments and inserting them into random genomic locations via a process called 3’ transduction. Approximately 70-100 different potentially active L1s are epigenetically silenced in normal tissues, but tend to be reactivated in cancers. We subsequently used PCR to validate these TTC28-L1 transductions, and compared their presence to single nucleotide variations (SNVs) and frameshift mutations in formalin-fixed, paraffin-embedded (FFPE) tumor tissues from different tumor sites for 4 ENOC and 3 CCOC cases. We found that these transduction events along with classical driver mutations were almost ubiquitous across the tumor sites, suggesting these L1 events likely occurred early in the malignant transformation of EAOCs.

We developed a low-input, probe-based capture assay to test the presence of TTC28-L1 transductions as an alternative method to performing WGS. Oligonucleotide probes tiling 1 kb downstream of active L1s are used to capture DNA fragments containing the transduced DNA, and the fragments are sequenced on the MiSeq next-generation sequencing platform. Analyses are performed using the Geneious software and the published bioinformatics tool Socrates, specific for detecting DNA fragments with split reads (fragments with ends aligning to different parts of the genome). We successfully validated the assay on 9 cases with WGS data: 7 EAOC cases with TTC28-L1 transductions and 2 EAOC cases without TTC28-L1 transductions. DNA extracted from frozen tumor and buffy coat (normal control) were used for each case, and FFPE tissues were used for selected cases. All reads containing the transduction events aligned to genomic coordinates corresponding to the WGS data.

While L1-mediated DNA transductions are often passenger events during tumorigenesis, our data suggest that they likely occur early in ovarian cancer tumorigenesis. Our data show that this probe-based capture assay provides an alternative method to WGS, and may be useful in detecting active 3’ transductions in novel cases to track the development of ovarian tumors.

Citation Format: Zhouchunyang Xia, Dawn Cochrane, Michael Anglesio, Winnie Yang, Miguel Alcaide, Tayyebeh Nazeran, Janine Senz, Amy Lum, Ali Bashashati, Yikan Wang, Ryan Morin, Sohrab Shah, David Huntsman. Capturing L1 retrotransposon-mediated DNA transductions in endometriosis associated ovarian cancers as a way to track tumor development. [abstract]. In: Proceedings of the AACR Conference: Addressing Critical Questions in Ovarian Cancer Research and Treatment; Oct 1-4, 2017; Pittsburgh, PA. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(15_Suppl):Abstract nr B22.

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.

Single-agent panobinostat for relapsed/refractory diffuse large B-cell lymphoma: clinical outcome and correlation with genomic data. A phase 2 study of the Fondazione Italiana Linfomi

We investigated panobinostat 40 mg three times weekly in 35 adult patients with relapsed/refractory (R/R) diffuse large B-cell lymphoma (DLBCL). Overall response rate and complete response were 17.1% and 11.4%, respectively. Median progression-free survival (PFS) and overall survival were 2.4 and 7.6 months, respectively. Calculated 12, 24 and 36 months PFS were 26%, 11% and 11%, respectively. Four patients who achieved a sustained CR, continued receiving panobinostat for an overall period of 44, 48, 50, 62 months. Thrombocytopenia grade 3 (5 patients) and 4 (24 patients) represented the main toxic effect, causing dose reduction or treatment suspension in 19 patients. Genomic analysis was unable to identify any relationship between mutations and response; TP53 mutation appeared not to impact the clinical outcome. Overall, panobinostat has a modest activity in R/R DLBCL patients, however it can induce very long lasting responses in some cases. Thrombocytopenia frequently limits the use of this agent.

The genomic landscape of two Burkitt lymphoma cases and derived cell lines: comparison between primary and relapse samples

Relapse occurs in 10-40% of Burkitt lymphoma (BL) patients that have completed intensive chemotherapy regimens and is typically fatal. While treatment-naive BL has been characterized, the genomic landscape of BL at the time of relapse (rBL) has never been reported. Here, we present a genomic characterization of two rBL patients. The diagnostic samples had mutations common in BL, including MYC and CCND3. Additional mutations were detected at relapse, affecting important pathways such as NFκB (IKBKB) and MEK/ERK (NRAS) signaling, glutamine metabolism (SIRT4), and RNA processing (ZFP36L2). Genes implicated in drug resistance were also mutated at relapse (TP53, BAX, ALDH3A1, APAF1, FANCI). This concurrent genomic profiling of samples obtained at diagnosis and relapse has revealed mutations not previously reported in this disease. The patient-derived cell lines will be made available and, along with their detailed genetics, will be a valuable resource to examine the role of specific mutations in therapeutic resistance.

Targeted error-suppressed quantification of circulating tumor DNA using semi-degenerate barcoded adapters and biotinylated baits

Ultrasensitive methods for rare allele detection are critical to leverage the full potential offered by liquid biopsies. Here, we describe a novel molecular barcoding method for the precise detection and quantification of circulating tumor DNA (ctDNA). The major benefits of our design include straightforward and cost-effective production of barcoded adapters to tag individual DNA molecules before PCR and sequencing, and better control over cross-contamination between experiments. We validated our approach in a cohort of 24 patients with a broad spectrum of cancer diagnoses by targeting and quantifying single-nucleotide variants (SNVs), indels and genomic rearrangements in plasma samples. By using personalized panels targeting a priori known mutations, we demonstrate comprehensive error-suppression capabilities for SNVs and detection thresholds for ctDNA below 0.1%. We also show that our semi-degenerate barcoded adapters hold promise for noninvasive genotyping in the absence of tumor biopsies and monitoring of minimal residual disease in longitudinal plasma samples. The benefits demonstrated here include broad applicability, flexibility, affordability and reproducibility in the research and clinical settings.

Investigating the Genetic Causes of Sudden Unexpected Death in Children Through Targeted Next-Generation Sequencing Analysis

Background—

Inherited arrhythmia syndromes are responsible for a significant portion of autopsy-negative sudden unexpected death (SUD) cases, but molecular autopsy used to identify potentially causal variants is not routinely included in SUD investigations. We collaborated with a medical examiner's office to assist in finding a diagnosis for their autopsy-negative child SUD cases.

Methods and Results—

191 child SUD cases (<5 years of age) were selected for analyses. Our next generation sequencing panel incorporated 38 inherited arrhythmia syndrome candidate genes and another 33 genes not previously investigated for variants that may underlie SUDY pathophysiology. Overall, we identified 11 potentially causal disease-associated variants in 12 cases, for an overall yield of 6.3%. We also identified 31 variants of uncertain significance in 36 cases and 16 novel variants predicted to be pathogenic in silico in 15 cases. The disease-associated variants were reported to the medical examiner to notify surviving relatives and recommend clinical assessment.

Conclusions—

We have identified variants that may assist in the diagnosis of at least 6.3% of autopsy-negative child SUD cases and reduce risk of future SUD in surviving relatives. We recommend a cautious approach to variant interpretation. We also suggest inclusion of cardiomyopathy genes as well as other candidate SUD genes in molecular autopsy analyses.

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.

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.

Genetic Landscapes of Relapsed and Refractory Diffuse Large B-Cell Lymphomas

PURPOSE

Relapsed or refractory diffuse large B-cell lymphoma (rrDLBCL) is fatal in 90% of patients, and yet little is known about its biology.

EXPERIMENTAL DESIGN

Using exome sequencing, we characterized the mutation profiles of 38 rrDLBCL biopsies obtained at the time of progression after immunochemotherapy. To identify genes that may be associated with relapse, we compared the mutation frequency in samples obtained at relapse to an unrelated cohort of 138 diagnostic DLBCLs and separately amplified specific mutations in their matched diagnostic samples to identify clonal expansions.

RESULTS

On the basis of a higher frequency at relapse and evidence for clonal selection, TP53, FOXO1, MLL3 (KMT2C), CCND3, NFKBIZ, and STAT6 emerged as top candidate genes implicated in therapeutic resistance. We observed individual examples of clonal expansions affecting genes whose mutations had not been previously associated with DLBCL including two regulators of NF-κB: NFKBIE and NFKBIZ We detected mutations that may be affect sensitivity to novel therapeutics, such as MYD88 and CD79B mutations, in 31% and 23% of patients with activated B-cell-type of rrDLBCL, respectively. We also identified recurrent STAT6 mutations affecting D419 in 36% of patients with the germinal center B (GCB) cell rrDLBCL. These were associated with activated JAK/STAT signaling, increased phospho-STAT6 protein expression and increased expression of STAT6 target genes.

CONCLUSIONS

This work improves our understanding of therapeutic resistance in rrDLBCL and has identified novel therapeutic opportunities especially for the high-risk patients with GCB-type rrDLBCL. Clin Cancer Res; 22(9); 2290-300. ©2015 AACR.

Determinants and short-term physiological consequences of PHA immune response in lesser kestrel nestlings

Individual immune responses are likely affected by genetic, physiological, and environmental determinants. We studied the determinants and short-term consequences of Phytohaemagglutinin (PHA) induced immune response, a commonly used immune challenge eliciting both innate and acquired immunity, on lesser kestrel (Falco naumanni) nestlings in semi-captivity conditions and with a homogeneous diet composition. We conducted a repeated measures analyses of a set of blood parameters (carotenoids, triglycerides, β-hydroxybutyrate, cholesterol, uric acid, urea, total proteins, and total antioxidant capacity), metabolic (resting metabolic rate), genotypic (MHC class II B heterozygosity), and biometric (body mass) variables. PHA challenge did not affect the studied physiological parameters on a short-term basis (<12 hr), except plasma concentrations of triglycerides and carotenoids, which decreased and increased, respectively. Uric acid was the only physiological parameter correlated with the PHA induced immune response (skin swelling), but the change of body mass, cholesterol, total antioxidant capacity, and triglycerides between sessions (i.e., post-pre treatment) were also positively correlated to PHA response. No relationships were detected between MHC gene heterozygosity or resting metabolic rate and PHA response. Our results indicate that PHA response in lesser kestrel nestlings growing in optimal conditions does not imply a severe energetic cost 12 hr after challenge, but is condition-dependent as a rapid mobilization of carotenoids and decrease of triglycerides is elicited on a short-term basis.

Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next-generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide-binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high-frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non-passerine bird.

Major histocompatibility complex class I evolution in songbirds: universal primers, rapid evolution and base compositional shifts in exon 3

Genes of the Major Histocompatibility Complex (MHC) have become an important marker for the investigation of adaptive genetic variation in vertebrates because of their critical role in pathogen resistance. However, despite significant advances in the last few years the characterization of MHC variation in non-model species still remains a challenging task due to the redundancy and high variation of this gene complex. Here we report the utility of a single pair of primers for the cross-amplification of the third exon of MHC class I genes, which encodes the more polymorphic half of the peptide-binding region (PBR), in oscine passerines (songbirds; Aves: Passeriformes), a group especially challenging for MHC characterization due to the presence of large and complex MHC multigene families. In our survey, although the primers failed to amplify exon 3 from two suboscine passerine birds, they amplified exon 3 of multiple MHC class I genes in all 16 species of oscine songbirds tested, yielding a total of 120 sequences. The 16 songbird species belong to 14 different families, primarily within the Passerida, but also in the Corvida. Using a conservative approach based on the analysis of cloned amplicons (n = 16) from each species, we found between 3 and 10 MHC sequences per individual. Each allele repertoire was highly divergent, with the overall number of polymorphic sites per species ranging from 33 to 108 (out of 264 sites) and the average number of nucleotide differences between alleles ranging from 14.67 to 43.67. Our survey in songbirds allowed us to compare macroevolutionary dynamics of exon 3 between songbirds and non-passerine birds. We found compelling evidence of positive selection acting specifically upon peptide-binding codons across birds, and we estimate the strength of diversifying selection in songbirds to be about twice that in non-passerines. Analysis using comparative methods suggest weaker evidence for a higher GC content in the 3rd codon position of exon 3 in non-passerine birds, a pattern that contrasts with among-clade GC patterns found in other avian studies and may suggests different mutational mechanisms. Our primers represent a useful tool for the characterization of functional and evolutionarily relevant MHC variation across the hyperdiverse songbirds.

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.

Feeding patterns of potential West Nile virus vectors in south-west Spain

Background

Mosquito feeding behaviour determines the degree of vector–host contact and may have a serious impact on the risk of West Nile virus (WNV) epidemics. Feeding behaviour also interacts with other biotic and abiotic factors that affect virus amplification and transmission.

Methodology/Principal Findings

We identified the origin of blood meals in five mosquito species from three different wetlands in SW Spain. All mosquito species analysed fed with different frequencies on birds, mammals and reptiles. Both ‘mosquito species’ and ‘locality’ explained a similar amount of variance in the occurrence of avian blood meals. However, ‘season of year’ was the main factor explaining the presence of human blood meals. The differences in diet resulted in a marked spatial heterogeneity in the estimated WNV transmission risk. Culex perexiguus, Cx. modestus and Cx. pipiens were the main mosquito species involved in WNV enzootic circulation since they feed mainly on birds, were abundant in a number of localities and had high vector competence. Cx. perexiguus may also be important for WNV transmission to horses, as are Cx. pipiens and Cx. theileri in transmission to humans. Estimates of the WNV transmission risk based on mosquito diet, abundance and vector competence matched the results of previous WNV monitoring programs in the area. Our sensitivity analyses suggested that mosquito diet, followed by mosquito abundance and vector competence, are all relevant factors in understanding virus amplification and transmission risk in the studied wild ecosystems. At some of the studied localities, the risk of enzootic circulation of WNV was relatively high, even if the risk of transmission to humans and horses was less.

Conclusions/Significance

Our results describe for first time the role of five WNV candidate vectors in SW Spain. Interspecific and local differences in mosquito diet composition has an important effect on the potential transmission risk of WNV to birds, horses and humans.

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.

Sampling strategies for accurate computational inferences of gametic phase across highly polymorphic major histocompatibility complex loci

Background

Genes of the Major Histocompatibility Complex (MHC) are very popular genetic markers among evolutionary biologists because of their potential role in pathogen confrontation and sexual selection. However, MHC genotyping still remains challenging and time-consuming in spite of substantial methodological advances. Although computational haplotype inference has brought into focus interesting alternatives, high heterozygosity, extensive genetic variation and population admixture are known to cause inaccuracies. We have investigated the role of sample size, genetic polymorphism and genetic structuring on the performance of the popular Bayesian PHASE algorithm. To cover this aim, we took advantage of a large database of known genotypes (using traditional laboratory-based techniques) at single MHC class I (N = 56 individuals and 50 alleles) and MHC class II B (N = 103 individuals and 62 alleles) loci in the lesser kestrel Falco naumanni.

Findings

Analyses carried out over real MHC genotypes showed that the accuracy of gametic phase reconstruction improved with sample size as a result of the reduction in the allele to individual ratio. We then simulated different data sets introducing variations in this parameter to define an optimal ratio.

Conclusions

Our results demonstrate a critical influence of the allele to individual ratio on PHASE performance. We found that a minimum allele to individual ratio (1:2) yielded 100% accuracy for both MHC loci. Sampling effort is therefore a crucial step to obtain reliable MHC haplotype reconstructions and must be accomplished accordingly to the degree of MHC polymorphism. We expect our findings provide a foothold into the design of straightforward and cost-effective genotyping strategies of those MHC loci from which locus-specific primers are available.

Rectal infections with chlamydia and gonorrhoea in women attending a multiethnic sexually transmitted diseases urban clinic

Sexually transmitted rectal infections with Chlamydia trachomatis (CT) and Neisseria gonorrhoeae (GC) have been well documented in men who have sex with men (MSM). Few studies have described infections in women who engage in anal intercourse. We performed testing for rectal infections in women who reported ano-receptive intercourse at the Miami Dade Health Department Sexually Transmitted Diseases (STD) clinic and report the prevalence and characteristics of women with rectal CT or GC infections. Our results revealed a prevalence of 17.5% for rectal chlamydia and 13.4% for rectal gonorrhoea. Urine-based screening alone would have missed 6% of rectal chlamydia infections and 35% of rectal gonorrhoea infections. Anal symptoms were reported in 12.5% of women with rectal chlamydia infections. The only associated factor identified was an age less than 28 years. We conclude that rectal screening for CT and GC should be included in STD prevention strategies, especially in the younger population.

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.

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.

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca

Background

Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts.

Findings

Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (β chain). Taking the pied flycatcher Ficedula hypoleuca as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species.

Conclusions

Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.

Interaction of an ordered mesoporous bioactive glass with osteoblasts, fibroblasts and lymphocytes, demonstrating its biocompatibility as a potential bone graft material

Ordered mesoporous 85SiO(2)-10CaO-5P(2)O(5) bioactive glass (MBG85) is an excellent candidate as a graft for bone tissue regeneration, owing to its excellent textured properties, structural characteristics and crystalline apatite rate formation. To assess MBG85 biocompatibility, different parameters have been evaluated (cell morphology, size/complexity, proliferation, viability, cell cycle, reactive oxygen species content, lactate dehydrogenase release) using human Saos-2 osteoblasts after treatment with either MBG85 extracts or 1% MBG85 directly added to cells. The osteoblast response to MBG85 was compared with L929 fibroblast behaviour after the same treatment. The high cell viability observed and the absence of signs of cell damage in both cell types demonstrates MBG85 biocompatibility. Only a cytostatic effect was observed through the reduction of cell proliferation, related with the initial Ca elution, whereas Si leaching did not result into any negative effect. In vitro lymphocytic proliferation analysis was also carried out with SR.D10 clone after treatment with either MBG85 extracts or culture supernatants of L929 fibroblasts previously treated with 1% MBG85 (cell-conditioned extracts). The absence of modification of in vitro T-cell response underlines the biocompatibility of MBG85 and its potential application in the field of bone and dental grafting.

Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals

Emerging infectious diseases represent a challenge for global economies and public health. About one fourth of the last pandemics have been originated by the spread of vector-borne pathogens. In this sense, the advent of modern molecular techniques has enhanced our capabilities to understand vector-host interactions and disease ecology. However, host identification protocols have poorly profited of international DNA barcoding initiatives and/or have focused exclusively on a limited array of vector species. Therefore, ascertaining the potential afforded by DNA barcoding tools in other vector-host systems of human and veterinary importance would represent a major advance in tracking pathogen life cycles and hosts. Here, we show the applicability of a novel and efficient molecular method for the identification of the vertebrate host's DNA contained in the midgut of blood-feeding arthropods. To this end, we designed a eukaryote-universal forward primer and a vertebrate-specific reverse primer to selectively amplify 758 base pairs (bp) of the vertebrate mitochondrial Cytochrome c Oxidase Subunit I (COI) gene. Our method was validated using both extensive sequence surveys from the public domain and Polymerase Chain Reaction (PCR) experiments carried out over specimens from different Classes of vertebrates (Mammalia, Aves, Reptilia and Amphibia) and invertebrate ectoparasites (Arachnida and Insecta). The analysis of mosquito, culicoid, phlebotomie, sucking bugs, and tick bloodmeals revealed up to 40 vertebrate hosts, including 23 avian, 16 mammalian and one reptilian species. Importantly, the inspection and analysis of direct sequencing electropherograms also assisted the resolving of mixed bloodmeals. We therefore provide a universal and high-throughput diagnostic tool for the study of the ecology of haematophagous invertebrates in relation to their vertebrate hosts. Such information is crucial to support the efficient management of initiatives aimed at reducing epidemiologic risks of arthropod vector-borne pathogens, a priority for public health.

Strong philopatry derived from capture-recapture records does not lead to fine-scale genetic differentiation in lesser kestrels

1. The integration of capture-recapture and molecular approaches can improve our understanding of the consequences of habitat fragmentation on population connectivity. Here we employed microsatellites to test dispersal hypotheses derived from intense and long-term ringing programmes of the lesser kestrel Falco naumanni in Western Europe. 2. Re-encounters of 1308 marked individuals in Spain have revealed that most first-time breeders settled within 10 km from their natal colony, with a negative association between dispersal and geographical distance. Although these findings would predict fine-scale spatial patterns of genetic differentiation, the genetic impact of rarely reported events concerning long-distance effective dispersal (> 100 km) is unknown. 3. First, we investigated a spatially structured and geographically isolated population located in north-eastern Spain, where capture-recapture records and genetic data could be appropriately compared over similar spatial and temporal scales. Spatial autocorrelation analyses (N = 174 nestlings from different broods) did not reveal either significant differences in average relatedness at any distance class nor decreased relatedness as a function of distance. At a broader spatial scale, Bayesian analysis of population structure (N = 432 nestlings) indicated panmixia across Western Europe. However, F(ST) comparisons between four geographically distinct populations indicated low but significant genetic differentiation. 4. Our genetic data would therefore challenge traditional assumptions associating philopatry with the emergence of fine-scale genetic structuring. This could be because even low levels of gene flow are enough to preclude the development of local genetic structure. Nevertheless, the analysis of a geographically isolated and small population from Southern France exemplifies a situation in which restricted dispersal has translated into weak but consistently significant genetic differentiation. 5. Relevant to conservation genetics and evolutionary biology, our results may lessen the genetic concerns derived from population fragmentation at relatively small geographical scales in species with apparently limited dispersal abilities, but raises concerns about increased genetic divergence in small and isolated demes.

Population fragmentation leads to isolation by distance but not genetic impoverishment in the philopatric Lesser Kestrel: a comparison with the widespread and sympatric Eurasian Kestrel

Population fragmentation is a widespread phenomenon usually associated with human activity. As a result of habitat transformation, the philopatric and steppe-specialist Lesser Kestrel Falco naumanni underwent a severe population decline during the last century that increased population fragmentation throughout its breeding range. In contrast, the ubiquitous Eurasian Kestrel Falco tinnunculus did not suffer such adverse effects, its breeding range still remaining rather continuous. Using microsatellites, we tested the effects of population fragmentation on large-scale spatial patterns of genetic differentiation and diversity by comparing these two sympatric and phylogenetically related species. Our results suggest that habitat fragmentation has increased genetic differentiation between Lesser Kestrel populations, following an isolation-by-distance pattern, while the population of Eurasian Kestrels is panmictic. Contrary to expectations, we did not detect significant evidence of reduced genetic variation or increased inbreeding in Lesser Kestrels. Although this study reports genetic differentiation in a species that has potential for long-distance dispersal but philopatry-limited gene flow, large enough effective population sizes and migration may have been sufficient to mitigate genetic depauperation. A serious reduction of genetic diversity in Lesser Kestrels would, therefore, only be expected after severe population bottlenecks following extreme geographic isolation.

Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.

Concurrent infection with Trichinella spiralis and other helminths in pigs

The aim of this study was to investigate possible influence of different helmintosis in the development of Trichinella spiralis in experimental infected pigs. Forty-two Iberian pigs were allocated to six groups. Three groups were single inoculated with Ascaris suum, Metastrongylus apri or T. spiralis, respectively. Two groups were co-infected with T. spiralis and A. suum or T. spiralis and M. apri, respectively, while the last group included uninfected control pigs. Clinical signs were only observed in pigs with single or concurrent M. apri infections, with more severe respiratory symptoms in pigs with mixed M. apri infection. The number of A. suum and M. apri lung larvae, intestinal larvae of A. suum and adult M. apri were reduced in pigs with mixed Trichinella infections compared to pigs with single infections. In contrast, the number of liver white spots was higher in pigs with mixed infections. While T. spiralis muscular larval burdens were increased in pigs concomitantly infected with M. apri, they were reduced in pigs concomitantly infected with A. suum, compared to pigs receiving single infections with either of these helminths. Pigs with single or mixed A. suum infections showed higher eosinophil levels compared to the remaining groups. IgGt, IgG1, IgG2 and IgM against T. spiralis antigen could not be detected in pigs with single Ascaris or Metastrongylus infections, indicating that no cross-antibodies were produced. IgGt, IgG1 and IgM antibodies were detected earlier and generally at higher levels in mixed T. spiralis infections compared to single T. spiralis infections. The results suggest that T. spiralis had a low synergistic interaction with M. apri in concomitantly infected pigs, and an antagonistic interaction in concurrent infection with A. suum.

Epidemiology of Oestrus ovis (Linneo, 1761) infestation in goats in Spain

This survey was conducted to determine the chronobiology and seroprevalence of nasal bot infestation (Oestrus ovis) in Spain and to identify the risk factors associated with this disease in caprine herds. A total of 1590 sera from adult goats were collected at random on 175 farms in southwestern Spain. Sera were tested by ELISA, using crude protein from second stage larvae as antigen. The mean seroprevalence was 46.04% and mean percentage of optical densities was 41.83. These data indicate a high prevalence of this parasite in the investigated areas. The serological survey revealed that goats managed at higher altitudes, at meridians latitudes and on farms with small herds had a smaller probability of infestation. Eighty goat heads, obtained from abattoirs in the central region of Spain, were collected and examined for nasal botflies from February to October 2002. O. ovis larval stages were recovered from the nasal-sinus cavities of 23 goats, reaching a prevalence of 34.94%. The mean larval burden was 3.9 larvae per infested head. No first instars were found during February and March, when the second instar reached its larger count. The third instar was observed in very small number during the whole period of study, with one peak occurring in July--August. These data show the existence of a favourable period for the development of larval instars of O. ovis in goats that starts in February and finishes in September.

Analysis of larval antigens of Oestrus ovis for the diagnosis of oestrosis by enzyme-linked immunosorbent assay

A serodiagnostic test for the diagnosis of infestation by the sheep nasal bot fly, Oestrus ovis (Linné) was examined. The enzyme-linked immunosorbent assay (ELISA) technique was used to analyze and compare the production of immunoglobulin G (IgG) antibodies against excretory-secretory products (ESP) and crude extract (CE) antigens from all the different larval stages of O. ovis in the sera of 276 adult sheep sampled in summer (n = 135) and winter (n = 141). ESP from first stage larvae was the most sensitive, coating antigen in winter and ESP from second stage larvae during summer. The most specific values were obtained by ESP against L1 in winter and by CE against L3 in summer. These results show that the stage of larval development has a significant impact on the humoral immune response over the course of a season. A significant correlation (P < 0.001) was found between the number of O. ovis larvae and the serum antibody levels using all differents antigens, except L3 CE. In Spain, where a long favourable period exists for the evolution and development of the different stage larvae between March and November, the ELISA test using L1 ESP antigen during winter and L2 ESP antigen in summer may be used for ovine oestrosis immunodiagnosis.