Long-range gene control and genetic disease

A twin analysis to estimate genetic and environmental factors contributing to variation in weighted gene co-expression network module eigengenes

Multivariate network-based analytic methods such as weighted gene co-expression network analysis are frequently applied to human and animal gene-expression data to estimate the first principal component of a module, or module eigengene (ME). MEs are interpreted as multivariate summaries of correlated gene-expression patterns and network connectivity across genes within a module. As such, they have the potential to elucidate the mechanisms by which molecular genomic variation contributes to individual differences in complex traits. Although increasingly used to test for associations between modules and complex traits, the genetic and environmental etiology of MEs has not been empirically established. It is unclear if, and to what degree, individual differences in blood-derived MEs reflect random variation versus familial aggregation arising from heritable or shared environmental influences. We used biometrical genetic analyses to estimate the contribution of genetic and environmental influences on MEs derived from blood lymphocytes collected on a sample of N = 661 older male twins from the Vietnam Era Twin Study of Aging (VETSA) whose mean age at assessment was 67.7 years (SD = 2.6 years, range = 62-74 years). Of the 26 detected MEs, 14 (56%) had statistically significant additive genetic variation with an average heritability of 44% (SD = 0.08, range = 35%-64%). Despite the relatively small sample size, this demonstration of significant family aggregation including estimates of heritability in 14 of the 26 MEs suggests that blood-based MEs are reliable and merit further exploration in terms of their associations with complex traits and diseases.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes^1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

iEnhancer-DLRA: identification of enhancers and their strengths by a self-attention fusion strategy for local and global features

Identification and classification of enhancers are highly significant because they play crucial roles in controlling gene transcription. Recently, several deep learning-based methods for identifying enhancers and their strengths have been developed. However, existing methods are usually limited because they use only local or only global features. The combination of local and global features is critical to further improve the prediction performance. In this work, we propose a novel deep learning-based method, called iEnhancer-DLRA, to identify enhancers and their strengths. iEnhancer-DLRA extracts local and multi-scale global features of sequences by using a residual convolutional network and two bidirectional long short-term memory networks. Then, a self-attention fusion strategy is proposed to deeply integrate these local and global features. The experimental results on the independent test dataset indicate that iEnhancer-DLRA performs better than nine existing state-of-the-art methods in both identification and classification of enhancers in almost all metrics. iEnhancer-DLRA achieves 13.8% (for identifying enhancers) and 12.6% (for classifying strengths) improvement in accuracy compared with the best existing state-of-the-art method. This is the first time that the accuracy of an enhancer identifier exceeds 0.9 and the accuracy of the enhancer classifier exceeds 0.8 on the independent test set. Moreover, iEnhancer-DLRA achieves superior predictive performance on the rice dataset compared with the state-of-the-art method RiceENN.

GRiNCH: simultaneous smoothing and detection of topological units of genome organization from sparse chromatin contact count matrices with matrix factorization

High-throughput chromosome conformation capture assays, such as Hi-C, have shown that the genome is organized into organizational units such as topologically associating domains (TADs), which can impact gene regulatory processes. The sparsity of Hi-C matrices poses a challenge for reliable detection of these units. We present GRiNCH, a constrained matrix-factorization-based approach for simultaneous smoothing and discovery of TADs from sparse contact count matrices. GRiNCH shows superior performance against seven TAD-calling methods and three smoothing methods. GRiNCH is applicable to multiple platforms including SPRITE and HiChIP and can predict novel boundary factors with potential roles in genome organization.

iEnhancer-EBLSTM: Identifying Enhancers and Strengths by Ensembles of Bidirectional Long Short-Term Memory

Enhancers are regulatory DNA sequences that could be bound by specific proteins named transcription factors (TFs). The interactions between enhancers and TFs regulate specific genes by increasing the target gene expression. Therefore, enhancer identification and classification have been a critical issue in the enhancer field. Unfortunately, so far there has been a lack of suitable methods to identify enhancers. Previous research has mainly focused on the features of the enhancer's function and interactions, which ignores the sequence information. As we know, the recurrent neural network (RNN) and long short-term memory (LSTM) models are currently the most common methods for processing time series data. LSTM is more suitable than RNN to address the DNA sequence. In this paper, we take the advantages of LSTM to build a method named iEnhancer-EBLSTM to identify enhancers. iEnhancer-ensembles of bidirectional LSTM (EBLSTM) consists of two steps. In the first step, we extract subsequences by sliding a 3-mer window along the DNA sequence as features. Second, EBLSTM model is used to identify enhancers from the candidate input sequences. We use the dataset from the study of Quang H et al. as the benchmarks. The experimental results from the datasets demonstrate the efficiency of our proposed model.

A holistic view of mouse enhancer architectures reveals analogous pleiotropic effects and correlation with human disease

BACKGROUND

Efforts to elucidate the function of enhancers in vivo are underway but their vast numbers alongside differing enhancer architectures make it difficult to determine their impact on gene activity. By systematically annotating multiple mouse tissues with super- and typical-enhancers, we have explored their relationship with gene function and phenotype.

RESULTS

Though super-enhancers drive high total- and tissue-specific expression of their associated genes, we find that typical-enhancers also contribute heavily to the tissue-specific expression landscape on account of their large numbers in the genome. Unexpectedly, we demonstrate that both enhancer types are preferentially associated with relevant 'tissue-type' phenotypes and exhibit no difference in phenotype effect size or pleiotropy. Modelling regulatory data alongside molecular data, we built a predictive model to infer gene-phenotype associations and use this model to predict potentially novel disease-associated genes.

CONCLUSION

Overall our findings reveal that differing enhancer architectures have a similar impact on mammalian phenotypes whilst harbouring differing cellular and expression effects. Together, our results systematically characterise enhancers with predicted phenotypic traits endorsing the role for both types of enhancers in human disease and disorders.

Rare and de novo duplications containing SHOX in clubfoot

Introduction

Congenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained.

Methods

Exome sequence data from 816 unrelated clubfoot cases and 2645 in-house controls were analysed using coverage data to identify rare CNVs. The precise size and location of duplications were then determined using high-density Affymetrix Cytoscan chromosomal microarray (CMA). Segregation in families and de novo status were determined using qantitative PCR.

Results

