Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.

     

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.     

The strength of selection on ultraconserved elements in the human genome

Ultraconserved elements are stretches of consecutive nucleotides that are perfectly conserved in multiple mammalian genomes. Although these sequences are identical in the reference human, mouse, and rat genomes, we identified numerous polymorphisms within these regions in the human population. To determine whether polymorphisms in ultraconserved elements affect fitness, we genotyped unrelated human DNA samples at loci within these sequences. For all single-nucleotide polymorphisms tested in ultraconserved regions, individuals homozygous for derived alleles (alleles that differ from the rodent reference genomes) were present, viable, and healthy. The distribution of allele frequencies in these samples argues against strong, ongoing selection as the force maintaining the conservation of these sequences. We then used two methods to determine the minimum level of selection required to generate these sequences. Despite the lack of fixed differences in these sequences between humans and rodents, the average level of selection on ultraconserved elements is less than that on essential genes. The strength of selection associated with ultraconserved elements suggests that mutations in these regions may have subtle phenotypic consequences that are not easily detected in the laboratory.     

Analysis of invariant sequences in 266 complete genomes

To date, the complete genome sequences of more than 250 organisms have been determined. This information can now be used to determine whether there exist any invariant sequences that are conserved among all organisms, from bacteria to plants, animals, and humans. The existence of invariant sequences would strongly suggest that these sequences have been inherited unchanged from the last common ancestor of all life, and that they have essential functions. We have developed a new software program to identify invariant sequences conserved among the currently sequenced genomes and applied this analysis to the complete genome sequences of 266 organisms. We have identified 3 invariant DNA sequences longer than or equal to 11 bp and 6 invariant amino acid sequences longer than or equal to 6 aa. The longest invariant DNA sequence, AAGTCGTACAAGGT (15 bp), was found in the 16S/18S rRNA gene. Two 8 aa sequences, GHVDHGKT in IF2 and EF-Tu and DTPGHVDF in EF-G, were the longest invariant amino acid sequences detected. These sequences could be essential elements from the genome of the last common ancestor and may have remained unchanged throughout evolution.     

S Elements: A Family of Tc1-like Transposons in the Genome of Drosophila Melanogaster

The S elements form a diverse family of long-inverted-repeat transposons within the genome of Drosophila melanogaster. These elements vary in size and sequence, the longest consisting of 1736 bp with 234-bp inverted terminal repeats. The longest open reading frame in an intact S element could encode a 345-amino acid polypeptide. This polypeptide is homologous to the transposases of the mariner-Tc1 superfamily of transposable elements. S elements are ubiquitous in D. melanogaster populations and also appear to be present in the genomes of two sibling species; however, they seem to be absent from 17 other Drosophila species that were examined. Within D. melanogaster strains, there are, on average, 37.4 cytologically detectable S elements per diploid genome. These elements are scattered throughout the chromosomes, but several sites in both the euchromatin and β heterochromatin are consistently occupied. The discovery of an S-element-insertion mutation and a reversion of this mutation indicates that S elements are at least occasionally mobile in the D. melanogaster genome. These elements seem to insert at an AT dinucleotide within a short palindrome and apparently duplicate that dinucleotide upon insertion.     

Diversification of the rice Waxy gene by insertion of mobile DNA elements into introns.

The waxy (wx) gene of Oryza glaberrima was cloned, and its nucleotide sequence was determined. A waxy mutant of O. glaberrima showing a glutinous phenotype was found to contain a substitution mutation generating a termination codon in the coding region of the wx gene. The Wx sequence of O. glaberrima was different from that of Oryza sativa by substitutions and insertions/deletions, among which only a few substitutions occurred in several exons not to severely alter the amino acid sequence of the Wx protein. The most striking difference observed in introns was a 139-bp deletion (or insertion) in intron 10 of O. glaberrima (or O. sativa). In O. sativa, 125 bp of the 139-bp sequence was flanked by direct repeats of a 14-bp sequence. A sequence homologous to the 125-bp sequence was found in the region preceding exon 2; this sequence was also flanked by direct repeats of another 14-bp sequence. This result and the observation that the 125-bp sequence was interspersed in rice genomes indicate that they are SINEs (short interspersed elements) in the plant system. We also identified a DNA sequence with long terminal inverted repeats in intron 13 of both O. glaberrima and O. sativa. This sequence was present in multiple copies in rice genomes, suggesting that it is a transposable element. These results obtained suggest that mobile DNA elements have diversified the rice Waxy gene by inserting into introns, each of which may originally have a length of about 100 bp.     

