Alu insertion polymorphisms for the study of human genomic diversity

A common Alu insertion in the 3'UTR of TMEM106B is associated with risk of dementia

INTRODUCTION

Sequence variants in TMEM106B have been associated with an increased risk of developing dementia.

METHODS

As part of our efforts to generate a set of mouse lines in which we replaced the mouse Tmem106b gene with a human TMEM106B gene comprised of either a risk or protective haplotype, we conducted an in-depth sequence analysis of these alleles. We also analyzed transcribed TMEM106B sequences using RNA-seq data (AD Knowledge portal) and full genome sequences (1000 Genomes).

RESULTS

We identified an AluYb8 insertion in the 3' untranslated region (3'UTR) of the TMEM106B risk haplotype. We found this AluYb8 insertion in every risk haplotype analyzed, but not in either protective haplotypes or in non-human primates.

DISCUSSION

We conclude that this risk haplotype arose early in human development with a single Alu-insertion event within a unique haplotype context. This AluYb8 element may act as a functional variant in conferring an increased risk of developing dementia.

HIGHLIGHTS

We conducted an in-depth sequence analysis of (1) a risk and (2) a protective haplotype of the human TMEM106B gene. We also analyzed transcribed TMEM106B sequences using RNA-seq data (AD Knowledge Portal) and full genome sequences (1000 Genomes). We identified an AluYb8 insertion in the 3' untranslated region (3'UTR) of the TMEM106B risk haplotype. We found this AluYb8 insertion in every risk haplotype analyzed, but not in either protective haplotypes or in non-human primates. This AluYb8 element may act as a functional variant in conferring an increased risk of developing dementia.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Polyamines and Their Metabolism: From the Maintenance of Physiological Homeostasis to the Mediation of Disease

The polyamines spermidine and spermine are positively charged aliphatic molecules. They are critical in the regulation of nucleic acid and protein structures, protein synthesis, protein and nucleic acid interactions, oxidative balance, and cell proliferation. Cellular polyamine levels are tightly controlled through their import, export, de novo synthesis, and catabolism. Enzymes and enzymatic cascades involved in polyamine metabolism have been well characterized. This knowledge has been used for the development of novel compounds for research and medical applications. Furthermore, studies have shown that disturbances in polyamine levels and their metabolic pathways, as a result of spontaneous mutations in patients, genetic engineering in mice or experimentally induced injuries in rodents, are associated with multiple maladaptive changes. The adverse effects of altered polyamine metabolism have also been demonstrated in in vitro models. These observations highlight the important role these molecules and their metabolism play in the maintenance of physiological normalcy and the mediation of injury. This review will attempt to cover the extensive and diverse knowledge of the biological role of polyamines and their metabolism in the maintenance of physiological homeostasis and the mediation of tissue injury.

Poly-Enrich: count-based methods for gene set enrichment testing with genomic regions

Gene set enrichment (GSE) testing enhances the biological interpretation of ChIP-seq data and other large sets of genomic regions. Our group has previously introduced two GSE methods for genomic regions: ChIP-Enrich for narrow regions and Broad-Enrich for broad regions. Here, we introduce Poly-Enrich, which has wider applicability, additional capabilities and models the number of peaks assigned to a gene using a generalized additive model with a negative binomial family to determine gene set enrichment, while adjusting for gene locus length. As opposed to ChIP-Enrich, Poly-Enrich works well even when nearly all genes have a peak, illustrated by using Poly-Enrich to characterize pathways and types of genic regions enriched with different families of repetitive elements. By comparing Poly-Enrich and ChIP-Enrich results with ENCODE ChIP-seq data, we found that the optimal test depends more on the pathway being regulated than on properties of the transcription factors. Using known transcription factor functions, we discovered clusters of related biological processes consistently better modeled with Poly-Enrich. This suggests that the regulation of certain processes may be modified by multiple binding events, better modeled by a count-based method. Our new hybrid method automatically uses the optimal method for each gene set, with correct FDR-adjustment.

True Homoplasy of Retrotransposon Insertions in Primates

How reliable are the presence/absence insertion patterns of the supposedly homoplasy-free retrotransposons, which were randomly inserted in the quasi infinite genomic space? To systematically examine this question in an up-to-date, multigenome comparison, we screened millions of primate transposed Alu SINE elements for incidences of homoplasious precise insertions and deletions. In genome-wide analyses, we identified and manually verified nine cases of precise parallel Alu insertions of apparently identical elements at orthologous positions in two ape lineages and twelve incidences of precise deletions of previously established SINEs. Correspondingly, eight precise parallel insertions and no exact deletions were detected in a comparison of lemuriform primate and human insertions spanning the range of primate diversity. With an overall frequency of homoplasious Alu insertions of only 0.01% (for human-chimpanzee-rhesus macaque) and 0.02-0.04% (for human-bushbaby-lemurs) and precise Alu deletions of 0.001-0.002% (for human-chimpanzee-rhesus macaque), real homoplasy is not considered to be a quantitatively relevant source of evolutionary noise. Thus, presence/absence patterns of Alu retrotransposons and, presumably, all LINE1-mobilized elements represent indeed the virtually homoplasy-free markers they are considered to be. Therefore, ancestral incomplete lineage sorting and hybridization remain the only serious sources of conflicting presence/absence patterns of retrotransposon insertions, and as such are detectable and quantifiable. [Homoplasy; precise deletions; precise parallel insertions; primates; retrotransposons.].

