Search

Abstract Genes of the major histocompatibility complex (MHC) are crucial for adaptive immunity in jawed vertebrates, and theory predicts that there should be mate choice for optimizing MHC constitution in the offspring. In a previous study, we demonstrated a non-random female choice of extra-pair males in the bluethroat (Luscinia svecica), yielding offspring that was closer to an intermediate MHC class II (MHCII) allele count than their within-pair halfsiblings. The present study tests whether social pairs with only within-pair young (WPY) in their brood, in the same study population, had a combined MHC-constitution closer to a presumed intermediate optimum, than social pairs with extra-pair young (EPY), with a corresponding pattern in their offspring. As expected, we found that WPY from pure WPY-broods were more MHC-optimal than WPY from mixed broods, but only in broods of young (second year) males. Correspondingly, there was a tendency for social pairs with only WPY in their brood to be more MHC-compatible than social pairs with EPY in their brood, when the male was young. Older bluethroat males have considerably larger testes than young males, and their higher sperm competitiveness could help them secure paternity in their own brood, also when they are not MHC-compatible. In other words, in the sexual conflict over paternity, females may be more likely to realise their preference for a MHC-compatible mate when paired to a young male. As a possible fitness indicator, immune responsiveness to an injected antigen (PHA) was elevated for offspring closer to “the golden mean” in MHCII allele count. Significance statement This study contributes to our understanding of MHC-based mate choice in extra-pair mating systems, by showing that female bluethroats (Luscinia svecica) with an MHCII-compatible social mate tend to have no extra-pair young in their brood, but only when the social male is young. This elucidates a possible sexual conflict, in which older social males are able to override female preferences and prevent other males from gaining paternity in their brood through higher sperm production. Studying systems in which extra-pair paternity occurs offers an insight into the genetic benefits of mate choice, as extra-pair males, in contrast to social males, generally contribute only sperm. Further, the strict and thorough genotyping scheme applied in this study enabled us to demonstrate a preference for “the golden mean” in MHC-diversity in a species with one of the highest MHC class II-diversity known to date.

Genes of the major histocompatibility complex (MHC) are crucial for adaptive immunity, and theory predicts that the optimal number of alleles is intermediate rather than maximal. Support for such an intermediate optimum in MHC allele count has mainly been demonstrated in studies on parasite load and reproductive traits. In this study we tested if an intermediate number of MHC class II (MHCII) alleles based on positively selected sites was correlated with an increased skin-swelling response to two consecutive injections of phytohaemagglutinin (PHA) in bluethroat Luscinia svecica nestlings. We found that nestlings that were closer to the population mean in their number of alleles had a higher response to the first injection, supporting an intermediate, optimal number of alleles. After the second PHA injection, we found that the skin swelling was associated with MHCII allele count only at the among-brood level, and with paternity (extra-pair offspring had a larger swelling than within-pair offspring). We therefore conclude that the PHA response is partly correlated with MHCII allele count, but that additional factors must be involved in both the PHA response and in female extra-pair mate choice in this species.

Growing fishing pressure worldwide has led to an increase in the exploitation of cephalopod products in fish markets where these taxa were not traditionally consumed. Squid catches have surged to meet that demand, yet for many species little is known about their role in food webs. Doryteuthis sanpaulensis is an important squid species in southeastern Brazilian fisheries. Despite many previous efforts at morphological analysis of its diet, few demersal and benthic species of fishes, crustaceans and mollusks have been identified to species level because the food is consumed in small digestible fragments. Here we used metabarcoding to analyze the diet of adult D. sanpaulensis caught as bycatch in the southern Brazilian sardine fishery. MOTUs generated from COI amplicons were assigned to taxa by matching against the NCBI nt database. Fishes constitute the majority of the diet of the analyzed samples. Considerable variability in the relative read abundance and frequency of occurrence of prey items across samples indicates the importance of increased sample sizes in analyses investigating ontogenetic, spatial, or temporal variation in diet. The results elucidate the rich diet of D. sanpaulensis off the Brazilian coast, and specifically that its varied diet includes more neritic diversity than previous studies have indicated.

