Utility of CYP2D6 copy number variants as prognostic biomarker in localized anal squamous cell carcinoma.

BACKGROUND: Anal squamous cell carcinoma (ASCC) is an infrequent tumor whose treatment has not changed since the 1970s. The aim of this study is the identification of biomarkers allowing personalized treatments and improvement of therapeutic outcomes. METHODS: Forty-six paraffin tumor samples from ASCC patients were analyzed by whole-exome sequencing. Copy number variants (CNVs) were identified and their relation to disease-free survival (DFS) was studied and validated in an independent retrospective cohort of 101 ASCC patients from the Multidisciplinary Spanish Digestive Cancer Group (GEMCAD). GEMCAD cohort proteomics allowed assessing the biological features of these tumors. RESULTS: On the discovery cohort, the median age was 61 years old, 50% were males, stages I/II/III: 3 (7%)/16 (35%)/27 (58%), respectively, median DFS was 33 months, and overall survival was 45 months. Twenty-nine genes whose duplication was related to DFS were identified. The most representative was duplications of the CYP2D locus, including CYP2D6, CYP2D7P, and CYP2D8P genes. Patients with CYP2D6 CNV had worse DFS at 5 years than those with two CYP2D6 copies (21% vs. 84%; p < .0002, hazard ratio [HR], 5.8; 95% confidence interval [CI], 2.7-24.9). In the GEMCAD validation cohort, patients with CYP2D6 CNV also had worse DFS at 5 years (56% vs. 87%; p = .02, HR = 3.6; 95% CI, 1.1-5.7). Mitochondria and mitochondrial cell-cycle proteins were overexpressed in patients with CYP2D6 CNV. CONCLUSIONS: Tumor CYP2D6 CNV identified patients with a significantly worse DFS at 5 years among localized ASCC patients treated with 5-fluorouracil, mitomycin C, and radiotherapy. Proteomics pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets for these high-risk patients. PLAIN LANGUAGE SUMMARY: Anal squamous cell carcinoma is an infrequent tumor whose treatment has not been changed since the 1970s. However, disease-free survival in late staged tumors is between 40% and 70%. The presence of an alteration in the number of copies of CYP2D6 gene is a biomarker of worse disease-free survival. The analysis of the proteins in these high-risk patients pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets. Therefore, the determination of the number of copies of CYP2D6 allows the identification of anal squamous carcinoma patients with a high-risk of relapse that could be redirected to a clinical trial. Additionally, this study may be useful to suggest new treatment strategies to increase current therapy efficacy.

Association of Dietary Patterns with MRI Markers of Hepatic Inflammation and Fibrosis in the MAST4HEALTH Study.

Whereas the etiology of non-alcoholic fatty liver disease (NAFLD) is complex, the role of nutrition as a causing and preventive factor is not fully explored. The aim of this study is to associate dietary patterns with magnetic resonance imaging (MRI) parameters in a European population (Greece, Italy, and Serbia) affected by NAFLD. For the first time, iron-corrected T1 (cT1), proton density fat fraction (PDFF), and the liver inflammation fibrosis score (LIF) were examined in relation to diet. A total of 97 obese patients with NAFLD from the MAST4HEALTH study were included in the analysis. A validated semi-quantitative food frequency questionnaire (FFQ) was used to assess the quality of diet and food combinations. Other variables investigated include anthropometric measurements, total type 2 diabetes risk, physical activity level (PAL), and smoking status. Principal component analysis (PCA) was performed to identify dietary patterns. Six dietary patterns were identified, namely "High-Sugar", "Prudent", "Western", "High-Fat and Salt", "Plant-Based", and "Low-Fat Dairy and Poultry". The "Western" pattern was positively associated with cT1 in the unadjusted model (beta: 0.020, p-value: 0.025) and even after adjusting for age, sex, body mass index (BMI), PAL, smoking, the center of the study, and the other five dietary patterns (beta: 0.024, p-value: 0.020). On the contrary, compared with low-intake patients, those with medium intake of the "Low-Fat Dairy and Poultry" pattern were associated with lower values of cT1, PDFF, and LIF. However, patients with a "Low-Fat Dairy and Poultry" dietary pattern were negatively associated with MRI parameters (cT1: beta: -0.052, p-value: 0.046, PDFF: beta: -0.448, p-value: 0.030, LIF: beta: -0.408, p-value: 0.025). Our findings indicate several associations between MRI parameters and dietary patterns in NAFLD patients, highlighting the importance of diet in NAFLD.

