A large deletion in a non-coding regulatory region leads to NFKB1 haploinsufficiency in two adult siblings

Mutations in NFkB pathway genes can cause inborn errors of immunity (IEI), with NFKB1 haploinsufficiency being a significant etiology for common variable immunodeficiency (CVID). Indeed, mutations in NFKB1 are found in 4 to 5% of in European and United States CVID cohorts, respectively; CVID representing almost ¼ of IEI patients in European countries registries. This case study presents a 49-year-old patient with respiratory infections, chronic diarrhea, immune thrombocytopenia, hypogammaglobulinemia, and secondary lymphoma. Comprehensive genetic analysis, including high-throughput sequencing of 300 IEI-related genes and copy number variation analysis, identified a critical 2.6-kb deletion spanning the first untranslated exon and its upstream region. The region's importance was confirmed through genetic markers indicative of enhancers and promoters. The deletion was also found in the patient's brother, who displayed similar but milder symptoms. Functional analysis supported haploinsufficiency with reduced mRNA and protein expression in both patients. This case underscores the significance of copy number variation (CNV) analysis and targeting noncoding exons within custom gene panels, emphasizing the broader genomic approaches needed in medical genetics.

Association of Structural Forms of 17q21.31 with the Risk of Progressive Supranuclear Palsy and MAPT Sub-haplotypes

Importance

The chromosome 17q21.31 region, containing a 900 Kb inversion that defines H1 and H2 haplotypes, represents the strongest genetic risk locus in progressive supranuclear palsy (PSP). In addition to H1 and H2, various structural forms of 17q21.31, characterized by the copy number of α, β, and γ duplications, have been identified. However, the specific effect of each structural form on the risk of PSP has never been evaluated in a large cohort study.

Objective

To assess the association of different structural forms of 17q.21.31, defined by the copy numbers of α, β, and γ duplications, with the risk of PSP and MAPT sub-haplotypes.

Design setting and participants

Utilizing whole genome sequencing data of 1,684 (1,386 autopsy confirmed) individuals with PSP and 2,392 control subjects, a case-control study was conducted to investigate the association of copy numbers of α, β, and γ duplications and structural forms of 17q21.31 with the risk of PSP. All study subjects were selected from the Alzheimer's Disease Sequencing Project (ADSP) Umbrella NG00067.v7. Data were analyzed between March 2022 and November 2023.

Main outcomes and measures

The main outcomes were the risk (odds ratios [ORs]) for PSP with 95% CIs. Risks for PSP were evaluated by logistic regression models.

Results

The copy numbers of α and β were associated with the risk of PSP only due to their correlation with H1 and H2, while the copy number of γ was independently associated with the increased risk of PSP. Each additional duplication of γ was associated with 1.10 (95% CI, 1.04-1.17; P = 0.0018) fold of increased risk of PSP when conditioning H1 and H2. For the H1 haplotype, addition γ duplications displayed a higher odds ratio for PSP: the odds ratio increases from 1.21 (95%CI 1.10-1.33, P = 5.47 × 10-5) for H1β1γ1 to 1.29 (95%CI 1.16-1.43, P = 1.35 × 10-6) for H1β1γ2, 1.45 (95%CI 1.27-1.65, P = 3.94 × 10-8) for H1β1γ3, and 1.57 (95%CI 1.10-2.26, P = 1.35 × 10-2) for H1β1γ4. Moreover, H1β1γ3 is in linkage disequilibrium with H1c (R2 = 0.31), a widely recognized MAPT sub-haplotype associated with increased risk of PSP. The proportion of MAPT sub-haplotypes associated with increased risk of PSP (i.e., H1c, H1d, H1g, H1o, and H1h) increased from 34% in H1β1γ1 to 77% in H1β1γ4.

Conclusions and relevance

This study revealed that the copy number of γ was associated with the risk of PSP independently from H1 and H2. The H1 haplotype with more γ duplications showed a higher odds ratio for PSP and were associated with MAPT sub-haplotypes with increased risk of PSP. These findings expand our understanding of how the complex structure at 17q21.31 affect the risk of PSP.

Establishment and characterization of a recurrent malignant peripheral nerve sheath tumor cell line: RsNF

Malignant peripheral nerve sheath tumor (MPNST) is a highly aggressive and recurrent soft tissue sarcoma. It most commonly occurs secondary to neurofibromatosis type I, and it has a 5-year survival rate of only 8-13%. To better study the tumor heterogeneity of MPNST and to develop diverse treatment options, more tumor-derived cell lines are needed to obtain richer biological information. Here, we established a primary cell line of relapsed MPNST RsNF cells derived from a patient diagnosed with NF1 and detected the presence of NF1 mutations and SUZ12 somatic mutations through whole-exome sequencing(WES). Through tumor molecular marker targeted sequencing and single-cell transcriptome sequencing, it was found that chromosome 7 copy number variation (CNV) was gained in this cell line, and ZNF804B, EGFR, etc., were overexpressed on chromosome 7. Therefore, RsNF cells can be used as a useful tool in NF1-associated MPNST genomic amplification studies and to develop new therapeutic strategies.

High overall copy number variation burden by genome-wide methylation profiling holds negative prognostic value in surgically treated pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive type of cancer with an overall 5-year survival of around 10 %. New prognostic tools to stratify patients are needed. Our main aim was to evaluate the prognostic value of overall copy number variation (CNV) burden in surgically treated PDAC. DNA extracted from 108 surgical PDAC specimens was examined to collect data on the genome-wide DNA methylation status of >850,000 CpG sites in promoter, gene body, and enhancer regions (Illumina Infinium Methylation EPIC BeadChip Kit). CNV profiles were obtained and all PDACs were stratified into one of three groups: Low, moderate, or high overall CNV burden. Tumors histologically showing a dominant conventional and/or tubulopapillary pattern in 60 %-100 % and 0-59 % were categorized as Group A and Group B as per Kalimuthu. We also performed targeted next-generation sequencing (NGS) and immunohistochemistry. High overall CNV burden held independent negative prognostic value with poor survival (HR 4.01 (95%CI 1.96-8.19), p = 0.00014) and was more frequent in Group B (p = 0.0003). Most frequent chromosomal arm-level aberrations were gains of 8q (29 %) and 1q (19 %) and losses of 17p (55 %), 18q (43 %), 6q (37 %), 9p (36 %), 6p (26 %), 19p (26 %), and 8p (25 %). Most frequent mutations found were in KRAS (95 %), TP53 (62 %), CDKN2A (24 %), SMAD4 (23 %), ATM (9 %), ARID1A (7 %), RNF43 (7 %), GNAS (6 %), and KDM6A (6 %). Group A PDACs showed more frequently KRAS variants other than Gly12Val and Gly12Asp (p = 0.012). Our data indicate that overall CNV burden using genome-wide methylation profiling may be a useful prognostic tool in surgically treated PDAC. Importantly, our approach, using data from genome-wide methylation profiling for analysis of overall CNV burden, can be performed on formalin-fixed and paraffin embedded PDAC tissues. Future studies should examine the prognostic value of overall CNV burden in unresectable PDAC.

Chromosome microarray analysis combined with karyotype analysis is a powerful tool for the detection in pregnant women with high-risk indicators

BACKGROUND

Karyotype analysis and fluorescence in situ hybridization (FISH) are commonly used for prenatal diagnosis, however they have many disadvantages. Chromosome microarray analysis (CMA) has the potential to overcome these disadvantages. This study aimed to evaluate the clinical value of CMA in the diagnosis of fetal chromosomal anomalies in southwest of China.

METHODS

A total of 3336 samples of amniotic fluid or umbilical cord blood from pregnant women with high-risk indicators at our center in southwest of China from June 2018 to January 2023 were included in the retrospective analysis. 3222 cases tested by CMA and karyotyping, 114 cases only tested by CMA.

RESULTS

3336 samples divided into 2911 cases with single and 425 cases with multiple high-risk indicators. The aneuploidy and pathogenic/likely pathogenic copy number variations (CNVs) of 2911 cases with single high-risk indicator were 4.43% (129/2911) and 2.44% (71/2911) respectively; the aneuploidy and pathogenic/likely pathogenic CNVs of 425 cases with multiple high-risk indicators were 6.82% (29/425) and 2.12% (9/425) respectively. The rate of aneuploidy increased significantly with pregnancy age or NT value. The detection rate of aneuploidy on cases with AMA combined NT ≥ 2.5 mm was significantly higher than that in cases only with AMA (p < 0.001); the detection rate of aneuploidy and pathogenic/likely pathogenic CNVs in cases with AMA combined NIPT high-risk were higher than that in cases only with AMA (p < 0.001, p < 0.05).

CONCLUSIONS

