Killer-cell Immunoglobulin-like Receptor gene linkage and copy number variation analysis by droplet digital PCR

KIR genotype and haplotype repertoire in Kuwaiti healthy donors, hematopoietic cell transplant recipients and healthy family members

The gene complex located on chromosome 19q13.4 encodes the Killer-cell Immunoglobulin-like Receptors (KIRs), which exhibit remarkable polymorphism in both gene content and sequences. Further, the repertoire of KIR genes varies within and between populations, creating a diverse pool of KIR genotypes. This study was carried out to characterize KIR genotypes and haplotypes among 379 Arab Kuwaiti individuals including 60 subjects from 20 trio families, 49 hematopoietic cell transplantation (HCT) recipients and 270 healthy Kuwaiti volunteer HCT donors. KIR Genotyping was performed by a combination of reverse sequence specific oligonucleotide probes (rSSO) and/or Real Time PCR. The frequencies of KIR genes in 270 healthy Kuwaiti volunteer donors were compared to previously reported frequencies in other populations. In addition, we compared the differences in KIR repertoire of patients and healthy donors to investigate the reproducibility of previously reported significant differences between patients with hematological malignancies and healthy donors. The observed frequencies in our cohort volunteer HCT donors was comparable to those reported in neighboring Arab populations. The activating genes KIR2DS1, KIR2DS5 and KIR3DS1 and the inhibitory gene KIR2DL5 were significantly more frequent in patients compared to healthy donors, however, none of the previously reported differences were reproducible in our Kuwaiti cohort. This report is the first description of KIR gene carrier frequency and haplotype characterization in a fairly large cohort of the Kuwaiti population, which may have implications in KIR based HCT donor selection strategies.

Current commercial dPCR platforms: technology and market review

Digital polymerase chain reaction (dPCR) technology has provided a new technique for molecular diagnostics, with superior advantages, such as higher sensitivity, precision, and specificity over quantitative real-time PCRs (qPCR). Eight companies have offered commercial dPCR instruments: Fluidigm Corporation, Bio-Rad, RainDance Technologies, Life Technologies, Qiagen, JN MedSys Clarity, Optolane, and Stilla Technologies Naica. This paper discusses the working principle of each offered dPCR device and compares the associated: technical aspects, usability, costs, and current applications of each dPCR device. Lastly, up-and-coming dPCR technologies are also presented, as anticipation of how the dPCR device landscape may likely morph in the next few years.

Development and Evaluation of a Droplet Digital PCR Assay for 8p23 β-Defensin Cluster Copy Number Determination

Introduction

It has been proposed that the copy number (CN) variation (CNV) in β-defensin genes (DEFB) on human chromosome 8p23 determines phenotypic differences in inflammatory diseases. However, no method for accurate and easy DEFB CN quantification is yet available.

Objective

Droplet digital polymerase chain reaction (ddPCR) is a novel method for CNV quantification and has been used for genes such as CCL4L, CCL3L1, AMY1, and HER2. However, to date, no ddPCR assay has been available for DEFB CN determination. In the present study, we aimed to develop and evaluate such a ddPCR assay.

Methods

The assay was designed using DEFB4 and RPP30 as the target and the reference gene, respectively. To evaluate the assay, 283 DNA samples with known CNs previously determined using the multiple ligation-dependent probe amplification (MLPA) method, the current gold standard, were used as standards. To discover the optimal DNA template amount, we tested 80 to 2.5 ng DNA by a serial of one to two dilutions of eight samples. To evaluate the reproducibility of the assay, 31 samples were repeated to calculate the intra- and inter-assay variations. To further validate the reliability of the assay, the CNs of all 283 samples were determined using ddPCR. To compare results with those using quantitative PCR (qPCR), DEFB CNs for 48 samples were determined using qPCR with the same primers and probes.

Results

