Cost of exome analysis in patients with intellectual disability: a micro-costing study in a French setting.

BACKGROUND: With the development of next generation sequencing technologies in France, exome sequencing (ES) has recently emerged as an opportunity to improve the diagnosis rate of patients presenting an intellectual disability (ID). To help French policy makers determine an adequate tariff for ES, we aimed to assess the unit cost per ES diagnostic test for ID from the preparation of the pre-analytical step until the report writing step and to identify its main cost drivers. METHODS: A micro-costing bottom-up approach was conducted for the year 2018 in a French setting as part of the DISSEQ study, a cost-effectiveness study funded by the Ministry of Health and performed in collaboration with the GAD (Génétique des Anomalies du Développement), a genetic team from the Dijon University Hospital, and a public sequencing platform, the Centre National de Recherche en Génomique Humaine (CNRGH). The analysis was conducted from the point of view of these two ES stakeholders. All of the resources (labor, equipment, disposables and reagents, reusable material) required to analyze blood samples were identified, collected and valued. Several sensitivity analyses were performed. RESULTS: The unit nominal cost per ES diagnostic test for ID was estimated to be €2,019.39. Labor represented 50.7% of the total cost. The analytical step (from the preparation of libraries to the analysis of sequences) represented 88% of the total cost. Sensitivity analyses suggested that a simultaneous price decrease of 20% for the capture kit and 50% for the sequencing support kit led to an estimation of €1,769 per ES diagnostic test for ID. CONCLUSION: This is the first estimation of ES cost to be done in the French setting of ID diagnosis. The estimation is especially influenced by the price of equipment kits, but more generally by the organization of the centers involved in the different steps of the analysis and the time period in which the study was conducted. This information can now be used to define an adequate tariff and assess the efficiency of ES. TRIAL REGISTRATION: ClinicalTrials.gov identifier NCT03287206 on September 19, 2017.

Genomic characterization of grapevine viruses N and O: novel vitiviruses from South Africa.

A survey was performed on a Vitis cultivar collection in Stellenbosch, South Africa. Metaviromes were generated for each cultivar, using an RNAtag-seq workflow. Analysis of assembled contigs indicated the presence of two putatively novel members of the genus Vitivirus, provisionally named "grapevine virus N" (GVN) and "grapevine virus O" (GVO). Comparisons of amino acid sequences showed that GVN and GVO are most closely related to grapevine virus G and grapevine virus E, respectively. The incidence of these novel viruses within the sampling site was low, with GVO and GVN associated with only five and two cultivars, respectively, of the 229 sampled.

Genomic and phylogenetic characterization of wheat yellows virus, a novel tenuivirus infecting wheat in South Africa.

During the 2019 winter wheat season, symptoms of severe chlorosis and stunting were observed on wheat in the irrigation production areas of South Africa. RNAtag-seq data were generated for seven samples from KwaZulu-Natal province and one from Limpopo. Analysis of assembled contigs indicated the presence of a putatively novel member of the genus Tenuivirus, tentatively named "wheat yellows virus" (WhYV). The genome is made up of four segments, which are 8952, 3451, 2338, and 2045 nucleotides in length and code for a total of seven ORFs. Phylogenies of each segment (nucleotide) and the polymerase gene (amino acid), as well as amino acid sequence comparisons of each gene product, showed that WhYV is most closely related to rice stripe virus.

Quantification of 1.5 T T1 and T2* Relaxation Times of Fetal Tissues in Uncomplicated Pregnancies.