The combined application of CMA and karyotyping were recommended in prenatal diagnosis for providing a scientific and accurate genetic diagnosis and improving the quality of prenatal genetic counseling.

Epilepsy plus blindness in microdeletion of GABRA1 and GABRG2 in mouse and human

OBJECTIVE

GABAA receptor subunit gene (GABR) mutations are significant causes of epilepsy, including syndromic epilepsy. This report for the first time, describes intractable epilepsy and blindness due to optic atrophy in our patient, who has a microdeletion of the GABRA1 and GABRG2 genes. We then characterized the molecular phenotypes and determined patho-mechanisms underlying the genotype-phenotype correlations in a mouse model who is haploinsufficient for both genes (Gabra1+/-/Gabrg2+/- mouse).

METHODS

Electroencephalography was conducted in both human and mice with the same gene loss. GABAA receptor expression was evaluated by biochemical and imaging approaches. Optic nerve atrophy was evaluated with fundus photography in human while electronic microscopy, visual evoked potential and electroretinography recordings were conducted in mice.

RESULTS

The patient has bilateral optical nerve atrophy. Mice displayed spontaneous seizures, reduced electroretinography oscillatory potential and reduced GABAA receptor α1, β2 and γ2 subunit expression in various brain regions. Electronic microscopy showed that mice also had optic nerve degeneration, as indicated by increased G-ratio, the ratio of the inner axonal diameter to the total outer diameter, suggesting impaired myelination of axons. More importantly, we identified that phenobarbital was the most effective anticonvulsant in mice and the patient's seizures were also controlled with phenobarbital after failing multiple anti-seizure drugs.

CONCLUSIONS

This study is the first report of haploinsufficiency of two GABR epilepsy genes and visual impairment due to altered axonal myelination and resultant optic nerve atrophy. The study suggests the far-reaching impact of GABR mutations and the translational significance of animal models with the same etiology.

Copy number variations (CNVs) and karyotyping analysis in males with azoospermia and oligospermia

BACKGROUND

Considering the essential roles that genetic factors play in azoospermia and oligospermia, this study aims to identify abnormal chromosomes using karyotyping and CNVs and elucidate the associated genes in patients.

METHODS

A total of 1157 azoospermia and oligospermia patients were recruited, of whom, 769 and 674 underwent next-generation sequencing (NGS) to identify CNVs and routine G-band karyotyping, respectively.

RESULTS

First, 286 patients were co-analyzed using CNV sequencing (CNV-seq) and karyotyping. Of the 725 and 432 patients with azoospermia and oligospermia, 33.8% and 48.9% had abnormal karyotypes and CNVs, respectively. In particular, 47,XXY accounted for 44.18% and 26.33% of abnormal karyotypes and CNVs, respectively, representing the most frequent genetic aberration in azoospermia and oligospermia patients. Nevertheless, big Y and small Y accounted for 7.46% and 16.67% of abnormal karyotypes, respectively. We also identified high-frequency CNVs-loci, such as Xp22.31 and 2p24.3, in azoospermia and oligospermia patients.

CONCLUSION

Sex chromosome and autosomal CNV loci, such as Xp22.31 and 2p24.3, as well as the associated genes, such as VCX and NACAP9, could be candidate spermatogenesis genes. The high-frequency abnormal karyotypes, CNV loci, and hot genes represent new targets for future research.

Detection of germline variants with pathogenic potential in 48 patients with familial colorectal cancer by using whole exome sequencing

BACKGROUND

Hereditary genetic mutations causing predisposition to colorectal cancer are accountable for approximately 30% of all colorectal cancer cases. However, only a small fraction of these are high penetrant mutations occurring in DNA mismatch repair genes, causing one of several types of familial colorectal cancer (CRC) syndromes. Most of the mutations are low-penetrant variants, contributing to an increased risk of familial colorectal cancer, and they are often found in additional genes and pathways not previously associated with CRC. The aim of this study was to identify such variants, both high-penetrant and low-penetrant ones.

METHODS

We performed whole exome sequencing on constitutional DNA extracted from blood of 48 patients suspected of familial colorectal cancer and used multiple in silico prediction tools and available literature-based evidence to detect and investigate genetic variants.

RESULTS

We identified several causative and some potentially causative germline variants in genes known for their association with colorectal cancer. In addition, we identified several variants in genes not typically included in relevant gene panels for colorectal cancer, including CFTR, PABPC1 and TYRO3, which may be associated with an increased risk for cancer.

CONCLUSIONS

Identification of variants in additional genes that potentially can be associated with familial colorectal cancer indicates a larger genetic spectrum of this disease, not limited only to mismatch repair genes. Usage of multiple in silico tools based on different methods and combined through a consensus approach increases the sensitivity of predictions and narrows down a large list of variants to the ones that are most likely to be significant.

Familial 5.29 Mb deletion in chromosome Xq22.1-q22.3 with a normal phenotype: a rare pedigree and literature review

BACKGROUND

Xq22.1-q22.3 deletion is a rare chromosome aberration. The purpose of this study was to identify the correlation between the phenotype and genotype of chromosome Xq22.1-q22.3 deletions.

METHODS

Chromosome aberrations were identified by copy number variation sequencing (CNV-seq) technology and karyotype analysis. Furthermore, we reviewed patients with Xq22.1-q22.3 deletions or a deletion partially overlapping this region to highlight the rare condition and analyse the genotype-phenotype correlations.

RESULTS

We described a female foetus who is the "proband" of a Chinese pedigree and carries a heterozygous 5.29 Mb deletion (GRCh37: chrX: 100,460,000-105,740,000) in chromosome Xq22.1-q22.3, which may affect 98 genes from DRP2 to NAP1L4P2. This deletion encompasses 7 known morbid genes: TIMM8A, BTK, GLA, HNRNPH2, GPRASP2, PLP1, and SERPINA7. In addition, the parents have a normal phenotype and are of normal intelligence. The paternal genotype is normal. The mother carries the same deletion in the X chromosome. These results indicate that the foetus inherited this CNV from her mother. Moreover, two more healthy female family members were identified to carry the same CNV deletion through pedigree analysis according to the next-generation sequencing (NGS) results. To our knowledge, this family is the first pedigree to have the largest reported deletion of Xq22.1-q22.3 but to have a normal phenotype with normal intelligence.

CONCLUSIONS

Our findings further improve the understanding of the genotype-phenotype correlations of chromosome Xq22.1-q22.3 deletions.This report may provide novel information for prenatal diagnosis and genetic counselling for patients who carry similar chromosome abnormalities.

Establishment and evaluation of digital PCR methods for HER2 copy number variation in breast cancer

Accurate measurement of human epidermal growth factor receptor 2 (HER2) copy number variation (CNV) is very important for guiding the tumor target therapy in breast cancer. Digital PCR (dPCR) is a sensitive and an absolute quantitative method, which can be used to detect HER2 CNV. Three HER2 exon-specific digital PCR assays along with three new reference genes assays (homo sapiens ribonuclease P RNA component H1 (RPPH1), glucose-6-phosphate isomerase (GPI), and chromosome 1 open reading frame 43 (C1ORF43), on different chromosomes) were established and validated by using standard reference material, 8 different cell lines and 110 clinical Formalin-fixed and paraffin-embedded (FFPE) samples. DPCR can achieve precise quantification of HER2 CNV by calculating the ratio of HER2/reference gene. The positive and negative coincidence rates were 98% (53/54) and 95% (53/56), respectively, compared with fluorescence in situ hybridization (FISH) diagnostic result 110 of FFPE samples. The common reference gene CEP17 used for FISH diagnostic was not suitable as single reference gene for HER2 CNV measurements by dPCR. The best practice of HER2 CNV determination by dPCR is to conduct the three duplex assays of H1 (HER2 exon 4) with the proposed three new reference genes, with a positive cut-off value of H1/RPPH1 ≥ 2.0 or H1/averaged reference gene ≥ 2.0. The proposed dPCR method in our study can accurately provide absolute copy number of HER2 and reference gene on an alternative chromosome, thus avoiding false negative caused by polysomy of chromosome 17. The improved molecular typing and diagnosis of breast cancer will better guide clinical medication.

Analysis of Copy Number Variation of DNA Repair/Damage Response Genes in Tumor Tissues

Cells experience increased genome instability through the course of disease development including cancer initiation and progression. Point mutations, insertion/deletions, translocations, and amplifications of both coding and noncoding regions all contribute to cancer phenotypes. Copy number variation (CNV), i.e., changes of the number of copies of nuclear DNA, occurs in the genome of even normal somatic cells. Studies to understand the effects of CNV on tumor development, especially aspects concerning tumor aggressiveness and the influence on outcomes of therapeutic modalities, have been reignited by the breakthrough technologies of the molecular genomics. This section discusses the significance of analyzing CNVs that cause simultaneous increase/decrease of clusters of genes, using the expression profile of XRCC1 with its neighbor genes LIG1, PNKP, and POLD1 as an example. Methods for CNV assay at the individual gene level on formalin-fixed, paraffin-embedded (FFPE) tissues using the NanoString nCounter technology will then be described.

