Identification of S71-related human endogenous retroviral sequences with full-length pol genes

Mechanistic regulation of HERV activation in tumors and implications for translational research in oncology

Transcription of distinct loci of human endogenous retroviruses (HERVs) and in some cases, translation of these transcripts have been consistently observed in many types of cancer. It is still debated whether HERV activation serves as a trigger for carcinogenesis or rather occurs as a consequence of epigenetic alterations and other molecular sequelae that characterize cellular transformation. Here we review the known molecular and epigenetic mechanisms of HERV activation in cancer cells as well as its potential contribution to carcinogenesis. Further, we describe the use of HERV expression in cancer diagnostic and characterize the potential of HERV-derived antigens to serve as novel targets for cancer immunotherapy. We believe this review, which summarizes both what is known as well as unknown in this rapidly developing field, will boost interest in research on the therapeutic potential of targeting HERV elements in tumors and the impact of HERV activation in oncogenesis.

Murine gammaretrovirus group G3 was not found in Swedish patients with myalgic encephalomyelitis/chronic fatigue syndrome and fibromyalgia

Background

The recent report of gammaretroviruses of probable murine origin in humans, called xenotropic murine retrovirus related virus (XMRV) and human murine leukemia virus related virus (HMRV), necessitated a bioinformatic search for this virus in genomes of the mouse and other vertebrates, and by PCR in humans.

Results

Three major groups of murine endogenous gammaretroviruses were identified. The third group encompassed both exogenous and endogenous Murine Leukemia Viruses (MLVs), and most XMRV/HMRV sequences reported from patients suffering from myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). Two sensitive real-time PCRs for this group were developed. The predicted and observed amplification range for these and three published XMRV/HMRV PCRs demonstrated conspicuous differences between some of them, partly explainable by a recombinatorial origin of XMRV. Three reverse transcription real-time PCRs (RTQPCRs), directed against conserved and not overlapping stretches of env, gag and integrase (INT) sequences of XMRV/HMRV were used on human samples. White blood cells from 78 patients suffering from ME/CFS, of which 30 patients also fulfilled the diagnostic criteria for fibromyalgia (ME/CFS/FM) and in 7 patients with fibromyalgia (FM) only, all from the Gothenburg area of Sweden. As controls we analyzed 168 sera from Uppsala blood donors. We controlled for presence and amplifiability of nucleic acid and for mouse DNA contamination. To score as positive, a sample had to react with several of the XMRV/HMRV PCRs. None of the samples gave PCR reactions which fulfilled the positivity criteria.

Conclusions

XMRV/HMRV like proviruses occur in the third murine gammaretrovirus group, characterized here. PCRs developed by us, and others, approximately cover this group, except for the INT RTQPCR, which is rather strictly XMRV specific. Using such PCRs, XMRV/HMRV could not be detected in PBMC and plasma samples from Swedish patients suffering from ME/CFS/FM, and in sera from Swedish blood donors.

Phylogeny-directed search for murine leukemia virus-like retroviruses in vertebrate genomes and in patients suffering from myalgic encephalomyelitis/chronic fatigue syndrome and prostate cancer

Gammaretrovirus-like sequences occur in most vertebrate genomes. Murine Leukemia Virus (MLV) like retroviruses (MLLVs) are a subset, which may be pathogenic and spread cross-species. Retroviruses highly similar to MLLVs (xenotropic murine retrovirus related virus (XMRV) and Human Mouse retrovirus-like RetroViruses (HMRVs)) reported from patients suffering from prostate cancer (PC) and myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) raise the possibility that also humans have been infected. Structurally intact, potentially infectious MLLVs occur in the genomes of some mammals, especially mouse. Mouse MLLVs contain three major groups. One, MERV G3, contained MLVs and XMRV/HMRV. Its presence in mouse DNA, and the abundance of xenotropic MLVs in biologicals, is a source of false positivity. Theoretically, XMRV/HMRV could be one of several MLLV transspecies infections. MLLV pathobiology and diversity indicate optimal strategies for investigating XMRV/HMRV in humans and raise ethical concerns. The alternatives that XMRV/HMRV may give a hard-to-detect “stealth” infection, or that XMRV/HMRV never reached humans, have to be considered.

Chemical induction of endogenous retrovirus particles from the vero cell line of African green monkeys

