Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis.

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges.

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Digital genotyping of macrosatellites and multicopy genes reveals novel biological functions associated with copy number variation of large tandem repeats.

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5-10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality 'finished' human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed "repeat induced gene silencing", which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.

Electron density shape analysis of a family of through-space and through-bond interactions.

A family of styrene derivatives has been used to study the effects of through-space and through-bond interactions on the local and global shapes of electron densities of complete molecules and a set of substituents on their central rings. Shape analysis methods which have been used extensively in the past for the study of molecular property-molecular shape correlations have shown that in these molecules a complementary role is played by the through-space and through-bond interactions. For each specific example, the dominance of either one of the two interactions can be identified and interpreted in terms of local shapes and the typical reactivities of the various substituents. Three levels of quantum chemical computational methods have been applied for these structures, including the B3LYP/cc-pVTZ level of density functional methodology, and the essential conclusions are the same for all three levels. The general approach is suggested as a tool for the identification of specific interaction types which are able to modify molecular electron densities. By separately influencing the through-space and through-bond components using polar groups and groups capable of conjugation, some fine-tuning of the overall effects becomes possible. The method described may contribute to an improved understanding and control of molecular properties involving complex interactions with a possible role in the emerging field of molecular design.

Expansion of a 12-kb VNTR containing the REXO1L1 gene cluster underlies the microscopically visible euchromatic variant of 8q21.2.

Copy number variants visible with the light microscope have been described as euchromatic variants (EVs) and EVs with extra G-light material at 8q21.2 have been reported only once before. We report four further patients with EVs of 8q21.2 ascertained for clinical (3) or reproductive reasons (1). Enhanced signal strength from two overlapping bacterial artificial chromosomes (BACs) and microarray analysis mapped the EV to a 284-kb interval in the reference genome. This interval consists of a sequence gap flanked by segmental duplications that contain the 12-kb components of one of the largest Variable Number Tandem Repeat arrays in the human genome. Using digital NanoString technology with a custom probe for the RNA exonuclease 1 homologue (S. cerevisiae)-like 1 (REXO1L1) gene within each 12-kb repeat, significantly enhanced diploid copy numbers of 270 and 265 were found in an EV family and a median diploid copy number of 166 copies in 216 controls. These 8q21.2 EVs are not thought to have clinical consequences as the phenotypes of the probands were inconsistent, those referred for reproductive reasons were otherwise phenotypically normal and the REXO1L1 gene has no known disease association. This EV was found in 4/3078 (1 in 770) consecutive referrals for chromosome analysis and needs to be distinguished from pathogenic imbalances of medial 8q. The REXO1L1 gene product is a marker of hepatitis C virus (HCV) infection and a possible association between REXO1L1 copy number and susceptibility to HCV infection, progression or response to treatment has not yet been excluded.

The octamer is the major form of CENP-A nucleosomes at human centromeres.

The centromere is the chromosomal locus that ensures fidelity in genome transmission at cell division. Centromere protein A (CENP-A) is a histone H3 variant that specifies centromere location independently of DNA sequence. Conflicting evidence has emerged regarding the histone composition and stoichiometry of CENP-A nucleosomes. Here we show that the predominant form of the CENP-A particle at human centromeres is an octameric nucleosome. CENP-A nucleosomes are very highly phased on α-satellite 171-base-pair monomers at normal centromeres and also display strong positioning at neocentromeres. At either type of functional centromere, CENP-A nucleosomes exhibit similar DNA-wrapping behavior, as do octameric CENP-A nucleosomes reconstituted with recombinant components, having looser DNA termini than those on conventional nucleosomes containing canonical histone H3. Thus, the fundamental unit of the chromatin that epigenetically specifies centromere location in mammals is an octameric nucleosome with loose termini.

Neocentric X-chromosome in a girl with Turner-like syndrome.

BACKGROUND: Neocentromeres are rare human chromosomal aberrations in which a new centromere has formed in a previously non-centromeric location. We report the finding of a structurally abnormal X chromosome with a neocentromere in a 15-year-old girl with clinical features suggestive of Turner syndrome, including short stature and primary amenorrhea. RESULT: G-banded chromosome analysis revealed a mosaic female karyotype involving two abnormal cell lines. One cell line (84% of analyzed metaphases) had a structurally abnormal X chromosome (duplication of the long arm and deletion of the short arm) and a normal X chromosome. The other cell line (16% of cells) exhibited monosomy X. C-banding studies were negative for the abnormal X chromosome. FISH analysis revealed lack of hybridization of the abnormal X chromosome with both the X centromere-specific probe and the "all human centromeres" probe, a pattern consistent with lack of the X chromosome endogenous centromere. A FISH study using an XIST gene probe revealed the presence of two XIST genes, one on each long arm of the iso(Xq), required for inactivation of the abnormal X chromosome. R-banding also demonstrated inactivation of the abnormal X chromosome. An assay for centromeric protein C (CENP-C) was positive on both the normal and the abnormal X chromosomes. The position of CENP-C in the abnormal X chromosome defined a neocentromere, which explains its mitotic stability. The karyotype is thus designated as 46,X,neo(X)(qter- > q12::q12- > q21.2- > neo- > q21.2- > qter)[42]/45,X[8], which is consistent with stigmata of Turner syndrome. The mother of this patient has a normal karyotype; however, the father was not available for study. CONCLUSION: To our knowledge, this is the first case of mosaic Turner syndrome involving an analphoid iso(Xq) chromosome with a proven neocentromere among 90 previously described cases with a proven neocentromere.

