Use of single nucleotide polymorphism-based mapping arrays to detect copy number changes and loss of heterozygosity in multiple myeloma

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Chromothripsis 18 in multiple myeloma patient with rapid extramedullary relapse

Background

Catastrophic chromosomal event known as chromothripsis was proven to be a significant hallmark of poor prognosis in several cancer diseases. While this phenomenon is very rare in among multiple myeloma (MM) patients, its presence in karyotype is associated with very poor prognosis.

Case presentation

In our case, we report a 62 year female patient with rapid progression of multiple myeloma (MM) into extramedullary disease and short overall survival (OS = 23 months). I-FISH investigation revealed presence of gain 1q21 and hyperdiploidy (+ 5,+ 9,+ 15) in 82% and 86%, respectively, while IgH rearrangements, del(17)(p13) and del(13)(q14) were evaluated as negative.Whole-genome profiling using array-CGH showed complex genomic changes including hyperdiploidy (+ 3,+ 5,+ 9,+ 11, + 15,+ 19), monosomy X, structural gains (1q21-1q23.1, 1q32-1q44, 16p13.13-16p11.2) and losses (1q23.1-1q32.1; 8p23.3-8p11.21) of genetic material and chromothripsis in chromosome 18 with 6 breakpoint areas. Next-generation sequencing showed a total of 338 variants with 1.8% (6/338) of pathological mutations in NRAS (c.181C > A; p.Gln61Lys) or variants of unknown significance in TP53, CUX1 and POU4F1.

Conclusions

Our findings suggest that presence of chromothripsis should be considered as another important genetic hallmark of poor prognosis in MM patients and utilization of genome-wide screening techniques such as array-CGH and NGS improves the clinical diagnostics of the disease.

Whole Exome Sequencing in Multiple Myeloma to Identify Somatic Single Nucleotide Variants and Key Translocations Involving Immunoglobulin Loci and MYC

Multiple myeloma is a malignancy of terminally differentiated plasma cells in the bone marrow. These plasma cells produce high levels of immunoglobulin which cause end-organ damage. Rearrangements within the immunoglobulin loci are a physiological part of B cell development, but these DNA level double-strand breaks may result in interchromosomal translocations. There are five main translocations involving the Ig loci: t(4;14) 12%, t(6;14) 1%, t(11;14) 15%, t(14;16) 3%, and t(14;20) 2%. These are primary events, found in all cells within the tumor clone and are associated with different prognosis. The t(4;14), t(14;16), and t(14;20) are associated with a poor prognosis, whereas the others are associated with a more favorable prognosis. Rearrangements at the MYC locus are also associated with a poor prognosis and increased expression of MYC. MYC rearrangements are frequent (25%) and involve interchromosomal translocations involving Ig loci or other partners, but also include intrachromosomal inversions, duplications and deletions. As such, the Ig and MYC loci are key players in the myeloma genome and including these in any genomic studies is key to understanding the relationship with other abnormalities. We have designed a custom capture of the Ig and MYC loci which can be added to exome or targeted captures to inform on these key events. This saves on performing additional tests to determine these events, which are generally mandatory for any genetic investigations in myeloma. This custom capture is also relevant to other B cell malignancies where MYC and Ig translocations occur.

Genome-wide screening of cytogenetic abnormalities in multiple myeloma patients using array-CGH technique: a Czech multicenter experience

