Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

Differences in the climatic adaptation of silver birch (Betula pendula) and downy birch (B. pubescens) in Finland based on male flowering phenology.

Male flowering was studied at the canopy level in 10 silver birch (Betula pendula Roth) stands from 8 localities and in 14 downy birch (B. pubescens Ehrh.) stands from 10 localities in Finland from 1963 to 1973. Distributions of cumulative pollen catches were compared to the normal Gaussian distribution. The basis for the timing of flowering was the 50 per cent point of the anthesis-fitted normal distribution. To eliminate effects of background pollen, only the central, normally distributed part of the cumulative distribution was used. Development up to the median point of the distribution was measured and tested in calendar days, in degree days (> 5 °C) and in period units. The count of each parameter began on and included March 19. Male flowering in silver birch occurred from late April to late June depending on latitude, and flowering in downy birch took place from early May to early July. The heat sums needed for male flowering varied in downy birch stands latitudinally but there was practically no latitudinal variation in heat sums needed for silver birch flowering. The amount of male flowering in stands of both birch species were found to have a large annual variation but without any clear periodicity. The between years pollen catch variation in stands of either birch species did not show any significant latitudinal correlation in contrast to Norway spruce stands. The period unit heat sum gave the most accurate forecast of the timing of flowering for 60 per cent of the silver birch stands and for 78.6 per cent of the for downy birch stands. Calendar days, however, gave the best forecast for silver birch in 25 per cent of the cases, while degree days gave the best forecast for downy birch in 21.4 per cent of the cases. Silver birch seems to have a local inclination for a more fixed flowering date compared to downy birch, which could mean a considerable photoperiodic influence on flowering time of silver birch. Silver birch and downy birch had different geographical correlations. Frequent hybridization of birch species occurs more often in northern Finland in than in more southern latitudes. The different timing in flowering caused increasing scatter in flowering times in the north, especially in the case of downy birch. The chance of simultaneous flowering of silver birch and downy birch so increased northwards due to a more variable climate and also higher altitudinal variations. Compared with conifers, the reproduction cycles of both birch species were found to be well protected from damage by frost.

Detection of clinical mastitis with the help of a thermal camera

Increasing dairy farm size and increase in automation in livestock production require that new methods are used to monitor animal health. In this study, a thermal camera was tested for its capacity to detect clinical mastitis. Mastitis was experimentally induced in 6 cows with 10 mu g of Escherichia coli lipopolysaccharide (LPS). The LPS was infused into the left forequarter of each cow, and the right forequarters served as controls. Clinical examination for systemic and local signs and sampling for indicators of inflammation in milk were carried out before morning and evening milking throughout the 5-d experimental period and more frequently on the challenge day. Thermal images of experimental and control quarters were taken at each sampling time from lateral and medial angles. The first signs of clinical mastitis were noted in all cows 2 h postchallenge and included changes in general appearance of the cows and local clinical signs in the affected udder quarter. Rectal temperature, milk somatic cell count, and electrical conductivity were increased 4 h postchallenge and milk N-acetyl-beta-D-glucosaminidase activity 8 h postchallenge. The thermal camera was successful in detecting the 1 to 1.5 degrees C temperature change on udder skin associated with clinical mastitis in all cows because temperature of the udder skin of the experimental and control quarters increased in line with the rectal temperature. Yet, local signs on the udder were seen before the rise in udder skin and body temperature. The udder represents a sensitive site for detection of any febrile disease using a noninvasive method. A thermal camera mounted in a milking or feeding parlor could detect temperature changes associated with clinical mastitis or other diseases in a dairy herd.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Impaired spermatogenesis in Finnish boars and bulls

