MysTR: an endogenous retrovirus family in mammals that is undergoing recent amplifications to unprecedented copy numbers

A large percentage of the repetitive elements in mammalian genomes are retroelements, which have been moved primarily by LINE-1 retrotransposons and endogenous retroviruses. Although LINE-1 elements have remained active throughout the mammalian radiation, specific groups of endogenous retroviruses generally remain active for comparatively shorter periods of time. Identification of an unusual extinction of LINE-1 activity in a group of South American rodents has opened a window for examination of the interplay in mammalian genomes between these ubiquitous retroelements. In the course of a search for any type of repetitive sequences whose copy numbers have substantially changed in Oryzomys palustris, a species that has lost LINE-1 activity, versus Sigmodon hispidus, a closely related species retaining LINE-1 activity, we have identified an endogenous retrovirus family differentially amplified in these two species. Analysis of three full-length, recently transposed copies, called mysTR elements, revealed gag, pro, and pol coding regions containing stop codons which may have accumulated either before or after retrotransposition. Isolation of related sequences in S. hispidus and the LINE-1 active outgroup species, Peromyscus maniculatus, by PCR of a pro-pol region has allowed determination of copy numbers in each species. Unusually high copy numbers of approximately 10,000 in O. palustris versus 1,000 in S. hispidus and 4,500 in the more distantly related P. maniculatus leave open the question of whether there is a connection between endogenous retrovirus activity and LINE-1 inactivity. Nevertheless, these independent expansions of mysTR represent recent amplifications of this endogenous retrovirus family to unprecedented levels.     

Replicon fusions promoted by the inverted repeats of Tn5. The right repeat is an insertion sequence

Members of three repetitive sequence families were isolated from recombinant λ-genome libraries, and were used to investigate sequence relationships within these families. Studies presented elsewhere show that members of all three repeat sequence families are transcribed tissue-specifically. The thermal stability of intrafamilial heteroduplexes was measured, and the extent of colinearity between related sequences was determined by restriction mapping, heteroduplex visualization, gel blot hybridization, and direct sequencing. One large and very divergent family, named 2108, was shown to consist of an assemblage of many small repeat sequence subfamilies. Each subfamily includes <40 members which are not contiguous in the genome but are very closely related colinear sequence elements several thousand nucleotides in length. The different 2108 subfamilies share only small sequence subelements, which in each subfamily occur in a different linear order and are surrounded by different sequences. A second divergent family consisting of short repetitive sequences, the 2109 family, includes many small internally homologous subfamilies as well. A third family, 2034, displays little internal sequence divergence and no apparent subfamily structure. The repeat sequence subfamilies may be biologically significant units of repetition. Thus specific 2108 subfamilies were shown to be evolutionary conserved to a remarkable degree. Highly homologous 2108 sequences were found shared among sea urchin species which diverged almost 200 million years ago, although only about 10% of the single copy DNA sequences of these species are now homologous enough to crossreact.     

The olfactory receptor gene superfamily: data mining, classification, and nomenclature

The vertebrate olfactory receptor (OR) subgenome harbors the largest known gene family, which has been expanded by the need to provide recognition capacity for millions of potential odorants. We implemented an automated procedure to identify all OR coding regions from published sequences. This led us to the identification of 831 OR coding regions (including pseudogenes) from 24 vertebrate species. The resulting dataset was subjected to neighbor-joining phylogenetic analysis and classified into 32 distinct families, 14 of which include only genes from tetrapodan species (Class II ORs). We also report here the first identification of OR sequences from a marsupial (koala) and a monotreme (platypus). Analysis of these OR sequences suggests that the ancestral mammal had a small OR repertoire, which expanded independently in all three mammalian subclasses. Classification of ``fish-like' (Class I) ORs indicates that some of these ancient ORs were maintained and even expanded in mammals. A nomenclature system for the OR gene superfamily is proposed, based on a divergence evolutionary model. The nomenclature consists of the root symbol `OR', followed by a family numeral, subfamily letter(s), and a numeral representing the individual gene within the subfamily. For example, OR3A1 is an OR gene of family 3, subfamily A, and OR7E12P is an OR pseudogene of family 7, subfamily E. The symbol is to be preceded by a species indicator. We have assigned the proposed nomenclature symbols for all 330 human OR genes in the database. A WWW tool for automated name assignment is provided. Received: / Accepted:     