Genotyping of classical major histocompatibility complex (MHC) genes is challenging when they are hypervariable and occur in multiple copies. In this study, we used several different approaches to genotype the moderately variable MHC class I exon 3 (MHCIe3) and the highly polymorphic MHC class II exon 2 (MHCIIβe2) in the bluethroat (Luscinia svecica). Two family groups (eight individuals) were sequenced in replicates at both markers using Ion Torrent technology with both a single‐ and a dual‐indexed primer structure. Additionally, MHCIIβe2 was sequenced on Illumina MiSeq. Allele calling was conducted by modifications of the pipeline developed by Sommer et al. (BMC Genomics, 14, 2013, 542) and the software AmpliSAS. While the different genotyping strategies gave largely consistent results for MHCIe3, with a maximum of eight alleles per individual, MHCIIβe2 was remarkably complex with a maximum of 56 MHCIIβe2 alleles called for one individual. Each genotyping strategy detected on average 50%–82% of all MHCIIβe2 alleles per individual, but dropouts were largely allele‐specific and consistent within families for each strategy. The discrepancies among approaches indicate PCR biases caused by the platform‐specific primer tails. Further, AmpliSAS called fewer alleles than the modified Sommer pipeline. Our results demonstrate that allelic dropout is a significant problem when genotyping the hypervariable MHCIIβe2. As these genotyping errors are largely nonrandom and method‐specific, we caution against comparing genotypes across different genotyping strategies. Nevertheless, we conclude that high‐throughput approaches provide a major advance in the challenging task of genotyping hypervariable MHC loci, even though they may not reveal the complete allelic repertoire.

Background With 25% prevalence of Crohn's disease, Familial GUCY2C diarrhea syndrome (FGDS) is a monogenic disorder potentially suited to study initiating factors in inflammatory bowel disease (IBD). We aimed to characterize the impact of an activating GUCY2C mutation on the gut microbiota in patients with FGDS controlling for Crohn's disease status and to determine whether changes share features with those observed in unrelated patients with IBD. Methods Bacterial DNA from fecal samples collected from patients with FGDS (N = 20), healthy relatives (N = 11), unrelated healthy individuals (N = 263), and IBD controls (N = 46) was subjected to sequencing of the V3-V4 region of the 16S rRNA gene to determine gut microbiota composition. Food frequency questionnaires were obtained from patients with FGDS and their relatives. Results Compared with healthy controls, FGDS displayed prominent changes in many microbial lineages including increase in Enterobacteriaceae, loss of Bifidobacterium and Faecalibacterium prausnitzii but an unchanged intraindividual (alpha) diversity. The depletion of F. prausnitzii is in line with what is typically observed in Crohn's disease. There was no significant difference in the dietary profile between the patients and related controls. The gut microbiota in related and unrelated healthy controls was also similar, suggesting that diet and familial factors do not explain the gut microbiota alterations in FGDS. Conclusions The findings support that the activating mutation in GUCY2C creates an intestinal environment with a major influence on the microbiota, which could contribute to the increased susceptibility to IBD in patients with FGDS.

Natural history museum collections represent a vast source of ancient and historical DNA samples from extinct taxa that can be utilized by high-throughput sequencing tools to reveal novel genetic and phylogenetic information about them. Here, we report on the successful sequencing of complete mitochondrial genome sequences (mitogenomes) from eleven extinct bird species, using de novo assembly of short sequences derived from toepad samples of degraded DNA from museum specimens. For two species (the Passenger Pigeon Ectopistes migratorius and the South Island Piopio Turnagra capensis), whole mitogenomes were already available from recent studies, whereas for five others (the Great Auk Pinguinis impennis, the Imperial Woodpecker Campehilus imperialis, the Huia Heteralocha acutirostris, the Kauai Oo Moho braccathus and the South Island Kokako Callaeas cinereus), there were partial mitochondrial sequences available for comparison. For all seven species, we found sequence similarities of >98%. For the remaining four species (the Kamao Myadestes myadestinus, the Paradise Parrot Psephotellus pulcherrimus, the Ou Psittirostra psittacea and the Lesser Akialoa Akialoa obscura), there was no sequence information available for comparison, so we conducted blast searches and phylogenetic analyses to determine their phylogenetic positions and identify their closest extant relatives. These mitogenomes will be valuable for future analyses of avian phylogenetics and illustrate the importance of museum collections as repositories for genomics resources. The final version of this research has been published in Molecular Ecology Resources. © 2017 Wiley