Tetrasomy 15q26: a distinct syndrome or Shprintzen-Goldberg syndrome phenocopy?

PURPOSE: The aim of this study was to characterize the clinical phenotype of patients with tetrasomy of the distal 15q chromosome in the form of a neocentric marker chromosome and to evaluate whether the phenotype represents a new clinical syndrome or is a phenocopy of Shprintzen-Goldberg syndrome. METHODS: We carried out comprehensive clinical evaluation of four patients who were identified with a supernumerary marker chromosome. The marker chromosome was characterized by G-banding, fluorescence in situ hybridization, single nucleotide polymorphism oligonucleotide microarray analysis, and immunofluorescence with antibodies to centromere protein C. RESULTS: The marker chromosomes were categorized as being neocentric with all showing tetrasomy for regions distal to 15q25 and the common region of overlap being 15q26→qter. CONCLUSION: Tetrasomy of 15q26 likely results in a distinct syndrome as the patients with tetrasomy 15q26 share a strikingly more consistent phenotype than do the patients with Shprintzen-Goldberg syndrome, who show remarkable clinical variation.

Comparisons of computational and experimental thermochemical properties of α-amino acids.

This study provides comprehensive benchmark calculations for the thermochemical properties of the common α-amino acids. Calculated properties include the proton affinity, gas-phase basicity, protonation entropy, ΔH°(acid), ΔG°(acid), and enthalpies of formation for the protonated and deprotonated α-amino acids. In order to determine the performance at various levels of theory, including density functional methods and composite methods, the calculated thermochemical properties are compared to experimental results. For all the common α-amino acids investigated, the thermochemical properties computed with the Gaussian-n theories were found to be quite consistent with each other in terms of mean absolute deviation from experiment. While all Gaussian-n theory values can serve as benchmarks, we focus on the G3MP2 values as it is the least resource-intensive of the Gaussian-n theories considered.

Duplication of the ZIC2 gene is not associated with holoprosencephaly.

Cytogenetic testing using genomic microarrays presents a clinical challenge when data regarding the phenotypic consequences of the genomic alteration are not available. We describe a chromosome 13q32.3 duplication discovered by microarray testing in a fetus with a prenatally detected apparently balanced de novo translocation 46,XY,t(2;13)(q37;q32). Microarray analysis on the fetal DNA showed duplications of 384 and 564 kb at the breakpoint regions on chromosomes 2q37.3 and 13q32.3, respectively. There were no disease-associated genes in the duplicated region on chromosome 2q37. The duplicated region on chromosome 13q contains the ZIC2 gene. Haploinsufficiency of ZIC2 is known to cause holoprosencephaly and other brain malformations. Studies in the mouse models have suggested that over expression of ZIC2 may also lead to brain malformations. Fetal MRI of the brain was normal and the family elected to continue the pregnancy. An apparently normal baby was born at term. At 3 months of age a physical exam showed no abnormalities and no developmental delay. This report shows that duplication of ZIC2 is not necessarily associated with brain malformations. We also describe the phenotype from four additional patients with duplications of the region of chromosome 13 containing ZIC2 and three previously described patients with supernumerary marker chromosomes derived from distal chromosome 13. None of the eight patients had holoprosencephaly or brain malformations, indicating that duplication of ZIC2 is not associated with brain anomalies. This information will be useful for counseling in other occurrences of this duplication identified by microarray.

Formation of novel CENP-A domains on tandem repetitive DNA and across chromosome breakpoints on human chromosome 8q21 neocentromeres.

