Identification and interruption of inheritance of familial cryptic translocations: A case report.

BACKGROUND: Cryptic translocations can be identified via genetic analysis of aborted tissues or malformed infants, but it is difficult to deduce the parental origins of the translocations. In the absence of such information, it is not easy to distinguish translocations from normal embryos during pre-implantation genetic testing, that seeks to block familial transmission of translocations. METHODS: Here, we present a new method that detects cryptic translocations and blocks familial transmission thereof. Whole-genome, low-coverage mate-pair sequencing (WGLMPS) revealed chromosome breakpoint sequences, and preimplantation genetic haplotyping (PGH) was then used to discard embryos with cryptic translocations. RESULTS: Cryptic translocations were found in all four families, and familial transmission was successfully blocked in one family. CONCLUSION: Whole-genome, low-coverage mate-pair sequencing combined with preimplantation genetic haplotyping methods powerfully and practically identify cryptic translocations and block familial transmissions.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization.

PURPOSE: To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples. METHODS: MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers. RESULTS: In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences. CONCLUSION: Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.

Molecular characterization of TCF3::PBX1 chromosomal breakpoints in acute lymphoblastic leukemia and their use for measurable residual disease assessment.

The translocation t(1;19)(q23;p13) with the resulting chimeric TCF3::PBX1 gene is the third most prevalent recurrent chromosomal translocation in acute lymphoblastic leukemia and accounts for 3-5% of cases. The molecular background of this translocation has been incompletely studied, especially in adult cases. We characterized the chromosomal breakpoints of 49 patients with TCF3::PBX1 and the corresponding reciprocal PBX1::TCF3 breakpoints in 15 cases at the molecular level, thus providing an extensive molecular overview of this translocation in a well-defined study patient population. Breakpoints were found to be remarkably clustered not only in TCF3 but also in PBX1. No association with DNA repeats or putative cryptic recombination signal sequence sites was observed. A simplified detection method for breakpoint identification was developed and the feasibility of patient-specific chromosomal break sites as molecular markers for detecting measurable residual disease (MRD) was explored. A highly sensitive generic real-time PCR for MRD assessment using these breakpoint sequences was established that could serve as a useful alternative to the classical method utilizing rearranged immune gene loci. This study provides the first extensive molecular data set on the chromosomal breakpoints of the t(1;19)/TCF3::PBX1 aberration in adult ALL. Based on the obtained data a generic MRD method was developed that has several theoretical advantages, including an on average higher sensitivity and a greater stability of the molecular marker in the course of disease.

Gene sequencing and result analysis of balanced translocation carriers by third-generation gene sequencing technology.

Because the total gene copy number remains constant and all genes are normally expressed, carriers of balanced chromosomal translocations usually have a normal phenotype but are able to produce many different types of gametes during meiosis, and unbalanced gametes lead to increased risks of infertility, recurrent spontaneous abortion, stillbirth, neonatal death or malformations and intellectual abnormalities in offspring. The key to balanced translocations lies in finding the breakpoints, but current genetic testing techniques are all short-read sequencing, with the disadvantage of procedural complexity and imprecision for precisely identifying the breakpoints. The latest third-generation sequencing technology overcomes these drawbacks and uses robust long-read sequencing to accurately and rapidly detect genome-wide information and identify breakpoint locations. In this paper, we performed whole genome long-read sequencing using an Oxford Nanopore sequencer to detect the breakpoints of 4 balanced chromosomal translocation carriers. The results showed that employing about ~ 10× coverage confirmed 6 of the 8 breakpoints, of which, 2 had microdeletions/insertions identified near the breakpoints and 4 had breakpoints that disrupted the normal gene structure and were simultaneously tested for genome-wide structural variation (SV). The results show that whole genome long-read sequencing is an efficient method for pinpointing translocation breakpoints and providing genome-wide information, which is essential for medical genetics and preimplantation genetic testing.

HPV integration generates a cellular super-enhancer which functions as ecDNA to regulate genome-wide transcription.