Oral microbiome in Proliferative Verrucous Leukoplakia exhibits loss of diversity and enrichment of pathogens.

OBJECTIVES: Oral microbiome plays an important role in oral diseases. Among them, proliferative verrucous leucoplakia (PVL) is an uncommon form of progressive multifocal leukoplakia with a worryingly rate of malignant transformation. Here, we aimed to characterize the oral microbiome of PVL patients and compare it with those of healthy controls. MATERIAL AND METHODS: Oral biopsies from ten PVL patients and five healthy individuals were obtained and used to compare their microbial communities. The sequence of the V3-V4 region of 16S rRNA gene was used as the taxonomic basis to estimate and analyze the composition and diversity of bacterial populations present in the samples. RESULTS: Our results show that the oral microbial composition and diversity are significantly different among PVL patients and healthy donors. The average number of observed operational taxonomic units (OTUs) was higher for healthy donors than for PVL, proving a loss of diversity in PVL. Several OTUs were found to be more abundant in either group. Among those that were significantly enriched in PVL patients, potential protumorigenic pathogens like Oribacterium sp. oral taxon 108, Campylobacter jejuni, uncultured Eubacterium sp., Tannerella, and Porphyromonas were identified. CONCLUSION: Oral microbiome dysbiosis was found in patients suffering from PVL. To the best of our knowledge, this is the first study investigating the oral microbiome alterations in PVL and, due to the limited number of participants, additional studies are needed. Oral microbiota-based biomarkers may be helpful in predicting the risks for the development of PVL.

Nutrigenetic Interactions Might Modulate the Antioxidant and Anti-Inflammatory Status in Mastiha-Supplemented Patients With NAFLD.

Non-alcoholic fatty liver disease (NAFLD) is the most common liver disease with no therapeutic consensus. Oxidation and inflammation are hallmarks in the progression of this complex disease, which also involves interactions between the genetic background and the environment. Mastiha is a natural nutritional supplement known to possess antioxidant and anti-inflammatory properties. This study investigated how a 6-month Mastiha supplementation (2.1 g/day) could impact the antioxidant and inflammatory status of patients with NAFLD, and whether genetic variants significantly mediate these effects. We recruited 98 patients with obesity (BMI ≥ 30 kg/m2) and NAFLD and randomly allocated them to either the Mastiha or the placebo group for 6 months. The anti-oxidative and inflammatory status was assessed at baseline and post-treatment. Genome-wide genetic data was also obtained from all participants, to investigate gene-by-Mastiha interactions. NAFLD patients with severe obesity (BMI > 35kg/m2) taking the Mastiha had significantly higher total antioxidant status (TAS) compared to the corresponding placebo group (P value=0.008). We did not observe any other significant change in the investigated biomarkers as a result of Mastiha supplementation alone. We identified several novel gene-by-Mastiha interaction associations with levels of cytokines and antioxidant biomarkers. Some of the identified genetic loci are implicated in the pathological pathways of NAFLD, including the lanosterol synthase gene (LSS) associated with glutathione peroxidase activity (Gpx) levels, the mitochondrial pyruvate carrier-1 gene (MPC1) and the sphingolipid transporter-1 gene (SPNS1) associated with hemoglobin levels, the transforming growth factor-beta-induced gene (TGFBI) and the micro-RNA 129-1 (MIR129-1) associated with IL-6 and the granzyme B gene (GZMB) associated with IL-10 levels. Within the MAST4HEALTH randomized clinical trial (NCT03135873, www.clinicaltrials.gov) Mastiha supplementation improved the TAS levels among NAFLD patients with severe obesity. We identified several novel genome-wide significant nutrigenetic interactions, influencing the antioxidant and inflammatory status in NAFLD. Clinical Trial Registration: ClinicalTrials.gov, identifier NCT03135873.

Description of the genetic variants identified in a cohort of patients diagnosed with localized anal squamous cell carcinoma and treated with panitumumab.