Chromosome Xp22.33 duplications involving SHOX were identified in 1.1% of cases (9/816) compared with 0.07% of in-house controls (2/2645) (p=7.98×10⁻⁵, OR=14.57) and 0.27% (38/13592) of Atherosclerosis Risk in Communities/the Wellcome Trust Case Control Consortium 2 controls (p=0.001, OR=3.97). CMA validation confirmed an overlapping 180.28 kb duplicated region that included SHOX exons as well as downstream non-coding regions. In four of six sporadic cases where DNA was available for unaffected parents, the duplication was de novo. The probability of four de novo mutations in SHOX by chance in a cohort of 450 sporadic clubfoot cases is 5.4×10^–10.

Conclusions

Microduplications of the pseudoautosomal chromosome Xp22.33 region (PAR1) containing SHOX and downstream enhancer elements occur in ~1% of patients with clubfoot. SHOX and regulatory regions have previously been implicated in skeletal dysplasia as well as idiopathic short stature, but have not yet been reported in clubfoot. SHOX duplications likely contribute to clubfoot pathogenesis by altering early limb development.

In silico prediction of high-resolution Hi-C interaction matrices

The three-dimensional (3D) organization of the genome plays an important role in gene regulation bringing distal sequence elements in 3D proximity to genes hundreds of kilobases away. Hi-C is a powerful genome-wide technique to study 3D genome organization. Owing to experimental costs, high resolution Hi-C datasets are limited to a few cell lines. Computational prediction of Hi-C counts can offer a scalable and inexpensive approach to examine 3D genome organization across multiple cellular contexts. Here we present HiC-Reg, an approach to predict contact counts from one-dimensional regulatory signals. HiC-Reg predictions identify topologically associating domains and significant interactions that are enriched for CCCTC-binding factor (CTCF) bidirectional motifs and interactions identified from complementary sources. CTCF and chromatin marks, especially repressive and elongation marks, are most important for HiC-Reg’s predictive performance. Taken together, HiC-Reg provides a powerful framework to generate high-resolution profiles of contact counts that can be used to study individual locus level interactions and higher-order organizational units of the genome.

Existing computational approaches to predict long-range regulatory interactions do not fully exploit high-resolution Hi-C datasets. Here the authors present a Random Forests regression-based approach to predict high-resolution Hi-C counts using one-dimensional regulatory genomic signals.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Versatility of multivalent orientation, inverted meiosis, and rescued fitness in holocentric chromosomal hybrids

Chromosomal rearrangements (e.g., fusions/fissions) have the potential to drive speciation. However, their accumulation in a population is generally viewed as unlikely, because chromosomal heterozygosity should lead to meiotic problems and aneuploid gametes. Canonical meiosis involves segregation of homologous chromosomes in meiosis I and sister chromatid segregation during meiosis II. In organisms with holocentric chromosomes, which are characterized by kinetic activity distributed along almost the entire chromosome length, this order may be inverted depending on their metaphase I orientation. Here we analyzed the evolutionary role of this intrinsic versatility of holocentric chromosomes, which is not available to monocentric ones, by studying F₁ to F₄ hybrids between two chromosomal races of the Wood White butterfly (Leptidea sinapis), separated by at least 24 chromosomal fusions/fissions. We found that these chromosomal rearrangements resulted in multiple meiotic multivalents, and, contrary to the theoretical prediction, the hybrids displayed relatively high reproductive fitness (42% of that of the control lines) and regular behavior of meiotic chromosomes. In the hybrids, we also discovered inverted meiosis, in which the first and critical stage of chromosome number reduction was replaced by the less risky stage of sister chromatid separation. We hypothesize that the ability to invert the order of the main meiotic events facilitates proper chromosome segregation and hence rescues fertility and viability in chromosomal hybrids, potentially promoting dynamic karyotype evolution and chromosomal speciation.

Structural variations causing inherited peripheral neuropathies: A paradigm for understanding genomic organization, chromatin interactions, and gene dysregulation

Inherited peripheral neuropathies (IPNs) are a clinically and genetically heterogeneous group of diseases affecting the motor and sensory peripheral nerves. IPNs have benefited from gene discovery and genetic diagnosis using next-generation sequencing with over 80 causative genes available for testing. Despite this success, up to 50% of cases remain genetically unsolved. In the absence of protein coding mutations, noncoding DNA or structural variation (SV) mutations are a possible explanation. The most common IPN, Charcot-Marie-Tooth neuropathy type 1A (CMT1A), is caused by a 1.5 Mb duplication causing trisomy of the dosage sensitive gene PMP22. Using genome sequencing, we recently identified two large genomic rearrangements causing IPN subtypes X-linked CMT (CMTX3) and distal hereditary motor neuropathy (DHMN1), thereby expanding the spectrum of SV mutations causing IPN. Understanding how newly discovered SVs can cause IPN may serve as a useful paradigm to examine the role of topologically associated domains (TADs), chromatin interactions, and gene dysregulation in disease. This review will describe the growing role of SV in the pathogenesis of IPN and the importance of considering this type of mutation in Mendelian diseases where protein coding mutations cannot be identified.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

Digging for known genetic mutations underlying inherited bone and cartilage characteristics and disorders in the dog and cat

Gene mapping projects for many traits in both dogs and cats have yielded new knowledge. Both researchers and the public alike have been fascinated by the inheritance of breed characteristic phenotypes and sporadic disorders. It has been proposed that selective breeding practices have on occasion generated alterations in structure that might be harmful. In this review, simply inherited disorders and characteristics affecting bone and cartilage for which a putative mutation is known are collected. A better understanding of the known inherited basis of skeletal conditions and disorders will assist veterinarians to improve their diagnoses and increase their effectiveness on advising clients on the prevention, management, prognosis and possible treatment of the conditions.

Functional monosomy of 6q27-qter and functional disomy of Xpter-p22.11 due to X;6 translocation with an atypical X-inactivation pattern

Pattern of X chromosome inactivation (XCI) is typically random in females. However, chromosomal rearrangements affecting the X chromosome can result in XCI skewing due to cell growth disadvantage. In case of an X;autosome translocation, this usually leads to an XCI pattern of 100:0 with the derivative X being the active one in the majority of females. A de novo balanced X;6 translocation [46,X,t(X;6)(p22.1;q27)] and a completely skewed XCI pattern (100:0) were detected in a female patient with microcephaly, cerebellar vermis hypoplasia, heart defect, and severe developmental delay. We mapped the breakpoint regions using fluorescence in situ hybridization and found the X-linked gene POLA1 to be disrupted. POLA1 codes for the catalytic subunit of the polymerase α-primase complex which is responsible for initiation of the DNA replication process; absence of POLA1 is probably incompatible with life. Consequently, by RBA banding we determined which of the X chromosomes was the active one in the patient. In all examined lymphocytes the wild-type X chromosome was active. We propose that completely skewed XCI favoring the normal X chromosome resulted from death of cells with an active derivative X that was caused by a non-functional POLA1 gene. In summary, we conclude that functional monosomy of 6q27-qter and functional disomy of Xpter-p22.11 are responsible for the clinical phenotype of the patient. This case demonstrates the importance of determining which one of the X chromosomes underwent inactivation to correlate clinical features of a female with an X;autosome translocation with the nature of the genetic alteration.

