Nimbus: a design-driven analyses suite for amplicon-based NGS data.

Motivation: PCR-based DNA enrichment followed by massively parallel sequencing is a straightforward and cost effective method to sequence genes up to high depth. The full potential of amplicon-based sequencing assays is currently not achieved as analysis methods do not take into account the source amplicons of the detected variants. Tracking the source amplicons has the potential to identify systematic biases, enhance variant calling and improve the designs of future assays. Results: We present Nimbus, a software suite for the analysis of amplicon-based sequencing data. Nimbus includes tools for data pre-processing, alignment, single nucleotide polymorphism (SNP), insertion and deletion calling, quality control and visualization. Nimbus can detect SNPs in its alignment seeds and reduces alignment issues by the usage of decoy amplicons. Tracking the amplicons throughout analysis allows easy and fast design optimization by amplicon performance comparison. It enables detection of probable false positive variants present in a single amplicon from real variants present in multiple amplicons and provides multiple sample visualization. Nimbus was tested using HaloPlex Exome datasets and outperforms other callers for low-frequency variants. The variants called by Nimbus were highly concordant between twin samples and SNP-arrays. The Nimbus suite provides an end-to-end solution for variant calling, design optimization and visualization of amplicon-derived next-generation sequencing datasets. Availability and implementation: https://github.com/erasmus-center-for-biomics/Nimbus. Supplementary information: Supplementary data are available at Bioinformatics online.

Patient survival in uveal melanoma is not affected by oncogenic mutations in GNAQ and GNA11.

BACKGROUND: Mutations in GNAQ and GNA11, encoding the oncogenic G-protein alpha subunit q and 11, respectively, occur frequently in the majority of uveal melanomas. METHODS: Exons 4 and 5 from GNAQ and GNA11 were amplified and sequenced from 92 ciliary body and choroidal melanomas. The mutation status was correlated with disease-free survival (DFS) and other parameters. RESULTS: None of the tumours harboured a GNAQ exon 4 mutation. A GNAQ mutation in exon 5 codon 209 was found in 46 out of 92 (50.0%) of the tumours. Only 1 out of 92 (1.1%) melanomas showed a mutation in GNA11 exon 4 codon 183, whereas 39 out of 92 (42.4%) harboured a mutation in exon 5 of GNA11 codon 209. Six tumours did not show any mutations in exons 4 and 5 of these genes. Univariate analyses showed no correlation between DFS and the mutation status. CONCLUSION: GNAQ and GNA11 mutations are, in equal matter, not associated with patient outcome.

Phenotype-genotype correlation in a familial IGF1R microdeletion case.

BACKGROUND: IGF1R (insulin-like growth factor 1 receptor) haploinsufficiency is a rare event causing difficulties in defining clear genotype-phenotype correlations, although short stature is its well established hallmark. Several pure 15q26 monosomies (n=22) have been described in the literature, including those with breakpoints proximal to the IGF1R gene. Clinical heterogeneity is characteristic for these mainly de novo telomeric deletions and is illustrated by the involvement of several different organ systems such as the heart, diaphragm, lungs, kidneys and limbs, besides growth failure in the patient's phenotype. The clinical variability in these patients could be explained by the haploinsufficiency of multiple genes besides the IGF1R gene. In comparison, the six different IGF1R mutations revealed to date exhibit some variance in their clinical features as well, probably because different parts of the downstream IGF1R signalling cascade were affected. METHODS AND RESULTS: Using the recently developed technique multiplex ligation dependent probe amplification (MLPA), a chromosome 15q26.3 microdeletion harbouring part of the IGF1R gene was identified in a Dutch family. This deletion segregated with short height in seven out of 14 relatives across three generations. Metaphase fluorescence in situ hybridisation (FISH) and Affymetrix 250k single nucleotide polymorphism (SNP) microarray were used to characterise the deletion into more detail and showed that exons 11-21 of the IGF1R and a small hypothetical protein (LOC 145814) were deleted. CONCLUSION: Clinical work-up of this newly identified family, which constitutes the smallest (0.095 Mb) pure 15q26.3 interstitial deletion to date, confirms that disruption of the IGF1R gene does not induce major organ malformation or severe mental retardation.