Squamous cell carcinoma is the most frequent histologic type of anal carcinoma. The standard of care since the 1970s has been a combination of 5-fluorouracil, mitomycin C, and radiotherapy. This treatment is very effective in T1/T2 tumors (achieving complete regression in 80-90% of tumors). However, in T3/T4 tumors, the 3-year relapse free survival rate is only 50%. The VITAL trial aimed to assess the efficacy and safety of panitumumab in combination with this standard treatment. In this study, 27 paraffin-embedded samples from the VITAL trial and 18 samples from patients from daily clinical practice were analyzed by whole-exome sequencing and the influence of the presence of genetic variants in the response to panitumumab was studied. Having a moderate- or high-impact genetic variant in PIK3CA seemed to be related to the response to panitumumab. Furthermore, copy number variants in FGFR3, GRB2 and JAK1 were also related to the response to panitumumab. These genetic alterations have also been studied in the cohort of patients from daily clinical practice (not treated with panitumumab) and they did not have a predictive value. Therefore, in this study, a collection of genetic alterations related to the response with panitumumab was described. These results could be useful for patient stratification in new anti-EGFR clinical trials.

Differentially methylated genes in proliferative verrucous leukoplakia reveal potential malignant biomarkers for oral squamous cell carcinoma.

OBJECTIVES: To explore the pathophysiology of proliferative verrucous leucoplakia (PVL) through a methylated DNA immunoprecipitation and high-throughput sequencing (MeDIP-seq) case-control study. MATERIALS AND METHODS: Oral biopsies from ten PVL patients and five healthy individuals were obtained and used to compare their epigenetic patterns. Network biology methods and integrative analyses of MeDIP-seq and RNAseq data were applied to investigate functional relations among differentially methylated genes (DMGs). The value of selected genes as malignant biomarkers was evaluated in a large cohort of oral squamous cell carcinoma (OSCC) patients from TCGA. RESULTS: A total of 4647 differentially methylated regions were found, with a prominent state of hypermethylation in PVL patients. At the gene level, differentially methylated regions (DMRs) covered 826 genes with distinct roles, including transcription factors and binding proteins with functions in cell adhesion, migration, proliferation, regulation of transcription, bone morphogenesis, and cell signalling. Network analysis revealed three major hubs, two of them collecting proteins related to the response of the patients to PVL and treatment and one hub collecting proteins related to PVL and cancer. The integrative analysis revealed 8 genes (ARTN, CD8A, GATA3, HOXD10, MYO7A, OSR2, PLCB1, and SPOCK2) significantly upregulated in PVL compared to control and 5 genes (ANKRD6, DLG2, GPX3, PITX2, and ZNF736) significantly downregulated. The status of de-regulation found for PVL patients was concordant with what was found for OSCC samples compared to normal adjacent tissue. CONCLUSION: Our findings show the potential of methylation markers in PVL and suggest novel OSCC diagnostic biomarkers which may boost the development of novel epigenetic-based therapies.

Immune expression profile identification in a group of proliferative verrucous leukoplakia patients: a pre-cancer niche for oral squamous cell carcinoma development.

OBJECTIVES: To explore the pathophysiology of proliferative verrucous leukoplakia, a rare oral disorder that exhibits high rates of recurrence and malignant transformation, through a RNAseq case-control study. MATERIAL AND METHODS: We obtained oral biopsies from 10 patients with verrucous leukoplakia lesions and from the mucosa of 5 healthy individuals for sequencing using RNAseq technology. Using bioinformatic methods, we investigated gene expression and enrichment differences between patients both with and without the disorder. We applied network biology methods to investigate functional relations among those genes that were differentially deregulated. RESULTS: We detected 140 differentially expressed genes with distinct roles in immune surveillance, tissue and organ morphogenesis, development, and organization. Of these 140 genes, 111 have been previously described as cancer expression biomarkers, being oral squamous cell carcinoma the most represented type of cancer among them. Of these 140 genes, 26 were prioritized for further investigation as biomarkers using larger sample sizes. CONCLUSIONS: The gene expression patterns of healthy and unhealthy patients differed in 140 genes whose deregulation has a functional impact on normal functioning of the immune system. This immune expression profile provides a plausible hypothesis to explain the transformation to oral squamous cell carcinoma observed in 6 of the 10 assayed cases. CLINICAL RELEVANCE: By determining the molecular bases of the proliferative verrucous leukoplakia disorder and identifying early biomarkers of malignancy, this can allow us to develop new treatment strategies.

Increased miR-7641 Levels in Peritoneal Hyalinizing Vasculopathy in Long-Term Peritoneal Dialysis Patients.