In a one-dimensional plot, the positive and negative droplets were clearly separated in both DEFB4 and RPP30 detection channels. In a two-dimensional plot, four populations of droplets were observed. The 20 ng template DNA proved optimal, with either high (80 ng) or low (10, 5, 2.5 ng) template input leading to ambiguous or inaccurate results. For the 31 standard samples, DEFB CNs were accurately determined with small intra- and inter-assay variations (coefficient of variation < 0.04 for both). In the validation cohort, ddPCR provided the correct CN for all 283 samples with high confidence. qPCR measurements for the 48 samples produced noisy data with high uncertainty and low accuracy.

Conclusions

ddPCR is an accurate, reproducible, easy-to-use, cheap, high-throughput method for DEFB CN determination. ddPCR could be applied to DEFB CN quantification in large-scale case–control studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40291-021-00546-2.

HLA Class I Molecules as Immune Checkpoints for NK Cell Alloreactivity and Anti-Viral Immunity in Kidney Transplantation

Natural killer (NK) cells are innate lymphocytes that can kill diseased- or virally-infected cells, mediate antibody dependent cytotoxicity and produce type I immune-associated cytokines upon activation. NK cells also contribute to the allo-immune response upon kidney transplantation either by promoting allograft rejection through lysis of cells of the transplanted organ or by promoting alloreactive T cells. In addition, they protect against viral infections upon transplantation which may be especially relevant in patients receiving high dose immune suppression. NK cell activation is tightly regulated through the integrated balance of signaling via inhibitory- and activating receptors. HLA class I molecules are critical regulators of NK cell activation through the interaction with inhibitory- as well as activating NK cell receptors, hence, HLA molecules act as critical immune checkpoints for NK cells. In the current review, we evaluate how NK cell alloreactivity and anti-viral immunity are regulated by NK cell receptors belonging to the KIR family and interacting with classical HLA class I molecules, or by NKG2A/C and LILRB1/KIR2DL4 engaging non-classical HLA-E or -G. In addition, we provide an overview of the methods to determine genetic variation in these receptors and their HLA ligands.

Long-Distance Phasing of a Tentative "Enhancer" Single-Nucleotide Polymorphism With CYP2D6 Star Allele Definitions

BACKGROUND

The CYP2D6 gene locus has been extensively studied over decades, yet a portion of variability in CYP2D6 activity cannot be explained by known sequence variations within the gene, copy number variation, or structural rearrangements. It was proposed that rs5758550, located 116 kb downstream of the CYP2D6 gene locus, increases gene expression and thus contributes to variability in CYP2D6 activity. This finding has, however, not been validated. The purpose of the study was to address a major technological barrier, i.e., experimentally linking rs5758550, also referred to as the "enhancer" single-nucleotide polymorphism (SNP), to CYP2D6 haplotypes >100 kb away. To overcome this challenge is essential to ultimately determine the contribution of the "enhancer" SNP to interindividual variability in CYP2D6 activity.

METHODS

A large ethnically mixed population sample (n=3,162) was computationally phased to determine linkage between the "enhancer" SNP and CYP2D6 haplotypes (or star alleles). To experimentally validate predicted linkages, DropPhase2D6, a digital droplet PCR (ddPCR)-based method was developed. 10X Genomics Linked-Reads were utilized as a proof of concept.

RESULTS