Genetic variation and population structure of five ethnic groups from Punjab, North-West India: Analysis of MHC class I polymorphic Alu insertions (POALINs)

Genetic variation and differentiation of five ethnic groups from Punjab, North-West India was characterized by analyzing data on polymorphic Alu insertions (POALINs) within the class I genomic region of major histocompatibility complex (MHC), which is completely non-existent in Indian population. The haplotype frequency, distribution and heterozygosity among these groups and their potential implications in molecular anthropology and evolutionary studies were also determined. A total of 479 unrelated healthy individuals representing five different ethnic groups: Banias, Brahmins, Khatri, Jat Sikhs and Scheduled Castes were genotyped for five MHC Alu elements (AluHG, AluMICB, AluHJ, AluTF and AluHF) using polymerase chain reaction (PCR). All the loci were found to be polymorphic among the studied populations. No significant deviation from Hardy-Weinberg equilibrium was observed, except for the AluHJ locus in Brahmins. The POALINs varied in allele frequency between 0.0260 and 0.4427. The average heterozygosity among the studied groups ranged from 0.1937 in Banias to 0.2666 in Jat Sikhs. The genetic differentiation among the studied groups was observed to be of the order of 0.01302. Single POALIN haplotypes were found to be more frequent than multiple POALIN haplotypes. The results of inter-population differentiations, haplotype frequencies, genetic distances, multidimensional scaling, phylogenetic and structure analyses indicated close genetic relationships between the five ethnic groups of Punjab, North-West India. Analyses of polymorphic Alu loci of MHC genomic region may represent reliable information about the ancestry, demographic history and geographic origins of the various human populations, facilitating better understanding of the evolutionary, forensic and epidemiological prospective.

AP-1 and the injury response of the GFAP gene

Increased GFAP gene expression is a common feature of CNS injury, resulting in its use as a reporter to investigate mechanisms producing gliosis. AP-1 transcription factors are among those proposed to participate in mediating the reactive response. Prior studies found a consensus AP-1 binding site in the GFAP promoter to be essential for activity of reporter constructs transfected into cultured cells, but to have little to no effect on basal transgene expression in mice. Since cultured astrocytes display some properties of reactive astrocytes, these findings suggested that AP-1 transcription factors are critical for the upregulation of GFAP in injury, but not for its resting level of expression. We have examined this possibility by comparing the injury response in mice of lacZ transgenes driven by human GFAP promoters that contain the wild-type AP-1 binding site to those in which the site is mutated. An intact AP-1 site was found critical for a GFAP promoter response to the three different injury models used: physical trauma produced by cryoinjury, seizures produced by kainic acid, and chronic gliosis produced in an Alexander disease model. An unexpected additional finding was that the responses of the lacZ transgenes driven by the wild-type promoters were substantially less than that of the endogenous mouse GFAP gene. This suggests that the GFAP gene has previously unrecognized injury-responsive elements that reside further upstream of the transcription start site than the 2.2 kb present in the GFAP promoter segments used here.

Apolipoprotein(a) is the Product of a Pseudogene: Implications for the Pathophysiology of Lipoprotein(a)

Apolipoprotein(a) [apo(a)] is an apolipoprotein unique to lipoprotein(a) [Lp(a)]. Although it has no known function, Lp(a) is a risk factor for accelerated atherothrombosis. We hypothesize that LPA, the gene which encodes apo(a), is a heretofore unrecognized unprocessed pseudogene created by duplication of PLG, the gene which encodes plasminogen. Unprocessed pseudogenes are genes which were created by duplication of functional genes and subsequently lost function after acquiring various mutations. This hypothesis explains many of the unusual features of Lp(a) and apo(a). Also, this hypothesis has implications for the therapy of elevated Lp(a) and atherothrombosis theory. Because apo(a) is functionless, the diseases associated with elevated levels of Lp(a) are due to its impact on blood viscosity.

Alu element insertion in PKLR gene as a novel cause of pyruvate kinase deficiency in Middle Eastern patients

Pyruvate kinase deficiency (PKD) is the most frequent red blood cell enzyme abnormality of the glycolytic pathway and the most common cause of hereditary nonspherocytic hemolytic anemia. Over 250 PKLR-gene mutations have been described, including missense/nonsense, splicing and regulatory mutations, small insertions, small and gross deletions, causing PKD and hemolytic anemia of variable severity. Alu retrotransposons are the most abundant mobile DNA sequences in the human genome, contributing to almost 11% of its mass. Alu insertions have been associated with a number of human diseases either by disrupting a coding region or a splice signal. Here, we report on two unrelated Middle Eastern patients, both born from consanguineous parents, with transfusion-dependent hemolytic anemia, where sequence analysis revealed a homozygous insertion of AluYb9 within exon 6 of the PKLR gene, causing precipitous decrease of PKLR RNA levels. This Alu element insertion consists a previously unrecognized mechanism underlying pathogenesis of PKD.

Sequence Analysis and Characterization of Active Human Alu Subfamilies Based on the 1000 Genomes Pilot Project

