The swan genome and transcriptome, it is not all black and white.

BACKGROUND: The Australian black swan (Cygnus atratus) is an iconic species with contrasting plumage to that of the closely related northern hemisphere white swans. The relative geographic isolation of the black swan may have resulted in a limited immune repertoire and increased susceptibility to infectious diseases, notably infectious diseases from which Australia has been largely shielded. Unlike mallard ducks and the mute swan (Cygnus olor), the black swan is extremely sensitive to highly pathogenic avian influenza. Understanding this susceptibility has been impaired by the absence of any available swan genome and transcriptome information. RESULTS: Here, we generate the first chromosome-length black and mute swan genomes annotated with transcriptome data, all using long-read based pipelines generated for vertebrate species. We use these genomes and transcriptomes to show that unlike other wild waterfowl, black swans lack an expanded immune gene repertoire, lack a key viral pattern-recognition receptor in endothelial cells and mount a poorly controlled inflammatory response to highly pathogenic avian influenza. We also implicate genetic differences in SLC45A2 gene in the iconic plumage of the black swan. CONCLUSION: Together, these data suggest that the immune system of the black swan is such that should any avian viral infection become established in its native habitat, the black swan would be in a significant peril.

Comparison of long-read methods for sequencing and assembly of a plant genome.

BACKGROUND: Sequencing technologies have advanced to the point where it is possible to generate high-accuracy, haplotype-resolved, chromosome-scale assemblies. Several long-read sequencing technologies are available, and a growing number of algorithms have been developed to assemble the reads generated by those technologies. When starting a new genome project, it is therefore challenging to select the most cost-effective sequencing technology, as well as the most appropriate software for assembly and polishing. It is thus important to benchmark different approaches applied to the same sample. RESULTS: Here, we report a comparison of 3 long-read sequencing technologies applied to the de novo assembly of a plant genome, Macadamia jansenii. We have generated sequencing data using Pacific Biosciences (Sequel I), Oxford Nanopore Technologies (PromethION), and BGI (single-tube Long Fragment Read) technologies for the same sample. Several assemblers were benchmarked in the assembly of Pacific Biosciences and Nanopore reads. Results obtained from combining long-read technologies or short-read and long-read technologies are also presented. The assemblies were compared for contiguity, base accuracy, and completeness, as well as sequencing costs and DNA material requirements. CONCLUSIONS: The 3 long-read technologies produced highly contiguous and complete genome assemblies of M. jansenii. At the time of sequencing, the cost associated with each method was significantly different, but continuous improvements in technologies have resulted in greater accuracy, increased throughput, and reduced costs. We propose updating this comparison regularly with reports on significant iterations of the sequencing technologies.

A High-Fat Diet Increases Influenza A Virus-Associated Cardiovascular Damage.

BACKGROUND: Influenza A virus (IAV) causes a wide range of extrarespiratory complications. However, the role of host factors in these complications of influenza virus infection remains to be defined. METHODS: Here, we sought to use transcriptional profiling, virology, histology, and echocardiograms to investigate the role of a high-fat diet in IAV-associated cardiac damage. RESULTS: Transcriptional profiling showed that, compared to their low-fat counterparts (LF mice), mice fed a high-fat diet (HF mice) had impairments in inflammatory signaling in the lung and heart after IAV infection. This was associated with increased viral titers in the heart, increased left ventricular mass, and thickening of the left ventricular wall in IAV-infected HF mice compared to both IAV-infected LF mice and uninfected HF mice. Retrospective analysis of clinical data revealed that cardiac complications were more common in patients with excess weight, an association which was significant in 2 out of 4 studies. CONCLUSIONS: Together, these data provide the first evidence that a high-fat diet may be a risk factor for the development of IAV-associated cardiovascular damage and emphasizes the need for further clinical research in this area.

Single-Cell Transcriptomic Analysis of Cardiac Differentiation from Human PSCs Reveals HOPX-Dependent Cardiomyocyte Maturation.