Characteristic recurrent copy number aberrations (CNAs) play a key role in multiple myeloma (MM) pathogenesis and have important prognostic significance for MM patients. Array-based comparative genomic hybridization (aCGH) provides a powerful tool for genome-wide classification of CNAs and thus should be implemented into MM routine diagnostics. We demonstrate the possibility of effective utilization of oligonucleotide-based aCGH in 91 MM patients. Chromosomal aberrations associated with effect on the prognosis of MM were initially evaluated by I-FISH and were found in 93.4% (85/91). Incidence of hyperdiploidy was 49.5% (45/91); del(13)(q14) was detected in 57.1% (52/91); gain(1)(q21) occurred in 58.2% (53/91); del(17)(p13) was observed in 15.4% (14/91); and t(4;14)(p16;q32) was found in 18.6% (16/86). Genome-wide screening using Agilent 44K aCGH microarrays revealed copy number alterations in 100% (91/91). Most common deletions were found at 13q (58.9%), 1p (39.6%), and 8p (31.1%), whereas gain of whole 1q was the most often duplicated region (50.6%). Furthermore, frequent homozygous deletions of genes playing important role in myeloma biology such as TRAF3, BIRC1/BIRC2, RB1, or CDKN2C were observed. Taken together, we demonstrated the utilization of aCGH technique in clinical diagnostics as powerful tool for identification of unbalanced genomic abnormalities with prognostic significance for MM patients.

DNA repair pathways in human multiple myeloma: role in oncogenesis and potential targets for treatment

Every day, cells are faced with thousands of DNA lesions, which have to be repaired to preserve cell survival and function. DNA repair is more or less accurate and could result in genomic instability and cancer. We review here the current knowledge of the links between molecular features, treatment, and DNA repair in multiple myeloma (MM), a disease characterized by the accumulation of malignant plasma cells producing a monoclonal immunoglobulin. Genetic instability and abnormalities are two hallmarks of MM cells and aberrant DNA repair pathways are involved in disease onset, primary translocations in MM cells, and MM progression. Two major drugs currently used to treat MM, the alkylating agent Melphalan and the proteasome inhibitor Bortezomib act directly on DNA repair pathways, which are involved in response to treatment and resistance. A better knowledge of DNA repair pathways in MM could help to target them, thus improving disease treatment.

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

SNP-based mapping arrays reveal high genomic complexity in monoclonal gammopathies, from MGUS to myeloma status

Genetic events mediating transformation from premalignant monoclonal gammopathies (MG) to multiple myeloma (MM) are unknown. To obtain a comprehensive genomic profile of MG from the early to late stages, we performed high-resolution analysis of purified plasma cells from 20 MGUS, 20 smoldering MM (SMM) and 34 MM by high-density 6.0 SNP array. A progressive increase in the incidence of copy number abnormalities (CNA) from MGUS to SMM and to MM (median 5, 7.5 and 12 per case, respectively) was observed (P=0.006). Gains on 1q, 3p, 6p, 9p, 11q, 19p, 19q and 21q along with 1p, 16q and 22q deletions were significantly less frequent in MGUS than in MM. Although 11q and 21q gains together with 16q and 22q deletions were apparently exclusive of MM status, we observed that these abnormalities were also present in minor subclones in MGUS. Overall, a total of 65 copy number-neutral LOH (CNN-LOH) were detected. Their frequency was higher in active MM than in the asymptomatic entities (P=0.047). A strong association between genetic lesions and fragile sites was also detected. In summary, our study shows an increasing genomic complexity from MGUS to MM and identifies new chromosomal regions involved in CNA and CNN-LOH.

Unique genetic profile of sporadic colorectal cancer liver metastasis versus primary tumors as defined by high-density single-nucleotide polymorphism arrays

Most genetic studies in colorectal carcinomas have focused on those abnormalities that are acquired by primary tumors, particularly in the transition from adenoma to carcinoma, whereas few studies have compared the genetic abnormalities of primary versus paired metastatic samples. In this study, we used high-density 500K single-nucleotide polymorphism arrays to map the overall genetic changes present in liver metastases (n=20) from untreated colorectal carcinoma patients studied at diagnosis versus their paired primary tumors (n=20). MLH1, MSH2 and MSH6 gene expression was measured in parallel by immunohistochemistry. Overall, metastatic tumors systematically contained those genetic abnormalities observed in the primary tumor sample from the same subject. However, liver metastases from many cases (up to 8 out of 20) showed acquisition of genetic aberrations that were not found in their paired primary tumors. These new metastatic aberrations mainly consisted of (1) an increased frequency of genetic lesions of chromosomes that have been associated with metastatic colorectal carcinoma (1p, 7p, 8q, 13q, 17p, 18q, 20q) and, more interestingly, (2) acquisition of new chromosomal abnormalities (eg, losses of chromosomes 4 and 10q and gains of chromosomes 5p and 6p). These genetic changes acquired by metastatic tumors may be associated with either the metastatic process and/or adaption of metastatic cells to the liver microenvironment. Further studies in larger series of patients are necessary to dissect the specific role of each of the altered genes and chromosomal regions in the metastatic spread of colorectal tumors.

