Non-invasive prenatal testing leading to a maternal diagnosis of Charcot-Marie-Tooth neuropathy.

Non-invasive prenatal testing (NIPT) is increasingly used in routine practice due to its high sensitivity and specificity in detecting fetal chromosomal anomalies. Several reports have highlighted that NIPT can diagnose previously unsuspected malignancy or benign copy number variation in the expectant mother. We report a case in which NIPT detected a duplication involving the 17p11.2-17p12 region with possible Potocki-Lupski syndrome in the fetus. However, on further questioning, the mother revealed that she had Charcot-Marie-Tooth neuropathy type IA (CMT1A) and thus using array CGH, we were able to confirm that the 17p duplication was maternal in origin, included only the typical CMT1A region and that the fetus had a normal chromosome complement. Although it may be rare for a mother to have a pathogenic chromosome duplication/deletion, with the expansion in scope of NIPT from classic trisomies to global chromosome coverage and monogenic conditions, more examples of fortuitous maternal diagnosis will certainly be forthcoming and this should be taken into account during pre-test genetic counseling.

A large genomic deletion leads to enhancer adoption by the lamin B1 gene: a second path to autosomal dominant adult-onset demyelinating leukodystrophy (ADLD).

Chromosomal rearrangements with duplication of the lamin B1 (LMNB1) gene underlie autosomal dominant adult-onset demyelinating leukodystrophy (ADLD), a rare neurological disorder in which overexpression of LMNB1 causes progressive central nervous system demyelination. However, we previously reported an ADLD family (ADLD-1-TO) without evidence of duplication or other mutation in LMNB1 despite linkage to the LMNB1 locus and lamin B1 overexpression. By custom array-CGH, we further investigated this family and report here that patients carry a large (∼660 kb) heterozygous deletion that begins 66 kb upstream of the LMNB1 promoter. Lamin B1 overexpression was confirmed in further ADLD-1-TO tissues and in a postmortem brain sample, where lamin B1 was increased in the frontal lobe. Through parallel studies, we investigated both loss of genetic material and chromosomal rearrangement as possible causes of LMNB1 overexpression, and found that ADLD-1-TO plausibly results from an enhancer adoption mechanism. The deletion eliminates a genome topological domain boundary, allowing normally forbidden interactions between at least three forebrain-directed enhancers and the LMNB1 promoter, in line with the observed mainly cerebral localization of lamin B1 overexpression and myelin degeneration. This second route to LMNB1 overexpression and ADLD is a new example of the relevance of regulatory landscape modifications in determining Mendelian phenotypes.

Domains of genome-wide gene expression dysregulation in Down's syndrome.

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

Simultaneous identification and prioritization of variants in familial, de novo, and somatic genetic disorders with VariantMaster.

There is increasing interest in clinical genetics pertaining to the utilization of high-throughput sequencing data for accurate diagnoses of monogenic diseases. Moreover, massive whole-exome sequencing of tumors has provided significant advances in the understanding of cancer development through the recognition of somatic driver variants. To improve the identification of the variants from HTS, we developed VariantMaster, an original program that accurately and efficiently extracts causative variants in familial and sporadic genetic diseases. The algorithm takes into account predicted variants (SNPs and indels) in affected individuals or tumor samples and utilizes the row (BAM) data to robustly estimate the conditional probability of segregation in a family, as well as the probability of it being de novo or somatic. In familial cases, various modes of inheritance are considered: X-linked, autosomal dominant, and recessive (homozygosity or compound heterozygosity). Moreover, VariantMaster integrates phenotypes and genotypes, and employs Annovar to produce additional information such as allelic frequencies in the general population and damaging scores to further reduce the number of putative variants. As a proof of concept, we successfully applied VariantMaster to identify (1) de novo mutations in a previously described data set, (2) causative variants in a rare Mendelian genetic disease, and (3) known and new "driver" mutations in previously reported cancer data sets. Our results demonstrate that VariantMaster is considerably more accurate in terms of precision and sensitivity compared with previously published algorithms.

A single-nucleotide substitution mutator phenotype revealed by exome sequencing of human colon adenomas.