Distant eQTLs and Non-coding Sequences Play Critical Roles in Regulating Gene Expression and Quantitative Trait Variation in Maize

A detailed understanding of genetic architecture of mRNA expression by millions of genetic variants is important for studying quantitative trait variation. In this study, we identified 1.25M SNPs with a minor allele frequency greater than 0.05 by combining reduced genome sequencing (GBS), high-density array technologies (600K), and previous deep RNA-sequencing data from 368 diverse inbred lines of maize. The balanced allelic frequencies and distributions in a relatively large and diverse natural panel helped to identify expression quantitative trait loci (eQTLs) associated with more than 18 000 genes (63.4% of tested genes). We found that distant eQTLs were more frequent (∼75% of all eQTLs) across the whole genome. Thirteen novel associated loci affecting maize kernel oil concentration were identified using the new dataset, among which one intergenic locus affected the kernel oil variation by controlling expression of three other known oil-related genes. Altogether, this study provides resources for expanding our understanding of cellular regulatory mechanisms of transcriptome variation and the landscape of functional variants within the maize genome, thereby enhancing the understanding of quantitative variations.

Decoding the role of regulatory element polymorphisms in complex disease

Genetic variation in gene regulatory elements contributes to diverse human diseases, ranging from rare and severe developmental defects to common and complex diseases such as obesity and diabetes. Early examples of regulatory mechanisms of human diseases involve large chromosomal rearrangements that change the regulatory connections within the genome. Single nucleotide variants in regulatory elements can also contribute to disease, potentially via demonstrated associations with changes in transcription factor binding, enhancer activity, post-translational histone modifications, long-range enhancer-promoter interactions, or RNA polymerase recruitment. Establishing causality between non-coding genetic variants, gene regulation, and disease has recently become more feasible with advances in genome-editing and epigenome-editing technologies. As establishing causal regulatory mechanisms of diseases becomes routine, functional annotation of target genes is likely to emerge as a major bottleneck for translation into patient benefits. In this review, we discuss the history and recent advances in understanding the regulatory mechanisms of human disease, and new challenges likely to be encountered once establishing those mechanisms becomes rote.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

Functional characterization of open chromatin in bidirectional promoters of rice

Bidirectional gene pairs tend to be highly coregulated and function in similar biological processes in eukaryotic genomes. Structural features and functional consequences of bidirectional promoters (BDPs) have received considerable attention among diverse species. However, the underlying mechanisms responsible for the bidirectional transcription and coexpression of BDPs remain poorly understood in plants. In this study, we integrated DNase-seq, RNA-seq, ChIP-seq and MNase-seq data and investigated the effect of physical DNase I hypersensitive site (DHS) positions on the transcription of rice BDPs. We found that the physical position of a DHS relative to the TSS of bidirectional gene pairs can affect the expression of the corresponding genes: the closer a DHS is to the TSS, the higher is the expression level of the genes. Most importantly, we observed that the distribution of DHSs plays a significant role in the regulation of transcription and the coexpression of gene pairs, which are possibly mediated by orchestrating the positioning of histone marks and canonical nucleosomes around BDPs. Our results demonstrate that the combined actions of chromatin structures with DHSs, which contain functional cis-elements for interaction with transcriptional machinery, may play an important role in the regulation of the bidirectional transcription or coexpression in rice BDPs. Our findings may help to enhance the understanding of DHSs in the regulation of bidirectional gene pairs.

Integrated small copy number variations and epigenome maps of disorders of sex development

Small copy number variations (CNVs) have typically not been analyzed or reported in clinical settings and hence have remained underrepresented in databases and the literature. Here, we focused our investigations on these small CNVs using chromosome microarray analysis (CMA) data previously obtained from patients with atypical characteristics or disorders of sex development (DSD). Using our customized CMA track targeting 334 genes involved in the development of urogenital and reproductive structures and a less stringent analysis filter, we uncovered small genes with recurrent and overlapping CNVs as small as 1 kb, and small regions of homozygosity (ROHs), imprinting and position effects. Detailed analysis of these high-resolution data revealed CNVs and ROHs involving structural and functional domains, repeat elements, active transcription sites and regulatory regions. Integration of these genomic data with DNA methylation, histone modification and predicted RNA expression profiles in normal testes and ovaries suggested spatiotemporal and tissue-specific gene regulation. This study emphasized a DSD-specific and gene-targeted CMA approach that uncovered previously unanalyzed or unreported small genes and CNVs, contributing to the growing resources on small CNVs and facilitating the narrowing of the genomic gap for identifying candidate genes or regions. This high-resolution analysis tool could improve the diagnostic utility of CMA, not only in patients with DSD but also in other clinical populations. These integrated data provided a better genomic-epigenomic landscape of DSD and greater opportunities for downstream research.

Three Dimensional Organization of Genome Might Have Guided the Dynamics of Gene Order Evolution in Eukaryotes

In eukaryotes, genes are nonrandomly organized into short gene-dense regions or "gene-clusters" interspersed by long gene-poor regions. How these gene-clusters have evolved is not entirely clear. Gene duplication may not account for all the gene-clusters since the genes in most of the clusters do not exhibit significant sequence similarity. In this study, using genome-wide data sets from budding yeast, fruit-fly, and human, we show that: 1) long-range evolutionary repositioning of genes strongly associate with their spatial proximity in the nucleus; 2) presence of evolutionary DNA break-points at involved loci hints at their susceptibility to undergo long-range genomic rearrangements; and 3) correlated epigenetic and transcriptional states of engaged genes highlight the underlying evolutionary constraints. The significance of observation 1, 2, and 3 are particularly stronger for the instances of inferred evolutionary gain, as compared with loss, of linear gene-clustering. These observations suggest that the long-range genomic rearrangements guided through 3D genome organization might have contributed to the evolution of gene order. We further hypothesize that the evolution of linear gene-clusters in eukaryotic genomes might have been mediated through spatial interactions among distant loci in order to optimize co-ordinated regulation of genes. We model this hypothesis through a heuristic model of gene-order evolution.