BACKGROUND: Despite its many advantages, experience with fetal magnetic resonance imaging (MRI) is limited, as is knowledge of how fetal tissue relaxation times change with gestational age (GA). Quantification of fetal tissue relaxation times as a function of GA provides insight into tissue changes during fetal development and facilitates comparison of images across time and subjects. This, therefore, can allow the determination of biophysical tissue parameters that may have clinical utility. PURPOSE: To demonstrate the feasibility of quantifying previously unknown T1 and T2* relaxation times of fetal tissues in uncomplicated pregnancies as a function of GA at 1.5 T. STUDY TYPE: Pilot. POPULATION: Nine women with singleton, uncomplicated pregnancies (28-38 weeks GA). FIELD STRENGTH/SEQUENCE: All participants underwent two iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL-IQ) acquisitions at different flip angles (6° and 20°) at 1.5 T. ASSESSMENT: Segmentations of the lungs, liver, spleen, kidneys, muscle, and adipose tissue (AT) were conducted using water-only images and proton density fat fraction maps. Driven equilibrium single pulse observation of T1 (DESPOT1 ) was used to quantify the mean water T1 of the lungs, intraabdominal organs, and muscle, and the mean water and lipid T1 of AT. IDEAL T2* maps were used to quantify the T2* values of the lungs, intraabdominal organs, and muscle. STATISTICAL TESTS: F-tests were performed to assess the T1 and T2* changes of each analyzed tissue as a function of GA. RESULTS: No tissue demonstrated a significant change in T1 as a function of GA (lungs [P = 0.89]; liver [P = 0.14]; spleen [P = 0.59]; kidneys [P = 0.97]; muscle [P = 0.22]; AT: water [P = 0.36] and lipid [P = 0.14]). Only the spleen and muscle T2* showed a significant decrease as a function of GA (lungs [P = 0.67); liver [P = 0.05]; spleen [P < 0.05]; kidneys [P = 0.70]; muscle [P < 0.05]). DATA CONCLUSION: These preliminary data suggest that the T1 of the investigated tissues is relatively stable over 28-38 weeks GA, while the T2* change in spleen and muscle decreases significantly in that period. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY STAGE: 2.

Novel viruses associated with plants of the family Amaryllidaceae in South Africa.

Nineteen samples from members of the plant genera Agapanthus, Clivia, Hippeastrum, and Scadoxus were collected from gardens in the Gauteng and Western Cape provinces of South Africa. The plants displayed highly variable symptoms of viral disease, including chlorosis, necrosis, streaking, and ringspot. RNAtag-seq was used to characterize the associated viral populations. Plants of the genus Agapanthus were found to be associated with three novel viruses from the families Caulimoviridae, Closteroviridae, and Betaflexiviridae; plants of the genus Clivia were associated with novel members of the families Potyviridae and Betaflexiviridae; and plants of the genus Scadoxus were associated with a novel member of the family Tospoviridae. Nerine latent virus was associated with plants of the genera Agapanthus, Clivia, and Hippeastrum, while hippeastrum mosaic virus was associated exclusively with a Hippeastrum cultivar.

Water-fat magnetic resonance imaging of adipose tissue compartments in the normal third trimester fetus.

BACKGROUND: Assessment of fetal adipose tissue gives information about the future metabolic health of an individual, with evidence that the development of this tissue has regional heterogeneity. OBJECTIVE: To assess differences in the proton density fat fraction (PDFF) between fetal adipose tissue compartments in the third trimester using water-fat magnetic resonance imaging (MRI). MATERIALS AND METHODS: Water-fat MRI was performed in a 1.5-T scanner. Fetal adipose tissue was segmented into cheeks, thorax, abdomen, upper arms, forearms, thighs and lower legs. PDFF and R2* values were measured in each compartment. RESULTS: Twenty-eight women with singleton pregnancies were imaged between 28 and 38 weeks of gestation. At 30 weeks' gestation (n=22), the PDFF was statistically different between the compartments (P<0.0001), with the highest PDFF in cheeks, followed by upper arms, thorax, thighs, forearms, lower legs and abdomen. There were no statistical differences in the rate of PDFF change with gestational age between the white adipose tissue compartments (P=0.97). Perirenal brown adipose tissue had a different PDFF and R2* compared to white adipose tissue, while the rate of R2* change did not significantly change with gestational age between white adipose tissue compartments (P=0.96). CONCLUSION: Fetal adipose tissue accumulates lipids at a similar rate in all white adipose tissue compartments. PDFF variances between the compartments suggest that accumulation begins at different gestational ages, starting with cheeks, followed by extremities, trunk and abdomen. Additionally, MRI was able to detect differences in the PDFF between fetal brown adipose tissue and white adipose tissue.

Skeletal cell YAP and TAZ combinatorially promote bone development.

