First Transcriptome Analysis of Hepatoblastoma in Brazil: Unraveling the Pivotal Role of Noncoding RNAs and Metabolic Pathways.

Hepatoblastoma stands as the most prevalent liver cancer in the pediatric population. Characterized by a low mutational burden, chromosomal and epigenetic alterations are key drivers of its tumorigenesis. Transcriptome analysis is a powerful tool for unraveling the molecular intricacies of hepatoblastoma, shedding light on the effects of genetic and epigenetic changes on gene expression. In this study conducted in Brazilian patients, an in-depth whole transcriptome analysis was performed on 14 primary hepatoblastomas, compared to control liver tissues. The analysis unveiled 1,492 differentially expressed genes (1,031 upregulated and 461 downregulated), including 920 protein-coding genes (62%). Upregulated biological processes were linked to cell differentiation, signaling, morphogenesis, and development, involving known hepatoblastoma-associated genes (DLK1, MEG3, HDAC2, TET1, HMGA2, DKK1, DKK4), alongside with novel findings (GYNG4, CDH3, and TNFRSF19). Downregulated processes predominantly centered around oxidation and metabolism, affecting amines, nicotinamides, and lipids, featuring novel discoveries like the repression of SYT7, TTC36, THRSP, CCND1, GCK and CAMK2B. Two genes, which displayed a concordant pattern of DNA methylation alteration in their promoter regions and dysregulation in the transcriptome, were further validated by RT-qPCR: the upregulated TNFRSF19, a key gene in the embryonic development, and the repressed THRSP, connected to lipid metabolism. Furthermore, based on protein-protein interaction analysis, we identified genes holding central positions in the network, such as HDAC2, CCND1, GCK, and CAMK2B, among others, that emerged as prime candidates warranting functional validation in future studies. Notably, a significant dysregulation of non-coding RNAs (ncRNAs), predominantly upregulated transcripts, was observed, with 42% of the top 50 highly expressed genes being ncRNAs. An integrative miRNA-mRNA analysis revealed crucial biological processes associated with metabolism, oxidation reactions of lipids and carbohydrates, and methylation-dependent chromatin silencing. In particular, four upregulated miRNAs (miR-186, miR-214, miR-377, and miR-494) played a pivotal role in the network, potentially targeting multiple protein-coding transcripts, including CCND1 and CAMK2B. In summary, our transcriptome analysis highlighted disrupted embryonic development as well as metabolic pathways, particularly those involving lipids, emphasizing the emerging role of ncRNAs as epigenetic regulators in hepatoblastomas. These findings provide insights into the complexity of the hepatoblastoma transcriptome and identify potential targets for future therapeutic interventions.

Genome-wide analysis and characterization of the LRR-RLK gene family provides insights into anthracnose resistance in common bean.

Anthracnose, caused by the hemibiotrophic fungus Colletotrichum lindemuthianum, is a damaging disease of common beans that can drastically reduce crop yield. The most effective strategy to manage anthracnose is the use of resistant cultivars. There are many resistance loci that have been identified, mapped and associated with markers in common bean chromosomes. The Leucine-rich repeat kinase receptor protein (LRR-RLK) family is a diverse group of transmembrane receptors, which potentially recognizes pathogen-associated molecular patterns and activates an immune response. In this study, we performed in silico analyses to identify, classify, and characterize common bean LRR-RLKs, also evaluating their expression profile in response to the infection by C. lindemuthianum. By analyzing the entire genome of Phaseolus vulgaris, we could identify and classify 230 LRR-RLKs into 15 different subfamilies. The analyses of gene structures, conserved domains and motifs suggest that LRR-RLKs from the same subfamily are consistent in their exon/intron organization and composition. LRR-RLK genes were found along the 11 chromosomes of the species, including regions of proximity with anthracnose resistance markers. By investigating the duplication events within the LRR-RLK family, we associated the importance of such a family with an expansion resulting from a strong stabilizing selection. Promoter analysis was also performed, highlighting cis-elements associated with the plant response to biotic stress. With regard to the expression pattern of LRR-RLKs in response to the infection by C. lindemuthianum, we could point out several differentially expressed genes in this subfamily, which were associated to specific molecular patterns of LRR-RLKs. Our work provides a broad analysis of the LRR-RLK family in P. vulgaris, allowing an in-depth structural and functional characterization of genes and proteins of this family. From specific expression patterns related to anthracnose response, we could infer a direct participation of RLK-LRR genes in the mechanisms of resistance to anthracnose, highlighting important subfamilies for further investigations.

Identification of a key protein set involved in Moniliophthora perniciosa necrotrophic mycelium and basidiocarp development.