Generating CRISPR/Cas9 Mediated Monoallelic Deletions to Study Enhancer Function in Mouse Embryonic Stem Cells

Enhancers control cell identity by regulating tissue-specific gene expression in a position and orientation independent manner. These enhancers are often located distally from the regulated gene in intergenic regions or even within the body of another gene. The position independent nature of enhancer activity makes it difficult to match enhancers with the genes they regulate. Deletion of an enhancer region provides direct evidence for enhancer activity and is the gold standard to reveal an enhancer's role in endogenous gene transcription. Conventional homologous recombination based deletion methods have been surpassed by recent advances in genome editing technology which enable rapid and precisely located changes to the genomes of numerous model organisms. CRISPR/Cas9 mediated genome editing can be used to manipulate the genome in many cell types and organisms rapidly and cost effectively, due to the ease with which Cas9 can be targeted to the genome by a guide RNA from a bespoke expression plasmid. Homozygous deletion of essential gene regulatory elements might lead to lethality or alter cellular phenotype whereas monoallelic deletion of transcriptional enhancers allows for the study of cis-regulation of gene expression without this confounding issue. Presented here is a protocol for CRISPR/Cas9 mediated deletion in F1 mouse embryonic stem (ES) cells (Mus musculus(129) x Mus castaneus). Monoallelic deletion, screening and expression analysis is facilitated by single nucleotide polymorphisms (SNP) between the two alleles which occur on average every 125 bp in these cells.

The segregation of different submicroscopic imbalances underlying the clinical variability associated with a familial karyotypically balanced translocation

Background

About 7 % of karyotypically balanced chromosomal rearrangements (BCRs) are associated with congenital anomalies due to gene or regulatory element disruption, and cryptic imbalances on rearranged chromosomes. Rare familial BCRs segregating with clinical features are a powerful source for the identifying of causative genes due to the presence of several affected carriers.

Case presentation

We report on a karyotypically balanced translocation t(2;22)(p13;q12.2) associated with variable learning disabilities, and craniofacial and hand dysmorphisms, detected in six individuals in a three-generation family. Combined a-CGH, FISH and mate-pair sequencing revealed a ten-break complex rearrangement, also involving chromosome 5. As the consequence of the segregation of the derivative chromosomes der(2), der(5) and der(22), different imbalances were present in affected and clinically normal family members, thus contributing to the clinical variability. A 6.64 Mb duplication of a 5q23.2-23.3 segment was the imbalance common to all affected individuals. Although LMNB1, implicated in adult-onset autosomal dominant leukodystrophy (ADLD) when overexpressed, was among the 18 duplicated genes, none of the adult carriers manifested ADLD, and LMNB1 overexpression was not detected in the two tested individuals, after qRT-PCR. The ectopic location of the extra copy of the LMBN1 gene on chromosome 22 might have negatively impacted its expression. In addition, two individuals presenting with more severe learning disabilities carried a 1.42 Mb 2p14 microdeletion, with three genes (CEP68, RAB1A and ACTR2),which are candidates for the intellectual impairment observed in the previously described 2p14p15 microdeletion syndrome, mapping to the minimal overlapping deleted segment. A 5p15.1 deletion, encompassing 1.47 Mb, also detected in the family, did not segregate with the clinical phenotype.

Conclusion

The disclosing of the complexity of an apparently simple two-break familial rearrangement illustrates the importance of reconstructing the precise structure of derivative chromosomes for establishing genotype-phenotype correlations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13039-015-0205-9) contains supplementary material, which is available to authorized users.

High throughput technologies for the functional discovery of mammalian enhancers: new approaches for understanding transcriptional regulatory network dynamics

Completion of the human and mouse genomes has inspired new initiatives to obtain a global understanding of the functional regulatory networks governing gene expression. Enhancers are primary regulatory DNA elements determining precise spatio- and temporal gene expression patterns, but the observation that they can function at any distance from the gene(s) they regulate has made their genome-wide characterization challenging. Since traditional, single reporter approaches would be unable to accomplish this enormous task, high throughput technologies for mapping chromatin features associated with enhancers have emerged as an effective surrogate for enhancer discovery. However, the last few years have witnessed the development of several new innovative approaches that can effectively screen for and discover enhancers based on their functional activation of transcription using massively parallel reporter systems. In addition to their application for genome annotation, these new high throughput functional approaches open new and exciting avenues for modeling gene regulatory networks.

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.

A novel acquired cryptic three-way translocation t(2;11;5)(p21.3;q13.5;q23.2) with a submicroscopic deletion at 11p14.3 in an adult with hypereosinophilic syndrome

Hypereosinophilic syndrome (HES) is a clinically and pathologically heterogeneous disease entity. It is characterized by persistent eosinophilia and organ damage after excluding other causes. Clonal eosinophilia is distinguished from idiopathic eosinophilia by the presence of histologic, cytogenetic, or molecular evidence of an underlying malignancy. There are two distinct subcategories of clonal eosinophilia: chronic eosinophilic leukemia, not otherwise specified and myeloid/lymphoid neoplasms with eosinophilia and mutations involving platelet-derived growth factor receptor α/β or fibroblast growth factor receptor 1. More than 50% of HES are without knowledge of underlying pathogenic molecular pathways. Here we examined a HES patient by oligo-based aCGH analysis and molecular cytogenetic methods. Examination for the common eosinophilia-related cytogenetic abnormalities involving the genes PDGFRA, PDGFRB, and FGFR1 together with BCR-ABL fusion gene was negative. Cytogenetic analysis and multi-color FISH analysis revealed a novel cryptic three-way translocation t(2;11;5)(p21.3;q13.5;q23.2). By oaCGH analysis we could not find any copy number changes related to the cytogenetic breakpoints but instead detected a 0.9Mb submicroscopic deletion at 11p14.3. The deleted region involved the 5'-upstream sequences and exons 1-4 of the LUZP2 gene, which encodes a leucine zipper protein. Analysis of surrogate germ-line cells revealed a normal result showing that the detected chromosomal aberrations were acquired. This is the first report on a HES patient associated with a novel complex three-way translocation t(2;11;5)(p21.3;q13.5;q23.2) and a submicroscopic deletion in chromosome band 11p14.3. The study also demonstrates the benefits of oligo-based aCGH analysis in detecting hidden disease related chromosomal abnormalities. The present findings provide additional clues to unravel important molecular pathways in HES to obtain the full spectrum of acquired chromosomal and genomic aberrations in this heterogeneous disease entity. As more cases become characterized this may eventually improve on classification and treatment options.