The functions of the paralogous transcriptional coactivators Yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) in bone are controversial. Each has been observed to promote or inhibit osteogenesis in vitro, with reports of both equivalent and divergent functions. Their combinatorial roles in bone physiology are unknown. We report that combinatorial YAP/TAZ deletion from skeletal lineage cells, using Osterix-Cre, caused an osteogenesis imperfecta-like phenotype with severity dependent on allele dose and greater phenotypic expressivity with homozygous TAZ vs. YAP ablation. YAP/TAZ deletion decreased bone accrual and reduced intrinsic bone material properties through impaired collagen content and organization. These structural and material defects produced spontaneous fractures, particularly in mice with homozygous TAZ deletion and caused neonatal lethality in dual homozygous knockouts. At the cellular level in vivo, YAP/TAZ ablation reduced osteoblast activity and increased osteoclast activity, in an allele dose-dependent manner, impairing bone accrual and remodeling. Transcriptionally, YAP/TAZ deletion and small-molecule inhibition of YAP/TAZ interaction with the transcriptional coeffector TEAD reduced osteogenic and collagen-related gene expression, both in vivo and in vitro. These data demonstrate that YAP and TAZ combinatorially promote bone development through regulation of osteoblast activity, matrix quality, and osteoclastic remodeling.-Kegelman, C. D., Mason, D. E., Dawahare, J. H., Horan, D. J., Vigil, G. D., Howard, S. S., Robling, A. G., Bellido, T. M., Boerckel, J. D. Skeletal cell YAP and TAZ combinatorially promote bone development.

Dermatological manifestations in cardiofaciocutaneous syndrome: a prospective multicentric study of 45 mutation-positive patients.

BACKGROUND: Data on dermatological manifestations of cardiofaciocutaneous syndrome (CFCS) remain heterogeneous and almost without expert dermatological classification. OBJECTIVES: To describe the dermatological manifestations of CFCS; to compare them with the literature findings; to assess those discriminating CFCS from other RASopathies, including Noonan syndrome (NS) and Costello syndrome (CS); and to test for dermatological phenotype-genotype correlations. METHODS: We performed a 4-year, large, prospective, multicentric, collaborative dermatological and genetic study. RESULTS: Forty-five patients were enrolled. Hair abnormalities were ubiquitous, including scarcity or absence of eyebrows and wavy or curly hair in 73% and 69% of patients, respectively. Keratosis pilaris (KP), ulerythema ophryogenes (UO), palmoplantar hyperkeratosis (PPHK) and multiple melanocytic naevi (MMN; over 50 naevi) were noted in 82%, 44%, 27% and 29% of patients, respectively. Scarcity or absence of eyebrows, association of UO and PPHK, diffuse KP and MMN best differentiated CFCS from NS and CS. Oral acitretin may be highly beneficial for therapeutic management of PPHK, whereas treatment of UO by topical sirolimus 1% failed. No significant dermatological phenotype-genotype correlation was determined. CONCLUSIONS: A thorough knowledge of CFCS skin manifestations would help in making a positive diagnosis and differentiating CFCS from CS and NS.

Dermatological manifestations in Noonan syndrome: a prospective multicentric study of 129 patients positive for mutation.

BACKGROUND: Data on dermatological manifestations of Noonan syndrome (NS) remain heterogeneous and are based on limited dermatological expertise. OBJECTIVES: To describe the dermatological manifestations of NS, compare them with the literature findings, and test for dermatological phenotype-genotype correlations with or without the presence of PTPN11 mutations. METHODS: We performed a large 4-year, prospective, multicentric, collaborative dermatological and genetic study. RESULTS: Overall, 129 patients with NS were enrolled, including 65 patients with PTPN11-NS, 34 patients with PTPN11-NS with multiple lentigines (NSML), and 30 patients with NS who had a mutation other than PTPN11. Easy bruising was the most frequent dermatological finding in PTPN11-NS, present in 53·8% of patients. Multiple lentigines and café-au-lait macules (n ≥ 3) were present in 94% and 80% of cases of NSML linked to specific mutations of PTPN11, respectively. Atypical forms of NSML could be associated with NS with RAF1 or NRAS mutations. In univariate analysis, patients without a PTPN11 mutation showed (i) a significantly higher frequency of keratinization disorders (P = 0·001), including keratosis pilaris (P = 0·005), ulerythema ophryogenes (P = 0·0001) and palmar and/or plantar hyperkeratosis (P = 0·06, trend association), and (ii) a significantly higher frequency of scarce scalp hair (P = 0·035) and scarce or absent eyelashes (P = 0·06, trend association) than those with PTPN11 mutations. CONCLUSIONS: The cutaneous phenotype of NS with a PTPN11 mutation is generally mild and nonspecific, whereas the absence of a PTPN11 mutation is associated with a high frequency of keratinization disorders and hair abnormalities.

