External validation of models to predict hepatocellular carcinoma in Hepatitis C Virus cured F3-F4 patients.

BACKGROUND & AIMS: Several hepatocellular carcinoma (HCC) risk-models have been developed to individualise patient surveillance following sustained viral response (SVR) in Hepatitis C Virus patients. Validation of these models in different cohorts is an important step to incorporate a more personalised risk assessment in clinical practice. We aimed at applying these models to stratify the risk in our patients and potentially determine cost-saving associated with individualised HCC risk-stratification screening strategy. METHODS: Patients with baseline F3-4 fibrosis treated with new oral direct-acting antivirals who had reached a SVR were regularly followed as part of the HCC surveillance strategy. Six models were applied: Pons, aMAP, Ioannou, HCC risk, Alonso and Semmler. Validation of the models was performed based on sensitivity and the proportion of patients labelled as "high risk". RESULTS: After excluding 557 with less than 3 fibrosis, 12 without SVR, 18 with a follow up (FU) <1 year, 17 transplant recipients, 16 lost to FU and 31 with HCC at time of antiviral therapy, our cohort consisted of 349 F3-4 SVR patients. Twenty-three patients (6.6%) developed HCC after a median FU of 5.12 years. The sensitivity of the different models varied between 0.17 (Semmler7noalcohol) and 1 (Alonso A and aMAP). The lowest proportion of high-risk patients corresponded to the Semmler-noalcohol model (5%). Sixty-three and 90% of the Alonso A and aMAP patients, respectively were labelled as high risk. The most reliable HCC risk-model applied to our cohort to predict HCC development is the Alonso model (based on fibrosis stage assessed by liver stiffness measurements or Fibrosis-4 index (FIB-4) at baseline and after 1 year, and albumin levels at 1 year) with a-100% sensitivity in detecting HCC among those at high risk and 63% labelled as high risk. The application of the model would have saved the cost of 1290 ultrasound no longer being performed in the 37% low-risk group. CONCLUSION: In our cohort, the Alonso A model allows the most reliable reduction in HCC screening resulting in safely stopping life-long monitoring in about a third of F3-F4 patients achieving SVR with DAAs.

Tannic Acid and Ethyl Gallate Potentialize Paclitaxel Effect on Microtubule Dynamics in Hep3B Cells.

Among broad-spectrum anticancer agents, paclitaxel (PTX) has proven to be one of the most effective against solid tumors for which more specific treatments are lacking. However, drawbacks such as neurotoxicity and the development of resistance reduce its therapeutic efficacy. Therefore, there is a need for compounds able to improve its activity by synergizing with it or potentiating its effect, thus reducing the doses required. We investigated the interaction between PTX and tannins, other compounds with anticancer activity known to act as repressors of several proteins involved in oncological pathways. We found that both tannic acid (TA) and ethyl gallate (EG) strongly potentiate the toxicity of PTX in Hep3B cells, suggesting their utility in combination therapy. We also found that AT and EG promote tubulin polymerization and enhance the effect of PTX on tubulin, suggesting a direct interaction with tubulin. Biochemical experiments confirmed that TA, but not EG, binds tubulin and potentiates the apparent binding affinity of PTX for the tubulin binding site. Furthermore, the molecular docking of TA to tubulin suggests that TA can bind to two different sites on tubulin, one at the PTX site and the second at the interface of α and ß-tubulin (cluster 2). The binding of TA to cluster 2 could explain the overstabilization in the tubulin + PTX combinatorial assay. Finally, we found that EG can inhibit PTX-induced expression of pAkt and pERK defensive protein kinases, which are involved in resistance to PXT, by limiting cell death (apoptosis) and favoring cell proliferation and cell cycle progression. Our results support that tannic acid and ethyl gallate are potential chemotherapeutic agents due to their potentiating effect on paclitaxel.

Rapid genome functional annotation pipeline anchored to the house sparrow (Passer domesticus, Linnaeus 1758) genome reannotation.

