Complex beta-satellite repeat structures and the expansion of the zinc finger gene cluster in 19p12

Rare Sequence Variation Underlying Suspected Familial Cerebral Small-Vessel Disease

BACKGROUND

Cerebral small-vessel disease (cSVD) is the leading monogenic cause of stroke. Despite genetic screening in routine diagnosis, many cases remain without a known causative variant. Using a cohort with suspected familial cSVD and whole-genome sequencing, we screened for variants in genes associated with monogenic cSVD and searched for novel variants associated with the disease.

METHODS AND RESULTS

Rare variants were identified in whole-genome sequencing data from the NBR (National Institute for Health Research BioResource Rare Disease) study. Pathogenic variants in known monogenic cSVD genes were identified. Gene-based burden tests and family analysis were performed to identify novel variants associated with familial cSVD. A total of 257 suspected cSVD cases (mean ± SD age, 56.2 ± 16.1 years), and 13 086 controls with other nonstroke diseases (5874 [44.9%] men) were studied. A total of 8.9% of the cases carried a variant in known cSVD genes. Excluding these known causes, 23.6% of unrelated subjects with cSVD carried predicted deleterious variants in the Genomics England gene panel, but no association was found with cSVD in burden tests. We identified potential associations with cSVD in noncoding genes, including RP4-568F9.3 (adjusted P = 7.1 × 10-25), RP3-466I7.1 (adjusted P = 8.9 × 10-16), and ZNF209P (adjusted P = 1.0 × 10-15), and matrisomal genes (adjusted P = 5.1 × 10-6), including FAM20C, INHA, LAMC1, and VWA5B2.

CONCLUSIONS

Predicted deleterious variants in known cSVD genes were present in 23.6% of unrelated cases with cSVD, but none of the genes were associated with the disease. Rare variants in noncoding and matrisomal genes could potentially contribute to cSVD development. These genes could play a role in tissue development and brain endothelial cell function. However, further studies are needed to confirm their pathophysiological roles.

Emergence of the ZNF675 Gene During Primate Evolution-Influenced Human Neurodevelopment Through Changing HES1 Autoregulation

In this study, we investigated recurrent copy number variations (CNVs) in the 19p12 locus, which are associated with neurodevelopmental disorders. The two genes in this locus, ZNF675 and ZNF681, arose via gene duplication in primates, and their presence in several pathological CNVs in the human population suggests that either or both of these genes are required for normal human brain development. ZNF675 and ZNF681 are members of the Krüppel-associated box zinc finger (KZNF) protein family, a class of transcriptional repressors important for epigenetic silencing of specific genomic regions. About 170 primate-specific KZNFs are present in the human genome. Although KZNFs are primarily associated with repressing retrotransposon-derived DNA, evidence is emerging that they can be co-opted for other gene regulatory processes. We show that genetic deletion of ZNF675 causes developmental defects in cortical organoids, and our data suggest that part of the observed neurodevelopmental phenotype is mediated by a gene regulatory role of ZNF675 on the promoter of the neurodevelopmental gene Hes family BHLH transcription factor 1 (HES1). We also find evidence for the recently evolved regulation of genes involved in neurological disorders, microcephalin 1 and sestrin 3. We show that ZNF675 interferes with HES1 auto-inhibition, a process essential for the maintenance of neural progenitors. As a striking example of how some KZNFs have integrated into preexisting gene expression networks, these findings suggest the emergence of ZNF675 has caused a change in the balance of HES1 autoregulation. The association of ZNF675 CNV with human developmental disorders and ZNF675-mediated regulation of neurodevelopmental genes suggests that it evolved into an important factor for human brain development.

Comparative Transcriptomics Data Profiling Reveals E2F Targets as an Important Biological Pathway Overexpressed in Intellectual Disability Disorder

Intellectual disability (ID) is an early childhood neurodevelopmental disorder that is characterized by impaired intellectual functioning and adaptive behavior. It is one of the major concerns in the field of neurodevelopmental disorders across the globe. Diversified approaches have been put forward to overcome this problem. Among all these approaches, high throughput transcriptomic analysis has taken an important dimension. The identification of genes causing ID rapidly increased over the past 3 to 5 years owing to the use of sophisticated high throughput sequencing platforms. Early monitoring and preventions are much important for such disorder as their progression occurs during fetal development. This study is an attempt to identify differentially expressed genes (DEGs) and upregulated biological processes involved in development of ID patients through comparative analysis of available transcriptomics data. A total of 7 transcriptomic studies were retrieved from National Center for Biotechnology Information (NCBI) and were subjected to quality check and trimming prior to alignment. The normalization and differential expression analysis were carried out using DESeq2 and EdgeR packages of Rstudio to identify DEGs in ID. In progression of the study, functional enrichment analysis of the results obtained from both DESeq2 and EdgeR was done using gene set enrichment analysis (GSEA) tool to identify major upregulated biological processes involved in ID. Our findings concluded that monitoring the level of E2F targets, estrogen, and genes related to oxidative phosphorylation, DNA repair, and glycolysis during the developmental stage of an individual can help in the early detection of ID disorder.

Genomic Tackling of Human Satellite DNA: Breaking Barriers through Time

(Peri)centromeric repetitive sequences and, more specifically, satellite DNA (satDNA) sequences, constitute a major human genomic component. SatDNA sequences can vary on a large number of features, including nucleotide composition, complexity, and abundance. Several satDNA families have been identified and characterized in the human genome through time, albeit at different speeds. Human satDNA families present a high degree of sub-variability, leading to the definition of various subfamilies with different organization and clustered localization. Evolution of satDNA analysis has enabled the progressive characterization of satDNA features. Despite recent advances in the sequencing of centromeric arrays, comprehensive genomic studies to assess their variability are still required to provide accurate and proportional representation of satDNA (peri)centromeric/acrocentric short arm sequences. Approaches combining multiple techniques have been successfully applied and seem to be the path to follow for generating integrated knowledge in the promising field of human satDNA biology.

The wide distribution and horizontal transfers of beta satellite DNA in eukaryotes

Beta satellite DNA (satDNA), also known as Sau3A sequences, are repetitive DNA sequences reported in human and primate genomes. It is previously thought that beta satDNAs originated in old world monkeys and bursted in great apes. In this study, we searched 7821 genome assemblies of 3767 eukaryotic species and found that beta satDNAs are widely distributed across eukaryotes. The four major branches of eukaryotes, animals, fungi, plants and Harosa/SAR, all have multiple clades containing beta satDNAs. These results were also confirmed by searching whole genome sequencing data (SRA) and PCR assay. Beta satDNA sequences were found in all the primate clades, as well as in Dermoptera and Scandentia, indicating that the beta satDNAs in primates might originate in the common ancestor of Primatomorpha or Euarchonta. In contrast, the widely patchy distribution of beta satDNAs across eukaryotes presents a typical scenario of multiple horizontal transfers.

Genome-Wide DNA Alterations in X-Irradiated Human Gingiva Fibroblasts