Human papillomavirus (HPV) integration is a critical step in cervical cancer development; however, the oncogenic mechanism at the genome-wide transcriptional level is still poorly understood. In this study, we employed integrative analysis on multi-omics data of six HPV-positive and three HPV-negative cell lines. Through HPV integration detection, super-enhancer (SE) identification, SE-associated gene expression and extrachromosomal DNA (ecDNA) investigation, we aimed to explore the genome-wide transcriptional influence of HPV integration. We identified seven high-ranking cellular SEs generated by HPV integration in total (the HPV breakpoint-induced cellular SEs, BP-cSEs), leading to intra-chromosomal and inter-chromosomal regulation of chromosomal genes. The pathway analysis revealed that the dysregulated chromosomal genes were correlated to cancer-related pathways. Importantly, we demonstrated that BP-cSEs existed in the HPV-human hybrid ecDNAs, explaining the above transcriptional alterations. Our results suggest that HPV integration generates cellular SEs that function as ecDNA to regulate unconstrained transcription, expanding the tumorigenic mechanism of HPV integration and providing insights for developing new diagnostic and therapeutic strategies.

Analysis of 2 men with t(8;22)(q13;q13) and t(8;14)(q13;q22) chromosomal translocation karyotypes.

Male infertility is a multifactorial condition that is closely associated with chromosomal abnormalities. Reciprocal chromosomal translocation (RCT) is a significant structural genetic abnormality. The specific mechanisms of forms of RCT affecting male infertility include the product of chromosomally unbalanced gametes, thereby disrupting the structure and function of important genes responsible for spermatogenesis. RCT breakpoints have been found to disrupt gene structure and function in many medical fields However, the relationship between RCT breakpoints and male infertility remains to be determined. The purpose of this study is to describe 2 male carriers of RCTs 46,XY,t(8;22)(q13;q13) and 46,XY,t(8;14)(q13;q22). Both patients were collected from the second hospital of Jilin University. Semen parameters were detected using the computer-aided semen analysis system. Cytogenetic analysis was performed using standard operating procedure. Related genes on chromosomal breakpoints were searched using Online Mendelian Inheritance in Man. One man had semen parameters within the normal range, but the couple was infertile after 5 years of marriage. The other man showed normal semen parameters, and his wife had experienced 2 spontaneous miscarriages. Using a literature search, the association between chromosome 22q13 breakpoint and fertility were investigated. The results suggest that physicians should focus on the clinical phenotype of the patients and the breakpoints of RCT in genetic counseling. An important gene related to human male infertility is clearly located in chromosome region 22q13, and its function is worthy of further study.

Cryptic MYC insertions in Burkitt lymphoma: New data and a review of the literature.

The occurrence of MYC-negative Burkitt lymphoma (BL) has been discussed for many years. The real frequency of the MYC insertion in MYC-negative BL is still unknown. Fine-needle aspiration biopsies of 108 consecutive patients with clinicopathologically suspected BL (suspBL) were evaluated by flow cytometry, classical cytogenetics, and fluorescence in situ hybridization (FISH). We found 12 cases (11%) without the MYC rearrangement by FISH with a MYC breakapart probe: two patients (1.9%) with cryptic MYC/IGH fusion (finally diagnosed as BL) and 10 patients (9.3%) with 11q gain/loss (finally diagnosed as Burkitt-like lymphoma with 11q aberration). The exact breakpoints of the cryptic MYC/IGH were investigated by next-generation sequencing. The MYC insertions' breakpoints were identified in PVT1 in the first case, and 42 kb upstream of 5'MYC in the second case. To date, a molecular characterization of the MYC insertion in BL has only been reported in one case. Detailed descriptions of our MYC insertions in a routinely and consecutively diagnosed suspBL cohort will contribute to resolving the issue of MYC negativity in BL. In our opinion, the presence of the MYC insertions in BL and other lymphomas might be underestimated, because routine genetic diagnostics are usually based on FISH only, without karyotyping.

Posterior Amorphous Corneal Dystrophy: New Chromosomal Breakpoints in the Small Leucine-Rich Proteoglycan-Coding Region.