Persistence of chronic myelocytic leukemia despite deletion of rearranged bcr/c-abl sequences in blast crisis.

We report on a Ph+ chronic myelocytic leukemia (CML) case, cytogenetically characterized by the occurrence of a second Philadelphia (Ph) chromosome in lymphoid blast crisis of T cell lineage. In situ hybridization analyses showed a deletion of translocated c-abl sequences, present on the Ph during chronic state, from both Ph in acute state. Moreover, Southern blot analyses of blastic cells exhibited a rearrangement within bcr, but a deletion of 5' bcr sequences, and Northern blots failed to detect the hybrid 8.5 kb bcr/c-abl transcript usually observed in Ph+ CML.

c-sis is translocated from chromosome 22 to chromosome 9 in chronic myelocytic leukemia.

By analysis of a series of somatic cell hybrids derived by fusion of either mouse or Chinese hamster cells with leukocytes from different chronic myelocytic leukemia (CML) patients or from normal donors, we have localized the human oncogene, c-sis, on the q11 to qter segment of chromosome 22 and demonstrated its translocation from chromosome 22 to chromosome 9 (q34) in CML.

Disturbed interaction of p21-rac with mutated p67-phox causes chronic granulomatous disease.

Chronic granulomatous disease (CGD) is characterized by the failure of phagocytic leukocytes to generate superoxide, needed for the intracellular killing of microorganisms. This is caused by mutations in any one of the four subunits of the nicotinamide adenine dinucleotide phosphate (NADPH) oxidase. In a rare, autosomal recessive form of CGD, a 67-kD cytosolic component of this enzyme (p67-phox) is missing. We here report on a patient with a mutation in the p67-phox gene that leads to expression of a nonfunctional p67-phox protein. The purified granulocytes of this patient failed to produce superoxide and contained about half of the normal amount of p67-phox. Analysis of the cDNA and genomic DNA of this patient showed that the patient is a compound heterozygote for a triplet nucleotide deletion in the p67-phox gene, predicting an in-frame deletion of lysine 58 in the p67-phox protein and a larger deletion of 11-13 kb in the other allele. Interestingly, the 58Lys deletion in p67-phox disrupts the interaction with p21-rac1, a ras-related protein involved in the activation of the NADPH oxidase. In contrast to normal neutrophils, in which p47-phox and p67-phox translocate to the plasma membrane upon cell activation, the cells of the patient did not show this translocation, indicating that an interaction between p67-phox and p21-rac1 is essential for translocation of these cytosolic proteins and activation of the NADPH oxidase. Moreover, this CGD patient represents the first case of disease caused by a disturbed binding of a ras-related protein to its target protein.

Pre- and postnatal diagnosis and outcome of fetuses and neonates with esophageal atresia and tracheoesophageal fistula.

OBJECTIVES: Clinical symptoms and ultrasound signs during pregnancy could suggest the presence of esophageal atresia (EA). However, most often EA is diagnosed postnatally. The aim of our study is to evaluate the course and outcome for prenatally and postnatally diagnosed EA. In addition, we studied the outcome of isolated versus nonisolated EA. METHODS: In a retrospective data analysis, ultrasound characteristics, maternal and neonatal variables as well as clinical outcome were compared for fetuses/neonates with prenatal (n = 30) or postnatal (n = 49) diagnosis of EA. Clinical outcome in terms of morbidity and mortality of isolated EA was compared with that of EA complicated by chromosomal or structural anomalies. RESULTS: Prenatally diagnosed children were born 2 weeks earlier than postnatally diagnosed children (36.4 weeks vs 38.2 weeks; P = 0.02). The former had higher mortality rates (30 vs 12%; P = 0.05) and more associated anomalies (80 vs 59%; P = 0.04). In both subsets, there was a high morbidity rate in the survivors (not significant). Nonisolated EA was associated with greater occurrence of polyhydramnios (53 vs 27%; P = 0.04) and higher mortality rate (28 vs 0%; P = 0.002). CONCLUSIONS: Mortality was significantly higher in prenatally diagnosed infants and in infants with additional congenital anomalies. Isolated EA is associated with good outcome.