Moniliophthora perniciosa is a hemibiotrophic fungus that causes witches' broom disease in cacao (Theobroma cacao L.). The biotrophic fungal phase initiates the disease and is characterized by a monokaryotic mycelium, while the necrotrophic phase is characterized by a dikaryotic mycelium and leads to necrosis of infected tissues. A study of the necrotrophic phase was conducted on bran-based solid medium, which is the only medium that enables basidiocarp and basidiospore production. Six different fungal developmental phases were observed according to the mycelium colour or the organ produced: white, yellow, pink, dark pink, primordium and basidiocarp. In this study, we identified notable proteins in each phase, particularly those accumulated prior to basidiocarp formation. Proteins were analysed by proteomics; 2-D gels showed 300-550 spots. Statistically differentially accumulated spots were sequenced by mass spectrometry and 259 proteins were identified and categorized into nine functional classes. Proteins related to energy metabolism, protein folding and morphogenesis that were potentially involved in primordium and basidiocarp formation were identified; these proteins may represent useful candidates for further analysis related to the spread and pathogenesis of this fungus. To the best of our knowledge, this report describes the first proteomic analysis of the developmental phases of Moniliophthora perniciosa.

Moniliophthora perniciosa development: key genes involved in stress-mediated cell wall organization and autophagy.

Moniliophthora perniciosa is a basidiomycete responsible for the witches' broom disease in cacao (Theobroma cacao L.). Chitin synthase (CHS), chitinase (CHIT) and autophagy (ATG) genes have been associated to stress response preceding the formation of basidiocarp. An analysis of literature mining, interactomics and gene expression was developed to identify the main proteins related to development, cell wall organization and autophagy in M. perniciosa. TORC2 complex elements were identified and were involved in the response to the nutrient starvation during the fungus development stages preceding the basidiocarp formation. This complex interacted with target proteins related to cell wall synthesis and to polarization and cell division (FKS1, CHS, CDC42, ROM2). Autolysis and autophagy processes were associated to CHIT2, ATG8 and to the TORC1 complex (TOR1 and KOG1), which is central in the upstream signalization of the stress response due to nutrient starvation and growth regulation. Other important elements that participate to steps preceding basidiocarp formation were also identified (KOG1, SSZ1, GDI1, FKS1, CCD10, CKS1, CDC42, RHO1, AVO1, BAG7). Similar gene expression patterns during fungus reproductive structure formation and when treated by rapamycin (a nutritional related-autophagy stress agent) were observed: cell division related-genes were repressed while those related to autolysis/autophagy were overexpressed.

Phylogenetic Origin of Primary and Secondary Metabolic Pathway Genes Revealed by C. maxima and C. reticulata Diagnostic SNPs.

Modern cultivated Citrus species and varieties result from interspecific hybridization between four ancestral taxa. Among them, Citrus maxima and Citrus reticulata, closely associated with the pummelo and mandarin horticultural groups, respectively, were particularly important as the progenitors of sour and sweet oranges (Citrus aurantium and Citrus sinensis), grapefruits (Citrus paradisi), and hybrid types resulting from modern breeding programs (tangors, tangelos, and orangelos). The differentiation between the four ancestral taxa and the phylogenomic structure of modern varieties widely drive the phenotypic diversity's organization. In particular, strong phenotypic differences exist in the coloration and sweetness and represent important criteria for breeders. In this context, focusing on the genes of the sugar, carotenoid, and chlorophyll biosynthesis pathways, the aim of this work was to develop a set of diagnostic single-nucleotide polymorphism (SNP) markers to distinguish the ancestral haplotypes of C. maxima and C. reticulata and to provide information at the intraspecific diversity level (within C. reticulata or C. maxima). In silico analysis allowed the identification of 3,347 SNPs from selected genes. Among them, 1,024 were detected as potential differentiation markers between C. reticulata and C. maxima. A total of 115 SNPs were successfully developed using a competitive PCR technology. Their transferability among all Citrus species and the true citrus genera was very good, with only 0.87% of missing data. The ancestral alleles of the SNPs were identified, and we validated the usefulness of the developed markers for tracing the ancestral haplotype in large germplasm collections and sexually recombined progeny issued from the C. reticulata/C. maxima admixture gene pool. These markers will pave the way for targeted association studies based on ancestral haplotypes.

Comparative proteomics of two citrus varieties in response to infection by the fungus Alternaria alternata.

Alternaria brown spot (ABS) is a disease caused by the necrotrophic fungus Alternaria alternata, which induces necrotic lesions on fruits and young leaves due to the production of the host-specific ACT toxin by the fungus. To better understand the citrus-A. alternata interaction and to identify putative resistance proteins, as well as the receptor of the ACT toxin, citrus plants susceptible ('Minneola' mandarin) and resistant ('Clemenules' tangor) to A. alternata, infected or not (control) with the pathogen were analyzed by proteomics. Protein changes were observed between citrus genotypes after infection, and 150 candidate proteins were obtained. A general scheme of the metabolic processes involved in susceptible and resistant citrus-A. alternata interactions was designed. Susceptible plants presented a high level of proteins involved in stress response at the final stages of the infection, whereas resistant plants presented high level of ROS proteins, metabolic proteins, and proteins involved in the immune system process. Proteins like ferredoxin and cyclophilin are specific to the susceptible variety and may be good candidates as fungal effector-interacting proteins. This is the first citrus-A. alternata proteomics analysis, which has allowed a better understanding of the molecular bases of the citrus response to ABS disease.