PURPOSE: The purpose of this study was to report the clinical features and describe the results obtained by multimodal corneal imaging of a patient with novel chromosomal breakpoints of the 12q21.33 locus. METHODS: This study was a case report and literature review. RESULTS: A 12-year-old girl presented with visual loss whose examination revealed a best-corrected visual acuity of 20/50 in her right eye and 20/35 in her left eye and corneal flattening and gray sheet-like opacities deep in the stroma. Anterior segment optical coherence tomography and ultrabiomicroscopy showed an evenly distributed hyperreflective line in the posterior stroma. Confocal microscopy revealed enlarged keratocytes and the presence of small reflective deposits from the pre-Descemet line to the endothelium. In addition, a 447-kb deletion that included the small leucine-rich proteoglycan-coding region in locus 12q21.33 was found. She was, therefore, diagnosed with PACD. CONCLUSIONS: PACD is a rare genetic disorder of the cornea characterized by gray sheet-like opacification of the posterior stroma in combination with corneal flattening. Confocal microscopy provides histologic segmentation of each corneal layer and shows the degree to which they are affected. New chromosomal breakpoints of a deletion in the small leucine-rich proteoglycan-coding region are hereby reported. PACD may be a contiguous gene syndrome, and further tests are required to identify the exact position responsible for the phenotypic variation.

Chromosome 10q-linked FSHD identifies DUX4 as principal disease gene.

BACKGROUND: Facioscapulohumeral dystrophy (FSHD) is an inherited muscular dystrophy clinically characterised by muscle weakness starting with the facial and upper extremity muscles. A disease model has been developed that postulates that failure in somatic repression of the transcription factor DUX4 embedded in the D4Z4 repeat on chromosome 4q causes FSHD. However, due to the position of the D4Z4 repeat close to the telomere and the complex genetic and epigenetic aetiology of FSHD, there is ongoing debate about the transcriptional deregulation of closely linked genes and their involvement in FSHD. METHOD: Detailed genetic characterisation and gene expression analysis of patients with clinically confirmed FSHD and control individuals. RESULTS: Identification of two FSHD families in which the disease is caused by repeat contraction and DUX4 expression from chromosome 10 due to a de novo D4Z4 repeat exchange between chromosomes 4 and 10. We show that the genetic lesion causal to FSHD in these families is physically separated from other candidate genes on chromosome 4. We demonstrate that muscle cell cultures from affected family members exhibit the characteristic molecular features of FSHD, including DUX4 and DUX4 target gene expression, without showing evidence for transcriptional deregulation of other chromosome 4-specific candidate genes. CONCLUSION: This study shows that in rare situations, FSHD can occur on chromosome 10 due to an interchromosomal rearrangement with the FSHD locus on chromosome 4q. These findings provide further evidence that DUX4 derepression is the dominant disease pathway for FSHD. Hence, therapeutic strategies should focus on DUX4 as the primary target.

A novel complex genomic rearrangement affecting the KCNJ2 regulatory region causes a variant of Cooks syndrome.

Cooks syndrome (CS) is an ultrarare limb malformation due to in tandem microduplications involving KCNJ2 and extending to the 5' regulatory element of SOX9. To date, six CS families were resolved at the molecular level. Subsequent studies explored the evolutionary and pathological complexities of the SOX9-KCNJ2/Sox9-Kcnj2 locus, and suggested a key role for the formation of novel topologically associating domain (TAD) by inter-TAD duplications in causing CS. Here, we report a unique case of CS associated with a de novo 1;17 translocation affecting the KCNJ2 locus. On chromosome 17, the breakpoint mapped between KCNJ16 and KCNJ2, and combined with a ~ 5 kb deletion in the 5' of KCNJ2. Based on available capture Hi-C data, the breakpoint on chromosome 17 separated KCNJ2 from a putative enhancer. Gene expression analysis demonstrated downregulation of KCNJ2 in both patient's blood cells and cultured skin fibroblasts. Our findings suggest that a complex rearrangement falling in the 5' of KCNJ2 may mimic the developmental consequences of in tandem duplications affecting the SOX9-KCNJ2/Sox9-Kcnj2 locus. This finding adds weight to the notion of an intricate role of gene regulatory regions and, presumably, the related three-dimensional chromatin structure in normal and abnormal human morphology.

