Characterization Of Copy Number Variations Research Articles

Abstract 1764: Long read sequencing allows comprehensive molecular profiling of complex karyotype acute myeloid leukemia (CKAML) with 5q deletions

Abstract Background Acute myeloid leukemia with complex karyotype (CKAML) does frequently show del(5q), deletion of chromosome 5q. Previous work in MDS and AML showed that CKAML form clusters distinguishing del(5q) cases with commonly retained regions (CRRs) in the telomeric (≥5q34) and centromeric (≤5q14.2) ends from cases with loss of the telomeric and/or centromeric regions. Interestingly, del(5q) patients with CRRs were shown to have less genomic lesions and were associated with a better prognosis than other 5q-deleted CKAML cases. Aim: Combining the advantages of two different long read sequencing technologies, Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), we aimed to more comprehensively characterize the molecular landscape of CKAML with del(5q). Our integrative workflow enabled an in-depth characterization of copy number variations (CNV), genomic rearrangements, transcriptomic and epigenomic changes. Methods: Low coverage whole genome sequencing (WGS) using ONT was performed on 209 patients with CKAML to identify del(5q) cases for further in-depth analysis. For 10 selected cases with known CRR status, genomic deep sequencing using PacBio Revio and ONT PromethION was performed, as well as PacBio Revio MAS-seq for fusion transcript detection. Results: A mean sequencing depth of 3-fold (range 0.4-17.3-fold) was obtained by low coverage WGS of the 209 CKAML cases. With a resolution of 0.1 Mbp n=134 del(5q) cases were identified in our CKAML cohort (64%) and commonly deleted regions on 5q could be resolved at 5q21.2 to 5q23.3 and 5q33.3 to 5q35 (Fisher test, FDR Benjamini-Hochberg, p<10−30). These regions separated del(5q) in n=52 cases showing CRR retention and n=82 showing partially no CRR, (either the centromeric or telomeric regions, n=74, or of both regions, n=8 involved). Subsequent co-occurrence analysis provided evidence of an association between the loss of telomeric 5q region with a 3p12-25 deletion. In-depth analysis (30-50 fold) of ten distinct del(5q) cases by PacBio and ONT revealed additional insights into the complex clonal composition of del(5q) cases. For example, one case showed two distinct subclones with 5q14.2-5q35.1 deletion and with translocation t(3;5)(p14;q14) respectively. Moreover, we could demonstrate that 5q-aberrations are far more complex than initially thought showing, e.g. inverted terminal duplications that resulted in fusion transcripts involving PDE4D and leading to its deregulated expression. Conclusion: Using our long read deep sequencing approach we could further delineate the structural complexity of del(5q) CKAML. While gene expression profiling, methylation analysis and correlation of findings with clinical outcome are ongoing and will be presented at the meeting, our results help to further pinpoint the biology underlying del(5q) and to further refine this heterogeneous CKAML subgroup. Citation Format: Felix Schick, Ekaterina Slonova, Eric Straeng, Sarah Huet, Maurizio Raso, Marlon Tilgner, Frederik Damm, Pierre Sujobert, Olga Blau, Lars Bullinger, Anna Dolnik. Long read sequencing allows comprehensive molecular profiling of complex karyotype acute myeloid leukemia (CKAML) with 5q deletions [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1764.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations.

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

Characterization of copy-number variants in a large cohort of patients with von Willebrand disease reveals a relationship between disrupted regions and disease type

BackgroundGenetic analysis for von Willebrand disease (VWD) commonly utilizes DNA sequencing to identify variants in the von Willebrand factor (VWF) gene; however, this technique cannot always detect copy-number variants (CNVs). Additional mapping of CNVs in patients with VWD is needed. ObjectivesThis study aimed to characterize CNVs in a large sample of VWF mutation-negative VWD patients. MethodsTo determine the role of CNVs in VWD, a VWF high-resolution comparative genomic hybridization array was custom-designed to avoid multiple sequence variations, repeated sequences, and the VWF pseudogene. This was performed on 204 mutation-negative subjects for whom clinical variables were also available. ResultsAmong the 204 patients, 7 unique CNVs were found, with a total of 24 CNVs (12%). Of the 7 unique CNVs, 1 was novel, 1 was found in a VWF database, and 5 were previously reported. All patients with type 1C VWD and a CNV had the same exon 33 and 34 in-frame deletion. Certain clinical variables were also significantly different between those with and without CNVs. ConclusionThe in-frame deletion in patients with type 1C VWD exactly matches the D4N module of the D4 domain, a region where mutations and deletions are known to affect clearance. We observed significantly higher VWF–to–ristocetin cofactor levels in patients with type 1C VWD and a CNV than in patients without a CNV, suggesting a relationship between CNVs and the increased clearance observed in patients with type 1C VWD. Glycoprotein IbM activity was significantly lower in patients with type 1 VWD and a CNV than in patients without a CNV, suggesting that platelet binding is more affected by CNVs than single base pair mutations. This work elucidates some of the underlying genetic mechanisms of CNVs in these patients.

Genome-Wide Discovery of Structural Variants Reveals Distinct Variant Dynamics for Two Closely Related Monilinia Species.

Structural variants (SVs) are variants with sizes bigger than 50 bp and capable of changing the size, copy number, location, orientation, and sequence content of genomic DNA. Although these variants have been proven to be extensive and involved in many evolutionary processes along the tree of life, there is still insufficient information on many fungal plant pathogens. In this study, the extent of SVs, as well as single-nucleotide polymorphisms (SNPs), has been determined for two prominent species of the Monilinia genus (the causal agents of brown rot disease in pome and stone fruits): Monilinia fructicola and Monilinia laxa for the first time. The genomes of M. fructicola were found to be more variant-rich in contrast to M. laxa based on the reference-based variant calling (with a total number of 266.618 and 190.599 SNPs and 1,540 and 918 SVs, respectively). The extent, as well as distribution of SVs, presented high conservation within the species and high diversity between the species. Investigation of potential functional effects of characterized variants revealed high potential relevance of SVs. Moreover, the detailed characterization of copy number variations (CNVs) for each isolate revealed that around 0.67% of M. fructicola genomes and 2.06% of M. laxa genomes are copy number variables. The variant catalog as well as distinct variant dynamics within and between the species presented in this study opens doors for many further research questions.

Copy Number Variations as Determinants of Colorectal Tumor Progression in Liquid Biopsies.

Over the years, increasing evidence has shown that copy number variations (CNVs) play an important role in the pathogenesis and prognosis of Colorectal Cancer (CRC). Colorectal adenomas are highly prevalent lesions, but only 5% of these adenomas ever progress to carcinoma. This review summarizes the different CNVs associated with adenoma-carcinoma CRC progression and with CRC staging. Characterization of CNVs in circulating free-RNA and in blood-derived exosomes augers well with the potential of using such assays for patient management and early detection of metastasis. To overcome the limitations related to tissue biopsies and tumor heterogeneity, using CNVs to characterize tumor-derived materials in biofluids provides less invasive sampling methods and a sample that collectively represents multiple tumor sites in heterogeneous samples. Liquid biopsies provide a source of circulating tumor DNA (ctDNA), circulating tumor cells (CTCs), tumor-derived exosomes (TDE), circulating free RNA, and non-coding RNA. This review provides an overview of the current diagnostic and predictive models from liquid biopsies.

Characterization of CNVs in a Large Cohort of VWD Patients Reveals Relationship between Disrupted Regions, Clinical Variables and Disease Type

BACKGROUND: Copy-number variants (CNVs) have been associated to all VWD types. Although many of the underlying pathological mechanisms have not been investigated, in-frame deletions have been shown to cause defects in VWF biosynthesis and secretion. While there are previous studies describing CNVs in VWD patients, these only discuss a small number of CNVs each, and are in varying populations, making it difficult to compare between populations as the differences could be population specific. METHODS: A VWF high-resolution aCGH array was custom designed using long oligonucleotide probes ~60 bp in length to avoid multiple sequence variations, repeated sequences and the pseudogene. Data were manually inspected and analyzed with CytoSure Interpret Software. A multiplex PCR approach was used to map deletion breakpoints and to confirm the presence and zygosity of the deletions in the index case as well as family members. CNVs were mapped and t-tests performed on clinical laboratory values within each VWD type for those with and without a CNV. aCGH was performed in two datasets on low-VWF and VWD index cases with no known pathogenic VWF sequence variant for a total of 205 patients. RESULTS: Of the 205 patients, CNVs were found in 25 (12%). No CNVs were found in low-VWF patients and the only duplication occurred in a Type 2 (n=1) patient. The remaining CNVs were deletions in VWD Type 1+1S (n=12), 1C (n=5), and 3 (n=7) patients. All Type 1C patients (n=5) had the same exon 33-34 deletion. All but 1 of the Type 1+1S patients (n=11) had the same exon 4-5 deletion, while the Type 3 patients (n=7) had deletions spread across the gene. Further, within the Type 1+1S patients, VWFpp/VWF:Ag ratio was significantly lower in those with a CNV than those without (two-tailed t-test, p=0.02). Lastly, ristocetin cofactor was significantly higher in the Type 1C patients with a CNV versus those without a CNV (two-tailed t-test, p=8.2x10-5). DISCUSSION: Here we have analyzed a large dataset of CNVs defined by aCGH in VWD patients with no known pathogenic VWF sequence variant. The location of the deletion in all Type 1C patients exactly matches the D4N module of the D4 domain of the mature VWF protein, a region where mutations and deletions are known to affect clearance. The D4 domain is the only D domain to also include an N module. Although previous work has shown that deletions of exons 33-34 potentially cause a mild reduction in VWF:Ag levels and only a moderate impact on secretion, we observed significantly higher VWF:RCo levels in those Type 1C patients with a CNV when compared to those without. As expected, Type 1C patients had a significantly higher VWFpp/VWF:Ag ratio than Type 1+1S patients overall (p=2.11x10-15). However, the Type 1+1S patient with a CNV had significantly lower VWFpp/VWF:Ag ratio than those without a CNV. This work elucidates some of the underlying genetic mechanisms of CNVs in these patients. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Identification of copy number variations and candidate genes for reproduction traits in Finnish pig populations

Animal breeding programs can be improved by genetic markers associated with production and reproduction traits.Reproduction traits are important for economic success of pig production and therefore development of genetictools for selection is of high interest to pig breeding. In this study our objective was to identify genomic regions associatedwith fertility traits in two Finnish pig breeds using large scale SNP genotyping and genome wide associationanalysis and characterization of copy number variations (CNV). Since CNVs are structural variations of the genomethey potentially have a large effect on gene expression and protein function. We analyzed 1265 genotyped boarsfor nine different reproduction traits and identified 46 CNV regions encompassing 13 genes. 11 of the CNV regionswere shared between the breeds, 20 were unique to the Finnish Yorkshire and 15 to the Finnish Landrace. The genomewide association (GWAS) analysis identified zero to five reproduction associated genomic regions per trait.Furthermore, we identified 23 genomic regions with 20 candidate genes associated with fertility traits using GWASanalysis. The identified CNV regions were compared against GWAS regions to detect candidate regions with an effecton reproduction traits. This study reports candidate genes and genomic regions within two Finnish pig breedsfor reproduction traits, which can be utilized in breeding programs.

Copy Number Variants in Two Northernmost Cattle Breeds Are Related to Their Adaptive Phenotypes.

Copy number variations (CNVs) are genomic structural variants with potential functional and evolutionary effects on phenotypes. In this study, we report the identification and characterization of CNVs from the whole-genome resequencing data of two northernmost cattle breeds from Russia: the Yakut and Kholmogory cattle and their phylogenetically most related breeds, Hanwoo and Holstein, respectively. Comparisons of the CNV regions (CNVRs) among the breeds led to the identification of breed-specific CNVRs shared by cold-adapted Kholmogory and Yakut cattle. An investigation of their overlap with genes, regulatory domains, conserved non-coding elements (CNEs), enhancers, and quantitative trait loci (QTLs) was performed to further explore breed-specific biology and adaptations. We found CNVRs enriched for gene ontology terms related to adaptation to environments in both the Kholmogory and Yakut breeds and related to thermoregulation specifically in Yakut cattle. Interestingly, the latter has also been supported when exploring the enrichment of breed-specific CNVRs in the regulatory domains and enhancers, CNEs, and QTLs implying the potential contribution of CNVR to the Yakut and Kholmogory cattle breeds’ adaptation to a harsh environment.

Copy Number Variation (CNV): A New Genomic Insight in Horses.

Simple SummaryThis study aimed to contribute to our knowledge of CNVs, a type of genomic marker in equines, by producing, for the first time, a fine-scale characterization of the CNV regions (CNVRs) in the Pura Raza Española horse breed. We found not only the existence of a unique pattern of genomic regions enriched in CNVs in the PRE in comparison with the data available from other breeds but also the incidence of CNVs across the entire genome. Since these regions could affect the structure and dose of the genes involved, we also performed a gene ontology analysis which revealed that most of the genes overlapping in CNVRs were related to the olfactory pathways and immune response.Copy number variations (CNVs) are a new-fangled source of genetic variation that can explain changes in the phenotypes in complex traits and diseases. In recent years, their study has increased in many livestock populations. However, the study and characterization of CNVs in equines is still very limited. Our study aimed to investigate the distribution pattern of CNVs, characterize CNV regions (CNVRs), and identify the biological pathways affected by CNVRs in the Pura Raza Española (PRE) breed. To achieve this, we analyzed high-density SNP genotyping data (670,804 markers) from a large cohort of 654 PRE horses. In total, we identified 19,902 CNV segments and 1007 CNV regions in the whole population. The length of the CNVs ranged from 1.024 kb to 4.55 Mb, while the percentage of the genome covered by CNVs was 4.4%. Interestingly, duplications were more abundant than deletions and mixed CNVRs. In addition, the distribution of CNVs across the chromosomes was not uniform, with ECA12 being the chromosome with the largest percentage of its genome covered (19.2%), while the highest numbers of CNVs were found in ECA20, ECA12, and ECA1. Our results showed that 71.4% of CNVRs contained genes involved in olfactory transduction, olfactory receptor activity, and immune response. Finally, 39.1% of the CNVs detected in our study were unique when compared with CNVRs identified in previous studies. To the best of our knowledge, this is the first attempt to reveal and characterize the CNV landscape in PRE horses, and it contributes to our knowledge of CNVs in equines, thus facilitating the understanding of genetic and phenotypic variations in the species. However, further research is still needed to confirm if the CNVs observed in the PRE are also linked to variations in the specific phenotypical differences in the breed.

Clinical Characterization of Copy Number Variants Associated With Neurodevelopmental Disorders in a Large-scale Multiancestry Biobank

Past studies identified rare copy number variants (CNVs) as risk factors for neurodevelopmental disorders (NDDs), including autism spectrum disorder and schizophrenia. However, the clinical characterization of NDD CNVs is understudied in population cohorts unselected for neuropsychiatric disorders and in cohorts of diverse ancestry. To identify individuals harboring NDD CNVs in a multiancestry biobank and to query their enrichment for select neuropsychiatric disorders as well as association with multiple medical disorders. In a series of phenotypic enrichment and association analyses, NDD CNVs were clinically characterized among 24 877 participants in the BioMe biobank, an electronic health record-linked biobank derived from the Mount Sinai Health System, New York, New York. Participants were recruited into the biobank since September 2007 across diverse ancestry and medical and neuropsychiatric specialties. For the current analyses, electronic health record data were analyzed from May 2004 through May 2019. NDD CNVs were identified using a consensus of 2 CNV calling algorithms, based on whole-exome sequencing and genotype array data, followed by novel in-silico clinical assessments. Of 24 877 participants, 14 586 (58.7%) were female; self-reported ancestry categories included 5965 (24.0%) who were of African ancestry, 7892 (31.7%) who were of European ancestry, and 8536 (34.3%) who were of Hispanic ancestry; and the mean (SD) age was 50.5 (17.3) years. Among 24 877 individuals, the prevalence of 64 NDD CNVs was 2.5% (n = 627), with prevalence varying by locus, corroborating the presence of some relatively highly prevalent NDD CNVs (eg, 15q11.2 deletion/duplication). An aggregate set of NDD CNVs were enriched for congenital disorders (odds ratio, 2.0; 95% CI, 1.1-3.5; P = .01) and major depressive disorder (odds ratio, 1.5; 95% CI, 1.1-2.0; P = .01). In a meta-analysis of medical diagnoses (n = 195 hierarchically clustered diagnostic codes), NDD CNVs were significantly associated with several medical outcomes, including essential hypertension (z score = 3.6; P = 2.8 × 10-4), kidney failure (z score = 3.3; P = 1.1 × 10-3), and obstructive sleep apnea (z score = 3.4; P = 8.1 × 10-4) and, in another analysis, morbid obesity (z score = 3.8; P = 1.3 × 10-4). Further, NDD CNVs were associated with increased body mass index in a multiancestry analysis (β = 0.19; 95% CI, 0.10-0.31; P = .003). For 36 common serum tests, there was no association with NDD CNVs. Clinical features of individuals harboring NDD CNVs were elucidated in a large-scale, multiancestry biobank, identifying enrichments for congenital disorders and major depressive disorder as well as associations with several medical outcomes, including hypertension, kidney failure, and obesity and obesity-related phenotypes, specifically obstructive sleep apnea and increased body mass index. The association between NDD CNVs and obesity outcomes indicate further potential pleiotropy of NDD CNVs beyond neurodevelopmental outcomes previously reported. Future clinical genetic investigations may lead to insights of at-risk individuals and therapeutic strategies targeting specific genetic variants. The importance of diverse inclusion within biobanks and considering the effect of rare genetic variants in a multiancestry context is evident.

KOMB: K-core based de novo characterization of copy number variation in microbiomes

Characterizing metagenomes via kmer-based, database-dependent taxonomic classification has yielded key insights into underlying microbiome dynamics. However, novel approaches are needed to track community dynamics and genomic flux within metagenomes, particularly in response to perturbations. We describe KOMB, a novel method for tracking genome level dynamics within microbiomes. KOMB utilizes K-core decomposition to identify Structural variations (SVs), specifically, population-level Copy Number Variation (CNV) within microbiomes. K-core decomposition partitions the graph into shells containing nodes of induced degree at least K, yielding reduced computational complexity compared to prior approaches. Through validation on a synthetic community, we show that KOMB recovers and profiles repetitive genomic regions in the sample. KOMB is shown to identify functionally-important regions in Human Microbiome Project datasets, and was used to analyze longitudinal data and identify keystone taxa in Fecal Microbiota Transplantation (FMT) samples. In summary, KOMB represents a novel graph-based, taxonomy-oblivious, and reference-free approach for tracking CNV within microbiomes. KOMB is open source and available for download athttps://gitlab.com/treangenlab/komb.

Comprehensive characterization of copy number variation (CNV) called from array, long- and short-read data

BackgroundSNP arrays, short- and long-read genome sequencing are genome-wide high-throughput technologies that may be used to assay copy number variants (CNVs) in a personal genome. Each of these technologies comes with its own limitations and biases, many of which are well-known, but not all of them are thoroughly quantified.ResultsWe assembled an ensemble of public datasets of published CNV calls and raw data for the well-studied Genome in a Bottle individual NA12878. This assembly represents a variety of methods and pipelines used for CNV calling from array, short- and long-read technologies. We then performed cross-technology comparisons regarding their ability to call CNVs. Different from other studies, we refrained from using the golden standard. Instead, we attempted to validate the CNV calls by the raw data of each technology.ConclusionsOur study confirms that long-read platforms enable recalling CNVs in genomic regions inaccessible to arrays or short reads. We also found that the reproducibility of a CNV by different pipelines within each technology is strongly linked to other CNV evidence measures. Importantly, the three technologies show distinct public database frequency profiles, which differ depending on what technology the database was built on.

Gene Variant of Barrier to Autointegration Factor 2 (Banf2w) Is Concordant with Female Determination in Cichlids.

Oreochromis fishes exhibit variability of sex-determination (SD) genes whose characterization contributes to understanding of the sex differentiation network, and to effective tilapia farming, which requires all-male culture. However, O. niloticus (On) amh is the only master-key regulator (MKR) of SD that has been mapped (XY/XX SD-system on LG23). In O. aureus (Oa), LG3 controls a WZ/ZZ SD-system that has recently been delimited to 9.2 Mbp, with an embedded interval rich with female-specific variation, harboring two paics genes and banf2. Developing genetic markers within this interval and using a hybrid Oa stock that demonstrates no recombination repression in LG3, we mapped the critical SD region to 235 Kbp on the orthologous On physical map (p < 1.5 × 10−26). DNA-seq assembly and peak-proportion analysis of variation based on Sanger chromatograms allowed the characterization of copy-number variation (CNV) of banf2. Oa males had three exons capable of encoding 90-amino-acid polypeptides, yet in Oa females, we found an extra copy with an 89-amino-acid polypeptide and three non-conservative amino acid substitutions, designated as banf2w. CNV analysis suggested the existence of two to five copies of banf2 in diploidic Cichlidae. Disrupting the Hardy–Weinberg equilibrium (p < 4.2 × 10−3), banf2w was concordant with female determination in Oa and in three cichlids with LG3 WZ/ZZ SD-systems (O. tanganicae, O. hornorum and Pelmatolapia mariae). Furthermore, exclusive RNA-seq expression in Oa females strengthened the candidacy of banf2w as the long-sought LG3 SD MKR. As banf genes mediate nuclear assembly, chromatin organization, gene expression and gonad development, banf2w may play a fundamental role inducing female nucleus formation that is essential for WZ/ZZ SD.

Characterization of Copy-Number Variations and Possible Candidate Genes in Recurrent Pregnancy Losses.

It is well established that embryonic chromosomal abnormalities (both in the number of chromosomes and the structure) account for 50% of early pregnancy losses. However, little is known regarding the potential differences in the incidence and distribution of chromosomal abnormalities between patients with sporadic abortion (SA) and recurrent pregnancy loss (RPL), let alone the role of submicroscopic copy-number variations (CNVs) in these cases. The aim of the present study was to systematically evaluate the role of embryonic chromosomal abnormalities and CNVs in the etiology of RPL compared with SA. Over a 3-year period, 1556 fresh products of conception (POCs) from miscarriage specimens were investigated using single nucleotide polymorphism array (SNP-array) and CNV sequencing (CNV-seq) in this study, along with further functional enrichment analysis. Chromosomal abnormalities were identified in 57.52% (895/1556) of all cases. Comparisons of the incidence and distributions of chromosomal abnormalities within the SA group and RPL group and within the different age groups were performed. Moreover, 346 CNVs in 173 cases were identified, including 272 duplications, 2 deletions and 72 duplications along with deletions. Duplications in 16q24.3 and 16p13.3 were significantly more frequent in RPL cases, and thereby considered to be associated with RPL. There were 213 genes and 131 signaling pathways identified as potential RPL candidate genes and signaling pathways, respectively, which were centered primarily on six functional categories. The results of the present study may improve our understanding of the etiologies of RPL and assist in the establishment of a population-based diagnostic panel of genetic markers for screening RPL amongst Chinese women.

Genome-wide characterization of copy number variations in the host genome in genetic resistance to Marek\u2019s disease using next generation sequencing

BackgroundMarek’s disease (MD) is a highly neoplastic disease primarily affecting chickens, and remains as a chronic infectious disease that threatens the poultry industry. Copy number variation (CNV) has been examined in many species and is recognized as a major source of genetic variation that directly contributes to phenotypic variation such as resistance to infectious diseases. Two highly inbred chicken lines, 63 (MD-resistant) and 72 (MD-susceptible), as well as their F1 generation and six recombinant congenic strains (RCSs) with varied susceptibility to MD, are considered as ideal models to identify the complex mechanisms of genetic and molecular resistance to MD.ResultsIn the present study, to unravel the potential genetic mechanisms underlying resistance to MD, we performed a genome-wide CNV detection using next generation sequencing on the inbred chicken lines with the assistance of CNVnator. As a result, a total of 1649 CNV regions (CNVRs) were successfully identified after merging all the nine datasets, of which 90 CNVRs were overlapped across all the chicken lines. Within these shared regions, 1360 harbored genes were identified. In addition, 55 and 44 CNVRs with 62 and 57 harbored genes were specifically identified in line 63 and 72, respectively. Bioinformatics analysis showed that the nearby genes were significantly enriched in 36 GO terms and 6 KEGG pathways including JAK/STAT signaling pathway. Ten CNVRs (nine deletions and one duplication) involved in 10 disease-related genes were selected for validation by using quantitative real-time PCR (qPCR), all of which were successfully confirmed. Finally, qPCR was also used to validate two deletion events in line 72 that were definitely normal in line 63. One high-confidence gene, IRF2 was identified as the most promising candidate gene underlying resistance and susceptibility to MD in view of its function and overlaps with data from previous study.ConclusionsOur findings provide valuable insights for understanding the genetic mechanism of resistance to MD and the identified gene and pathway could be considered as the subject of further functional characterization.

Genome-wide characterization of copy number variations in diffuse large B-cell lymphoma with implications in targeted therapy.

Diffuse large B-cell lymphoma (DLBCL) is the aggressive form of haematological malignancies with relapse/refractory in ~ 40% of cases. It mostly develops due to accumulation of various genetic and epigenetic variations that contribute to its aggressiveness. Though large-scale structural alterations have been reported in DLBCL, their functional role in pathogenesis and as potential targets for therapy is not yet well understood. In this study we performed detection and analysis of copy number variations (CNVs) in 11 human DLBCL cell lines (4 activated B-cell–like [ABC] and 7 germinal-centre B-cell–like [GCB]), that serve as model systems for DLBCL cancer cell biology. Significant heterogeneity observed in CNV profiles of these cell lines and poor prognosis associated with ABC subtype indicates the importance of individualized screening for diagnostic and prognostic targets. Functional analysis of key cancer genes exhibiting copy alterations across the cell lines revealed activation/disruption of ten potentially targetable immuno-oncogenic pathways. Genome guided in silico therapy that putatively target these pathways is elucidated. Based on our analysis, five CNV-genes associated with worst survival prognosis are proposed as potential prognostic markers of DLBCL.

Characterization of Copy Number Variations in Oral Cavity Squamous Cell Carcinoma Reveals a Novel Role for MLLT3 in Cell Invasiveness.

DNA copy number variations (CNVs) are a hallmark of cancer, and the current study aimed to demonstrate the profile of the CNVs for oral cavity squamous cell carcinoma (OSCC) and elucidate the clinicopathological associations and molecular mechanisms of a potential marker derived from CNVs, mixed-lineage leukemia translocated to chromosome 3 protein (MLLT3), in OSCC carcinogenesis. CNVs in 37 OSCC tissue specimens were analyzed using a high-resolution microarray, the OncoScan array. Gene expression was analyzed by real-time polymerase chain reaction in 127 OSCC and normal tissue samples. Cell function assays included cell cycle, migration, invasion and chromatin immunoprecipitation assays. We found a novel copy number amplified region, chromosome 9p, encompassing MLLT3 via the comparison of our data set with six other OSCC genome-wide CNV data sets. MLLT3 overexpression was associated with poorer overall survival in patients with OSCC (p = .048). MLLT3 knockdown reduced cell migration and invasion. The reduced invasion ability in MLLT3-knockdown cells was rescued with double knockdown of MLLT3 and CBP/p300-interacting transactivator with ED rich carboxy-terminal domain 4 (CITED4; 21.0% vs. 61.5%). Knockdown of MLLT3 impaired disruptor of telomeric silencing-1-like (Dot1L)-associated hypermethylation in the promoter of the tumor suppressor, CITED4 (p < .001), and hence dysregulated HIF-1α-mediated genes (TWIST, MMP1, MMP2, VIM, and CDH1) in OSCC cells. We identified unique CNVs in tumors of Taiwanese patients with OSCC. Notably, MLLT3 overexpression is related to the poorer prognosis of patients with OSCC and is required for Dot1L-mediated transcriptional repression of CITED4, leading to dysregulation of HIF-1α-mediated genes. This article reports unique copy number variations in oral cavity squamous cell carcinoma (OSCC) tumors of Taiwanese patients. Notably, MLLT3 overexpression is related to the poorer prognosis of patients with OSCC and is required for Dot1L-mediated transcriptional repression of CITED4, leading to dysregulation of HIF-1α-mediated genes.

Characterization of copy number variants in a large multibreed population of beef and dairy cattle using high-density single nucleotide polymorphism genotype data.

Copy number variants (CNVs) are a form of genomic variation that changes the structure of the genome through deletion or duplication of stretches of DNA. The objective of the present study was to characterize CNVs in a large multibreed population of beef and dairy bulls. The CNVs were called on the autosomes of 5,551 cattle from 22 different beef and dairy breeds, using 2 freely available software suites, QuantiSNP and PennCNV. All CNVs were classified into either deletions or duplications. The median concordance between PennCNV and QuantiSNP, per animal, was 18.5% for deletions and 0% for duplications. The low concordance rate between PennCNV and QuantiSNP indicated that neither algorithm, by itself, could identify all CNVs in the population. In total, PennCNV and QuantiSNP collectively identified 747,129 deletions and 432,523 duplications; 80.2% of all duplications and 69.1% of all deletions were present only once in the population. Only 0.154% of all CNVs identified were present in more than 50 animals in the population. The distribution of the percentage of the autosomes that were composed of deletions, per animal, was positively skewed, as was the distribution for the percentage of the autosomes that were composed of duplications, per animal. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of deletions were 0.019%, 0.037%, and 0.201%, respectively. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of duplications were 0.013%, 0.028%, and 0.076%, respectively. The distributions of the number of deletions and duplications per animal were both positively skewed. The interquartile range for the number of deletions per animal in the population was between 16 and 117, whereas for duplications it was between 8 and 23. Per animal, there tended to be twice as many deletions as duplications. The distribution of the length of deletions was positively skewed, as was the distribution of the length of duplications. The interquartile range for the length of deletions in the population was between 25 and 101 kb, and for duplications the interquartile range was between 46 and 235 kb. Per animal, duplications tended to be twice as long as deletions. This study provides a description of the characteristics and distribution of CNVs in a large multibreed population of beef and dairy cattle.

A population genomic characterization of copy number variation in the opportunistic fungal pathogen Aspergillus fumigatus.

Aspergillus fumigatus is a potentially deadly opportunistic fungal pathogen. Molecular studies have shaped our understanding of the genes, proteins, and molecules that contribute to A. fumigatus pathogenicity, but few studies have characterized genome-wide patterns of genetic variation at the population level. Of A. fumigatus genomic studies to-date, most focus mainly on single nucleotide polymorphisms and large structural variants, while overlooking the contribution of copy number variation (CNV). CNV is a class of small structural variation defined as loci that vary in their number of copies between individuals due to duplication, gain, or deletion. CNV can influence phenotype, including fungal virulence. In the present study, we characterized the population genomic patterns of CNV in a diverse collection of 71 A. fumigatus isolates using publicly available sequencing data. We used genome-wide single nucleotide polymorphisms to infer the population structure of these isolates and identified three populations consisting of at least 8 isolates. We then computationally predicted genome-wide CNV profiles for each isolate and conducted analyses at the species-, population-, and individual levels. Our results suggest that CNV contributes to genetic variation in A. fumigatus, with ~10% of the genome being CN variable. Our analysis indicates that CNV is non-randomly distributed across the A. fumigatus genome, and is overrepresented in subtelomeric regions. Analysis of gene ontology categories in genes that overlapped CN variants revealed an enrichment of genes related to transposable element and secondary metabolism functions. We further identified 72 loci containing 33 genes that showed divergent copy number profiles between the three A. fumigatus populations. Many of these genes encode proteins that interact with the cell surface or are involved in pathogenicity. Our results suggest that CNV is an important source of genetic variation that could account for some of the phenotypic differences between A. fumigatus populations and isolates.

Global characterization of copy number variants in epilepsy patients from whole genome sequencing.

Epilepsy will affect nearly 3% of people at some point during their lifetime. Previous copy number variants (CNVs) studies of epilepsy have used array-based technology and were restricted to the detection of large or exonic events. In contrast, whole-genome sequencing (WGS) has the potential to more comprehensively profile CNVs but existing analytic methods suffer from limited accuracy. We show that this is in part due to the non-uniformity of read coverage, even after intra-sample normalization. To improve on this, we developed PopSV, an algorithm that uses multiple samples to control for technical variation and enables the robust detection of CNVs. Using WGS and PopSV, we performed a comprehensive characterization of CNVs in 198 individuals affected with epilepsy and 301 controls. For both large and small variants, we found an enrichment of rare exonic events in epilepsy patients, especially in genes with predicted loss-of-function intolerance. Notably, this genome-wide survey also revealed an enrichment of rare non-coding CNVs near previously known epilepsy genes. This enrichment was strongest for non-coding CNVs located within 100 Kbp of an epilepsy gene and in regions associated with changes in the gene expression, such as expression QTLs or DNase I hypersensitive sites. Finally, we report on 21 potentially damaging events that could be associated with known or new candidate epilepsy genes. Our results suggest that comprehensive sequence-based profiling of CNVs could help explain a larger fraction of epilepsy cases.