Maintenance of breeding efficiency and high semen quality is essential for reproductive success in farm animals. Early recognition of possible inheritable factors causing infertility requires constant attention. This thesis focuses on describing different manifestations of impaired spermatogenesis, their impact on fertility and partly also their incidence in populations. The reasons for spermatogenic failure are various. An interruption of germ cell differentiation, spermatogenic arrest, can lead to infertility. The incidence of azoospermia was investigated in the 1996 2005 survey of Finnish AI and farm breeding boars. We focused on the diagnosis, testicular morphometry and the possible reasons for the condition. The incidence of azoospermia was significantly higher in Yorkshire boars than in the Landrace breed. The most common diagnosis in Yorkshire boars was germ cell arrest at the primary spermatocyte level. The second most frequent diagnosis in Yorkshire boars was segmental aplasia of the Wolffian ducts with idiopathic epididymal obstruction. Other reasons for azoospermia were infrequent. In the second study we investigated the incidence of two relatively well-defined specific sperm defects in Finnish Yorkshire and Landrace boars during the same survey, the immotile short-tail sperm (ISTS) defect and the knobbed acrosome (KA) defect. In the Finnish Yorkshire boars the inherited ISTS defect, and the probably inherited KA defect, were important causes of infertility during 1996 2005. The ISTS defect was found in 7.6% and the KA defect in 0.8% of the Yorkshire boars. No Landrace boars were diagnosed with either of these two defects. In the third study we described a new sterilizing sperm defect in an oligoasthenoterazoospermic bull. Because of its morphological characteristics this defect was termed the multinuclear-multiflagellar sperm (MNMFS) defect. The number of Sertoli cells in the seminiferous tubuli was highly increased in the MNMFS bull compared with the number in normal bulls. In the following two studies we used a combined approach of fluorescence in situ hybridization (FISH), flow cytometry and morphometric studies to provide information on the cytogenetic background of macrocephalic bull spermatozoa. We described cellular features of diploid spermatozoa and compared the failures in the first and second meiotic divisions. In the last study we describe how the transplantation of testicular cells was used to determine whether spermatogonia derived from donor animals are able to colonize and produce motile spermatozoa in immune-competent unrelated boars suffering the ISTS defect. Transplantation resulted in complete focal spermatogenesis, indicated by the appearance of motile spermatozoa and confirmed by genotyping.

The many ways to delimit species: hairs, genes and surface chemistry

Species identification forms the basis for understanding the diversity of the living world, but it is also a prerequisite for understanding many evolutionary patterns and processes. The most promising approach for correctly delimiting and identifying species is to integrate many types of information in the same study. Our aim was to test how cuticular hydro- carbons, traditional morphometrics, genetic polymorphisms in nuclear markers (allozymes and DNA microsatellites) and DNA barcoding (partial mitochondrial COI gene) perform in delimiting species. As an example, we used two closely related Formica ants, F. fusca and F. lemani, sampled from a sympatric population in the northern part of their distribu- tion. Morphological characters vary and overlap in different parts of their distribution areas, but cuticular hydrocarbons include a strong taxonomic signal and our aim is to test the degree to which morphological and genetic data correspond to the chemical data. In the morphological analysis, species were best separated by the combined number of hairs on pro- notum and mesonotum, but individual workers overlapped in hair numbers, as previously noted by several authors. Nests of the two species were separated but not clustered according to species in a Principal Component Analysis made on nuclear genetic data. However, model-based Bayesian clustering resulted in perfect separation of the species and gave no indication of hybridization. Furthermore, F. lemani and F. fusca did not share any mitochondrial haplotypes, and the species were perfectly separated in a phylogenetic tree. We conclude that F. fusca and F. lemani are valid species that can be separated in our study area relatively well with all methods employed. However, the unusually small genetic differen- tiation in nuclear markers (FST = 0.12) shows that they are closely related, and occasional hybridization between F. fusca and F. lemani cannot be ruled out.

Suolistopatogeenisten Escherichia coli -bakteerien esiintyvyys ripulipotilailla ja raa’oissa lihatuotteissa Burkina Fasossa sekä enterohemorraagisten E. coli -kantojen karakterisointi