Fertility problems in men carrying a translocation involved in breakpoints on chromosome 17p13: A retrospective, observational study.

Male infertility is a multifactorial reproductive disorder. The effect of genetic factors on male infertility has been the focus of research. Although a variety of genetic techniques are applied to male infertility in clinical practice, karyotype analysis remains a powerful and inexpensive technology. Reciprocal chromosomal translocation (RCT) is closely related to male infertility, but the clinical phenotypes of RCT carriers are varied, and the underlying pathological mechanism is unclear. Some studies suggest that RCT breakpoints disrupt the structure and function of important genes responsible for spermatogenesis. Several breakpoints of chromosome 17 are related to important genes, which can lead to spermatogenic failure. This study aimed to identify the clinical features of 3 men with translocation karyotypes involving breakpoints on chromosome 17p13. Semen analysis and cytogenetic analysis were performed with informed consent. Gene ontology analysis was performed for 60 pathogenic genes on chromosome band 17p13. Cytogenetic analysis showed that the karyotypes were 46, XY, t(6;17) (p21;p13), 46,XY,t(10;17)(q11.2;p13), and 46, XY, t(17;20) (p13;q13), respectively. Relevant studies and genes on breakpoints on chromosome 17p13 were searched for using PubMed. Fourteen reported cases of the same karyotype were reviewed. The results suggest that although chromosome 17 is closely related to spermatogenic failure, the clinical phenotypes of RCT carriers with involvement of 17p13 breakpoints are varied. The important genes involved in the breakpoint were analyzed. The results of molecular functions suggested that these targets genes on chromosome band 17p13 were mostly involved in microfilament motor activity, ATPase activity. These results suggested that the translocation chromosome and breakpoint analysis should be considered in the clinical assessment of the patients. Physicians should be aware of these in genetic counseling. These breakpoints and the function of related genes require further study.

Fine Breakpoint Mapping by Genome Sequencing Reveals the First Large X Inversion Disrupting the NHS Gene in a Patient with Syndromic Cataracts.

Inversions are structural variants that are generally balanced. However, they could lead to gene disruptions or have positional effects leading to diseases. Mutations in the NHS gene cause Nance-Horan syndrome, an X-linked disorder characterised by congenital cataracts and dental anomalies. Here, we aimed to characterise a balanced pericentric inversion X(p22q27), maternally inherited, in a child with syndromic bilateral cataracts by breakpoint mapping using whole-genome sequencing (WGS). 30× Illumina paired-end WGS was performed in the proband, and breakpoints were confirmed by Sanger sequencing. EdU assays and FISH analysis were used to assess skewed X-inactivation patterns. RNA expression of involved genes in the breakpoint boundaries was evaluated by droplet-digital PCR. We defined the breakpoint position of the inversion at Xp22.13, with a 15 bp deletion, disrupting the unusually large intron 1 of the canonical NHS isoform, and also perturbing topologically-associated domains (TADs). Moreover, a microhomology region of 5 bp was found on both sides. RNA analysis confirmed null and reduced NHS expression in the proband and his unaffected mother, respectively. In conclusion, we report the first chromosomal inversion disrupting NHS, fine-mapped by WGS. Our data expand the clinical spectrum and the pathogenic mechanisms underlying the NHS defects.

Marfan Syndrome Caused by Disruption of the FBN1 Gene due to A Reciprocal Chromosome Translocation.

Marfan syndrome (MFS) is a hereditary connective tissue disease caused by heterozygous mutations in the fibrillin-1 gene (FBN1) located on chromosome 15q21.1. A complex chromosomal rearrangement leading to MFS has only been reported in one case so far. We report on a mother and daughter with marfanoid habitus and no pathogenic variant in the FBN1 gene after next generation sequencing (NGS) analysis, both showing a cytogenetically reciprocal balanced translocation between chromosomes 2 and 15. By means of fluorescence in situ hybridization of Bacterial artificial chromosome (BAC) clones from the breakpoint area on chromosome 15 the breakpoint was narrowed down to a region of approximately 110 kb in FBN1. With the help of optical genome mapping (OGM), the translocation breakpoints were further refined on chromosomes 2 and 15. Sequencing of the regions affected by the translocation identified the breakpoint of chromosome 2 as well as the breakpoint of chromosome 15 in the FBN1 gene leading to its disruption. To our knowledge, this is the first report of patients with typical clinical features of MFS showing a cytogenetically reciprocal translocation involving the FBN1 gene. Our case highlights the importance of structural genome variants as an underlying cause of monogenic diseases and the useful clinical application of OGM in the elucidation of structural variants.