MED13L-related intellectual disability: involvement of missense variants and delineation of the phenotype.

Molecular anomalies in MED13L, leading to haploinsufficiency, have been reported in patients with moderate to severe intellectual disability (ID) and distinct facial features, with or without congenital heart defects. Phenotype of the patients was referred to "MED13L haploinsufficiency syndrome." Missense variants in MED13L were already previously described to cause the MED13L-related syndrome, but only in a limited number of patients. Here we report 36 patients with MED13L molecular anomaly, recruited through an international collaboration between centers of expertise for developmental anomalies. All patients presented with intellectual disability and severe language impairment. Hypotonia, ataxia, and recognizable facial gestalt were frequent findings, but not congenital heart defects. We identified seven de novo missense variations, in addition to protein-truncating variants and intragenic deletions. Missense variants clustered in two mutation hot-spots, i.e., exons 15-17 and 25-31. We found that patients carrying missense mutations had more frequently epilepsy and showed a more severe phenotype. This study ascertains missense variations in MED13L as a cause for MED13L-related intellectual disability and improves the clinical delineation of the condition.

Wiedemann-Steiner syndrome as a major cause of syndromic intellectual disability: A study of 33 French cases.

Wiedemann-Steiner syndrome (WSS) is a rare syndromic condition in which intellectual disability (ID) is associated with hypertrichosis cubiti, short stature, and characteristic facies. Following the identification of the causative gene (KMT2A) in 2012, only 31 cases of WSS have been described precisely in the literature. We report on 33 French individuals with a KMT2A mutation confirmed by targeted gene sequencing, high-throughput sequencing or exome sequencing. Patients' molecular and clinical features were recorded and compared with the literature data. On the molecular level, we found 29 novel mutations. We observed autosomal dominant transmission of WSS in 3 families and mosaicism in one family. Clinically, we observed a broad phenotypic spectrum with regard to ID (mild to severe), the facies (typical or not of WSS) and associated malformations (bone, cerebral, renal, cardiac and ophthalmological anomalies). Hypertrichosis cubiti that was supposed to be pathognomonic in the literature was found only in 61% of our cases. This is the largest series of WSS cases yet described to date. A majority of patients exhibited suggestive features, but others were less characteristic, only identified by molecular diagnosis. The prevalence of WSS was higher than expected in patients with ID, suggesting than KMT2A is a major gene in ID.

TBL1XR1 mutations in Pierpont syndrome are not restricted to the recurrent p.Tyr446Cys mutation.

Pierpont syndrome is a rare and sporadic syndrome, including developmental delay, facial characteristics, and abnormal extremities. Recently, a recurrent de novo TBL1XR1 variant (c.1337A > G; p.Tyr446Cys) has been identified in eight patients by whole-exome sequencing. A dominant-negative effect of this mutation is strongly suspected, since patients with TBL1XR1 deletion and other variants predicting loss of function do not share the same phenotype. We report two patients with typical Pierpont-like syndrome features. Exome sequencing allowed identifying a de novo heterozygous missense TBL1XR1 variant in both patients, different from those already reported: p.Cys325Tyr and p.Tyr446His. The localization of these mutations and clinical features of Pierpont-like syndrome suggest that their functional consequences are comparable with the recurrent mutation previously described, and provided additional data to understand molecular mechanisms of TBL1XR1 anomalies.

Three-dimensional computer-assisted dissection of pancreatic lymphatic anatomy on human fetuses: a step toward automatic image alignment.

PURPOSE: Pancreatic cancer is the fourth cause of death by cancer worldwide. Lymph node (LN) involvement is known to be the main prognostic factor. However, lymphatic anatomy is complex and only partially characterized. The aim of the study was to study the pancreatic lymphatic system using computer-assisted anatomic dissection (CAAD) technique and also to update CAAD technique by automatizing slice alignment. METHODS: We dissected three human fetuses aged from 18 to 34 WA. 5-µm serial sections of duodeno-pancreas and spleen blocks were stained (hematoxylin-eosin, hematoxylin of Mayer and Masson trichrome), scanned, aligned and modeled in three dimensions. RESULTS: We observed a rich, diffuse but not systematized lymphatic network in the peri-pancreatic region. There was an equal distribution of LNs between the cephalic and body-tail portions. The lymphatic vascularization appeared in continuity from the celiac trunk to the distal ends of its hepatic and splenic arterial branches parallel to the nerve ramifications of the celiac plexus. We also observed a continuity between the drainage of the pancreatic head and the para-aortic region posteriorly. CONCLUSION: In view of the wealth of peri-pancreatic LNs, the number of LNs to harvest could be increased to improve nodal staging and prognostic evaluation. Pancreatic anatomy as described does not seem to be compatible with the sentinel LN procedure in pancreatic surgery. Finally, we are now able to offer an alternative to manual alignment with a semi-automated alignment.