Suolistopatogeeniset Escherichia coli -bakteerit eli ripulikolit aiheuttavat ihmisellä suolistoinfektioita. Kuten normaalimikrobiston E. coli -bakteerit, ne esiintyvät ihmisen lisäksi muiden nisäkkäiden, etenkin märehtijöiden, ja lintujen suolistossa. Lisäksi ne voivat esiintyä maaperässä ja vesistöissä. Ihminen voi saada tartunnan eläinperäisten elintarvikkeiden välityksellä tai juomalla eläinten tai ihmisen ulosteilla saastunutta vettä. Ripulikolit voidaan jakaa ainakin viiteen ryhmään perustuen niiden erilaisiin virulenssiominaisuuksiin: enteropatogeeninen E. coli (EPEC), enterotoksigeeninen E. coli (ETEC), enterohemorraaginen E. coli (EHEC), enteroinvasiivinen E. coli (EIEC) ja enteroaggregatiivinen E. coli (EAEC). EPEC aiheuttaa etenkin kehitysmaissa pikkulapsille ripulia. ETEC aiheuttaa turistiripulia ja vastasyntyneiden ripulia kehitysmaissa. EHEC aiheuttaa ripulia tai veriripulia, joka voi varsinkin pienillä lapsilla johtaa hemolyyttis-ureemiseen oireyhtymään (HUS) ja munuaisten vaurioitumiseen. EIEC aiheuttaa Shigellan kaltaista ripulia, joka voi olla veristä. EAEC on yhdistetty lähinnä pitkittyneisiin ripuleihin. Tutkimuksessa selvitettiin suolistopatogeenisten E. coli -bakteerien esiintyvyyttä Burkina Fasossa, josta ei ole saatavilla aikaisempaa tietoa ripulikolien esiintymisestä ihmisissä ja elintarvikkeissa. Ulostenäytteitä otettiin ripulia sairastavilta alle viisivuotiailta lapsilta maaseudulta kahdesta kylästä, Boromosta ja Gourcysta, ja maan pääkaupungista Ouagadougousta (110 näytettä). Lihanäytteitä (kanaa, nautaa, lammasta ja naudan suolta, jota käytetään ihmisravinnoksi) otettiin Ouagadougoun toreilla myytävistä kypsentämättömistä lihoista (120 näytettä). Näytteistä saadut bakteerisekaviljelmät tutkittiin monialukkeisella PCR-menetelmällä, joka tunnistaa viiden ripulikoliryhmän virulenssigeenejä. Lisäksi lihanäytteistä eristettiin 20 EHEC-kantaa shigatoksiinin stx-geenin havaitsemiseen perustuvalla pesäkehybridisaatiolla ja PCR-seulonnalla, ja karakterisoitiin mahdollisten virulenssiominaisuuksien selvittämiseksi. Tutkimus osoitti, että ripulikolien aiheuttamat suolistoinfektiot ovat yleisiä ripulia sairastavilla pikkulapsilla Burkina Fasossa. Ulostenäytteistä 59 % oli positiivisia. Useimmiten lapsilla esiintyi EAEC- (32 %) ETEC- (31 %) ja EPEC-patoryhmiä (20 %). EIEC- (2 %) ja EHEC-patoryhmiä (1 %) esiintyi vähän. Myös useamman patoryhmän sekainfektiot olivat yleisiä (24 %). Eri paikkakuntien välillä oli tilastollisesti merkitseviä eroja ripulikolien esiintymisessä. Gourcyssa ripulikoleja esiintyi useammin kuin Ouagadougoussa ja Boromossa. Tutkimuksessa kävi ilmi, että Ouagadougoun toreilla myytävissä lihoissa on paljon ripulikoleja. Lihanäytteistä 43 % oli positiivisia. Yleisimmin lihoissa esiintyi EHEC (28 %), EPEC (20 %), ETEC (8 %) ja EAEC (5 %). EIEC-ryhmää ei havaittu lihoissa. Myös useamman patoryhmän sekakontaminaatioita löytyi (17 %) lihoista. Ripulikolien esiintyvyydessä eri lihojen välillä ei ollut tilastollisesti merkitseviä eroja, kun tarkasteltiin kaikkia patoryhmiä yhdessä. Eri patoryhmien esiintyvyyttä tarkasteltaessa EHEC-patoryhmää ei esiintynyt ollenkaan kanassa ja ero oli tilastollisesti merkitsevä muihin lihoihin verrattuna. Lihoista eristetyt 20 EHEC-kantaa kuuluivat 14 eri serotyyppiin, joista osa on aikaisemmin eristetty suolistoinfektioihin ja HUSoireyhtymään sairastuneilta ihmisiltä. Kaikki kannat olivat stx1-positiivisia ja puolella oli lisäksi stx2-geeni, jota pidetään shigatoksiinin virulentimpana muotona. Kahdelta EHEC-kannalta löytyi myös ETECpatoryhmän lämpöstabiilin enterotoksiini Ia:n geeni eli kannat olivat kahden patoryhmän välimuotoa ja osoitus geenien siirtymisestä eri patoryhmien välillä. Vaikka nuorimmat näytteen antaneet lapsipotilaat tuskin söivät lihaa, sen voidaan ajatella silti olevan edustava näyte lasten elinympäristöstä, sillä lasten ruoka valmistetaan usein samoissa oloissa, joissa raakaa lihaa käsitellään. Saastunut liha voi siten olla pikkulasten ripulikoli-infektioiden aiheuttaja.

Microbial identification by detection of ligation probes on DNA microarray