Adult acampomelic campomelic dysplasia and disorders of sex development due to a reciprocal translocation involving chromosome 17q24.3 upstream of the SOX9 gene.

Balanced chromosomal rearrangements with a breakpoint located upstream of the sex determining region Y-box 9 (SOX9) gene on chromosome 17q24.3 are associated with skeletal abnormalities, campomelic dysplasia (CMPD), or acampomelic campomelic dysplasia (ACMPD). We report on a female patient with a reciprocal translocation of t (11; 17) (p15.4; q24.3), who was diagnosed with acampomelic campomelic dysplasia. The 34-year-old Japanese patient presented with distinct skeletal abnormalities, profound intellectual disability, and female phenotype despite the presence of Y chromosome and the sex determining region Y (SRY) gene. Her menarche started at 33 years and 4 months after hormone therapy of estrogen therapy followed by estrogen progesterone therapy. By conducting whole genome sequencing followed by Sanger sequencing validation, we determined the precise breakpoint positions of the reciprocal translocation, one of which was located 203 kb upstream of the SOX9 gene. Considering the phenotypic variations previously reported among the CMPD/ACMPD patients with a chromosomal translocation in the vicinity of SOX9, the identified translocation was concluded to be responsible for all major phenotypes observed in the patient.

Determination of COL1A1-PDGFB breakpoints by next-generation sequencing in the molecular diagnosis of dermatofibrosarcoma protuberans.

OBJECTIVE: In most cases, dermatofibrosarcoma protuberans (DFSP) is characterized by the chromosomal translocation t (17; 22) (q22; q13) that leads to a fusion of collagen type 1 alpha 1 (COL1A1) and platelet-derived growth factor beta chain (PDGFB). Recently, next-generation sequencing (NGS) has been reported to detect fusion transcripts in some malignancies. Therefore, the present study aimed to evaluate the utility of the targeted NGS in detecting the COL1A1-PDGFB fusion in patients with DFSP. METHODS: We designed a targeted DNA capture panel to tile along the fusion regions, including exon, intron, and untranslated regions of the COL1A1 and PDGFB. A cohort of 18 DNA samples extracted from formalin-fixed, paraffin-embedded tissues was used to evaluate the targeted NGS. The results were compared with that of fluorescence in situ hybridization (FISH). RESULTS: The COL1A1-PDGFB fusion was identified in 13 of 18 cases (72.2%) by targeted NGS assay. PDGFB breakpoints were constantly found in exon 2, while breakpoints in COL1A1 varied from exon 15 to 46. Of these 18 cases assayed by FISH, 12 (66.7%) exhibited COL1A1-PDGFB fusion signals. One case (P9), which was FISH-negative, was demonstrated with the fusion by targeted NGS and validated by PCR and Sanger sequencing. The targeted NGS results showed a high concordance with the results of the FISH assay (94.4%). CONCLUSION: Our study reported a targeted NGS assay for detecting the breakpoints of the COL1A1-PDGFB fusion gene, which can be implemented in diagnosing patients with DFSP.

R-loops and regulatory changes in chronologically ageing fission yeast cells drive non-random patterns of genome rearrangements.