Oncogene-induced DNA replication stress is thought to drive genomic instability in cancer. In particular, replication stress can explain the high prevalence of focal genomic deletions mapping within very large genes in human tumors. However, the origin of single-nucleotide substitutions (SNS) in nonfamilial cancers is strongly debated. Some argue that cancers have a mutator phenotype, whereas others argue that the normal DNA replication error rates are sufficient to explain the number of observed SNSs. Here, we sequenced the exomes of 24, mostly precancerous, colon polyps. Analysis of the sequences revealed mutations in the APC, CTNNB1, and BRAF genes as the presumptive cancer-initiating events and many passenger SNSs. We used the number of SNSs in the various lesions to calculate mutation rates for normal colon and adenomas and found that colon adenomas exhibit a mutator phenotype. Interestingly, the SNSs in the adenomas mapped more often than expected within very large genes, where focal deletions in response to DNA replication stress also map. We propose that single-stranded DNA generated in response to oncogene-induced replication stress compromises the repair of deaminated cytosines and other damaged bases, leading to the observed SNS mutator phenotype.

Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma.

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.

Chromosome conformation capture uncovers potential genome-wide interactions between human conserved non-coding sequences.

Comparative analyses of various mammalian genomes have identified numerous conserved non-coding (CNC) DNA elements that display striking conservation among species, suggesting that they have maintained specific functions throughout evolution. CNC function remains poorly understood, although recent studies have identified a role in gene regulation. We hypothesized that the identification of genomic loci that interact physically with CNCs would provide information on their functions. We have used circular chromosome conformation capture (4C) to characterize interactions of 10 CNCs from human chromosome 21 in K562 cells. The data provide evidence that CNCs are capable of interacting with loci that are enriched for CNCs. The number of trans interactions varies among CNCs; some show interactions with many loci, while others interact with few. Some of the tested CNCs are capable of driving the expression of a reporter gene in the mouse embryo, and associate with the oligodendrocyte genes OLIG1 and OLIG2. Our results underscore the power of chromosome conformation capture for the identification of targets of functional DNA elements and raise the possibility that CNCs exert their functions by physical association with defined genomic regions enriched in CNCs. These CNC-CNC interactions may in part explain their stringent conservation as a group of regulatory sequences.

A systematic enhancer screen using lentivector transgenesis identifies conserved and non-conserved functional elements at the Olig1 and Olig2 locus.

Finding sequences that control expression of genes is central to understanding genome function. Previous studies have used evolutionary conservation as an indicator of regulatory potential. Here, we present a method for the unbiased in vivo screen of putative enhancers in large DNA regions, using the mouse as a model. We cloned a library of 142 overlapping fragments from a 200 kb-long murine BAC in a lentiviral vector expressing LacZ from a minimal promoter, and used the resulting vectors to infect fertilized murine oocytes. LacZ staining of E11 embryos obtained by first using the vectors in pools and then testing individual candidates led to the identification of 3 enhancers, only one of which shows significant evolutionary conservation. In situ hybridization and 3C/4C experiments suggest that this enhancer, which is active in the neural tube and posterior diencephalon, influences the expression of the Olig1 and/or Olig2 genes. This work provides a new approach for the large-scale in vivo screening of transcriptional regulatory sequences, and further demonstrates that evolutionary conservation alone seems too limiting a criterion for the identification of enhancers.

Detection of genomic variation by selection of a 9 mb DNA region and high throughput sequencing.

Detection of the rare polymorphisms and causative mutations of genetic diseases in a targeted genomic area has become a major goal in order to understand genomic and phenotypic variability. We have interrogated repeat-masked regions of 8.9 Mb on human chromosomes 21 (7.8 Mb) and 7 (1.1 Mb) from an individual from the International HapMap Project (NA12872). We have optimized a method of genomic selection for high throughput sequencing. Microarray-based selection and sequencing resulted in 260-fold enrichment, with 41% of reads mapping to the target region. 83% of SNPs in the targeted region had at least 4-fold sequence coverage and 54% at least 15-fold. When assaying HapMap SNPs in NA12872, our sequence genotypes are 91.3% concordant in regions with coverage > or = 4-fold, and 97.9% concordant in regions with coverage > or = 15-fold. About 81% of the SNPs recovered with both thresholds are listed in dbSNP. We observed that regions with low sequence coverage occur in close proximity to low-complexity DNA. Validation experiments using Sanger sequencing were performed for 46 SNPs with 15-20 fold coverage, with a confirmation rate of 96%, suggesting that DNA selection provides an accurate and cost-effective method for identifying rare genomic variants.

In vitro whole-genome analysis identifies a susceptibility locus for HIV-1.

