A computational framework discovers new copy number variants with functional importance

Copy Number Variations and Gene Mutations Identified by Multiplex Ligation-Dependent Probe Amplification in Romanian Chronic Lymphocytic Leukemia Patients

Chronic lymphocytic leukemia (CLL) is known for its wide-ranging clinical and genetic diversity. The study aimed to assess the associations between copy number variations (CNVs) and various biological and clinical features, as well as the survival rates of CLL patients and to evaluate the effectiveness of the multiplex ligation-dependent probe amplification (MLPA) technique in CLL patients.DNA was extracted from 110 patients, and MLPA was performed. Mutations in NOTCH1, SF3B1, and MYD88 were also analyzed. A total of 52 patients showed at least one CNV, 26 had at least one somatic mutation, and 10 presented both, CNVs, and somatic mutations. The most commonly identified CNVs were del(114.3), del(11q22.3), and dup(12q23.2). Other CNVs identified included del(17p13.1), del(14q32.33), dup(10q23.31), and del(19p13.2). One patient was identified with concomitant trisomy 12, 13, and 19. NOTCH1 and SF3B1 mutations were found in 13 patients each, either alone or in combination with other mutations or CNVs, while MYD88 mutation was identified in one patient. Forty-two patients had normal results. Associations between the investigated CNVs and gene mutations and patients' overall survival were found. The presence of NOTCH1 and SF3B1 mutations or the combination of NOTCH1 mutation and CNVs significantly influenced the survival of patients with CLL. Both mutations are frequently associated with different CNVs. Del(13q) is associated with the longest survival rate, while the shortest survival is found in patients with del(17p). Even if MLPA has constraints, it may be used as the primary routine analysis in patients with CLL.

A Mild PUM1 Mutation Is Associated with Adult-Onset Ataxia, whereas Haploinsufficiency Causes Developmental Delay and Seizures

Certain mutations can cause proteins to accumulate in neurons, leading to neurodegeneration. We recently showed, however, that upregulation of a wild-type protein, Ataxin1, caused by haploinsufficiency of its repressor, the RNA-binding protein Pumilio1 (PUM1), also causes neurodegeneration in mice. We therefore searched for human patients with PUM1 mutations. We identified eleven individuals with either PUM1 deletions or de novo missense variants who suffer a developmental syndrome (Pumilio1-associated developmental disability, ataxia, and seizure; PADDAS). We also identified a milder missense mutation in a family with adult-onset ataxia with incomplete penetrance (Pumilio1-related cerebellar ataxia, PRCA). Studies in patient-derived cells revealed that the missense mutations reduced PUM1 protein levels by ∼25% in the adult-onset cases and by ∼50% in the infantile-onset cases; levels of known PUM1 targets increased accordingly. Changes in protein levels thus track with phenotypic severity, and identifying posttranscriptional modulators of protein expression should identify new candidate disease genes.

A multiplex ligation‑dependent probe amplification‑based next‑generation sequencing approach for the detection of copy number variations in the human genome

The aim of the present study was to describe a multiplex ligation‑dependent probe amplification (MLPA)‑based next‑generation sequencing (NGS) assay that exhibited a significantly higher efficiency in detecting copy number variations (CNVs) and known single‑nucleotide variants, compared with traditional MLPA. MLPA polymerase chain reaction products were used to construct a library with indexed adapters, which was subsequently tested on an NGS platform, and the resulting data were analyzed by a series of analytical software. The reads from each probe reflected genetic variations in the target regions, and fragment differentiation was based on the specific base composition of the sequences, rather than fragment length, which was determined by capillary electrophoresis. The results of this approach were not only consistent with the MLPA results following capillary electrophoresis, but also coincided with the CNV results from the single‑nucleotide polymorphism array chip. This method allowed high‑throughput screening for the number of fragments and samples by integrating additional indices for detection. Furthermore, this technology precisely and accurately performed large‑scale detection and quantification of DNA variations, thereby serving as an effective and sensitive method for diagnosing genetic disorders caused by CNVs and known single‑nucleotide variations. Notably, MLPA‑NGS circumvents the problems associated with the inaccuracies of NGS in CNV detection due to the use of target sequence capture.