The house sparrow (Passer domesticus) is a valuable avian model for studying evolutionary genetics, development, neurobiology, physiology, behavior, and ecology, both in laboratory and field-based settings. The current annotation of the P. domesticus genome available at the Ensembl Rapid Release site is primarily focused on gene set building and lacks functional information. In this study, we present the first comprehensive functional reannotation of the P. domesticus genome using intestinal Illumina RNA sequencing (RNA-Seq) libraries. Our revised annotation provides an expanded view of the genome, encompassing 38592 transcripts compared to the current 23574 transcripts in Ensembl. We also predicted 14717 protein-coding genes, achieving 96.4% completeness for Passeriformes lineage BUSCOs. A substantial improvement in this reannotation is the accurate delineation of untranslated region (UTR) sequences. We identified 82.7% and 93.8% of the transcripts containing 5'- and 3'-UTRs, respectively. These UTR annotations are crucial for understanding post-transcriptional regulatory processes. Our findings underscore the advantages of incorporating additional specific RNA-Seq data into genome annotation, particularly when leveraging fast and efficient data processing capabilities. This functional reannotation enhances our understanding of the P. domesticus genome, providing valuable resources for future investigations in various research fields.

Liquidhope: methylome and genomic profiling from very limited quantities of plasma-derived DNA.

Analysis of the methylome of tumor cell-free deoxyribonucleic acid (DNA; cfDNA) has emerged as a powerful non-invasive technique for cancer subtyping and prognosis. However, its application is frequently hampered by the quality and total cfDNA yield. Here, we demonstrate the feasibility of very low-input cfDNA for whole-methylome and copy-number profiling studies using enzymatic conversion of unmethylated cysteines [enzymatic methyl-seq (EM-seq)] to better preserve DNA integrity. We created a model for predicting genomic subtyping and prognosis with high accuracy. We validated our tool by comparing whole-genome CpG sequencing with in situ cohorts generated with bisulfite conversion and array hybridization, demonstrating that, despite the different techniques and sample origins, information on cfDNA methylation is comparable with in situ cohorts. Our findings support use of liquid biopsy followed by EM-seq to assess methylome of cancer patients, enabling validation in external cohorts. This advance is particularly relevant for rare cancers like neuroblastomas where liquid-biopsy volume is restricted by ethical regulations in pediatric patients.

Online polarization and cross-fertilization in multi-cleavage societies: the case of Spain.

The impact of the social media (SM) has been seen on the one hand as the cause of large exacerbation of negative messages, responsible for massively harmful societal phenomenon against democracies. On the other hand, recent studies have begun to look at how these online channels were able to provide a new impulse in human communication. The novelty of our work resides on analysing several axes of polarizations related to different societal topics. We believe our approach to reflect a more complex society, differing from the recent literature, which has considered a unique left-right dichotomic cleavage. Our methodology consists of extracting topics from the priority themes of the SM debate, using BERT language processing techniques and TF-IDF model. Our results show situation of social media interactions in a multidimensional space does exist. We highlight how social media behaviours, polarization and cross-fertilization differ as upon concrete topics. We argue therefore the 'mega-identity partisanship' which differentiate US online users in two different spaces cannot be extended for the rest of countries taking as first evidence the case of Spain. Further research should extend our conclusions for a possible generalization.

Opportunities Lost? Evolutionary Causes and Ecological Consequences of the Absence of Trehalose Digestion in Birds.

AbstractTrehalose is a nonreducing disaccharide that is a primary storage and energy source in prokaryotes, yeasts, fungi, and invertebrates. Vertebrates digest trehalose with the intestinal brush border membrane (BBM) enzyme trehalase. Intestinal trehalase activity is reported to be either very low or absent in several bird species. We assayed trehalase activity in 19 avian species, used proteomic analysis to quantify its abundance in the intestinal BBM, and used analyses of available genomes to detect the presence of the gene that codes for trehalase (Treh). We found no intestinal trehalase activity in birds, trehalase was absent from the proteome of their intestinal BBM, and the gene coding for trehalase was absent in their genomes. Surveys of available transcriptomes support the hypothesis that Treh is absent in birds. The trehalase gene was found in the same conserved syntenic block within the genome of all vertebrates surveyed except birds. Our analysis suggests that Treh was lost in an inversion followed by a reinsertion of a large gene block. This event appears to have taken place after the split between crocodiles and birds and dinosaurs. Birds are unable to digest a common dietary sugar like trehalose because their ancestor lost the trehalase gene. The loss of this gene seems to represent an ecological cost, as insectivorous birds seem to be unable to digest a carbohydrate present in their prey. We also speculate that the paucity of mycophagy in birds is due to the presence of large amounts of this sugar in fungal tissues.