The goal of the 1000 Genomes Consortium is to characterize human genome structural variation (SV), including forms of copy number variations such as deletions, duplications, and insertions. Mobile element insertions, particularly Alu elements, are major contributors to genomic SV among humans. During the pilot phase of the project we experimentally validated 645 (611 intergenic and 34 exon targeted) polymorphic "young" Alu insertion events, absent from the human reference genome. Here, we report high resolution sequencing of 343 (322 unique) recent Alu insertion events, along with their respective target site duplications, precise genomic breakpoint coordinates, subfamily assignment, percent divergence, and estimated A-rich tail lengths. All the sequenced Alu loci were derived from the AluY lineage with no evidence of retrotransposition activity involving older Alu families (e.g., AluJ and AluS). AluYa5 is currently the most active Alu subfamily in the human lineage, followed by AluYb8, and many others including three newly identified subfamilies we have termed AluYb7a3, AluYb8b1, and AluYa4a1. This report provides the structural details of 322 unique Alu variants from individual human genomes collectively adding about 100 kb of genomic variation. Many Alu subfamilies are currently active in human populations, including a surprising level of AluY retrotransposition. Human Alu subfamilies exhibit continuous evolution with potential drivers sprouting new Alu lineages.

Unity in diversity: an overview of the genomic anthropology of India

CONTEXT

India is considered a treasure for geneticists and evolutionary biologists due to its vast human diversity, consisting of more than 4500 anthropologically well-defined populations (castes, tribes and religious groups). Each population differs in terms of endogamy, language, culture, physical features, geographic and climatic position and genetic architecture. These factors contributed to India-specific genetic variations which may be responsible for various common diseases in India and its migratory populations. As a result, interpretations of the origins and affinities of Indian populations as well as health and disease conditions require complex and sophisticated genetic analysis. Evidence of ancient human dispersals and settlements is preserved in the genome of Indian inhabitants and this has been extensively analysed in conventional and genomic analyses.

OBJECTIVE AND METHODS

Using genomic analyses of STRs and Alu on a set of populations, this study estimates the level and extent of genetic variation and its implications.

RESULTS

The results show that Indian populations have a higher level of unique genetic diversity which is structured by many social processes and geographical attributes of the country.

CONCLUSION

This overview highlights the need to study the anthropological structure and evolutionary history of Indian populations while designing genomic and epigenomic investigations.

Assembly and characterization of novel Alu inserts detected from next-generation sequencing data

Repetitive elements generally, and Alu inserts specifically are a large contributor to the recent evolution of the human genome. By assembling the sequences of novel Alu inserts using their respective subfamily consensus sequences as references, we found an exponential decay in the Alu subfamily call enrichment with increased number of sequence variants (Pearson correlation [Formula: see text], [Formula: see text]). By mapping the sequences of these inserts to a human reference genome, we infer the reference Alu sources of a subset of the novel Alus, of which 85% were previously shown to be active. We also evaluate relationships between the loci of the novel inserts and their inferred sources.

The impact of transposable elements in genome evolution and genetic instability and their implications in various diseases

Approximately 45% of the human genome is comprised of transposable elements (TEs). Results from the Human Genome Project have emphasized the biological importance of TEs. Many studies have revealed that TEs are not simply "junk" DNA, but rather, they play various roles in processes, including genome evolution, gene expression regulation, genetic instability, and cancer disposition. The effects of TE insertion in the genome varies from negligible to disease conditions. For the past two decades, many studies have shown that TEs are the causative factors of various genetic disorders and cancer. TEs are a subject of interest worldwide, not only in terms of their clinical aspects but also in basic research, such as evolutionary tracking. Although active TEs contribute to genetic instability and disease states, non-long terminal repeat transposons are well studied, and their roles in these processes have been confirmed. In this review, we will give an overview of the importance of TEs in studying genome evolution and genetic instability, and we suggest that further in-depth studies on the mechanisms related to these phenomena will be useful for both evolutionary tracking and clinical diagnostics.

Study of Modern Human Evolution via Comparative Analysis with the Neanderthal Genome

Many other human species appeared in evolution in the last 6 million years that have not been able to survive to modern times and are broadly known as archaic humans, as opposed to the extant modern humans. It has always been considered fascinating to compare the modern human genome with that of archaic humans to identify modern human-specific sequence variants and figure out those that made modern humans different from their predecessors or cousin species. Neanderthals are the latest humans to become extinct, and many factors made them the best representatives of archaic humans. Even though a number of comparisons have been made sporadically between Neanderthals and modern humans, mostly following a candidate gene approach, the major breakthrough took place with the sequencing of the Neanderthal genome. The initial genome-wide comparison, based on the first draft of the Neanderthal genome, has generated some interesting inferences regarding variations in functional elements that are not shared by the two species and the debated admixture question. However, there are certain other genetic elements that were not included or included at a smaller scale in those studies, and they should be compared comprehensively to better understand the molecular make-up of modern humans and their phenotypic characteristics. Besides briefly discussing the important outcomes of the comparative analyses made so far between modern humans and Neanderthals, we propose that future comparative studies may include retrotransposons, pseudogenes, and conserved non-coding regions, all of which might have played significant roles during the evolution of modern humans.

Association and differentiation of MHC class I and II polymorphic Alu insertions and HLA-A, -B, -C and -DRB1 alleles in the Chinese Han population