While ionizing radiation (IR) is a powerful tool in medical diagnostics, nuclear medicine, and radiology, it also is a serious threat to the integrity of genetic material. Mutagenic effects of IR to the human genome have long been the subject of research, yet still comparatively little is known about the genome-wide effects of IR exposure on the DNA-sequence level. In this study, we employed high throughput sequencing technologies to investigate IR-induced DNA alterations in human gingiva fibroblasts (HGF) that were acutely exposed to 0.5, 2, and 10 Gy of 240 kV X-radiation followed by repair times of 16 h or 7 days before whole-genome sequencing (WGS). Our analysis of the obtained WGS datasets revealed patterns of IR-induced variant (SNV and InDel) accumulation across the genome, within chromosomes as well as around the borders of topologically associating domains (TADs). Chromosome 19 consistently accumulated the highest SNVs and InDels events. Translocations showed variable patterns but with recurrent chromosomes of origin (e.g., Chr7 and Chr16). IR-induced InDels showed a relative increase in number relative to SNVs and a characteristic signature with respect to the frequency of triplet deletions in areas without repetitive or microhomology features. Overall experimental conditions and datasets the majority of SNVs per genome had no or little predicted functional impact with a maximum of 62, showing damaging potential. A dose-dependent effect of IR was surprisingly not apparent. We also observed a significant reduction in transition/transversion (Ti/Tv) ratios for IR-dependent SNVs, which could point to a contribution of the mismatch repair (MMR) system that strongly favors the repair of transitions over transversions, to the IR-induced DNA-damage response in human cells. Taken together, our results show the presence of distinguishable characteristic patterns of IR-induced DNA-alterations on a genome-wide level and implicate DNA-repair mechanisms in the formation of these signatures.

GDF15: A Hormone Conveying Somatic Distress to the Brain

GDF15 has recently gained scientific and translational prominence with the discovery that its receptor is a GFRAL-RET heterodimer of which GFRAL is expressed solely in the hindbrain. Activation of this receptor results in reduced food intake and loss of body weight and is perceived and recalled by animals as aversive. This information encourages a revised interpretation of the large body of previous research on the protein. GDF15 can be secreted by a wide variety of cell types in response to a broad range of stressors. We propose that central sensing of GDF15 via GFRAL-RET activation results in behaviors that facilitate the reduction of exposure to a noxious stimulus. The human trophoblast appears to have hijacked this signal, producing large amounts of GDF15 from early pregnancy. We speculate that this encourages avoidance of potential teratogens in pregnancy. Circulating GDF15 levels are elevated in a range of human disease states, including various forms of cachexia, and GDF15-GFRAL antagonism is emerging as a therapeutic strategy for anorexia/cachexia syndromes. Metformin elevates circulating GDF15 chronically in humans and the weight loss caused by this drug appears to be dependent on the rise in GDF15. This supports the concept that chronic activation of the GDF15-GFRAL axis has efficacy as an antiobesity agent. In this review, we examine the science of GDF15 since its identification in 1997 with our interpretation of this body of work now being assisted by a clear understanding of its highly selective central site of action.

Instability of Short Arm of Acrocentric Chromosomes: Lesson from Non-Acrocentric Satellited Chromosomes. Report of 24 Unrelated Cases

Satellited non-acrocentric autosomal chromosomes (ps–qs-chromosomes) are the result of an interchange between sub- or telomeric regions of autosomes and the p arm of acrocentrics. The sequence homology at the rearrangement breakpoints appears to be, among others, the most frequent mechanism generating these variant chromosomes. The unbalanced carriers of this type of translocation may or may not display phenotypic abnormalities. With the aim to understand the causative mechanism, we revised all the ps–qs-chromosomes identified in five medical genetics laboratories, which used the same procedures for karyotype analysis, reporting 24 unrelated cases involving eight chromosomes. In conclusion, we observed three different scenarios: true translocation, benign variant and complex rearrangement. The detection of translocation partners is essential to evaluate possible euchromatic unbalances and to infer their effect on phenotype. Moreover, we emphasize the importance to perform both, molecular and conventional cytogenetics methods, to better understand the behavior of our genome.

On the relations of phase separation and Hi-C maps to epigenetics

The relationship between compartmentalization of the genome and epigenetics is long and hoary. In 1928, Heitz defined heterochromatin as the largest differentiated chromatin compartment in eukaryotic nuclei. Müller's discovery of position-effect variegation in 1930 went on to show that heterochromatin is a cytologically visible state of heritable (epigenetic) gene repression. Current insights into compartmentalization have come from a high-throughput top-down approach where contact frequency (Hi-C) maps revealed the presence of compartmental domains that segregate the genome into heterochromatin and euchromatin. It has been argued that the compartmentalization seen in Hi-C maps is owing to the physiochemical process of phase separation. Oddly, the insights provided by these experimental and conceptual advances have remained largely silent on how Hi-C maps and phase separation relate to epigenetics. Addressing this issue directly in mammals, we have made use of a bottom-up approach starting with the hallmarks of constitutive heterochromatin, heterochromatin protein 1 (HP1) and its binding partner the H3K9me2/3 determinant of the histone code. They are key epigenetic regulators in eukaryotes. Both hallmarks are also found outside mammalian constitutive heterochromatin as constituents of larger (0.1-5 Mb) heterochromatin-like domains and smaller (less than 100 kb) complexes. The well-documented ability of HP1 proteins to function as bridges between H3K9me2/3-marked nucleosomes contributes to polymer-polymer phase separation that packages epigenetically heritable chromatin states during interphase. Contacts mediated by HP1 'bridging' are likely to have been detected in Hi-C maps, as evidenced by the B4 heterochromatic subcompartment that emerges from contacts between large KRAB-ZNF heterochromatin-like domains. Further, mutational analyses have revealed a finer, innate, compartmentalization in Hi-C experiments that probably reflect contacts involving smaller domains/complexes. Proteins that bridge (modified) DNA and histones in nucleosomal fibres-where the HP1-H3K9me2/3 interaction represents the most evolutionarily conserved paradigm-could drive and generate the fundamental compartmentalization of the interphase nucleus. This has implications for the mechanism(s) that maintains cellular identity, be it a terminally differentiated fibroblast or a pluripotent embryonic stem cell.

Telomere length, ZNF208 genetic variants and risk of chronic obstructive pulmonary disease in the Hainan Li population

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a disease characterized by airflow limitation. It is not completely reversible and shows progressive development. ZNF208 rs8105767 affects telomere length, although the impact of telomere on COPD is still controversial. In the present study, we aimed to explore the impact of the ZNF208 gene polymorphism on telomere length and also that of telomere length on COPD in the Hainan Li population.

METHODS

In total, 270 COPD patients and 288 controls were recruited. Telomere length was measured by a quantitative real-time polymerase chain reaction. Five single nucleotide polymorphisms in ZNF208 were selected and genotyping was performed using MassARRAY software (Agena Bioscience Co. Ltd, San Diego, CA, USA). Differences in telomere length among the subjects with three genotypes of related genes were assessed using analysis of variance. Unconditional logistic regression was used to calculate odds ratios (OR) as the indicator of association between telomere length and COPD risk.

RESULTS