In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Chromatin topology defines estradiol-primed progesterone receptor and PAX2 binding in endometrial cancer cells.

Estrogen (E2) and Progesterone (Pg), via their specific receptors (ERalpha and PR), are major determinants in the development and progression of endometrial carcinomas, However, their precise mechanism of action and the role of other transcription factors involved are not entirely clear. Using Ishikawa endometrial cancer cells, we report that E2 treatment exposes a set of progestin-dependent PR binding sites which include both E2 and progestin target genes. ChIP-seq results from hormone-treated cells revealed a non-random distribution of PAX2 binding in the vicinity of these estrogen-promoted PR sites. Altered expression of hormone regulated genes in PAX2 knockdown cells suggests a role for PAX2 in fine-tuning ERalpha and PR interplay in transcriptional regulation. Analysis of long-range interactions by Hi-C coupled with ATAC-seq data showed that these regions, that we call 'progestin control regions' (PgCRs), exhibited an open chromatin state even before hormone exposure and were non-randomly associated with regulated genes. Nearly 20% of genes potentially influenced by PgCRs were found to be altered during progression of endometrial cancer. Our findings suggest that endometrial response to progestins in differentiated endometrial tumor cells results in part from binding of PR together with PAX2 to accessible chromatin regions. What maintains these regions open remains to be studied.

A set of accessible enhancers enables the initial response of breast cancer cells to physiological progestin concentrations.

Here, we report that in T47D breast cancer cells 50 pM progestin is sufficient to activate cell cycle entry and the progesterone gene expression program. At this concentration, equivalent to the progesterone blood levels found around the menopause, progesterone receptor (PR) binds only to 2800 genomic sites, which are accessible to ATAC cleavage prior to hormone exposure. These highly accessible sites (HAs) are surrounded by well-organized nucleosomes and exhibit breast enhancer features, including estrogen receptor alpha (ERα), higher FOXA1 and BRD4 (bromodomain containing 4) occupancy. Although HAs are enriched in RAD21 and CTCF, PR binding is the driving force for the most robust interactions with hormone-regulated genes. HAs show higher frequency of 3D contacts among themselves than with other PR binding sites, indicating colocalization in similar compartments. Gene regulation via HAs is independent of classical coregulators and ATP-activated remodelers, relying mainly on MAP kinase activation that enables PR nuclear engagement. HAs are also preferentially occupied by PR and ERα in breast cancer xenografts derived from MCF-7 cells as well as from patients, indicating their potential usefulness as targets for therapeutic intervention.

Monitoring Sexually Transmitted Infections in Cervicovaginal Exfoliative Samples in Mexican Women.

BACKGROUND: Globally, Sexually Transmitted Infections (STIs) are a major cause of morbidity in sexually active individuals, having complications in reproduction health and quality of life. In concordance with the Sustainable Development Goals (SDG), the study aimed to investigate the prevalence of Candida spp., Ureaplasma spp., Trichomonas vaginalis, Neisseria gonorrhoeae, Chlamydia trachomatis, HSV, and Mycoplasma spp. from cervicovaginal samples and to correlate them with the gynecological history of the patients. METHODS: Our analytical, prospective, and cross-sectional study included 377 women who participated in a reproductive health campaign during 2015-2016. Anthropometric and gynecological variables were obtained. Cervicovaginal specimens were collected and analyzed with a multiplex in-house PCR to detect Candida spp., Ureaplasma spp., Trichomonas vaginalis, Neisseria gonorrhoeae, HSV, Mycoplasma spp., and Chlamydia trachomatis. RESULTS: The positive cases were 175/377 (46.4%) to at least one of the microorganisms. The most frequent pathogen detected in this population was Ureaplasma spp. (n = 111, 29.4%), followed by Mycoplasma spp. (n = 56, 14.9%) and Candida spp. (n = 47, 12.5%); 33.7% of the positive cases were single infections, whereas 12.7% had coinfection. The multiplex PCR assay was designed targeting nucleotide sequences. CONCLUSIONS: Our data demonstrated that monitoring STIs among asymptomatic patients will encourage target programs to be more precisely and effectively implemented, as well as make these programs more affordable, to benefit society by decreasing the prevalence of STIs.

