Detection of pathogenic copy number variants in children with idiopathic intellectual disability using 500 K SNP array genomic hybridization.

BACKGROUND: Array genomic hybridization is being used clinically to detect pathogenic copy number variants in children with intellectual disability and other birth defects. However, there is no agreement regarding the kind of array, the distribution of probes across the genome, or the resolution that is most appropriate for clinical use. RESULTS: We performed 500 K Affymetrix GeneChip array genomic hybridization in 100 idiopathic intellectual disability trios, each comprised of a child with intellectual disability of unknown cause and both unaffected parents. We found pathogenic genomic imbalance in 16 of these 100 individuals with idiopathic intellectual disability. In comparison, we had found pathogenic genomic imbalance in 11 of 100 children with idiopathic intellectual disability in a previous cohort who had been studied by 100 K GeneChip array genomic hybridization. Among 54 intellectual disability trios selected from the previous cohort who were re-tested with 500 K GeneChip array genomic hybridization, we identified all 10 previously-detected pathogenic genomic alterations and at least one additional pathogenic copy number variant that had not been detected with 100 K GeneChip array genomic hybridization. Many benign copy number variants, including one that was de novo, were also detected with 500 K array genomic hybridization, but it was possible to distinguish the benign and pathogenic copy number variants with confidence in all but 3 (1.9%) of the 154 intellectual disability trios studied. CONCLUSION: Affymetrix GeneChip 500 K array genomic hybridization detected pathogenic genomic imbalance in 10 of 10 patients with idiopathic developmental disability in whom 100 K GeneChip array genomic hybridization had found genomic imbalance, 1 of 44 patients in whom 100 K GeneChip array genomic hybridization had found no abnormality, and 16 of 100 patients who had not previously been tested. Effective clinical interpretation of these studies requires considerable skill and experience.

Novel deletions of 14q11.2 associated with developmental delay, cognitive impairment and similar minor anomalies in three children.

METHODS AND RESULTS: We identified de novo submicroscopic chromosome 14q11.2 deletions in two children with idiopathic developmental delay and cognitive impairment. Vancouver patient 5566 has a approximately 200 kb deletion and Vancouver patient 8326 has a approximately 1.6 Mb deletion. The Database of Chromosomal Imbalance and Phenotype in Humans using Ensembl Resources (DECIPHER) revealed a third patient with idiopathic developmental delay and cognitive impairment, DECIPHER patient 126, who has a approximately 1.1 Mb deletion of 14q11.2. The deletion of patient 5566 overlaps that of patient 126 and both of these deletions lie entirely within that of patient 8326. All three children have similar dysmorphic features, including widely-spaced eyes, short nose with flat nasal bridge, long philtrum, prominent Cupid's bow of the upper lip, full lower lip and similar auricular anomalies. CONCLUSION: The minimal common deletion region on chromosome 14q11.2 is only approximately 35 kb (from 20.897 to 20.932, University of California at Santa Cruz (UCSC) Genome Browser; build hg18, March 2006) and includes only two genes, SUPT16H and CHD8, which are good candidate genes for the phenotypes. The non-recurrent breakpoints of these patients, the presence of normal copy number variants in the region and the local genomic structure support the notion that this region has reduced stability.

Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data.

BACKGROUND: Genomic deletions and duplications are important in the pathogenesis of diseases, such as cancer and mental retardation, and have recently been shown to occur frequently in unaffected individuals as polymorphisms. Affymetrix GeneChip whole genome sampling analysis (WGSA) combined with 100 K single nucleotide polymorphism (SNP) genotyping arrays is one of several microarray-based approaches that are now being used to detect such structural genomic changes. The popularity of this technology and its associated open source data format have resulted in the development of an increasing number of software packages for the analysis of copy number changes using these SNP arrays. RESULTS: We evaluated four publicly available software packages for high throughput copy number analysis using synthetic and empirical 100 K SNP array data sets, the latter obtained from 107 mental retardation (MR) patients and their unaffected parents and siblings. We evaluated the software with regards to overall suitability for high-throughput 100 K SNP array data analysis, as well as effectiveness of normalization, scaling with various reference sets and feature extraction, as well as true and false positive rates of genomic copy number variant (CNV) detection. CONCLUSION: We observed considerable variation among the numbers and types of candidate CNVs detected by different analysis approaches, and found that multiple programs were needed to find all real aberrations in our test set. The frequency of false positive deletions was substantial, but could be greatly reduced by using the SNP genotype information to confirm loss of heterozygosity.

Strain-specific telomere length revealed by single telomere length analysis in Caenorhabditis elegans.

Terminal restriction fragment analysis is the only method currently available for measuring telomere length in Caenorhabditis elegans. Its limitations include low sensitivity and interference by the presence of interstitial telomeric sequences in the C.elegans genome. Here we report the adaptation of single telomere length analysis (STELA) to measure the length of telomeric repeats on the left arm of chromosome V in C.elegans. This highly sensitive PCR-based method allows telomere length measurement from as few as a single worm. The application of STELA to eight wild-type C.elegans strains revealed considerable strain-specific differences in telomere length. Within individual strains, short outlying telomeres were observed that were clearly distinct from the bulk telomere length distributions, suggesting that processes other than end-replication losses and telomerase-mediated lengthening may generate telomere length heterogeneity in C.elegans. The utility of this method was further demonstrated by the characterization of telomere shortening in mrt-2 mutants. We conclude that STELA appears to be a valuable tool for studying telomere biology in C.elegans.

Elevated mutant frequencies and predominance of G:C to A:T transition mutations in Msh6(-/-) small intestinal epithelium.