Peritoneal hyalinizing vasculopathy (PHV) represents the cornerstone of long-term peritoneal dialysis (PD), and especially characterizes patients associated with encapsulating peritoneal sclerosis. However, the mechanisms of PHV development remain unknown. A cross sectional study was performed in 100 non-selected peritoneal biopsies of PD patients. Clinical data were collected and lesions were evaluated by immunohistochemistry. In selected biopsies a microRNA (miRNA)-sequencing analysis was performed. Only fifteen patients (15%) showed PHV at different degrees. PHV prevalence was significantly lower among patients using PD fluids containing low glucose degradation products (GDP) (5.9% vs. 24.5%), angiotensin converting enzyme inhibitors (ACEIs) (7.5% vs. 23.4%), statins (6.5% vs. 22.6%) or presenting residual renal function, suggesting the existence of several PHV protective factors. Peritoneal biopsies from PHV samples showed loss of endothelial markers and induction of mesenchymal proteins, associated with collagen IV accumulation and wide reduplication of the basement membrane. Moreover, co-expression of endothelial and mesenchymal markers, as well as TGF-ß1/Smad3 signaling activation were found in PHV biopsies. These findings suggest that an endothelial-to-mesenchymal transition (EndMT) process was taking place. Additionally, significantly higher levels of miR-7641 were observed in severe PHV compared to non-PHV peritoneal biopsies. Peritoneal damage by GDPs induce miRNA deregulation and an EndMT process in submesothelial vessels, which could contribute to collagen IV accumulation and PHV.

Comprehensive Characterization of the Mutational Landscape in Localized Anal Squamous Cell Carcinoma.

Anal squamous cell carcinoma (ASCC) is a rare neoplasm. Chemoradiotherapy is the standard of care, with no therapeutic advances achieved over the past three decades. Thus, a deeper molecular characterization of this disease is still necessary. We analyzed 46 paraffin-embedded tumor samples from patients diagnosed with primary ASCC by exome sequencing. A bioinformatics approach focused in the identification of high-impact genetic variants, which may act as drivers of oncogenesis, was performed. The relation between genetics variants and prognosis was also studied. The list of high-impact genetic variants was unique for each patient. However, the pathways in which these genes are involved are well-known hallmarks of cancer, such as angiogenesis or immune pathways. Additionally, we determined that genetic variants in BRCA2, ZNF750, FAM208B, ZNF599, and ZC3H13 genes are related with poor disease-free survival in ASCC. This may help to stratify the patient's prognosis and open new avenues for potential therapeutic intervention. In conclusion, sequencing of ASCC clinical samples appears an encouraging tool for the molecular portrait of this disease.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification.

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

Analysis of the Transcriptome of the Red Seaweed Grateloupia imbricata with Emphasis on Reproductive Potential.

Grateloupia imbricata is an intertidal marine seaweed and candidate model organism for both industry and academic research, owing to its ability to produce raw materials such as carrageenan. Here we report on the transcriptome of G. imbricata with the aim of providing new insights into the metabolic pathways and other functional pathways related to the reproduction of Grateloupia species. Next-generation sequencing was carried out with subsequent de novo assembly and annotation using state-of-the-art bioinformatic protocols. The results show the presence of transcripts required for the uptake of glycerol, which is a specific carbon source for in vitro culture of G. imbricata and nucleotide sequences that are involved in polyamine-based biosynthesis, polyamine degradation, and metabolism of jasmonates and ethylene. Polyamines, ethylene and methyl jasmonate are plant growth regulators that elicit the development and maturation of cystocarps and the release of spores from seaweeds. Our results will inform studies of the mechanisms that control polysaccharide accumulation, cystocarp formation and spore release. Moreover, our transcriptome information clarifies aspects of red seaweed carposporogenesis with potential benefits for enhancing reproduction.

Differential distribution and enrichment of non-coding RNAs in exosomes from normal and Cancer-associated fibroblasts in colorectal cancer.

Exosome production from cancer-associated fibroblasts seems to be an important driver of tumor progression. We report the first in-depth biotype characterization of ncRNAs, analyzed by Next Generation Sequencing and Bioinformatics, expressed in established primary human normal and cancer-associated fibroblasts (CAFs) from cancer and normal mucosa tissues from 9 colorectal cancer patients, and/or packaged in their derived exosomes. Differential representation and enrichment analyses based on these ncRNAs revealed a significant number of differences between the ncRNA content of exosomes and the expression patterns of the normal and cancer-associated fibroblast cells. ncRNA regulatory elements are specifically packaged in CAF-derived exosomes, supporting a specific cross-talk between CAFs and colon cancer cells and/or other stromal cells, mediated by exosomes. These sncRNAs are potential biomarkers present in cancer-associated fibroblast-derived exosomes, which should thereby contribute to developing new non-invasive diagnostic, prognostic and predictive methods for clinical applications in management of cancer patients.