In order to investigate the polymorphism of Alu insertions (POALINs) in the HLA region, we genotyped ten Alu loci (AluMICB, AluTF, AluHJ, AluHG, AluHF in the HLA class I region and AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10 in the HLA class II region) to determine their allele frequencies and associations with the HLA-A, HLA-B, HLA-C and HLA-DRB1 genes in the Chinese Han population. Our results showed the ten-loci POALINs varied in frequency between 0.003 and 0.425. By comparing the data of the ten-loci POALIN in Chinese Han with Japanese and Caucasian data, marked differences were observed between the three ethnic groups at the allelic or haplotypic levels. Each POALIN was in significant linkage disequilibrium with a variety of HLA-A, -B, -C and -DRB1 alleles, and was associated with a variety of HLA-A, -B, -C and -DRB1 allele in Chinese Han. This comparative study of multilocus POALINs in the HLA class I and II regions of the Chinese Han population shows that POALINs alone or as haplotypes together with the HLA class I and II alleles are informative genetic markers for the identification of HLA class I and II allele and variations, such as crossing over events within the same and/or different populations.

Identification of three new Alu Yb subfamilies by source tracking of recently integrated Alu Yb elements

BACKGROUND

Alu elements are the most abundant mobile elements in the human genome, with over 1 million copies and constituting more than 10% of the genome. The majority of these Alu elements were inserted into the primate genome 35 to 60 million years ago, but certain subfamilies of Alu elements are relatively very new and suspected to be still evolving. We attempted to trace the source/master copies of all human-specific members of the Alu Yb lineage using a computational approach by clustering similar Yb elements and constructing an evolutionary relation among the members of a cluster.

RESULTS

We discovered that one copy of Yb8 at 10p14 is the source of several active Yb8 copies, which retrotransposed to generate 712 copies or 54% of all human-specific Yb8 elements. We detected eight other Yb8 elements that had generated ten or more copies, potentially acting as 'stealth drivers'. One Yb8 element at 14q32.31 seemed to act as the source copy for all Yb9 elements tested, having producing 13 active Yb9 elements, and subsequently generated a total of 131 full-length copies. We identified and characterized three new subclasses of Yb elements: Yb8a1, Yb10 and Yb11. Their copy numbers in the reference genome are 75, 8 and 16. We analysed personal genome data from the 1000 Genome Project and detected an additional 6 Yb8a1, 3 Yb10 and 15 Yb11 copies outside the reference genome. Our analysis indicates that the Yb8a1 subfamily has a similar age to Yb9 (1.93 million years and 2.15 million years, respectively), while Yb10 and Yb11 evolved only 1.4 and 0.71 million years ago, suggesting a linear evolutionary path from Yb8a1 to Yb10 and then to Yb11. Our preliminary data indicate that members in Yb10 and Yb11 are mostly polymorphic, indicating their young age.

CONCLUSIONS

Our findings suggest that the Yb lineage is still evolving with new subfamilies being formed. Due to their very young age and the high rate of being polymorphic, insertions from these young subfamilies are very useful genetic markers for studying human population genetics and migration patterns, and the trend for mobile element insertions in the human genome.

Genome-wide analysis of short interspersed nuclear elements SINES revealed high sequence conservation, gene association and retrotranspositional activity in wheat

Short interspersed nuclear elements (SINEs) are non-autonomous non-LTR retroelements that are present in most eukaryotic species. While SINEs have been intensively investigated in humans and other animal systems, they are poorly studied in plants, especially in wheat (Triticum aestivum). We used quantitative PCR of various wheat species to determine the copy number of a wheat SINE family, termed Au SINE, combined with computer-assisted analyses of the publicly available 454 pyrosequencing database of T. aestivum. In addition, we utilized site-specific PCR on 57 Au SINE insertions, transposon methylation display and transposon display on newly formed wheat polyploids to assess retrotranspositional activity, epigenetic status and genetic rearrangements in Au SINE, respectively. We retrieved 3706 different insertions of Au SINE from the 454 pyrosequencing database of T. aestivum, and found that most of the elements are inserted in A/T-rich regions, while approximately 38% of the insertions are associated with transcribed regions, including known wheat genes. We observed typical retrotransposition of Au SINE in the second generation of a newly formed wheat allohexaploid, and massive hypermethylation in CCGG sites surrounding Au SINE in the third generation. Finally, we observed huge differences in the copy numbers in diploid Triticum and Aegilops species, and a significant increase in the copy numbers in natural wheat polyploids, but no significant increase in the copy number of Au SINE in the first four generations for two of three newly formed allopolyploid species used in this study. Our data indicate that SINEs may play a prominent role in the genomic evolution of wheat through stress-induced activation.

[The functional analysis of polymorphic insertions of Alu retroelements at acute lymphoblastic leukemia]

Human genome variability observed within patient cohorts is considered as a goal of functional genomics essential for personalized medicine progress. In the current research we implement functional analysis of 31 polymorphic Alu insertions located within gene introns for individual genomes of patients with acute lymphoblastic leukemia (ALL). As a result we demonstrated a decrease of the primary transcripts content for 21 Alu-containing alleles. The most strong inhibitory effect of 10 Alu insertions was observed in both mononuclear blood cells of healthy donors and B-lymphoblasts of ALL patients. Allele frequencies of three Alu insertions that are located in MEF2C (two of them) and TAX1BP1 genes significantly differ (p-value 0.027. 0.052. 0.014 accordingly) between cohorts of healthy donors and ALL patients. Prolong influence of the Alu insertions on intracellular content of mature mRNA was studied for corresponding allele of TARBP1 gene.

Why share data? Lessons learned from the fMRIDC

