Identification of differentially expressed genomic repeats in primary hepatocellular carcinoma and their potential links to biological processes and survival

Hepatocellular carcinoma (HCC) is one of the deadliest cancers. Research on HCC so far primarily focused on genes and provided limited information on genomic repeats, which constitute more than half of the human genome and contribute to genomic stability. In line with this, repeat dysregulation was significantly shown to be pathological in various cancers and other diseases. In this study, we aimed to determine the full repeat expression profile of HCC for the first time. We utilised two independent RNA-seq datasets obtained from primary HCC tumours with matched normal tissues of 20 and 17 HCC patients, respectively. We quantified repeat expressions and analysed their differential expression. We also identified repeats that are cooperatively expressed with genes by constructing a gene coexpression network. Our results indicated that HCC tumours in both datasets harbour 24 differentially expressed repeats and even more elements were coexpressed with genes involved in various metabolic pathways. We discovered that two L1 elements (L1M3b, L1M3de) were downregulated and a handful of HERV subfamily repeats (HERV-Fc1-int, HERV3-int, HERVE_a-int, HERVK11D-int, HERVK14C-int, HERVL18-int) were upregulated with the exception of HERV1_LTRc, which was downregulated. Various LTR elements (LTR32, LTR9, LTR4, LTR52-int, LTR70) and MER elements (MER11C, MER11D, MER57C1, MER9a1, MER74C) were implicated along with few other subtypes including Charlie12, MLT2A2, Tigger15a, Tigger 17b. The only satellite repeat differentially expressed in both datasets was GSATII, whose expression was upregulated in 33 (>90%) out of 37 patients. Notably, GSATII expression correlated with HCC survival genes. Elements discovered here promise future studies to be considered for biomarker and HCC therapy research. The coexpression pattern of the GSATII satellite with HCC survival genes and the fact that it has been upregulated in the vast majority of patients make this repeat particularly stand out for HCC.     

Nucleotide sequence of the rat gamma-crystallin gene region and comparison with an orthologous human region

The sequences of a 51-kb region containing the cluster of five rat gamma-crystallin-coding genes (CRYG) and of a 7-kb region surrounding the sixth rat CRYG gene were determined. Approximately 78% of the total sequence represents intergenic DNA. We also sequenced 22 kb of DNA from the human CRYG gene cluster. All CRYG genes are associated with CpG-rich regions. The sequence similarity between the human and rat gene regions drops sharply (to 65%) in intronic and 3'-flanking regions but decreases only gradually in the 5'-flanking region. Highly conserved regions (greater than 80%) are found as far upstream as 1.5 kb. Overall intergenic distances are conserved. The human region contains much more repetitive DNA (24% vs. 10%) but less simple-sequence (sps) DNA (0.7% vs. 4%) than the rat region. Almost all repeats and spsDNA elements are located in the intergenic region. The location of repetitive and spsDNA differs between the orthologous regions and these elements were probably inserted after the evolutionary separation of rat and man. The Alu repeats in man and the B3 repeats in the rat are close copies of their respective consensus sequences and bordered by virtually perfect repeats. In contrast, the B1 and B2 repeats in the rat have diverged considerably from the consensus sequence and the surrounding direct repeats are usually imperfect. Thus the dispersion of the B1 and B2 repeats in the rat probably preceded that of the B3 repeats. Within the rat genomic region the spacing of Z-DNA elements is surprisingly regular, they are located about 12 kb apart. A search for putative matrix-associated regions suggests that the rat CRYG gene cluster is organized into two chromosomal domains.     

Structure of the human CRFB4 gene: Comparison with its IFNAR neighbor

The cytokine receptor family consists of a growing number of structurally and evolutionarily related transmembrane receptors. CRFB4 and IFNAR are two of the most similar members of this family. They are encoded by two neighboring genes on both human chromosome 21 and murine chromosome 16. The sequence of the human CRFB4 gene was determined from the first exon to the last intron. The nature of the repetitive sequences present in the introns was analyzed and compared with those present in the human IFNAR gene. This analysis leads to considerations of the antiquity of the duplication that gave rise to both genes from a common ancestor. A pseudogene for USF has been identified in the IFNAR gene and a new definition for the repetitive sequence MER37 is proposed. The polymorphism associated with two CA repeats present in the CRFB4 gene is described.The nucleotide sequence reported in this paper has been deposited to GenBank with accession numbers U08988 and U12021 Correspondence to: G. Lutfalla     

The beta globin gene cluster of the prosimian primate Galago crassicaudatus: nucleotide sequence determination of the 41-kb cluster and comparative sequence analyses.

