Probe-specific mixed-model approach to detect copy number differences using multiplex ligation-dependent probe amplification (MLPA)

Background: MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions. In this paper, we propose a method for the normalization procedure based on a non-linear mixed-model, as well as a new approach for determining the statistical significance of altered probes based on linear mixed-model. This method establishes a threshold by using different tolerance intervals that accommodates the specific random error variability observed in each test sample.Results: Through simulation studies we have shown that our proposed method outperforms two existing methods that are based on simple threshold rules or iterative regression. We have illustrated the method using a controlled MLPA assay in which targeted regions are variable in copy number in individuals suffering from different disorders such as Prader-Willi, DiGeorge or Autism showing the best performace.Conclusion: Using the proposed mixed-model, we are able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed.

Copy number variation modifies expression profiles over developmental time

Abstract : A preliminary understanding of the phenotypic effect of copy number variation (CNV) of DNA segments is emerging. These rearrangements were shown to influence, in a somewhat dose-dependent manner, the expression of genes mapping within them. They were also shown to modify the expression of genes located on their Hanks, sometimes at great distance. Here, we demonstrate by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, 'we find that some brain- expressed genes mapping within CNVS appear to be under compensatory loops only at specific time-points, indicating that the effect of CNVS on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time-course expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but its timing of expression as well. Résume : Nous commençons à comprendre les effets phénotypiques liés aux séquences d'ADN qui changent de nombre de copies d'un individu a l'autre. Des travaux précédents ont montré que ces variante de nombre de copies (CNVS) avaient une influence sur l'expression non seulement des gènes se trouvant dans le réarrangement, mais aussi sur ceux se trouvant à une certaine distance. Le présent travail étudie ces effets à différents stades du développement de la souris allant d'un embryon de deux semaines à la souris adulte. Nous avons observé que certains gènes exprimés dans le cerveau semblent soumis à un contrôle plus strict a certaines étapes du développement suggèrent que l'effet des CNVs est modulé différemment au cours de la vie. Notre travail sur trois souches différentes de souris a permis de montrer que les gènes ayant un profil d'expression différent dans le temps entre souches sont enrichis en gènes se trouvant dans des CNVs. Ceci nous amène à penser que les CNVs ont, non seulement une influence sur le niveau d'expression des gènes, mais aussi sur les moments durant lesquels ils seront exprimés. Résumé pour un large public : De nombreuses maladies sont dues soit a un gain (on parle alors de duplication) soit à une perte de matériel génétique (il s'agit dune délétion). Bien que les recherches visant à identifier les mécanismes moléculaires liés à ces réarrangements de notre génome progressent continuellement, la plupart des causes des maladies génétiques restent à élucider. Certaines parties de notre génome sont présentes en un nombre de copies qui diffère d'un individu à l'autre sans pour autant provoquer une ou des maladies. Ces segments d'ADN qui varient en nombre sont appelés Copy Number Variant (CNVs). Ils couvrent environ 12% de notre matériel génétique. Des études menées sur différents modèles animaux ont montré que les CNVs avaient une influence aussi bien sur les gènes qui sont a l'intérieur des CNVs que sur ceux qui sont dans leur voisinage. Ce travail étudie l'effet des CNVs à travers différents stades du développement de la souris. Nous avons démontré que les segments d'ADN qui varient en nombre de copies ont des effets variables selon le stade auxquels ils sont mesurés. Ainsi, les CNVs ont non seulement un impact sur l'expression des gènes présents dans ces régions et dans leur voisinage, mais influencent également leurs profils d'expression au cours du temps.

Genetic factors of obesity and eating disorders: Copy number variations (CNV) involving SH2B1