Phasing predicted that the "enhancer" SNP can occur on numerous CYP2D6 haplotypes including CYP2D6*1, *2, *5, and *41 and suggested that linkage is incomplete, i.e., a portion of these alleles do not have the "enhancer" SNP. Phasing also revealed differences among the European and African ancestry data sets regarding the proportion of alleles with and without the "enhancer" SNP. DropPhase2D6 was utilized to confirm or refute the predicted "enhancer" SNP location for individual samples, e.g., of n=3 samples genotyped as *1/*41, rs5758550 was on the *41 allele of two samples and on the *1 allele of one sample. Our findings highlight that the location of the "enhancer" SNP must not be assigned by "default." Furthermore, linkage between the "enhancer" SNP and CYP2D6 star allele haplotypes was confirmed with 10X Genomics technology.

CONCLUSIONS

Since the "enhancer" SNP can be present on a portion of normal, decreased, or no function alleles, the phase of the "enhancer" SNP must be considered when investigating the impact of the "enhancer" SNP on CYP2D6 activity.

Haplotype-Based Analysis of KIR-Gene Profiles in a South European Population-Distribution of Standard and Variant Haplotypes, and Identification of Novel Recombinant Structures

Inhibitory Killer-cell Immunoglobulin-like Receptors (KIR) specific for HLA class I molecules enable human natural killer cells to monitor altered antigen presentation in pathogen-infected and tumor cells. KIR genes display extensive copy-number variation and allelic polymorphism. They organize in a series of variable arrangements, designated KIR haplotypes, which derive from duplications of ancestral genes and sequence diversification through point mutation and unequal crossing-over events. Genomic studies have established the organization of multiple KIR haplotypes—many of them are fixed in most human populations, whereas variants of those have less certain distributions. Whilst KIR-gene diversity of many populations and ethnicities has been explored superficially (frequencies of individual genes and presence/absence profiles), less abundant are in-depth analyses of how such diversity emerges from KIR-haplotype structures. We characterize here the genetic diversity of KIR in a sample of 414 Spanish individuals. Using a parsimonious approach, we manage to explain all 38 observed KIR-gene profiles by homo- or heterozygous combinations of six fixed centromeric and telomeric motifs; of six variant gene arrangements characterized previously by us and others; and of two novel haplotypes never detected before in Caucasoids. Associated to the latter haplotypes, we also identified the novel transcribed KIR2DL5B^*0020202 allele, and a chimeric KIR2DS2/KIR2DL3 gene (designated KIR2DL3^*033) that challenges current criteria for classification and nomenclature of KIR genes and haplotypes.

Conjunctival Microbiome-Host Responses Are Associated With Impaired Epithelial Cell Health in Both Early and Late Stages of Trachoma

Background: Trachoma, a neglected tropical disease, is the leading infectious cause of blindness and visual impairment worldwide. Host responses to ocular chlamydial infection resulting in chronic inflammation and expansion of non-chlamydial bacteria are hypothesized risk factors for development of active trachoma and conjunctival scarring.

Methods: Ocular swabs from trachoma endemic populations in The Gambia were selected from archived samples for 16S sequencing and host conjunctival gene expression. We recruited children with active trachoma and adults with conjunctival scarring, alongside corresponding matched controls.

Findings: In children, active trachoma was not associated with significant changes in the ocular microbiome. Haemophilus enrichment was associated with antimicrobial responses but not linked to active trachoma. Adults with scarring trachoma had a reduced ocular bacterial diversity compared to controls, with increased relative abundance of Corynebacterium. Increased abundance of Corynebacterium in scarring disease was associated with innate immune responses to the microbiota, dominated by altered mucin expression and increased matrix adhesion.

Interpretation: In the absence of current Chlamydia trachomatis infection, changes in the ocular microbiome associate with differential expression of antimicrobial and inflammatory genes that impair epithelial cell health. In scarring trachoma, expansion of non-pathogenic bacteria such as Corynebacterium and innate responses are coincident, warranting further investigation of this relationship. Comparisons between active and scarring trachoma supported the relative absence of type-2 interferon responses in scarring, whilst highlighting a common suppression of re-epithelialization with altered epithelial and bacterial adhesion, likely contributing to development of scarring pathology.

Killer Ig-Like Receptors (KIRs): Their Role in NK Cell Modulation and Developments Leading to Their Clinical Exploitation

Natural killer (NK) cells contribute to the first line of defense against viruses and to the control of tumor growth and metastasis spread. The discovery of HLA class I specific inhibitory receptors, primarily of killer Ig-like receptors (KIRs), and of activating receptors has been fundamental to unravel NK cell function and the molecular mechanisms of tumor cell killing. Stemmed from the seminal discoveries in early '90s, in which Alessandro Moretta was the major actor, an extraordinary amount of research on KIR specificity, genetics, polymorphism, and repertoire has followed. These basic notions on NK cells and their receptors have been successfully translated to clinical applications, primarily to the haploidentical hematopoietic stem cell transplantation to cure otherwise fatal leukemia in patients with no HLA compatible donors. The finding that NK cells may express the PD-1 inhibitory checkpoint, particularly in cancer patients, may allow understanding how anti-PD-1 therapy could function also in case of HLA class I^neg tumors, usually susceptible to NK-mediated killing. This, together with the synergy of therapeutic anti-checkpoint monoclonal antibodies, including those directed against NKG2A or KIRs, emerging in recent or ongoing studies, opened new solid perspectives in cancer therapy.

Determination of copy number and circularization ratio of Tn916-Tn1545 family of conjugative transposons in oral streptococci by droplet digital PCR

Background: Tn916 and Tn1545 are paradigms of a large family of related, broad host range, conjugative transposons that are widely distributed in bacteria and contribute to the spread of antibiotic resistance genes (ARGs). Variation in the copy number (CN) of Tn916-Tn1545 elements and the circularization ratio (CR) may play an important role in propagation of ARGs carried by these elements.

Objectives and Design: In this study, the CN and CR of Tn916-Tn1545 elements in oral streptococci were determined using droplet digital PCR (ddPCR). In addition, we investigated the influence of tetracycline on the CR of Tn916-Tn1545 elements.

Results: The ddPCR assay designed in this study is a reliable way to rapidly determine CN and CR of Tn916-Tn1545 elements.

Conclusions: Our data also suggest that Tn916-Tn1545 elements are generally stable without selective pressure in the clinical oral Streptococcus strains investigated in this study.

Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

dPCR: A Technology Review

Digital Polymerase Chain Reaction (dPCR) is a novel method for the absolute quantification of target nucleic acids. Quantification by dPCR hinges on the fact that the random distribution of molecules in many partitions follows a Poisson distribution. Each partition acts as an individual PCR microreactor and partitions containing amplified target sequences are detected by fluorescence. The proportion of PCR-positive partitions suffices to determine the concentration of the target sequence without a need for calibration. Advances in microfluidics enabled the current revolution of digital quantification by providing efficient partitioning methods. In this review, we compare the fundamental concepts behind the quantification of nucleic acids by dPCR and quantitative real-time PCR (qPCR). We detail the underlying statistics of dPCR and explain how it defines its precision and performance metrics. We review the different microfluidic digital PCR formats, present their underlying physical principles, and analyze the technological evolution of dPCR platforms. We present the novel multiplexing strategies enabled by dPCR and examine how isothermal amplification could be an alternative to PCR in digital assays. Finally, we determine whether the theoretical advantages of dPCR over qPCR hold true by perusing studies that directly compare assays implemented with both methods.

Normalization of cell associated antiretroviral drug concentrations with a novel RPP30 droplet digital PCR assay

Quantification of antiretroviral (ARV) drug concentrations in peripheral blood mononuclear cells (PBMCs) and tissue isolated mononuclear cells (TIMCs) from lymph node (LNMC) and rectum (RMC) is an important measure of bio-distribution. Normalization of drug concentrations is critical to represent tissue drug concentrations and to analyze both intra-individual and inter-individual variability in drug distribution. However, a molecular method to normalize intracellular drug concentrations in PBMCs and TIMCs methanol extracts is currently unavailable. In this study, a novel droplet digital PCR (ddPCR) assay was designed to amplify RPP30 gene sequence conserved in human and non-human primates (NHP). Genomic DNA (gDNA) isolated from 70 percent methanol embedded PBMCs and TIMCs was used as ddPCR template to quantitate precise RPP30 copies to derive cell counts. The novel molecular method quantitated RPP30 copies in human and rhesus macaque gDNA templates with greater accuracy and precision than qPCR. RPP30 ddPCR derived cell counts are strongly correlated with automated cytometer based cell counts in PBMC (R = 0.90, p = 0.001 and n = 20); LNMC (R = 0.85 p = 0.0001 and n = 22) and RMC (R = 0.92, p = 0.0001 and n = 20) and achieved comparable normalized drug concentrations. Therefore, the RPP30 ddPCR assay is an important normalization method in drug bio-distribution and pharmacokinetic studies in humans and NHPs.

A scalable self-priming fractal branching microchannel net chip for digital PCR

As an absolute quantification method at the single-molecule level, digital PCR has been widely used in many bioresearch fields, such as next generation sequencing, single cell analysis, gene editing detection and so on. However, existing digital PCR methods still have some disadvantages, including high cost, sample loss, and complicated operation. In this work, we develop an exquisite scalable self-priming fractal branching microchannel net digital PCR chip. This chip with a special design inspired by natural fractal-tree systems has an even distribution and 100% compartmentalization of the sample without any sample loss, which is not available in existing chip-based digital PCR methods. A special 10 nm nano-waterproof layer was created to prevent the solution from evaporating. A vacuum pre-packaging method called self-priming reagent introduction is used to passively drive the reagent flow into the microchannel nets, so that this chip can realize sequential reagent loading and isolation within a couple of minutes, which is very suitable for point-of-care detection. When the number of positive microwells stays in the range of 100 to 4000, the relative uncertainty is below 5%, which means that one panel can detect an average of 101 to 15 374 molecules by the Poisson distribution. This chip is proved to have an excellent ability for single molecule detection and quantification of low expression of hHF-MSC stem cell markers. Due to its potential for high throughput, high density, low cost, lack of sample and reagent loss, self-priming even compartmentalization and simple operation, we envision that this device will significantly expand and extend the application range of digital PCR involving rare samples, liquid biopsy detection and point-of-care detection with higher sensitivity and accuracy.

Copy Number Variation Analysis by Droplet Digital PCR

The health impact of many copy number variants in our genome remains still largely to be discovered. Detecting and genotyping this often complex variation presents a technical challenge. Here we describe a 96-well format droplet digital PCR (ddPCR) protocol for genotyping a common copy variant in the human haptoglobin gene. ddPCR allows for high-throughput and accurate quantitation of gene copy numbers.

Deciphering the killer-cell immunoglobulin-like receptor system at super-resolution for natural killer and T-cell biology

Killer-cell immunoglobulin-like receptors (KIRs) are components of two fundamental biological systems essential for human health and survival. First, they contribute to host immune responses, both innate and adaptive, through their expression by natural killer cells and T cells. Second, KIR play a key role in regulating placentation, and hence reproductive success. Analogous to the diversity of their human leucocyte antigen class I ligands, KIR are extremely polymorphic. In this review, we describe recent developments, fuelled by methodological advances, that are helping to decipher the KIR system in terms of haplotypes, polymorphisms, expression patterns and their ligand interactions. These developments are delivering deeper insight into the relevance of KIR in immune system function, evolution and disease.

qKAT: a high-throughput qPCR method for KIR gene copy number and haplotype determination

Killer cell immunoglobulin-like receptors (KIRs), expressed on natural killer cells and T cells, have considerable biomedical relevance playing significant roles in immunity, pregnancy and transplantation. The KIR locus is one of the most complex and polymorphic regions of the human genome. Extensive sequence homology and copy number variation makes KIRs technically laborious and expensive to type. To aid the investigation of KIRs in human disease we developed a high-throughput, multiplex real-time polymerase chain reaction method to determine gene copy number for each KIR locus. We used reference DNA samples to validate the accuracy and a cohort of 1698 individuals to evaluate capability for precise copy number discrimination. The method provides improved information and identifies KIR haplotype alterations that were not previously visible using other approaches.

Surveillance of cell and tissue perturbation by receptors in the LRC

The human leukocyte receptor complex (LRC) encompasses several sets of genes with a common evolutionary origin and which form a branch of the immunoglobulin superfamily (IgSF). Comparisons of LRC genes both within and between species calls for a high degree of plasticity. The drive for this unprecedented level of variation is not known, but it relates in part to interaction of several LRC products with polymorphic human leukocyte antigen (HLA) class I molecules. However, the range of other proposed ligands for LRC products indicates a dynamic set of receptors that have adapted to detect target molecules relating to numerous cellular pathways. Several receptors in the complex bind a molecular signature in collagenous ligands. Others detect a variety of motifs relating to pathogens in addition to cellular stress, attesting to the opportunistic versatility of LRC receptors.

Identification of genomic aberrations in hemangioblastoma by droplet digital PCR and SNP microarray highlights novel candidate genes and pathways for pathogenesis

Background

The genetic mechanisms underlying hemangioblastoma development are still largely unknown. We used high-resolution single nucleotide polymorphism microarrays and droplet digital PCR analysis to detect copy number variations (CNVs) in total of 45 hemangioblastoma tumors.

Results

We identified 94 CNVs with a median of 18 CNVs per sample. The most frequently gained regions were on chromosomes 1 (p36.32) and 7 (p11.2). These regions contain the EGFR and PRDM16 genes. Recurrent losses were located at chromosome 12 (q24.13), which includes the gene PTPN11.

Conclusions

Our findings provide the first high-resolution genome-wide view of chromosomal changes in hemangioblastoma and identify 23 candidate genes: EGFR, PRDM16, PTPN11, HOXD11, HOXD13, FLT3, PTCH, FGFR1, FOXP1, GPC3, HOXC13, HOXC11, MKL1, CHEK2, IRF4, GPHN, IKZF1, RB1, HOXA9, and micro RNA, such as hsa-mir-196a-2 for hemangioblastoma pathogenesis. Furthermore, our data implicate that cell proliferation and angiogenesis promoting pathways may be involved in the molecular pathogenesis of hemangioblastoma.

Genomic-assisted phylogenetic analysis and marker development for next generation soybean cyst nematode resistance breeding

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is a serious soybean pest. The use of resistant cultivars is an effective approach for preventing yield loss. In this study, 19,652 publicly available soybean accessions that were previously genotyped with the SoySNP50K iSelect BeadChip were used to evaluate the phylogenetic diversity of SCN resistance genes Rhg1 and Rhg4 in an attempt to identify novel sources of resistance. The sequence information of soybean lines was utilized to develop KASPar (KBioscience Competitive Allele-Specific PCR) assays from single nucleotide polymorphisms (SNPs) of Rhg1, Rhg4, and other novel quantitative trait loci (QTL). These markers were used to genotype a diverse set of 95 soybean germplasm lines and three recombinant inbred line (RIL) populations. SNP markers from the Rhg1 gene were able to differentiate copy number variation (CNV), such as resistant-high copy (PI 88788-type), low copy (Peking-type), and susceptible-single copy (Williams 82) numbers. Similarly, markers for the Rhg4 gene were able to detect Peking-type (resistance) genotypes. The phylogenetic information of SCN resistance loci from a large set of soybean accessions and the gene/QTL specific markers that were developed in this study will accelerate SCN resistance breeding programs.

Association of Killer Cell Immunoglobulin- Like Receptor Genes in Iranian Patients with Rheumatoid Arthritis

Objectives

Rheumatoid arthritis (RA) is a chronic inflammatory disorder characterized by persistent synovitis, ultimately leading to cartilage and bone degeneration. Natural Killer cells and CD28 null T-cells are suspected as role players in RA pathogenesis. These cells are similar in feature and function, as they both exert their cytotoxic effect via Killer Cell Immunoglobulin- Like Receptors (KIR) on their surface. KIR genes have either an inhibitory or activating effect depending on their intracytoplasmic structure. Herein we genotyped 16 KIR genes, 3 pseudo genes and 6 HLA class І genes as their corresponding ligands in RA patients and control subjects.

Methods

In this case-control study, KIR and HLA genes were genotyped in 400 RA patients and 372 matched healthy controls using sequence-specific primers (SSP-PCR). Differences in the frequency of genes and haplotypes were determined by χ² test.

Results

KIR2DL2, 2DL5a, 2DL5b and activating KIR: KIR2DS5 and 3DS1 were all protective against RA. KIR2DL5 removal from a full Inhibitory KIR haplotype converted the mild protection (OR = 0.56) to a powerful predisposition to RA (OR = 16.47). Inhibitory haplotype No. 7 comprising KIR2DL5 in the absence of KIR2DL1 and KIR2DL3 confers a 14-fold protective effect against RA.

Conclusion

Individuals carrying the inhibitory KIR haplotype No. 6 have a high potential risk for developing RA.

Identification of Novel Fusion Genes in Testicular Germ Cell Tumors

Testicular germ cell tumors (TGCT) are the most frequently diagnosed solid tumors in young men ages 15 to 44 years. Embryonal carcinomas (EC) comprise a subset of TGCTs that exhibit pluripotent characteristics similar to embryonic stem (ES) cells, but the genetic drivers underlying malignant transformation of ECs are unknown. To elucidate the abnormal genetic events potentially contributing to TGCT malignancy, such as the existence of fusion genes or aberrant fusion transcript expression, we performed RNA sequencing of EC cell lines and their nonmalignant ES cell line counterparts. We identified eight novel fusion transcripts and one gene with alternative promoter usage, ETV6. Four out of nine transcripts were found recurrently expressed in an extended panel of primary TGCTs and additional EC cell lines, but not in normal parenchyma of the testis, implying tumor-specific expression. Two of the recurrent transcripts involved an intrachromosomal fusion between RCC1 and HENMT1 located 80 Mbp apart and an interchromosomal fusion between RCC1 and ABHD12B. RCC1-ABHD12B and the ETV6 transcript variant were found to be preferentially expressed in the more undifferentiated TGCT subtypes. In vitro differentiation of the NTERA2 EC cell line resulted in significantly reduced expression of both fusion transcripts involving RCC1 and the ETV6 transcript variant, indicating that they are markers of pluripotency in a malignant setting. In conclusion, we identified eight novel fusion transcripts that, to our knowledge, are the first fusion genes described in TGCT and may therefore potentially serve as genomic biomarkers of malignant progression.

Imputation of KIR Types from SNP Variation Data

Large population studies of immune system genes are essential for characterizing their role in diseases, including autoimmune conditions. Of key interest are a group of genes encoding the killer cell immunoglobulin-like receptors (KIRs), which have known and hypothesized roles in autoimmune diseases, resistance to viruses, reproductive conditions, and cancer. These genes are highly polymorphic, which makes typing expensive and time consuming. Consequently, despite their importance, KIRs have been little studied in large cohorts. Statistical imputation methods developed for other complex loci (e.g., human leukocyte antigen [HLA]) on the basis of SNP data provide an inexpensive high-throughput alternative to direct laboratory typing of these loci and have enabled important findings and insights for many diseases. We present KIR^∗IMP, a method for imputation of KIR copy number. We show that KIR^∗IMP is highly accurate and thus allows the study of KIRs in large cohorts and enables detailed investigation of the role of KIRs in human disease.

Complex and multi-allelic copy number variation in human disease

Hundreds of copy number variants are complex and multi-allelic, in that they have many structural alleles and have rearranged multiple times in the ancestors who contributed chromosomes to current humans. Not only are the relationships of these multi-allelic CNVs (mCNVs) to phenotypes generally unknown, but many mCNVs have not yet been described at the basic levels—alleles, allele frequencies, structural features—that support genetic investigation. To date, most reported disease associations to these variants have been ascertained through candidate gene studies. However, only a few associations have reached the level of acceptance defined by durable replications in many cohorts. This likely stems from longstanding challenges in making precise molecular measurements of the alleles individuals have at these loci. However, approaches for mCNV analysis are improving quickly, and some of the unique characteristics of mCNVs may assist future association studies. Their various structural alleles are likely to have different magnitudes of effect, creating a natural allelic series of growing phenotypic impact and giving investigators a set of natural predictions and testable hypotheses about the extent to which each allele of an mCNV predisposes to a phenotype. Also, mCNVs’ low-to-modest correlation to individual single-nucleotide polymorphisms (SNPs) may make it easier to distinguish between mCNVs and nearby SNPs as the drivers of an association signal, and perhaps, make it possible to preliminarily screen candidate loci, or the entire genome, for the many mCNV–disease relationships that remain to be discovered.

A droplet digital PCR detection method for rare L1 insertions in tumors

Background

The active human mobile element, long interspersed element 1 (L1) currently populates human genomes in excess of 500,000 copies per haploid genome. Through its mobility via a process called target primed reverse transcription (TPRT), L1 mobilization has resulted in over 100 de novo cases of human disease and has recently been associated with various cancer types. Large advances in high-throughput sequencing (HTS) technology have allowed for an increased understanding of the role of L1 in human cancer; however, researchers are still limited by the ability to validate potentially rare L1 insertion events detected by HTS that may occur in only a small fraction of tumor cells. Additionally, HTS detection of rare events varies greatly as a function of read depth, and new tools for de novo element discovery are needed to fill in gaps created by HTS.

Results

We have employed droplet digital PCR (ddPCR) to detect rare L1 loci in mosaic human genomes. Our assay allows for the detection of L1 insertions as rare as one cell in every 10,000.

Conclusions

ddPCR represents a robust method to be used alongside HTS techniques for detecting, validating and quantitating rare L1 insertion events in tumors and other tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-014-0030-4) contains supplementary material, which is available to authorized users.

Germ-line transmitted, chromosomally integrated HHV-6 and classical Hodgkin lymphoma

A unique feature of both human herpesvirus 6A and B (HHV-6A and B) among human herpesviruses is their ability to integrate into chromosomal telomeres. In some individuals integrated viral genomes are present in the germ-line and result in the vertical transmission of HHV-6; however, little is known about the disease associations of germ-line transmitted, chromosomally integrated HHV-6 (ciHHV-6). Recent publications suggest that HHV-6 is associated with classical Hodgkin lymphoma (cHL). Here we examine the prevalence of ciHHV-6 in 936 cases of cHL and 563 controls by screening with a duplex TaqMan assay and confirming with droplet digital PCR. ciHHV-6 was detected in 10/563 (1.8%) controls and in all but one individual the virus was HHV-6B. Amongst cases 16/936 (1.7%) harboured ciHHV-6, thus demonstrating no association between ciHHV-6 and risk of cHL.

Considerations for digital PCR as an accurate molecular diagnostic tool

BACKGROUND

Digital PCR (dPCR) is an increasingly popular manifestation of PCR that offers a number of unique advantages when applied to preclinical research, particularly when used to detect rare mutations and in the precise quantification of nucleic acids. As is common with many new research methods, the application of dPCR to potential clinical scenarios is also being increasingly described.

CONTENT

This review addresses some of the factors that need to be considered in the application of dPCR. Compared to real-time quantitative PCR (qPCR), dPCR clearly has the potential to offer more sensitive and considerably more reproducible clinical methods that could lend themselves to diagnostic, prognostic, and predictive tests. But for this to be realized the technology will need to be further developed to reduce cost and simplify application. Concomitantly the preclinical research will need be reported with a comprehensive understanding of the associated errors. dPCR benefits from a far more predictable variance than qPCR but is as susceptible to upstream errors associated with factors like sampling and extraction. dPCR can also suffer systematic bias, particularly leading to underestimation, and internal positive controls are likely to be as important for dPCR as they are for qPCR, especially when reporting the absence of a sequence.

SUMMARY

In this review we highlight some of the considerations that may be needed when applying dPCR and discuss sources of error. The factors discussed here aim to assist in the translation of dPCR to diagnostic, predictive, or prognostic applications.