Advances in large-scale analysis of human genomic variability provide unprecedented opportunities to study the genetic basis of susceptibility to infectious agents. We report here the use of an in vitro system for the identification of a locus on HSA8q24.3 associated with cellular susceptibility to HIV-1. This locus was mapped through quantitative linkage analysis using cell lines from multigeneration families, validated in vitro, and followed up by two independent association studies in HIV-positive individuals. Single nucleotide polymorphism rs2572886, which is associated with cellular susceptibility to HIV-1 in lymphoblastoid B cells and in primary T cells, was also associated with accelerated disease progression in one of two cohorts of HIV-1-infected patients. Biological analysis suggests a role of the rs2572886 region in the regulation of the LY6 family of glycosyl-phosphatidyl-inositol (GPI)-anchored proteins. Genetic analysis of in vitro cellular phenotypes provides an attractive approach for the discovery of susceptibility loci to infectious agents.

Genomewide analysis of nucleosome density histone acetylation and HDAC function in fission yeast.

We have conducted a genomewide investigation into the enzymatic specificity, expression profiles, and binding locations of four histone deacetylases (HDACs), representing the three different phylogenetic classes in fission yeast (Schizosaccharomyces pombe). By directly comparing nucleosome density, histone acetylation patterns and HDAC binding in both intergenic and coding regions with gene expression profiles, we found that Sir2 (class III) and Hos2 (class I) have a role in preventing histone loss; Clr6 (class I) is the principal enzyme in promoter-localized repression. Hos2 has an unexpected role in promoting high expression of growth-related genes by deacetylating H4K16Ac in their open reading frames. Clr3 (class II) acts cooperatively with Sir2 throughout the genome, including the silent regions: rDNA, centromeres, mat2/3 and telomeres. The most significant acetylation sites are H3K14Ac for Clr3 and H3K9Ac for Sir2 at their genomic targets. Clr3 also affects subtelomeric regions which contain clustered stress- and meiosis-induced genes. Thus, this combined genomic approach has uncovered different roles for fission yeast HDACs at the silent regions in repression and activation of gene expression.

Genomewide histone acetylation microarrays.

Histone acetylation and methylation are important regulators of gene activity. Chromatin immunoprecipitation (ChIP or ChrIP) has made it possible to examine not only the state of histone acetylation at a gene but also that of histone methylation and may soon be extended to other histone modifications such as phosphorylation and ubiquitination. In principle such studies are possible as long as an antibody is available to the particular histone modification. Once a target gene is identified it is instructive to see the effect of mutating putative enzymes responsible for the modification to determine how a particular enzyme is responsible for altering chromatin of that gene. Although specific target genes have been studied that contain such modifications recent technical advances have made it possible to study histone modifications genomewide. This not only allows for alternate views of particular paradigms to be investigated, but also uncovers chromosomal patterns of histone modification that would be missed in analyzing individual genes. We describe here an approach to rapidly study histone modifications genomewide by combining chromatin immunoprecipitation and DNA microarrays.

Genome-wide binding map of the histone deacetylase Rpd3 in yeast.

We describe the genome-wide distribution of the histone deacetylase and repressor Rpd3 and its associated proteins Ume1 and Ume6 in Saccharomyces cerevisiae. Using a new cross-linking protocol, we found that Rpd3 binds upstream of many individual genes and upstream of members of gene classes with similar functions in anabolic processes. In addition, Rpd3 is preferentially associated with promoters that direct high transcriptional activity. We also found that Rpd3 was absent from large sub-telomeric domains. We show by co-immunoprecipitation and by the high similarity of their binding maps that Ume1 interacts with Rpd3. In contrast, despite the known role of Ume6 in Rpd3 recruitment, only a limited number of the genes targeted by Rpd3 are also enriched for (or targeted by) Ume6. This suggests that Rpd3 is brought to many promoters by alternative recruiters, some of which may bind the putative cis-regulatory DNA elements that we have identified in sets of Rpd3 target genes. Finally, we show that comparing the genome-wide pattern of Rpd3 binding with gene expression and histone acetylation in the rpd3 Delta mutant strain reveals new sites of Rpd3 function.

Microarray deacetylation maps determine genome-wide functions for yeast histone deacetylases.

Yeast contains a family of five related histone deacetylases (HDACs) whose functions are known at few genes. Therefore, we used chromatin immunoprecipitation and intergenic microarrays to generate genome-wide HDAC enzyme activity maps. Rpd3 and Hda1 deacetylate mainly distinct promoters and gene classes where they are recruited largely by novel mechanisms. Hda1 also deacetylates subtelomeric domains containing normally repressed genes that are used instead for gluconeogenesis, growth on carbon sources other than glucose, and adverse growth conditions. These domains have certain features of heterochromatin but are distinct from subtelomeric heterochromatin repressed by the deacetylase Sir2. Finally, Hos1/Hos3 and Hos2 preferentially affect ribosomal DNA and ribosomal protein genes, respectively. Thus, acetylation microarrays uncover the "division of labor" for yeast histone deacetylases.