The DNA mismatch repair (MMR) system is primarily responsible for purging newly synthesized DNA of errors incurred during semi-conservative replication. Lesion recognition is initially carried out by one of two heterodimeric protein complexes, MutS(alpha) or MutS(beta). While the former, comprised of MSH2 and MSH6, recognizes mispairs as well as short (1-2 nucleotide) insertions/deletions (IDLs), the latter, made up of MSH2 and MSH3, is primarily responsible for recognizing 2-6 nucleotide IDLs. As most of the functional information on these heterodimers is derived from in vitro studies, it was of interest to study the in vivo consequences of a lack of MutS(alpha). To this end, Big Blue( trade mark ) mice, that carry a lacI(+) transgenic lambda shuttle-phage mutational reporter, were crossed with Msh6(-/-) mice to evaluate the specific contribution of MutS(alpha) to genome integrity. Consistent with the importance of MutS(alpha) in lesion surveillance, small intestine epithelial cell DNA derived from lacI(+) Msh6(-/-) mice exhibited striking increases (average of 41-fold) in spontaneous mutant frequencies. Furthermore, the lacI gene mutation spectrum was dominated by G:C to A:T transitions, highlighting the critical importance of the MutS(alpha) complex in suppressing this frequently observed type of spontaneous mutation.

LongSAGE profiling of nine human embryonic stem cell lines.

To facilitate discovery of novel human embryonic stem cell (ESC) transcripts, we generated 2.5 million LongSAGE tags from 9 human ESC lines. Analysis of this data revealed that ESCs express proportionately more RNA binding proteins compared with terminally differentiated cells, and identified novel ESC transcripts, at least one of which may represent a marker of the pluripotent state.

Oligonucleotide microarray analysis of genomic imbalance in children with mental retardation.

The cause of mental retardation in one-third to one-half of all affected individuals is unknown. Microscopically detectable chromosomal abnormalities are the most frequently recognized cause, but gain or loss of chromosomal segments that are too small to be seen by conventional cytogenetic analysis has been found to be another important cause. Array-based methods offer a practical means of performing a high-resolution survey of the entire genome for submicroscopic copy-number variants. We studied 100 children with idiopathic mental retardation and normal results of standard chromosomal analysis, by use of whole-genome sampling analysis with Affymetrix GeneChip Human Mapping 100K arrays. We found de novo deletions as small as 178 kb in eight cases, de novo duplications as small as 1.1 Mb in two cases, and unsuspected mosaic trisomy 9 in another case. This technology can detect at least twice as many potentially pathogenic de novo copy-number variants as conventional cytogenetic analysis can in people with mental retardation.

Effect of TERT and ATM on gene expression profiles in human fibroblasts.

Telomeres protect chromosomes from degradation, end-to-end fusion, and illegitimate recombination. Loss of telomeres may lead to cell death or senescence or may cause genomic instability, leading to tumor formation. Expression of human telomerase reverse transcriptase (TERT) in human fibroblast cells elongates their telomeres and extends their lifespan. Ataxia telangiectasia mutated (ATM) deficiency in A-T human fibroblasts results in accelerated telomere shortening, abnormal cell-cycle response to DNA damage, and early senescence. Gene expression profiling was performed by serial analysis of gene expression (SAGE) on BJ normal human skin fibroblasts, A-T cells, and BJ and A-T cells transduced with TERT cDNA and expressing telomerase activity. In the four SAGE libraries, 36,921 unique SAGE tags were detected. Pairwise comparisons between the libraries showed differential expression levels of 1%-8% of the tags. Transcripts affected by both TERT and ATM were identified according to expression patterns, making them good candidates for further studies of pathways affected by both TERT and ATM. These include MT2A, P4HB, LGALS1, CFL1, LDHA, S100A10, EIF3S8, RANBP9, and SEC63. These genes are involved in apoptosis or processes related to cell growth, and most have been found to be deregulated in cancer. Our results have provided further insight into the roles of TERT and ATM by identifying genes likely to be involved in their function. Supplementary material for this article can be found on the Genes, Chromosomes and Cancer website at http://www.interscience.wiley.com/jpages/1045-2257/suppmat/index.html.

Systematic recovery and analysis of full-ORF human cDNA clones.

The Mammalian Gene Collection (MGC) consortium (http://mgc.nci.nih.gov) seeks to establish publicly available collections of full-ORF cDNAs for several organisms of significance to biomedical research, including human. To date over 15,200 human cDNA clones containing full-length open reading frames (ORFs) have been identified via systematic expressed sequence tag (EST) analysis of a diverse set of cDNA libraries; however, further systematic EST analysis is no longer an efficient method for identifying new cDNAs. As part of our involvement in the MGC program, we have developed a scalable method for targeted recovery of cDNA clones to facilitate recovery of genes absent from the MGC collection. First, cDNA is synthesized from various RNAs, followed by polymerase chain reaction (PCR) amplification of transcripts in 96-well plates using gene-specific primer pairs flanking the ORFs. Amplicons are cloned into a sequencing vector, and full-length sequences are obtained. Sequences are processed and assembled using Phred and Phrap, and analyzed using Consed and a number of bioinformatics methods we have developed. Sequences are compared with the Reference Sequence (RefSeq) database, and validation of sequence discrepancies is attempted using other sequence databases including dbEST and dbSNP. Clones with identical sequence to RefSeq or containing only validated changes will become part of the MGC human gene collection. Clones containing novel splice variants or polymorphisms have also been identified. Our approach to clone recovery, applied at large scale, has the potential to recover many and possibly most of the genes absent from the MGC collection.