Context : It is now clearly shown that genetic factors in association with environment play a key role in obesity and eating disorders. This project studies the clinical symptoms and molecular abnormalities in patients carrying a strong hereditary predisposition to obesity and eating behavior disorders. We have previously published the association between the 16:29.5-30.1 deletion and a very penetrant form of morbid obesity and macrocephaly. We have also demonstrated the association between the reciprocal 16:29.5-30.1 duplication and underweight and small head circumference. These 2 studies demonstrate that gene dosage of one or several genes in this region regulates BMI as well as brain growth. At present, there are no data pointing towards particular candidate genes. We are currently investigating a second non-overlapping recurrent CNV encompassing SH2B1, upstream of the aforementioned rearrangement. SNPs in this gene have been associated with BMI in GWAS studies and mice models confirmed this association. Bokuchova et al have reported an association between deletions encompassing this gene and severe early onset obesity, as well as insulin resistance. We are currently collecting and analyzing data to fully characterize the phenotype and the transcriptional patterns associated with this rearrangement. Aims : 1. Identify carriers of any CNVs in the greater 16p11.2 region (between 16:28MB and 32MB) in the EGG consortium. 2. Perform association studies between SNPs in the greater 16p11.2 region (16:28-32MB) and anthropometric measures with adjusted "locus-wide significance", to identify or prioritize candidate genes potentially driving the association observed in patients with the CNVs (and thus worthy of further validation and sequencing). 3. Explore associations between GSV genome-wide and brain volume. 4. Explore relationship between brain volumes (whole brain and regional for those who underwent brain MRI), head circumference and BMI. 5. Extrapolate this procedure to other regions covered by the Metabochip. Methods : - Examine and collect clinical informations, as well as molecular informations in these patients. - Analysis of MRI data in children and adults with BMI > 2SD. Compare changes to MRI data obtained in patients with monogenic forms of obesity (data from Lausanne study) and to underweight (BMI<-2SD) individuals from EGG. - Test whether opposite extremes of the phenotypic distribution may be highly informative Expected results : This is a highly focused study, pertaining to approximately 1 0/00 of the human genome. Yet it is clear that if successful, the lessons learned from this study could be extrapolated to other segments of the genome and would need validation and replication by additional studies. Altogether they will contribute to further explore the missing heritability and point to etiologic genes and pathways underlying these important health burdens.

Copy number variation modifies expression time courses.

A preliminary understanding into the phenotypic effect of DNA segment copy number variation (CNV) is emerging. These rearrangements were demonstrated to influence, in a somewhat dose-dependent manner, the expression of genes that map within them. They were also shown to modify the expression of genes located on their flanks and sometimes those at a great distance from their boundary. Here we demonstrate, by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, we find that some brain-expressed genes mapping within CNVs appear to be under compensatory loops only at specific time points, indicating that the effect of CNVs on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time courses of expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but also the timing of its expression.

DNA Topoisomerase I Gene Copy Number and mRNA Expression Assessed as Predictive Biomarkers for Adjuvant Irinotecan in Stage II/III Colon Cancer.

PURPOSE: Prospective-retrospective assessment of theTOP1gene copy number andTOP1mRNA expression as predictive biomarkers for adjuvant irinotecan in stage II/III colon cancer. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded tissue microarrays were obtained from an adjuvant colon cancer trial (PETACC3) where patients were randomized to 5-fluorouracil/folinic acid with or without additional irinotecan.TOP1copy number status was analyzed by fluorescencein situhybridization (FISH) using aTOP1/CEN20 dual-probe combination.TOP1mRNA data were available from previous analyses. RESULTS: TOP1FISH and follow-up data were obtained from 534 patients.TOP1gain was identified in 27% using a single-probe enumeration strategy (≥4TOP1signals per cell) and in 31% when defined by aTOP1/CEN20 ratio ≥ 1.5. The effect of additional irinotecan was not dependent onTOP1FISH status.TOP1mRNA data were available from 580 patients with stage III disease. Benefit of irinotecan was restricted to patients characterized byTOP1mRNA expression ≥ third quartile (RFS: HRadjusted, 0.59;P= 0.09; OS: HRadjusted, 0.44;P= 0.03). The treatment byTOP1mRNA interaction was not statistically significant, but in exploratory multivariable fractional polynomial interaction analysis, increasingTOP1mRNA values appeared to be associated with increasing benefit of irinotecan. CONCLUSIONS: In contrast to theTOP1copy number, a trend was demonstrated for a predictive property ofTOP1mRNA expression. On the basis ofTOP1mRNA, it might be possible to identify a subgroup of patients where an irinotecan doublet is a clinically relevant option in the adjuvant setting of colon cancer.Clin Cancer Res; 22(7); 1621-31. ©2015 AACR.

Copy number variations of colorectal cancer by whole exome sequencing data