Cardiac differentiation of human pluripotent stem cells (hPSCs) requires orchestration of dynamic gene regulatory networks during stepwise fate transitions but often generates immature cell types that do not fully recapitulate properties of their adult counterparts, suggesting incomplete activation of key transcriptional networks. We performed extensive single-cell transcriptomic analyses to map fate choices and gene expression programs during cardiac differentiation of hPSCs and identified strategies to improve in vitro cardiomyocyte differentiation. Utilizing genetic gain- and loss-of-function approaches, we found that hypertrophic signaling is not effectively activated during monolayer-based cardiac differentiation, thereby preventing expression of HOPX and its activation of downstream genes that govern late stages of cardiomyocyte maturation. This study therefore provides a key transcriptional roadmap of in vitro cardiac differentiation at single-cell resolution, revealing fundamental mechanisms underlying heart development and differentiation of hPSC-derived cardiomyocytes.

Single-cell RNA-seq of human induced pluripotent stem cells reveals cellular heterogeneity and cell state transitions between subpopulations.

Heterogeneity of cell states represented in pluripotent cultures has not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC-CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method and, through this, identified four subpopulations distinguishable on the basis of their pluripotent state, including a core pluripotent population (48.3%), proliferative (47.8%), early primed for differentiation (2.8%), and late primed for differentiation (1.1%). For each subpopulation, we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets composed of 165 unique genes that denote the specific pluripotency states; using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to threefold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations and support our conclusions with results from two orthogonal pseudotime trajectory methods.

Targeted Next-Generation Sequencing for Detecting MLL Gene Fusions in Leukemia.

Mixed lineage leukemia (MLL) gene rearrangements characterize approximately 70% of infant and 10% of adult and therapy-related leukemia. Conventional clinical diagnostics, including cytogenetics and fluorescence in situ hybridization (FISH) fail to detect MLL translocation partner genes (TPG) in many patients. Long-distance inverse (LDI)-PCR, the "gold standard" technique that is used to characterize MLL breakpoints, is laborious and requires a large input of genomic DNA (gDNA). To overcome the limitations of current techniques, a targeted next-generation sequencing (NGS) approach that requires low RNA input was tested. Anchored multiplex PCR-based enrichment (AMP-E) was used to rapidly identify a broad range of MLL fusions in patient specimens. Libraries generated using Archer FusionPlex Heme and Myeloid panels were sequenced using the Illumina platform. Diagnostic specimens (n = 39) from pediatric leukemia patients were tested with AMP-E and validated by LDI-PCR. In concordance with LDI-PCR, the AMP-E method successfully identified TPGs without prior knowledge. AMP-E identified 10 different MLL fusions in the 39 samples. Only two specimens were discordant; AMP-E successfully identified a MLL-MLLT1 fusion where LDI-PCR had failed to determine the breakpoint, whereas a MLL-MLLT3 fusion was not detected by AMP-E due to low expression of the fusion transcript. Sensitivity assays demonstrated that AMP-E can detect MLL-AFF1 in MV4-11 cell dilutions of 10-7 and transcripts down to 0.005 copies/ng.Implications: This study demonstrates a NGS methodology with improved sensitivity compared with current diagnostic methods for MLL-rearranged leukemia. Furthermore, this assay rapidly and reliably identifies MLL partner genes and patient-specific fusion sequences that could be used for monitoring minimal residual disease. Mol Cancer Res; 16(2); 279-85. ©2017 AACR.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours.

This corrects the article DOI: 10.1038/nature21063.

Whole-genome landscape of pancreatic neuroendocrine tumours.

The diagnosis of pancreatic neuroendocrine tumours (PanNETs) is increasing owing to more sensitive detection methods, and this increase is creating challenges for clinical management. We performed whole-genome sequencing of 102 primary PanNETs and defined the genomic events that characterize their pathogenesis. Here we describe the mutational signatures they harbour, including a deficiency in G:C > T:A base excision repair due to inactivation of MUTYH, which encodes a DNA glycosylase. Clinically sporadic PanNETs contain a larger-than-expected proportion of germline mutations, including previously unreported mutations in the DNA repair genes MUTYH, CHEK2 and BRCA2. Together with mutations in MEN1 and VHL, these mutations occur in 17% of patients. Somatic mutations, including point mutations and gene fusions, were commonly found in genes involved in four main pathways: chromatin remodelling, DNA damage repair, activation of mTOR signalling (including previously undescribed EWSR1 gene fusions), and telomere maintenance. In addition, our gene expression analyses identified a subgroup of tumours associated with hypoxia and HIF signalling.