The nucleotide sequence of the beta globin gene cluster of the prosimian Galago crassicaudatus has been determined. A total sequence spanning 41,101 bp contains and links together previously published sequences of the five galago beta-like globin genes (5'-epsilon-gamma-psi eta-delta-beta-3'). A computer-aided search for middle interspersed repetitive sequences identified 10 LINE (L1) elements, including a 5' truncated repeat that is orthologous to the full-length L1 element found in the human epsilon-gamma intergenic region. SINE elements that were identified included one Alu type I repeat, four Alu type II repeats, and two methionine tRNA-derived Monomer (type III) elements. Alu type II and Monomer sequences are unique to the galago genome. Structural analyses of the cluster sequence reveals that it is relatively A+T rich (about 62%) and regions with high G+C content are associated primarily with globin coding regions. Comparative analyses with the beta globin cluster sequences of human, rabbit, and mouse reveal extensive sequence homologies in their genic regions, but only human, galago, and rabbit sequences share extensive intergenic sequence homologies. Divergence analyses of aligned intergenic and flanking sequences from orthologous human, galago, and rabbit sequences show a gradation in the rate of nucleotide sequence evolution along the cluster where sequences 5' of the epsilon globin gene region show the least sequence divergence and sequences just 5' of the beta globin gene region show the greatest sequence divergence.     

MER1, a yeast gene required for chromosome pairing and genetic recombination, is induced in meiosis.

The yeast MER1 gene is required for the production of viable meiotic products and for meiotic recombination. Cytological analysis of chromosome spreads from a mer1 mutant indicates that the MER1 gene product is also required for normal chromosome pairing. mer1 strains make axial elements, precursors to the synaptonemal complex; however, the chromosomes in most nuclei do not become fully synapsed. The DNA sequence of the MER1 coding region was determined; the MER1 open reading frame encodes a 270-amino-acid protein with a molecular mass of 31.1 kilodaltons. The MER1 protein shows limited sequence similarity to calmodulin. Expression of the MER1 gene was examined by RNA blot hybridization analysis and through the construction and analysis of mer1::lacZ fusion genes. Expression of the MER1 gene is meiotically induced and required the IME1 gene product. Thus, expression of the MER1 gene early in meiosis is required for proper chromosome pairing and meiotic recombination.     

Human transposon-like elements insert at a preferred target site: evidence for a retrovirally mediated process.

Members of the human transposon-like family of repetitive sequences (called THE 1 repeats) like many other repetitive DNA sequences are flanked by short direct repeats. Comparison of the base sequences of twelve examples of these flanking direct repeats indicates that THE 1 repeats insert into a preferred genomic target site. In one case, we have identified the sequence of an empty site into which a THE 1 element inserted. The sequence of this empty site and sequences of truncated THE 1 LTRs are consistent with a retroviral mechanism for the insertion of THE 1 elements. Truncated transposon structures illustrate for the first time that intermediate structures of retrotransposition may also be integrated into the genome.     

Association of some potential hormone response elements in human genes with the Alu family repeats.

Short interspersed repeats of the Alu family located in promoters of some human genes contain high-affinity binding sites for thyroid hormone receptor, retinoic acid receptor and estrogen receptor. The standard binding sites for the receptors represent variants of duplicated AGGTCA motif with different spacing and orientation (direct, DR, or inverted, IR), and Alu sequences were found to have functional DR-4, DR-2 or variant IR-3/IR-17 elements. In this study we analyzed distribution and abundance of the elements in a set of human genomic sequences from GenBank and their association with Alu repeats. Our results indicate that a major fraction of potentially active DR-4, DR-2 and variant IR-3/IR-17 elements in the genes is located within Alu repeats. Alu-associated DR-2 elements are conserved in primate evolution. However, very few Alu have potential DR-3 glucocorticoid-response elements. Gel-shift experiments with the probe (AUB) corresponding to the consensus Alu sequence just upstream of the RNA polymerase III promoter B-box and containing duplicated AGGTCA motif indicate that the probe interacts in a sequence-specific manner with human nuclear proteins which bind to standard IR-0, DR-1, DR-4 or DR-5 elements. The AUB sequence was also able to promote thyroid hormone-dependent trans-activation of a reporter gene. The results support the view that Alu retroposons played an important role in evolution of regulation of the primate gene expression by nuclear hormone receptors.     

Identification of a new, abundant superfamily of mammalian LTR-transposons.

A new superfamily of mammalian transposable genetic elements is described with an estimated 40,000 to 100,000 members in both primate and rodent genomes. Sequences known before as MT, ORR-1, MstII, MER15 and MER18 are shown to represent (part of) the long terminal repeats of retrotransposon-like elements related to THE1 in humans. These transposons have structural similarities to retroviruses. However, the putative product of a 1350 base pair open reading frame detected in the consensus internal sequence of THE1 does not resemble retroviral proteins. The elements are named 'Mammalian apparent LTR-retrotransposons' (MaLRs). The internal sequence is usually found to be excised. Their presence in rodents, artiodactyls, lagomorphs, and primates, the divergence of the individual elements from their consensus, and the existence of a probably orthologous element in mouse and man suggest that the first MaLRs were distributed before the radiation of eutherian mammals 80-100 million years ago. MaLRs may prove to be very helpful in determining the evolutionary branching pattern of mammalian orders and suborders.     

Invasion of the human albumin-alpha-fetoprotein gene family by Alu, Kpn, and two novel repetitive DNA elements