Association of HSF1 gene copy number variation with growth traits in the Ashidan yak

Copy Number Variation (CNV) is the major manner for the variation of genome structure, which is associated with numerous important traits. The heat shock factor 1 (HSF1) gene is a stress response transcriptional regulator. It participates in the heat shock response, simultaneously participated in the development of tissue. The objective of this research was to explore the influence of CNV of the HSF1 gene on the growth traits of the Ashidan yak. In this study, the growth traits (withers height, body weight, chest girth, and body length) of 274 Ashidan yaks were divided into four stages (6, 12, 18, and 30 months old). Moreover, quantitative polymerase chain reaction (qPCR) was exploited for determining the HSF1 gene relative expression level, and SPSS software was utilized for the statistical analysis. The outcomes indicated that HSF1-CNV was significantly associated with body length (p < 0.05) and was extremely significant associated with withers height (p < 0.01) of 18-month-old Ashidan yaks. Besides, the HSF1 relative expression in heart and muscle was higher than that existed in other tissues (p < 0.01). The outcomes suggested that the CNV of HSF1 gene would affect the growth and development of the Ashidan yak, which is conducive to the early breeding of yak.

Comparative analyses of copy number variations between swamp and river buffalo

Domestic buffalo is an important livestock in the tropical and sub-tropical region, including two types: swamp and river buffalo. The swamp buffalo is mainly used as draft animal, while the river buffalo is raised for milk production. In this study, based on the new high-quality buffalo reference genome UOA_WB_1, we firstly investigated the copy number variants in buffalo using whole-genome Illumina sequencing. A total of 3,734 CNV regions (CNVRs) were detected in 106 buffalo population with a total length of 23,429,066 bp, corresponding to ∼ 0.88% of the water buffalo genome (UOA_WB_1). Our results revealed a clear population differentiation in CNV between swamp and river buffalo. In addition, a total of 667 highly differentiated CNVRs (covering 886 genes) were detected between river and swamp buffalo population. We detected a set of CNVR-overlapping genes associated with exercise, immunity, nerve, and milk trait which exhibited different copy numbers between swamp and river buffalo population. This study provides valuable genome variation resources for buffalo and would contribute to understanding the genetic differences between swamp and river buffalo.

Copy number variation (CNV) of the AHR gene in the Ashidan yak and its association with growth traits

Copy number variation (CNV) is a principal genomic structure variation affecting the gene expression through the dose-effect and change of gene regulatory region. It plays an important role in regulating the various complex traits of vertebrates. The aromatic hydrocarbon receptor (AHR) is a member of ligand-dependent transcription factors which belong to the alkaline helix-loop-helix PASS family. It is used as a conservative environmental sensor during biological evolution. This study, tracked the growth data (body weight, withers height, body length, chest girth) of 332 yaks in four stages (6, 12, 18, and 30 months) were tracked. The CNV of the yaks was analyzed using real-time quantitative PCR, and the correlation between CNV of AHR and yak growth traits was analyzed using the SPSS and R software. The AHR gene expression profiles were assessed in different tissues of the 18-month-old yak. The statistical analysis indicated the AHR-CNV of the Ashidan yak to significantly correlate with the body length (P < 0.05), and was found to be correlated with the withers height at 18 months old (P < 0.01) with extreme significance. To sum up, this study for the first time discussed the relationship between AHR-CNV and the growth traits of the Ashidan yak. The results indicated that the AHR gene might become a new molecular marker in the breeding yak.

Identification of a novel CNV at the EYA4 gene in a Chinese family with autosomal dominant nonsyndromic hearing loss

BACKGROUND

Hereditary hearing loss is a heterogeneous class of disorders that exhibits various patterns of inheritance and involves many genes. Variants in the EYA4 gene in DFNA10 are known to lead to postlingual, progressive, autosomal dominant nonsyndromic hereditary hearing loss.

PATIENTS AND METHODS

We collected a four-generation Chinese family with autosomal-dominant nonsyndromic hearing loss (ADNSHL). We applied targeted next-generation sequencing (TNGS) in three patients of this pedigree and whole-genome sequencing (WGS) in the proband. The intrafamilial cosegregation of the variant and the deafness phenotype were confirmed by PCR, gap-PCR and Sanger sequencing.

RESULTS

A novel CNV deletion at 6q23 in exons 8-11 of the EYA4 gene with a 10 bp insertion was identified by TNGS and WGS and segregated with the ADNSHL phenotypes.

CONCLUSIONS

Our results expanded the variant spectrum and genotype‒phenotype correlation of the EYA4 gene and autosomal dominant nonsyndromic hereditary hearing loss in Chinese Han individuals. WGS is an accurate and effective method for verifying the genomic features of CNVs.

Expanding the application of non-invasive prenatal testing in the detection of foetal chromosomal copy number variations

Purpose

The aim of this study was to assess the detection efficiency and clinical application value of non-invasive prenatal testing (NIPT) for foetal copy number variants (CNVs) in clinical samples from 39,002 prospective cases.

Methods

A total of 39,002 pregnant women who received NIPT by next-generation sequencing (NGS) with a sequencing depth of 6 M reads in our centre from January 2018 to April 2020 were enrolled. Chromosomal microarray analysis (CMA) was further used to diagnose suspected chromosomal aneuploidies and chromosomal microdeletion/microduplication for consistency assessment.

Results

A total of 473 pregnancies (1.213%) were positive for clinically significant foetal chromosome abnormalities by NIPT. This group comprised 99 trisomy 21 (T21, 0.254%), 30 trisomy 18 (T18, 0.077%), 25 trisomy 13 (T13, 0.064%), 155 sex chromosome aneuploidy (SCA, 0.398%), 69 rare trisomy (0.177%), and 95 microdeletion/microduplication syndrome (MMS, 0.244%) cases. Based on follow-up tests, the positive predictive values (PPVs) for the T21, T18, T13, SCA, rare trisomy, and MMS cases were calculated to be 88.89%, 53.33%, 20.00%, 40.22%, 4.88%, and 49.02%, respectively. In addition, the PPVs of CNVs of < 5 Mb, 5–10 Mb, and > 10 Mb were 54.55%, 38.46%, and 40.00%, respectively. Among the 95 cases with suspected CNVs, 25 were diagnosed as true positive and 26 cases as false positive; follow-up prenatal diagnosis by CMA was not performed for 44 cases. Moreover, among the 25 true positive cases, 10 were pathogenic, 3 were likely pathogenic, and 12 were of uncertain significance.

Conclusion

NIPT is not only suitable for screening T21, T18, T13, and SCA but also has potential significance for CNV detection. As combined with ultrasound, extended NIPT is effective for screening MMS. However, NIPT should not be recommended for whole-chromosome aneuploidy screening.

Detecting copy number variation in next generation sequencing data from diagnostic gene panels

Background

Detection of copy number variation (CNV) in genes associated with disease is important in genetic diagnostics, and next generation sequencing (NGS) technology provides data that can be used for CNV detection. However, CNV detection based on NGS data is in general not often used in diagnostic labs as the data analysis is challenging, especially with data from targeted gene panels. Wet lab methods like MLPA (MRC Holland) are widely used, but are expensive, time consuming and have gene-specific limitations. Our aim has been to develop a bioinformatic tool for CNV detection from NGS data in medical genetic diagnostic samples.

Results

Our computational pipeline for detection of CNVs in NGS data from targeted gene panels utilizes coverage depth of the captured regions and calculates a copy number ratio score for each region. This is computed by comparing the mean coverage of the sample with the mean coverage of the same region in other samples, defined as a pool. The pipeline selects pools for comparison dynamically from previously sequenced samples, using the pool with an average coverage depth that is nearest to the one of the samples. A sliding window-based approach is used to analyze each region, where length of sliding window and sliding distance can be chosen dynamically to increase or decrease the resolution. This helps in detecting CNVs in small or partial exons. With this pipeline we have correctly identified the CNVs in 36 positive control samples, with sensitivity of 100% and specificity of 91%. We have detected whole gene level deletion/duplication, single/multi exonic level deletion/duplication, partial exonic deletion and mosaic deletion. Since its implementation in mid-2018 it has proven its diagnostic value with more than 45 CNV findings in routine tests.

Conclusions

With this pipeline as part of our diagnostic practices it is now possible to detect partial, single or multi-exonic, and intragenic CNVs in all genes in our target panel. This has helped our diagnostic lab to expand the portfolio of genes where we offer CNV detection, which previously was limited by the availability of MLPA kits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01059-x.