Hypermutation In Pancreatic Cancer.

Pancreatic cancer is molecularly diverse, with few effective therapies. Increased mutation burden and defective DNA repair are associated with response to immune checkpoint inhibitors in several other cancer types. We interrogated 385 pancreatic cancer genomes to define hypermutation and its causes. Mutational signatures inferring defects in DNA repair were enriched in those with the highest mutation burdens. Mismatch repair deficiency was identified in 1% of tumors harboring different mechanisms of somatic inactivation of MLH1 and MSH2. Defining mutation load in individual pancreatic cancers and the optimal assay for patient selection may inform clinical trial design for immunotherapy in pancreatic cancer.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus).

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake's body length and nestle within the rib cavity. Despite the snake's notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

Assessing the genetic diversity of Cu resistance in mine tailings through high-throughput recovery of full-length copA genes.

Characterizing the genetic diversity of microbial copper (Cu) resistance at the community level remains challenging, mainly due to the polymorphism of the core functional gene copA. In this study, a local BLASTN method using a copA database built in this study was developed to recover full-length putative copA sequences from an assembled tailings metagenome; these sequences were then screened for potentially functioning CopA using conserved metal-binding motifs, inferred by evolutionary trace analysis of CopA sequences from known Cu resistant microorganisms. In total, 99 putative copA sequences were recovered from the tailings metagenome, out of which 70 were found with high potential to be functioning in Cu resistance. Phylogenetic analysis of selected copA sequences detected in the tailings metagenome showed that topology of the copA phylogeny is largely congruent with that of the 16S-based phylogeny of the tailings microbial community obtained in our previous study, indicating that the development of copA diversity in the tailings might be mainly through vertical descent with few lateral gene transfer events. The method established here can be used to explore copA (and potentially other metal resistance genes) diversity in any metagenome and has the potential to exhaust the full-length gene sequences for downstream analyses.

Integrated genomic and transcriptomic analysis of human brain metastases identifies alterations of potential clinical significance.

Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2) = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.

Recommendations for Accurate Resolution of Gene and Isoform Allele-Specific Expression in RNA-Seq Data.

Genetic variation modulates gene expression transcriptionally or post-transcriptionally, and can profoundly alter an individual's phenotype. Measuring allelic differential expression at heterozygous loci within an individual, a phenomenon called allele-specific expression (ASE), can assist in identifying such factors. Massively parallel DNA and RNA sequencing and advances in bioinformatic methodologies provide an outstanding opportunity to measure ASE genome-wide. In this study, matched DNA and RNA sequencing, genotyping arrays and computationally phased haplotypes were integrated to comprehensively and conservatively quantify ASE in a single human brain and liver tissue sample. We describe a methodological evaluation and assessment of common bioinformatic steps for ASE quantification, and recommend a robust approach to accurately measure SNP, gene and isoform ASE through the use of personalized haplotype genome alignment, strict alignment quality control and intragenic SNP aggregation. Our results indicate that accurate ASE quantification requires careful bioinformatic analyses and is adversely affected by sample specific alignment confounders and random sampling even at moderate sequence depths. We identified multiple known and several novel ASE genes in liver, including WDR72, DSP and UBD, as well as genes that contained ASE SNPs with imbalance direction discordant with haplotype phase, explainable by annotated transcript structure, suggesting isoform derived ASE. The methods evaluated in this study will be of use to researchers performing highly conservative quantification of ASE, and the genes and isoforms identified as ASE of interest to researchers studying those loci.

Genomic catastrophes frequently arise in esophageal adenocarcinoma and drive tumorigenesis.