Neuroimaging and the discipline of cognitive neuroscience have grown together in lock-step with each pushing the other toward an improved ability to explore and examine brain function and form. However successful neuroimaging and the examination of cognitive processes may seem today, the culture of data sharing in these fields remains underdeveloped. In this article, we discuss our own experience in the development of the fMRI Data Center (fMRIDC) - a large-scale effort to gather, curate, and openly share the complete data sets from published research articles of brain activation studies using fMRI. We outline the fMRIDC effort's beginnings, how it operated, note some of the sociological reactions we received, and provide several examples of prominent new studies performed using data drawn from the archive. Finally, we provide comment on what considerations are needed for successful neuroimaging databasing and data sharing as existing and emerging efforts take the next steps in archiving and disseminating the field's valuable and irreplaceable data.

Development of multilocus putatively neutral DNA markers in the X-chromosome for population genetic studies in humans

BACKGROUND

It has now been well documented that the type (coding, non-coding) and location (nuclear, mitochondrial etc.) of genetic markers heavily influence evolutionary inferences; realistic assumptions can be drawn if multiple putatively neutral DNA fragments spread across the genome are used.

AIM

To infer human population history, Single Nucleotide Polymorphisms (SNPs), located in the non-coding regions of different genes in the X-chromosome have been developed as 'putatively neutral markers'.

SUBJECTS AND METHODS

A population sample consisting of 16 male individuals from the western part of India was utilized for sequencing eight DNA fragments located in introns of three genes (Duchenne muscular dystrophy, Factor IX and Pyruvate dehydrogenase E1 sub-unit) on the human X-chromosome. PCR amplification and DNA sequencing confirmed the polymorphic status of all the fragments.

RESULTS

Twenty nine SNPs were found to be segregating in the Western Indian population samples. Using these SNPs the nucleotide diversity and demographic parameters of the Western Indian population were estimated. Several tests of neutrality ascertained that all eight fragments evolve putatively neutrally. Further, linkage disequilibrium analyses confirmed this fact.

CONCLUSION

All eight DNA fragments seem to bear the characteristics to be considered as 'putatively neutral genetic markers' and thus, could be utilized for inference of human population and demographic histories.

Marker utility of miniature inverted-repeat transposable elements for wheat biodiversity and evolution

Transposable elements (TEs) account for up to 80% of the wheat genome and are considered one of the main drivers of wheat genome evolution. However, the contribution of TEs to the divergence and evolution of wheat genomes is not fully understood. In this study, we have developed 55 miniature inverted-repeat transposable element (MITE) markers that are based on the presence/absence of an element, with over 60% of these 55 MITE insertions associated with wheat genes. We then applied these markers to assess genetic diversity among Triticum and Aegilops species, including diploid (AA, BB and DD genomes), tetraploid (BBAA genome) and hexaploid (BBAADD genome) species. While 18.2% of the MITE markers showed similar insertions in all species indicating that those are fossil insertions, 81.8% of the markers showed polymorphic insertions among species, subspecies, and accessions. Furthermore, a phylogenetic analysis based on MITE markers revealed that species were clustered based on genus, genome composition, and ploidy level, while 47.13% genetic divergence was observed between the two main clusters, diploids versus polyploids. In addition, we provide evidence for MITE dynamics in wild emmer populations. The use of MITEs as evolutionary markers might shed more light on the origin of the B-genome of polyploid wheat.

The human genome: a multifractal analysis

BACKGROUND

Several studies have shown that genomes can be studied via a multifractal formalism. Recently, we used a multifractal approach to study the genetic information content of the Caenorhabditis elegans genome. Here we investigate the possibility that the human genome shows a similar behavior to that observed in the nematode.

RESULTS

We report here multifractality in the human genome sequence. This behavior correlates strongly on the presence of Alu elements and to a lesser extent on CpG islands and (G+C) content. In contrast, no or low relationship was found for LINE, MIR, MER, LTRs elements and DNA regions poor in genetic information. Gene function, cluster of orthologous genes, metabolic pathways, and exons tended to increase their frequencies with ranges of multifractality and large gene families were located in genomic regions with varied multifractality. Additionally, a multifractal map and classification for human chromosomes are proposed.

CONCLUSIONS

Based on these findings, we propose a descriptive non-linear model for the structure of the human genome, with some biological implications. This model reveals 1) a multifractal regionalization where many regions coexist that are far from equilibrium and 2) this non-linear organization has significant molecular and medical genetic implications for understanding the role of Alu elements in genome stability and structure of the human genome. Given the role of Alu sequences in gene regulation, genetic diseases, human genetic diversity, adaptation and phylogenetic analyses, these quantifications are especially useful.

[Retroposons in modern human genome evolution]

The ascertainment of the rates and driving forces of human genome evolution along with the genetic diversity of populations or separate population groups remains a topical problem of fundamental and applied genomics. According to the results of comparative analysis, the most numerous human genome structure peculiarities are connected with the distribution of mobile genetic retroelements - LTR, LINE1, SVA, and Alu repeats. Due to the wide distribution in different genome loci, conversed retropositional activity, and the retroelements regulatory potential, let us regard them as one of the significant evolutionary driving forces and the source of human genome variability. In the current review, we summarize published data and recent results of our research aimed at the analysis of the evolutionary impact of the young retroelements group on the function and variability of the human genome. We examine modern approaches of the polygenomic identification of polymorphic retroelements inserts. Using an original Internet resource, we analyze special features of the genomic polymorphic inserts of Alu repeats. We thoroughly characterize the strategy of large-scale functional analysis of polymorphic retroelement inserts. The presented results confirm the hypothesis of the roles of retroelements as active cis regulatory elements that are able to modulate surrounding genes.