Molecular, clinical and neuropsychological study in 31 patients with Kabuki syndrome and KMT2D mutations.

Kabuki syndrome (KS-OMIM 147920) is a rare developmental disease characterized by the association of multiple congenital anomalies and intellectual disability. This study aimed to investigate intellectual performance in children with KS and link the performance to several clinical features and molecular data. We recruited 31 children with KMT2D mutations who were 6 to 16 years old. They all completed the Weschler Intelligence Scale for Children, fourth edition. We calculated all indexes: the Full Scale Intellectual Quotient (FSIQ), Verbal Comprehension Index (VCI), Perceptive Reasoning Index (PRI), Processing Speed Index (PSI), and Working Memory Index (WMI). In addition, molecular data and several clinical symptoms were studied. FSIQ and VCI scores were 10 points lower for patients with a truncating mutation than other types of mutations. In addition, scores for FSIQ, VCI and PRI were lower for children with visual impairment than normal vision. We also identified a discrepancy in indexes characterized by high WMI and VCI and low PRI and PSI. We emphasize the importance of early identification and intensive care of visual disorders in patients with KS and recommend individual assessment of intellectual profile.

Saturation-compensated measurements for fluorescence lifetime imaging microscopy.

Fluorophore saturation is the key factor limiting the speed and excitation range of fluorescence lifetime imaging microscopy (FLIM). For example, fluorophore saturation causes incorrect lifetime measurements when using conventional frequency-domain FLIM at high excitation powers. In this Letter, we present an analytical theoretical description of this error and present a method for compensating for this error in order to extract correct lifetime measurements in the limit of fluorophore saturation. We perform a series of simulations and experiments to validate our methods. The simulations and experiments show a 13.2× and a 2.6× increase in excitation range, respectively. The presented method is based on algorithms that can be easily applied to existing FLIM setups.

Disruptions in a cluster of computationally identified enhancers near FOXC1 and GMDS may influence brain development.

Regulatory elements are more evolutionarily conserved and provide a larger mutational target than coding regions of the human genome, suggesting that mutations in non-coding regions contribute significantly to development and disease. Using a computational approach to predict gene regulatory enhancers, we found that many known and predicted embryonic enhancers cluster in genomic loci harboring development-associated genes. One of the densest clusters of predicted enhancers in the human genome is near the genes GMDS and FOXC1. GMDS encodes a short-chain mannose dehydrogenase enzyme involved in the regulation of hindbrain neural migration, and FOXC1 encodes a developmental transcription factor required for brain, heart, and eye development. We experimentally validate four novel enhancers in this locus and demonstrate that these enhancers show consistent activity during embryonic development in domains that overlap with the expression of FOXC1 and GMDS. These four enhancers contain binding motifs for several transcription factors, including the ZIC family of transcription factors. Removal of the ZIC binding sites significantly alters enhancer activity in three of these enhancers, reducing expression in the eye, hindbrain, and limb, suggesting a mechanism whereby ZIC family members may transcriptionally regulate FOXC1 and/or GMDS expression. Our findings uncover novel enhancer regions that may control transcription in a topological domain important for embryonic development.

Super-sensitivity multiphoton frequency-domain fluorescence lifetime imaging microscopy.

We present a series of experiments that demonstrate a super-sensitive chemical imaging technique based on multiphoton frequency-domain fluorescence lifetime imaging microscopy (MPM-FD-FLIM) that shows a 2× improvement in imaging speed compared to the theoretical limit of conventional MPM-FD-FLIM. Additionally, this technique produces unprecedented sensitivity over a large range of fluorescence lifetimes. These results are achieved through simple modifications to data analysis in a conventional MPM-FD-FLIM microscope and are based on an analytical model describing the signal-to-noise ratio (SNR) of a MPM-FD-FLIM system [J. Opt. Soc. Am. A33, B1 (2016)]. Here we experimentally validate this model.