Relative telomere length in the COPD group and control group was 0.66 ± 0.47 and 1.44 ± 0.89, respectively. We grouped according to a median of 0.8284 for telomere length and observed that the risk of COPD for individuals with a telomere length less than 0.8284 is 2.92 times that for individuals with a telomere length longer than 0.8284 (OR = 2.92, 95% confidence interval = 2.01-4.25, p = 1.91 × 10^-8 ). Subjects carrying the G allele of rs2188972 had a longer telomere length. Subjects carrying the carrying the CA genotype of rs8103163 and AC genotype of rs7248488 had a longer telomere length compared to wild-type individuals.

CONCLUSIONS

Shorter telomeres increase COPD risk and the ZNF208 polymorphism affects telomere length in COPD patients.

Single nucleotide polymorphisms in ZNF208 are associated with increased risk for HBV in Chinese people

Single nucleotide polymorphisms (SNPs) in ZNF208 may be associated with susceptibility to Hepatitis B virus (HBV). In the current study, we analyzed the association between ZNF208 SNPs and risk of HBV in 242 HBV patients and 300 healthy subjects from the Xi'an area of Chinese Han Population. Of the five SNPs examined, rs2188971 (OR: 1.36, 95% CI: 1.04-1.76, P = 0.022), rs8103163 (OR: 1.40, 95% CI: 1.08-1.82, P = 0.010) and rs7248488 (OR: 1.38, 95% CI: 1.07-1.79, P = 0.014) were correlated with HBV susceptibility based on Chi-square tests. After the P-values were adjusted by Bonferroni correction, there only rs8103163 (P = 0.050) was slightly with increased HBV risk. Additionally, haplotype A_rs2188972T_rs2188971A_rs8103163A_rs7248488 (OR = 1.42; 95% C I, 1.10-1.85; P = 0.008) was found to increase susceptibility of suffering from HBV. These findings suggest that ZNF208 polymorphisms may contribute to the development of HBV.

Association between TNIP1, MPHOSPH6 and ZNF208 genetic polymorphisms and the coronary artery disease risk in Chinese Han population

Introduction

Coronary artery disease (CAD) is a common disease and among the leading cause of death in the general population. Inherited factors are involved in the pathogenesis of CAD.

Aims

Our study examined whether SNPs in TNIP1, MPHOSPH6, ZNF208 to be associated with CAD risk in a Chinese Han population. We recruited 596 CAD patients, 603 controls and genotyping fifteen SNPs using Sequenom MassARRAY. For association analysis between TNIP1, MPHOSPH6 and ZNF208 and CAD was determined by Odds ratios (ORs) with 95% confidence intervals (CIs) using Logistic Regression.

Results

The results indicated in allel model, the rs960709 in TNIP1 was associated with CAD risk (OR = 0.78, 95%CI = 0.65-0.94, P=0.010). The genetic model results showed that the rs960709 (A/G) polymorphism was associated with the risk of developing CAD in codominant, Dominant and Log-additive. The rs1056654 A/A allele and CAD patients compared to the healthy controls in recessive model (OR = 0.55, 95%CI = 0.34-0.90; P = 0.018). We also found that three SNPS in ZNF208 associated with CAD, respectively, rs2188971, rs8103163 and rs7248488.

Linkage disequilibrium (LD) and haplotype analyses of the SNPs found that the CTA haplotype (rs1056675, rs1056654, rs11859599) and rs2188972A/rs2188971T/rs8103163A/rs7248488A (ATAA) were associated with CAD.

Conclusion

In conclusion, the present study provided evidence that SNPs in the TNIP1, ZNF208 and MPHOSPH6 were associated with CAD in Chinese Han population. It is possible that these SNPs are CAD risk factors and these data can provide.

Meta-analysis of five genome-wide association studies identifies multiple new loci associated with testicular germ cell tumor

The international Testicular Cancer Consortium (TECAC) combined five published genome-wide association studies of testicular germ cell tumor (TGCT; 3,558 cases and 13,970 controls) to identify new susceptibility loci. We conducted a fixed-effects meta-analysis, including, to our knowledge, the first analysis of the X chromosome. Eight new loci mapping to 2q14.2, 3q26.2, 4q35.2, 7q36.3, 10q26.13, 15q21.3, 15q22.31, and Xq28 achieved genome-wide significance (P < 5 × 10-8). Most loci harbor biologically plausible candidate genes. We refined previously reported associations at 9p24.3 and 19p12 by identifying one and three additional independent SNPs, respectively. In aggregate, the 39 independent markers identified to date explain 37% of father-to-son familial risk, 8% of which can be attributed to the 12 new signals reported here. Our findings substantially increase the number of known TGCT susceptibility alleles, move the field closer to a comprehensive understanding of the underlying genetic architecture of TGCT, and provide further clues to the etiology of TGCT.

ZNF208 polymorphisms associated with ischemic stroke in a southern Chinese Han population

BACKGROUND

Ischemic stroke is one of the most common diseases with a high burden of neurological deficits, disability and death. Zinc finger protein 208 (ZNF208) was found to be involved in coronary heart disease, although little information is available about its association with ischemic stroke. We performed the present case-control study to clarify the association between single-nucleotide polymorphisms (SNPs) within ZNF208 and the risk of ischemic stroke in a southern Chinese Han population.

METHODS

A total of 799 subjects (400 cases and 399 healthy controls) were enrolled in the present study. Five SNPs within ZNF208 gene were selected and genotyped using Sequenom MassARRY technology (Sequenom, Inc., San Diego, CA, USA). Data management and statistical analyses were conducted using Sequenom Typer, version 4.0, and a chi-squared test, as well as unconditional logistic regression.

RESULTS

Statistical results showed that three variants were associated with the risk of ischemic stroke under allele models (rs2188971, rs2188972, rs8103163 and rs7248488). The variant rs2188972 was also associated with the risk of ischemic stroke in a recessive model after adjustment for age and sex. Haplotype analysis suggested that a significant difference existed between the A_rs2188972 T_rs2188971 A_rs8103163 A_rs7248488 haplotype and the risk of ischemic stroke, although this disappeared after adjustment for sex and age.

CONCLUSIONS

The results obtained in the present study indicate a potential association between ZNF208 variants and the risk of ischemic risk in a southern Chinese Han population.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

Association between genetic risk score for telomere length and risk of breast cancer

PURPOSE

While leukocyte telomere length (TL) has been associated with breast cancer risk, limited information is available regarding the role of genetically-determined TL on breast cancer risk. We investigated whether aggregated TL-associated variants are associated with the risk of breast cancer in 2,865 breast cancer cases and 2,285 controls from the Shanghai Breast Cancer Genetics Study.

METHODS