Endogenous retroviral sequences are present in high copy numbers in the genomes of all species and may be expressed as RNAs; however, the majority are defective for virus production. Although virus has been isolated from various Old World monkey and New World monkey species, there has been no report of endogenous retroviruses produced from African green monkey (AGM) tissues or cell lines. We have recently developed a stepwise approach for evaluating the presence of latent viruses by chemical induction (Khan et al., Biologicals 37:196-201, 2009). Based upon this strategy, optimum conditions were determined for investigating the presence of inducible, endogenous retroviruses in the AGM-derived Vero cell line. Low-level reverse transcriptase activity was produced with 5-azacytidine (AzaC) and with 5'-iodo-2'-deoxyuridine (IUdR); none was detected with sodium butyrate. Nucleotide sequence analysis of PCR-amplified fragments from the gag, pol, and env regions of RNAs, prepared from ultracentrifuged pellets of filtered supernatants, indicated that endogenous retrovirus particles related to simian endogenous type D betaretrovirus (SERV) sequences and baboon endogenous virus type C gammaretrovirus (BaEV) sequences were induced by AzaC, whereas SERV sequences were also induced by IUdR. Additionally, sequence heterogeneity was seen in the RNAs of SERV- and BaEV-related particles. Infectivity analysis of drug-treated AGM Vero cells showed no virus replication in cell lines known to be susceptible to type D simian retroviruses (SRVs) and to BaEV. The results indicated that multiple, inducible endogenous retrovirus loci are present in the AGM genome that can encode noninfectious, viruslike particles.

(Epi)genomics and neurodevelopment in schizophrenia: monozygotic twins discordant for schizophrenia augment the search for disease-related (epi)genomic alterations

Unlike stunning breakthroughs in the identification of genes for Mendelian disorders during the last three decades, gene identification in most complex disorders has been full of twists and turns and little progress. Doing more of the same will not guarantee success. The lessons learned argue for a need to reconsider genetic models that are appropriate for the disorder in question along with an interdisciplinary, systematic approach using genomic methods that have now become possible. We will use schizophrenia as an example to review the genetic progress to date that has been disappointing. We will argue that the causation of this complex disease may involve heterogeneous genomic changes of major effect. We will provide three approaches, retroviral transpositions, methylation, and copy number variations, to test this hypothesis. We will present arguments to suggest that such experiments will be most effective if undertaken on monozygotic twins. It will include our experience with associated experiments on the monozygotic twins discordant for schizophrenia. The results support that (epi)genomic changes of major effect, rather than accumulation of mutations of small effect, underlie the causation of this complex disease. More important, this experimental strategy will be an effective strategy for studies on other complex (behavioural) disorders as well.

Expression and phylogenetic analyses of human endogenous retrovirus HC2 belonging to the HERV-T family in human tissues and cancer cells

Recently, a new HERV-T family, representative of the HERV-S71 and HERV-HC2 family, was identified using a screen for envelope genes and a computer-assisted database search. Here, we investigate expression of pol fragments of HERV-HC2 belonging to the HERV-T family in various human tissues and cancer cells. The pol gene was expressed in nearly all human tissues examined and in all cancer cell lines. Expression analyses suggest that the pol gene of HERV-HC2 family is more actively transcribed in human cancer cells than in normal tissues, suggesting a functional role during carcinogenesis. Phylogenetic analysis of the HERV-HC2 pol family revealed three groups (I, II, and III) generated through evolutionary divergence during primate evolution, indicating that they were integrated into primate genomes approximately 56 million years (MY) ago and have evolved at a rate of 0.2% nucleotide differences per MY. Our data might contribute to an understanding of the information on the transcriptional and pathological potential of the HERV-T family in human disease, including cancer.

The distribution of pol containing human endogenous retroviruses in non-human primates

Few human endogenous retroviruses (HERVs) have been extensively studied in non-human primates. Such investigations have demonstrated that several element classes are primate unique, contain members with important biological function, are conserved in specific primate lineages, and have in some cases expanded in copy number. We have examined multiple sub-families of all major groups of HERVs using a DNA microarray based on the reverse transcriptase (RT) domain of the viral polymerase gene (pol). The microarray was used to investigate the distribution of HERVs in non-human primates with particular focus on the differences between New World monkeys (NWMs) and other anthropoids. This is the first study examining most HERV families in multiple non-human primate DNAs using a uniform and sensitive method and suggests that major differences exist between primate groups. The results indicate that a major invasion and expansion of pol containing HERVs occurred after the platyrrhine (NWM) lineage separated from the catarrhines (Old World Monkeys and apes).