Aberrant repair of DNA double-strand breaks can recombine distant chromosomal breakpoints. Chromosomal rearrangements compromise genome function and are a hallmark of ageing. Rearrangements are challenging to detect in non-dividing cell populations, because they reflect individually rare, heterogeneous events. The genomic distribution of de novo rearrangements in non-dividing cells, and their dynamics during ageing, remain therefore poorly characterized. Studies of genomic instability during ageing have focussed on mitochondrial DNA, small genetic variants, or proliferating cells. To characterize genome rearrangements during cellular ageing in non-dividing cells, we interrogated a single diagnostic measure, DNA breakpoint junctions, using Schizosaccharomyces pombe as a model system. Aberrant DNA junctions that accumulated with age were associated with microhomology sequences and R-loops. Global hotspots for age-associated breakpoint formation were evident near telomeric genes and linked to remote breakpoints elsewhere in the genome, including the mitochondrial chromosome. Formation of breakpoint junctions at global hotspots was inhibited by the Sir2 histone deacetylase and might be triggered by an age-dependent de-repression of chromatin silencing. An unexpected mechanism of genomic instability may cause more local hotspots: age-associated reduction in an RNA-binding protein triggering R-loops at target loci. This result suggests that biological processes other than transcription or replication can drive genome rearrangements. Notably, we detected similar signatures of genome rearrangements that accumulated in old brain cells of humans. These findings provide insights into the unique patterns and possible mechanisms of genome rearrangements in non-dividing cells, which can be promoted by ageing-related changes in gene-regulatory proteins.

Genomic Chaos (Multiple Copy Number Variations and Structural Reorganization) Detected in Two Prenatal Cases.

The use of new technologies in the routine diagnosis of constitutional abnormalities, such as high-resolution chromosomal microarray and next-generation sequencing, has unmasked new mechanisms for generating structural variation of the human genome. For example, complex chromosome rearrangements can originate by a chromosome catastrophe phenomenon in which numerous genomic rearrangements are apparently acquired in a single catastrophic event. This phenomenon is named chromoanagenesis (from the Greek "chromo" for chromosome and "anagenesis" for rebirth). Herein, we report 2 cases of genomic chaos detected at prenatal diagnosis. The terms "chromothripsis" and "chromoanasynthesis" and the challenge of genetic counseling are discussed.

GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing.

GRIDSS2 is the first structural variant caller to explicitly report single breakends-breakpoints in which only one side can be unambiguously determined. By treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints, GRIDSS2 can explain 47% of somatic centromere copy number changes using single breakends to non-centromere sequence. On a cohort of 3782 deeply sequenced metastatic cancers, GRIDSS2 achieves an unprecedented 3.1% false negative rate and 3.3% false discovery rate and identifies a novel 32-100 bp duplication signature. GRIDSS2 simplifies complex rearrangement interpretation through phasing of structural variants with 16% of somatic calls phasable using paired-end sequencing.

Breakpoint characterization of a rare alpha0 -thalassemia deletion using targeted locus amplification on genomic DNA.

INTRODUCTION: The high-sequence homology of the α-globin-gene cluster is responsible for microhomology-mediated recombination events during meiosis, resulting in a high density of deletion breakpoints within a 10 kb region. Commonly used deletion detection methods, such as multiplex ligation-dependent probe amplification (MLPA) and Southern blot, cannot exactly define the breakpoints. This typically requires long-range PCR, which is not always successful. Targeted locus amplification (TLA) is a targeted enrichment method that can be used to sequence up to 70 kb of neighboring DNA sequences without prior knowledge about the target site. METHODS: Genomic DNA (gDNA) TLA is a technique that folds isolated DNA, ensuring that adjacent loci are in a close spatial proximity. Subsequent digestion and religation form DNA circles that are amplified using fragment-specific inverse primers, creating a library that is suitable for Illumina sequencing. RESULTS: Here, we describe the characterization of a rare 16 771 bp deletion, utilizing gDNA TLA with a single inverse PCR primer set on one end of the breakpoint. Primers for breakpoint PCR were designed to confirm the deletion breakpoints and were consequently used to characterize the same deletion in 10 additional carriers sharing comparable hematologic data and similar MLPA results. CONCLUSIONS: The gDNA TLA technology was successfully used to identify deletion breakpoints within the alpha-globin cluster. The deletion was described only once in an earlier study as the --gb , but as it was not registered correctly in the available databases, it was not initially recognized as such.