Association between genetic subgroups of pancreatic ductal adenocarcinoma defined by high density 500 K SNP-arrays and tumor histopathology

The specific genes and genetic pathways associated with pancreatic ductal adenocarcinoma are still largely unknown partially due to the low resolution of the techniques applied so far to their study. Here we used high-density 500 K single nucleotide polymorphism (SNP)-arrays to define those chromosomal regions which most commonly harbour copy number (CN) alterations and loss of heterozygozity (LOH) in a series of 20 PDAC tumors and we correlated the corresponding genetic profiles with the most relevant clinical and histopathological features of the disease. Overall our results showed that primary PDAC frequently display (>70%) extensive gains of chromosomes 1q, 7q, 8q and 20q, together with losses of chromosomes 1p, 9p, 12q, 17p and 18q, such chromosomal regions harboring multiple cancer- and PDAC-associated genes. Interestingly, these alterations clustered into two distinct genetic profiles characterized by gains of the 2q14.2, 3q22.1, 5q32, 10q26.13, 10q26.3, 11q13.1, 11q13.3, 11q13.4, 16q24.1, 16q24.3, 22q13.1, 22q13.31 and 22q13.32 chromosomal regions (group 1; n = 9) versus gains at 1q21.1 and losses of the 1p36.11, 6q25.2, 9p22.1, 9p24.3, 17p13.3 and Xp22.33 chromosomal regions (group 2; n = 11). From the clinical and histopathological point of view, group 1 cases were associated with smaller and well/moderately-differentiated grade I/II PDAC tumors, whereas and group 2 PDAC displayed a larger size and they mainly consisted of poorly-differentiated grade III carcinomas. These findings confirm the cytogenetic complexity and heterogenity of PDAC and provide evidence for the association between tumor cytogenetics and its histopathological features. In addition, we also show that the altered regions identified harbor multiple cancer associate genes that deserve further investigation to determine their relevance in the pathogenesis of PDAC.

Mapping of genetic abnormalities of primary tumours from metastatic CRC by high-resolution SNP arrays

Background

For years, the genetics of metastatic colorectal cancer (CRC) have been studied using a variety of techniques. However, most of the approaches employed so far have a relatively limited resolution which hampers detailed characterization of the common recurrent chromosomal breakpoints as well as the identification of small regions carrying genetic changes and the genes involved in them.

Methodology/Principal Findings

Here we applied 500K SNP arrays to map the most common chromosomal lesions present at diagnosis in a series of 23 primary tumours from sporadic CRC patients who had developed liver metastasis. Overall our results confirm that the genetic profile of metastatic CRC is defined by imbalanced gains of chromosomes 7, 8q, 11q, 13q, 20q and X together with losses of the 1p, 8p, 17p and 18q chromosome regions. In addition, SNP-array studies allowed the identification of small (<1.3 Mb) and extensive/large (>1.5 Mb) altered DNA sequences, many of which contain cancer genes known to be involved in CRC and the metastatic process. Detailed characterization of the breakpoint regions for the altered chromosomes showed four recurrent breakpoints at chromosomes 1p12, 8p12, 17p11.2 and 20p12.1; interestingly, the most frequently observed recurrent chromosomal breakpoint was localized at 17p11.2 and systematically targeted the FAM27L gene, whose role in CRC deserves further investigations.

Conclusions/Significance

In summary, in the present study we provide a detailed map of the genetic abnormalities of primary tumours from metastatic CRC patients, which confirm and extend on previous observations as regards the identification of genes potentially involved in development of CRC and the metastatic process.