Using the ENCODE Resource for Functional Annotation of Genetic Variants

This article illustrates the use of the Encyclopedia of DNA Elements (ENCODE) resource to generate or refine hypotheses from genomic data on disease and other phenotypic traits. First, the goals and history of ENCODE and related epigenomics projects are reviewed. Second, the rationale for ENCODE and the major data types used by ENCODE are briefly described, as are some standard heuristics for their interpretation. Third, the use of the ENCODE resource is examined. Standard use cases for ENCODE, accessing the ENCODE resource, and accessing data from related projects are discussed. Although the focus of this article is the use of ENCODE data, some of the same approaches can be used with data from other projects.

Evaluation of copy number variation and gene expression in neurofibromatosis type-1-associated malignant peripheral nerve sheath tumours

BACKGROUND

Neurofibromatosis type-1 (NF1) is a complex neurogenetic disorder characterised by the development of benign and malignant tumours of the peripheral nerve sheath (MPNSTs). Whilst biallelic NF1 gene inactivation contributes to benign tumour formation, additional cellular changes in gene structure and/or expression are required to induce malignant transformation. Although few molecular profiling studies have been performed on the process of progression of pre-existing plexiform neurofibromas to MPNSTs, the integrated analysis of copy number alterations (CNAs) and gene expression is likely to be key to understanding the molecular mechanisms underlying NF1-MPNST tumorigenesis. In a pilot study, we employed this approach to identify genes differentially expressed between benign and malignant NF1 tumours.

RESULTS

SPP1 (osteopontin) was the most differentially expressed gene (85-fold increase in expression), compared to benign plexiform neurofibromas. Short hairpin RNA (shRNA) knockdown of SPP1 in NF1-MPNST cells reduced tumour spheroid size, wound healing and invasion in four different MPNST cell lines. Seventy-six genes were found to exhibit concordance between CNA and gene expression level.

CONCLUSIONS

Pathway analysis of these genes suggested that glutathione metabolism and Wnt signalling may be specifically involved in NF1-MPNST development. SPP1 is associated with malignant transformation in NF1-associated MPNSTs and could prove to be an important target for therapeutic intervention.

Chromosome boundary elements and regulation of heterochromatin spreading

Chromatin is generally classified as euchromatin or heterochromatin, each with distinct histone modifications, compaction levels, and gene expression patterns. Although the proper formation of heterochromatin is essential for maintaining genome integrity and regulating gene expression, heterochromatin can also spread into neighboring regions in a sequence-independent manner, leading to the inactivation of genes. Because the distance of heterochromatin spreading is stochastic, the formation of boundaries, which block the spreading of heterochromatin, is critical for maintaining stable gene expression patterns. Here we review the current understanding of the mechanisms underlying heterochromatin spreading and boundary formation.

Structure-function relationships in human testis-determining factor SRY: an aromatic buttress underlies the specific DNA-bending surface of a high mobility group (HMG) box

Human testis determination is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in SRY cause 46 XY gonadal dysgenesis with female somatic phenotype (Swyer syndrome) and confer a high risk of malignancy (gonadoblastoma). Such mutations cluster in the SRY high mobility group (HMG) box, a conserved motif of specific DNA binding and bending. To explore structure-function relationships, we constructed all possible substitutions at a site of clinical mutation (W70L). Our studies thus focused on a core aromatic residue (position 15 of the consensus HMG box) that is invariant among SRY-related HMG box transcription factors (the SOX family) and conserved as aromatic (Phe or Tyr) among other sequence-specific boxes. In a yeast one-hybrid system sensitive to specific SRY-DNA binding, the variant domains exhibited reduced (Phe and Tyr) or absent activity (the remaining 17 substitutions). Representative nonpolar variants with partial or absent activity (Tyr, Phe, Leu, and Ala in order of decreasing side-chain volume) were chosen for study in vitro and in mammalian cell culture. The clinical mutation (Leu) was found to markedly impair multiple biochemical and cellular activities as respectively probed through the following: (i) in vitro assays of specific DNA binding and protein stability, and (ii) cell culture-based assays of proteosomal degradation, nuclear import, enhancer DNA occupancy, and SRY-dependent transcriptional activation. Surprisingly, however, DNA bending is robust to this or the related Ala substitution that profoundly impairs box stability. Together, our findings demonstrate that the folding, trafficking, and gene-regulatory function of SRY requires an invariant aromatic "buttress" beneath its specific DNA-bending surface.

A high resolution map of mammalian X chromosome fragile regions assessed by large-scale comparative genomics

Chromosomal evolution involves multiple changes at structural and numerical levels. These changes, which are related to the variation of the gene number and their location, can be tracked by the identification of syntenic blocks (SB). First reports proposed that ~180-280 SB might be shared by mouse and human species. More recently, further studies including additional genomes have identified up to ~1,400 SB during the evolution of eutherian species. A considerable number of studies regarding the X chromosome's structure and evolution have been undertaken because of its extraordinary biological impact on reproductive fitness and speciation. Some have identified evolutionary breakpoint regions and fragile sites at specific locations in the human X chromosome. However, mapping these regions to date has involved using low-to-moderate resolution techniques. Such scenario might be related to underestimating their total number and giving an inaccurate location. The present study included using a combination of bioinformatics methods for identifying, at base-pair level, chromosomal rearrangements occurring during X chromosome evolution in 13 mammalian species. A comparative technique using four different algorithms was used for optimizing the detection of hotspot regions in the human X chromosome. We identified a significant interspecific variation in SB size which was related to genetic information gain regarding the human X chromosome. We found that human hotspot regions were enriched by LINE-1 and Alu transposable elements, which may have led to intraspecific chromosome rearrangement events. New fragile regions located in the human X chromosome have also been postulated. We estimate that the high resolution map of X chromosome fragile sites presented here constitutes useful data concerning future studies on mammalian evolution and human disease.

Function-based identification of mammalian enhancers using site-specific integration

The accurate and comprehensive identification of functional regulatory sequences in mammalian genomes remains a major challenge. Here we describe site-specific integration fluorescence-activated cell sorting followed by sequencing (SIF-seq), an unbiased, medium-throughput functional assay for the discovery of distant-acting enhancers. Targeted single-copy genomic integration into pluripotent cells, reporter assays and flow cytometry are coupled with high-throughput DNA sequencing to enable parallel screening of large numbers of DNA sequences. By functionally interrogating >500 kilobases (kb) of mouse and human sequence in mouse embryonic stem cells for enhancer activity we identified enhancers at pluripotency loci including NANOG. In in vitro-differentiated cardiomyocytes and neural progenitor cells, we identified cardiac enhancers and neuronal enhancers, respectively. SIF-seq is a powerful and flexible method for de novo functional identification of mammalian enhancers in a potentially wide variety of cell types.