CTCF chromatin residence time controls three-dimensional genome organization, gene expression and DNA methylation in pluripotent cells.

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Chromosome compartments on the inactive X guide TAD formation independently of transcription during X-reactivation.

A hallmark of chromosome organization is the partition into transcriptionally active A and repressed B compartments, and into topologically associating domains (TADs). Both structures were regarded to be absent from the inactive mouse X chromosome, but to be re-established with transcriptional reactivation and chromatin opening during X-reactivation. Here, we combine a tailor-made mouse iPSC reprogramming system and high-resolution Hi-C to produce a time course combining gene reactivation, chromatin opening and chromosome topology during X-reactivation. Contrary to previous observations, we observe A/B-like compartments on the inactive X harbouring multiple subcompartments. While partial X-reactivation initiates within a compartment rich in X-inactivation escapees, it then occurs rapidly along the chromosome, concomitant with downregulation of Xist. Importantly, we find that TAD formation precedes transcription and initiates from Xist-poor compartments. Here, we show that TAD formation and transcriptional reactivation are causally independent during X-reactivation while establishing Xist as a common denominator.

MyoD induces ARTD1 and nucleoplasmic poly-ADP-ribosylation during fibroblast to myoblast transdifferentiation.

While protein ADP-ribosylation was reported to regulate differentiation and dedifferentiation, it has so far not been studied during transdifferentiation. Here, we found that MyoD-induced transdifferentiation of fibroblasts to myoblasts promotes the expression of the ADP-ribosyltransferase ARTD1. Comprehensive analysis of the genome architecture by Hi-C and RNA-seq analysis during transdifferentiation indicated that ARTD1 locally contributed to A/B compartmentalization and coregulated a subset of MyoD target genes that were however not sufficient to alter transdifferentiation. Surprisingly, the expression of ARTD1 was accompanied by the continuous synthesis of nuclear ADP ribosylation that was neither dependent on the cell cycle nor induced by DNA damage. Conversely to the H2O2-induced ADP-ribosylation, the MyoD-dependent ADP-ribosylation was not associated to chromatin but rather localized to the nucleoplasm. Together, these data describe a MyoD-induced nucleoplasmic ADP-ribosylation that is observed particularly during transdifferentiation and thus potentially expands the plethora of cellular processes associated with ADP-ribosylation.

Intratumoral immunosuppression profiles in 11q-deleted neuroblastomas provide new potential therapeutic targets.

High-risk neuroblastoma (NB) patients with 11q deletion frequently undergo late but consecutive relapse cycles with fatal outcome. To date, no actionable targets to improve current multimodal treatment have been identified. We analyzed immune microenvironment and genetic profiles of high-risk NB correlating with 11q immune status. We show in two independent cohorts that 11q-deleted NB exhibits various immune inhibitory mechanisms, including increased CD4+ resting T cells and M2 macrophages, higher expression of programmed death-ligand 1, interleukin-10, transforming growth factor-beta-1, and indoleamine 2,3-dioxygenase 1 (P < 0.05), and also higher chromosomal breakages (P ≤ 0.02) and hemizygosity of immunosuppressive miRNAs than MYCN-amplified and other 11q-nondeleted high-risk NB. We also analyzed benefits of maintenance treatment in 83 high-risk stage M NB patients focusing on 11q status, either with standard anti-GD2 immunotherapy (n = 50) or previous retinoic acid-based therapy alone (n = 33). Immunotherapy associated with higher EFS (50 vs. 30, P = 0.028) and OS (72 vs. 52, P = 0.047) at 3 years in the overall population. Despite benefits from standard anti-GD2 immunotherapy in high-risk NB patients, those with 11q deletion still face poor outcome. This NB subgroup displays intratumoral immune suppression profiles, revealing a potential therapeutic strategy with combination immunotherapy to circumvent this immune checkpoint blockade.