Integrative analysis of prognostic long non-coding RNAs with copy number variation in bladder cancer

Copy number variations (CNVs), which can affect the role of long non-coding RNAs (lncRNAs), are important genetic changes seen in some malignant tumors. We analyzed lncRNAs with CNV to explore the relationship between lncRNAs and prognosis in bladder cancer (BLCA). Messenger RNA (mRNA) expression levels, DNA methylation, and DNA copy number data of 408 BLCA patients were subjected to integrative bioinformatics analysis. Cluster analysis was performed to obtain different subtypes and differently expressed lncRNAs and coding genes. Weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression gene and lncRNA modules. CNV-associated lncRNA data and their influence on cancer prognosis were assessed with Kaplan-Meier survival curve. Multi-omics integration analysis revealed five prognostic lncRNAs with CNV, namely NR2F1-AS1, LINC01138, THUMPD3-AS1, LOC101928489,and TMEM147-AS1,and a risk-score signature related to overall survival in BLCA was identified. Moreover, validated results in another independent Gene Expression Omnibus (GEO) dataset, GSE31684, were consistent with these results. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the mitogen-activated protein kinase (MAPK) signaling pathway, focal adhesion pathway, and Janus kinase-signal transducers and activators of transcription (JAK-STAT) signaling pathway were enriched in a high-risk score pattern, suggesting that imbalance in these pathways is closely related to tumor development. We revealed the prognosis-related lncRNAs by analyzing the expression profiles of lncRNAs and CNVs, which can be used as prognostic biomarkers for BLCA.

Population differentiated copy number variation of Bos taurus, Bos indicus and their African hybrids

Background

CNV comprises a large proportion in cattle genome and is associated with various traits. However, there were few population-scale comparison studies on cattle CNV.

Results

Here, autosome-wide CNVs were called by read depth of NGS alignment result and copy number variation regions (CNVRs) defined from 102 Eurasian taurine (EAT) of 14 breeds, 28 Asian indicine (ASI) of 6 breeds, 22 African taurine (AFT) of 2 breeds, and 184 African humped cattle (AFH) of 17 breeds. The copy number of every CNVRs were compared between populations and CNVRs with population differentiated copy numbers were sorted out using the pairwise statistics V_ST and Kruskal-Wallis test. Three hundred sixty-two of CNVRs were significantly differentiated in both statistics and 313 genes were located on the population differentiated CNVRs.

Conclusion

For some of these genes, the averages of copy numbers were also different between populations and these may be candidate genes under selection. These include olfactory receptors, pathogen-resistance, parasite-resistance, heat tolerance and productivity related genes. Furthermore, breed- and individual-level comparison was performed using the presence or copy number of the autosomal CNVRs. Our findings were based on identification of CNVs from short Illumina reads of 336 individuals and 39 breeds, which to our knowledge is the largest dataset for this type of analysis and revealed important CNVs that may play a role in cattle adaption to various environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07808-7.

Performances of NIPT for copy number variations at different sequencing depths using the semiconductor sequencing platform

Objective

To evaluate the performance of noninvasive prenatal testing (NIPT) and NIPT-PLUS for the detection of genome-wide microdeletion and microduplication syndromes (MMSs) at different sequencing depths. The NIPT sequencing depth was 0.15X, and the data volume was 3 million reads; the NIPT-PLUS sequencing depth was 0.4X, and the data volume was 8 million reads.

Methods

A cohort of 50,679 pregnancies was recruited. A total of 42,969 patients opted for NIPT, and 7710 patients opted for NIPT-PLUS. All high-risk cases were advised to undergo invasive prenatal diagnosis and were followed up.

Results

A total of 373 cases had a high risk of a copy number variation (CNV) as predicted by NIPT and NIPT-PLUS: NIPT predicted 250 high-risk CNVs and NIPT-PLUS predicted 123. NIPT-PLUS increased the detection rate by 1.02% (0.58% vs 1.60%, p < 0.001). A total of 291 cases accepted noninvasive prenatal diagnosis, with 197 cases of NIPT and 94 cases of NIPT-PLUS. The PPV of CNV > 10 Mb for NIPT-PLUS was significantly higher than that for NIPT (p = 0.02). The total PPV of NIPT-PLUS was 12.56% higher than that of NIPT (43.61% vs 30.96%, p = 0.03).

Conclusion

NIPT-PLUS had a better performance in detecting CNVs in terms of the total detection rate and total PPV. However, great care must be taken in presenting results and providing appropriate counseling to patients when deeper sequencing is performed in clinical practice.

CNVs in the 22q11.2 Chromosomal Region Should Be an Early Suspect in Infants with Congenital Cardiac Disease

BACKGROUND

Congenital heart disease (CHD) is the most common congenital malformation, it is frequently found as an isolated defect, and the etiology is not completely understood. Although most of the cases have multifactorial causes, they can also be secondary to chromosomal abnormalities, monogenic diseases, microduplications or microdeletions, among others. Copy number variations (CNVs) at 22q11.2 are associated with a variety of symptoms including CHD, thymic aplasia, and developmental and behavioral manifestations. We tested CNVs in the 22q11.2 chromosomal region by MLPA in a cohort of Colombian patients with isolated CHD to establish the frequency of these CNVs in the cohort.

METHODS

CNVs analysis of 22q11.2 by MLPA were performed in 32 patients with apparently isolate CHD during the neonatal period. Participants were enrolled from different hospitals in Bogotá, and they underwent a clinical assessment by a cardiologist and a clinical geneticist.

RESULTS

CNVs in the 22q11.2 chromosomal region were found in 7 patients (21.9%). The typical deletion was found in 6 patients (18.75%) and atypical 1.5 Mb duplication was found in 1 patient (3.1%).

CONCLUSIONS

CNVs in 22q11.2 is a common finding in patients presenting with isolated congenital cardiac disease, therefore these patients should be tested early despite the absence of other clinical manifestations. MLPA is a very useful molecular method and provides an accurate diagnosis.

Re-evaluation of single nucleotide variants and identification of structural variants in a cohort of 45 sudden unexplained death cases

Sudden unexplained death (SUD) takes up a considerable part in overall sudden death cases, especially in adolescents and young adults. During the past decade, many channelopathy- and cardiomyopathy-associated single nucleotide variants (SNVs) have been identified in SUD studies by means of postmortem molecular autopsy, yet the number of cases that remain inconclusive is still high. Recent studies had suggested that structural variants (SVs) might play an important role in SUD, but there is no consensus on the impact of SVs on inherited cardiac diseases. In this study, we searched for potentially pathogenic SVs in 244 genes associated with cardiac diseases. Whole-exome sequencing and appropriate data analysis were performed in 45 SUD cases. Re-analysis of the exome data according to the current ACMG guidelines identified 14 pathogenic or likely pathogenic variants in 10 (22.2%) out of the 45 SUD cases, whereof 2 (4.4%) individuals had variants with likely functional effects in the channelopathy-associated genes SCN5A and TRDN and 1 (2.2%) individual in the cardiomyopathy-associated gene DTNA. In addition, 18 structural variants (SVs) were identified in 15 out of the 45 individuals. Two SVs with likely functional impairment were found in the coding regions of PDSS2 and TRPM4 in 2 SUD cases (4.4%). Both were identified as heterozygous deletions, which were confirmed by multiplex ligation-dependent probe amplification. In conclusion, our findings support that SVs could contribute to the pathology of the sudden death event in some of the cases and therefore should be investigated on a routine basis in suspected SUD cases.

Cotton pan-genome retrieves the lost sequences and genes during domestication and selection

BACKGROUND

Millennia of directional human selection has reshaped the genomic architecture of cultivated cotton relative to wild counterparts, but we have limited understanding of the selective retention and fractionation of genomic components.

RESULTS

We construct a comprehensive genomic variome based on 1961 cottons and identify 456 Mb and 357 Mb of sequence with domestication and improvement selection signals and 162 loci, 84 of which are novel, including 47 loci associated with 16 agronomic traits. Using pan-genome analyses, we identify 32,569 and 8851 non-reference genes lost from Gossypium hirsutum and Gossypium barbadense reference genomes respectively, of which 38.2% (39,278) and 14.2% (11,359) of genes exhibit presence/absence variation (PAV). We document the landscape of PAV selection accompanied by asymmetric gene gain and loss and identify 124 PAVs linked to favorable fiber quality and yield loci.

CONCLUSIONS

This variation repertoire points to genomic divergence during cotton domestication and improvement, which informs the characterization of favorable gene alleles for improved breeding practice using a pan-genome-based approach.

Y chromosome structural variation in infertile men detected by targeted next-generation sequencing

PURPOSE