Microbes in natural and artificial environments as well as in the human body are a key part of the functional properties of these complex systems. The presence or absence of certain microbial taxa is a correlate of functional status like risk of disease or course of metabolic processes of a microbial community. As microbes are highly diverse and mostly notcultivable, molecular markers like gene sequences are a potential basis for detection and identification of key types. The goal of this thesis was to study molecular methods for identification of microbial DNA in order to develop a tool for analysis of environmental and clinical DNA samples. Particular emphasis was placed on specificity of detection which is a major challenge when analyzing complex microbial communities. The approach taken in this study was the application and optimization of enzymatic ligation of DNA probes coupled with microarray read-out for high-throughput microbial profiling. The results show that fungal phylotypes and human papillomavirus genotypes could be accurately identified from pools of PCR amplicons generated from purified sample DNA. Approximately 1 ng/μl of sample DNA was needed for representative PCR amplification as measured by comparisons between clone sequencing and microarray. A minimum of 0,25 amol/μl of PCR amplicons was detectable from amongst 5 ng/μl of background DNA, suggesting that the detection limit of the test comprising of ligation reaction followed by microarray read-out was approximately 0,04%. Detection from sample DNA directly was shown to be feasible with probes forming a circular molecule upon ligation followed by PCR amplification of the probe. In this approach, the minimum detectable relative amount of target genome was found to be 1% of all genomes in the sample as estimated from 454 deep sequencing results. Signal-to-noise of contact printed microarrays could be improved by using an internal microarray hybridization control oligonucleotide probe together with a computational algorithm. The algorithm was based on identification of a bias in the microarray data and correction of the bias as shown by simulated and real data. The results further suggest semiquantitative detection to be possible by ligation detection, allowing estimation of target abundance in a sample. However, in practise, comprehensive sequence information of full length rRNA genes is needed to support probe design with complex samples. This study shows that DNA microarray has the potential for an accurate microbial diagnostic platform to take advantage of increasing sequence data and to replace traditional, less efficient methods that still dominate routine testing in laboratories. The data suggests that ligation reaction based microarray assay can be optimized to a degree that allows good signal-tonoise and semiquantitative detection.

Molecular diagnosis of developmental disorders

ABSTRACT Idiopathic developmental disorders (DDs) affect ~1% of the population worldwide. This being a considerable amount, efforts are being made to elucidate the disease mechanisms. One or several genetic factors cause 30-40% of DDs, and only 10% are caused by environmental factors. The remaining 50% of DD patients go undiagnosed, mostly due to a lack of diagnostic techniques. The cause in most undiagnosed cases is though to be a genetic factor or a combination of genetic and environmental factors. Despite the surge of new technologies entering the market, their implementation into diagnostic laboratories is hampered by costs, lack of information about the expected diagnostic yield, and the wide range of selection. This study evaluates new microarray methods in diagnosing idiopathic DDs, providing information about their added diagnostic value. Study I analysed 150 patients by array comparative genomic hybridization (array CGH, 44K and 244K), with a subsequent 18% diagnostic yield. These results are supported by other studies, indicating an enourmous added diagnostic value of array CGH, compared with conventional cytogenetic analysis. Nevertheless, 80% of the patients remained undiagnosed in Study I. In an effort to diagnose more patients, in Study IV the resolution was increased from 8.9 Kb of the 244K CGH array to 0.7 Kb, by using a single-nucleotide polymorphism (SNP) array. However, no additional pathogenic changes were detected in the 35 patients assessed, and thus, for diagnostic purposes, an array platform with ca 9 Kb resolution appears adequate. The recent vast increase in reports of detected aberrations and associated phenotypes has enabled characterization of several new syndromes first based on a common aberration and thereafter by delineation of common clinical characteristics. In Study II, a familial deletion at 9q22.2q22.32 with variable penetrance was described. Despite several reports of aberrations in the adjacent area at 9q associated with Gorlin syndrome, the patients in this family had a unique phenotype and did not present with the syndrome. In Study III, a familial duplication of chromosome 6p22.2 was described. The duplication caused increased expression of an important enzyme of the γ-aminobutyric acid (GABA) degradation pathway, causing oxidative stress of the brain, and thus, very likely, the mild mental retardation of these patients. These two case studies attempted to pinpoint candidate genes and to resolve the pathogenic mechanism causing the clinical characteristics of the patients. Presenting rare genetic and clinical findings to the international science and medical community enables interpretation of similar findings in other patients. The added value of molecular karyotyping in patients with idiopathic DD is evident. As a first line of testing, arrays with a median resolution of at least 9 Kb should be considered and further characterization of detected aberrations undertaken when possible. Diagnostic whole-exome sequencing may be the best option for patients who remain undiagnosed after high-resolution array analysis.