Genetic and epigenetic variations contributed by Alu retrotransposition

<p>Abstract</p> <p>Background</p> <p><it>De novo </it>retrotransposition of Alu elements has been recognized as a major driver for insertion polymorphisms in human populations. In this study, we exploited Alu-anchored bisulfite PCR libraries to identify evolutionarily recent Alu element insertions, and to investigate their genetic and epigenetic variation.</p> <p>Results</p> <p>A total of 327 putatively recent Alu insertions were identified, altogether represented by 1,762 sequence reads. Nearly all such <it>de novo </it>retrotransposition events (316/327) were novel. Forty-seven out of forty-nine randomly selected events, corresponding to nineteen genomic loci, were sequence-verified. Alu element insertions remained hemizygous in one or more individuals in sixteen of the nineteen genomic loci. The Alu elements were found to be enriched for young Alu families with characteristic sequence features, such as the presence of a longer poly(A) tail. In addition, we documented the occurrence of a duplication of the AT-rich target site in their immediate flanking sequences, a hallmark of retrotransposition. Furthermore, we found the sequence motif (TT/AAAA) that is recognized by the ORF2P protein encoded by LINE-1 in their 5'-flanking regions, consistent with the fact that Alu retrotransposition is facilitated by LINE-1 elements. While most of these Alu elements were heavily methylated, we identified an Alu localized 1.5 kb downstream of TOMM5 that exhibited a completely unmethylated left arm. Interestingly, we observed differential methylation of its immediate 5' and 3' flanking CpG dinucleotides, in concordance with the unmethylated and methylated statuses of its internal 5' and 3' sequences, respectively. Importantly, TOMM5's CpG island and the 3 Alu repeats and 1 MIR element localized upstream of this newly inserted Alu were also found to be unmethylated. Methylation analyses of two additional genomic loci revealed no methylation differences in CpG dinucleotides flanking the Alu insertion sites in the two homologous chromosomes, irrespective of the presence or absence of the insertion.</p> <p>Conclusions</p> <p>We anticipate that the combination of methodologies utilized in this study, which included repeat-anchored bisulfite PCR sequencing and the computational analysis pipeline herein reported, will prove invaluable for the generation of genetic and epigenetic variation maps.</p

Additions to the Mycosphaerella complex

Species in the present study were compared based on their morphology, growth characteristics in culture, and DNA sequences of the nuclear ribosomal RNA gene operon (including ITS1, ITS2, 5.8S nrDNA and the first 900 bp of the 28S nrDNA) for all species and partial actin and translation elongation factor 1-alpha gene sequences for Cladosporium species. New species of Mycosphaerella (Mycosphaerellaceae) introduced in this study include M. cerastiicola (on Cerastium semidecandrum, The Netherlands), and M. etlingerae (on Etlingera elatior, Hawaii). Mycosphaerella holualoana is newly reported on Hedychium coronarium (Hawaii). Epitypes are also designated for Hendersonia persooniae, the basionym of Camarosporula persooniae, and for Sphaerella agapanthi, the basionym of Teratosphaeria agapanthi comb. nov. (Teratosphaeriaceae) on Agapathus umbellatus from South Africa. The latter pathogen is also newly recorded from A. umbellatus in Europe (Portugal). Furthermore, two sexual species of Cladosporium (Davidiellaceae) are described, namely C. grevilleae (on Grevillea sp., Australia), and C. silenes (on Silene maritima, UK). Finally, the phylogenetic position of two genera are newly confirmed, namely Camarosporula (based on C. persooniae, teleomorph Anthracostroma persooniae), which is a leaf pathogen of Persoonia spp. in Australia, belongs to the Teratosphaeriaceae, and Sphaerulina (based on S. myriadea), which occurs on leaves of Fagaceae (Carpinus, Castanopsis, Fagus, Quercus), and belongs to the Mycosphaerellaceae

Identification of Widespread Ultra-Edited Human RNAs

Adenosine-to-inosine modification of RNA molecules (A-to-I RNA editing) is an important mechanism that increases transciptome diversity. It occurs when a genomically encoded adenosine (A) is converted to an inosine (I) by ADAR proteins. Sequencing reactions read inosine as guanosine (G); therefore, current methods to detect A-to-I editing sites align RNA sequences to their corresponding DNA regions and identify A-to-G mismatches. However, such methods perform poorly on RNAs that underwent extensive editing (“ultra”-editing), as the large number of mismatches obscures the genomic origin of these RNAs. Therefore, only a few anecdotal ultra-edited RNAs have been discovered so far. Here we introduce and apply a novel computational method to identify ultra-edited RNAs. We detected 760 ESTs containing 15,646 editing sites (more than 20 sites per EST, on average), of which 13,668 are novel. Ultra-edited RNAs exhibit the known sequence motif of ADARs and tend to localize in sense strand Alu elements. Compared to sites of mild editing, ultra-editing occurs primarily in Alu-rich regions, where potential base pairing with neighboring, inverted Alus creates particularly long double-stranded RNA structures. Ultra-editing sites are underrepresented in old Alu subfamilies, tend to be non-conserved, and avoid exons, suggesting that ultra-editing is usually deleterious. A possible biological function of ultra-editing could be mediated by non-canonical splicing and cleavage of the RNA near the editing sites