To provide a validated method to identify copy number variation (CNV) in regions of the Y chromosome of infertile men by next-generation sequencing (NGS).

METHODS

Semen analysis was used to determine the quality of semen and diagnose infertility. Deletion of the azoospermia factor (AZF) region in the Y chromosome was detected by a routine sequence-tagged-site PCR (STS-PCR) method. We then used the NGS method to detect CNV in the AZF region, including deletions and duplications.

RESULTS

A total of 326 samples from male infertility patients, family members, and sperm donors were studied between January 2011 and May 2017. AZF microdeletions were detected in 120 patients by STS-PCR, and these results were consistent with the results from NGS. In addition, of the 160 patients and male family members who had no microdeletions detected by STS-PCR, 51 cases were found to exhibit Y chromosome structural variations by the NGS method (31.88%, 51/160). No microdeletions were found in 46 donors by STS-PCR, but the NGS method revealed 11 of these donors (23.91%, 11/46) carried structural variations, which were mainly in the AZFc region, including partial deletions and duplications.

CONCLUSION

The established NGS method can replace the conventional STS-PCR method to detect Y chromosome microdeletions. The NGS method can detect CNV, such as partial deletion or duplication, and provide details of the abnormal range and size of variations.

A molecular marker of milk composition traits in NCAM2 gene of Chinese Holstein

The aim of this study was to detect the novel copy number variation (CNV) locus of NCAM2 gene in Chinese Holstein, and to analyze the effect of the novel CNV locus in NCAM2 gene on milk composition traits. The novel CNV locus of NCAM2 gene in 310 Chinese Holstein was detected by real-time quantitative fluorescent PCR (qPCR) and association analysis was performed between the novel CNV locus in NCAM2 gene and milk composition traits in Chinese Holstein. There are three CNV types of NCAM2 gene in Chinese Holstein: gain (increased copy number), median (normal copy number) and loss (deleted copy number). Statistical analysis revealed that there was a significant association between CNV types and milk fat rate (p < 0.05). Moreover, we also discovered that the milk production and milk protein rate of gain type is higher than that of loss type, but that of mediate type is lower than that of loss type. However, in terms of somatic cell score, loss type is higher than that of gain type, but that of mediate type is lower than that of gain type. These observations suggested that gain type can be used as a candidate molecular genetic marker of milk fat rate.HighlightsThe CNVs of the NCAM2 gene were detected and validated in Chinese Holstein.The type of CNV was successfully implemented using qPCR.The statistical analysis indicated that the CNV of the NCAM2 gene are significantly associated with milk fat rate.

Identification and utilization of copy number information for correcting Hi-C contact map of cancer cell lines

Background

Hi-C and its variant techniques have been developed to capture the spatial organization of chromatin. Normalization of Hi-C contact map is essential for accurate modeling and interpretation of high-throughput chromatin conformation capture (3C) experiments. Hi-C correction tools were originally developed to normalize systematic biases of karyotypically normal cell lines. However, a vast majority of available Hi-C datasets are derived from cancer cell lines that carry multi-level DNA copy number variations (CNVs). CNV regions display over- or under-representation of interaction frequencies compared to CN-neutral regions. Therefore, it is necessary to remove CNV-driven bias from chromatin interaction data of cancer cell lines to generate a euploid-equivalent contact map.

Results

We developed the HiCNAtra framework to compute high-resolution CNV profiles from Hi-C or 3C-seq data of cancer cell lines and to correct chromatin contact maps from systematic biases including CNV-associated bias. First, we introduce a novel ‘entire-fragment’ counting method for better estimation of the read depth (RD) signal from Hi-C reads that recapitulates the whole-genome sequencing (WGS)-derived coverage signal. Second, HiCNAtra employs a multimodal-based hierarchical CNV calling approach, which outperformed OneD and HiNT tools, to accurately identify CNVs of cancer cell lines. Third, incorporating CNV information with other systematic biases, HiCNAtra simultaneously estimates the contribution of each bias and explicitly corrects the interaction matrix using Poisson regression. HiCNAtra normalization abolishes CNV-induced artifacts from the contact map generating a heatmap with homogeneous signal. When benchmarked against OneD, CAIC, and ICE methods using MCF7 cancer cell line, HiCNAtra-corrected heatmap achieves the least 1D signal variation without deforming the inherent chromatin interaction signal. Additionally, HiCNAtra-corrected contact frequencies have minimum correlations with each of the systematic bias sources compared to OneD’s explicit method. Visual inspection of CNV profiles and contact maps of cancer cell lines reveals that HiCNAtra is the most robust Hi-C correction tool for ameliorating CNV-induced bias.

Conclusions

HiCNAtra is a Hi-C-based computational tool that provides an analytical and visualization framework for DNA copy number profiling and chromatin contact map correction of karyotypically abnormal cell lines. HiCNAtra is an open-source software implemented in MATLAB and is available at https://github.com/AISKhalil/HiCNAtra.

Exome sequencing identifies frequent genomic loss of TET1 in IDH-wild-type glioblastoma

Glioblastoma (GBM) is the most common and malignant brain tumor in adults. Genomic and epigenomic alterations of multiple cancer-driving genes are frequent in GBM. To identify molecular alterations associated with epigenetic aberrations, we performed whole exome sequencing-based analysis of DNA copy number variations in 55 adult patients with IDH-wild-type GBM. Beside mutations in common GBM driver genes such as TERTp (76%), TP53 (22%) and PTEN (20%), 67% of patients were affected by amplifications of genes associated with RTK/Rb/p53 cell signaling, including EGFR (45%), CDK4 (13%), and MDM2/4 (both 7%). The minimal deleted region at chromosome 10 was detected at the DNA demethylase TET1 (93%), mainly due to a loss-of-heterozygosity of complete chromosome 10 (53%) or by a mono-allelic microdeletion at 10q21.3 (7%). In addition, bi-allelic TET1 deletions, detected in 18 patients (33%), frequently co-occurred with EGFR amplification and were associated with low levels of TET1 mRNA expression, pointing at loss of TET1 activity. Bi-allelic TET1 loss was not associated with global concentrations of 5-hydroxymethylcytosine, indicating a site-specific effect of TET1 for DNA (de)methylation. Focal amplification of EGFR positively correlated with overall mutational burden, tumor size, and poor long-term survival. Bi-allelic TET1 loss was not an independent prognostic factor, but significantly associated with poor survival in patients with concomitant EGFR amplification. Rates of genomic TET1 deletion were significantly lower in a cohort of IDH1-mutated patients. Despite the relevance of TET1 for DNA demethylation and as potential therapeutic target, a frequent genomic loss of TET1 has not previously been reported in GBM.

Comparative analyses of copy number variations between Bos taurus and Bos indicus

Background

Bos taurus and Bos indicus are two main sub-species of cattle. However, the differential copy number variations (CNVs) between them are not yet well studied.

Results

Based on the new high-quality cattle reference genome ARS-UCD1.2, we identified 13,234 non-redundant CNV regions (CNVRs) from 73 animals of 10 cattle breeds (4 Bos taurus and 6 Bos indicus), by integrating three detection strategies. While 6990 CNVRs (52.82%) were shared by Bos taurus and Bos indicus, large CNV differences were discovered between them and these differences could be used to successfully separate animals into two subspecies. We found that 2212 and 538 genes uniquely overlapped with either indicine-specific CNVRs and or taurine-specific CNVRs, respectively. Based on F_ST, we detected 16 candidate lineage-differential CNV segments (top 0.1%) under selection, which overlapped with eight genes (CTNNA1, ENSBTAG00000004415, PKN2, BMPER, PDE1C, DNAJC18, MUSK, and PLCXD3). Moreover, we obtained 1.74 Mbp indicine-specific sequences, which could only be mapped on the Bos indicus reference genome UOA_Brahman_1. We found these sequences and their associated genes were related to heat resistance, lipid and ATP metabolic process, and muscle development under selection. We further analyzed and validated the top significant lineage-differential CNV. This CNV overlapped genes related to muscle cell differentiation, which might be generated from a retropseudogene of CTH but was deleted along Bos indicus lineage.

Conclusions

This study presents a genome wide CNV comparison between Bos taurus and Bos indicus. It supplied essential genome diversity information for understanding of adaptation and phenotype differences between the Bos taurus and Bos indicus populations.

Copy number variation analysis reveals variants associated with milk production traits in dairy goats