The Longest (A+T) and (G+C) Blocks in the Human and Other Genomes

Abstract

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16^th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are embeeded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

The longest (A+T) and (G+C) blocks in the human and other genomes

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are imbedded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

Widespread occurrence of small inversions in the chloroplast genomes of land plants

Large inversions are well characterized in the chloroplast genomes of land plants. In contrast, reports of small inversions are rare and involve limited plant groups. In this study, we report the widespread occurrence of small inversions ranging from 5 to 50 bp in fully and partially sequenced chloroplast genomes of both monocots and dicots. We found that small inversions were much more common than large inversions. The small inversions were scattered over the chloroplast genome including the IR, SSC, and LSC regions. Several small inversions were uncovered in chloroplast genomes even though they shared the same overall gene order. The majority of these small inversions were located within 100 bp downstream of the 3' ends of genes. All had inverted repeat sequences, ranging from 11 to 24 bp, at their ends. Such small inversions form stem-loop hairpin structures that usually have the function of stabilizing the corresponding mRNA molecules. Intra-molecular recombination between the inverted sequences in the stem-forming regions are responsible for generating flip-flop orientations of the loops. The presence of two different orientations of the stem-loop in the trnL-F noncoding region of a single species of Jasminum elegans suggests that a short inversion can be generated within a short period of time. Small inversions of non-coding sequences may influence sequence alignment and character interpretation in phylogeny reconstructions, as shown in nine species of Jasminum. Many small inversions may have been generated by parallel or back mutation events during chloroplast genome evolution. Our data indicate that caution is needed when using chloroplast non-coding sequences for phylogenetic analysis.     

Analysis of the human <Emphasis Type="Italic">Alu</Emphasis> Ye lineage

Background  

Alu elements are short (~300 bp) interspersed elements that amplify in primate genomes through a process termed retroposition. The expansion of these elements has had a significant impact on the structure and function of primate genomes. Approximately 10 % of the mass of the human genome is comprised of Alu elements, making them the most abundant short interspersed element (SINE) in our genome. The majority of Alu amplification occurred early in primate evolution, and the current rate of Alu retroposition is at least 100 fold slower than the peak of amplification that occurred 30–50 million years ago. Alu elements are therefore a rich source of inter- and intra-species primate genomic variation.     

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.     

Size-dependent palindrome-induced intrachromosomal recombination in yeast

Palindromic and quasi-palindromic sequences are important DNA motifs found in various cis-acting genetic elements, but are also known to provoke different types of genetic alterations. The instability of such motifs is clearly size-related and depends on their potential to adopt secondary structures known as hairpins and cruciforms. Here we studied the influence of palindrome size on recombination between two directly repeated copies of the yeast CYC1 gene leading to the loss of the intervening sequence (“pop-out” recombination). We show that palindromes inserted either within one copy or between the two copies of the CYC1 gene become recombinogenic only when they attain a certain critical size and we estimate this critical size to be about 70 bp. With the longest palindrome used in this study (150 bp) we observed a more than 20-fold increase in the pop-out recombination. In the sae2/com1 mutant the palindrome-stimulated recombination was completely abolished. Suppression of palindrome recombinogenicity may be crucial for the maintenance of genetic stability in organisms containing a significant number of large palindromes in their genomes, like humans.     

Internal deletions of transposable elements: the case of Lemi elements

Mobile elements using a “cut and paste” mechanism of transposition (Class II) are frequently prone to internal deletions and the question of the origin of these copies remains elusive. In this study, we looked for copies belonging to the Lemi Family (Tc1-mariner-IS630 SuperFamily) in the plant genomes, and copies within internal deletions were analyzed in detail. Lemi elements are found exclusively in Eudicots, and more than half of the copies have been deleted. All deletions occur between microhomologies (direct repeats from 2 to 13 bp). Copies less than 500 bp long, similar to MITEs, are frequent. These copies seem to result from large deletions occurring between microhomologies present within a region of 300 bp at both extremities of the element. These regions are particularly A/T rich, compared to the internal part of the element, which increases the probability of observing short direct repeats. Most of the molecular mechanisms responsible for double strand break repair are able to induce deletions between microhomologies during the repair process. This could be a quick way to reduce the population of active copies within a genome and, more generally, to reduce the overall activity of the element after it has entered a naive genome.     

A highly conserved, small LTR retrotransposon that preferentially targets genes in grass genomes