Genome of wild olive and the evolution of oil biosynthesis.

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Gut Bacterial Community of the Xylophagous Cockroaches Cryptocercus punctulatus and Parasphaeria boleiriana.

Cryptocercus punctulatus and Parasphaeria boleiriana are two distantly related xylophagous and subsocial cockroaches. Cryptocercus is related to termites. Xylophagous cockroaches and termites are excellent model organisms for studying the symbiotic relationship between the insect and their microbiota. In this study, high-throughput 454 pyrosequencing of 16S rRNA was used to investigate the diversity of metagenomic gut communities of C. punctulatus and P. boleiriana, and thereby to identify possible shifts in symbiont allegiances during cockroaches evolution. Our results revealed that the hindgut prokaryotic communities of both xylophagous cockroaches are dominated by members of four Bacteria phyla: Bacteroidetes, Firmicutes, Proteobacteria, and Actinobacteria. Other identified phyla were Spirochaetes, Planctomycetes, candidatus Saccharibacteria (formerly TM7), and Acidobacteria, each of which represented 1-2% of the total population detected. Community similarity based on phylogenetic relatedness by unweighted UniFrac analyses indicated that the composition of the bacterial community in the two species was significantly different (P < 0.05). Phylogenetic analysis based on the characterized clusters of Bacteroidetes, Spirochaetes, and Deltaproteobacteria showed that many OTUs present in both cockroach species clustered with sequences previously described in termites and other cockroaches, but not with those from other animals or environments. These results suggest that, during their evolution, those cockroaches conserved several bacterial communities from the microbiota of a common ancestor. The ecological stability of those microbial communities may imply the important functional role for the survival of the host of providing nutrients in appropriate quantities and balance.

Non-specific symptomatogy presenting in a patient with Gaucher's disease creating a challenging diagnosis: a case report.

Gaucher's disease is the most common of the lysosomal storage diseases; however, with a current worldwide incidence of 1/75,000, it is still a rare occurrence. We present a case of Gaucher's disease type 1 in a Hispanic patient, the first incidence of this specific subtype of Gaucher's disease to be reported in Puerto Rico. Due to the patient's uncharacteristic presentation of negative initial bone marrow biopsy followed by the findings of Gaucher cells on splenic tissue examination, we investigate the specificity and sensitivity of various diagnostic methods to the disease as evidenced in the current literature. Physician awareness of an effective diagnostic approach as illustrated by this report is critical for the early diagnosis of Gaucher's disease so as to prevent disease progression and increased patient morbidity.

The Gypsy Database (GyDB) of mobile genetic elements: release 2.0.

This article introduces the second release of the Gypsy Database of Mobile Genetic Elements (GyDB 2.0): a research project devoted to the evolutionary dynamics of viruses and transposable elements based on their phylogenetic classification (per lineage and protein domain). The Gypsy Database (GyDB) is a long-term project that is continuously progressing, and that owing to the high molecular diversity of mobile elements requires to be completed in several stages. GyDB 2.0 has been powered with a wiki to allow other researchers participate in the project. The current database stage and scope are long terminal repeats (LTR) retroelements and relatives. GyDB 2.0 is an update based on the analysis of Ty3/Gypsy, Retroviridae, Ty1/Copia and Bel/Pao LTR retroelements and the Caulimoviridae pararetroviruses of plants. Among other features, in terms of the aforementioned topics, this update adds: (i) a variety of descriptions and reviews distributed in multiple web pages; (ii) protein-based phylogenies, where phylogenetic levels are assigned to distinct classified elements; (iii) a collection of multiple alignments, lineage-specific hidden Markov models and consensus sequences, called GyDB collection; (iv) updated RefSeq databases and BLAST and HMM servers to facilitate sequence characterization of new LTR retroelement and caulimovirus queries; and (v) a bibliographic server. GyDB 2.0 is available at http://gydb.org.

Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees.

