LINE-1 elements and X chromosome inactivation: a function for "junk" DNA?

Polyamine Dysregulation and Nucleolar Disruption in Alzheimer's Disease

A hypothesis of Alzheimer's disease etiology is proposed describing how cellular stress induces excessive polyamine synthesis and recycling which can disrupt nucleoli. Polyamines are essential in nucleolar functions, such as RNA folding and ribonucleoprotein assembly. Changes in the nucleolar pool of anionic RNA and cationic polyamines acting as counterions can cause significant nucleolar dynamics. Polyamine synthesis reduces S-adenosylmethionine which, at low levels, triggers tau phosphorylation. Also, polyamine recycling reduces acetyl-CoA needed for acetylcholine, which is low in Alzheimer's disease. Extraordinary nucleolar expansion and/or contraction can disrupt epigenetic control in peri-nucleolar chromatin, such as chromosome 14 with the presenilin-1 gene; chromosome 21 with the amyloid precursor protein gene; chromosome 17 with the tau gene; chromosome 19 with the APOE4 gene; and the inactive X chromosome (Xi; aka "nucleolar satellite") with normally silent spermine synthase (polyamine synthesis) and spermidine/spermine-N1-acetyltransferase (polyamine recycling) alleles. Chromosomes 17, 19 and the Xi have high concentrations of Alu elements which can be transcribed by RNA polymerase III if positioned nucleosomes are displaced from the Alu elements. A sudden flood of Alu RNA transcripts can competitively bind nucleolin which is usually bound to Alu sequences in structural RNAs that stabilize the nucleolar heterochromatic shell. This Alu competition leads to loss of nucleolar integrity with leaking of nucleolar polyamines that cause aggregation of phosphorylated tau. The hypothesis was developed with key word searches (e.g., PubMed) using relevant terms (e.g., Alzheimer's, lupus, nucleolin) based on a systems biology approach and exploring autoimmune disease tautology, gaining synergistic insights from other diseases.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

Graves’ Disease in a Young Patient With Turner’s Syndrome: The Genetic Association

INTRODUCTION

Autoimmune diseases occur more often in females, suggesting a key role for the X chromosome. Curiously, individuals with Turner syndrome (TS), with fewer copies of X-linked genes, are prone to develop autoimmune conditions. Hashimoto's thyroiditis (HT) is described with a relatively high frequency in patients with TS while the association with Graves' disease (GD) is rare. Here we report a rare case of TS with GD in a young patient.

METHOD

A 14-year-old girl presented with hyperthyroid symptoms and eye signs that developed over the past six months. She had somatic stigmata of TS. TS was diagnosed by karyotyping (45,XO/46,XX del Xq22) and GD was diagnosed by a thyroid function test and the presence of autoantibodies. She was treated effectively with carbimazole for GD. Estrogen replacement therapy was also initiated to induce the development of secondary sex characteristics.

CONCLUSION

X chromosome inactivation, an epigenetic process that establishes and maintains dosage compensation of X-linked genes, is especially vulnerable to disruption and may contribute to an autoimmune disease process. The occurrence of autoimmune diseases in patients with TS is discussed with regard to possible abnormalities in X-linked dosage compensation.

Accumulation of retrotransposons contributes to W chromosome differentiation in the willow beauty Peribatodes rhomboidaria (Lepidoptera: Geometridae)

The W chromosome of Lepidoptera is typically gene-poor, repeat-rich and composed of heterochromatin. Pioneering studies investigating this chromosome reported an abundance of mobile elements. However, the actual composition of the W chromosome varies greatly between species, as repeatedly demonstrated by comparative genomic hybridization (CGH) or fluorescence in situ hybridization (FISH). Here we present an analysis of repeats on the W chromosome in the willow beauty, Peribatodes rhomboidaria (Geometridae), a species in which CGH predicted an abundance of W-enriched or W-specific sequences. Indeed, comparative analysis of male and female genomes using RepeatExplorer identified ten putative W chromosome-enriched repeats, most of which are LTR or LINE mobile elements. We analysed the two most abundant: PRW LINE-like and PRW Bel-Pao. The results of FISH mapping and bioinformatic analysis confirmed their enrichment on the W chromosome, supporting the hypothesis that mobile elements are the driving force of W chromosome differentiation in Lepidoptera. As the W chromosome is highly underrepresented in chromosome-level genome assemblies of Lepidoptera, this recently introduced approach, combining bioinformatic comparative genome analysis with molecular cytogenetics, provides an elegant tool for studying this elusive and rapidly evolving part of the genome.

X-linked genes exhibit miR6891-5p-regulated skewing in Sjögren's syndrome

Many autoimmune diseases exhibit a strikingly increased prevalence in females, with primary Sjögren's syndrome (pSS) being the most female-predominant example. However, the molecular basis underlying the female-bias in pSS remains elusive. To address this knowledge gap, we performed genome-wide, allele-specific profiling of minor salivary gland-derived mesenchymal stromal cells (MSCs) from pSS patients and control subjects, and detected major differences in the regulation of X-linked genes. In control female MSCs, X-linked genes were expressed from both paternal and maternal X chromosomes with a median paternal ratio of ~ 0.5. However, in pSS female MSCs, X-linked genes exhibited preferential expression from one of the two X chromosomes. Concomitantly, pSS MSCs showed decrease in XIST levels and reorganization of H3K27me3+ foci in the nucleus. Moreover, the HLA-locus-expressed miRNA miR6891-5p was decreased in pSS MSCs. miR6891-5p inhibition in control MSCs caused XIST dysregulation, ectopic silencing, and allelic skewing. Allelic skewing was accompanied by the mislocation of protein products encoded by the skewed genes, which was recapitulated by XIST and miR6891-5p disruption in control MSCs. Our data reveal X skewing as a molecular hallmark of pSS and highlight the importance of restoring X-chromosomal allelic balance for pSS treatment. KEY MESSAGES: X-linked genes exhibit skewing in primary Sjögren's syndrome (pSS). X skewing in pSS associates with alterations in H3K27me3 deposition. pSS MSCs show decreased levels of miR6891-5p, a HLA-expressed miRNA. miR6891-5p inhibition causes H3K27me3 dysregulation and allelic skewing.

Off-Target Effect of Activation of NF-κB by HIV Latency Reversal Agents on Transposable Elements Expression

Many drugs have been evaluated to reactivate HIV-1 from cellular reservoirs, but the off-target effects of these latency reversal agents (LRA) remain poorly defined. Transposable elements (TEs) are reactivated during HIV-1 infection, but studies of potential off-target drug effects on TE expression have been limited. We analyzed the differential expression of TEs induced by canonical and non-canonical NF-κB signaling. We evaluated the effect of PKC agonists (Bryostatin and Ingenol B) on the expression of TEs in memory CD4+ T cells. Ingenol B induced 38 differentially expressed TEs (17 HERV (45%) and 21 L1 (55%)). Interestingly, TE expression in effector memory CD4+ T cells was more affected by Bryostatin compared to other memory T-cell subsets, with 121 (107 upregulated and 14 downregulated) differentially expressed (DE) TEs. Of these, 31% (n = 37) were HERVs, and 69% (n = 84) were LINE-1 (L1). AZD5582 induced 753 DE TEs (406 HERV (54%) and 347 L1 (46%)). Together, our findings show that canonical and non-canonical NF-κB signaling activation leads to retroelement expressions as an off-target effect. Furthermore, our data highlights the importance of exploring the interaction between LRAs and the expression of retroelements in the context of HIV-1 eradication strategies.

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

A comprehensive analysis of gorilla-specific LINE-1 retrotransposons

BACKGROUND

Long interspersed element-1 (LINE-1 or L1) is the most abundant retrotransposons in the primate genome. They have approximately 520,000 copies and make up ~ 17% of the primate genome. Full-length L1s can mobilize to a new genomic location using their enzymatic machinery. Gorilla is the second closest species to humans after the chimpanzee, and human-gorilla split 7-12 million years ago. The gorilla genome provides an opportunity to explore primate origins and evolution.

OBJECTIVE

L1s have contributed to genome diversity and variations during primate evolution. This study aimed to identify gorilla-specific L1s using a more recent version of the gorilla reference genome (Mar. 2016 GSMRT3/gorGor5).

METHODS

We collected gorilla-specific L1 candidates through computational analysis and manual inspection. L1Xplorer was used to identify whether full-length gorilla-specific L1s were intact. In addition, to determine the level of sequence conservation between intact fulllength gorilla-specific L1s, two ORFs of intact L1s were aligned with the L1PA2 consensus sequence.

RESULTS

2002 gorilla-specific L1 candidates were identified through computational analysis. Among them, we manually inspected 1,883 gorilla-specific L1s, among which most of them belong to the L1PA2 subfamily and 12 were intact L1s that could influence genomic variations in the gorilla genome. Interestingly, the 12 intact full-length gorilla-specific L1s have 14 highly conserved nonsynonymous mutations, including 6 mutations and 8 mutations in ORF1 and ORF2, respectively. In comparison to the intact full-length chimpanzee-specific L1s and human-specific hot-L1s, two of these in ORF1 (L256F and E293G) were shown as gorilla-specific nonsynonymous mutations.

CONCLUSION

The gorilla-specific L1s may have had significantly affected the gorilla genome to compose a genome different form that of other primates during primate evolution.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

XIST RNA: a window into the broader role of RNA in nuclear chromosome architecture

XIST RNA triggers the transformation of an active X chromosome into a condensed, inactive Barr body and therefore provides a unique window into transitions of higher-order chromosome architecture. Despite recent progress, how XIST RNA localizes and interacts with the X chromosome remains poorly understood. Genetic engineering of XIST into a trisomic autosome demonstrates remarkable capacity of XIST RNA to localize and comprehensively silence that autosome. Thus, XIST does not require X chromosome-specific sequences but operates on mechanisms available genome-wide. Prior results suggested XIST localization is controlled by attachment to the insoluble nuclear scaffold. Our recent work affirms that scaffold attachment factor A (SAF-A) is involved in anchoring XIST, but argues against the view that SAF-A provides a unimolecular bridge between RNA and the chromosome. Rather, we suggest that a complex meshwork of architectural proteins interact with XIST RNA. Parallel work studying the territory of actively transcribed chromosomes suggests that repeat-rich RNA 'coats' euchromatin and may impact chromosome architecture in a manner opposite of XIST A model is discussed whereby RNA may not just recruit histone modifications, but more directly impact higher-order chromatin condensation via interaction with architectural proteins of the nucleus.This article is part of the themed issue 'X-chromosome inactivation: a tribute to Mary Lyon'.

Chromosome mapping of a Tc1-like transposon in species of the catfish Ancistrus

The Tc1 mariner element is widely distributed among organisms and have been already described in different species of fish. The genus Ancistrus (Kner, 1854) has 68 nominal species and is part of an interesting taxonomic and cytogenetic group, as well as presenting a variation of chromosome number, ranging from 2n=34 to 54 chromosomes, and the existence of simple and multiple sex chromosome system and the occurrence of chromosomal polymorphisms involving chromosomes that carry the nucleolus organizer region. In this study, a repetitive element by restriction enzyme, from Ancistrus sp.1 "Flecha" was isolated, which showed similarity with a transposable element Tc1-mariner. Its chromosomal location is distributed in heterochromatic regions and along the chromosomal arms of all specimens covered in this study, confirming the pattern dispersed of this element found in other studies carried out with other species. Thus, this result reinforces the hypothesis that the sequence AnDraI is really a dispersed element isolated. As this isolated sequence showed the same pattern in all species which have different sex chromosomes systems, including in all sex chromosomes, we could know that it is not involved in sex chromosome differentiation.

Sex-Biased Transcriptome of Schistosoma mansoni: Host-Parasite Interaction, Genetic Determinants and Epigenetic Regulators Are Associated with Sexual Differentiation

Background

Among more than 20,000 species of hermaphroditic trematodes, Schistosomatidae are unusual since they have evolved gonochorism. In schistosomes, sex is determined by a female heterogametic system, but phenotypic sexual dimorphism appears only after infection of the vertebrate definitive host. The completion of gonad maturation occurs even later, after pairing. To date, the molecular mechanisms that trigger the sexual differentiation in these species remain unknown, and in vivo studies on the developing schistosomulum stages are lacking. To study the molecular basis of sex determination and sexual differentiation in schistosomes, we investigated the whole transcriptome of the human parasite Schistosoma mansoni in a stage- and sex-comparative manner.

Methodology/ Principal Findings

We performed a RNA-seq on males and females for five developmental stages: cercariae larvae, three in vivo schistosomulum stages and adults. We detected 7,168 genes differentially expressed between sexes in at least one of the developmental stages, and 4,065 of them were functionally annotated. Transcriptome data were completed with H3K27me3 histone modification analysis using ChIP-Seq before (in cercariae) and after (in adults) the phenotypic sexual dimorphism appearance. In this paper we present (i) candidate determinants of the sexual differentiation, (ii) sex-biased players of the interaction with the vertebrate host, and (iii) different dynamic of the H3K27me3 histone mark between sexes as an illustration of sex-biased epigenetic landscapes.

Conclusions/ Significance

Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, genetic determinants and epigenetic regulators. Our results suggest that such combination could lead to the optimized sexual dimorphism of this parasitic species. As S. mansoni is pathogenic for humans, this study represents a promising source of therapeutic targets, providing not only data on the parasite development in interaction with its vertebrate host, but also new insights on its reproductive function.

Parasitic flatworms include more than 20,000 species that are classically hermaphrodites. Among them, the roughly hundred species of Schistosomatidae are intriguing because they are gonochoric. Schistosomes are responsible of the second most important parasitic disease worldwide, and eggs are the main cause of the inflammatory symptoms. Thus, studying the sexual reproduction mechanisms of schistosomes is of particular interest for drug development. Schistosome’s sex is genetically determined by the presence of sex chromosomes: ZZ in males or ZW in females. There is, however, no phenotypic dimorphism in the larval stages: sexual dimorphism appears only in the vertebrate host. In order to understand the molecular mechanisms underlying phenotypic sexual dimorphism, we performed a transcriptome analysis (RNA-Seq) in five different stages of the parasite lifecycle as well as a chromatin status analysis (ChIP-Seq) in the non-differentiated stage cercariae and in the adult differentiated stage, for males and females separately. Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, developmental pathways and epigenetic regulators. Our sex-comparative approach provides therefore new potential therapeutic targets to affect development and sexual reproduction of parasite.

LINE-1 DNA methylation: A potential forensic marker for discriminating monozygotic twins

Discriminating individuals within a pair of monozygotic (MZ) twins using genetic markers remains unresolved. This inability causes problems in criminal or paternity cases involving MZ twins as suspects or alleged fathers. Our previous study showed DNA methylation differences in interspersed repeat sequences such as Alu and LINE-1 within pairs of newborn MZ twins. To further evaluate the possible value of LINE-1 DNA methylation for discriminating MZ twins, this study investigated the LINE-1 DNA methylation of a large number of twins. We collected blood samples and buccal cell samples from 119 pairs of MZ and 57 pairs of dizygotic (DZ) twins. Genomic DNA was extracted and LINE-1 methylation level was detected using bisulfite pyrosequencing. The mean methylation level of the three CpG sites in the blood sample among the 176 unrelated individuals was 76.60% and 70.08% in buccal samples. This difference was significant, indicating the tissue specificity of LINE-1 DNA methylation. Among 119 pairs of MZ twins, 15 pairs could be discriminated according to the difference of CpG methylation level between them, which accounted for 12.61% of total number of MZ pairs. As for DZ twins, 10 pairs had significant differences between two individuals, which accounted for 17.54% of the total 57 DZ pairs. In conclusion, there are global DNA methylation differences within some healthy concordant monozygotic (MZ) twin pairs. LINE-1 DNA methylation might be a potential marker for helping to discriminate individuals within MZ twin pairs, and the tissue specificity must be considered in practice.

LINEs in mice: features, families, and potential roles in early development

Approximately half of the mammalian genome is composed of repetitive elements, including LINE-1 (L1) elements. Because of their potential ability to transpose and integrate into other regions of the genome, their activation represents a threat to genome stability. Molecular pathways have emerged to tightly regulate and repress their transcriptional activity, including DNA methylation, histone modifications, and RNA pathways. It has become evident that Line-L1 elements are evolutionary diverse and dedicated repression pathways have been recently uncovered that discriminate between evolutionary old and young elements, with RNA-directed silencing mechanisms playing a prominent role. During periods of epigenetic reprogramming in development, specific classes of repetitive elements are upregulated, presumably due to the loss of most heterochromatic marks in this process. While we have learnt a lot on the molecular mechanisms that regulate Line-L1 expression over the last years, it is still unclear whether reactivation of Line-L1 after fertilization serves a functional purpose or it is a simple side effect of reprogramming.

Understanding the Complex Circuitry of lncRNAs at the X-inactivation Center and Its Implications in Disease Conditions

Balanced gene expression is a high priority in order to maintain optimal functioning since alterations and variations could result in acute consequences. X chromosome inactivation (X-inactivation) is one such strategy utilized by mammalian species to silence the extra X chromosome in females to uphold a similar level of expression between the two sexes. A functionally versatile class of molecules called long noncoding RNA (lncRNA) has emerged as key regulators of gene expression and plays important roles during development. An lncRNA that is indispensable for X-inactivation is X-inactive specific transcript (Xist), which induces a repressive epigenetic landscape and creates the inactive X chromosome (Xi). With recent advents in the field of X-inactivation, novel positive and negative lncRNA regulators of Xist such as Jpx and Tsix, respectively, have broadened the regulatory network of X-inactivation. Xist expression failure or dysregulation has been implicated in producing developmental anomalies and disease states. Subsequently, reactivation of the Xi at a later stage of development has also been associated with certain tumors. With the recent influx of information about lncRNA biology and advancements in methods to probe lncRNA, we can now attempt to understand this complex network of Xist regulation in development and disease. It has become clear that the presence of an extra set of genes could be fatal for the organism. Only by understanding the precise ways in which lncRNAs function can treatments be developed to bring aberrations under control. This chapter summarizes our current understanding and knowledge with regard to how lncRNAs are orchestrated at the X-inactivation center (Xic), with a special focus on how genetic diseases come about as a consequence of lncRNA dysregulation.

Quadruplex-forming DNA sequences spread by retrotransposons may serve as genome regulators

Transposable elements (TEs) are ubiquitous genome inhabitants in eukaryotes. Increasing evidence shows that TEs are involved in regulatory networks of eukaryotic cells and contribute to genome evolution. Recently, we reported that many plant long-terminal repeat (LTR) retrotransposons contain DNA quadruplex-forming sequences at precise positions inside their LTRs and that quadruplexes are better preserved in evolutionary younger elements. As quadruplexes can modulate molecular processes, quadruplexes found at specific distances upstream and downstream from the endogenous TE promoter can affect transcription of the element. Moreover, quadruplexes found in solo LTRs, as well as in 3' ends of 5'-truncated copies of LINE-1 elements, can affect expression of neighboring genes. Here, we propose that this way retrotransposons can serve as vehicles for spread of DNA quadruplexes. Quadruplexes can thus fulfill a dual regulatory role-to influence both the retrotransposons carrying them and the neighboring host genes, e.g., by direct effect on transcription or by modifying the local chromatin state. Additionally, four-stranded DNA structures may serve as hotspots for recombination-based genome rearrangements.

Plasma S-adenosylmethionine, DNMT polymorphisms, and peripheral blood LINE-1 methylation among healthy Chinese adults in Singapore

BACKGROUND

Global hypomethylation of repetitive DNA sequences is believed to occur early in tumorigenesis. There is a great interest in identifying factors that contribute to global DNA hypomethylation and associated cancer risk. We tested the hypothesis that plasma S-adenosylmethionine (SAM) level alone or in combination with genetic variation in DNA methyltransferases (DNMT1, DNMT3A and DNMT3B) was associated with global DNA methylation extent at long interspersed nucleotide element-1 (LINE-1) sequences.

METHODS

Plasma SAM level and LINE-1 DNA methylation index were measured using stored blood samples collected from 440 healthy Singaporean Chinese adults during 1994-1999. Genetic polymorphisms of 13 loci in DNMT1, DNMT3A and DNMT3B were determined.

RESULTS

LINE-1 methylation index was significantly higher in men than in women (p = 0.001). LINE-1 methylation index was positively associated with plasma SAM levels (p ≤ 0.01), with a plateau at approximately 78% of LINE-1 methylation index (55 nmol/L plasma SAM) in men and 77% methylation index (50 nmol/L plasma SAM) in women. In men only, the T allele of DNMT1 rs21124724 was associated with a statistically significantly higher LINE-1 methylation index (ptrend = 0.001). The DNMT1 rs2114724 genotype modified the association between plasma SAM and LINE-1 methylation index at low levels of plasma SAM in men.

CONCLUSIONS

Circulating SAM level was associated with LINE-1 methylation status among healthy Chinese adults. The DNMT1 genetic polymorphism may exert a modifying effect on the association between SAM and LINE-1 methylation status in men, especially when plasma SAM level is low. Our findings support a link between plasma SAM and global DNA methylation status at LINE-1 sequences.

Balancing the dose in the mouse

Organisms that use a chromosomal basis of sex determination have a problem of gene inequality. In the mouse, this dimorphism is evident by the presence of two X-chromosomes in females, while males have a single X- and a single Y-chromosome. To balance this disparity, one of the two female X-chromosomes is transcriptionally silenced to neutralize the gene dose with the XY male. Dosage compensation in mammals is known as X-chromosome inactivation (XCI) and is a crucial early developmental process. XCI is an example of epigenetics: a phenotype resulting in changes on a chromosome without a change in nucleic acid sequence. Studies in mouse embryology and genetics have answered many questions about the process of balancing the dose. In this chapter, I highlight how the mouse dosage compensates the gene disparity between XX females and XY males in a crucial epigenetic process called X-chromosome inactivation (XCI).

Impact of retrotransposons in pluripotent stem cells

Retrotransposons, which constitute approximately 40% of the human genome, have the capacity to 'jump' across the genome. Their mobility contributes to oncogenesis, evolution, and genomic plasticity of the host genome. Induced pluripotent stem cells as well as embryonic stem cells are more susceptible than differentiated cells to genomic aberrations including insertion, deletion and duplication. Recent studies have revealed specific behaviors of retrotransposons in pluripotent cells. Here, we review recent progress in understanding retrotransposons and provide a perspective on the relationship between retrotransposons and genomic variation in pluripotent stem cells.

Molecular cytogenetic characterization of the Amazon River dolphin Inia geoffrensis

Classical and molecular cytogenetic (18S rDNA, telomeric sequence, and LINE-1 retrotransposon probes) studies were carried out to contribute to an understanding of the organization of repeated DNA elements in the Amazon River dolphin (boto, Inia geoffrensis). Twenty-seven specimens were examined, each presenting 2n = 44 chromosomes, the karyotype formula 12m + 14sm + 6st + 10t + XX/XY, and fundamental number (FN) = 74. C-positive heterochromatin was observed in terminal and interstitial positions, with the occurrence of polymorphism. Interstitial telomeric sequences were not observed. The nucleolar organizer region (NOR) was located at a single site on a smallest autosomal pair. LINE-1 was preferentially distributed in the euchromatin regions, with the greatest accumulation on the X chromosome. Although the karyotype structure in cetaceans is considered to be conserved, the boto karyotype demonstrated significant variations in its formula, heterochromatin distribution, and the location of the NOR compared to other cetacean species. These results contribute to knowledge of the chromosome organization in boto and to a better understanding of karyoevolution in cetaceans.

Spreading of X chromosome inactivation via a hierarchy of defined Polycomb stations

X chromosome inactivation (XCI) achieves dosage balance in mammals by repressing one of two X chromosomes in females. During XCI, the long noncoding Xist RNA and Polycomb proteins spread along the inactive X (Xi) to initiate chromosome-wide silencing. Although inactivation is known to commence at the X-inactivation center (Xic), how it propagates remains unknown. Here, we examine allele-specific binding of Polycomb repressive complex 2 (PRC2) and chromatin composition during XCI and generate a chromosome-wide profile of Xi and Xa (active X) at nucleosome-resolution. Initially, Polycomb proteins are localized to ∼150 strong sites along the X and concentrated predominantly within bivalent domains coinciding with CpG islands (“canonical sites”). As XCI proceeds, ∼4000 noncanonical sites are recruited, most of which are intergenic, nonbivalent, and lack CpG islands. Polycomb sites are depleted of LINE repeats but enriched for SINEs and simple repeats. Noncanonical sites cluster around the ∼150 strong sites, and their H3K27me3 levels reflect a graded concentration originating from strong sites. This suggests that PRC2 and H3K27 methylation spread along a gradient unique to XCI. We propose that XCI is governed by a hierarchy of defined Polycomb stations that spread H3K27 methylation in cis.

Autoimmune diseases and polyamines

Genetics and environmental factors have important roles in autoimmune diseases but neither has given us sufficient understanding of these mysterious diseases. Therefore, we are now looking closer at epigenetics, an interface between genetics and environmental factors. Epigenetics can be defined as reversible heritable changes to chromatin that can alter gene expression without altering the gene's DNA sequence. Methylation of DNA and histones are primary means of epigenetic control. By adding methyl groups to DNA and histones, it can limit accessibility of the underlying gene thereby altering the amount of gene expression. The methyl group is derived from an essential molecule in the cell, S-adenosylmethionine (SAM). However, a group of small molecules called polyamines also require SAM for their synthesis. Polyamines are essential for many cellular functions and polyamine activity is increased in many autoimmune diseases. Presented here is the "polyamine hypothesis" in which increased polyamine synthesis competes with cellular methylation (epigenetic control) for SAM. It is proposed that increased polyamine activity can cause disruption of cellular methylation, which can lead to abnormal expression of previously sequestered genes and disruption of other methylation-dependent cellular processes.

Methylation of L1Hs promoters is lower on the inactive X, has a tendency of being higher on autosomes in smaller genomes and shows inter-individual variability at some loci

LINE-1 repeats account for ∼17% of the human genome. Little is known about their individual methylation patterns, because their repetitive, almost identical sequences make them difficult to be individually targeted. Here, we used bisulfite conversion to study methylation at individual LINE-1 repeats. The loci studied included 39 X-linked loci and 5 autosomal loci. On the X chromosome in women, we found statistically significant less methylation at almost all L1Hs compared with men. Methylation at L1P and L1M did not correlate with the inactivation status of the host DNA, while the majority of L1Hs that were possible to be studied lie in inactivated regions. To investigate whether the male–female differences at L1Hs on the X are linked to the inactivation process itself rather than to a mere influence of gender, we analyzed six of the L1Hs loci on the X chromosome in Turners and Klinefelters which have female and male phenotype, respectively, but with reversed number of X chromosomes. We could confirm that all samples with two X chromosomes are hypomethylated at the L1Hs loci. Therefore, the inactive X is hypomethylated at L1Hs; the latter could play an exclusive role in the X chromosome inactivation process. At autosomal L1Hs, methylation levels showed a correlation tendency between methylation level and genome size, with higher methylation observed at most loci in individuals with one X chromosome and the lowest in XXY individuals. In summary, loci-specific LINE-1 methylation levels show considerable plasticity and depend on genomic position and constitution.

The demoiselle of X-inactivation: 50 years old and as trendy and mesmerising as ever

In humans, sexual dimorphism is associated with the presence of two X chromosomes in the female, whereas males possess only one X and a small and largely degenerate Y chromosome. How do men cope with having only a single X chromosome given that virtually all other chromosomal monosomies are lethal? Ironically, or even typically many might say, women and more generally female mammals contribute most to the job by shutting down one of their two X chromosomes at random. This phenomenon, called X-inactivation, was originally described some 50 years ago by Mary Lyon and has captivated an increasing number of scientists ever since. The fascination arose in part from the realisation that the inactive X corresponded to a dense heterochromatin mass called the “Barr body” whose number varied with the number of Xs within the nucleus and from the many intellectual questions that this raised: How does the cell count the X chromosomes in the nucleus and inactivate all Xs except one? What kind of molecular mechanisms are able to trigger such a profound, chromosome-wide metamorphosis? When is X-inactivation initiated? How is it transmitted to daughter cells and how is it reset during gametogenesis? This review retraces some of the crucial findings, which have led to our current understanding of a biological process that was initially considered as an exception completely distinct from conventional regulatory systems but is now viewed as a paradigm “par excellence” for epigenetic regulation.

The single active X in human cells: evolutionary tinkering personified

All mammals compensate for sex differences in numbers of X chromosomes by transcribing only a single X chromosome in cells of both sexes; however, they differ from one another in the details of the compensatory mechanisms. These species variations result from chance mutations, species differences in the staging of developmental events, and interactions between events that occur concurrently. Such variations, which have only recently been appreciated, do not interfere with the strategy of establishing a single active X, but they influence how it is carried out. In an overview of X dosage compensation in human cells, I point out the evolutionary variations. I also argue that it is the single active X that is chosen, rather than inactive ones. Further, I suggest that the initial events in the process-those that precede silencing of future inactive X chromosomes-include randomly choosing the future active X, most likely by repressing its XIST locus.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

X chromosome inactivation and autoimmunity

Autoimmune diseases appear to have multiple contributing factors including genetics, epigenetics, environmental factors, and aging. The predominance of females among patients with autoimmune diseases suggests possible involvement of the X chromosome and X chromosome inactivation. X chromosome inactivation is an epigenetic event resulting in multiple levels of control for modulation of the expression of X-linked genes in normal female cells such that there remains only one active X chromosome in the cell. The extent of this control is unique among the chromosomes and has the potential for problems when regulation is disrupted. Here we discuss the X chromosome inactivation process and how the X chromosome and X chromosome inactivation may be involved in development of autoimmune disorders.

Repeat performance: how do genome packaging and regulation depend on simple sequence repeats?

Non-coding DNA has consistently increased during evolution of higher eukaryotes. Since the number of genes has remained relatively static during the evolution of complex organisms, it is believed that increased degree of sophisticated regulation of genes has contributed to the increased complexity. A higher proportion of non-coding DNA, including repeats, is likely to provide more complex regulatory potential. Here, we propose that repeats play a regulatory role by contributing to the packaging of the genome during cellular differentiation. Repeats, and in particular the simple sequence repeats, are proposed to serve as landmarks that can target regulatory mechanisms to a large number of genomic sites with the help of very few factors and regulate the linked loci in a coordinated manner. Repeats may, therefore, function as common target sites for regulatory mechanisms involved in the packaging and dynamic compartmentalization of the chromatin into active and inactive regions during cellular differentiation.

X chromosome inactivation and embryonic stem cells

X chromosome inactivation (XCI) is a process required to equalize the dosage of X-encoded genes between female and male cells. XCI is initiated very early during female embryonic development or upon differentiation of female embryonic stem (ES) cells and results in inactivation of one X chromosome in every female somatic cell. The regulation of XCI involves factors that also play a crucial role in ES cell maintenance and differentiation and the XCI process therefore provides a beautiful paradigm to study ES cell biology. In this chapter we describe the important cis and trans acting regulators of XCI and introduce the models that have been postulated to explain initiation of XCI in female cells only. We also discuss the proteins involved in the establishment of the inactive X chromosome and describe the different chromatin modifications associated with the inactivation process. Finally, we describe the potential of mouse and human ES and induced pluripotent stem (iPS) cells as model systems to study the XCI process.

Effects of L1-ORF2 fragments on green fluorescent protein gene expression

The retrotransposon known as long interspersed nuclear element-1 (L1) is 6 kb long, although most L1s in mammalian and other eukaryotic cells are truncated. L1 contains two open reading frames, ORF1 and ORF2, that code for an RNA-binding protein and a protein with endonuclease and reverse transcriptase activities, respectively. In this work, we examined the effects of full length L1-ORF2 and ORF2 fragments on green fluorescent protein gene (GFP) expression when inserted into the pEGFP-C1 vector downstream of GFP. All of the ORF2 fragments in sense orientation inhibited GFP expression more than when in antisense orientation, which suggests that small ORF2 fragments contribute to the distinct inhibitory effects of this ORF on gene expression. These results provide the first evidence that different 280-bp fragments have distinct effects on the termination of gene transcription, and that when inserted in the antisense direction, fragment 280-9 (the 3' end fragment of ORF2) induces premature termination of transcription that is consistent with the effect of ORF2.

A surrogate approach to study the evolution of noncoding DNA elements that organize eukaryotic genomes

Comparative genomics provides a facile way to address issues of evolutionary constraint acting on different elements of the genome. However, several important DNA elements have not reaped the benefits of this new approach. Some have proved intractable to current day sequencing technology. These include centromeric and heterochromatic DNA, which are essential for chromosome segregation as well as gene regulation, but the highly repetitive nature of the DNA sequences in these regions make them difficult to assemble into longer contigs. Other sequences, like dosage compensation X chromosomal sites, origins of DNA replication, or heterochromatic sequences that encode piwi-associated RNAs, have proved difficult to study because they do not have recognizable DNA features that allow them to be described functionally or computationally. We have employed an alternate approach to the direct study of these DNA elements. By using proteins that specifically bind these noncoding DNAs as surrogates, we can indirectly assay the evolutionary constraints acting on these important DNA elements. We review the impact that such "surrogate strategies" have had on our understanding of the evolutionary constraints shaping centromeres, origins of DNA replication, and dosage compensation X chromosomal sites. These have begun to reveal that in contrast to the view that such structural DNA elements are either highly constrained (under purifying selection) or free to drift (under neutral evolution), some of them may instead be shaped by adaptive evolution and genetic conflicts (these are not mutually exclusive). These insights also help to explain why the same elements (e.g., centromeres and replication origins), which are so complex in some eukaryotic genomes, can be simple and well defined in other where similar conflicts do not exist.

The role of LINEs and CpG islands in dosage compensation on the chicken Z chromosome

Most avian Z genes are expressed more highly in ZZ males than ZW females, suggesting that chromosome-wide mechanisms of dosage compensation have not evolved. Nevertheless, a small percentage of Z genes are expressed at similar levels in males and females, an indication that a yet unidentified mechanism compensates for the sex difference in copy number. Primary DNA sequences are thought to have a role in determining chromosome gene inactivation status on the mammalian X chromosome. However, it is currently unknown whether primary DNA sequences also mediate chicken Z gene compensation status. Using a combination of chicken DNA sequences and Z gene compensation profiles of 310 genes, we explored the relationship between Z gene compensation status and primary DNA sequence features. Statistical analysis of different Z chromosomal features revealed that long interspersed nuclear elements (LINEs) and CpG islands are enriched on the Z chromosome compared with 329 other DNA features. Linear support vector machine (SVM) classifiers, using primary DNA sequences, correctly predict the Z compensation status for >60% of all Z-linked genes. CpG islands appear to be the most accurate classifier and alone can correctly predict compensation of 63% of Z genes. We also show that LINE CR1 elements are enriched 2.7-fold on the chicken Z chromosome compared with autosomes and that chicken chromosomal length is highly correlated with percentage LINE content. However, the position of LINE elements is not significantly associated with dosage compensation status of Z genes. We also find a trend for a higher proportion of CpG islands in the region of the Z chromosome with the fewest dosage-compensated genes compared with the region containing the greatest concentration of compensated genes. Comparison between chicken and platypus genomes shows that LINE elements are not enriched on sex chromosomes in platypus, indicating that LINE accumulation is not a feature of all sex chromosomes. Our results suggest that CpG islands are not randomly distributed on the Z chromosome and may influence Z gene dosage compensation status.

The role of repetitive DNA in structure and evolution of sex chromosomes in plants

Eukaryotic genomes contain a large proportion of repetitive DNA sequences, mostly transposable elements (TEs) and tandem repeats. These repetitive sequences often colonize specific chromosomal (Y or W chromosomes, B chromosomes) or subchromosomal (telomeres, centromeres) niches. Sex chromosomes, especially non-recombining regions of the Y chromosome, are subject to different evolutionary forces compared with autosomes. In non-recombining regions of the Y chromosome repetitive DNA sequences are accumulated, representing a dominant and early process forming the Y chromosome, probably before genes start to degenerate. Here we review the occurrence and role of repetitive DNA in Y chromosome evolution in various species with a focus on dioecious plants. We also discuss the potential link between recombination and transposition in shaping genomes.

Quantitative variation of LINE-1 sequences in five species and three subspecies of the subgenus Mus and in five Robertsonian races of Mus musculus domesticus

The quantitative variation of a conserved region of the LINE-1 ORF2 sequence was determined in eight species and subspecies of the subgenus Mus (M. m. domesticus, M. m. musculus, M. m. castaneus, M. spicilegus, M. spretus, M. cervicolor, M. cookii, M. caroli) and five Robertsonian races of M. m. domesticus. No differences in LINE-1 ORF2 content were found between all acrocentric or Robertsonian chromosome races, whereas the quantitative variation of the LINE-1 ORF2 sequences detected among the eight taxa partly matches with the clades into which the subgenus is divided. An accumulation of LINE-1 ORF2 elements likely occurred during the evolution of the subgenus. Within the Asiatic clade, M. cervicolor, cookii, and caroli show a low quantity of LINE-1 sequences, also detected within the Palearctic clade in M. m. castaneus and M. spretus, representing perhaps the ancestral condition within the subgenus. On the other hand, M. m. domesticus, M. m. musculus and M. spicilegus showed a high content of LINE-1 ORF2 sequences. Comparison between the chromosomal hybridization pattern of M. m. domesticus, which possesses the highest content, and M. spicilegus did not show any difference in the LINE-1 ORF2 distribution, suggesting that the quantitative variation of this sequence family did not involve chromosome restructuring or a preferential chromosome accumulation, during the evolution of M. m. domesticus.

Retroelements (LINEs and SINEs) in vole genomes: differential distribution in the constitutive heterochromatin

The chromosomal distribution of mobile genetic elements is scarcely known in Arvicolinae species, but could be of relevance to understand the origin and complex evolution of the sex chromosome heterochromatin. In this work we cloned two retrotransposon sequences, L1 and SINE-B1, from the genome of Chionomys nivalis and investigated their chromosomal distribution on several arvicoline species. Our results demonstrate first that both retroelements are the most abundant repeated DNA sequences in the genome of these species. L1 elements, in most species, are highly accumulated in the sex chromosomes compared to the autosomes. This favoured L1 insertion could have played an important role in the origin of the enlarged heterochromatic blocks existing in the sex chromosomes of some Microtus species. Also, we propose that L1 accumulation on the X heterochromatin could have been the consequence of different, independent and rapid amplification processes acting in each species. SINE elements, however, were completely lacking from the constitutive heterochromatin, either in autosomes or in the heterochromatic blocks of sex chromosomes. These data could indicate that some SINE elements are incompatible with the formation of heterochromatic complexes and hence are necessarily missing from the constitutive heterochromatin.

Identification and characterization of cryptic SHOX intragenic deletions in three Japanese patients with Léri-Weill dyschondrosteosis

Although short-stature homeobox-containing gene (SHOX ) haploinsufficiency is responsible for Léri-Weill dyschondrosteosis (LWD), the molecular defect has not been identified in approximately 20% of Japanese LWD patients. Furthermore, although high prevalence of microdeletions affecting SHOX is primarily ascribed to the presence of repeat sequences such as Alu elements around SHOX, it remains to be determined whether microdeletions are actually mediated by repeat sequences. We performed multiple ligation probe amplification (MLPA) assay in six Japanese LWD patients with apparently normal SHOX, followed by fluorescent in situ hybridization (FISH) analysis and sequencing for polymerase chain reaction (PCR) products encompassing the deletion junctions in patients with abnormal MLPA patterns. Consequently, heterozygous intragenic deletions were identified in three cases, i.e., a 5,906-bp deletion involving exons 4-5 in case 1, a 5,594-bp deletion involving exons 4-6a in case 2, and a 50,199-bp deletion involving exons 4-6b in case 3. The deletion breakpoints of cases 1 and 2 were present in nonrepeat sequences, whereas those of case 3 resided within Alu elements. The results suggest that cryptic SHOX intragenic deletions account for a small fraction of LWD and that microdeletions affecting SHOX can be generated by repeat-sequence-mediated aberrant recombinations and by nonhomologous end joining.

Species-specific positive selection of the male-specific lethal complex that participates in dosage compensation in Drosophila

In many taxa, males and females have unequal ratios of sex chromosomes to autosomes, which has resulted in the invention of diverse mechanisms to equilibrate gene expression between the sexes (dosage compensation). Failure to compensate for sex chromosome dosage results in male lethality in Drosophila. In Drosophila, a male-specific lethal (MSL) complex of proteins and noncoding RNAs binds to hundreds of sites on the single male X chromosome and up-regulates gene expression. Here we use population genetics of two closely related Drosophila species to show that adaptive evolution has occurred in all five protein-coding genes of the MSL complex. This positive selection is asymmetric between closely related species, with a very strong signature apparent in Drosophila melanogaster but not in Drosophila simulans. In particular, the MSL1 and MSL2 proteins have undergone dramatic positive selection in D. melanogaster, in domains previously shown to be responsible for their specific targeting to the X chromosome. This signature of positive selection at an essential protein-DNA interface of the complex is unexpected and suggests that X chromosomal MSL-binding DNA segments may themselves be changing rapidly. This highly asymmetric, rapid evolution of the MSL genes further suggests that misregulated dosage compensation may represent one of the underlying causes of male hybrid inviability in Drosophila, wherein the fate of hybrid males depends on which species' X chromosome is inherited.

Sequence and characterization of the Ig heavy chain constant and partial variable region of the mouse strain 129S1

Although the entire mouse genome has been sequenced, there remain challenges concerning the elucidation of particular complex and polymorphic genomic loci. In the murine Igh locus, different haplotypes exist in different inbred mouse strains. For example, the Igh(b) haplotype sequence of the Mouse Genome Project strain C57BL/6 differs considerably from the Igh(a) haplotype of BALB/c, which has been widely used in the analyses of Ab responses. We have sequenced and annotated the 3' half of the Igh(a) locus of 129S1/SvImJ, covering the C(H) region and approximately half of the V(H) region. This sequence comprises 128 V(H) genes, of which 49 are judged to be functional. The comparison of the Igh(a) sequence with the homologous Igh(b) region from C57BL/6 revealed two major expansions in the germline repertoire of Igh(a). In addition, we found smaller haplotype-specific differences like the duplication of five V(H) genes in the Igh(a) locus. We generated a V(H) allele table by comparing the individual V(H) genes of both haplotypes. Surprisingly, the number and position of D(H) genes in the 129S1 strain differs not only from the sequence of C57BL/6 but also from the map published for BALB/c. Taken together, the contiguous genomic sequence of the 3' part of the Igh(a) locus allows a detailed view of the recent evolution of this highly dynamic locus in the mouse.

FluoMEP: a new genotyping method combining the advantages of randomly amplified polymorphic DNA and amplified fragment length polymorphism

PCR-based identification of differences between two unknown genomes often requires complex manipulation of the templates prior to amplification and/or gel electrophoretic separation of a large number of samples with manual methods. Here, we describe a new genotyping method, called fluorescent motif enhanced polymorphism (fluoMEP). The fluoMEP method is based on random amplified polymorphic DNA (RAPD) assay, but combines the advantages of the large collection of unlabelled 10mer primers (ca. 5000) from commercial sources and the power of the automated CE devices used for the detection of amplified fragment length polymorphism (AFLP) patterns. The link between these two components is provided by a fluorescently labeled "common primer" that is used in a two-primer PCR together with an unlabeled RAPD primer. By using the same "common primer" and a series of RAPD primers, DNA templates can be screened quickly and effectively for polymorphisms. Our manuscript describes the optimization of the method and its characterization on different templates. We demonstrate by using several different approaches that the addition of the "common primer" to the PCR changes the profile of amplified fragments, allowing for screening various parts of the genome with the same set of unlabeled primers. We also present an in silico analysis of the genomic localization of fragments amplified by a RAPD primer with two different "common primers" and alone.

Turning junk into gold: domestication of transposable elements and the creation of new genes in eukaryotes

Autonomous transposable elements, generally considered as junk and selfish, encode transposition proteins that can bind, copy, break, join or degrade nucleic acids as well as process or interact with other proteins. Such a repertoire of activities might be of interest for the host cell. There is indeed substantial evidence that mobile DNA can serve as a dynamic reservoir for new cellular functions. Transposable element genes encoding transposase, integrase, reverse transcriptase as well as structural and envelope proteins have been repeatedly recruited by their host during evolution in most eukaryotic lineages. Such domesticated sequences protect us against infections, are necessary for our reproduction, allow the replication of our chromosomes and control cell proliferation and death; others are essential for plant development. Many new candidates for domesticated sequences have been revealed by sequencing projects. Their functional analysis will uncover new aspects of evolutionary alchemy, the turning of junk into gold within genomes.

Genetic resources, genome mapping and evolutionary genomics of the pig (Sus scrofa)

The pig, a representative of the artiodactyla clade, is one of the first animals domesticated, and has become an important agriculture animal as one of the major human nutritional sources of animal based protein. The pig is also a valuable biomedical model organism for human health. The pig's importance to human health and nutrition is reflected in the decision to sequence its genome (3X). As an animal species with its wild ancestors present in the world, the pig provides a unique opportunity for tracing mammalian evolutionary history and defining signatures of selection resulting from both domestication and natural selection. Completion of the pig genome sequencing project will have significant impacts on both agriculture and human health. Following the pig whole genome sequence drafts, along with large-scale polymorphism data, it will be possible to conduct genome sweeps using association mapping, and identify signatures of selection. Here, we provide a description of the pig genome sequencing project and perspectives on utilizing genomic technologies to exploit pig genome evolution and the molecular basis for phenotypic traits for improving pig production and health.

Classifying genomic sequences by sequence feature analysis

Traditional sequence analysis depends on sequence alignment. In this study, we analyzed various functional regions of the human genome based on sequence features, including word frequency, dinucleotide relative abundance, and base-base correlation. We analyzed the human chromosome 22 and classified the upstream, exon, intron, downstream, and intergenic regions by principal component analysis and discriminant analysis of these features. The results show that we could classify the functional regions of genome based on sequence feature and discriminant analysis.