Role of endogenous retroviruses in autoimmune diseases

Endogenous retroviruses (ERVs) correspond to the integrated proviral form of infectious retroviruses that are trapped within the genome by mutations. Endogenous retroviruses represent a key molecular link between the host genome and infectious viral particles. Proteins encoded by ERVs are recognized by antiviral immune responses and become targets of autoreactivity. Activation of ERVs, such as human ERV-K or a human T-cell lymphotropic virus-related endogenous sequence, may also mediate pathogenicity of Epstein-Barr virus. Endogenous retrovirus peptides can directly regulate immune responses. Thus, molecular mimicry and immunomodulation by ERVs may account for self-reactivity and abnormal T- and B-cell functions in autoimmune disorders.

Molecular characterization of a 2.7-kb, 12q13-specific, retroviral-related sequence isolated by RDA from monozygotic twin pairs discordant for schizophrenia

This report deals with the molecular characterization of a representational difference analysis (RDA)-derived sequence (SZRV-2, GenBank accession No. AF135486; Genome Database accession Nos. 7692183 and 7501402) from three monozygotic twin pairs discordant for schizophrenia (MZD). The results suggest that it is a primate-specific, heavily methylated, and placentally expressed (-7-kb mRNA) endogenous retroviral-related (ERV) sequence of the human genome. We have mapped this sequence to 12q13 using two SZRV-2 positive BAC clones (4K11 (Genome Survey Sequence Database No. 1752076; GenBank accession No. AZ301773) and 501H16) by fluorescence in situ hybridization. End sequencing of the 4K11 BAC clone has allowed identification of nearby genes from the human genome database at NCBI that may be of interest in schizophrenia research. These include viral-related sequences (potential hot spots for insertions), developmental, channel, and signal transduction genes, as well as genes affecting expression of certain receptors in neurons. Furthermore, when used as a probe on Southern blots, SZRV-2 detected no difference between schizophrenia patients from southwestern Ontario and their matched controls. However, it identified aberrant methylation in one of the eight patients and none of the 21 unaffected controls. Although additional experiments will be required to establish the significance, if any, of SZRV-2 methylation in the complex etiology of schizophrenia, molecular results included offer a novel insight into the role of retroviral-related sequences in the origin, organization, and regulation of the human genome.

Isolation and phylogeny of new endogenous retroviral sequences belonging to the HERV-F family

A new human endogenous retroviral family (HERV-F) has been identified from human chromosome 7q31.1-q31.3 that was identical to the XA34 cDNA clone isolated from a human glioma cDNA library with an ERV-9 env probe. We investigated pol gene sequences of the HERV-F family from a human monochromosomal DNA panel and analyzed these with HERV-F. The pol gene sequences of the HERV-F family were detected on chromosomes 3, 6, 7, 10, 11, 14, 19, 20, X, and Y as examined by PCR. Thirty-six pol gene sequences identified from the human chromosomes have a high degree of sequence similarity (80-99%) with that of the HERV-F. Phylogenetic analysis of pol gene sequences distinctively showed four groups, indicating that the HERV-F family could be amplified at least four times after the original integration into the human genome or represent integration events separately during hominid evolution. One clone (HFY-3) on chromosome Y shared 100% sequence identity with a clone (HF19-2) on chromosome 19, and a clone (HF20-6) on chromosome 20 suggests either a recent retrotransposition or a chromosomal translocation. The history of endogenous retroviral sequences may contribute to an understanding of evolutionary change in human genomes.

Retrotransposable L1 elements expressed in rheumatoid arthritis synovial tissue: association with genomic DNA hypomethylation and influence on gene expression

OBJECTIVE

Rheumatoid arthritis (RA) is characterized by a progressive destruction of joints by invasive synovial fibroblasts (SF). We searched for retroviral sequences in RA synovial fluid pellets, identified a sequence similar to that of open reading frame 2 (ORF2)/L1 retrotransposable elements, explored the expression of L1 in RA synovial tissues and cultured RA SF, and investigated the link to genomic DNA hypomethylation and the influence of functional L1 on gene expression.