Copy number variation (CNV) is a major type of genomic structural variation. We investigated their impacts on goat dairy traits using the CaprineSNP50 array. From 120 samples of five dairy goat breeds, we totally identified 42 CNVs ranging from 56,044 bp to 4,337,625 bp. We found significant associations between two CNVs (CNV5 and CNV25) and two milk production traits (mean of milk fat yield and mean of milk protein yield) after false discovery rate (FDR) correction (P < 0.05). CNV5 overlaps the ADAMTS20 gene, which is involved in the differentiation of mammary cell and plays a crucial role in lactogenic activity of bovine mammary epithelial cells. CNV25 overlaps with PAPPA2, which has been found to be associated with bovine reproduction and milk production traits. Our results revealed that CNVs overlapped with ADAMTS20 and PAPPA2 could be involved in goat dairy traits and function as candidate markers for further genetic selection.

Mitochondrial dysfunction: A hidden trigger of autism?

Autism is a heterogeneous neurodevelopmental and neuropsychiatric disorder with no precise etiology. Deficits in cognitive functions uncover at early stages and are known to have an environmental and genetic basis. Since autism is multifaceted and also linked with other comorbidities associated with various organs, there is a possibility that there may be a fundamental cellular process responsible for this. These reasons place mitochondria at the point of interest as it is involved in multiple cellular processes predominantly involving metabolism. Mitochondria encoded genes were taken into consideration lately because it is inherited maternally, has its own genome and also functions the time of embryo development. Various researches have linked mitochondrial mishaps like oxidative stress, ROS production and mt-DNA copy number variations to autism. Despite dramatic advances in autism research worldwide, the studies focusing on mitochondrial dysfunction in autism is rather minimal, especially in India. India, owing to its rich diversity, may be able to contribute significantly to autism research. It is vital to urge more studies in this domain as it may help to completely understand the basics of the condition apart from a genetic standpoint. This review focuses on the worldwide and Indian scenario of autism research; mitochondrial abnormalities in autism and possible therapeutic approaches to combat it.

The MUC6/AP2A2 Locus and Its Relevance to Alzheimer's Disease: A Review

We recently reported evidence of Alzheimer's disease (AD)-linked genetic variation within the mucin 6 (MUC6) gene on chromosome 11p, nearby the adaptor-related protein complex 2 subunit alpha 2 (AP2A2) gene. This locus has interesting features related to human genomics and clinical research. MUC6 gene variants have been reported to potentially influence viral-including herpesvirus-immunity and the gut microbiome. Within the MUC6 gene is a unique variable number of tandem repeat (VNTR) region. We discovered an association between MUC6 VNTR repeat expansion and AD pathologic severity, particularly tau proteinopathy. Here, we review the relevant literature. The AD-linked VNTR polymorphism may also influence AP2A2 gene expression. AP2A2 encodes a polypeptide component of the adaptor protein complex, AP-2, which is involved in clathrin-coated vesicle function and was previously implicated in AD pathogenesis. To provide background information, we describe some key knowledge gaps in AD genetics research. The "missing/hidden heritability problem" of AD is highlighted. Extensive portions of the human genome, including the MUC6 VNTR, have not been thoroughly evaluated due to limitations of existing high-throughput sequencing technology. We present and discuss additional data, along with cautionary considerations, relevant to the hypothesis that MUC6 repeat expansion influences AD pathogenesis.

CNVmap: A Method and Software To Detect and Map Copy Number Variants from Segregation Data

Single nucleotide polymorphisms (SNPs) are used widely for detecting quantitative trait loci, or for searching for causal variants of diseases. Nevertheless, structural variations such as copy-number variants (CNVs) represent a large part of natural genetic diversity, and contribute significantly to trait variation. Numerous methods and softwares based on different technologies (amplicons, CGH, tiling, or SNP arrays, or sequencing) have already been developed to detect CNVs, but they bypass a wealth of information such as genotyping data from segregating populations, produced, e.g., for QTL mapping. Here, we propose an original method to both detect and genetically map CNVs using mapping panels. Specifically, we exploit the apparent heterozygous state of duplicated loci: peaks in appropriately defined genome-wide allelic profiles provide highly specific signatures that identify the nature and position of the CNVs. Our original method and software can detect and map automatically up to 33 different predefined types of CNVs based on segregation data only. We validate this approach on simulated and experimental biparental mapping panels in two maize populations and one wheat population. Most of the events found correspond to having just one extra copy in one of the parental lines, but the corresponding allelic value can be that of either parent. We also find cases with two or more additional copies, especially in wheat, where these copies locate to homeologues. More generally, our computational tool can be used to give additional value, at no cost, to many datasets produced over the past decade from genetic mapping panels.

Non-invasive prenatal testing reveals copy number variations related to pregnancy complications

Background

Pregnancy complications could lead to maternal and fetal morbidity and mortality. Early diagnosing and managing complications have been associated with good outcomes. The placenta was an important organ for development of pregnancy complications. Thus, non-invasive prenatal testing technologies could detect genetic variations, such as aneuploidies and sub-chromosomal copy number variations, reflecting defective placenta by maternal plasma cffDNAs. Maternal cffDNAs had been proved to derive from trophoblast cells of placenta.

Results

In order to find out the relationship between genetic variations and pregnancy complications, we reviewed NIPT results for subchromosomal copy number variations in a cohort of 3890 pregnancies without complications and 441 pregnancies with pregnancy complications including gestational diabetes mellitus (GDM), pregnancy-induced hypertension (PIH), preterm prelabor rupture of membranes (PPROM) and placenta implantation abnormalities (PIA). For GDMs, we identified three CNV regions containing some members of alpha- and beta-defensins, such as DEFA1, DEFA3, DEFB1. For PIHs, we found three duplication and one deletion region including Pcdhα, Pcdhβ, and Pcdhγ, known as protocadherins, which were complicated by hypertensive disorders. For PPROMs and PIAs, we identified one and two CNV regions, respectively. SFTPA2, SFTPD and SFTPA1, belonging to surfactant protein, was considered to moderated the inflammatory activation within the fetal extra-embryonic compartment, associated to duration of preterm prelabor rupture of fetal membranes, while MEF2C and TM6SF1 could be involved in trophoblast invasion and differentiation.

Conclusions

Our findings gave a clue to correlation between genetic variations of maternal cell-free DNAs and pregnancy complications.

Electronic supplementary material

The online version of this article (10.1186/s13039-019-0451-3) contains supplementary material, which is available to authorized users.

Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins

BACKGROUND

Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.

RESULTS

We identified 1043 CNV regions (CNVRs) from array comparative genomic hybridization (CGH) data of 47 Holstein bulls. Using a probe-based CNV association approach, we detected 87 CNVRs significantly (Bonferroni-corrected P value < 0.05) associated with at least one out of 41 complex traits. Within them, 39 CNVRs were simultaneously associated with at least 2 complex traits. Notably, 24 CNVRs were markedly related to daughter pregnancy rate (DPR). For example, CNVR661 containing CYP4A11 and CNVR213 containing CTR9, respectively, were associated with DPR and other traits related to reproduction, production, and body conformation. CNVR758 was also significantly related to DPR, with a nearby gene CAPZA3, encoding one of F-actin-capping proteins which play a role in determining sperm architecture and male fertility. We corroborated these CNVRs by examining their overlapped quantitative trait loci and comparing with previously published CNV results.

CONCLUSION

To our knowledge, this is one of the first genome-wide association studies based on CNVs called by array CGH in Holstein cattle. Our results contribute substantial information about the potential CNV impacts on reproduction, health, production, and body conformation traits, which lay the foundation for incorporating CNV into the future dairy cattle breeding program.

Single Color Multiplexed ddPCR Copy Number Measurements and Single Nucleotide Variant Genotyping

Droplet digital PCR (ddPCR) allows for accurate quantification of genetic events such as copy number variation and single nucleotide variants. Probe-based assays represent the current "gold-standard" for detection and quantification of these genetic events. Here, we introduce a cost-effective single color ddPCR assay that allows for single genome resolution quantification of copy number and single nucleotide variation.

Analysis of copy number variation in the NEDD4L gene potentially implicated in body height in the Japanese population

Recently it has been recognized that a considerable number of copy number variations (CNVs) are associated with diseases and other complex human traits. In our previous study, we developed a simple quantitative real-time PCR (Q-PCR) method for analysis of CNV copy number, which had the advantage of obviating the need for reference DNA with a known copy number. Using DNA samples obtained from 231 Japanese individuals, we applied this method for analyzing the copy number of a candidate CNV associated with body height, located in the neural precursor cell expressed, developmentally down-regulated 4-like, E3 ubiquitin protein ligase (NEDD4L) gene. In addition, the appropriateness of the results was evaluated and confirmed by quantification of amplicons with an Agilent 2100 Bioanalyzer. The NEDD4L gene encodes a member of the Nedd4 family of HECT domain E3 ubiquitin ligases. The target CNV located in the intron has been found to be significantly associated with height variation in Chinese. However, it remains unknown whether such an association exists in other populations, including Japanese. Analysis of the correlations between copy number and body height using ANOVA revealed no statistically significant correlations in Japanese.