Laboratory methods for the analysis of primate mobile elements

Mobile elements represent a unique and powerful set of tools for understanding the variation in a genome. Methods exist not only to utilize the polymorphisms among and within taxa to various ends but also to investigate the mechanism through which mobilization occurs. The number of methods to accomplish these ends is ever growing. Here, we present several protocols designed to assay mobile element-based variation within and among individual genomes.

Detecting HIV-1 integration by repetitive-sampling Alu-gag PCR

In this review, we compare four assays that are currently used to measure HIV integration and discuss their strengths and weaknesses. We then outline advances that have been made toward development of a more robust, more sensitive, quantitative HIV integration assay suitable for clinical use. The assay that we have developed uses repetitive-sampling Alu-gag PCR. The detailed protocol describes our assay step-by-step, the creation of an integration standard cell line and accompanying standard curve, as well as the quantitation of integration and calculation of associated error estimates. Finally, we speculate on fundamental, unresolved issues in HIV latency that can be addressed by measuring HIV integration.

Analysis of newly identified low copy AluYj subfamily

Human specific AluY elements were investigated by comparative analysis between human chromosome 21 and chimpanzee chromosome 22. Human specific AluY element was identified on human chromosome 21q22 (accession no. AL163282), and then that was a new member of AluYj subfamily. From the bioinformatic analysis, AluYj subfamily was investigated in human whole genome using AluYj4 consensus sequence (accession no. AL163282). Thirteen members of the AluYj4 elements (4 diagnostic mutations) and eight members of the AluYj3 elements (3 diagnostic mutations) were identified with distinct diagnostic mutation from AluY consensus sequence. The results of the molecular clock calculation of non-CpG region substitution indicated that, AluYj4 elements (2.1 million years old) may be proliferated more recent time than AluYj3 elements (14.1 million years old). For the verification of recent insertion time, four of AluYj4 elements (ch2-AC017101, ch10-AC044786, ch12-AC007656 and ch21-AL163282) from human chromosomes 2, 10, 12, 21 were analyzed by PCR amplification using various human and primate DNA samples. Though, no polymorphism was detected in human population, we identified the new AluYj4 subfamily as the human specific elements.

Similarity searching using BLAST

Similarity searches are an essential component of most bioinformatic applications. They form the bases of structural motif identification, gene identification, and insights into functional associations. With the rapid increase in the available genetic data through a wide variety of databases, similarity searches are an essential tool for accessing these data in an informative and productive way. In this chapter, we provide an overview of similarity searching approaches, related databases, and parameter options to achieve the best results for a variety of applications. We then provide a worked example and some notes for consideration.

Systems approaches to preventing transplanted cell death in cardiac repair

Stem cell transplantation may repair the injured heart, but tissue regeneration is limited by death of transplanted cells. Most cell death occurs in the first few days post-transplantation, likely from a combination of ischemia, anoikis and inflammation. Interventions known to enhance transplanted cell survival include heat shock, over-expressing anti-apoptotic proteins, free radical scavengers, anti-inflammatory therapy and co-delivery of extracellular matrix molecules. Combinatorial use of such interventions markedly enhances graft cell survival, but death still remains a significant problem. We review these challenges to cardiac cell transplantation and present an approach to systematically address them. Most anti-death studies use histology to assess engraftment, which is time- and labor-intensive. To increase throughput, we developed two biochemical approaches to follow graft viability in the mouse heart. The first relies on LacZ enzymatic activity to track genetically modified cells, and the second quantifies human genomic DNA content using repetitive Alu sequences. Both show linear relationships between input cell number and biochemical signal, but require correction for the time lag between cell death and loss of signal. Once optimized, they permit detection of as few as 1 graft cell in 40,000 host cells. Pro-survival effects measured biochemically at three days predict long-term histological engraftment benefits. These methods permitted identification of carbamylated erythropoietin (CEPO) as a pro-survival factor for human embryonic stem cell-derived cardiomyocyte grafts. CEPO's effects were additive to heat shock, implying independent survival pathways. This system should permit combinatorial approaches to enhance graft viability in a fraction of the time required for conventional histology.

Mobile DNA elements in primate and human evolution

Roughly 50% of the primate genome consists of mobile, repetitive DNA sequences such as Alu and LINE1 elements. The causes and evolutionary consequences of mobile element insertion, which have received considerable attention during the past decade, are reviewed in this article. Because of their unique mutational mechanisms, these elements are highly useful for answering phylogenetic questions. We demonstrate how they have been used to help resolve a number of questions in primate phylogeny, including the human-chimpanzee-gorilla trichotomy and New World primate phylogeny. Alu and LINE1 element insertion polymorphisms have also been analyzed in human populations to test hypotheses about human evolution and population affinities and to address forensic issues. Finally, these elements have had impacts on the genome itself. We review how they have influenced fundamental ongoing processes like nonhomologous recombination, genomic deletion, and X chromosome inactivation.

How many populations set foot through the Patagonian door? Genetic composition of the current population of Bahía Blanca (Argentina) based on data from 19 Alu polymorphisms