Oncogenic GNAQ mutations are not correlated with disease-free survival in uveal melanoma.

BACKGROUND: Recently, oncogenic G protein alpha subunit q (GNAQ) mutations have been described in about 50% of uveal melanomas and in the blue nevi of the skin. METHODS: GNAQ exon 5 was amplified from 75 ciliary body and choroidal melanoma DNAs and sequenced directly. GNAQ mutation status was correlated with disease-free survival (DFS), as well as other clinical and histopathological factors, and with chromosomal variations detected by FISH and CGH. RESULTS: Of the 75 tumour DNA samples analysed, 40 (53.3%) harboured oncogenic mutations in GNAQ codon 209. Univariate and multivariate analysis showed that GNAQ mutation status was not significantly correlated with DFS. CONCLUSION: The GNAQ mutation status is not suitable to predict DFS. However, the high frequency of GNAQ mutations may render it a promising target for therapeutic intervention.

TACI mutations and disease susceptibility in patients with common variable immunodeficiency.

The most prevalent primary immunodeficiency is common variable immunodeficiency (CVID). Mutations have been described in four genes, ICOS, CD19, BAFF-R and TNFRSF13B (encoding TACI), together associated with 10-15% of CVID cases. We investigated a family with CVID and identified the heterozygous C104R TNFRSF13B mutation in two of the three index-children with CVID, a mother with selective immunoglobulin A deficiency, a mother with recurrent infections and a healthy grandfather. Remarkably, we did not find the TNFRSF13B mutation in the third index-child with CVID, despite his hypogammaglobulinaemia and decreased response to unconjugated pneumococcal vaccine. This family illustrates that TNFRSF13B mutations induce disease susceptibility rather than cause disease directly. Apparently, other genetic or environmental factors, still to be identified, contributed to the development of CVID in this family. Consequently, TNFRSF13B mutations must be interpreted with caution in the clinical setting.

Genetic and environmental factors in the etiology of esophageal atresia and/or tracheoesophageal fistula: an overview of the current concepts.

Esophageal atresia and/or tracheoesophageal fistula (EA/TEF) are severe congenital anomalies. Although recent years have brought significant improvement in clinical treatment, our understanding of the etiology of these defects is lagging. Many genes and genetic pathways have been implicated in the development of EA/TEF, but only a few genes have been shown to be involved in humans, in animals, or in both. Extrapolating data from animal models to humans is not always straightforward. Environmental factors may also carry a risk, but the mechanisms are yet to be elucidated. This review gives an overview of the current state of knowledge about both genetic and environmental risk factors in the etiology of EA/TEF.

Mapping of 5q35 chromosomal rearrangements within a genomically unstable region.

BACKGROUND: Recent molecular studies of breakpoints of recurrent chromosome rearrangements revealed the role of genomic architecture in their formation. In particular, segmental duplications representing blocks of >1 kb with >90% sequence homology were shown to mediate non-allelic homologous recombination (NAHR). However, the occurrence of the majority of newly detected submicroscopic imbalances cannot be explained by the presence of segmental duplications. Therefore, further studies are needed to investigate whether architectural features other than segmental duplications mediate these rearrangements. METHODS: We analysed a series of patients with breakpoints clustering within chromosome band 5q35. Using high density arrays and subsequent quantitative polymerase chain reaction (qPCR), we characterised the breakpoints of four interstitial deletions (including one associated with an unbalanced paracentric inversion), a duplication and a familial reciprocal t(5;18)(q35;q22) translocation. RESULTS AND CONCLUSION: Five of the breakpoints were located within an interval of approximately 265 kb encompassing the RANBP17 and TLX3 genes. This region is also targeted by the recurrent cryptic t(5;14)(q35;q32) translocation, which occurs in approximately 20% of childhood T cell acute lymphoblastic leukaemia (T-ALL). In silico analysis indicated the architectural features most likely to contribute to the genomic instability of this region, which was supported by our molecular data. Of further interest, in two patients and the familial translocation, the delineated breakpoint regions encompassed highly homologous LINEs (long interspersed nuclear elements), suggesting that NAHR between these LINEs may have mediated these rearrangements.