Genomic alterations in gastric cancers discovered via whole-exome sequencing

Background

Gastric cancer (GC) ranks the second in mortality rate among all cancers. Metastases account for most of the deaths in GC patients. Yet our understanding of GC and its metastasis mechanism is still very limited.

Methods

We performed 20 whole-exome sequencing (WES) on 5 typical metastatic gastric adenocarcinoma (GAC) patients with lymph node metastasis. We compared both the primary tumors to their metastatic lymph nodes, and a specific analysis pipeline was used to detect single nucleotide variants (SNVs), small insertions/deletions (indels) and copy number variants (CNVs).

Results

(1) We confirmed 30 candidate mutations in both primary and lymph nodes tissues, and other 7 only in primary tumors. (2) Copy number gains were observed in a large section of 17q12–21, as well as copy number losses in regions containing CDKN2A and CDKN2B in both primary and lymph nodes tissues.

Conclusions

Our results provide preliminary insights in the molecular mechanisms of GC initiation, development, and metastatic progression. These results need to be validated through large-scale studies.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-5097-8) contains supplementary material, which is available to authorized users.

Chromosomal microarray analysis of patients with Duane retraction syndrome

PURPOSE

Duane retraction syndrome (DS) is a rare congenital strabismus with genetic heterogeneity. The genetic causes of DS are not always of monogenic origin; various chromosomal copy number variations (CNVs) have also been reported. The objective of our study was to characterize the CNVs, including gains and losses detected by high-resolution chromosomal microarray in patients with DS.

METHODS

Twenty patients with DS were investigated using high-resolution chromosomal microarray analysis (CMA) (Affymetrix CytoScan Array 750 K). Conventional cytogenetic analysis was also performed.

RESULTS

All samples revealed normal karyotype by cytogenetic analysis. However, in all our patients, multiple CNVs, including gains and losses, were detected using the high-resolution CMA method. Chromosomal loci 1q21.2, 2p11.2-q11.1, 2q21.1-q21.2, 4p16.1, 7p11.2-q11.21, 14q32.33, 17p11.2-q11.1 and 20p11.1-q11.21 were the most frequently affected regions.

CONCLUSIONS

This study emphasized that CNVs in several chromosomal regions are known to be involved in DS. We also underscore the genetic heterogeneity of DS. Our suggestion is that genes located in the most frequently affected regions should be focused on in the following candidate gene studies.

Clinical analysis of germline copy number variation in DMD using a non-conjugate hierarchical Bayesian model

Background

Detection of copy number variants (CNVs) is an important aspect of clinical testing for several disorders, including Duchenne muscular dystrophy, and is often performed using multiplex ligation-dependent probe amplification (MLPA). However, since many genetic carrier screens depend instead on next-generation sequencing (NGS) for wider discovery of small variants, they often do not include CNV analysis. Moreover, most computational techniques developed to detect CNVs from exome sequencing data are not suitable for carrier screening, as they require matched normals, very large cohorts, or extensive gene panels.

Methods

We present a computational software package, geneCNV (http://github.com/vkozareva/geneCNV), which can identify exon-level CNVs using exome sequencing data from only a few genes. The tool relies on a hierarchical parametric model trained on a small cohort of reference samples.

Results

Using geneCNV, we accurately inferred heterozygous CNVs in the DMD gene across a cohort of 15 test subjects. These results were validated against MLPA, the current standard for clinical CNV analysis in DMD. We also benchmarked the tool’s performance against other computational techniques and found comparable or improved CNV detection in DMD using data from panels ranging from 4,000 genes to as few as 8 genes.

Conclusions

geneCNV allows for the creation of cost-effective screening panels by allowing NGS sequencing approaches to generate results equivalent to bespoke genotyping assays like MLPA. By using a parametric model to detect CNVs, it also fulfills regulatory requirements to define a reference range for a genetic test. It is freely available and can be incorporated into any Illumina sequencing pipeline to create clinical assays for detection of exon duplications and deletions.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0404-4) contains supplementary material, which is available to authorized users.

Genetic and epigenetic changes in clonal descendants of irradiated human fibroblasts

To study delayed genetic and epigenetic radiation effects, which may trigger radiation-induced carcinogenesis, we have established single-cell clones from irradiated and non-irradiated primary human fibroblasts. Stable clones were endowed with the same karyotype in all analyzed metaphases after 20 population doublings (PDs), whereas unstable clones displayed mosaics of normal and abnormal karyotypes. To account for variation in radiation sensitivity, all experiments were performed with two different fibroblast strains. After a single X-ray dose of 2 Gy more than half of the irradiated clones exhibited radiation-induced genome instability (RIGI). Irradiated clones displayed an increased rate of loss of chromosome Y (LOY) and copy number variations (CNVs), compared to controls. CNV breakpoints clustered in specific chromosome regions, in particular 3p14.2 and 7q11.21, coinciding with common fragile sites. CNVs affecting the FHIT gene in FRA3B were observed in independent unstable clones and may drive RIGI. Bisulfite pyrosequencing of control clones and the respective primary culture revealed global hypomethylation of ALU, LINE-1, and alpha-satellite repeats as well as rDNA hypermethylation during in vitro ageing. Irradiated clones showed further reduced ALU and alpha-satellite methylation and increased rDNA methylation, compared to controls. Methylation arrays identified several hundred differentially methylated genes and several enriched pathways associated with in vitro ageing. Methylation changes in 259 genes and the MAP kinase signaling pathway were associated with delayed radiation effects (after 20 PDs). Collectively, our results suggest that both genetic (LOY and CNVs) and epigenetic changes occur in the progeny of exposed cells that were not damaged directly by irradiation, likely contributing to radiation-induced carcinogenesis. We did not observe epigenetic differences between stable and unstable irradiated clones. The fact that the DNA methylation (DNAm) age of clones derived from the same primary culture varied greatly suggests that DNAm age of a single cell (represented by a clone) can be quite different from the DNAm age of a tissue. We propose that DNAm age reflects the emergent property of a large number of individual cells whose respective DNAm ages can be highly variable.

Genome-wide copy number variant analysis reveals variants associated with 10 diverse production traits in Holstein cattle

BACKGROUND

Copy number variation (CNV) is an important type of genetic variation contributing to phenotypic differences among mammals and may serve as an alternative molecular marker to single nucleotide polymorphism (SNP) for genome-wide association study (GWAS). Recently, GWAS analysis using CNV has been applied in livestock, although few studies have focused on Holstein cattle.

RESULTS

We describe 191 CNV detected using intensity data from over 700,000 SNP genotypes generated with the BovineHD Genotyping BeadChip (Illumina, San Diego, CA) in 528 Holstein cows. The CNV were used for GWAS analysis of 10 important production traits of 473 cattle related to feed intake, milk quality, and female fertility, as well as 2 composite traits of net merit and productive life. In total, we detected 57 CNV associated (P < 0.05 after false discovery rate correction) with at least one of the 10 phenotypes. Focusing on feed efficiency and intake-related phenotypes of residual feed intake and dry matter intake, we detected a single CNV associated with both traits which overlaps a predicted olfactory receptor gene OR2A2 (LOC787786). Additionally, 2 CNV within the RXFP4 (relaxin/insulin like family peptide receptor 4) and 2 additional olfactory receptor gene regions, respectively, were associated with residual feed intake. The RXFP4 gene encodes a receptor for an orexigenic peptide, insulin-like peptide 5 produced by intestinal L cells, which is expressed by enteric neurons. Olfactory receptors are critical for transmitting the effects of odorants, contributing to the sense of smell, and have been implicated in participating in appetite regulation.

CONCLUSIONS

Our results identify CNV for genomic evaluation in Holstein cattle, and provide candidate genes, such as RXFP4, contributing to variation in feed efficiency and feed intake-related traits. These results indicate potential novel targets for manipulating feed intake-related traits of livestock.

Gene expression patterns of chicken neuregulin 3 in association with copy number variation and frameshift deletion

Background

Neuregulin 3 (NRG3) plays a key role in central nervous system development and is a strong candidate for human mental disorders. Thus, genetic variation in NRG3 may have some impact on a variety of phenotypes in non-mammalian vertebrates. Recently, genome-wide screening for short insertions and deletions in chicken (Gallus gallus) genomes has provided useful information about structural variation in functionally important genes. NRG3 is one such gene that has a putative frameshift deletion in exon 2, resulting in premature termination of translation. Our aims were to characterize the structure of chicken NRG3 and to compare expression patterns between NRG3 isoforms.

Results