The city of Bahía Blanca occupies a strategic place in Argentina south of the Pampean region in the north-east corner of the Patagonia. Since 1828, this city has been the historical and political border between Amerindian lands in the south, and the lands of European colonists. Nowadays, Bahía Blanca is an urban population mainly composed by descendents of immigrants from Spain and other European countries with apparently low admixture with Amerindians. In view of the unexpectedly high Amerindian admixture levels (about 46.7%) suggested by mtDNA data, and protein markers (19.5%), we analyzed a set of 19 Alu polymorphisms (18 autosomal, 1 of Chromosome Y) in a well-documented genealogical sample from Bahía Blanca. The genotyped sample was made up of 119 unrelated healthy individuals whose birth place and grandparent origins were fully documented. According to available genealogical records, the total sample has been subdivided into two groups: Bahía Blanca Original (64 individuals with all 4 gandparents born in Argentina) and Bahía Blanca Mix (55 individuals with one to three grandparents born out of Argentina). Allele frequencies and gene diversity values in Bahía Blanca fit well into the European ranges. Population relationships have been tested for 8 Alu markers, whose variation has been described in several Amerindian and European samples. Reynolds genetic distances underline the significant genetic similarity of Bahía Blanca to Europeans (mean distance 0.044) and their differentiation from Amerindians (0.146). Interestingly enough, when the general sample is divided, Bahía Blanca Original appears slightly closer to Amerindians (0.127) in contrast to Bahía Blanca Mix (0.161). Furthermore, the genetic relationships depicted through a principal components analysis emphasize the relative similarity of Bahía Blanca Original to Amerindians. A thorough knowledge of the sample origins has allowed us to make a subtle distinction of the genetic composition of Bahía Blanca.

Web-based analysis for student-generated complex genetic profiles

A simple, rapid method for generating complex genetic profiles using Alu-based markers was recently developed for students primarily at the undergraduate level to learn more about forensics and paternity analysis. On the basis of the Cold Spring Harbor Allele Server, which provides an excellent tool for analyzing a single Alu variant, we present a new web-based system for analyzing several genetic loci, including Hardy-Weinberg equilibrium, genetic drift, and Fst genetic-distance calculations, as well as analyzing eight loci profiles simultaneously for forensic purposes. By analyzing several loci, students can determine more precisely the relatedness of populations as well as develop a greater appreciation for the use of DNA markers in forensic analysis by concurrently assessing the frequencies of genotypes for eight genetic loci.

SINEs of speciation: tracking lineages with retroposons

The value of short interspersed elements (SINEs) for diagnosing common ancestry is being expanded to examine the differential sorting of lineages through the course of speciation events. Because most SINEs are neutral markers of identical descent, are not precisely excised from the genome and have a known ancestral condition, they are advantageous for reconciling gene trees and species trees with minimal phylogenetic error. A population perspective on SINE evolution combined with coalescence theory provides a context for investigating the phenomenon of ancestral polymorphism and its role in producing incongruent SINE insertion patterns among multiple loci. Studies of human Alu repeats demonstrate the value of young polymorphic SINEs for assessing human genomic diversity and tracking ancient demographics of human populations, whereas incongruent insertion patterns revealed by older fixed SINE loci, such as those in African cichlid fishes, contain information that might help identify ancient radiations that are otherwise obscured by accumulated mutations in sequence data. Here, we review the utility of retroposons for inferring common ancestry, discuss limits to the method, and clarify confusion by providing examples from the literature that illustrate how discordant multi-locus insertion patterns of retroelements can indicate lineage-sorting events that should not be misinterpreted as phylogenetic noise.

Hypomorphic mutations in the gene encoding a key Fanconi anemia protein, FANCD2, sustain a significant group of FA-D2 patients with severe phenotype

FANCD2 is an evolutionarily conserved Fanconi anemia (FA) gene that plays a key role in DNA double-strand-type damage responses. Using complementation assays and immunoblotting, a consortium of American and European groups assigned 29 patients with FA from 23 families and 4 additional unrelated patients to complementation group FA-D2. This amounts to 3%-6% of FA-affected patients registered in various data sets. Malformations are frequent in FA-D2 patients, and hematological manifestations appear earlier and progress more rapidly when compared with all other patients combined (FA-non-D2) in the International Fanconi Anemia Registry. FANCD2 is flanked by two pseudogenes. Mutation analysis revealed the expected total of 66 mutated alleles, 34 of which result in aberrant splicing patterns. Many mutations are recurrent and have ethnic associations and shared allelic haplotypes. There were no biallelic null mutations; residual FANCD2 protein of both isotypes was observed in all available patient cell lines. These analyses suggest that, unlike the knockout mouse model, total absence of FANCD2 does not exist in FA-D2 patients, because of constraints on viable combinations of FANCD2 mutations. Although hypomorphic mutations arie involved, clinically, these patients have a relatively severe form of FA.

Mobile element-based forensic genomics

Mobile elements are commonly referred to as selfish repetitive DNA sequences. However, mobile elements represent a unique and underutilized group of molecular markers. Several of their characteristics make them ideally suited for use as tools in forensic genomic applications. These include their nature as essentially homoplasy-free characters, they are identical by descent, the ancestral state of any insertion is known to be the absence of the element, and many mobile element insertions are lineage specific. In this review, we provide an overview of mobile element biology and describe the application of certain mobile elements, especially the SINEs and other retrotransposons, to forensic genomics. These tools include quantitative species-specific DNA detection, analysis of complex biomaterials, and the inference of geographic origin of human DNA samples.