Can we improve outcome of congenital diaphragmatic hernia?

This review gives an overview of the disease spectrum of congenital diaphragmatic hernia (CDH). Etiological factors, prenatal predictors of survival, new treatment strategies and long-term morbidity are described. Early recognition of problems and improvement of treatment strategies in CDH patients may increase survival and prevent secondary morbidity. Multidisciplinary healthcare is necessary to improve healthcare for CDH patients. Absence of international therapy guidelines, lack of evidence of many therapeutic modalities and the relative low number of CDH patients calls for cooperation between centers with an expertise in the treatment of CDH patients. The international CDH Euro-Consortium is an example of such a collaborative network, which enhances exchange of knowledge, future research and development of treatment protocols.

Germline variant in MSX1 identified in a Dutch family with clustering of Barrett's esophagus and esophageal adenocarcinoma.

The vast majority of esophageal adenocarcinoma cases are sporadic and caused by somatic mutations. However, over the last decades several families have been identified with clustering of Barrett's esophagus and esophageal adenocarcinoma. This observation suggests that one or more hereditary factors may play a role in the initiation of Barrett's esophagus and esophageal adenocarcinoma in these families. A Dutch family with clustering of Barrett's esophagus and esophageal adenocarcinoma was identified. Normal DNA obtained from the proband diagnosed with Barrett's esophagus was analyzed with SNP array and exome sequencing. A custom-made panel consisting of potential germline variants was verified in the normal DNA of the affected family members. In addition, the respective tumors were analyzed for somatic loss of the wild type allele or the presence of an inactivating somatic mutation in the wild type allele. Exome sequencing revealed 244 candidate variants in the normal DNA of the proband, of which 212 variants were verified successfully. After the normal DNA of the affected family members was analyzed for the presence of the 212 potential germline variants and subsequently the respective tumors, only one potential germline variant in MSX1 (chr4: 4861985 T > G, c.359T > G, p.V120G, NM_002448) showed loss of the wild type allele in the tumor DNAs of the affected family members. A germline variant in MSX1 was identified in a Dutch family with clustering of Barrett's esophagus and esophageal adenocarcinoma. This finding indicates that the germline defect in MSX1 may be associated with Barrett's esophagus and cancer in this particular family.

Targeted disruption of the cell-cycle checkpoint gene ATR leads to early embryonic lethality in mice.

Checkpoints of DNA integrity are conserved throughout evolution, as are the kinases ATM (Ataxia Telangiectasia mutated) and ATR (Ataxia- and Rad-related), which are related to phosphatidylinositol (PI) 3-kinase [1] [2] [3]. The ATM gene is not essential, but mutations lead to ataxia telangiectasia (AT), a pleiotropic disorder characterised by radiation sensitivity and cellular checkpoint defects in response to ionising radiation [4] [5] [6]. The ATR gene has not been associated with human syndromes and, structurally, is more closely related to the canonical yeast checkpoint genes rad3(Sp) and MEC1(Sc) [7] [8]. ATR has been implicated in the response to ultraviolet (UV) radiation and blocks to DNA synthesis [8] [9] [10] [11], and may phosphorylate p53 [12] [13], suggesting that ATM and ATR may have similar and, perhaps, complementary roles in cell-cycle control after DNA damage. Here, we report that targeted inactivation of ATR in mice by disruption of the kinase domain leads to early embryonic lethality before embryonic day 8.5 (E8.5). Heterozygous mice were fertile and had no aberrant phenotype, despite a lower ATR mRNA level. No increase was observed in the sensitivity of ATR(+/-) embryonic stem (ES) cells to a variety of DNA-damaging agents. Attempts to target the remaining wild-type ATR allele in heterozygous ATR(+/-) ES cells failed, supporting the idea that loss of both alleles of the ATR gene, even at the ES-cell level, is lethal. Thus, in contrast to the closely related checkpoint gene ATM, ATR has an essential function in early mammalian development.

