Human-specific subfamilies of HERV-K (HML-2) long terminal repeats: three master genes were active simultaneously during branching of hominoid lineages

Using 40 known human-specific LTR sequences, we have derived a consensus sequence for an evolutionary young HERV-K (HML-2) LTR family, which was named the HS family. In the human genome the HS family is represented by approximately 150-160 LTR sequences, 90% of them being human-specific (hs). The family can be subdivided into two subfamilies differing in five linked nucleotide substitutions: HS-a and HS-b of 5.8 and 10.3 Myr evolutionary ages, respectively. The HS-b subfamily members were transpositionally active both before the divergence of the human and chimpanzee ancestor lineages and after it in both lineages. The HS-a subfamily comprises only hs LTRs. These and other data strongly suggest that at least three "master genes" of HERV-K (HML-2) LTRs were active in the human ancestor lineage after the human-chimpanzee divergence. We also found hs HERV-K (HML-2) LTRs integrations in introns of 12 human genes and identified 13 new hs HERV-K (HML-2) LTRs.     

Subfamilies and nearest-neighbour dendrogram for the LTRs of human endogenous retroviruses HERV-K mapped on human chromosome 19: physical neighbourhood does not correlate with identity level

Sequences of 45 long terminal repeats (LTRs) of the human endogenous retroviruses HERV-K family, precisely mapped by us earlier on human chromosome 19, were determined and a nearest-neighbour dendrogram was constructed. No correlation was observed between the degree of identity of the LTR pairs and their relative positions on the chromosome. Thus, sequences of distantly located LTRs, even positioned on different chromosome arms, could be highly similar to each other, whereas those of closely located LTRs could differ significantly. We conclude that the LTRs have randomly transposed across the chromosome in the course of evolution. The alignment of the LTR sequences allowed us to assign most of the LTRs to two major subfamilies. The LTRs belonging to the first subfamily (LTR-I) are characterised by higher intrasubfamily sequence divergence than those of the second subfamily (LTR-II). The two subfamilies are easily distinguished by the presence of characteristic deletions/insertions in the LTR sequences. The higher divergence of the first subfamily members suggests that their propagation started at earlier stages of evolution, probably soon after the insertion of their ancestral sequence into the primate genome. In turn, each of the subfamilies includes several distinct branches with various degrees of intragroup divergence and with characteristic diagnostic features, suggesting that the members of the branches represent amplified copies of particular master genes which had appeared at different periods of evolution. The sequences of the LTRs demonstrate a characteristic distribution of conservative and variable regions, indicating that the LTRs might have some sequence-dependent functions in the primate genome. Received: 11 August 1997 / Accepted: 22 September 1997     

Characterization of the intragenomic spread of the human endogenous retrovirus family HERV-W

This study examines the intragenomic spread of the human endogenous retrovirus family HERV-W from insertions present within the draft sequence of the human genome. Identification of shared diagnostic differences and phylogenetic analyses revealed the existence of three main subfamilies. The average divergence between sequences for each of the subfamilies suggests that most of the HERV-W elements were inserted within the genome during a short period of evolutionary time. Each one of the subfamilies consists of two types of insertions, the expected proviral sequences and other sequences resembling the structure of processed retrogenes. These HERV-W retrosequences extend from the R region of the 5' long-terminal repeat (LTR) to the R region of the 3' LTR (as viral genomic RNAs), end in poly(A) 3' tails, and are flanked by direct repeats longer than the proviral integrations. Furthermore, several of the HERV-W retrosequences are 5'-truncated at different sites. I suggest the involvement of the L1 machinery in these integrations and discuss the characteristic features of the evolutionary history of HERV-W, with emphasis on the putative impact of HERV-W retrosequence integrations on the mammalian genome.     

Species specificity of macaque rhadinovirus glycoprotein B sequences