Oesophageal adenocarcinoma (EAC) incidence is rapidly increasing in Western countries. A better understanding of EAC underpins efforts to improve early detection and treatment outcomes. While large EAC exome sequencing efforts to date have found recurrent loss-of-function mutations, oncogenic driving events have been underrepresented. Here we use a combination of whole-genome sequencing (WGS) and single-nucleotide polymorphism-array profiling to show that genomic catastrophes are frequent in EAC, with almost a third (32%, n=40/123) undergoing chromothriptic events. WGS of 22 EAC cases show that catastrophes may lead to oncogene amplification through chromothripsis-derived double-minute chromosome formation (MYC and MDM2) or breakage-fusion-bridge (KRAS, MDM2 and RFC3). Telomere shortening is more prominent in EACs bearing localized complex rearrangements. Mutational signature analysis also confirms that extreme genomic instability in EAC can be driven by somatic BRCA2 mutations. These findings suggest that genomic catastrophes have a significant role in the malignant transformation of EAC.

A workflow to increase verification rate of chromosomal structural rearrangements using high-throughput next-generation sequencing.

Somatic rearrangements, which are commonly found in human cancer genomes, contribute to the progression and maintenance of cancers. Conventionally, the verification of somatic rearrangements comprises many manual steps and Sanger sequencing. This is labor intensive when verifying a large number of rearrangements in a large cohort. To increase the verification throughput, we devised a high-throughput workflow that utilizes benchtop next-generation sequencing and in-house bioinformatics tools to link the laboratory processes. In the proposed workflow, primers are automatically designed. PCR and an optional gel electrophoresis step to confirm the somatic nature of the rearrangements are performed. PCR products of somatic events are pooled for Ion Torrent PGM and/or Illumina MiSeq sequencing, the resulting sequence reads are assembled into consensus contigs by a consensus assembler, and an automated BLAT is used to resolve the breakpoints to base level. We compared sequences and breakpoints of verified somatic rearrangements between the conventional and high-throughput workflow. The results showed that next-generation sequencing methods are comparable to conventional Sanger sequencing. The identified breakpoints obtained from next-generation sequencing methods were highly accurate and reproducible. Furthermore, the proposed workflow allows hundreds of events to be processed in a shorter time frame compared with the conventional workflow.

The first mouse mutants of D14Abb1e (Fam208a) show that it is critical for early development.

An ENU mutagenesis screen to identify novel epigenetic modifiers was established in mice carrying a multi-copy GFP transgene, which is expressed in a variegated manner in erythrocytes and is highly sensitive to epigenetic silencing. The screen has produced mouse mutants of both known modifiers of epigenetic state, such as Dnmt1 and Smarca5, and novel modifiers, such as Smchd1 and Rlf. Here we report two mouse lines generated from the screen, MommeD6 and MommeD20, with point mutations in D14Abb1e. These are the first mouse mutants of D14Abb1e (also known as Fam208a), a gene about which little is known. Heterozygous intercrosses show that homozygous mutants from both the MommeD6 and MommeD20 lines are not viable beyond gastrulation, demonstrating an important role for D14Abb1e in development. We demonstrate that haploinsufficiency for D14Abb1e effects transgene expression at the RNA level. Analysis of the predicted D14Abb1e protein sequence reveals that it contains putative nuclear localisation signals and a domain of unknown function, DUF3715. Our studies reveal that D14Abb1e is localised to the nucleus and is expressed in skin and testes.

Somatic point mutation calling in low cellularity tumors.