Endogenous human centromeres form on megabase-sized arrays of tandemly repeated alpha satellite DNA. Human neocentromeres form epigenetically at ectopic sites devoid of alpha satellite DNA and permit analysis of centromeric DNA and chromatin organization. In this study, we present molecular cytogenetic and CENP-A chromatin immunoprecipitation (ChIP) on CHIP analyses of two neocentromeres that have formed in chromosome band 8q21 each with a unique DNA and CENP-A chromatin configuration. The first neocentromere was found on a neodicentric chromosome 8 with an inactivated endogenous centromere, where the centromeric activity and CENP-A domain were repositioned to band 8q21 on a large tandemly repeated DNA. This is the first example of a neocentromere forming on repetitive DNA, as all other mapped neocentromeres have formed on single copy DNA. Quantitative fluorescent in situ hybridization (FISH) analysis showed a 60% reduction in the alpha satellite array size at the inactive centromere compared to the active centromere on the normal chromosome 8. This neodicentric chromosome may provide insight into centromere inactivation and the role of tandem DNA in centromere structure. The second neocentromere was found on a neocentric ring chromosome that contained the 8q21 tandemly repeated DNA, although the neocentromere was localized to a different genomic region. Interestingly, this neocentromere is composed of two distinct CENP-A domains in bands 8q21 and 8q24, which are brought into closer proximity on the ring chromosome. This neocentromere suggests that chromosomal rearrangement and DNA breakage may be involved in neocentromere formation. These novel examples provide insight into the formation and structure of human neocentromeres.

Atoms and bonds in molecules from radial densities.

Radial densities are explored as an alternative method for partitioning the molecular density into atomic regions and bonding regions. The radial densities for atoms in molecules are similar to those of an isolated atom. The method may also provide an alternative to Bragg-Slater radii.

Large inverted repeats within Xp11.2 are present at the breakpoints of isodicentric X chromosomes in Turner syndrome.

Turner syndrome (TS) results from whole or partial monosomy X and is mediated by haploinsufficiency of genes that normally escape X-inactivation. Although a 45,X karyotype is observed in half of all TS cases, the most frequent variant TS karyotype includes the isodicentric X chromosome alone [46,X,idic(X)(p11)] or as a mosaic [46,X,idic(X)(p11)/45,X]. Given the mechanism of idic(X)(p11) rearrangement is poorly understood and breakpoint sequence information is unknown, this study sought to investigate the molecular mechanism of idic(X)(p11) formation by determining their precise breakpoint intervals. Karyotype analysis and fluorescence in situ hybridization mapping of eight idic(X)(p11) cell lines and three unbalanced Xp11.2 translocation lines identified the majority of breakpoints within a 5 Mb region, from approximately 53 to 58 Mb, in Xp11.1-p11.22, clustering into four regions. To further refine the breakpoints, a high-resolution oligonucleotide microarray (average of approximately 350 bp) was designed and array-based comparative genomic hybridization (aCGH) was performed on all 11 idic(X)(p11) and Xp11.2 translocation lines. aCGH analyses identified all breakpoint regions, including an idic(X)(p11) line with two potential breakpoints, one breakpoint shared between two idic(X)(p11) lines and two Xp translocations that shared breakpoints with idic(X)(p11) lines. Four of the breakpoint regions included large inverted repeats composed of repetitive gene clusters and segmental duplications, which corresponded to regions of copy-number variation. These data indicate that the rearrangement sites on Xp11.2 that lead to isodicentric chromosome formation and translocations are probably not random and suggest that the complex repetitive architecture of this region predisposes it to rearrangements, some of which are recurrent.

A paucity of heterochromatin at functional human neocentromeres.

BACKGROUND: Centromeres are responsible for the proper segregation of replicated chromatids during cell division. Neocentromeres are fully functional ectopic human centromeres that form on low-copy DNA sequences and permit analysis of centromere structure in relation to the underlying DNA sequence. Such structural analysis is not possible at endogenous centromeres because of the large amounts of repetitive alpha satellite DNA present. RESULTS: High-resolution chromatin immunoprecipitation (ChIP) on CHIP (microarray) analysis of three independent neocentromeres from chromosome 13q revealed that each neocentromere contained approximately 100 kb of centromere protein (CENP)-A in a two-domain organization. Additional CENP-A domains were observed in the vicinity of neocentromeres, coinciding with CpG islands at the 5' end of genes. Analysis of histone H3 dimethylated at lysine 4 (H3K4me2) revealed small domains at each neocentromere. However, these domains of H3K4me2 were also found in the equivalent non-neocentric chromosomes. A surprisingly minimal (approximately 15 kb) heterochromatin domain was observed at one of the neocentromeres, which formed in an unusual transposon-free region distal to the CENP-A domains. Another neocentromere showed a distinct absence of nearby significant domains of heterochromatin. A subtle defect in centromere cohesion detected at these neocentromeres may be due to the paucity of heterochromatin domains. CONCLUSIONS: This high-resolution mapping suggests that H3K4me2 does not seem sufficiently abundant to play a structural role at neocentromeres, as proposed for endogenous centromeres. Large domains of heterochromatin also do not appear necessary for centromere function. Thus, this study provides important insight into the structural requirements of human centromere function.