Chromosome 5q33 deletions associated with congenital heart defects

Congenital heart defects (CHD) are present in over 1% of all newborns and are the leading cause of birth-defect-related deaths in the United States. We describe two male subjects with CHD, one with an atrial septal defect, a ventricular septal defect, and pulmonary artery stenosis; and the other with tetralogy of Fallot and a right aortic arch, who carry partially overlapping, de novo deletions of chromosome 5q33. The maximum region of overlap between these deletions encompasses HAND1 and SAP30L, two genes that have previously been shown to play a role in cardiac development. HAND1 encodes a basic helix-loop-helix transcription factor. Cardiac-specific ablation of Hand1 in mice causes septal, valvular, and outflow tract defects. SAP30L, its paralog SAP30, and other SAP proteins form part of a multi-subunit complex involved in transcriptional regulation via histone deacetylation. Morpholino knockdown of sap30L in zebrafish, which do not have a distinct sap30 gene, leads to cardiac hypoplasia and cardiac insufficiency. We subsequently identified two other individuals with chromosomal deletions involving HAND1 and SAP30L in whom cardiac-related medical problems were not described. These observations suggest that haploinsufficiency of HAND1 and/or SAP30L may contribute to the development of CHD, although the contribution of other genes on chromosome 5q33 cannot be excluded. Our findings also suggest that the penetrance of CHD associated with 5q33 deletions is incomplete and may be influenced by other genetic, environmental or stochastic factors. © 2016 Wiley Periodicals, Inc.

NUDT21-spanning CNVs lead to neuropsychiatric disease and altered MeCP2 abundance via alternative polyadenylation

The brain is sensitive to the dose of MeCP2 such that small fluctuations in protein quantity lead to neuropsychiatric disease. Despite the importance of MeCP2 levels to brain function, little is known about its regulation. In this study, we report eleven individuals with neuropsychiatric disease and copy-number variations spanning NUDT21, which encodes a subunit of pre-mRNA cleavage factor Im. Investigations of MECP2 mRNA and protein abundance in patient-derived lymphoblastoid cells from one NUDT21 deletion and three duplication cases show that NUDT21 regulates MeCP2 protein quantity. Elevated NUDT21 increases usage of the distal polyadenylation site in the MECP2 3′ UTR, resulting in an enrichment of inefficiently translated long mRNA isoforms. Furthermore, normalization of NUDT21 via siRNA-mediated knockdown in duplication patient lymphoblasts restores MeCP2 to normal levels. Ultimately, we identify NUDT21 as a novel candidate for intellectual disability and neuropsychiatric disease, and elucidate a mechanism of pathogenesis by MeCP2 dysregulation via altered alternative polyadenylation.

DOI:http://dx.doi.org/10.7554/eLife.10782.001

The X-chromosome carries a number of genes that are involved in a child's intellectual development. One of these genes encodes a protein called MeCP2, which is important for brain function after birth. Mutations in the MECP2 gene cause a disorder known as Rett syndrome. At around 18 months of age, affected children begin to lose the cognitive and motor skills that they had previously acquired. Individuals with extra copies of this gene also show cognitive impairments. For both diseases, individuals with levels of the MeCP2 protein that are the most different from those found in healthy individuals also show the most severe symptoms.

To produce the protein that is encoded by a particular gene, enzymes inside the cell must first make a copy of that gene using a molecule called messenger ribonucleic acid (or mRNA). This mRNA is then used as a template to assemble the protein itself. In the case of MECP2, two different mRNA templates are produced: a long version and a short version. A gene called NUDT21 makes a protein that regulates whether the long or short version of MECP2 mRNA is made.