Does Selection against Transcriptional Interference Shape Retroelement-Free Regions in Mammalian Genomes?

BACKGROUND: Eukaryotic genomes are scattered with retroelements that proliferate through retrotransposition. Although retroelements make up around 40 percent of the human genome, large regions are found to be completely devoid of retroelements. This has been hypothesised to be a result of genomic regions being intolerant to insertions of retroelements. The inadvertent transcriptional activity of retroelements may affect neighbouring genes, which in turn could be detrimental to an organism. We speculate that such retroelement transcription, or transcriptional interference, is a contributing factor in generating and maintaining retroelement-free regions in the human genome. METHODOLOGY/PRINCIPAL FINDINGS: Based on the known transcriptional properties of retroelements, we expect long interspersed elements (LINEs) to be able to display a high degree of transcriptional interference. In contrast, we expect short interspersed elements (SINEs) to display very low levels of transcriptional interference. We find that genomic regions devoid of long interspersed elements (LINEs) are enriched for protein-coding genes, but that this is not the case for regions devoid of short interspersed elements (SINEs). This is expected if genes are subject to selection against transcriptional interference. We do not find microRNAs to be associated with genomic regions devoid of either SINEs or LINEs. We further observe an increased relative activity of genes overlapping LINE-free regions during early embryogenesis, where activity of LINEs has been identified previously. CONCLUSIONS/SIGNIFICANCE: Our observations are consistent with the notion that selection against transcriptional interference has contributed to the maintenance and/or generation of retroelement-free regions in the human genome

Involvement of JNK-mediated pathway in EGF-mediated protection against paclitaxel-induced apoptosis in SiHa human cervical cancer cells

We investigated the signalling pathways by which epidermal growth factor (EGF) modulates paclitaxel-induced apoptosis in SiHa human cervical cancer cells. SiHa cells exposed to paclitaxel underwent apoptosis, which was strongly inhibited by EGF. This inhibition of apoptosis by EGF was not altered by pharmacological blockade of phosphatidylinositol 3′-OH kinase (PI-3K) with the PI-3K specific inhibitor LY294002 or blockade of the mitogen-activated protein kinase (MAPK) kinase (MEK) with the MEK specific inhibitor PD98059, or by transfection of the cells with PI-3K or MEK dominant-negative expression vectors. EGF did not stimulate PI-3K/Akt, MEK/MAPK, or p38 MAPK activity in SiHa cells but did transiently activate the c-Jun NH2-terminal kinase (JNK). Co-exposure of SiHa cells to SB202190 at concentrations that inhibit JNK abolished the protective effect of EGF on SiHa cells against paclitaxel-induced apoptosis. Our findings indicate that the JNK signaling pathway plays an important role in EGF-mediated protection from paclitaxel-induced apoptosis in SiHa cells. © 2001 Cancer Research Campaign http://www.bjcancer.co

Positive Selection within the Schizophrenia-Associated GABA(A) Receptor β(2) Gene

The gamma-aminobutyric acid type-A (GABA(A)) receptor plays a major role in inhibitory neurotransmissions. Intronic SNPs and haplotypes in GABRB2, the gene for GABA(A) receptor β(2) subunit, are associated with schizophrenia and correlated with the expression of two alternatively spliced β(2) isoforms. In the present study, using chimpanzee as an ancestral reference, high frequencies were observed for the derived (D) alleles of the four SNPs rs6556547, rs187269, rs1816071 and rs1816072 in GABRB2, suggesting the occurrence of positive selection for these derived alleles. Coalescence-based simulation showed that the population frequency spectra and the frequencies of H56, the haplotype having all four D alleles, significantly deviated from neutral-evolution expectation in various demographic models. Haplotypes containing the derived allele of rs1816072 displayed significantly less diversity compared to haplotypes containing its ancestral allele, further supporting positive selection. The variations in DD-genotype frequencies in five human populations provided a snapshot of the evolutionary history, which suggested that the positive selections of the D alleles are recent and likely ongoing. The divergence between the DD-genotype profiles of schizophrenic and control samples pointed to the schizophrenia-relevance of positive selections, with the schizophrenic samples showing weakened selections compared to the controls. These DD-genotypes were previously found to increase the expression of β(2), especially its long isoform. Electrophysiological analysis showed that this long β(2) isoform favored by the positive selections is more sensitive than the short isoform to the inhibition of GABA(A) receptor function by energy depletion. These findings represent the first demonstration of positive selection in a schizophrenia-associated gene