Translocation of oncogene c-sis from chromosome 22 to chromosome 11 in a Ewing sarcoma-derived cell line.

Somatic cell hybrids, obtained after fusion of translocation (11;22)-positive Ewing sarcoma cells and Chinese hamster fibroblasts, were assayed for the presence of immunoglobulin C lambda, Philadelphia chromosome breakpoint cluster region, and c-sis oncogene sequences. It was found that c-sis was translocated from chromosome 22 to chromosome 11 in the Ewing sarcoma cells used, indicating that the breakpoint must be proximal to this locus. Moreover, we found that the chromosome 22-linked C lambda and breakpoint cluster region sequences are not translocated. This result confirms an earlier cytogenetic observation that the Ewing sarcoma-associated breakpoint in chromosome 22 is distal to those observed in translocation (8;22)-positive Burkitt lymphoma and in Philadelphia chromosome-positive chronic myeloid leukemia.

The chronic myelocytic cell line K562 contains a breakpoint in bcr and produces a chimeric bcr/c-abl transcript.

In the DNAs of all Ph1-positive chronic myelocytic leukemia patients studied to date, a breakpoint on chromosome 22 (the Ph1 chromosome) can be demonstrated with a probe from the bcr (breakpoint cluster region). Although the K562 cell line was established from cells of a chronic myelocytic leukemia patient, we have been unable to detect the Ph1 chromosome by cytogenetic means. Employing a probe from the 5' region of bcr, we have cloned an amplified Ph1 breakpoint fragment from K562. This demonstrates that K562 contains multiple remnants of a Ph1 chromosome with a breakpoint within bcr and thus may serve as a model system for the study of Ph1-positive chronic myelocytic leukemia at a molecular level. The isolation of bcr cDNA sequences shows that parts of bcr encode a protein. Employing K562, we demonstrate the presence of an abnormally sized mRNA species hybridizing to c-abl and to a bcr cDNA probe, indicating the possible consequence of the Ph1 translocation on a transcriptional level in chronic myelocytic leukemia. The isolation and sequencing of a cDNA containing the breakpoint area of this mRNA provide further evidence for its chimeric structure. Cloning of large stretches of chromosomal DNA flanking bcr and c-abl sequences in K562 and identification of the exons participating in the formation of the chimeric mRNA shows that a splice of at least 99 kilobases is made to fuse the 3' bcr exon to the 5' c-abl exon. Furthermore two chimeric cDNAs were isolated containing chromosome 9 sequences that map 43.5 kilobases downstream from the K562 breakpoint. These chromosome 9 sequences neither hybridize to the 8.5-kilobase chimeric c-abl mRNA nor to normal c-abl mRNAs in Hela cells and probably represent incorrect splicing products present in the K562 cell line.

Breakpoint heterogeneity in t(10;11) translocation in AML-M4/M5 resulting in fusion of AF10 and MLL is resolved by fluorescent in situ hybridization analysis.

Ten AML-M4/M5 patients' samples containing a t(10;11) translocation, but with different cytogenetic breakpoints on chromosome 11q (11q13-23), were studied by G- and R-banding and fluorescent in situ hybridization. Southern blotting analysis, studied in five patients, revealed a rearranged MLL gene. Reverse transcription-PCR analysis carried out in six patients showed a 5' MLL-3' AF-10 fusion transcript. Fluorescent in situ hybridization studies suggested that in 8 of 10 patients, the rearrangement/fusion transcript resulted from an inversion of a part of 11q (q13q23) translocated to 10p12. In the other two patients, it is assumed that an inversion/translocation has occurred of a part of 10p to the der(11). The results suggest that the orientation of the AF-10 gene on 10p is 5' telomeric and 3' centromeric. This is the first example of opposite-oriented genes being involved in translocation to yield fusion transcripts.