Depending on the presence or absence of the 2-bp deletion in chicken NRG3, 3 breeds (red junglefowl [RJF], Boris Brown [BB], and Hinai-jidori [HJ]) were genotyped using flanking primers. In the commercial breeds (BB and HJ), approximately 45% of individuals had at least one exon 2 allele with the 2-bp deletion, whereas there was no deletion allele in RJF. The lack of a homozygous mutant indicated the existence of duplicated NRG3 segments in the chicken genome. Indeed, highly conserved elements consisting of exon 1, intron 1, exon 2, and part of intron 2 were found in the reference RJF genome, and quantitative PCR detected copy number variation (CNV) between breeds as well as between individuals. The copy number of conserved elements was significantly higher in chicks harboring the 2-bp deletion in exon 2. We identified 7 novel transcript variants using total mRNA isolated from the amygdala. Novel isoforms were found to lack the exon 2 cassette, which probably harbored the premature termination codon. The relative transcription levels of the newly identified isoforms were almost the same between chick groups with and without the 2-bp deletion, while chicks with the deletion showed significant suppression of the expression of previously reported isoforms.

Conclusions

A putative frameshift deletion and CNV in chicken NRG3 are structural mutations that occurred before the establishment of commercial chicken lines. Our results further suggest that the putative frameshift deletion in exon 2 may potentially affect the expression level of particular isoforms of chicken NRG3.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0537-z) contains supplementary material, which is available to authorized users.

Copy-number variation of housekeeping gene rpl13a in rat strains selected for nervous system excitability

We evaluated copy number variation (CNV) for four genes in rat strains differing in nervous system excitability. rpl13a copy number is significantly reduced in hippocampus and bone marrow in rats with a high excitability threshold and stress. The observed phenomenon may be associated with a role for rpl13a in lipid metabolism.

Simple screening method for copy number variations associated with physical features

Recent studies of copy number variations (CNVs) associated with physical features, such as body mass index, body height or bone length, have suggested that such CNVs could serve as markers in forensic cases involving unidentified individuals. However, the process of cataloging CNVs has been slow because of the cumbersome nature and low reliability of the procedures involved. Here we describe a simple quantitative real-time PCR (Q-PCR) method for screening of medicolegally useful CNVs, which does not require reference DNA with known copy number. The first step is to prepare a chimeric plasmid vector including one copy each of the single-copy gene-specific sequence as the internal standard, and the target CNV-specific sequence. To assess the validity of this new method, we analyzed CNVs in the LTBP1 and ETV6 gene regions, both of which are candidate CNVs associated with body height. The PCR efficiencies for the single-copy (reference) gene and the target CNV were similar, indicating that quantitation was reliable. Furthermore, simulated analysis of the LTBP1 CNV using mock samples prepared by mixing vectors in varying proportions showed that this analytical method allowed correct determination of the LTBP1 copy number. These results demonstrated that our simple method has considerable potential for screening of trait-related CNVs that would be useful for forensic casework.

Array-Based Comparative Genomic Hybridization (aCGH)

Copy number variations (CNVs) in the genomes have been suggested to play important roles in human evolution, genetic diversity, and disease susceptibility. A number of assays have been developed for the detection of CNVs, including fluorescent in situ hybridization (FISH), array-based comparative genomic hybridization (aCGH), PCR-based assays, and next-generation sequencing (NGS). In this chapter, we describe a microarray method that has been used for the detection of genome-wide CNVs, loss of heterozygosity (LOH), and uniparental disomy (UPD) associated with constitutional and neoplastic disorders.

Analyzing Copy Number Variation Using Pulsed-Field Gel Electrophoresis: Providing a Genetic Diagnosis for FSHD1

The myopathy facioscapulohumeral muscular dystrophy type 1 (FSHD1) is caused by copy number variation of the D4Z4 macrosatellite repeat on chromosome 4. In unaffected individuals the number of 3.3 kb D4Z4 units varies between 8 and 100, whereas 1-10 units are seen in FSHD1 cases. A homologous and heterogenous D4Z4 array can be found on chromosome 10q, but contractions of this array are typically not associated with FSHD. Discriminating between the chromosome 4 and chromosome 10 D4Z4 arrays, as well as determining the array size, requires the use of pulsed-field gel electrophoresis, Southern blotting, and the isolation of high-quality DNA.

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.

Genome-wide CNV analysis reveals variants associated with growth traits in Bos indicus

BACKGROUND

Apart from single nucleotide polymorphism (SNP), copy number variation (CNV) is another important type of genetic variation, which may affect growth traits and play key roles for the production of beef cattle. To date, no genome-wide association study (GWAS) for CNV and body traits in beef cattle has been reported, so the present study aimed to investigate this type of association in one of the most important cattle subspecies: Bos indicus (Nellore breed).

RESULTS

We have used intensity data from over 700,000 SNP probes across the bovine genome to detect common CNVs in a sample of 2230 Nellore cattle, and performed GWAS between the detected CNVs and nine growth traits. After filtering for frequency and length, a total of 231 CNVs ranging from 894 bp to 4,855,088 bp were kept and tested as predictors for each growth trait using linear regression analysis with principal components correction. There were 49 significant associations identified among 17 CNVs and seven body traits after false discovery rate correction (P < 0.05). Among the 17 CNVs, three were significant or marginally significant for all the traits. We have compared the locations of associated CNVs with quantitative trait locus and the RefGene database, and found two sets of 9 CNVs overlapping with either known QTLs or genes, respectively. The gene overlapping with CNV100, KCNJ12, is a functional candidate for muscle development and plays critical roles in muscling traits.

CONCLUSION

This study presents the first CNV-based GWAS of growth traits using high density SNP microarray data in cattle. We detected 17 CNVs significantly associated with seven growth traits and one of them (CNV100) may be involved in growth traits through KCNJ12.

Multiallelic copy number variation in the complement component 4A (C4A) gene is associated with late-stage age-related macular degeneration (AMD)

BACKGROUND

Age-related macular degeneration (AMD) is the leading cause of vision loss in Western societies with a strong genetic component. Candidate gene studies as well as genome-wide association studies strongly implicated genetic variations in complement genes to be involved in disease risk. So far, no association of AMD with complement component 4 (C4) was reported probably due to the complex nature of the C4 locus on chromosome 6.

METHODS

We used multiplex ligation-dependent probe amplification (MLPA) to determine the copy number of the C4 gene as well as of both relevant isoforms, C4A and C4B, and assessed their association with AMD using logistic regression models.

RESULTS

Here, we report on the analysis of 2645 individuals (1536 probands and 1109 unaffected controls), across three different centers, for multiallelic copy number variation (CNV) at the C4 locus. We find strong statistical significance for association of increased copy number of C4A (OR 0.81 (0.73; 0.89);P = 4.4 × 10(-5)), with the effect most pronounced in individuals over 78 years (OR 0.67 (0.55; 0.81)) and females (OR 0.77 (0.68; 0.87)). Furthermore, this association is independent of known AMD-associated risk variants in the nearby CFB/C2 locus, particularly in females and in individuals over 78 years.

CONCLUSIONS

Our data strengthen the notion that complement dysregulation plays a crucial role in AMD etiology, an important finding for early intervention strategies and future therapeutics. In addition, for the first time, we provide evidence that multiallelic CNVs are associated with AMD pathology.

Array-based molecular karyotyping in fetuses with isolated brain malformations identifies disease-causing CNVs

Background

The overall birth prevalence for congenital malformations of the central nervous system (CNS) among Europeans may be as high as 1 in 100 live births. The etiological factors remain largely unknown. The aim of this study was to detect causative copy number variations (CNVs) in fetuses of terminated pregnancies with prenatally detected isolated brain malformations.

Methods

Array-based molecular karyotyping was performed in a cohort of 35 terminated fetuses with isolated CNS malformations. Identified putative disease-causing CNVs were confirmed using quantitative polymerase chain reaction or multiplex ligation-dependent probe amplification.

Results

Based on their de novo occurrence and/or their established association with congenital brain malformations, we detected five disease-causing CNVs in four fetuses involving chromosomal regions 6p25.1-6p25.3 (FOXC1), 6q27, 16p12.3, Xp22.2-Xp22.32 (MID1), and Xp22.32-Xp22.33. Furthermore, we detected a probably disease-causing CNV involving chromosomal region 3p26.3 in one fetus, and in addition, we detected 12 CNVs in nine fetuses of unknown clinical significance. All CNVs except for two were absent in 1307 healthy in-house controls (frequency <0.0008). Each of the two CNVs present in in-house controls was present only once (frequency = 0.0008). Furthermore, our data suggests the involvement of CNTN6 and KLHL15 in the etiology of agenesis of the corpus callosum, the involvement of RASD1 and PTPRD in Dandy-Walker malformation, and the involvement of ERMARD in ventriculomegaly.

Conclusions

Our study suggests that CNVs play an important role in the etiology of isolated brain malformations.