Multi-Locus Phylogeographic and Population Genetic Analysis of Anolis carolinensis: Historical Demography of a Genomic Model Species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N = 62 to 152) and one mitochondrial locus (N = 226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, Bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ∼2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range

High Differentiation among Eight Villages in a Secluded Area of Sardinia Revealed by Genome-Wide High Density SNPs Analysis

To better design association studies for complex traits in isolated populations it's important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not “a priori” be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of sub-populations and stratification within genetic isolates

Nuclear versus mitochondrial DNA: evidence for hybridization in colobine monkeys

<p>Abstract</p> <p>Background</p> <p>Colobine monkeys constitute a diverse group of primates with major radiations in Africa and Asia. However, phylogenetic relationships among genera are under debate, and recent molecular studies with incomplete taxon-sampling revealed discordant gene trees. To solve the evolutionary history of colobine genera and to determine causes for possible gene tree incongruences, we combined presence/absence analysis of mobile elements with autosomal, X chromosomal, Y chromosomal and mitochondrial sequence data from all recognized colobine genera.</p> <p>Results</p> <p>Gene tree topologies and divergence age estimates derived from different markers were similar, but differed in placing <it>Piliocolobus/Procolobus </it>and langur genera among colobines. Although insufficient data, homoplasy and incomplete lineage sorting might all have contributed to the discordance among gene trees, hybridization is favored as the main cause of the observed discordance. We propose that African colobines are paraphyletic, but might later have experienced female introgression from <it>Piliocolobus</it>/<it>Procolobus </it>into <it>Colobus</it>. In the late Miocene, colobines invaded Eurasia and diversified into several lineages. Among Asian colobines, <it>Semnopithecus </it>diverged first, indicating langur paraphyly. However, unidirectional gene flow from <it>Semnopithecus </it>into <it>Trachypithecus </it>via male introgression followed by nuclear swamping might have occurred until the earliest Pleistocene.</p> <p>Conclusions</p> <p>Overall, our study provides the most comprehensive view on colobine evolution to date and emphasizes that analyses of various molecular markers, such as mobile elements and sequence data from multiple loci, are crucial to better understand evolutionary relationships and to trace hybridization events. Our results also suggest that sex-specific dispersal patterns, promoted by a respective social organization of the species involved, can result in different hybridization scenarios.</p

Nuclear Receptor HNF4α Binding Sequences are Widespread in Alu Repeats

<p>Abstract</p> <p>Background</p> <p>Alu repeats, which account for ~10% of the human genome, were originally considered to be junk DNA. Recent studies, however, suggest that they may contain transcription factor binding sites and hence possibly play a role in regulating gene expression.</p> <p>Results</p> <p>Here, we show that binding sites for a highly conserved member of the nuclear receptor superfamily of ligand-dependent transcription factors, hepatocyte nuclear factor 4alpha (HNF4α, NR2A1), are highly prevalent in Alu repeats. We employ high throughput protein binding microarrays (PBMs) to show that HNF4α binds > 66 unique sequences in Alu repeats that are present in ~1.2 million locations in the human genome. We use chromatin immunoprecipitation (ChIP) to demonstrate that HNF4α binds Alu elements in the promoters of target genes (<it>ABCC3, APOA4, APOM, ATPIF1, CANX, FEMT1A, GSTM4, IL32, IP6K2, PRLR, PRODH2, SOCS2, TTR</it>) and luciferase assays to show that at least some of those Alu elements can modulate HNF4α-mediated transactivation <it>in vivo </it>(<it>APOM, PRODH2, TTR, APOA4</it>). HNF4α-Alu elements are enriched in promoters of genes involved in RNA processing and a sizeable fraction are in regions of accessible chromatin. Comparative genomics analysis suggests that there may have been a gain in HNF4α binding sites in Alu elements during evolution and that non Alu repeats, such as Tiggers, also contain HNF4α sites.</p> <p>Conclusions</p> <p>Our findings suggest that HNF4α, in addition to regulating gene expression via high affinity binding sites, may also modulate transcription via low affinity sites in Alu repeats.</p