All members of the Herpesviridae family contain sequences for a highly conserved glycoprotein B (gB) gene. We investigated the phylogenetic relationships of gB sequences from eight independent rhadinovirus isolates obtained from three species: rhesus (Macaca mulatta), cynomologus (Macaca fasicularis), and pig-tailed (Macaca nemestrina) macaques. Samples were derived from monkeys housed at four separate facilities. Analysis of these eight independent gB sequences revealed five regions of heterogeneity within the 823- to 829-amino-acid polypeptides: residues 1 to 65, 120 to 185, 255 to 300, 352 to 393, and 412 to 457. The remaining regions of gB were highly conserved among the different macaque isolates. Overall divergence among these gene sequences ranged from 0.1 to 7.2% at the amino acid level. Phylogenetic trees constructed with our macaque rhadinovirus gB sequences and those derived from additional subfamilies or genera (alpha, beta, gamma-1, and gamma-2) revealed that the macaque gB sequences branched with other gamma-2 herpesvirus gB sequences and that within the gamma-2 genera, the macaque gB sequences clustered as a distinct branch. The eight macaque rhadinovirus gB sequences were all approximately equidistant from Kaposi sarcoma-associated herpesvirus (KSHV) gB sequences and had a shorter evolutionary distance to KSHV gB sequences than to any other herpesvirus, including the gamma-2 herpesvirus saimiri (HVS) of New World squirrel monkeys. The macaque gB sequences did not cluster according to the facility of origin, but did cluster according to the species of origin, displaying less intraspecies divergence (0.1 to 2.9%) than interspecies divergence (3.3 to 7.2%). These results demonstrate a close relatedness of rhadinovirus isolates from different macaque species.     

Diverse long terminal repeats are associated with murine retroviruslike (VL30) elements.

The VL30 family is a retroviruslike gene family with no apparent nucleic acid homology to any known retrovirus. Over 100 copies of VL30 DNA elements are dispersed throughout the mouse genome. Sequence analysis of the VL30 long terminal repeat (LTR) units showed that, whereas the LTR units of a given VL30 DNA element were almost identical, the LTR units associated with distinct members of the family were very different from one another. Comparison of the LTR sequences possessed by two particular VL30 DNA elements revealed a pattern of extensively homologous DNA segments adjacent to only distantly related DNA sequences. With the aid of sub-LTR probes, it was shown that a certain LTR is composed of both U5 sequences that are abundantly present in all species of the genus Mus and a U3 region detected only in Mus musculus. In addition, we isolated a VL30 DNA element in which the LTR units were replaced by the LTR units of an apparently novel retroviruslike family. These findings suggest that recombinations have played a role in generating the diverse population of VL30-associated LTRs.     

Subfamily divergence in the multigene family of ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcS) in Triticeae and its relatives

To investigate genetic mechanisms acting on multigene family in plants, we analyzed sequence variation in the rbcS gene of 13 species of Triticeae and one species each of related tribes (Bromeae and Aveneae). A total of 36 rbcS genes were analyzed. Based on dimorphism in the length of intron, the rbcSs of investigated species were classified into two subfamilies A and B. The difference in intron length was caused by an indel of about 200 bp in the middle of the intron. The two subfamilies of rbcS were present in the three tribes, indicating that the divergence of rbcS subfamilies occurred before the split of these tribes. Generally, variation between the two subfamilies of rbcS was larger than that within subfamily, but these two measures were about the same at the tribe level. This result suggested that divergence of the subfamilies of rbcS occurred at about the same time of tribe diversification. The level of nucleotide variation in the exon region between subfamilies was reduced in the Triticeae, but clear change was not detected in the intron sequence. This result suggested that the exon sequences between subfamilies of rbcS were homogenized without affecting the intron sequence in the Triticeae lineage.     

Comparison of mRNA and protein levels of four members of the protein 4.1 family: the type II brain 4.1/4.1B/KIAA0987 is the most predominant member of the protein 4.1 family in rat brain