Defining functional DNA elements in the human genome

With the completion of the human genome sequence, attention turned to identifying and annotating its functional DNA elements. As a complement to genetic and comparative genomics approaches, the Encyclopedia of DNA Elements Project was launched to contribute maps of RNA transcripts, transcriptional regulator binding sites, and chromatin states in many cell types. The resulting genome-wide data reveal sites of biochemical activity with high positional resolution and cell type specificity that facilitate studies of gene regulation and interpretation of noncoding variants associated with human disease. However, the biochemically active regions cover a much larger fraction of the genome than do evolutionarily conserved regions, raising the question of whether nonconserved but biochemically active regions are truly functional. Here, we review the strengths and limitations of biochemical, evolutionary, and genetic approaches for defining functional DNA segments, potential sources for the observed differences in estimated genomic coverage, and the biological implications of these discrepancies. We also analyze the relationship between signal intensity, genomic coverage, and evolutionary conservation. Our results reinforce the principle that each approach provides complementary information and that we need to use combinations of all three to elucidate genome function in human biology and disease.

Targeting of Slc25a21 is associated with orofacial defects and otitis media due to disrupted expression of a neighbouring gene

Homozygosity for Slc25a21^{tm1a(KOMP)Wtsi} results in mice exhibiting orofacial abnormalities, alterations in carpal and rugae structures, hearing impairment and inflammation in the middle ear. In humans it has been hypothesised that the 2-oxoadipate mitochondrial carrier coded by SLC25A21 may be involved in the disease 2-oxoadipate acidaemia. Unexpectedly, no 2-oxoadipate acidaemia-like symptoms were observed in animals homozygous for Slc25a21^{tm1a(KOMP)Wtsi} despite confirmation that this allele reduces Slc25a21 expression by 71.3%. To study the complete knockout, an allelic series was generated using the loxP and FRT sites typical of a Knockout Mouse Project allele. After removal of the critical exon and neomycin selection cassette, Slc25a21 knockout mice homozygous for the Slc25a21^{tm1b(KOMP)Wtsi} and Slc25a21^{tm1d(KOMP)Wtsi} alleles were phenotypically indistinguishable from wild-type. This led us to explore the genomic environment of Slc25a21 and to discover that expression of Pax9, located 3′ of the target gene, was reduced in homozygous Slc25a21^{tm1a(KOMP)Wtsi} mice. We hypothesize that the presence of the selection cassette is the cause of the down regulation of Pax9 observed. The phenotypes we observed in homozygous Slc25a21^{tm1a(KOMP)Wtsi} mice were broadly consistent with a hypomorphic Pax9 allele with the exception of otitis media and hearing impairment which may be a novel consequence of Pax9 down regulation. We explore the ramifications associated with this particular targeted mutation and emphasise the need to interpret phenotypes taking into consideration all potential underlying genetic mechanisms.

Megabladder mouse model of congenital obstructive nephropathy: genetic etiology and renal adaptation

Congenital obstructive nephropathy remains one of the leading causes of chronic renal failure in children. The direct link between obstructed urine flow and abnormal renal development and subsequent dysfunction represents a central paradigm of urogenital pathogenesis that has far-reaching clinical implications. Even so, a number of diagnostic, prognostic, and therapeutic quandaries still exist in the management of congenital obstructive nephropathy. Studies in our laboratory have characterized a unique mutant mouse line that develops in utero megabladder, variable hydronephrosis, and progressive renal failure. Megabladder mice represent a valuable functional model for the study of congenital obstructive nephropathy. Recent studies have begun to shed light on the genetic etiology of mgb (-/-) mice as well as the molecular pathways controlling disease progression in these animals.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish

Many important model organisms for biomedical and evolutionary research have sequenced genomes, but occupy a phylogenetically isolated position, evolutionarily distant from other sequenced genomes. This phylogenetic isolation is exemplified for zebrafish, a vertebrate model for cis-regulation, development and human disease, whose evolutionary distance to all other currently sequenced fish exceeds the distance between human and chicken. Such large distances make it difficult to align genomes and use them for comparative analysis beyond gene-focused questions. In particular, detecting conserved non-genic elements (CNEs) as promising cis-regulatory elements with biological importance is challenging. Here, we develop a general comparative genomics framework to align isolated genomes and to comprehensively detect CNEs. Our approach integrates highly sensitive and quality-controlled local alignments and uses alignment transitivity and ancestral reconstruction to bridge large evolutionary distances. We apply our framework to zebrafish and demonstrate substantially improved CNE detection and quality compared with previous sets. Our zebrafish CNE set comprises 54 533 CNEs, of which 11 792 (22%) are conserved to human or mouse. Our zebrafish CNEs (http://zebrafish.stanford.edu) are highly enriched in known enhancers and extend existing experimental (ChIP-Seq) sets. The same framework can now be applied to the isolated genomes of frog, amphioxus, Caenorhabditis elegans and many others.

Enhancers: five essential questions

It is estimated that the human genome contains hundreds of thousands of enhancers, so understanding these gene-regulatory elements is a crucial goal. Several fundamental questions need to be addressed about enhancers, such as how do we identify them all, how do they work, and how do they contribute to disease and evolution? Five prominent researchers in this field look at how much we know already and what needs to be done to answer these questions.

Functional anatomy of distant-acting mammalian enhancers

Transcriptional enhancers are a major class of functional element embedded in the vast non-coding portion of the human genome. Acting over large genomic distances, enhancers play critical roles in the tissue and cell type-specific regulation of genes, and there is mounting evidence that they contribute to the aetiology of many human diseases. Methods for genome-wide mapping of enhancer regions are now available, but the functional architecture contained within human enhancer elements remains unclear. Here, we review recent approaches aimed at understanding the functional anatomy of individual enhancer elements, using systematic qualitative and quantitative assessments of mammalian enhancer variants in cultured cells and in vivo. These studies provide direct insight into common architectural characteristics of enhancers including the presence of multiple transcription factor-binding sites and the mixture of both transcriptionally activating and repressing domains within the same enhancer. Despite such progress in understanding the functional composition of enhancers, the inherent complexities of enhancer anatomy continue to limit our ability to predict the impact of sequence changes on in vivo enhancer function. While providing an initial glimpse into the mutability of mammalian enhancers, these observations highlight the continued need for experimental enhancer assessment as genome sequencing becomes routine in the clinic.

Long-range transcriptional regulation of breast cancer genes

Breast cancer is a major health problem and understanding the genetic basis of this disease is crucial for predicting risk and developing effective targeted therapeutics. Several breast cancer predisposing genes have been identified, but mutations in the coding regions of these genes only accounts for a small proportion of risk. Research now suggests that combinations of multiple non-coding changes in breast cancer susceptibility genes, which cause moderate alterations in gene expression, will be responsible for the remaining inherited risk. These non-coding changes will include variants in proximal and distal transcriptional and post-transcriptional regulatory elements and may affect the levels and function of trans-acting factors, including proteins and RNAs, which act on these elements. Somatic changes in such elements and factors have also been associated with breast cancer progression. With the recent advent of techniques allowing the detection of long-range DNA interactions spanning the human genome, it has become increasingly clear that long-range regulatory elements constitute an important mechanism for gene regulation. Recent studies have identified several such elements that are important for regulating genes involved in breast cancer, raising the possibility that defects in these sequences may contribute to breast cancer predisposition and progression. In this review, we discuss the emerging functions of cis-regulatory elements and a subset of trans-acting factors in breast tumorigenesis. We also discuss some recent progress in our understanding of how dysregulation in these transcriptional components may contribute to breast cancer, and the potential implications for molecular diagnosis, prognosis prediction, and the treatment of this disease.

Transcription factor binding at enhancers: shaping a genomic regulatory landscape in flux

The mammalian genome is packed tightly in the nucleus of the cell. This packing is primarily facilitated by histone proteins and results in an ordered organization of the genome in chromosome territories that can be roughly divided in heterochromatic and euchromatic domains. On top of this organization several distinct gene regulatory elements on the same chromosome or other chromosomes are thought to dynamically communicate via chromatin looping. Advances in genome-wide technologies have revealed the existence of a plethora of these regulatory elements in various eukaryotic genomes. These regulatory elements are defined by particular in vitro assays as promoters, enhancers, insulators, and boundary elements. However, recent studies indicate that the in vivo distinction between these elements is often less strict. Regulatory elements are bound by a mixture of common and lineage-specific transcription factors which mediate the long-range interactions between these elements. Inappropriate modulation of the binding of these transcription factors can alter the interactions between regulatory elements, which in turn leads to aberrant gene expression with disease as an ultimate consequence. Here we discuss the bi-modal behavior of regulatory elements that act in cis (with a focus on enhancers), how their activity is modulated by transcription factor binding and the effect this has on gene regulation.

Human genetics and genomics a decade after the release of the draft sequence of the human genome

Substantial progress has been made in human genetics and genomics research over the past ten years since the publication of the draft sequence of the human genome in 2001. Findings emanating directly from the Human Genome Project, together with those from follow-on studies, have had an enormous impact on our understanding of the architecture and function of the human genome. Major developments have been made in cataloguing genetic variation, the International HapMap Project, and with respect to advances in genotyping technologies. These developments are vital for the emergence of genome-wide association studies in the investigation of complex diseases and traits. In parallel, the advent of high-throughput sequencing technologies has ushered in the 'personal genome sequencing' era for both normal and cancer genomes, and made possible large-scale genome sequencing studies such as the 1000 Genomes Project and the International Cancer Genome Consortium. The high-throughput sequencing and sequence-capture technologies are also providing new opportunities to study Mendelian disorders through exome sequencing and whole-genome sequencing. This paper reviews these major developments in human genetics and genomics over the past decade.

Spatial partitioning of the regulatory landscape of the X-inactivation centre

In eukaryotes transcriptional regulation often involves multiple long-range elements and is influenced by the genomic environment. A prime example of this concerns the mouse X-inactivation centre (Xic), which orchestrates the initiation of X-chromosome inactivation (XCI) by controlling the expression of the non-protein-coding Xist transcript. The extent of Xic sequences required for the proper regulation of Xist remains unknown. Here we use chromosome conformation capture carbon-copy (5C) and super-resolution microscopy to analyse the spatial organization of a 4.5-megabases (Mb) region including Xist. We discover a series of discrete 200-kilobase to 1 Mb topologically associating domains (TADs), present both before and after cell differentiation and on the active and inactive X. TADs align with, but do not rely on, several domain-wide features of the epigenome, such as H3K27me3 or H3K9me2 blocks and lamina-associated domains. TADs also align with coordinately regulated gene clusters. Disruption of a TAD boundary causes ectopic chromosomal contacts and long-range transcriptional misregulation. The Xist/Tsix sense/antisense unit illustrates how TADs enable the spatial segregation of oppositely regulated chromosomal neighbourhoods, with the respective promoters of Xist and Tsix lying in adjacent TADs, each containing their known positive regulators. We identify a novel distal regulatory region of Tsix within its TAD, which produces a long intervening RNA, Linx. In addition to uncovering a new principle of cis-regulatory architecture of mammalian chromosomes, our study sets the stage for the full genetic dissection of the X-inactivation centre.

A genetic approach to the transcriptional regulation of Hox gene clusters

The evolution of vertebrate genomes was accompanied by an astounding increase in the complexity of their regulatory modalities. Genetic redundancy resulting from large-scale genome duplications at the base of the chordate tree was repeatedly exploited by the functional redeployment of paralogous genes via innovations in their regulatory circuits. As a paradigm of such regulatory evolution, we have extensively studied those control mechanisms at work in-cis over vertebrate Hox gene clusters. Here, we review the portfolio of genetic strategies that have been developed to tackle the intricate relationship between genomic topography and the transcriptional activities in this gene family, and we describe some of the mechanistic insights we gained by using the HoxD cluster as an example. We discuss the high heuristic value of this system in our general understanding of how changes in transcriptional regulation can diversify gene function and thereby fuel morphological evolution.

[Cis-ruptions of highly conserved non-coding genomic elements distant from the SOX9 gene in the Pierre Robin sequence]

Major developmental genes, exhibiting complex expression patterns, are often embedded within a genic desert particularly rich in regions, which though non-coding are highly conserved. The developmental expression of these genes in many areas requires coordinated regulation in time and space, which is orchestrated by some of these conserved non-coding regions, acting as transcriptional regulators. SOX9 is an essential gene for many developmental processes, such as chondrogenesis, migration and differentiation of neural crest cells and testis development. In agreement with these major expression areas, SOX9 haploinsufficiency, linked to alterations in coding sequence, leads to a polymorphic malformation syndrome - campomelic dysplasia - whose major symptoms are a bone anomaly, a Pierre Robin sequence, and a sexual differentiation anomaly (Disorder of Sex Development, DSD). SOX9 is located in a  ~2.5 Mb gene desert extremely rich in conserved sequences. We have used the SOX9 locus and campomelic dysplasia as a model to show that one or several endophenotypes within a complex syndrome may arise from a tissue-specific deregulation of a major developmental gene transcription. Our work has focused on one of these endophenotypes, SPR, characterized by the triad micro- and/or retrognathy, glossoptosis and cleft palate. Here we report in detail how we identified alterations (translocations, deletions, point mutations) in non-coding regions, located far away (more than 1.2 Mb) upstream and downstream of SOX9, in clustered or sporadic SPR cases.

Histone crosstalk directed by H2B ubiquitination is required for chromatin boundary integrity

Genomic maps of chromatin modifications have provided evidence for the partitioning of genomes into domains of distinct chromatin states, which assist coordinated gene regulation. The maintenance of chromatin domain integrity can require the setting of boundaries. The HS4 insulator element marks the 3′ boundary of a heterochromatin region located upstream of the chicken β-globin gene cluster. Here we show that HS4 recruits the E3 ligase RNF20/BRE1A to mediate H2B mono-ubiquitination (H2Bub1) at this insulator. Knockdown experiments show that RNF20 is required for H2Bub1 and processive H3K4 methylation. Depletion of RNF20 results in a collapse of the active histone modification signature at the HS4 chromatin boundary, where H2Bub1, H3K4 methylation, and hyperacetylation of H3, H4, and H2A.Z are rapidly lost. A remarkably similar set of events occurs at the HSA/HSB regulatory elements of the FOLR1 gene, which mark the 5′ boundary of the same heterochromatin region. We find that persistent H2Bub1 at the HSA/HSB and HS4 elements is required for chromatin boundary integrity. The loss of boundary function leads to the sequential spreading of H3K9me2, H3K9me3, and H4K20me3 over the entire 50 kb FOLR1 and β-globin region and silencing of FOLR1 expression. These findings show that the HSA/HSB and HS4 boundary elements direct a cascade of active histone modifications that defend the FOLR1 and β-globin gene loci from the pervasive encroachment of an adjacent heterochromatin domain. We propose that many gene loci employ H2Bub1-dependent boundaries to prevent heterochromatin spreading.

The transcription of genes in eukaryotes occurs within the context of chromatin, a complex of DNA, histone proteins, and regulatory factors. Whole-genome profiling of chromatin proteins and histones that are post-translationally modified has revealed that genomes are organized into domains of distinct chromatin states that coordinate gene regulation. The integrity of chromatin domains can require the setting of their boundaries. DNA sequences known as chromatin insulator or boundary elements can establish boundaries between transcriptionally permissive and repressive chromatin domains. We have studied two chromatin boundary elements that flank a condensed chromatin region located between the chicken FOLR1 and β-globin genes, respectively. These elements recruit enzymes that mediate the ubiquitination of histone H2B. Histone H2B ubiquitination directs a cascade of so-called “active” histone modification events that favor chromatin accessibility. We observe a striking collapse of the active histone modification signature at both chromatin boundaries following the depletion of ubiquitinated H2B. This loss of boundary function leads to the comprehensive spreading of repressive chromatin over the entire FOLR1 and β-globin gene region, resulting in gene silencing. We propose that chromatin boundaries at many gene loci employ H2B ubiquitination to restrict the encroachment of repressive chromatin.

The core FOXG1 syndrome phenotype consists of postnatal microcephaly, severe mental retardation, absent language, dyskinesia, and corpus callosum hypogenesis

BACKGROUND

Submicroscopic deletions in 14q12 spanning FOXG1 or intragenic mutations have been reported in patients with a developmental disorder described as a congenital variant of Rett syndrome. This study aimed to further characterise and delineate the phenotype of FOXG1 mutation positive patients.

METHOD

The study mapped the breakpoints of a 2;14 translocation by fluorescence in situ hybridisation and analysed three chromosome rearrangements in 14q12 by cytogenetic analysis and/or array comparative genomic hybridisation. The FOXG1 gene was sequenced in 210 patients, including 129 patients with unexplained developmental disorders and 81 MECP2 mutation negative individuals.

RESULTS

One known mutation, seen in two patients, and nine novel mutations of FOXG1 including two deletions, two chromosome rearrangements disrupting or displacing putative cis-regulatory elements from FOXG1, and seven sequence changes, are reported. Analysis of 11 patients in this study, and a further 15 patients reported in the literature, demonstrates a complex constellation of features including mild postnatal growth deficiency, severe postnatal microcephaly, severe mental retardation with absent language development, deficient social reciprocity resembling autism, combined stereotypies and frank dyskinesias, epilepsy, poor sleep patterns, irritability in infancy, unexplained episodes of crying, recurrent aspiration, and gastro-oesophageal reflux. Brain imaging studies reveal simplified gyral pattern and reduced white matter volume in the frontal lobes, corpus callosum hypogenesis, and variable mild frontal pachgyria.

CONCLUSIONS

These findings have significantly expanded the number of FOXG1 mutations and identified two affecting possible cis-regulatory elements. While the phenotype of the patients overlaps both classic and congenital Rett syndrome, extensive clinical evaluation demonstrates a distinctive and clinically recognisable phenotype which the authors suggest designating as the FOXG1 syndrome.

Revealing the human mutome

The number of known mutations in human nuclear genes, underlying or associated with human inherited disease, has now exceeded 100,000 in more than 3700 different genes (Human Gene Mutation Database). However, for a variety of reasons, this figure is likely to represent only a small proportion of the clinically relevant genetic variants that remain to be identified in the human genome (the 'mutome'). With the advent of next-generation sequencing, we are currently witnessing a revolution in medical genetics. In particular, whole-genome sequencing (WGS) has the potential to identify all disease-causing or disease-associated DNA variants in a given individual. Here, we use examples of recent advances in our understanding of mutational/pathogenic mechanisms to guide our thinking about possible locations outwith gene-coding sequences for those disease-causing or disease-associated variants that are likely so often to have been overlooked because of the inadequacy of current mutation screening protocols. Such considerations are important not only for improving mutation-screening strategies but also for enhancing the interpretation of findings derived from genome-wide association studies, whole-exome sequencing and WGS. An improved understanding of the human mutome will not only lead to the development of improved diagnostic testing procedures but should also improve our understanding of human genome biology.