LTR retrotransposons are often the most abundant components of plant genomes and can impact gene and genome evolution. Most reported LTR retrotransposons are large elements (>4 kb) and are most often found in heterochromatic (gene poor) regions. We report the smallest LTR retrotransposon found to date, only 292 bp. The element is found in rice, maize, sorghum and other grass genomes, which indicates that it was present in the ancestor of grass species, at least 50-80 MYA. Estimated insertion times, comparisons between sequenced rice lines, and mRNA data indicate that this element may still be active in some genomes. Unlike other LTR retrotransposons, the small LTR retrotransposons (SMARTs) are distributed throughout the genomes and are often located within or near genes with insertion patterns similar to MITEs (miniature inverted repeat transposable elements). Our data suggests that insertions of SMARTs into or near genes can, in a few instances, alter both gene structures and gene expression. Further evidence for a role in regulating gene expression, SMART-specific small RNAs (sRNAs) were identified that may be involved in gene regulation. Thus, SMARTs may have played an important role in genome evolution and genic innovation and may provide a valuable tool for gene tagging systems in grass.     

Recent amplification of miniature inverted-repeat transposable elements in the vector mosquito Culex pipiens: characterization of the Mimo family

We describe a new family of repetitive elements, named Mimo, from the mosquito Culex pipiens. Structural characteristics of these elements fit well with those of miniature inverted-repeat transposable elements (MITEs), which are ubiquitous and highly abundant in plant genomes. The occurrence of Mimo in C. pipiens provides new evidence that MITEs are not restricted to plant genomes, but may be widespread in arthropods as well. The copy number of Mimo elements in C. pipiens (1000 copies in a 540 Mb genome) supports the hypothesis that there is a positive correlation between genome size and the magnitude of MITE proliferation. In contrast to most MITE families described so far, members of the Mimo family share a high sequence conservation, which may reflect a recent amplification history in this species. In addition, we found that Mimo elements are a frequent nest for other MITE-like elements, suggesting that multiple and successive MITE transposition events have occurred very recently in the C. pipiens genome. Despite evidence for recent mobility of these MITEs, no element has been found to encode a protein; therefore, we do not know how they have transposed and have spread in the genome. However, some sequence similarities in terminal inverted-repeats suggest a possible filiation of some of these mosquito MITEs with pogo-like DNA transposons.     

Structural organization of the human gastrointestinal glutathione peroxidase (GPX2) promoter and 3'-nontranscribed region: transcriptional response to exogenous redox agents

We describe a new family of repetitive elements, named Mimo, from the mosquito Culex pipiens. Structural characteristics of these elements fit well with those of miniature inverted-repeat transposable elements (MITEs), which are ubiquitous and highly abundant in plant genomes. The occurrence of Mimo in C. pipiens provides new evidence that MITEs are not restricted to plant genomes, but may be widespread in arthropods as well. The copy number of Mimo elements in C. pipiens (1000 copies in a 540 Mb genome) supports the hypothesis that there is a positive correlation between genome size and the magnitude of MITE proliferation. In contrast to most MITE families described so far, members of the Mimo family share a high sequence conservation, which may reflect a recent amplification history in this species. In addition, we found that Mimo elements are a frequent nest for other MITE-like elements, suggesting that multiple and successive MITE transposition events have occurred very recently in the C. pipiens genome. Despite evidence for recent mobility of these MITEs, no element has been found to encode a protein; therefore, we do not know how they have transposed and have spread in the genome. However, some sequence similarities in terminal inverted-repeats suggest a possible filiation of some of these mosquito MITEs with pogo-like DNA transposons.     

Deletion of ultraconserved elements yields viable mice

Ultraconserved elements have been suggested to retain extended perfect sequence identity between the human, mouse, and rat genomes due to essential functional properties. To investigate the necessities of these elements in vivo, we removed four noncoding ultraconserved elements (ranging in length from 222 to 731 base pairs) from the mouse genome. To maximize the likelihood of observing a phenotype, we chose to delete elements that function as enhancers in a mouse transgenic assay and that are near genes that exhibit marked phenotypes both when completely inactivated in the mouse and when their expression is altered due to other genomic modifications. Remarkably, all four resulting lines of mice lacking these ultraconserved elements were viable and fertile, and failed to reveal any critical abnormalities when assayed for a variety of phenotypes including growth, longevity, pathology, and metabolism. In addition, more targeted screens, informed by the abnormalities observed in mice in which genes in proximity to the investigated elements had been altered, also failed to reveal notable abnormalities. These results, while not inclusive of all the possible phenotypic impact of the deleted sequences, indicate that extreme sequence constraint does not necessarily reflect crucial functions required for viability.