Fusion of the homeobox gene HLXB9 and the ETV6 gene in infant acute myeloid leukemias with the t(7;12)(q36;p13).

Recently, we and others reported a recurrent t(7;12)(q36;p13) found in myeloid malignancies in children < or =18 months of age and associated with a poor prognosis. Fluorescence in situ hybridization studies mapped the 12p13 breakpoint to the first intron of ETV6 and narrowed down the region of 7q36 involved. By using the sequences made public recently by the Human Genome Project, two candidate genes in 7q36 were identified: the homeobox gene HLXB9 and c7orf3, a gene with unknown function. Reverse transcription-PCR of two cases with t(7;12), using primers for c7orf3 and ETV6, was negative. However, reverse transcription-PCR for HLXB9-ETV6 demonstrated alternative splicing; the two major bands corresponded to fusion of exon 1 of HLXB9 to exons 2 and 3, respectively, of ETV6. The reciprocal ETV6-HLXB9 transcript was not detected. It remains to be elucidated if the leukemic phenotype is attributable to the formation of the HLXB9-ETV6 fusion protein, which includes the helix-loop-helix and E26 transformation-specific DNA binding domains of ETV6 or to the disruption of the normal ETV6 protein.

High expression of ets-1 gene in human thymocytes and immature T leukemic cells.

The cellular ets-1 gene homologous to the 5' region of the v-ets sequence of the E26 retrovirus codes for a 6.8-kb mRNA that is translated into a 51-kDa protein in human cells. A survey of mRNA from human tissues showed the thymus as the tissue with the highest level of ets-1 transcription, within other hematopoietic organs and tissues, including spleen, fetal liver, lymph nodes, bone marrow, and peripheral lymphocytes exhibiting low or undetectable levels of hybridization. A high level of ets-1 expression was found in murine thymocyte mRNA as well. Investigation of the ets-1 expression levels in human leukemic samples showed that primary malignant T cells (T-ALLs), corresponding to intrathymic stages of maturation, have a much higher level of ets-1 mRNA than malignant T lymphoid cells with a mature phenotype, such as adult T cell leukemias (ATLs). T-ALLs were also higher in ets-1 expression than the other lymphoid (pre-T-ALL, c-ALL, pre-B-ALL) malignant cells analyzed. Insignificant amounts of the specific ets-1 mRNA were detected in several acute myeloid leukemias representing various degrees of maturation. The elevated ets-1 mRNA in thymocytes suggests a biological role for the ets-1 product in these cells that could be explored to investigate ets-1 function. Finally, the exhibited expression of ets-1 in lymphoid cells and absence from malignant myeloid cells makes it a candidate marker for phenotyping human hematopoietic tumors.

Development of an interphase fluorescent in situ hybridization (FISH) test to detect t(8;21) in AML patients.

The translocation (8;21) is a chromosome abnormality associated with acute myeloid leukemia (AML). As a consequence of the translocation the AML1 (CBFA2) gene in the 21q22 region is fused to the ETO(CDR,MTG8) gene in the 8q22 region, resulting in one transcriptionally active gene on the 8q- derivative chromosome. In this report we demonstrate the use of a highly specific dual-colour FISH method for the detection of t(8;21) on interphase cells. Genomic probes able to detect the chimeric AML1/ETO gene on the 8q- derivative chromosome were assayed on both normal and leukemic bone marrow and peripheral blood samples. Cut-off values were established by independent analysis of 15 bone marrow specimens negative for the translocation. The cut-off value of positive nuclei was determined to be 2% and the cut-off value for both positive nuclei and nuclei of uncertain classification, 4%. Persistence of cells above these cut-off values was interpreted as persistence of the mutated clone. A total of 36 samples at different disease stages were tested. Interphase cytogenetics detected the translocation at the onset and relapse in the BM or the PB of 14 AML patients with t(8;21). The technique appears to be an alternative tool to both conventional cytogenetics and reverse transcription polymerase chain reaction (RT-PCR) for the monitoring of disease during patients' follow-up. By enabling the analysis of individual cells, interphase FISH is ideal for clonality studies both for clinical and experimental applications.