Somatic mutation calling from next-generation sequencing data remains a challenge due to the difficulties of distinguishing true somatic events from artifacts arising from PCR, sequencing errors or mis-mapping. Tumor cellularity or purity, sub-clonality and copy number changes also confound the identification of true somatic events against a background of germline variants. We have developed a heuristic strategy and software (http://www.qcmg.org/bioinformatics/qsnp/) for somatic mutation calling in samples with low tumor content and we show the superior sensitivity and precision of our approach using a previously sequenced cell line, a series of tumor/normal admixtures, and 3,253 putative somatic SNVs verified on an orthogonal platform.

Thermal stress responses in the bacterial biosphere of the Great Barrier Reef sponge, Rhopaloeides odorabile.

Marine sponges are diverse, abundant and provide a crucial coupling point between benthic and pelagic habitats due to their high filtration rates. They also harbour extensive microbial communities, with many microbial phylotypes found exclusively in sponge hosts and not in the seawater or surrounding environment, i.e. so-called sponge-specific clusters (SCs) or sponge- and coral-specific clusters (SCCs). We employed DNA (16S rRNA gene) and RNA (16S rRNA)-based amplicon pyrosequencing to investigate the effects of sublethal thermal stress on the bacterial biosphere of the Great Barrier Reef sponge Rhopaloeides odorabile. A total of 8381 operational taxonomic units (OTUs) (97% sequence similarity) were identified, affiliated with 32 bacterial phyla from seawater samples, 23 bacterial phyla from sponge DNA extracts and 18 bacterial phyla from sponge RNA extracts. Sublethal thermal stress (31°C) had no effect on the present and/or active portions of the R. odorabile bacterial community but a shift in the bacterial assemblage was observed in necrotic sponges. Over two-thirds of DNA and RNA sequences could be assigned to previously defined SCs/SCCs in healthy sponges whereas only 12% of reads from necrotic sponges could be assigned to SCs/SCCs. A rapid decline in host health over a 1°C temperature increment suggests that sponges such as R. odorabile may be highly vulnerable to the effects of global climate change.

Reduced levels of two modifiers of epigenetic gene silencing, Dnmt3a and Trim28, cause increased phenotypic noise.

BACKGROUND: Inbred individuals reared in controlled environments display considerable variance in many complex traits but the underlying cause of this intangible variation has been an enigma. Here we show that two modifiers of epigenetic gene silencing play a critical role in the process. RESULTS: Inbred mice heterozygous for a null mutation in DNA methyltransferase 3a (Dnmt3a) or tripartite motif protein 28 (Trim28) show greater coefficients of variance in body weight than their wild-type littermates. Trim28 mutants additionally develop metabolic syndrome and abnormal behavior with incomplete penetrance. Genome-wide gene expression analyses identified 284 significantly dysregulated genes in Trim28 heterozygote mutants compared to wild-type mice, with Mas1, which encodes a G-protein coupled receptor implicated in lipid metabolism, showing the greatest average change in expression (7.8-fold higher in mutants). This gene also showed highly variable expression between mutant individuals. CONCLUSIONS: These studies provide a molecular explanation of developmental noise in whole organisms and suggest that faithful epigenetic control of transcription is central to suppressing deleterious levels of phenotypic variation. These findings have broad implications for understanding the mechanisms underlying sporadic and complex disease in humans.

A genome-wide screen for modifiers of transgene variegation identifies genes with critical roles in development.

BACKGROUND: Some years ago we established an N-ethyl-N-nitrosourea screen for modifiers of transgene variegation in the mouse and a preliminary description of the first six mutant lines, named MommeD1-D6, has been published. We have reported the underlying genes in three cases: MommeD1 is a mutation in SMC hinge domain containing 1 (Smchd1), a novel modifier of epigenetic gene silencing; MommeD2 is a mutation in DNA methyltransferase 1 (Dnmt1); and MommeD4 is a mutation in Smarca 5 (Snf2h), a known chromatin remodeler. The identification of Dnmt1 and Smarca5 attest to the effectiveness of the screen design. RESULTS: We have now extended the screen and have identified four new modifiers, MommeD7-D10. Here we show that all ten MommeDs link to unique sites in the genome, that homozygosity for the mutations is associated with severe developmental abnormalities and that heterozygosity results in phenotypic abnormalities and reduced reproductive fitness in some cases. In addition, we have now identified the underlying genes for MommeD5 and MommeD10. MommeD5 is a mutation in Hdac1, which encodes histone deacetylase 1, and MommeD10 is a mutation in Baz1b (also known as Williams syndrome transcription factor), which encodes a transcription factor containing a PHD-type zinc finger and a bromodomain. We show that reduction in the level of Baz1b in the mouse results in craniofacial features reminiscent of Williams syndrome. CONCLUSIONS: These results demonstrate the importance of dosage-dependent epigenetic reprogramming in the development of the embryo and the power of the screen to provide mouse models to study this process.