Colorectal cancer (CRC) is the third most common cancer and the fourth leading cause of cancer death worldwide. About 85% of the cases of CRC are known to have chromosomal instability, an allelic imbalance at several chromosomal loci, and chromosome amplification and translocation. The aim of this study is to determine the recurrent copy number variant (CNV) regions present in stage II of CRC through whole exome sequencing, a rapidly developing targeted next-generation sequencing (NGS) technology that provides an accurate alternative approach for accessing genomic variations. 42 normal-tumor paired samples were sequenced by Illumina Genome Analyzer. Data was analyzed with Varscan2 and segmentation was performed with R package R-GADA. Summary of the segments across all samples was performed and the result was overlapped with DEG data of the same samples from a previous study in the group1. Major and more recurrent segments of CNV were: gain of chromosome 7pq(13%), 13q(31%) and 20q(75%) and loss of 8p(25%), 17p(23%), and 18pq(27%). This results are coincident with the known literature of CNV in CRC or other cancers, but our methodology should be validated by array comparative genomic hybridisation (aCGH) profiling, which is currently the gold standard for genetic diagnosis of CNV.

The Tnt1 family member Retrosol copy number and structure disclose retrotransposon diversification in different Solanum species

Eukaryotic genome expansion/retraction caused by LTR-retrotransposon activity is dependent on the expression of full length copies to trigger efficient transposition and recombination-driven events. The Tnt1 family of retrotransposons has served as a model to evaluate the diversity among closely related elements within Solanaceae species and found that members of the family vary mainly in their U3 region of the long terminal repeats (LTRs). Recovery of a full length genomic copy of Retrosol was performed through a PCR-based approach from wild potato, Solanum oplocense. Further characterization focusing on both LTR sequences of the amplified copy allowed estimating an approximate insertion time at 2 million years ago thus supporting the occurrence of transposition cycles after genus divergence. Copy number of Tnt1-like elements in Solanum species were determined through genomic quantitative PCR whereby results sustain that Retrosol in Solanum species is a low copy number retrotransposon (1-4 copies) while Retrolyc1 has an intermediate copy number (38 copies) in S. peruvianum. Comparative analysis of retrotransposon content revealed no correlation between genome size or ploidy level and Retrosol copy number. The tetraploid cultivated potato with a cellular genome size of 1,715 Mbp harbours similar copy number per monoploid genome than other diploid Solanum species (613-884 Mbp). Conversely, S. peruvianum genome (1,125 Mbp) has a higher copy number. These results point towards a lineage specific dynamic flux regarding the history of amplification/activity of Tnt1-like elements in the genome of Solanum species.

Recurrent copy number gains of ACVR1 and corresponding transcript overexpression are associated with survival in head and neck squamous cell carcinomas

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Copy number of P elements, KP/full-sized P element ratio and their relationships with environmental factors in Brazilian Drosophila melanogaster populations

The P transposable element copy numbers and the KP/full-sized P element ratios were determined in eight Brazilian strains of Drosophila melanogaster. Strains from tropical regions showed lower overall P element copy numbers than did strains from temperate regions. Variable numbers of full-sized and defective elements were detected, but the full-sized P and KP elements were the predominant classes of elements in all strains. The full-sized P and KP element ratios were calculated and compared with latitude. The northernmost and southernmost Brazilian strains showed fewer full-sized elements than KP elements per genome, and the strains from less extreme latitudes had many more full-sized P than KP elements. However, no clinal variation was observed. Strains from different localities, previously classified as having P cytotype, displayed a higher or a lower proportion of KP elements than of full-sized P elements, as well as an equal number of the two element types, showing that the same phenotype may be produced by different underlying genomic components of the P-M system.

Characterization of hobo element copy number and integrity in Brasilian populations of Drosophila melanogaster

We have determined the copy number and the presence of full-size hobo transposable elements in eight Brasilian strains of Drosophila melanogaster. Genomic DNA was digested with AvaII and XhoI restriction enzymes, respectively, and probed with a 963 bp sequence of the hobo element. Variable numbers of full-sized and defective elements were detected in all strains. The range of the copy number was 22.13 +/- 4.52. Blots showed the presence of a 2.6 kb fragment, corresponding to the complete element, in all strains exception of one and the 1.0 kb sequence, correponding to the Th1 and Th2 repressor elements. There was neither association among copy numbers of hobo elements and latitude nor the mean annual temperatures in the original geographical region of each strain.