Prevalence and correlation of human papillomavirus genotypes with clinical factors in cervical samples from Mexican women.

IMPACT STATEMENT: We are submitting data regarding the prevalence and type distribution of the HPV infection and the risk factors associated with it, which may provide a valuable reference to reinforce screening strategies, and to maintain HPV genotype surveillance in Mexico. We discuss the overall prevalence of HPV infection as detected in normal cytological samples stratified by age, different types of infection, and oncogenic capacity. One of the most important findings was that common HPV genotypes detected in healthy women were the genotype numbers: 6, 31, 16, and 56, likewise, smoking and having a history of more than three sexual partners over their lifetime, represented the main risk factors in this study. Furthermore, we found a low frequency of cytological abnormalities and CIN 1-3 in women with HR-HPV.

Oxygen, reactive oxygen species and developmental redox networks: Evo-Devo Evil-Devils?

Molecular oxygen (O2), reactive oxygen species (ROS), and associated redox networks are cornerstones of aerobic life. These molecules and networks have gained recognition as fundamental players in mechanisms that regulate the development of multicellular organisms. First, we present a brief review in which we provide a historical description of some relevant discoveries that led to this recognition. We also discuss the fact that, despite its abundance in nature, oxygen is a limiting factor, and its high availability variation impacted the evolution of adaptive mechanisms to guarantee the proper development of diverse species under such extreme environments. Finally, some examples of when oxygen and ROS were identified as relevant for the control of developmental processes are discussed. We take into account not only the current knowledge on animal redox developmental biology, but also briefly discuss potential scenarios on the origin and evolution of redox developmental mechanisms and the importance of the ever-changing environment.

Adaptation of intestinal epithelial hydrolysis and absorption of dietary carbohydrate and protein in mammals and birds.

The small intestine of mammals and birds exhibits fascinating variation across taxa, body size, and life history features such as locomotion and diet. In the intestine's brush border membrane (BBM), hydrolases are more abundant than transporters in both mammals and birds, but there are differences among the groups in abundance of certain hydrolases and possibly in transporters. For example, mammals express two α-glucosidases, sucrase-isomaltase (SI) and maltase glucoamylase (MGAM), whereas songbirds we studied have only SI, and the chicken expresses SI plus another α-glucosidase that functions similarly to MGAM but is not a true ortholog. For intestinal absorption of sugars and amino acids, small fliers rely on a paracellular pathway to a greater extent than do nonflying mammals, which rely more on transporters. Possibly having evolved in fliers as compensation for lower intestinal nominal surface area (NSA), the fliers' reliance on paracellular absorption is supported by their greater villous surface enlargement that leads to more (per cm2 NSA) tight junctions and greater clearance of passively absorbed compounds. To match digestive capacity to nutrient load, a positive relationship is often observed between dietary intake of macronutrients and intestinal activity of the enzymes and transporters of their respective constituents. In enterocytes, rapid, fine-tuned adjustment to high dietary carbohydrate and protein involves rapid, specific correlated increase in activity and abundance of hydrolases and transporters in the BBM and increases in their mRNA.

Rapid and parallel changes in activity and mRNA of intestinal peptidase to match altered dietary protein levels in juvenile house sparrows (Passer domesticus).