Cytogenetic and molecular evaluation of centromere-associated DNA sequences from a marsupial (Macropodidae: Macropus rufogriseus) X chromosome

The constitution of the centromeric portions of the sex chromosomes of the red-necked wallaby, Macropus rufogriseus (family Macropodidae, subfamily Macropodinae), was investigated to develop an overview of the sequence composition of centromeres in a marsupial genome that harbors large amounts of centric and pericentric heterochromatin. The large, C-band-positive centromeric region of the X chromosome was microdissected and the isolated DNA was microcloned. Further sequence and cytogenetic analyses of three representative clones show that all chromosomes in this species carry a 178-bp satellite sequence containing a CENP-B DNA binding domain (CENP-B box) shown herein to selectively bind marsupial CENP-B protein. Two other repeats isolated in this study localize specifically to the sex chromosomes yet differ in copy number and intrachromosomal distribution. Immunocytohistochemistry assays with anti-CENP-E, anti-CREST, anti-CENP-B, and anti-trimethyl-H3K9 antibodies defined a restricted point localization of the outer kinetochore at the functional centromere within an enlarged pericentric and heterochromatic region. The distribution of these repeated sequences within the karyotype of this species, coupled with the apparent high copy number of these sequences, indicates a capacity for retention of large amounts of centromere-associated DNA in the genome of M. rufogriseus.     

Divergence of unique DNA sequences of bivalve mollusks from Mytilinae subfamily (Bivalvia mytilidae)

The DNA-DNA hybridization method was used to determine the divergence degree of unique sequences of five bivalvia species belonging to the Mytilinae subfamily. The matrix of delta Tm values for heteroduplexes of unique sequences was found which made it possible to define three phylogenetic branches within the subfamily. Under non-stringent hybridization conditions (55 degrees C, 0.5M PB) the divergence between species of any two branches was about 14% of nucleotide substitutions. The hybridization of [3H] unique sequences of Mytilinae with DNA fragments of Modiolus modiolus, a representative of the closest relative Modiolinae subfamily, showed that the divergence rate of unique sequences in two phylogenetic lines of the Mytilus genusis higher than in the line of the Crenomytilus genus I0.35-0.23% and 0.1% of nucleotide substitutions per one million years, respectively). According to the matrix of delta Tm values for five species of Mytilinae and Modiolus modiolus, a phylogenetic tree was built reflecting the differences between the divergence rates in different branches.     

Satellite-DNA diversification and the evolution of major lineages in Cardueae (Carduoideae Asteraceae)

In a previous work, we characterized the HinfI satellite DNA family in the subtribe Centaureinae (Cardueae) demonstrating that a “library” of eight HinfI subfamilies would exist in the common ancestor of all Centaureinae, which were differentially amplified in different lineages. Now, we extend our study by analyzing a total of 219 additional repeats from fifteen species belonging to Carlininae, Echinopsinae and Carduinae, and comparing them to those of Centaureinae. Most HinfI sequences belonged to the subfamily II, although a few sequences of other subfamilies were detected in some species. Additionally, a new subfamily characteristic of several Carduinae species was discovered. Although phylogenetic trees grouped sequences by subfamily affinity instead of species provenance, when comparing repeats of the same subfamily, the degree of divergence between any pair of sequences was related to the evolutionary distance between the species compared in most cases. Exceptions were in comparisons between sequences of some Centaureinae species, and between sequences of some Carduinae species and those of Centaureinae. Our results demonstrate that: (1) At least nine HinfI subfamilies would exist in the common ancestor of Cardueae, each one differentially amplified in different lineages; (2) After differential spreading, sequences of each subfamily evolved concertedly through molecular drive, resulting in the gradual divergence of repeats between different species; (3) The rate to which concerted evolution occurred was different between lineages according to the evolutionary history of each one.     