Gennarino, Alcott et al. have now discovered that people with too many, or too few, copies of the NUDT21 gene have intellectual disabilities and altered levels of MeCP2 protein. Specifically, individuals with extra copies of NUDT21—and thus higher levels of the corresponding protein—produce more of the long MECP2 mRNA. The production of proteins from this long mRNA is less efficient than from the short mRNA; therefore, these individuals have lower levels of MeCP2 protein. The opposite is true for individuals who lack a copy of the NUDT21 gene.

To confirm these data, Gennarino, Alcott et al. grew cells in the laboratory from patients with extra copies of the NUDT21 gene and found that reducing the production of its protein returned the levels of the MeCP2 protein back to normal. These findings show that alterations in the NUDT21 gene cause changes in the level of MeCP2 protein in cells and leads to neuropsychiatric diseases.

DOI:http://dx.doi.org/10.7554/eLife.10782.002

In-silico identification and functional validation of allele-dependent AR enhancers

Androgen Receptor (AR) and Estrogen Receptors (ERs) are key nuclear receptors that can cooperate in orchestrating gene expression programs in multiple tissues and diseases, targeting binding elements in promoters and distant enhancers. We report the unbiased identification of enhancer elements bound by AR and ER-α whose activity can be allele-specific depending on the status of nearby Single Nucleotide Polymorphisms (SNP). ENCODE data were computationally mined to nominate genomic loci with: (i) chromatin signature of enhancer activity from activation histone marks, (ii) binding evidence by AR and ER-α, (iii) presence of a SNP. Forty-one loci were identified and two, on 1q21.3 and 13q34, selected for characterization by gene reporter, Chromatin immunoprecipitation (ChIP) and RT-qPCR assays in breast (MCF7) and prostate (PC-3) cancer-derived cell lines. We observed allele-specific enhancer activity, responsiveness to ligand-bound AR, and potentially influence on the transcription of closely located genes (RAB20, ING1, ARHGEF7, ADAM15). The 1q21.3 variant, rs2242193, showed impact on AR binding in MCF7 cells that are heterozygous for the SNP. Our unbiased genome-wide search proved to be an efficient methodology to discover new functional polymorphic regulatory regions (PRR) potentially acting as risk modifiers in hormone-driven cancers and overall nominated SNPs in PRR across 136 transcription factors.

Genetic predisposition to prostate cancer: Update and future perspectives

OBJECTIVE

Prostate cancer is the second most frequent cancer in men worldwide and kills over 250,000 men worldwide every year. Prostate cancer is a heterogeneous disease at the clinical and the molecular level. The Scandinavian Twin Registry Study demonstrated that in contrast to most malignancies where environment was the overriding influence, heritable factors account for more than fifty percent of prostate cancers.

METHODS AND MATERIALS

We review the literature on prostate cancer risk variants (rare and common) including SNPs and Copy Number Variants (CNVs) and discuss the potential implications of significant variants for prostate cancer patient care.

RESULTS

The search for prostate cancer susceptibility genes has included both family-based studies and case-control studies utilizing a variety of approaches from array-based to sequencing-based studies. A major challenge is to identify genetic variants associated with more aggressive, potentially lethal prostate cancer and to understand their role in the progression of the disease.

CONCLUSION

Future risk models useful in the clinical setting will likely incorporate several risk loci rather than single variants and may be dependent on an individual patient's ethnic background.

Copy number variation detection using next generation sequencing read counts

Background

A copy number variation (CNV) is a difference between genotypes in the number of copies of a genomic region. Next generation sequencing (NGS) technologies provide sensitive and accurate tools for detecting genomic variations that include CNVs. However, statistical approaches for CNV identification using NGS are limited. We propose a new methodology for detecting CNVs using NGS data. This method (henceforth denoted by m-HMM) is based on a hidden Markov model with emission probabilities that are governed by mixture distributions. We use the Expectation-Maximization (EM) algorithm to estimate the parameters in the model.