Although dietary flexibility in digestive enzyme activity (i.e. reaction rate) is widespread in vertebrates, mechanisms are poorly understood. When laboratory rats are switched to a higher protein diet, the activities of apical intestinal peptidases increase within 15 h, in some cases by rapid increase in enzyme transcription followed by rapid translation and translocation to the intestine's apical, brush-border membrane (BBM). Focusing on aminopeptidase-N (APN), we studied intestinal digestive enzyme flexibility in birds, relying on activity and mRNA data from the same animals. Our model was nestling house sparrows (Passer domesticus), already known to modulate intestinal peptidase activity when switching between lower and higher protein diets. Twenty-four hours after a switch from an adequate, lower protein diet to a higher protein diet, APN activity was increased in both whole intestinal tissue homogenates and in isolated BBM, but not at 12 h post-diet switch. Twenty-four hours after a reverse switch back to the lower protein diet, APN activity was decreased, but not at 12 h post-diet switch. Changes in APN activity in both diet switch experiments were associated with parallel changes in APN mRNA. Although transcriptional changes seem to be an important mechanism underlying dietary modulation of intestinal peptidase in both nestling house sparrows and laboratory rodents, the time course for modulation in nestlings seemed slower (taking approximately twice as long) compared with laboratory rodents. It may be ecologically advantageous if nestlings biochemically restructure their gut in response to a sustained increase in insects and protein intake rather than one or a few lucky insect meals.

Dietary adaptation to high starch involves increased relative abundance of sucrase-isomaltase and its mRNA in nestling house sparrows.

Dietary flexibility in digestive enzyme activity is widespread in vertebrates but mechanisms are poorly understood. When laboratory rats are switched to a higher carbohydrate diet, the activities of the apical intestinal α-glucosidases (AGs) increase within 6-12 h, mainly by rapid increase in enzyme transcription, followed by rapid translation and translocation to the intestine's apical, brush-border membrane (BBM). We performed the first unified study of the overall process in birds, relying on activity, proteomic, and transcriptomic data from the same animals. Our avian model was nestling house sparrows (Passer domesticus), which switch naturally from a low-starch insect diet to a higher starch seed diet and in whom the protein sucrase-isomaltase (SI) is responsible for all maltase and sucrase intestinal activities. Twenty-four hours after the switch to a high-starch diet, SI activity was increased but not at 12 h post diet switch. SI was the only hydrolase increased in the BBM, and its relative abundance and activity were positively correlated. Twenty-four hours after a reverse switch back to the lower starch diet, SI activity was decreased but not at 12 h post diet switch. Parallel changes in SI mRNA relative abundance were associated with the changes in SI activity in both diet-switch experiments, but our data also revealed an apparent diurnal rhythm in SI mRNA. This is the first demonstration that birds may rely on rapid increase in abundance of SI and its mRNA when adjusting to high-starch diet. Although the mechanisms underlying dietary induction of intestinal enzymes seem similar in nestling house sparrows and laboratory rodents, the time course for modulation in nestlings seemed half as fast compared with laboratory rodents. Before undertaking modulation, an opportunistic forager facing limited resources might rely on more extensive or prolonged environmental sampling, because the redesign of the intestine's hydrolytic capacity shortly after just one or a few meals of a new substrate might be a costly mistake.

Multi-organ transcriptomic landscape of Ambystoma velasci metamorphosis.

Metamorphosis is a postembryonic developmental process that involves morphophysiological and behavioral changes, allowing organisms to adapt into a novel environment. In some amphibians, aquatic organisms undergo metamorphosis to adapt in a terrestrial environment. In this process, these organisms experience major changes in their circulatory, respiratory, digestive, excretory and reproductive systems. We performed a transcriptional global analysis of heart, lung and gills during diverse stages of Ambystoma velasci to investigate its metamorphosis. In our analyses, we identified eight gene clusters for each organ, according to the expression patterns of differentially expressed genes. We found 4064 differentially expressed genes in the heart, 4107 in the lung and 8265 in the gills. Among the differentially expressed genes in the heart, we observed genes involved in the differentiation of cardiomyocytes in the interatrial zone, vasculogenesis and in the maturation of coronary vessels. In the lung, we found genes differentially expressed related to angiogenesis, alveolarization and synthesis of the surfactant protein. In the case of the gills, the most prominent biological processes identified are degradation of extracellular matrix, apoptosis and keratin production. Our study sheds light on the transcriptional responses and the pathways modulation involved in the transformation of the facultative metamorphic salamander A. velasci in an organ-specific manner.