The human gene for the fourth member of the protein 4.1 family, KIAA0987, was recently identified by comprehensive cDNA analysis. To further characterize the corresponding gene and its product in rats, we cloned and sequenced rat KIAA0987 cDNA. RNA blot analyses revealed that the rat KIAA0987 gene was abundantly expressed only in the brain, kidney, and testis. Although we have previously reported that the third member of the protein 4.1 family, the KIAA0338 gene product, is predominantly expressed in rat brain, and thus was named brain 4.1, quantitative RNA blot analyses indicated that KIAA0987 should be called something other than brain 4.1 because the level of KIAA0987 mRNA was found to be of the same order as that of KIAA0338 mRNA. Our quantitative immunoblot analysis showed that the most predominant member of the protein 4.1 family at the protein level was the product of the KIAA0987 gene, not that of the KIAA0338 gene. Taking these results together, we consider it reasonable to name the KIAA0338 and KIAA0987 gene products 'type I brain 4.1' and 'type II brain 4.1,' respectively, because these two products were found to be more prominently produced in rat brain than the other two members of the protein 4.1 family, erythroid 4.1 and 4.1G.     

Characterization of the repetitive sequences in a 200-kb region around the rice waxy locus: diversity of transposable elements and presence of veiled repetitive sequences

Repetitive genomic sequences might have various structural features and properties distinct from those of the known transposable elements (TE). Here, the content and properties of the repetitive sequences present in a 200-kb region around the rice waxy locus were analyzed using the available rice genomic database. In our previous Southern blotting analysis, 70% of the segments in this region showed smeared patterns, but according to the present database analysis, the proportion of repetitive sequences in this region was only 15%. The repetitive segments in this 200-kb region comprised 75 repetitive sequences that we classified into 46 subfamilies: 21 subfamilies were known TEs or repetitive sequences and 25 subfamilies consisted of newly identified TEs or novel types of repetitive sequences. The region contains no long terminal repeat (LTR) retrotransposable elements, but miniature inverted repeat transposable elements (MITEs) constituted a major class among the elements identified. These MITEs showed remarkable structural divergence: 12 elements were found to be new members of known MITE superfamilies, while five elements had novel terminal structures, and did not belong to any known TE families. Interestingly, about 10% of the repetitive sequences, including virus-like sequences did not have any of the usual characteristics of TEs, suggesting that a certain proportion of repetitive sequences that might not share the transpositional mechanisms of known elements are dispersed in the compact rice genome.     

人X染色体含有一个黑色素瘤抗原基因亚家族

肿瘤相关基因的研究是肿瘤基因形成学说的核心内容。肿瘤相关基因家族的研究则是其中的重点和难点,从4－6月孕龄人胎肝cDNA文库中克隆到一个黑色素瘤抗原基因亚家族,称为MAGE－D亚家族,其成员包括3个直系同源体（人MAGE－D1、大鼠SNERG－1和小鼠DLXIN－1）和2个旁系同源体（人MAGE－D和人KIAA1114）。该家族的3个人类成员均定位于染色体Xp11.21-p11.23,同时具有独特的基因组结构。分子进化树分析表明,该家族与已知MAGE－A、－B和－C3个亚家族之间具有明显的进化上分歧。该亚家族的发现为研究肿瘤相关基因新功能提供了重要线索。     

The evolutionary analysis of the Tnt1 retrotransposon in Nicotiana species reveals the high variability of its regulatory sequences

We studied the evolution of the tobacco Tnt1 retrotransposon by analyzing Tnt1 partial sequences containing both coding domains and U3 regulatory sequences obtained from a number of Nicotiana species. We detected three different subfamilies of Tnt1 elements, Tnt1A, Tnt1B, and Tnt1C, that differ completely in their U3 regions but share conserved flanking coding and LTR regions. U3 divergence between the three subfamilies is found in the region that contains the regulatory sequences that control the expression of the well-characterized Tnt1-94 element. This suggests that expression of the three Tnt1 subfamilies might be differently regulated. The three Tnt1 subfamilies were present in the Nicotiana genome at the time of species divergence, but have evolved independently since then in the different genomes. Each Tnt1 subfamily seems to have conserved its ability to transpose in a limited and different number of Nicotiana species. Our results illustrate the high variability of Tnt1 regulatory sequences. We propose that this high sequence variability could allow these elements to evolve regulatory mechanisms in order to optimize their coexistence with their host genome.      