Evolutionary dynamics of satellite DNA family PIM357 in species of the genus Pimelia (Tenebrionidae,Coleoptera)

A large number of repeats of a satellite DNA (stDNA) family have been cloned and sequenced from species and populations of the genus Pimelia (Tenebrionidae, Coleoptera). The beetles were collected in the Canary Islands, Morocco, the Iberian Peninsula, and the Balearic Islands in order to analyze the evolutionary forces and processes acting on abundant stDNAs conserved at the genus level. This repetitive family is composed of an abundant A-T-rich stDNA, with basic units of 357 bp. All the sequences obtained showed similarity to the 22 repeat units of the PIM357 stDNA family described previously for six Iberian Pimelia species (Pons et al. 1997 ). An analysis based on similarity shows the presence of three different groups of sequences clearly in accordance with their geographical origin. One is composed of satellite sequences from Iberian and Balearic species, a second group from the Moroccan taxa, whereas the third one is from the Pimelia species endemic to the Canary Islands. The latter group shows higher nucleotide diversities for their stDNA sequences and a lack of relationship between transition stages to fixation and sequence divergence. Phylogeographic data of Canarian Pimelia show that the PIM357 stDNA family has persisted for more than 8 Myr and could probably be traced to the origin of the lineage. The data suggest that distinct demographic and phylogenetic patterns related to the colonization of the volcanic Canarian island chain account for particular evolutionary dynamics of the repeat DNA family in this group.     

Identification,characterization, and utilization of single copy genes in 29 angiosperm genomes

Background

Single copy genes are common across angiosperm genomes. With the sufficiently high quality sequenced genomes, the identification of large-scale single copy genes among multiple species is possible. Although some characteristics have been reported, our study provides novel insights into single copy genes.

Results

We identified single copy genes across 29 angiosperm genomes. A significant negative correlation was found between the number of duplicate blocks and the number of single copy genes. We found that a considerable number of single copy genes are located in organelles, showing a preference for binding and catalytic activity. The analysis of effective number of codons (Nc) illustrates that single copy genes have a stronger codon bias than non-single copy genes in eudicots. The relative high expression level of single copy genes was partially confirmed by the RNA-seq data, rather than the Codon Adaptation Index (CAI). Unlike in most other species, a strongly negatively correlation occurs between Nc and GC3 among single copy genes in grass genomes. When compared to all non-single copy genes, single copy genes indicate more conservation (as indicated by Ka and Ks values). But our alternative splicing (AS) results reveal that selective constraints are weaker in single copy genes than in low copy family genes (1–10 in-paralogs) and stronger than high copy family genes (>10 in-paralogs). Using concatenated shared single copy genes, we obtained a well-resolved phylogenetic tree. With the addition of intron sequences, the branch support is improved, but striking incongruences are also evident. Therefore, it is noteworthy that inclusion of intron sequences seems more appropriate for the phylogenetic reconstruction at lower taxonomic levels.

Conclusions

Our analysis provides insight into the evolutionary characteristics of single copy genes across 29 angiosperm genomes. The results suggest that there are key differences in evolutionary constraints between single copy genes and non-single copy genes. And to some extent, these evolutionary constraints show some species-specific differences, especially between eudicots and monocots. Our preliminary evidence also suggests that the concatenated shared single copy genes are well suited for use in resolving phylogenetic relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-504) contains supplementary material, which is available to authorized users.     

Genus-specific localization of the TaiI family of tandem-repetitive sequences in either the centromeric or subtelomeric regions in Triticeae species (Poaceae) and its evolution in wheat.