Different adaptation strategies of two citrus scion/rootstock combinations in response to drought stress.

Scion/rootstock interaction is important for plant development and for breeding programs. In this context, polyploid rootstocks presented several advantages, mainly in relation to biotic and abiotic stresses. Here we analyzed the response to drought of two different scion/rootstock combinations presenting different polyploidy: the diploid (2x) and autotetraploid (4x) Rangpur lime (Citrus limonia, Osbeck) rootstocks grafted with 2x Valencia Delta sweet orange (Citrus sinensis) scions, named V/2xRL and V/4xRL, respectively. Based on previous gene expression data, we developed an interactomic approach to identify proteins involved in V/2xRL and V/4xRL response to drought. A main interactomic network containing 3,830 nodes and 97,652 edges was built from V/2xRL and V/4xRL data. Exclusive proteins of the V/2xRL and V/4xRL networks (2,056 and 1,001, respectively), as well as common to both networks (773) were identified. Functional clusters were obtained and two models of drought stress response for the V/2xRL and V/4xRL genotypes were designed. Even if the V/2xRL plant implement some tolerance mechanisms, the global plant response to drought was rapid and quickly exhaustive resulting in a general tendency to dehydration avoidance, which presented some advantage in short and strong drought stress conditions, but which, in long terms, does not allow the plant survival. At the contrary, the V/4xRL plants presented a response which strong impacts on development but that present some advantages in case of prolonged drought. Finally, some specific proteins, which presented high centrality on interactomic analysis were identified as good candidates for subsequent functional analysis of citrus genes related to drought response, as well as be good markers of one or another physiological mechanism implemented by the plants.

The pathogenesis-related protein PR-4b from Theobroma cacao presents RNase activity, Ca(2+) and Mg(2+) dependent-DNase activity and antifungal action on Moniliophthora perniciosa.

BACKGROUND: The production and accumulation of pathogenesis-related proteins (PR proteins) in plants in response to biotic or abiotic stresses is well known and is considered as a crucial mechanism for plant defense. A pathogenesis-related protein 4 cDNA was identified from a cacao-Moniliophthora perniciosa interaction cDNA library and named TcPR-4b. RESULTS: TcPR-4b presents a Barwin domain with six conserved cysteine residues, but lacks the chitin-binding site. Molecular modeling of TcPR-4b confirmed the importance of the cysteine residues to maintain the protein structure, and of several conserved amino acids for the catalytic activity. In the cacao genome, TcPR-4b belonged to a small multigene family organized mainly on chromosome 5. TcPR-4b RT-qPCR analysis in resistant and susceptible cacao plants infected by M. perniciosa showed an increase of expression at 48 hours after infection (hai) in both cacao genotypes. After the initial stage (24-72 hai), the TcPR-4b expression was observed at all times in the resistant genotypes, while in the susceptible one the expression was concentrated at the final stages of infection (45-90 days after infection). The recombinant TcPR-4b protein showed RNase, and bivalent ions dependent-DNase activity, but no chitinase activity. Moreover, TcPR-4b presented antifungal action against M. perniciosa, and the reduction of M. perniciosa survival was related to ROS production in fungal hyphae. CONCLUSION: To our knowledge, this is the first report of a PR-4 showing simultaneously RNase, DNase and antifungal properties, but no chitinase activity. Moreover, we showed that the antifungal activity of TcPR-4b is directly related to RNase function. In cacao, TcPR-4b nuclease activities may be related to the establishment and maintenance of resistance, and to the PCD mechanism, in resistant and susceptible cacao genotypes, respectively.

Polygalacturonases from Moniliophthora perniciosa are regulated by fermentable carbon sources and possible post-translational modifications.

We report the first molecular and in silico analysis of Monilophthora perniciosa polygalacturonases (PGs). Three MpPG genes (MpPG1, MpPG2 and MpPG3) were identified and analyzed at transcriptional level, by RT-qPCR, in dikaryotic M. perniciosa mycelium grown on solid-bran based medium and on liquid medium supplemented with different fermentable and non-fermentable carbon sources. The MpPG genes presented different expression patterns suggesting different individual regulation. However, all are mainly regulated by fermentable carbon sources (galactose and mannose). The integrated analysis of PG gene expression and systems biology (using MpG1 and MpG2 orthologs in Neurospora crassa, named NCU06961 and NCU02369, respectively) allowed identifying some possible mechanism of protein regulation during the necrotrophic fungal phase. MpPG1-NCU06961 and MpPG2-NCU02369 directly or indirectly interacted with central and highly connected proteins involved in protein synthesis and protein regulation associated to post-translational modifications, in cell wall metabolism, and in cellular metabolism related to energy production. This analysis also allowed the identification of key proteins for further studies of M. perniciosa development and/or for disease management, such as MpPG2, a pectin methylesterase, an acetolactate synthase and the small ubiquitin-like modifier SMT3-like.