Objective Gut microbiota could influence gut, as well as hepatic and biliary immune responses. We therefore thoroughly characterised the gut microbiota in primary sclerosing cholangitis (PSC) compared with healthy controls (HC) and patients with ulcerative colitis without liver disease. Design We prospectively collected 543 stool samples. After a stringent exclusion process, bacterial DNA was submitted for 16S rRNA gene sequencing. PSC and HC were randomised to an exploration panel or a validation panel, and only significant results (p

Natural killer T (NKT) cells are activated by lipid antigens presented by CD1d molecules and represent a major lymphocyte subset of the liver. NOD.c3c4 mice spontaneously develop biliary inflammation in extra- and intrahepatic bile ducts. We demonstrated by flow cytometry that invariant NKT (iNKT) cells were more abundant in the thymus, spleen, and liver of NOD.c3c4 mice compared to NOD mice. iNKT cells in NOD.c3c4 mice displayed an activated phenotype. Further, NOD and NOD.Cd1d-/- mice were irradiated and injected with NOD.c3c4 bone marrow, and injection of NOD.c3c4 bone marrow resulted in biliary infiltrates independently of CD1d expression in recipient mice. Activation or blocking of NKT cells with α-galactosylceramide or anti-CD1d antibody injections did not affect the biliary phenotype of NOD.c3c4 mice. NOD.c3c4.Cd1d-/- mice were generated by crossing NOD.Cd1d-/- mice onto a NOD.c3c4 background. NOD.c3c4.Cd1d-/- and NOD.c3c4 mice developed the same extent of biliary disease. This study demonstrates that iNKT cells are more abundant and activated in the NOD.c3c4 model. The portal inflammation of NOD.c3c4 mice can be transferred to irradiated recipients, which suggests an immune-driven disease. Our findings imply that NKT cells can potentially participate in the biliary inflammation, but are not the primary drivers of disease in NOD.c3c4 mice.

Previous studies based on single mitochondrial markers have shown that the common raven (Corvus corax) consists of two highly diverged lineages that are hypothesised to have undergone speciation reversal upon secondary contact. Furthermore, common ravens are paraphyletic with respect to the Chihuahuan raven (C. cryptoleucus) based on mitochondrial DNA (mtDNA). Here we explore the causes of mtDNA paraphyly by sequencing whole mitochondrial genomes of 12 common ravens from across the Northern Hemisphere, in addition to three Chihuahuan ravens and one closely related brown-necked raven (C. ruficollis) using a long-range PCR protocol. Our raven mitogenomes ranged between 16925-16928 bp in length. GC content varied from 43.3% to 43.8% and the 13 protein coding genes, two rRNAs and 22 tRNAs followed a standard avian mitochondrial arrangement. The overall divergence between the two common raven clades was 3% (range 0.3-5.8% in 16 regions including the protein coding genes, rRNAs and the control region). Phylogenies constructed from whole mitogenomes recovered the previously found mitochondrial sister relationship between the common raven California clade and the Chihuahuan raven (overall divergence 1.1%), which strengthens the hypothesis that mtDNA paraphyly in the common raven results from speciation reversal of previously distinct Holarctic and California lineages.

Background One of the biggest challenges in avian taxonomy is the delimitation of allopatric species because their reproductive incompatibility cannot be directly studied in the wild. Instead, reproductive incompatibility has to be inferred from multiple, divergent character sets that indicate a low likelihood of allopatric populations amalgamating upon secondary contact. A set of quantitative criteria for species delimitation has been developed for avian taxonomy. Results Here, we report a broad multi-trait comparison of the two insular subspecies of the Blue Chaffinch Fringilla teydea, endemic to the pine forests of Tenerife (ssp. teydea) and Gran Canaria (ssp. polatzeki) in the Canary Islands. We found that the two taxa were reciprocally monophyletic in their whole mitogenomes and two Z chromosome introns. The genetic distance in mitogenomes indicates around 1 Mya of allopatric evolution. There were diagnostic differences in body morphometrics, song and plumage reflectance spectra, whose combined divergence score (=11) exceeds the threshold level (=7) set for species delimitation by Tobias et al. (Ibis 152:724–746, 2010). Moreover, we found a marked divergence in sperm lengths with little range overlap. Relatively long sperm with low intra- and intermale CV compared to other passerines suggest a mating system with high levels of sperm competition (extrapair paternity) in these taxa. Conclusion The large and diagnostic divergences in multiple functional traits qualify for species rank, i.e., Tenerife Blue Chaffinch (Fringilla teydea) and Gran Canaria Blue Chaffinch (Fringilla polatzeki). We encourage a wider use of sperm traits in avian taxonomy because sperm divergences might signal reproductive incompatibility at the postcopulatory prezygotic stage, especially in species with sperm competition.

Background Post-operative liver regeneration may contribute to tumor recurrence. There is a theoretical need for an adjuvant therapy that can suppress tumor growth without adversely affecting post-operative liver regeneration. Objective To evaluate the effect of RAF inhibitor Sorafenib on cell viability and proliferation of hepatoma cells and hepatocytes in vitro and in an in vivo rat model. Methods Cell viability, DNA synthesis, and RAF/MAPK kinase activity in the primary hepatocyte and hepatoma cell lines were investigated after Sorafenib exposure. Sequence analysis of the B-RAF gene in hepatic cells was determined. Tumor markers were compared within the rats after 70% hepatectomy with or without daily oral gavages of Sorafenib. Liver regeneration was assessed by liver function tests and proliferation markers. Results Primary hepatocytes showed higher cell viability, proliferation rate, and stronger RAF/MAPK kinase activity compared with hepatoma cell lines. The in vivo tumor volumes, size, and metastases were significantly decreased (P 0.05). B-RAF V600E mutation was not detected neither in the hepatic cells nor untransformed hepatocytes. Conclusions The RAF targeted inhibitor can reduce tumor growth without retarding liver regeneration in this experiment. J. Surg. Oncol. 2013;107:393–401. © 2012 Wiley Periodicals, Inc.

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide-binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.

Allelic variants of fucosyltransferases 2 and 3 (FUT2/3) influence serum levels of CA19-9, a screening parameter commonly used for detection of biliary malignancy in PSC. We aimed at improving diagnostic accuracy of CA19-9 by determining the impact of FUT2/3 genotypes.CA19-9 levels were measured in 433 PSC patients, 41 of whom had biliary malignancy. Genotypes for FUT3 and FUT2 were used to assign patients to one of three groups: A, no FUT3 activity regardless of FUT2 activity; B, both FUT2 and FUT3 activity and C, no FUT2 activity without loss of FUT3 activity. Group-specific cut-off values were determined by Youden's index.The median CA19-9 values of cancer-free patients were significantly different (p0.001) in Groups A (2.0U/ml), B (17.0U/ml), and C (37.0U/ml). Biliary malignancy patients in Groups B and C had significantly higher CA19-9 values than cancer-free patients (p0.001). The optimal cut-off, as determined by ROC analysis, for all patients was 88.5U/ml. Optimal cut-off values in Groups A, B, and C were 4.0U/ml, 74.5U/ml, and 106.8U/ml, respectively. Use of these values improved sensitivity of CA19-9 in Groups B and C. Further, use of group-dependent cut-off values with 90% sensitivity resulted in a 42.9% reduction of false positive results.Use of FUT2/3 genotype-dependent cut-off values for CA19-9 improved sensitivity and reduced the number of false positive results.

Passerine birds show large interspecific variation in extrapair paternity rates. There is accumulating evidence that such promiscuous behavior is driven by indirect, genetic benefits to females. Sexual selection theory distinguishes between two types of genetic benefits, additive and nonadditive effects, mediated by preferences for good and compatible genes, respectively. Good genes preferences should imply directional selection and mating skew among males, and thus reduced genetic diversity in the population. In contrast, compatible genes preferences should give balancing selection that retains genetic diversity. Here, we test how well these predictions fit with patterns of variation in genetic diversity and promiscuity levels among passerine birds. We found that more promiscuous species had higher nucleotide diversity at autosomal introns, but not at Z-chromosome introns. We also found that major histocompatibility complex (MHC) class IIB alleles had higher sequence diversity, and therefore should recognize a broader spectrum of pathogens, in more promiscuous species. Our results suggest that female promiscuity targets a multitude of autosomal genes for their nonadditive, compatibility benefits. Also, as immunity genes seem to be of particular importance, we hypothesize that interspecific variation in female promiscuity among passerine birds has arisen in response to the strength of pathogen-mediated selection.

Post-operative liver regeneration may contribute to tumor recurrence. There is a theoretical need for an adjuvant therapy that can suppress tumor growth without adversely affecting post-operative liver regeneration.To evaluate the effect of RAF inhibitor Sorafenib on cell viability and proliferation of hepatoma cells and hepatocytes in vitro and in an in vivo rat model.Cell viability, DNA synthesis, and RAF/MAPK kinase activity in the primary hepatocyte and hepatoma cell lines were investigated after Sorafenib exposure. Sequence analysis of the B-RAF gene in hepatic cells was determined. Tumor markers were compared within the rats after 70% hepatectomy with or without daily oral gavages of Sorafenib. Liver regeneration was assessed by liver function tests and proliferation markers.Primary hepatocytes showed higher cell viability, proliferation rate, and stronger RAF/MAPK kinase activity compared with hepatoma cell lines. The in vivo tumor volumes, size, and metastases were significantly decreased (P 0.05) whereas no significant change in liver regeneration related to Sorafenib exposure was found (P 0.05). B-RAF V600E mutation was not detected neither in the hepatic cells nor untransformed hepatocytes.The RAF targeted inhibitor can reduce tumor growth without retarding liver regeneration in this experiment.

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., his study was funded by the Spanish Ministry of Education, Science and Culture (projects CGL2006-07481/BOS and CGL2009-10652/BOS), which also granted the PhD fellowship to DC (I3P-BDP2005) and the Natural History Museum, University of Oslo (PhD fellowship to JAA).

Background Microsatellites are frequently used genetic markers in a wide range of applications, primarily due to their high length polymorphism levels that can easily be genotyped by fragment length analysis. However, the mode of microsatellite evolution is yet not fully understood, and the role of interrupting motifs for the stability of microsatellites remains to be explored in more detail. Here we present a sequence analysis of mutation events and a description of the structure of repeated regions in the hypervariable, pentanucleotide microsatellite locus HrU10 in barn swallows (Hirundo rustica) and tree swallows (Tachycineta bicolor). Results In a large-scale parentage analysis in barn swallows and tree swallows, broods were screened for mutations at the HrU10 locus. In 41 cases in the barn swallows and 15 cases in the tree swallows, mutations corresponding to the loss or gain of one or two repeat units were detected. The parent and mutant offspring alleles were sequenced for 33 of these instances (26 in barn swallows and 7 in tree swallows). Replication slippage was considered the most likely mutational process. We tested the hypothesis that HrU10, a microsatellite with a wide allele size range, has an increased probability of introductions of interruptive motifs (IMs) with increasing length of the repeated region. Indeed, the number and length of the IMs was strongly positively correlated with the total length of the microsatellite. However, there was no significant correlation with the length of the longest stretch of perfectly repeated units, indicating a threshold level for the maximum length of perfectly repeated pentanucleotide motifs in stable HrU10 alleles. The combination of sequence and pedigree data revealed that 15 barn swallow mutations (58%) produced alleles that were size homoplasic to other alleles in the data set. Conclusion Our results give further insights into the mode of microsatellite evolution, and support the assumption of increased slippage rate with increased microsatellite length and a stabilizing effect of interrupting motifs for microsatellite regions consisting of perfect repeats. In addition, the observed extent of size homoplasy may impose a general caution against using hypervariable microsatellites in genetic diversity measures when alleles are identified by fragment length analysis only.

The human leukocyte antigen (HLA) complex contains the most polymorphic genes in the human genome. The classical HLA class I and II genes define the specificity of adaptive immune responses. Genetic variation at the HLA genes is associated with susceptibility to autoimmune and infectious diseases and plays a major role in transplantation medicine and immunology. Currently, the HLA genes are characterized using Sanger- or next-generation sequencing (NGS) of a limited amplicon repertoire or labeled oligonucleotides for allele-specific sequences. High-quality NGS-based methods are in proprietary use and not publicly available. Here, we introduce the first highly automated open-kit/open-source HLA-typing method for NGS. The method employs in-solution targeted capturing of the classical class I (HLA-A, HLA-B, HLA-C) and class II HLA genes (HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). The calling algorithm allows for highly confident allele-calling to three-field resolution (cDNA nucleotide variants). The method was validated on 357 commercially available DNA samples with known HLA alleles obtained by classical typing. Our results showed on average an accurate allele call rate of 0.99 in a fully automated manner, identifying also errors in the reference data. Finally, our method provides the flexibility to add further enrichment target regions.