The TaiI family sequences are classified as tandem repetitive DNA sequences present in the genome of tribe Triticeae, and are localized in the centromeric regions of common wheat, but in the subtelomeric heterochromatic regions of Leymus racemosus and related species. In this study, we investigated the chromosomal distribution of TaiI family sequences in other Triticeae species. The results demonstrated a centromeric localization in genera Triticum and Aegilops and subtelomeric localization in other genera, thus showing a genus-dependent localization of TaiI family sequences in one or the other region. The copy numbers of TaiI family sequences in species in the same genus varied greatly, whether in the centromeric or subtelomeric regions (depending on genus). We also examined the evolution of TaiI family sequences during polyploidization of hexaploid common wheat. A comparison of chromosomal locations of the major TaiI family signals in common wheat and in its ancestral species suggested that the centromeric TaiI family sequences in common wheat were inherited from its ancestors with little modification, whereas a mixed origin for the B genome of common wheat was indicated.     

Evolution of the Xenopus piggyBac transposon family TxpB: domesticated and untamed strategies of transposon subfamilies

A new family, termed TxpB, of DNA transposons belonging to the piggyBac superfamily was found in 3 Xenopus species (Xenopus tropicalis, Xenopus laevis, and Xenopus borealis). Two TxpB subfamilies of Kobuta and Uribo1 were found in all the 3 species, and another subfamily termed Uribo2 was found in X. tropicalis. Molecular phylogenetic analyses of their open reading frames (ORFs) revealed that TxpB transposons have been maintained for over 100 Myr. Both the Uribo1 and the Uribo2 ORFs were present as multiple copies in each genome, and some of them were framed by terminal inverted repeat sequences. In contrast, all the Kobuta ORFs were present as a single copy in each genome and exhibited high evolutionary conservation, suggesting domestication of Kobuta genes by the host. Genomic insertion polymorphisms of the Uribo1 and Uribo2 transposons (nonautonomous type) were observed in a single species of X. tropicalis, indicating recent transposition events. Transfection experiments in cell culture revealed that an expression vector construct for the intact Uribo2 ORF caused precise excision of a nonautonomous Uribo2 element from the target vector construct but that for the Kobuta ORF did not. The present results support our viewpoint that some Uribo2 members are naturally active autonomous transposons, whereas Kobuta members may be domesticated by hosts.     

Protein tyrosine phosphatase domains from the protochordate Styela plicata

Protein tyrosine phosphorylation is an important regulatory mechanisms in cell physiology. While the protein tyrosine kinase (PTKase) family has been extensively studied, only six protein tyrosine phosphatases (PTPases) have been described. By Southern blot analysis, genomic DNA from several different phyla were found to cross-hybridize with a cDNA probe encoding the human leukocyte-common antigen (LCA; CD45) PTPase domains. To pursue this observation further, total mRNA from the protochordate Styela plicata was used as a tempalte to copy and amplify, using polymerase chain reaction (PCR) technology, PTPase domains. Twenty-seven distinct sequences were identified that contain hallmark residues of PTPases; two of these are similar to described mammalian PTPases. Southern blot analysis indicates that at least one other Styela sequence is highly conserved in a variety of phyla. Seven of the Styela domains have significant similarity to each other, indicating a subfamily of PTPases. However, most of the sequences are disparate. A comparison of the 27 Styela sequences with the ten known PTPase domain sequences reveals that only three residues are absolutely conserved and identifies regions that are highly divergent. The data indicate that the PTPase family will be equally as large and diverse as the PTKases. The extent and diversity of the PTPase family suggests that these enzymes are, in their own right, important regulators of cell behavior.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M37986-M38041.     

Genome-wide identification and analysis of membrane-bound O-acyltransferase (MBOAT) gene family in plants