Analysis of the largest tandemly repeated DNA families in the human genome.

BACKGROUND: Tandemly Repeated DNA represents a large portion of the human genome, and accounts for a significant amount of copy number variation. Here we present a genome wide analysis of the largest tandem repeats found in the human genome sequence. RESULTS: Using Tandem Repeats Finder (TRF), tandem repeat arrays greater than 10 kb in total size were identified, and classified into simple sequence e.g. GAATG, classical satellites e.g. alpha satellite DNA, and locus specific VNTR arrays. Analysis of these large sequenced regions revealed that several "simple sequence" arrays actually showed complex domain and/or higher order repeat organization. Using additional methods, we further identified a total of 96 additional arrays with tandem repeat units greater than 2 kb (the detection limit of TRF), 53 of which contained genes or repeated exons. The overall size of an array of tandem 12 kb repeats which spanned a gap on chromosome 8 was found to be 600 kb to 1.7 Mbp in size, representing one of the largest non-centromeric arrays characterized. Several novel megasatellite tandem DNA families were observed that are characterized by repeating patterns of interspersed transposable elements that have expanded presumably by unequal crossing over. One of these families is found on 11 different chromosomes in >25 arrays, and represents one of the largest most widespread megasatellite DNA families. CONCLUSION: This study represents the most comprehensive genome wide analysis of large tandem repeats in the human genome, and will serve as an important resource towards understanding the organization and copy number variation of these complex DNA families.

Analysis of transposon interruptions suggests selection for L1 elements on the X chromosome.

It has been hypothesised that the massive accumulation of L1 transposable elements on the X chromosome is due to their function in X inactivation, and that the accumulation of Alu elements near genes is adaptive. We tested the possible selective advantage of these two transposable element (TE) families with a novel method, interruption analysis. In mammalian genomes, a large number of TEs interrupt other TEs due to the high overall abundance and age of repeats, and these interruptions can be used to test whether TEs are selectively neutral. Interruptions of TEs, which are beneficial for the host, are expected to be deleterious and underrepresented compared with neutral ones. We found that L1 elements in the regions of the X chromosome that contain the majority of the inactivated genes are significantly less frequently interrupted than on the autosomes, while L1s near genes that escape inactivation are interrupted with higher frequency, supporting the hypothesis that L1s on the X chromosome play a role in its inactivation. In addition, we show that TEs are less frequently interrupted in introns than in intergenic regions, probably due to selection against the expansion of introns, but the insertion pattern of Alus is comparable to other repeats.

The mouse X chromosome is enriched for multicopy testis genes showing postmeiotic expression.

According to the prevailing view, mammalian X chromosomes are enriched in spermatogenesis genes expressed before meiosis and deficient in spermatogenesis genes expressed after meiosis. The paucity of postmeiotic genes on the X chromosome has been interpreted as a consequence of meiotic sex chromosome inactivation (MSCI)--the complete silencing of genes on the XY bivalent at meiotic prophase. Recent studies have concluded that MSCI-initiated silencing persists beyond meiosis and that most genes on the X chromosome remain repressed in round spermatids. Here, we report that 33 multicopy gene families, representing approximately 273 mouse X-linked genes, are expressed in the testis and that this expression is predominantly in postmeiotic cells. RNA FISH and microarray analysis show that the maintenance of X chromosome postmeiotic repression is incomplete. Furthermore, X-linked multicopy genes exhibit a similar degree of expression as autosomal genes. Thus, not only is the mouse X chromosome enriched for spermatogenesis genes functioning before meiosis, but in addition, approximately 18% of mouse X-linked genes are expressed in postmeiotic cells.

Biased distributions and decay of long interspersed nuclear elements in the chicken genome.

The genomes of birds are much smaller than mammalian genomes, and transposable elements (TEs) make up only 10% of the chicken genome, compared with the 45% of the human genome. To study the mechanisms that constrain the copy numbers of TEs, and as a consequence the genome size of birds, we analyzed the distributions of LINEs (CR1's) and SINEs (MIRs) on the chicken autosomes and Z chromosome. We show that (1) CR1 repeats are longest on the Z chromosome and their length is negatively correlated with the local GC content; (2) the decay of CR1 elements is highly biased, and the 5'-ends of the insertions are lost much faster than their 3'-ends; (3) the GC distribution of CR1 repeats shows a bimodal pattern with repeats enriched in both AT-rich and GC-rich regions of the genome, but the CR1 families show large differences in their GC distribution; and (4) the few MIRs in the chicken are most abundant in regions with intermediate GC content. Our results indicate that the primary mechanism that removes repeats from the chicken genome is ectopic exchange and that the low abundance of repeats in avian genomes is likely to be the consequence of their high recombination rates.