Results

A simulation study demonstrates that our proposed m-HMM approach has greater power for detecting copy number gains and losses relative to existing methods. Furthermore, application of our m-HMM to DNA sequencing data from the two maize inbred lines B73 and Mo17 to identify CNVs that may play a role in creating phenotypic differences between these inbred lines provides results concordant with previous array-based efforts to identify CNVs.

Conclusions

The new m-HMM method is a powerful and practical approach for identifying CNVs from NGS data.

Short copy number variations potentially associated with tonic immobility responses in newly hatched chicks

Introduction

Tonic immobility (TI) is fear-induced freezing that animals may undergo when confronted by a threat. It is principally observed in prey species as defence mechanisms. In our preliminary research, we detected large inter-individual variations in the frequency and duration of freezing behavior among newly hatched domestic chicks (Gallus gallus). In this study we aim to identify the copy number variations (CNVs) in the genome of chicks as genetic candidates that underlie the behavioral plasticity to fearful stimuli.

Methods

A total of 110 domestic chicks were used for an association study between TI responses and copy number polymorphisms. Array comparative genomic hybridization (aCGH) was conducted between chicks with high and low TI scores using an Agilent 4×180 custom microarray. We specifically focused on 3 genomic regions (>60 Mb) of chromosome 1 where previous quantitative trait loci (QTL) analysis showed significant F-values for fearful responses.

Results

ACGH successfully detected short CNVs within the regions overlapping 3 QTL peaks. Eleven of these identified loci were validated by real-time quantitative polymerase chain reaction (qPCR) as copy number polymorphisms. Although there wkas no significant p value in the correlation analysis between TI scores and the relative copy number within each breed, several CNV loci showed significant differences in the relative copy number between 2 breeds of chicken (White Leghorn and Nagoya) which had different quantitative characteristics of fear-induced responses.

Conclusion

Our data shows the potential CNVs that may be responsible for innate fear response in domestic chicks.

Variants at IRX4 as prostate cancer expression quantitative trait loci

Genome-wide association studies (GWAS) have identified numerous prostate cancer-associated risk loci. Some variants at these loci may be regulatory and influence expression of nearby genes. Such loci are known as cis-expression quantitative trait loci (cis-eQTL). As cis-eQTLs are highly tissue-specific, we asked if GWAS-identified prostate cancer risk loci are cis-eQTLs in human prostate tumor tissues. We investigated 50 prostate cancer samples for their genotype at 59 prostate cancer risk-associated single-nucleotide polymorphisms (SNPs) and performed cis-eQTL analysis of transcripts from paired primary tumors within two megabase windows. We tested 586 transcript-genotype associations, of which 27 were significant (false discovery rate ≤10%). An equivalent eQTL analysis of the same prostate cancer risk loci in lymphoblastoid cell lines did not result in any significant associations. The top-ranked cis-eQTL involved the IRX4 (Iroquois homeobox protein 4) transcript and rs12653946, tagged by rs10866528 in our study (P=4.91 × 10(-5)). Replication studies, linkage disequilibrium, and imputation analyses highlight population specificity at this locus. We independently validated IRX4 as a potential prostate cancer risk gene through cis-eQTL analysis of prostate cancer risk variants. Cis-eQTL analysis in relevant tissues, even with a small sample size, can be a powerful method to expedite functional follow-up of GWAS.

Impact of constitutional copy number variants on biological pathway evolution

Background

Inherited Copy Number Variants (CNVs) can modulate the expression levels of individual genes. However, little is known about how CNVs alter biological pathways and how this varies across different populations. To trace potential evolutionary changes of well-described biological pathways, we jointly queried the genomes and the transcriptomes of a collection of individuals with Caucasian, Asian or Yoruban descent combining high-resolution array and sequencing data.

Results