METHODS

RA synovial fluid pellets were screened by reverse transcriptase-polymerase chain reaction (RT-PCR) using degenerated pol primers. The sequences were identified by GenBank search. Riboprobes to ORF2/L1 and galectin-3 and antibodies to the ORF1/L1-related p40 protein were used for in situ hybridization and immunohistochemistry of synovial tissues and cultured RA SF. Real-time quantitative RT-PCR was used for detecting ORF1 messenger RNA (mRNA). Since DNA hypomethylation occurs in inflammatory diseases, we incubated cells with the methylation inhibitor 5-aza-2'-deoxycytidine (5-azaC) and compared RA SF and osteoarthritis (OA) SF. L1-negative RA SF were transfected with the functional L1.2 construct, and differential gene expression was analyzed by subtractive hybridization combined with nested PCR.

RESULTS

RNA sequences similar to those of ORF2/L1 retrotransposable elements, THE1 transposon, human endogenous retrovirus (ERV)-E, human ERV-HC2, and gibbon ape leukemia virus pol genes were isolated from different RA synovial fluid pellets. In RA synovial tissues, ORF2/L1 transcripts were detected in the sublining layer and at sites of cartilage and bone destruction. Galectin-3 mRNA and L1-related ORF1/ p40 protein showed similar expression patterns. In contrast, OA synovial tissues in situ and cultures in vitro were negative. Real-time quantitative RT-PCR confirmed the presence of ORF1 mRNA in cultured RA SF (30-300-fold the amount in normal SF), demonstrating the existence of a nondegenerated and functional L1 element. In vitro, the majority of RA SF expressed ORF2/L1 mRNA. After incubation of SF with 5-azaC, L1 mRNA appeared in a time- and dose-dependent manner. Compared with OA SF, RA SF were more sensitive to 5-azaC. After transfection of RA SF with a functional L1.2 element, human stress-activated protein kinase 2 delta (SAPK2delta [or SAPK4]), met protooncogene, and galectin-3 binding protein genes were differentially expressed. The transcription of the SAPK2delta gene, favored also by DNA hypomethylation in vitro, was confirmed in RA synovial tissues.

CONCLUSION

Taken together, these data suggest that L1 elements and SAPK2delta pathways play a role in the activation of RA SF.

HERV-IP-T47D, a novel type C-related human endogenous retroviral sequence derived from T47D particles

A new type C retrovirus-related endogenous pol sequence (ERV-FTD) found to be occasionally copackaged in retrovirus-like particles released by the human mammary carcinoma cell line T47D was used to screen a human genomic library (Seifarth W, Skladny H, Krieg-Schneider F, Reichert A, Hehlmann R, and Leib-Mösch C: J Virol 1995;69:6408-6416). The DNA sequence of one full-length clone now reveals a human endogenous proviral sequence (HERV) of 4190 bp in length comprising a 5' LTR (489 bp) and regions with 37 and 74% overall amino acid homology to RTVL-Ia gag and pol genes, respectively. About 35 related elements were found to be distributed on all human chromosomes except 16, 17, and Y. Sequence comparisons with Mo-MuLV and various type C-related HERVs suggest that despite a proline primer-binding site this novel HERV element, now named HERV-IP-T47D, can be assigned to one family together with known HERV-I elements. Phylogenetic analyses of 5 proviral and 25 solitary LTR sequences confirmed the existence of two distinct but closely related subgroups of the HERV-IP superfamily in the primate genome. In contrast to most known HERV-families, the evolutionary age of HERV-IP elements dates back prior to the divergence of New and Old World monkeys. Despite their old age, members of the HERV-IP family are still transcriptionally active and were found to be highly expressed in specific human tissues such as liver and kidney.

Molecular characterization of a MSRV-like sequence identified by RDA from monozygotic twin pairs discordant for schizophrenia