Genome-Wide Characterization of Endogenous Retroviruses in the Bat Myotis lucifugus Reveals Recent and Diverse Infections

Bats are increasingly recognized as reservoir species for a variety of zoonotic viruses that pose severe threats to human health. While many RNA viruses have been identified in bats, little is known about bat retroviruses. Endogenous retroviruses (ERVs) represent genomic fossils of past retroviral infections and, thus, can inform us on the diversity and history of retroviruses that have infected a species lineage. Here, we took advantage of the availability of a high-quality genome assembly for the little brown bat, Myotis lucifugus, to systematically identify and analyze ERVs in this species. We mined an initial set of 362 potentially complete proviruses from the three main classes of ERVs, which were further resolved into 13 major families and 86 subfamilies by phylogenetic analysis. Consensus or representative sequences for each of the 86 subfamilies were then merged to the Repbase collection of known ERV/long terminal repeat (LTR) elements to annotate the retroviral complement of the bat genome. The results show that nearly 5% of the genome assembly is occupied by ERV-derived sequences, a quantity comparable to findings for other eutherian mammals. About one-fourth of these sequences belong to subfamilies newly identified in this study. Using two independent methods, intraelement LTR divergence and analysis of orthologous loci in two other bat species, we found that the vast majority of the potentially complete proviruses identified in M. lucifugus were integrated in the last ∼25 million years. All three major ERV classes include recently integrated proviruses, suggesting that a wide diversity of retroviruses is still circulating in Myotis bats.     

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.     

利用iPBS技术克隆牡丹反转录转座子LTR序列

利用iPBS方法从西北牡丹(Paeonia suffruticosa)品种红绣球和中原牡丹品种洛阳红中扩增出相应片段,经回收、克隆及测序,获得了12条来自牡丹LTR类反转录转座子的LTR序列,并用相关生物信息学软件对序列进行分析。结果表明,这些核苷酸序列表现出较高的异质性,主要表现为缺失突变,序列长度变化范围为313–894 bp,同源性从31.1%–65.8%不等。将其氨基酸序列与已登录的不同植物LTR类反转录转座子LTR氨基酸序列进行聚类分析,结果显示与某些植物相应序列具有较高的同源性,表明可能存在LTR类反转录转座子的横向传递关系。根据克隆出的LTR序列设计SSAP引物,对牡丹29个品种进行了SSAP分子标记分析,结果显示具丰富的多态性。实验验证了用iPBS技术分离牡丹LTR序列的适用性,并为牡丹种质资源评价提供了新的技术手段。     

Evolution of thesrc-related protein tyrosine kinases

A phylogenetic analysis ofsrc-related protein tyrosine kinases (PTKs) showed that one group of these genes is quite ancient in the animals, its divergence predating the divergence of the diploblast and triploblast phyla. Three other major groupings of genes were found to predate the divergence of protostome and deuterostome phyla. Most knownsrc-related PTKs of mammals were found to belong to five well-differentiated families: srcA, srcB, abl, csk, and tec. One srcA gene (fyn) has an alternatively spliced seventh exon which shows a different pattern of relationship from the remainder of the gene; this suggests that this exon may have been derived by a recombinational event with another gene, perhaps one related tofgr. The recently published claim that mammalian members of this family expressed in the nervous system evolve more slowly at nonsynonymous nucleotide sites than do those expressed in the immune system was not supported by an analysis of 13 pairs of human and mouse orthologues. Rather, T-cell-specificsrc-related PTKs were found to have higher rates of nonsynonymous substitution than were those having broader expression. This effect was particularly marked in the peptide binding site of the SH2 domain. While the SH2 binding site was highly conserved among paralogous mammalian members of the srcA and srcB subfamilies, no such effect was seen in the comparison of paralogous members of the csk and tec subfamilies. This suggests that, while the peptide binding function of SH2 is conserved within both srcA and srcB subfamilies, paralogous members of the csk and tec subfamilies have diverged functionally with respect to peptide recognition by SH2.     

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.