We implemented an enrichment analysis of pathways accounting for CNVs and genes sizes and detected significant enrichment not only in signal transduction and extracellular biological processes, but also in metabolism pathways. Upon the estimation of CNV population differentiation (CNVs with different polymorphism frequencies across populations), we evaluated that 22% of the pathways contain at least one gene that is proximal to a CNV (CNV-gene pair) that shows significant population differentiation. The majority of these CNV-gene pairs belong to signal transduction pathways and 6% of the CNV-gene pairs show statistical association between the copy number states and the transcript levels.

Conclusions

The analysis suggested possible examples of positive selection within individual populations including NF-kB, MAPK signaling pathways, and Alu/L1 retrotransposition factors. Altogether, our results suggest that constitutional CNVs may modulate subtle pathway changes through specific pathway enzymes, which may become fixed in some populations.

Identification of functionally active, low frequency copy number variants at 15q21.3 and 12q21.31 associated with prostate cancer risk

Copy number variants (CNVs) are a recently recognized class of human germ line polymorphisms and are associated with a variety of human diseases, including cancer. Because of the strong genetic influence on prostate cancer, we sought to identify functionally active CNVs associated with susceptibility of this cancer type. We queried low-frequency biallelic CNVs from 1,903 men of Caucasian origin enrolled in the Tyrol Prostate Specific Antigen Screening Cohort and discovered two CNVs strongly associated with prostate cancer risk. The first risk locus (P = 7.7 × 10(-4), odds ratio = 2.78) maps to 15q21.3 and overlaps a noncoding enhancer element that contains multiple activator protein 1 (AP-1) transcription factor binding sites. Chromosome conformation capture (Hi-C) data suggested direct cis-interactions with distant genes. The second risk locus (P = 2.6 × 10(-3), odds ratio = 4.8) maps to the α-1,3-mannosyl-glycoprotein 4-β-N-acetylglucosaminyltransferase C (MGAT4C) gene on 12q21.31. In vitro cell-line assays found this gene to significantly modulate cell proliferation and migration in both benign and cancer prostate cells. Furthermore, MGAT4C was significantly overexpressed in metastatic versus localized prostate cancer. These two risk associations were replicated in an independent PSA-screened cohort of 800 men (15q21.3, combined P = 0.006; 12q21.31, combined P = 0.026). These findings establish noncoding and coding germ line CNVs as significant risk factors for prostate cancer susceptibility and implicate their role in disease development and progression.

Plasticity of the myelination genomic fabric

This study aimed to quantify the influence of the astrocyte proximity on myelination genomic fabric (MYE) of oligodendrocytes, defined as the most interconnected and stably expressed gene web responsible for myelination. Such quantitation is important to evaluate whether astrocyte signaling may contribute to demyelination when impaired and remyelination when properly restored. For this, we compared changes in the gene expression profiles of immortalized precursor oligodendrocytes (Oli-neu), stimulated to differentiate by the proximity of nontouching astrocytes or treatment with db-cAMP. In a previous paper, we reported that the astrocyte proximity upregulated or turned-on a large number of myelination genes and substantially enriched the Ca(2+)-signaling and cytokine receptor regulatory networks of MYE in Oli-neu cells. Here, we introduce the "transcriptomic distance" to evaluate fabric remodeling and "pair-wise relevance" to identify the most influential gene pairs. Together with the prominence gene analysis used to select and rank the fabric genes, these novel analytical tools provide a comprehensively quantitative view of the physio/pathological transformations of the transcriptomic programs of myelinating cells. Applied to our data, the analyses revealed not only that the astrocyte neighborhood is a substantially more powerful regulator of myelination than the differentiating treatment but also the molecular mechanisms of the two differentiating paradigms are different. By inducing a profound remodeling of MYE and regulatory transcriptomic networks, the astrocyte-oligodendrocyte intercommunication may be considered as a major player in both pathophysiology and therapy of neurodegenerative diseases related to myelination.