The human albumin-alpha-fetoprotein genomic domain contains 13 repetitive DNA elements randomly distributed throughout the symmetrical structures of these genes. These repeated sequences are located at different sites within the two genes. The human albumin gene contains five Alu elements within four of its 14 intervening sequences. Two of these repeats are located in intron 2, and the remaining three are located in introns 7, 8, and 11. The human alpha-fetoprotein gene contains three of these Alu elements, one in intron 4 and the remaining two in the 3'-untranslated region. In addition, the human alpha-fetoprotein gene contains a Kpn repeat and two classes of novel repeats that are absent from the human albumin gene. Six of the Alu elements within the two genes are bound by short direct repeats that harbor five base substitutions in 120 possible positions (60 bp times 2 termini). The absence of Alu repeats from analogous positions in rodents indicates that these repeats invaded the albumin-alpha-fetoprotein domain less than 85 Myr ago (the time of mammalian radiation). Furthermore, considering the conservation of terminal repeats flanking the Alu sequences of the albumin-alpha-fetoprotein domain (0.042 changes per site), we submit that the average time of Alu insertion into this gene family could have been as recently as 15-30 Myr ago.     

Nucleotide sequence analysis of avian retroviruses: structural similarities with transposable elements

Integrated retroviral genomes are flanked by direct repeats of sequences derived from the termini of the viral RNA genome. These sequences are designated long terminal repeats (LTRs). We have determined and analyzed the nucleotide sequence of the LTRs from several exogenous and endogenous avian retroviruses. These LTRs possess several structural similarities with eukaryotic and prokaryotic transposable elements: 1) inverted complementary repeats at the termini, 2) deletions of sequences adjacent to the LTR, 3) small duplications of host sequences flanking the integrated provirus, and 4) sequence homologies with transposable and other genetic elements. These observations suggest that LTRs function in the integration and perhaps transposition of retrovirus genomes. Evidence exists for the presence of a strong promoter sequence within the LTR. The retroviral LTR also contains a "Hogness box" up-stream of the capping site and a poly(A) signal. These features suggest an additional role for the LTR in the regulation of gene expression.     

Long terminal repeat-like elements flank a human immunoglobulin epsilon pseudogene that lacks introns

There are at least three immunoglobulin epsilon genes (C epsilon 1, C epsilon 2, and C epsilon 3) in the human genome. The nucleotide sequences of the expressed epsilon gene (C epsilon 1) and one (C epsilon 3) of the two epsilon pseudogenes were compared. The results show that the C epsilon 3 gene lacks the three intervening sequences entirely and has a 31-base A-rich sequence 16 bases 3' to the putative poly(A) addition signal, indicating that the C epsilon 3 gene is a processed gene. The C epsilon 3 gene sequence is homologous to the five separate DNA segments of the C epsilon 1 gene; namely, a segment in the 5'-flanking region (100 bases) and four exons, which are interrupted by a spacer region or intervening sequences. Long terminal repeat (LTR)-like sequences which contain TATAAA and AATAAA sequences as well as terminal inverted repeats are present in both 5'- and 3'-flanking regions. The 5' and 3' LTR-like sequences do not, however, constitute a direct repeat, unlike transposable elements of eukaryotes and retroviruses. The 3' LTR-like sequence is repetitive in the human genome, but is not homologous to the Alu family DNA. Models for the evolutionary origin of the processed gene flanked by the LTR-like sequences are discussed. The C epsilon 3 gene has a new open frame which codes potentially for an unknown protein of 292 amino acid residues.     

Junctions between genes in the haptoglobin gene cluster of primates.

To investigate the nature of the recombination that generated the haptoglobin three-gene cluster in Old World primates, we sequenced the region between the second gene (HPR) and the third gene (HPP) in chimpanzees (15 kb), as well as the region 3' to the cluster in humans (14 kb). Comparison to the previously sequenced human haptoglobin (HP) and HPR genes showed that the junction point between HP and HPR in humans (junction 1) was not identical to the junction point between the HPR and HPP genes of the chimpanzee (junction 2). An Alu sequence was found at each junction, but both Alu sequences lacked short direct repeats of the flanking genomic DNA. The lack of direct repeats implies that both junction Alu sequences are the products of recombination between different Alu elements. In addition, other insertion and deletion events are clustered in the regions near the junction Alu sequences. The observation that Alu sequences define the junctions between genes in the haptoglobin gene cluster emphasizes the importance of Alu sequences in the evolution of multigene families.     

The human THE-LTR(O) and MstII interspersed repeats are subfamilies of a single widely distributed highly variable repeat family.

Fifteen examples of the transposon-like human element (THE) LTR and thirteen examples of the MstII interspersed repeat are aligned to generate new consensus sequences for these human repetitive elements. The consensus sequences of these elements are very similar, indicating that they compose subfamilies of a single human interspersed repetitive sequence family. Members of this highly polymorphic repeat family have been mapped to at least 11 chromosomes. Seven examples of the THE internal sequence are also aligned to generate a new consensus sequence for this element. Estimates of the abundance of this repetitive sequence family, derived from both hybridization analysis and frequency of occurrence in GenBank, indicate that THE-LTR/MstII sequences are present every 100-3000 kb in human DNA. The widespread occurrence of members of this family makes them useful landmarks, like Alu, L1, and (GT)n repeats, for physical and genetic mapping of human DNA.