Copy number variation of individual cattle genomes using next-generation sequencing

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising similar to 55.6-Mbp sequence-476 of which (similar to 38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (similar to 52%, chi(2) test; P-value <0.05). In contrast, genes related to pathogen- and parasite-resistance, such as CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Análise do gene TUBB3 e de Copy Number Variations em pacientes com Polimicrogiria

A polimicrogiria (PMG) é uma malformação do córtex cerebral causada por falhas no seu desenvolvimento, caracterizando-se por um número excessivo de pequenos giros e laminação anormal, dando à superfície cortical uma aparência irregular e grosseira. A gravidade de suas manifestações clínicas se relaciona diretamente com a extensão da malformação e das regiões cerebrais afetadas, sendo que a presença de lesões bilaterais ou unilaterais extensas indica um pior prognóstico. Uma das síndromes de polimicrogiria mais freqüentes e, conseqüentemente, mais bem descritas clinicamente, é a polimicrogiria perisylviana bilateral (PPB). Essa forma de PMG atinge a região que tange a fenda Sylviana, podendo apresentar-se tanto unilateralmente quanto em ambos os hemisférios. Vários genes têm sido relacionados a diferentes formas de polimicrogiria, são eles AFF2,TUBA1A, TUBB2B e TUBA8, SRPX2 e WDR62. Estes genes já foram estudados pelo nosso grupo de pesquisa em um grupo de pacientes compostos de casos familiares e esporádicos, acometidos em sua maioria pela forma perisylviana de PMG. Nenhuma variante deletéria foi identificada nestes genes. Recentemente um novo gene foi implicado na etiologia molecular das PMG, o TUBB3. O gene em questão pertence à mesma família de TUBA1A, TUBB2B e TUBA8 e codifica uma proteína de ligação aos microtúbulos, tendo importante papel na formação do fuso. Além deste gene, também tem sido descritas alterações genômicas, denominadas de Copy Number Variations (CNV), estas variações estruturais tem sido associadas com diversos distúrbios neurológicos, que vão desde transtornos psiquiátricos até malformações do córtex cerebral como a PMG. Desta forma, o objetivo deste trabalho foi analisar a existência de alterações de ponto deletérias no gene TUBB3 em pacientes com PMG e também, o envolvimento de CNVsna etiologia deste tipo de malformação ...

Gene copy number and differential gene expression in haploid and diploid males of the stingless bee, Melipona quadrifasciata

Complementary sex determination in Hymenoptera implies that heterozygosity at the sex locus leads to the development of diploid females, whereas hemizygosity results in haploid males. Diploid males can arise through inbreeding. In social species, these pose a double burden on colony fitness, from significant reduction in its worker force and through being less viable and fertile than haploid males. Apart from being "misfits", diploid males are of interest to assess molecular correlates for possibly ploidy-related bionomic differences. Herein, we generated suppression subtractive cDNA libraries from newly emerged haploid and diploid males of the stingless bee Melipona quadrifasciata to enrich for differentially expressed genes. Gene Ontology classification revealed that in haploid males more DEGs were related to stress responsiveness, biosynthetic processes, reproductive processes and spermatogenesis, whereas in diploid ones differentially expressed genes were associated with cellular organization, nervous system development and amino acid transport were prevalent. Furthermore, both libraries contained over 40 % ESTs representing possibly novel transcripts. Quantitative RT-PCR analyses confirmed the differential expression of a representative DEG set in newly emerged males. Several muscle formation and energy metabolism-related genes were under-expressed in diploid males. On including 5-day-old males in the analysis, changes in transcript abundance during sexual maturation were revealed.

MSX2 copy number increase and craniosynostosis: copy number variation detected by array comparative genomic hybridization

Fleury Group

Single-Nucleotide Polymorphism and Copy Number Variation of the Multidrug Resistance-1 Locus of Plasmodium vivax: Local and Global Patterns

Emerging resistance to chloroquine (CQ) poses a major challenge for Plasmodium vivax malaria control, and nucleotide substitutions and copy number variation in the P. vivax multidrug resistance 1 (pvmdr-1) locus, which encodes a digestive vacuole membrane transporter, may modulate this phenotype. We describe patterns of genetic variation in pvmdr-1 alleles from Acre and Amazonas in northwestern Brazil, and compare then with those reported in other malaria-endemic regions. The pvmdr-1 mutation Y976F, which is associated with CQ resistance in Southeast Asia and Oceania, remains rare in northwestern Brazil (1.8%) and its prevalence mirrors that of CO resistance worldwide. Gene amplification of pvmdr-1, which is associated with mefloquine resistance but increased susceptibility to CO, remains relatively rare in northwestern Brazil (0.9%) and globally (< 4%), but became common (> 10%) in Tak Province, Thailand, possibly because of drug-mediated selection. The global database we have assembled provides a baseline for further studies of genetic variation in pvmdr-1 and drug resistance in P. vivax malaria.