Chediak-Higashi syndrome with early developmental delay resulting from paternal heterodisomy of chromosome 1

Chediak-Higashi syndrome (CHS) is a rare autosomal recessive disease characterized by variable oculocutaneous albinism, immunodeficiency, mild bleeding diathesis, and an accelerated lymphoproliferative state. Abnormal lysosome-related organelle membrane function leads to the accumulation of large intracellular vesicles in several cell types, including granulocytes, melanocytes, and platelets. This report describes a severe case of CHS resulting from paternal heterodisomy of chromosome 1, causing homozygosity for the most distal nonsense mutation (p.E3668X, exon 50) reported to date in the LYST/CHS1 gene. The mutation is located in the WD40 region of the CHS1 protein. The patient's fibroblasts expressed no detectable CHS1. Besides manifesting the classical CHS findings, the patient exhibited hypotonia and global developmental delays, raising concerns about other effects of heterodisomy. An interstitial 747 kb duplication on 6q14.2-6q14.3 was identified in the propositus and paternal samples by comparative genomic hybridization. SNP genotyping revealed no additional whole chromosome or segmental isodisomic regions or other dosage variations near the crossover breakpoints on chromosome 1. Unmasking of a separate autosomal recessive cause of developmental delay, or an additive effect of the paternal heterodisomy, could underlie the severity of the phenotype in this patient.

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome

PURPOSE

Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes.

EXPERIMENTAL DESIGN

Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome.

RESULTS

We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets.

CONCLUSIONS

We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.

A SNP microarray and FISH-based procedure to detect allelic imbalances in multiple myeloma: an integrated genomics approach reveals a wide gene dosage effect

Multiple myeloma (MM) is characterized by marked genomic heterogeneity. Beyond structural rearrangements, a relevant role in its biology is represented by allelic imbalances leading to significant variations in ploidy status. To elucidate better the genomic complexity of MM, we analyzed a panel of 45 patients using combined FISH and microarray approaches. We firstly generated genome-wide profiles of 41 MMs and four plasma cell leukemias, using a self-developed procedure to infer exact local copy numbers (CNs) for each sample. Our analysis allowed the identification of a significant fraction of patients showing near-tetraploidy. Furthermore, a conventional hierarchical clustering analysis showed that near-tetraploidy, 1q gain, hyperdiploidy, and recursive deletions at 1p and chromosomes 13, 14, and 22 were the main aberrations driving samples grouping. Moreover, mapping information was integrated with gene expression profiles of the tumor samples. A multiclass analysis of transcriptional profiles characterizing the different clusters showed marked gene-dosage effects, particularly concerning 1q transcripts; this finding was also confirmed by a nonparametric analysis between normalized gene expression levels and local CN variations (1027 highly-significant correlated genes). Finally, we identified several loci in which gene expression correlated with the occurrence of loss of heterozygosity. Our results provide insights into the composite network linking genome structure and transcriptional features in MM.

Genome-wide Mapping of Copy Number Variations Using SNP Arrays

The availability of high-density single nucleotide polymorphism (SNP) microarrays in recent years has proven to be a great step forward in the context of global analysis of genomic abnormalities in disease. SNP arrays offer great robustness, high resolution and the possibility to detect a variety of different genomic copy number variations such as submicroscopic deletions, amplifications, loss of heterozygosity and uniparental disomy. Moreover, they can be used to perform genome wide association studies. Therefore, SNP arrays harbor several advancements over traditional molecular methods to analyze genomic aberrations, such as cytogenetic analyses, fluorescence in situ hybridization or comparative genomic hybridization methods. Until now, SNP arrays have exclusively been used in experimental research and have enabled seminal new discoveries in many fields by identifying common genomic lesions underlying specific diseases, especially cancer. However, it is foreseeable that SNP arrays will also take up a position in routine diagnostic processes in the future. This review focuses on technical principles of the SNP array technology and their utilization to detect submicroscopic genomic and polymorphic markers associated with disease.