Six genetic variants, identified through a genome-wide association study (GWAS) of TL in European-ancestry participants, were included in the study. A separate sample [n = 1,536, from the Shanghai Women's Health Study (SWHS), for whom information on both phenotypical leukocyte TL and genetic information was collected] was used to evaluate the association of six variants with TL in Asians. Three genetic risk scores (GRSs), based on the number of alleles associated with shorter TL that each individual carries for the six variants, were derived for the study: un-weighted, internally weighted (from the SWHS), and externally weighted (from the European-ancestry GWAS study), and evaluated for their association with breast cancer risk by applying logistic regression analysis.

RESULTS

Both internally and externally weighted GRSs were significantly associated with a decreased risk of breast cancer (OR 0.83, 95 % CI 0.72-0.95 and OR 0.84, 95 % CI 0.74-0.96, respectively, for tertile 3 vs. tertile 1). Non-genetic risk factors for breast cancer (i.e., age, years of menstruation/reproduction, oral contraceptive usage, and BMI) did not modify the association between GRSs and the risk of breast cancer.

CONCLUSION

Our results suggest that short TL, determined by genetic factors, may be associated with a reduced susceptibility to breast cancer.

The octopus genome and the evolution of cephalopod neural and morphological novelties

Coleoid cephalopods (octopus, squid, and cuttlefish) are active, resourceful predators with a rich behavioral repertoire¹. They have the largest nervous systems among the invertebrates² and present other striking morphological innovations including camera-like eyes, prehensile arms, a highly derived early embryogenesis, and the most sophisticated adaptive coloration system among all animals^1,3. To investigate the molecular bases of cephalopod brain and body innovations we sequenced the genome and multiple transcriptomes of the California two-spot octopus, Octopus bimaculoides. We found no evidence for hypothesized whole genome duplications in the octopus lineage^4–6. The core developmental and neuronal gene repertoire of the octopus is broadly similar to that found across invertebrate bilaterians, except for massive expansions in two gene families formerly thought to be uniquely enlarged in vertebrates: the protocadherins, which regulate neuronal development, and the C2H2 superfamily of zinc finger transcription factors. Extensive mRNA editing generates transcript and protein diversity in genes involved in neural excitability, as previously described⁷, as well as in genes participating in a broad range of other cellular functions. We identified hundreds of cephalopod-specific genes, many of which showed elevated expression levels in such specialized structures as the skin, the suckers, and the nervous system. Finally, we found evidence for large-scale genomic rearrangements that are closely associated with transposable element expansions. Our analysis suggests that substantial expansion of a handful of gene families, along with extensive remodeling of genome linkage and repetitive content, played a critical role in the evolution of cephalopod morphological innovations, including their large and complex nervous systems.

G9a/GLP-dependent histone H3K9me2 patterning during human hematopoietic stem cell lineage commitment

G9a and GLP are conserved protein methyltransferases that play key roles during mammalian development through mono- and dimethylation of histone H3 Lys 9 (H3K9me1/2), modifications associated with transcriptional repression. During embryogenesis, large H3K9me2 chromatin territories arise that have been proposed to reinforce lineage choice by affecting high-order chromatin structure. Here we report that in adult human hematopoietic stem and progenitor cells (HSPCs), H3K9me2 chromatin territories are absent in primitive cells and are formed de novo during lineage commitment. In committed HSPCs, G9a/GLP activity nucleates H3K9me2 marks at CpG islands and other genomic sites within genic regions, which then spread across most genic regions during differentiation. Immunofluorescence assays revealed the emergence of H3K9me2 nuclear speckles in committed HSPCs, consistent with progressive marking. Moreover, gene expression analysis indicated that G9a/GLP activity suppresses promiscuous transcription of lineage-affiliated genes and certain gene clusters, suggestive of regulation of HSPC chromatin structure. Remarkably, HSPCs continuously treated with UNC0638, a G9a/GLP small molecular inhibitor, better retain stem cell-like phenotypes and function during in vitro expansion. These results suggest that G9a/GLP activity promotes progressive H3K9me2 patterning during HSPC lineage specification and that its inhibition delays HSPC lineage commitment. They also inform clinical manipulation of donor-derived HSPCs.

Expression, tandem repeat copy number variation and stability of four macrosatellite arrays in the human genome

Background

Macrosatellites are some of the largest variable number tandem repeats in the human genome, but what role these unusual sequences perform is unknown. Their importance to human health is clearly demonstrated by the 4q35 macrosatellite D4Z4 that is associated with the onset of the muscle degenerative disease facioscapulohumeral muscular dystrophy. Nevertheless, many other macrosatellite arrays in the human genome remain poorly characterized.

Results

Here we describe the organization, tandem repeat copy number variation, transmission stability and expression of four macrosatellite arrays in the human genome: the TAF11-Like array located on chromosomes 5p15.1, the SST1 arrays on 4q28.3 and 19q13.12, the PRR20 array located on chromosome 13q21.1, and the ZAV array at 9q32. All are polymorphic macrosatellite arrays that at least for TAF11-Like and SST1 show evidence of meiotic instability. With the exception of the SST1 array that is ubiquitously expressed, all are expressed at high levels in the testis and to a lesser extent in the brain.

Conclusions

Our results extend the number of characterized macrosatellite arrays in the human genome and provide the foundation for formulation of hypotheses to begin assessing their functional role in the human genome.

Complex evolution of a Y-chromosomal double homeobox 4 (DUX4)-related gene family in hominoids

The human Y chromosome carries four human Y-chromosomal euchromatin/heterochromatin transition regions, all of which are characterized by the presence of interchromosomal segmental duplications. The Yq11.1/Yq11.21 transition region harbours a peculiar segment composed of an imperfectly organized tandem-repeat structure encoding four members of the double homeobox (DUX) gene family. By comparative fluorescence in situ hybridization (FISH) analysis we have documented the primary appearance of Y-chromosomal DUX genes (DUXY) on the gibbon Y chromosome. The major amplification and dispersal of DUXY paralogs occurred after the gibbon and hominid lineages had diverged. Orthologous DUXY loci of human and chimpanzee show a highly similar structural organization. Sequence alignment survey, phylogenetic reconstruction and recombination detection analyses of human and chimpanzee DUXY genes revealed the existence of all copies in a common ancestor. Comparative analysis of the circumjacent beta-satellites indicated that DUXY genes and beta-satellites evolved in concert. However, evolutionary forces acting on DUXY genes may have induced amino acid sequence differences in the orthologous chimpanzee and human DUXY open reading frames (ORFs). The acquisition of complete ORFs in human copies might relate to evolutionary advantageous functions indicating neo-functionalization. We propose an evolutionary scenario in which an ancestral tandem array DUX gene cassette transposed to the hominoid Y chromosome followed by lineage-specific chromosomal rearrangements paved the way for a species-specific evolution of the Y-chromosomal members of a large highly diverged homeobox gene family.

Modulation of gene expression in U251 glioblastoma cells by binding of mutant p53 R273H to intronic and intergenic sequences

Missense point mutations in the TP53 gene are frequent genetic alterations in human tumor tissue and cell lines derived thereof. Mutant p53 (mutp53) proteins have lost sequence-specific DNA binding, but have retained the ability to interact in a structure-selective manner with non-B DNA and to act as regulators of transcription. To identify functional binding sites of mutp53, we established a small library of genomic sequences bound by p53^R273H in U251 human glioblastoma cells using chromatin immunoprecipitation (ChIP). Mutp53 binding to isolated DNA fragments confirmed the specificity of the ChIP. The mutp53 bound DNA sequences are rich in repetitive DNA elements, which are dispersed over non-coding DNA regions. Stable down-regulation of mutp53 expression strongly suggested that mutp53 binding to genomic DNA is functional. We identified the PPARGC1A and FRMD5 genes as p53^R273H targets regulated by binding to intronic and intra-genic sequences. We propose a model that attributes the oncogenic functions of mutp53 to its ability to interact with intronic and intergenic non-B DNA sequences and modulate gene transcription via re-organization of chromatin.

Origin of primate orphan genes: a comparative genomics approach

Genomes contain a large number of genes that do not have recognizable homologues in other species and that are likely to be involved in important species-specific adaptive processes. The origin of many such "orphan" genes remains unknown. Here we present the first systematic study of the characteristics and mechanisms of formation of primate-specific orphan genes. We determine that codon usage values for most orphan genes fall within the bulk of the codon usage distribution of bona fide human proteins, supporting their current protein-coding annotation. We also show that primate orphan genes display distinctive features in relation to genes of wider phylogenetic distribution: higher tissue specificity, more rapid evolution, and shorter peptide size. We estimate that around 24% are highly divergent members of mammalian protein families. Interestingly, around 53% of the orphan genes contain sequences derived from transposable elements (TEs) and are mostly located in primate-specific genomic regions. This indicates frequent recruitment of TEs as part of novel genes. Finally, we also obtain evidence that a small fraction of primate orphan genes, around 5.5%, might have originated de novo from mammalian noncoding genomic regions.

Evolution of C2H2-zinc finger genes and subfamilies in mammals: species-specific duplication and loss of clusters, genes and effector domains

Background

C2H2 zinc finger genes (C2H2-ZNF) constitute the largest class of transcription factors in humans and one of the largest gene families in mammals. Often arranged in clusters in the genome, these genes are thought to have undergone a massive expansion in vertebrates, primarily by tandem duplication. However, this view is based on limited datasets restricted to a single chromosome or a specific subset of genes belonging to the large KRAB domain-containing C2H2-ZNF subfamily.

Results

Here, we present the first comprehensive study of the evolution of the C2H2-ZNF family in mammals. We assembled the complete repertoire of human C2H2-ZNF genes (718 in total), about 70% of which are organized into 81 clusters across all chromosomes. Based on an analysis of their N-terminal effector domains, we identified two new C2H2-ZNF subfamilies encoding genes with a SET or a HOMEO domain. We searched for the syntenic counterparts of the human clusters in other mammals for which complete gene data are available: chimpanzee, mouse, rat and dog. Cross-species comparisons show a large variation in the numbers of C2H2-ZNF genes within homologous mammalian clusters, suggesting differential patterns of evolution. Phylogenetic analysis of selected clusters reveals that the disparity in C2H2-ZNF gene repertoires across mammals not only originates from differential gene duplication but also from gene loss. Further, we discovered variations among orthologs in the number of zinc finger motifs and association of the effector domains, the latter often undergoing sequence degeneration. Combined with phylogenetic studies, physical maps and an analysis of the exon-intron organization of genes from the SCAN and KRAB domains-containing subfamilies, this result suggests that the SCAN subfamily emerged first, followed by the SCAN-KRAB and finally by the KRAB subfamily.

Conclusion

Our results are in agreement with the "birth and death hypothesis" for the evolution of C2H2-ZNF genes, but also show that this hypothesis alone cannot explain the considerable evolutionary variation within the subfamilies of these genes in mammals. We, therefore, propose a new model involving the interdependent evolution of C2H2-ZNF gene subfamilies.

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

Effects of recombination rate on human endogenous retrovirus fixation and persistence

Endogenous retroviruses (ERVs) result from germ line infections by exogenous retroviruses. They can proliferate within the genome of their host species until they are either inactivated by mutation or removed by recombinational deletion. ERVs belong to a diverse group of mobile genetic elements collectively termed transposable elements (TEs). Numerous studies have attempted to elucidate the factors determining the genomic distribution and persistence of TEs. Here we show that, within humans, gene density and not recombination rate correlates with fixation of endogenous retroviruses, whereas the local recombination rate determines their persistence in a full-length state. Recombination does not appear to influence fixation either via the ectopic exchange model or by indirect models based on the efficacy of selection. We propose a model linking rates of meiotic recombination to the probability of recombinational deletion to explain the effect of recombination rate on persistence. Chromosomes 19 and Y are exceptions, possessing more elements than other regions, and we suggest this is due to low gene density and elevated rates of human ERV integration in males for chromosome Y and segmental duplication for chromosome 19.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5'-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Human heterochromatin proteins form large domains containing KRAB-ZNF genes

Heterochromatin is important for gene regulation and chromosome structure, but the genes that are occupied by heterochromatin proteins in the mammalian genome are largely unknown. We have adapted the DamID method to systematically identify target genes of the heterochromatin proteins HP1 and SUV39H1 in human and mouse cells. Unexpectedly, we found that CBX1 (formerly HP1beta) and SUV39H1 bind to genes encoding KRAB domain containing zinc finger (KRAB-ZNF) transcriptional repressors. These genes constitute one of the largest gene families and are organized in clusters in the human genome. Preference of CBX1 for this gene family was observed in both human and mouse cells. High-resolution mapping on human chromosome 19 revealed that CBX1 coats large domains 0.1-4 Mb in size, which coincide with the position of KRAB-ZNF gene clusters. These domains show an intricate CBX1 binding pattern: While CBX1 is globally elevated throughout the domains, it is absent from the promoters and binds more strongly to the 3' ends of KRAB-ZNF genes. KRAB-ZNF domains contain large numbers of LINE elements, which may contribute to CBX1 recruitment. These results uncover a surprising link between heterochromatin and a large family of regulatory genes in mammals. We suggest a role for heterochromatin in the evolution of the KRAB-ZNF gene family.

Chromosome-specific DNA repeat probes

In research as well as in clinical applications, fluorescence in situ hybridization (FISH) has gained increasing popularity as a highly sensitive technique to study cytogenetic changes. Today, hundreds of commercially available DNA probes serve the basic needs of the biomedical research community. Widespread applications, however, are often limited by the lack of appropriately labeled, specific nucleic acid probes. We describe two approaches for an expeditious preparation of chromosome-specific DNAs and the subsequent probe labeling with reporter molecules of choice. The described techniques allow the preparation of highly specific DNA repeat probes suitable for enumeration of chromosomes in interphase cell nuclei or tissue sections. In addition, there is no need for chromosome enrichment by flow cytometry and sorting or molecular cloning. Our PCR-based method uses either bacterial artificial chromosomes or human genomic DNA as templates with alpha-satellite-specific primers. Here we demonstrate the production of fluorochrome-labeled DNA repeat probes specific for human chromosomes 17 and 18 in just a few days without the need for highly specialized equipment and without the limitation to only a few fluorochrome labels.

A novel Krüppel related factor consisting of only a KRAB domain is expressed in the murine trigeminal ganglion

The largest family of zinc-finger (ZnF) transcription factors is that containing the Krüppel-associated box, or KRAB domain. The amino-terminal KRAB domain of these proteins functions as a transcriptional repressor with the downstream ZnF motifs providing DNA-binding specificity. Here we report the identification and characterization of a novel murine Krüppel-related factor (KLF), MIF1, which contains a KRAB domain but lacks a ZnF motif. Western blot analysis identified MIF1-like proteins in the murine trigeminal ganglion (TG) and immunostaining localized these proteins primarily to the cytoplasm of TG neuronal cell bodies. In situ hybridization for Mif1 transcripts confirms the selective expression of Mif1 in TG neurons. Consistent with the non-nuclear localization of MIF1 we could detect no transcriptional repressor activity of the MIF1 protein. However MIF1 appears to be capable of interacting with the co-repressor TIF1beta and exhibits transcription repressor activity when fused to yeast GAL4 binding domain protein. Genomic analysis of Mif1 sequence suggests that the Mif1 transcript may result from splicing of a longer KRAB-ZnF containing transcript.

Evolutionary expansion and divergence in the ZNF91 subfamily of primate-specific zinc finger genes

Most genes are conserved in mammals, but certain gene families have acquired large numbers of lineage-specific loci through repeated rounds of gene duplication, divergence, and loss that have continued in each mammalian group. One such family encodes KRAB-zinc finger (KRAB-ZNF) proteins, which function as transcriptional repressors. One particular subfamily of KRAB-ZNF genes, including ZNF91, has expanded specifically in primates to comprise more than 110 loci in the human genome. Genes of the ZNF91 subfamily reside in large gene clusters near centromeric regions of human chromosomes 19 and 7 with smaller clusters or isolated copies in other locations. Phylogenetic analysis indicates that many of these genes arose before the split between the New and Old World monkeys, but the ZNF91 subfamily has continued to expand and diversify throughout the evolution of apes and humans. Paralogous loci are distinguished by divergence within their zinc finger arrays, indicating selection for proteins with different regulatory targets. In addition, many loci produce multiple alternatively spliced transcripts encoding proteins that may serve separate and perhaps even opposing regulatory roles because of the modular motif structure of KRAB-ZNF genes. The tissue-specific expression patterns and rapid structural divergence of ZNF91 subfamily genes suggest a role in determining gene expression differences between species and the evolution of novel primate traits.

A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors

Krüppel-type zinc finger (ZNF) motifs are prevalent components of transcription factor proteins in all eukaryotes. KRAB-ZNF proteins, in which a potent repressor domain is attached to a tandem array of DNA-binding zinc-finger motifs, are specific to tetrapod vertebrates and represent the largest class of ZNF proteins in mammals. To define the full repertoire of human KRAB-ZNF proteins, we searched the genome sequence for key motifs and then constructed and manually curated gene models incorporating those sequences. The resulting gene catalog contains 423 KRAB-ZNF protein-coding loci, yielding alternative transcripts that altogether predict at least 742 structurally distinct proteins. Active rounds of segmental duplication, involving single genes or larger regions and including both tandem and distributed duplication events, have driven the expansion of this mammalian gene family. Comparisons between the human genes and ZNF loci mined from the draft mouse, dog, and chimpanzee genomes not only identified 103 KRAB-ZNF genes that are conserved in mammals but also highlighted a substantial level of lineage-specific change; at least 136 KRAB-ZNF coding genes are primate specific, including many recent duplicates. KRAB-ZNF genes are widely expressed and clustered genes are typically not coregulated, indicating that paralogs have evolved to fill roles in many different biological processes. To facilitate further study, we have developed a Web-based public resource with access to gene models, sequences, and other data, including visualization tools to provide genomic context and interaction with other public data sets.

Identification of a novel zinc finger protein gene (ZNF298) in the GAP2 of human chromosome 21q

We have isolated a novel zinc finger protein gene, designated ZNF298, as a candidate gene for a particular phenotype of Down syndrome or bipolar affective disorder (BPAD) which maps to human chromosome 21q22.3. ZNF298 gene consists of 25 exons spanning approximately 80kb in a direction from the telomere to centromere. There are four kinds of transcripts that harbor three types of 3' UTR. These four transcripts (ZNF298a, ZNF298b, ZNF298c, and ZNF298d) contain putative open reading frames encoding 1178, 1198, 555, and 515 amino acids, respectively. ZNF298 gene was ubiquitously expressed in various tissues at very low level. The protein motif analysis revealed that ZNF298 proteins contain a SET [Su(var)3-9, Enhancer-of-zeste, Trithorax] domain, multiple C2H2-type zinc finger (ZnF_C2H2) domains, several nuclear localization signals (NLSs), and PEST sequences. Nuclear localization of ZNF298 protein was confirmed by transfection of expression vector of GFP-tagged protein into two human cell lines. Interestingly, this gene crosses over a clone gap (GAP2) remaining in the band 21q22.3. We obtained the DNA fragments corresponding to GAP2 using ZNF298 cDNA sequence as anchor primers for PCR and determined its genomic DNA sequence.

Clustering of genes coding for DNA binding proteins in a region of atypical evolution of the human genome

Comparison of the human and mouse genomes has revealed that significant variations in evolutionary rates exist among genomic regions and that a large part of this variation is interchromosomal. We confirm in this work, using a large collection of introns, that human chromosome 19 is the one that shows the highest divergence with respect to mouse. To search for other differences among chromosomes, we examine the distribution of gene functions in human and mouse chromosomes using the Gene Ontology definitions. We found by correspondence analysis that among the strongest clusterings of gene functions in human chromosomes is a group of genes coding for DNA binding proteins in chromosome 19. Interestingly, chromosome 19 also has a very high GC content, a feature that has been proposed to promote an opening of the chromatin, thereby facilitating binding of proteins to the DNA helix. In the mouse genome, however, a similar aggregation of genes coding for DNA binding proteins and high GC content cannot be found. This suggests that the distribution of genes coding for DNA binding proteins and the variations of the chromatin accessibility to these proteins are different in the human and mouse genomes. It is likely that the overall high synonymous and intron rates in chromosome 19 are a by-product of the high GC content of this chromosome.

An ORFeome-based analysis of human transcription factor genes and the construction of a microarray to interrogate their expression

Transcription factors (TFs) are essential regulators of gene expression, and mutated TF genes have been shown to cause numerous human genetic diseases. Yet to date, no single, comprehensive database of human TFs exists. In this work, we describe the collection of an essentially complete set of TF genes from one depiction of the human ORFeome, and the design of a microarray to interrogate their expression. Taking 1468 known TFs from TRANSFAC, InterPro, and FlyBase, we used this seed set to search the ScriptSure human transcriptome database for additional genes. ScriptSure's genome-anchored transcript clusters allowed us to work with a nonredundant high-quality representation of the human transcriptome. We used a high-stringency similarity search by using BLASTN, and a protein motif search of the human ORFeome by using hidden Markov models of DNA-binding domains known to occur exclusively or primarily in TFs. Four hundred ninety-four additional TF genes were identified in the overlap between the two searches, bringing our estimate of the total number of human TFs to 1962. Zinc finger genes are by far the most abundant family (762 members), followed by homeobox (199 members) and basic helix-loop-helix genes (117 members). We designed a microarray of 50-mer oligonucleotide probes targeted to a unique region of the coding sequence of each gene. We have successfully used this microarray to interrogate TF gene expression in species as diverse as chickens and mice, as well as in humans.

Sry associates with the heterochromatin protein 1 complex by interacting with a KRAB domain protein

In mammals, the SRY/Sry gene on the Y chromosome is necessary and sufficient for a bipotential gonad to develop into a testis, regardless of its chromosomal sex. The SRY/Sry gene encodes a protein that belongs to a high-mobility-group (HMG) box protein family and that has been postulated to modulate the expression of genes necessary for male gonadal differentiation. Using a yeast two-hybrid screen, we identified a novel protein containing only a Krüppel-associated box (KRAB) domain, which is hereafter named KRAB-O (KRAB Only), as an SRY-interacting protein. The KRAB-O protein is encoded by an alternatively spliced transcript from the Zfp208 locus that also produces another transcript coding for a KRAB-zinc finger protein, ZFP208. The interaction of the mouse SRY with KRAB-O was further confirmed by glutathione S-transferase pull-down assay and coimmunoprecipitation in transfected COS7 cells. The KRAB-O interaction domain in both the human and mouse SRY was mapped to the bridge region outside the HMG box. Indirect immunofluorescence and confocal microscopy show that the mouse SRY colocalizes with KRAB-O in nuclear dots in transiently transfected COS7 cells and primary fetal mouse gonadal cells. Using similar approaches, we demonstrate that KRAB-O interacts directly with KAP1 (KRAB-associated protein 1), the obligatory corepressor for KRAB domain proteins. Furthermore, we show that the mouse SRY is associated indirectly with KAP1 and heterochromatin protein 1 (HP1) through its interaction with KRAB-O, suggesting that the mouse SRY could use the KRAB-KAP1-HP1 organized transcriptional regulatory complex to regulate its yet-to-be-identified downstream target genes.

ECR Browser: a tool for visualizing and accessing data from comparisons of multiple vertebrate genomes

With an increasing number of vertebrate genomes being sequenced in draft or finished form, unique opportunities for decoding the language of DNA sequence through comparative genome alignments have arisen. However, novel tools and strategies are required to accommodate this large volume of genomic information and to facilitate the transfer of predictions generated by comparative sequence alignment to researchers focused on experimental annotation of genome function. Here, we present the ECR Browser, a tool that provides easy and dynamic access to whole genome alignments of human, mouse, rat and fish sequences. This web-based tool (http://ecrbrowser.dcode.org) provides the starting point for discovery of novel genes, identification of distant gene regulatory elements and prediction of transcription factor binding sites. The genome alignment portal of the ECR Browser also permits fast and automated alignments of any user-submitted sequence to the genome of choice. The interconnection of the ECR Browser with other DNA sequence analysis tools creates a unique portal for studying and exploring vertebrate genomes.

N/A

N/A

The DNA sequence and biology of human chromosome 19

Chromosome 19 has the highest gene density of all human chromosomes, more than double the genome-wide average. The large clustered gene families, corresponding high G + C content, CpG islands and density of repetitive DNA indicate a chromosome rich in biological and evolutionary significance. Here we describe 55.8 million base pairs of highly accurate finished sequence representing 99.9% of the euchromatin portion of the chromosome. Manual curation of gene loci reveals 1,461 protein-coding genes and 321 pseudogenes. Among these are genes directly implicated in mendelian disorders, including familial hypercholesterolaemia and insulin-resistant diabetes. Nearly one-quarter of these genes belong to tandemly arranged families, encompassing more than 25% of the chromosome. Comparative analyses show a fascinating picture of conservation and divergence, revealing large blocks of gene orthology with rodents, scattered regions with more recent gene family expansions and deletions, and segments of coding and non-coding conservation with the distant fish species Takifugu.

Duplicate, decouple, disperse: the evolutionary transience of human centromeric regions

Human centromeric regions are enriched for segmental duplications, which elsewhere in the genome precipitate both genetic disease and gene formation. Molecular cytogenetic analyses of primate chromosomes have established that centromeres frequently move without altering the surrounding gene order. Recently, the positions of two ancestral centromeres have been mapped to regions of the human genome that are both rich in segmental duplications and are associated with duplication-based clinical phenotypes. This suggests a model for the evolution of euchromatic segmental duplication families involving the localised elevation of recombination rates within the duplication-rich heterochromatin of recently inactivated centromeres, and raises the possibility that the distribution of duplication/deletion syndromes within our genome has been heavily influenced by such events. The relaxation of the heterochromatin environment that must accompany centromere inactivation would also increase the transcriptional activity within previously pericentromeric DNA, increasing the likelihood of chimaeric gene creation through pericentromeric-directed duplication events.

Sp1- and Krüppel-like transcription factors

Sp1-like proteins and Krüppel-like factors (KLFs) are highly related zinc-finger proteins that are important components of the eukaryotic cellular transcriptional machinery. By regulating the expression of a large number of genes that have GC-rich promoters, Sp1-like/KLF transcription regulators may take part in virtually all facets of cellular function, including cell proliferation, apoptosis, differentiation, and neoplastic transformation. Individual members of the Sp1-like/KLF family can function as activators or repressors depending on which promoter they bind and the coregulators with which they interact. A long-standing research aim has been to define the mechanisms by which Sp1-like factors and KLFs regulate gene expression and cellular function in a cell- and promoter-specific manner. Most members of this family have been identified in mammals, with at least 21 Sp1-like/KLF proteins encoded in the human genome, and members are also found in frogs, worms and flies. Sp1-like/KLF proteins have highly conserved carboxy-terminal zinc-finger domains that function in DNA binding. The amino terminus, containing the transcription activation domain, can vary significantly between family members.

TRAP: Tandem Repeat Assembly Program produces improved shotgun assemblies of repetitive sequences

The software commonly used for assembly of shotgun sequence data has several limitations. One such limitation becomes obvious when repetitive sequences are encountered. Shotgun assembly is a difficult task, even for non-repetitive regions, but the use of quality assessments of the data and efficient matching algorithms have made it possible to assemble most sequences efficiently. In the case of highly repetitive sequences, however, these algorithms fail to distinguish between sequencing errors and single base differences in regions containing nearly identical repeats. None of the currently available fragment assembly programs are able to correctly assemble highly similar repetitive data, and we, therefore, present a novel shotgun assembly program, Tandem Repeat Assembly Program (TRAP). The main feature of this program is the ability to separate long repetitive regions from each other by distinguishing single base substitutions as well as insertions/deletions from sequencing errors. This is accomplished by using a novel multiple-alignment based analysis method. Since repeats are a common complication in most sequencing projects, this software should be of use for the whole sequencing community.

KRAB zinc finger proteins: an analysis of the molecular mechanisms governing their increase in numbers and complexity during evolution

Krüppel-related zinc finger proteins, with 564 members in the human genome, probably constitute the largest individual family of transcription factors in mammals. Approximately 30% of these proteins carry a potent repressor domain called the Krüppel associated box (KRAB). Depending on the structure of the KRAB domain, these proteins have been further divided into three subfamilies (A + B, A + b, and A only). In addition, some KRAB zinc finger proteins contain another conserved motif called SCAN. To study their molecular evolution, an extensive comparative analysis of a large panel of KRAB zinc finger genes was performed. The results show that both the KRAB A + b and the KRAB A subfamilies have their origin in a single member or a few closely related members of the KRAB A + B family. The KRAB A + B family is also the most prevalent among the KRAB zinc finger genes. Furthermore, we show that internal duplications of individual zinc finger motifs or blocks of several zinc finger motifs have occurred quite frequently within this gene family. However, zinc finger motifs are also frequently lost from the open reading frame, either by functional inactivation by point mutations or by the introduction of a stop codon. The introduction of a stop codon causes the exclusion of part of the zinc finger region from the coding region and the formation of graveyards of degenerate zinc finger motifs in the 3'-untranslated region of these genes. Earlier reports have shown that duplications of zinc finger genes commonly occur throughout evolution. We show that there is a relatively low degree of sequence conservation of the zinc finger motifs after these duplications. In many cases this may cause altered binding specificities of the transcription factors encoded by these genes. The repetitive nature of the zinc finger region and the structural flexibility within the zinc finger motif make these proteins highly adaptable. These factors may have been of major importance for their massive expansion in both number and complexity during metazoan evolution.

Characterization of ZNF333, a novel double KRAB domain containing zinc finger gene on human chromosome 19p13.1

ZNF333 is a novel human KRAB-zinc finger protein gene on chromosome 19p13.1 encompassing 14 exons. ZNF333 is highly expressed in heart and encodes a 665 amino acid protein that contains a rare combination of double KRAB-domains, each consisting of a classical KRAB-A and a highly divergent KRAB-B box at the N-terminus. ZNF333 further contains 10 C2H2 zinc finger motifs at the C-terminus.

Human-specific duplication and mosaic transcripts: the recent paralogous structure of chromosome 22

In recent decades, comparative chromosomal banding, chromosome painting, and gene-order studies have shown strong conservation of gross chromosome structure and gene order in mammals. However, findings from the human genome sequence suggest an unprecedented degree of recent (<35 million years ago) segmental duplication. This dynamism of segmental duplications has important implications in disease and evolution. Here we present a chromosome-wide view of the structure and evolution of the most highly homologous duplications (> or = 1 kb and > or = 90%) on chromosome 22. Overall, 10.8% (3.7/33.8 Mb) of chromosome 22 is duplicated, with an average sequence identity of 95.4%. To organize the duplications into tractable units, intron-exon structure and well-defined duplication boundaries were used to define 78 duplicated modules (minimally shared evolutionary segments) with 157 copies on chromosome 22. Analysis of these modules provides evidence for the creation or modification of 11 novel transcripts. Comparative FISH analyses of human, chimpanzee, gorilla, orangutan, and macaque reveal qualitative and quantitative differences in the distribution of these duplications--consistent with their recent origin. Several duplications appear to be human specific, including a approximately 400-kb duplication (99.4%-99.8% sequence identity) that transposed from chromosome 14 to the most proximal pericentromeric region of chromosome 22. Experimental and in silico data further support a pericentromeric gradient of duplications where the most recent duplications transpose adjacent to the centromere. Taken together, these data suggest that segmental duplications have been an ongoing process of primate genome evolution, contributing to recent gene innovation and the dynamic transformation of genome architecture within and among closely related species.

The myosin light chain kinase gene is not duplicated in mouse: partial structure and chromosomal localization of Mylk

The gene encoding myosin light chain kinase (MYLK) is duplicated on human chromosome 3 (HSA3; 3p13;3q21) and on a chromosome with conserved synteny to HSA3 in most non-human primate species. In human, the functional copy resides on 3q21, whereas the 3p13 site contains a pseudogene. To trace the origin of the duplication, we characterized the mouse gene Mylk. A single sequence corresponding to the functional Mylk was detected. We sequenced a 180-kb bacterial artificial chromosome clone containing the 24 first exons of Mylk; the complete mouse gene is expected to span >200 kb. Comparisons with the draft of the human genome revealed that the sequence and structure of MYLK are conserved in mammals. Fluorescence in situ hybridization (FISH) analysis indicated that the mouse gene localizes to a single site on chromosome 16B4-B5, a region with conserved synteny with HSA3q. Our study provides information on both the structure and the evolution of MYLK in mammals and suggests that it was duplicated after the divergence of rodents and primates.

Homology-driven assembly of a sequence-ready mouse BAC contig map spanning regions related to the 46-Mb gene-rich euchromatic segments of human chromosome 19

Draft sequence derived from the 46-Mb gene-rich euchromatic portion of human chromosome 19 (HSA19) was utilized to generate a sequence-ready physical map spanning homologous regions of mouse chromosomes. Sequence similarity searches with the human sequence identified more than 1000 individual orthologous mouse genes from which 382 overgo probes were developed for hybridization. Using human gene order and spacing as a model, these probes were used to isolate and assemble bacterial artificial chromosome (BAC) clone contigs spanning homologous mouse regions. Each contig was verified, extended, and joined to neighboring contigs by restriction enzyme fingerprinting analysis. Approximately 3000 mouse BACs were analyzed and assembled into 44 contigs with a combined length of 41.4 Mb. These BAC contigs, covering 90% of HSA19-related mouse DNA, are distributed throughout 15 homology segments derived from different regions of mouse chromosomes 7, 8, 9, 10, and 17. The alignment of the HSA19 map with the ordered mouse BAC contigs revealed a number of structural differences in several overtly conserved homologous regions and more precisely defined the borders of the known regions of HSA19-syntenic homology. Our results demonstrate that given a human draft sequence, BAC contig maps can be constructed quickly for comparative sequencing without the need for preestablished mouse-specific genetic or physical markers and indicate that similar strategies can be applied with equal success to genomes of other vertebrate species.

Insertion of beta-satellite repeats identifies a transmembrane protease causing both congenital and childhood onset autosomal recessive deafness

Approximately 50% of childhood deafness is caused by mutations in specific genes. Autosomal recessive loci account for approximately 80% of nonsyndromic genetic deafness. Here we report the identification of a new transmembrane serine protease (TMPRSS3; also known as ECHOS1) expressed in many tissues, including fetal cochlea, which is mutated in the families used to describe both the DFNB10 and DFNB8 loci. An 8-bp deletion and insertion of 18 monomeric (approximately 68-bp) beta-satellite repeat units, normally present in tandem arrays of up to several hundred kilobases on the short arms of acrocentric chromosomes, causes congenital deafness (DFNB10). A mutation in a splice-acceptor site, resulting in a 4-bp insertion in the mRNA and a frameshift, was detected in childhood onset deafness (DFNB8). This is the first description of beta-satellite insertion into an active gene resulting in a pathogenic state, and the first description of a protease involved in hearing loss.

The family of Caenorhabditis elegans tyrosine kinase receptors: similarities and differences with mammalian receptors

Transmembrane receptors with tyrosine kinase activity (RTK) constitute a superfamily of proteins present in all metazoans that is associated with the control and regulation of cellular processes. They have been the focus of numerous studies and are a good subject for comparative analyses of multigene families in different species aimed at understanding metazoan evolution. The sequence of the genome of the nematode worm Caenorhabditis elegans is available. This offers a good opportunity to study the superfamily of nematode RTKs in its entirety and to compare it with its mammalian counterpart. We show that the C. elegans RTKs constitute various groups with different phylogenetic relationships with mammalian RTKs. A group of four RTKs show structural similarity with the three mammalian receptors for the vascular endothelial growth factors. Another group comprises RTKs with a short extracellular region, a feature not known in mammals; the genes encoding these RTKs are clustered on chromosome II with other gene families, including genes encoding chitinase-like proteins. Most of the C. elegans RTKs have no direct orthologous relationship with any mammalian RTK, providing an illustration of the importance of the separate evolution of the different phyla.