BACKGROUND: Sequencing projects have allowed diverse retroviruses and LTR retrotransposons from different eukaryotic organisms to be characterized. It is known that retroviruses and other retro-transcribing viruses evolve from LTR retrotransposons and that this whole system clusters into five families: Ty3/Gypsy, Retroviridae, Ty1/Copia, Bel/Pao and Caulimoviridae. Phylogenetic analyses usually show that these split into multiple distinct lineages but what is yet to be understood is how deep evolution occurred in this system. RESULTS: We combined phylogenetic and graph analyses to investigate the history of LTR retroelements both as a tree and as a network. We used 268 non-redundant LTR retroelements, many of them introduced for the first time in this work, to elucidate all possible LTR retroelement phylogenetic patterns. These were superimposed over the tree of eukaryotes to investigate the dynamics of the system, at distinct evolutionary times. Next, we investigated phenotypic features such as duplication and variability of amino acid motifs, and several differences in genomic ORF organization. Using this information we characterized eight reticulate evolution markers to construct phenotypic network models. CONCLUSION: The evolutionary history of LTR retroelements can be traced as a time-evolving network that depends on phylogenetic patterns, epigenetic host-factors and phenotypic plasticity. The Ty1/Copia and the Ty3/Gypsy families represent the oldest patterns in this network that we found mimics eukaryotic macroevolution. The emergence of the Bel/Pao, Retroviridae and Caulimoviridae families in this network can be related with distinct inflations of the Ty3/Gypsy family, at distinct evolutionary times. This suggests that Ty3/Gypsy ancestors diversified much more than their Ty1/Copia counterparts, at distinct geological eras. Consistent with the principle of preferential attachment, the connectivities among phenotypic markers, taken as network-represented combinations, are power-law distributed. This evidences an inflationary mode of evolution where the system diversity; 1) expands continuously alternating vertical and gradual processes of phylogenetic divergence with episodes of modular, saltatory and reticulate evolution; 2) is governed by the intrinsic capability of distinct LTR retroelement host-communities to self-organize their phenotypes according to emergent laws characteristic of complex systems. REVIEWERS: This article was reviewed by Eugene V. Koonin, Eric Bapteste, and Enmanuelle Lerat (nominated by King Jordan).

Relationships of gag-pol diversity between Ty3/Gypsy and Retroviridae LTR retroelements and the three kings hypothesis.

BACKGROUND: The origin of vertebrate retroviruses (Retroviridae) is yet to be thoroughly investigated, but due to their similarity and identical gag-pol (and env) genome structure, it is accepted that they evolve from Ty3/Gypsy LTR retroelements the retrotransposons and retroviruses of plants, fungi and animals. These 2 groups of LTR retroelements code for 3 proteins rarely studied due to the high variability - gag polyprotein, protease and GPY/F module. In relation to 3 previously proposed Retroviridae classes I, II and II, investigation of the above proteins conclusively uncovers important insights regarding the ancient history of Ty3/Gypsy and Retroviridae LTR retroelements. RESULTS: We performed a comprehensive study of 120 non-redundant Ty3/Gypsy and Retroviridae LTR retroelements. Phylogenetic reconstruction inferred based on the concatenated analysis of the gag and pol polyproteins shows a robust phylogenetic signal regarding the clustering of OTUs. Evaluation of gag and pol polyproteins separately yields discordant information. While pol signal supports the traditional perspective (2 monophyletic groups), gag polyprotein describes an alternative scenario where each Retroviridae class can be distantly related with one or more Ty3/Gypsy lineages. We investigated more in depth this evidence through comparative analyses performed based on the gag polyprotein, the protease and the GPY/F module. Our results indicate that contrary to the traditional monophyletic view of the origin of vertebrate retroviruses, the Retroviridae class I is a molecular fossil, preserving features that were probably predominant among Ty3/Gypsy ancestors predating the split of plants, fungi and animals. In contrast, classes II and III maintain other phenotypes that emerged more recently during Ty3/Gypsy evolution. CONCLUSION: The 3 Retroviridae classes I, II and III exhibit phenotypic differences that delineate a network never before reported between Ty3/Gypsy and Retroviridae LTR retroelements. This new scenario reveals how the diversity of vertebrate retroviruses is polyphyletically recurrent into the Ty3/Gypsy evolution, i.e. older than previously thought. The simplest hypothesis to explain this finding is that classes I, II and III trace back to at least 3 Ty3/Gypsy ancestors that emerged at different evolutionary times prior to protostomes-deuterostomes divergence. We have called this "the three kings hypothesis" concerning the origin of vertebrate retroviruses.