Retroviral-related amplicons were used in modified RDA to identify four sequences from affected members of three pairs of monozygotic twins discordant for schizophrenia. One sequence (schizophrenia associated retrovirus, SZRV-1, GenBank Accession No. AF135487) is characterized here. It is similar to two known sequences of retroviral origin: multiple sclerosis-associated retrovirus, MSRV (GenBank Accession No. AF009668), and ERV-9 (GenBank Accession No. S77575). It is present in multiple copies in the human genome and has been localized to six different chromosomal sites. A zooblot shows that this multicopy sequence is predominant in the primate lineage and present in rhesus monkeys and humans. SZRV-1 is expressed as a 9-kb RNA band in the placenta. This could offer support to the hypothesis that retroviral sequences transposing during fetal growth may alter neurodevelopmental genes and cause diseases, although its direct involvement in the causation of schizophrenia remains to be established.

Isolation and phylogeny of endogenous retrovirus sequences belonging to the HERV-W family in primates

An investigation was undertaken of primate pol gene sequences from a novel endogenous retrovirus family, ERV-W, related to a new human endogenous retrovirus family (HERV-W) that includes multiple sclerosis-associated retrovirus (MSRV) sequences identified in particles recovered from monocyte cultures from patients with multiple sclerosis. The pol gene sequences of the ERV-W family were detected in hominoids and Old World monkeys, but not in New World monkeys, whereas ERV-W long terminal repeat-like elements were detected in all primates (hominoids, Old World monkeys and New World monkeys). Thirty-two pol gene sequences from hominoids and Old World monkeys showed a high degree of sequence identity to MSRV and other HERV-W sequences. Phylogenetic analysis indicated close relationships of pol gene sequences across primate species. The analysis suggests that the ERV-W family has evolved independently but in constrained patterns ('parallel evolution') in different primate species, including man. The ratio of synonymous to non-synonymous substitutions indicated that negative selective pressure is acting on CHW1-1 from chimpanzee, HBW6-6 from baboon and HWX5 from man, sequences that have no disruption by point mutation or insertions/deletions. Therefore, these pol gene sequences could be associated with an active provirus in primates. The findings indicate that the ERV-W family has continued to evolve in the course of the primate radiation and may include members with a capacity to influence gene function and possibly cause disease.

Characterization of human MMTV-like (HML) elements similar to a sequence that was highly expressed in a human breast cancer: further definition of the HML-6 group

Previously, we found a retroviral sequence, HML-6.2BC1, to be expressed at high levels in a multifocal ductal breast cancer from a 41-year-old woman who also developed ovarian carcinoma. The sequence of a human genomic clone (HML-6.28) selected by high-stringency hybridization with HML-6.2BC1 is reported here. It was 99% identical to HML-6.2BC1 and gave the same restriction fragments as total DNA. HML-6.28 is a 4.7-kb provirus with a 5'LTR, truncated in RT. Data from two similar genomic clones and sequences found in GenBank are also reported. Overlaps between them gave a rather complete picture of the HML-6.2BC1-like human endogenous retroviral elements. Work with somatic cell hybrids and FISH localized HML-6.28 to chromosome 6, band p21, close to the MHC region. The causal role of HML-6.28 in breast cancer remains unclear. Nevertheless, the ca. 20 Myr old HML-6 sequences enabled the definition of common and unique features of type A, B, and D (ABD) retroviruses. In Gag, HML-6 has no intervening sequences between matrix and capsid proteins, unlike extant exogenous ABD viruses, possibly an ancestral feature. Alignment of the dUTPase showed it to be present in all ABD viruses, but gave a phylogenetic tree different from trees made from other ABD genes, indicating a distinct phylogeny of dUTPase. A conserved 24-mer sequence in the amino terminus of some ABD envelope genes suggested a conserved function.

Comparative analyses of LTRs of the ERV-H family of primate-specific retrovirus-like elements isolated from marmoset, African green monkey, and man

We have isolated 8 different long terminal repeat (LTR) sequences of the ERV-H family of endogenous retrovirus-like elements from human chromosome 18, 9 from African green monkey, and 28 from marmoset. Human ERV-H LTRs have been divided into three types designated Type I, Type Ia, and Type II. Comparative analyses of the 45 isolated LTRs and 60 human ERV-H LTRs enabled a further subdivision into 13 subtypes. Type I elements were widely distributed in all three species. Their average evolutionary age (40 MYr), estimated by a consensus sequence approach, suggests that they first expanded in the genomes at the time New- and Old World monkeys diverged. The occurence of some very old Type I sequences indicate that ERV-H elements may have integrated even before prosimians and primates diverged. Type Ia and - II elements were found in both monkey species. Promoter active Type I and Type Ia LTRs were found while Type II LTRs were inactive. Promoter active Type I LTRs generally contained a functional GC/GT box immediately 3' to the TATA box, providing strong binding of Sp1 family proteins, while the highly promoter active Type Ia element H6 contained synergistically acting Sp1 binding sites located in the U3 enhancer region.