Investigation of signal-to-noise ratio in frequency-domain multiphoton fluorescence lifetime imaging microscopy.

Multiphoton microscopy (MPM) combined with fluorescence lifetime imaging microscopy (FLIM) has enabled three-dimensional quantitative molecular microscopy in vivo. The signal-to-noise ratio (SNR), and thus the imaging rate of MPM-FLIM, which is fundamentally limited by the shot noise and fluorescence saturation, has not been quantitatively studied yet. In this paper, we investigate the SNR performance of the frequency-domain (FD) MPM-FLIM with two figures of merit: the photon economy in the limit of shot noise, and the normalized SNR in the limit of saturation. The theoretical results and Monte Carlo simulations find that two-photon FD-FLIM requires 50% fewer photons to achieve the same SNR as conventional one-photon FLIM. We also analytically show that the MPM-FD-FLIM can exploit the DC and higher harmonic components generated by nonlinear optical mixing of the excitation light to improve SNR, reducing the required number of photons by an additional 50%. Finally, the effect of fluorophore saturation on the experimental SNR performance is discussed.

Search for ReCQL4 mutations in 39 patients genotyped for suspected Rothmund-Thomson/Baller-Gerold syndromes.

Three overlapping conditions, namely Rothmund-Thomson (RTS), Baller-Gerold (BGS) and RAPADILINO syndromes, have been attributed to RECQL4 mutations. Differential diagnoses depend on the clinical presentation, but the numbers of known genes remain low, leading to the widespread prescription of RECQL4 sequencing. The aim of our study was therefore to determine the best clinical indicators for the presence of RECQL4 mutations in a series of 39 patients referred for RECQL4 molecular analysis and belonging to the RTS (27 cases) and BGS (12 cases) spectrum. One or two deleterious RECQL4 mutations were found in 10/27 patients referred for RTS diagnosis. Clinical and molecular reevaluation led to a different diagnosis in 7/17 negative cases, including Clericuzio-type poikiloderma with neutropenia, hereditary sclerosing poikiloderma, and craniosynostosis/anal anomalies/porokeratosis. No RECQL4 mutations were found in the BGS group without poikiloderma, confirming that RECQL4 sequencing was not indicated in this phenotype. One chromosomal abnormality and one TWIST mutation was found in this cohort. This study highlights the search for differential diagnoses before the prescription of RECQL4 sequencing in this clinically heterogeneous group. The combination of clinically defined subgroups and next-generation sequencing will hopefully bring to light new molecular bases of syndromes with poikiloderma, as well as BGS without poikiloderma.

Integrating diverse datasets improves developmental enhancer prediction.

Gene-regulatory enhancers have been identified using various approaches, including evolutionary conservation, regulatory protein binding, chromatin modifications, and DNA sequence motifs. To integrate these different approaches, we developed EnhancerFinder, a two-step method for distinguishing developmental enhancers from the genomic background and then predicting their tissue specificity. EnhancerFinder uses a multiple kernel learning approach to integrate DNA sequence motifs, evolutionary patterns, and diverse functional genomics datasets from a variety of cell types. In contrast with prediction approaches that define enhancers based on histone marks or p300 sites from a single cell line, we trained EnhancerFinder on hundreds of experimentally verified human developmental enhancers from the VISTA Enhancer Browser. We comprehensively evaluated EnhancerFinder using cross validation and found that our integrative method improves the identification of enhancers over approaches that consider a single type of data, such as sequence motifs, evolutionary conservation, or the binding of enhancer-associated proteins. We find that VISTA enhancers active in embryonic heart are easier to identify than enhancers active in several other embryonic tissues, likely due to their uniquely high GC content. We applied EnhancerFinder to the entire human genome and predicted 84,301 developmental enhancers and their tissue specificity. These predictions provide specific functional annotations for large amounts of human non-coding DNA, and are significantly enriched near genes with annotated roles in their predicted tissues and lead SNPs from genome-wide association studies. We demonstrate the utility of EnhancerFinder predictions through in vivo validation of novel embryonic gene regulatory enhancers from three developmental transcription factor loci. Our genome-wide developmental enhancer predictions are freely available as a UCSC Genome Browser track, which we hope will enable researchers to further investigate questions in developmental biology.