Identification and characterization of novel polymorphic LINE-1 insertions through comparison of two human genome sequence assemblies

Mobile elements represent a relatively new class of markers for the study of human evolution. Long interspersed elements (LINEs) belong to a group of retrotransposons comprising approximately 21% of the human genome. Young LINE-1 (L1) elements that have integrated recently into the human genome can be polymorphic for insertion presence/absence in different human populations at particular chromosomal locations. To identify putative novel L1 insertion polymorphisms, we computationally compared two draft assemblies of the whole human genome (Public and Celera Human Genome assemblies). We identified a total of 148 potential polymorphic L1 insertion loci, among which 73 were candidates for novel polymorphic loci. Based on additional analyses we selected 34 loci for further experimental studies. PCR-based assays and DNA sequence analysis were performed for these 34 loci in 80 unrelated individuals from four diverse human populations: African-American, Asian, Caucasian, and South American. All but two of the selected loci were confirmed as polymorphic in our human population panel. Approximately 47% of the analyzed loci integrated into other repetitive elements, most commonly older L1s. One of the insertions was accompanied by a BC200 sequence. Collectively, these mobile elements represent a valuable source of genomic polymorphism for the study of human population genetics. Our results also suggest that the exhaustive identification of L1 insertion polymorphisms is far from complete, and new whole genome sequences are valuable sources for finding novel retrotransposon insertion polymorphisms.

Identification of a unique Alu-based polymorphism and its use in human population studies

Alu elements represent a family of short interspersed DNA elements (SINEs) found in primate genomes. These are members of a group of transposable elements that integrate into the genome by the process of retrotransposition. Recent integrations of Alu elements within the human genome have generated presence/absence variants useful as DNA markers in human population studies as well as in forensic and paternity analyses. Besides the ease of use, this type of marker is unique because the absence of the Alu represents the ancestral form. We have identified an Alu-based polymorphism that consists of four alleles in which we can predict the evolutionary order. Additionally, we have developed a simple PCR plus restriction endonuclease assay to readily distinguish the four alleles. We have thus far analyzed DNA from a small set of samples comprising ten different ethnic groups. The three populations of African descent exhibited a relatively low frequency of the absence allele in contrast to the other populations, as well as being the only populations in which all four alleles were identified. One presence allele was not found in both European Caucasian and South American populations that were sampled, whereas a different presence allele was not observed among the sampled Asian populations. Additionally, the four-allele system identified variations among populations not observed by simply scoring as presence/absence variants. Therefore, extending beyond the two-allele dimorphic Alu system further elucidates population variations. These features afford this marker as a unique tool in the study of both global and regional analyses of human populations.

Recently integrated Alu retrotransposons are essentially neutral residents of the human genome

Alu elements represent the largest family of human mobile elements in copy number. A controversial issue with implications for both Alu biology and human genome evolution is whether selective pressures are affecting Alu elements on a large scale. To address this issue, we analyzed the genomic distribution of the three youngest known human Alu subfamilies (Ya5a2, Ya8 and Yb9) in conjunction with their insertion polymorphism status in the human population, since selection can only act on polymorphic elements. Our results indicate that: (i) polymorphic and fixed recently integrated Alu elements are found in genomic regions whose GC contents are statistically indistinguishable, and (ii) recently integrated Alu elements are inserted randomly, regardless of the GC content of the surrounding genomic DNA. These results provide strong evidence that recently integrated "young" Alu elements are not subject to positive or negative selection on a large scale. Therefore, young Alu elements can be regarded as essentially neutral residents of the human genome. These results also imply that selective processes specifically targeting Alu elements can be ruled out as explanations for the accumulation of Alu elements in GC-rich regions of the human genome.

Whole genome computational comparative genomics: A fruitful approach for ascertaining Alu insertion polymorphisms

Alu elements are the most active and predominant type of short interspersed elements (SINEs) in the human genome. Recently inserted polymorphic (for presence/absence) Alu elements contribute to genome diversity among different human populations, and they are useful genetic markers for population genetic studies. The objective of this study is to identify polymorphic Alu insertions through an in silico comparative genomics approach and to analyze their distribution pattern throughout the human genome. By computationally comparing the public and Celera sequence assemblies of the human genome, we identified a total of 800 polymorphic Alu elements. We used polymerase chain reaction-based assays to screen a randomly selected set of 16 of these 800 Alu insertion polymorphisms using a human diversity panel to demonstrate the efficiency of our approach. Based on sequence analysis of the 800 Alu polymorphisms, we report three new Alu subfamilies, Ya3, Ya4b, and Yb11, with Yb11 being the smallest known Alu subfamily. Analysis of retrotransposition activity revealed Yb11, Ya8, Ya5, Yb9, and Yb8 as the most active Alu subfamilies and the maintenance of a very low level of retrotransposition activity or recent gene conversion events involving S subfamilies. The 800 polymorphic Alu insertions are characterized by the presence of target site duplications (TSDs) and longer than average polyA-tail length. Their pre-integration sites largely follow an extended "NT-AARA" motif. Among chromosomes, the density of Alu insertion polymorphisms is positively correlated with the Alu-site availability and is inversely correlated with the densities of older Alu elements and genes.