Identification of endogenous retroviral sequences based on modular organization: proviral structure at the SSAV1 locus

The current genome sequencing projects reveal megabases of unknown genomic sequences. About 1% of these sequences can be expected to be of retroviral origin. These are often severely deleted or mutated. Therefore, identification of the retroviral origin of these sequences can be very difficult due to the absence of convincing overall sequence similarity. There are also many copies of solo-LTRs (long terminal repeats) distributed throughout genomic sequences. LTR and envelope sequences in general are among the most divergent parts of the retroviral genome and thus especially hard to detect in mutated endogenous sequences. We took advantage of the fact that these retroviral sections contain short highly conserved sequence regions providing retroviral hallmarks even after loss of overall similarity. We defined several sequence elements and peptide motifs within LTR and Env sequences and used these elements to construct models for LTRs and Env proteins of mammalian C-type retroviruses. We then used this strategy to identify successfully the hitherto missing LTRs and an env-like region in the S71 human retroviral sequence. Our approach provides a new strategy for identifying remotely related retroviral sequences in genomic DNA (especially human DNA), of potential significance for the interpretation of genomic sequences obtained from the current large-scale sequencing projects.

Retrovirus-like particles released from the human breast cancer cell line T47-D display type B- and C-related endogenous retroviral sequences

The human mammary carcinoma cell line T47-D releases retrovirus-like particles of type B morphology in a steroid-dependent manner (I. Keydar, T. Ohno, R. Nayak, R. Sweet, F. Simoni, F. Weiss, S. Karby, R. Mesa-Tejada, and S. Spiegelman, Proc. Natl. Acad. Sci. USA 81:4188-4192, 1984). Furthermore, reverse transcriptase (RT) activity is found to be associated with particle preparations. Using a set of degenerate primers derived from a conserved region of retroviral pol genes, we repeatedly amplified three different retroviral sequences (MLN, FRD, and FTD) from purified T47-D particles in several RT-PCR experiments. Screening of a human genomic library and Southern blot analysis revealed that these sequences are of endogenous origin. ERV-MLN represents a multicopy family of human endogenous retroviral elements (HERVs) with two closely related copies and up to 20 more distantly related members. In contrast, ERV-FRD and ERV-FTD comprise only one copy and five to seven related elements per haploid human genome. DNA sequence analysis of the proviral pol region of ERV-MLN revealed an uninterrupted stretch of 241 amino acids that shows 65% identity with the RT of the type B-related HERV designated HERV-K10. ERV-FRD and ERV-FTD are defective type C-related HERVs. The pol gene of ERV-FRD displays a nucleotide homology of 54% to the gibbon ape leukemia virus, and the pol gene of ERV-FTD is about 67% homologous to members of the RTVL-I family of HERVs. Our results thus indicate that the retroviral particles released by the breast cancer cell line T47-D are probably generated by complementation of several endogenous proviruses and can package retroviral transcripts of different origins.

Evolution and biological significance of human retroelements

Retroelements comprise a substantial portion of the human genome. Their large number and ubiquitous distribution has led scientists to speculate about their evolutionary origin and their biological functions. Human endogenous retroviruses and their retrotransposon relatives represent a reservoir of possibly pathogenic retroviral genes that may be activated spontaneously or by environmental conditions. They can act as insertion mutagens and activate or inactivate cellular genes, or may be involved in chromosome aberrations by recombination of related elements on different chromosomal locations. Retroviral gene products themselves may also be pathogenic and, for example, could be implicated in the development of tumors and autoimmune diseases. On the other hand, endogenous retroviral elements and nonviral retroposons are thought to have played an important role in shaping the genomes of vertebrates by intracellular transposition events and by generating hot spots of recombination. In the course of time, some of these elements have acquired cellular functions, such as, for instance, in the regulation of gene expression. Therefore, the role of human endogenous retroviruses and retroposons in biological processes is currently a subject of great interest.