Membrane bound O-acyl transferase (MBOAT) family is composed of gene members encoding a variety of acyltransferase enzymes, which play important roles in plant acyl lipid metabolism. Here, we present the first genome-enabled identification and analysis of MBOAT gene models in plants. In total, we identified 136 plant MBOAT sequences from 14 plant species with complete genomes. Phylogenetic relationship analyses suggested the plant MBOAT gene models fell into four major groups, two of which likely encode enzymes of diacylglycerol acyltransferase 1 (DGAT1) and lysophospholipid acyltransferase (LPLAT), respectively, with one–three copies of paralogs present in each of the most plant species. A group of gene sequences, which are homologous to Saccharomyces cerevisiae glycerol uptake proteins (GUP), was identified in plants; copy numbers were conserved, with only one copy represented in each of the most plant species; analyses showed that residues essential for acyltransferases were more prone to be conserved than vertebrate orthologs. Among four groups, one was inferred to emerge in land plants and experience a rapid expansion in genomes of angiosperms, which suggested their important roles in adaptation of plants in lands. Sequence and phylogeny analyses indicated that genes in all four groups encode enzymes with acyltransferases. Comprehensive sequence identification of MBOAT family members and investigation into classification provide a complete picture of the MBOAT gene family in plants, and could shed light into enzymatic functions of different MBOAT genes in plants.     

Molecular structure and chromosomal localization of major repetitive DNA families in the chickpea (Cicer arietinum L.) genome

Three major repetitive DNA sequences were isolated from a genomic library of chickpea (Cicer arietinum L.) and characterized with respect to their genomic organization and chromosomal localization. All repetitive elements are genus-specific and mostly located in the AT-rich pericentric heterochromatin. Two families are organized as satellite DNAs with repeat lengths of 162–168 bp (CaSat1) and 100 bp (CaSat2). CaSat1 is mainly located adjacent to the 18S rDNA clusters on chromosomes A and B, whereas CaSat2 is a major component of the pericentric heterochromatin on all chromosomes. The high abundance of these sequences in closely related species of the genus Cicer as well as their variation in structure and copy number among the annual species provide useful tools for taxonomic studies. The retrotransposon-like sequences of the third family (CaRep) display a more complex organization and are represented by two independent sets of clones (CaRep1 and CaRep2) with homology to different regions of Ty3-gypsy-like retrotransposons. They are distributed over the pericentric heterochromatin block on all chromosomes with extensions into euchromatic regions. Conserved structures within different crossability groups of related Cicer species suggest independent amplification or transposition events during the evolution of the annual species of the genus.     

Properties of satellite DNA of three plant species of the subtribe Citrinae

The properties of stDNA of three species of the subtribe Citrinae have been investigated. The buoyant density of the main component is 1.693 g/cm3, that of satellite component is 1.712 g/cm3 and 1.715 g/cm3. The differential melting curves of satellite components reveal two melting zones. Some of stDNAs are melted within a broad temperature range, while others--at higher temperatures as a narrow peak. The reassociation kinetics suggest that 50-54% of stDNA are a fast reassociating fraction with the length of repeating sequences of 0.8-1.1 x 10(2) base pairs. Based on the values of Tm and buoyant density the 5-methylcytosine content in stDNA was calculated and was found equal to 20-35%. Using equilibrium ultracentrifugation in the actinomycin D--CsCl density gradient the stDNAs of the subtribe Citrinae were separated into constituent components.     

Characterization and Evolution of mariner Elements from Closely Related Species of Fruit Flies (Diptera: Tephritidae)

Mariner elements were amplified using the polymerase chain reaction from two species of tephritid flies, Ceratitis rosa and Trirhithrum coffeae. The sequences were ∼1.3 kb in length. None of these elements appeared to be functional, as in every case the open reading frame (ORF) was disrupted by the presence of frameshifts or stop codons. These elements, Crmar1 and Tcmar1, are very similar to the Ccmar1 element previously amplified from the closely related tephritid species C. capitata and are members of the mellifera subfamily of mariner elements. The phylogeny and pattern of divergence of these elements were examined in relation to the phylogeny of the host species. It is highly probable that the elements were present in the ancestral lineage prior to the divergence of the three species. The copy numbers of the elements within each species are very different, ranging from about 10 in T. coffeae to 5,000 in C. rosa. The possible mechanisms which determine the copy number of an element in the host genome are discussed. Received: 25 April 1997 / Accepted: 31 July 1997