Combining Genetic and Transcriptomic Approaches to Identify Transporter-Coding Genes as Likely Responsible for a Repeatable Salt Tolerance QTL in Citrus

The excessive accumulation of chloride (Cl-) in leaves due to salinity is frequently related to decreased yield in citrus. Two salt tolerance experiments to detect quantitative trait loci (QTLs) for leaf concentrations of Cl-, Na+, and other traits using the same reference progeny derived from the salt-tolerant Cleopatra mandarin (Citrus reshni) and the disease-resistant donor Poncirus trifoliata were performed with the aim to identify repeatable QTLs that regulate leaf Cl- (and/or Na+) exclusion across independent experiments in citrus, as well as potential candidate genes involved. A repeatable QTL controlling leaf Cl- was detected in chromosome 6 (LCl-6), where 23 potential candidate genes coding for transporters were identified using the C. clementina genome as reference. Transcriptomic analysis revealed two important candidate genes coding for a member of the nitrate transporter 1/peptide transporter family (NPF5.9) and a major facilitator superfamily (MFS) protein. Cell wall biosynthesis- and secondary metabolism-related processes appeared to play a significant role in differential gene expression in LCl-6. Six likely gene candidates were mapped in LCl-6, showing conserved synteny in C. reshni. In conclusion, markers to select beneficial Cleopatra mandarin alleles of likely candidate genes in LCl-6 to improve salt tolerance in citrus rootstock breeding programs are provided.

Circulating mitochondrial genes detect acute cardiac allograft rejection: Role of the mitochondrial calcium uniporter complex

Acute rejection after heart transplantation increases the risk of chronic dysfunction. Disturbances in mitochondrial function may play a contributory role, however, the relationship between histological signs of rejection in the human transplanted heart and expression levels of circulating mitochondrial genes, such as the mitochondrial Ca²⁺ uniporter (MCU) complex, remains unexplored. We conducted an RNA-sequencing analysis to identify altered mitochondrial genes in serum and to evaluate their diagnostic accuracy for rejection episodes. We included 40 consecutive samples from transplant recipients undergoing routine endomyocardial biopsies. In total, 112 mitochondrial genes were identified in the serum of posttransplant patients, of which 28 were differentially expressed in patients with acute rejection (p < .05). Considering the receiver operating characteristic analysis with an area under the curve (AUC) >0.900 to discriminate patients with moderate or severe degrees of rejection, we found that the MCU system showed a strong capability for detection: MCU (AUC = 0.944, p < .0001), MCU/MCUR1 ratio (AUC = 0.972, p < .0001), MCU/MCUB ratio (AUC = 0.970, p < .0001), and MCU/MICU1 ratio (AUC = 0.970, p < .0001). Mitochondrial alterations are reflected in peripheral blood and are capable of discriminating between patients with allograft rejection and those not experiencing rejection with excellent accuracy. The dysregulation of the MCU complex was found to be the most relevant finding.

Transcriptional Response in the Digestive Gland of the King Scallop (Pecten maximus) After the Injection of Domoic Acid

Some diatom species of the genus Pseudo-nitzschia produce the toxin domoic acid. The depuration rate of domoic acid in Pecten maximus is very low; for this reason, king scallops generally contain high levels of domoic acid in their tissues. A transcriptomic approach was used to identify the genes differentially expressed in the P. maximus digestive gland after the injection of domoic acid. The differential expression analysis found 535 differentially expressed genes (226 up-regulated and 309 down-regulated). Protein–protein interaction networks obtained with the up-regulated genes were enriched in gene ontology terms, such as vesicle-mediated transport, response to stress, signal transduction, immune system process, RNA metabolic process, and autophagy, while networks obtained with the down-regulated genes were enriched in gene ontology terms, such as response to stress, immune system process, ribosome biogenesis, signal transduction, and mRNA processing. Genes that code for cytochrome P450 enzymes, glutathione S-transferase theta-1, glutamine synthase, pyrroline-5-carboxylate reductase 2, and sodium- and chloride-dependent glycine transporter 1 were among the up-regulated genes. Therefore, a stress response at the level of gene expression, that could be caused by the domoic acid injection, was evidenced by the alteration of several biological, cellular, and molecular processes.

Allelic Variation of MYB10 Is the Major Force Controlling Natural Variation in Skin and Flesh Color in Strawberry (Fragaria spp.) Fruit

The fruits of diploid and octoploid strawberry (Fragaria spp) show substantial natural variation in color due to distinct anthocyanin accumulation and distribution patterns. Anthocyanin biosynthesis is controlled by a clade of R2R3 MYB transcription factors, among which MYB10 is the main activator in strawberry fruit. Here, we show that mutations in MYB10 cause most of the variation in anthocyanin accumulation and distribution observed in diploid woodland strawberry (F. vesca) and octoploid cultivated strawberry (F ×ananassa). Using a mapping-by-sequencing approach, we identified a gypsy-transposon in MYB10 that truncates the protein and knocks out anthocyanin biosynthesis in a white-fruited F. vesca ecotype. Two additional loss-of-function mutations in MYB10 were identified among geographically diverse white-fruited F. vesca ecotypes. Genetic and transcriptomic analyses of octoploid Fragaria spp revealed that FaMYB10-2, one of three MYB10 homoeologs identified, regulates anthocyanin biosynthesis in developing fruit. Furthermore, independent mutations in MYB10-2 are the underlying cause of natural variation in fruit skin and flesh color in octoploid strawberry. We identified a CACTA-like transposon (FaEnSpm-2) insertion in the MYB10-2 promoter of red-fleshed accessions that was associated with enhanced expression. Our findings suggest that cis-regulatory elements in FaEnSpm-2 are responsible for enhanced MYB10-2 expression and anthocyanin biosynthesis in strawberry fruit flesh.

Combination of phenotype and polygenic risk score in breast cancer risk evaluation in the Spanish population: a case -control study

Background

In recent years, the identification of genetic and phenotypic biomarkers of cancer for prevention, early diagnosis and patient stratification has been a main objective of research in the field. Different multivariable models that use biomarkers have been proposed for the evaluation of individual risk of developing breast cancer.

Methods

This is a case control study based on a population-based cohort. We describe and evaluate a multivariable model that incorporates 92 Single-nucleotide polymorphisms (SNPs) (Supplementary Table S1) and five different phenotypic variables and which was employed in a Spanish population of 642 healthy women and 455 breast cancer patients.

Results

Our model allowed us to stratify two groups: high and low risk of developing breast cancer. The 9th decile included 1% of controls vs 9% of cases, with an odds ratio (OR) of 12.9 and a p-value of 3.43E-07. The first decile presented an inverse proportion: 1% of cases and 9% of controls, with an OR of 0.097 and a p-value of 1.86E-08.

Conclusions

These results indicate the capacity of our multivariable model to stratify women according to their risk of developing breast cancer. The major limitation of our analysis is the small cohort size. However, despite the limitations, the results of our analysis provide proof of concept in a poorly studied population, and opens up the possibility of using this method in the routine screening of the Spanish population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-020-07584-9.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-020-07584-9.

Prefoldin-like Bud27 influences the transcription of ribosomal components and ribosome biogenesis in Saccharomyces cerevisiae

Understanding the functional connection that occurs for the three nuclear RNA polymerases to synthesize ribosome components during the ribosome biogenesis process has been the focal point of extensive research. To preserve correct homeostasis on the production of ribosomal components, cells might require the existence of proteins that target a common subunit of these RNA polymerases to impact their respective activities. This work describes how the yeast prefoldin-like Bud27 protein, which physically interacts with the Rpb5 common subunit of the three RNA polymerases, is able to modulate the transcription mediated by the RNA polymerase I, likely by influencing transcription elongation, the transcription of the RNA polymerase III, and the processing of ribosomal RNA. Bud27 also regulates both RNA polymerase II-dependent transcription of ribosomal proteins and ribosome biogenesis regulon genes, likely by occupying their DNA ORFs, and the processing of the corresponding mRNAs. With RNA polymerase II, this association occurs in a transcription rate-dependent manner. Our data also indicate that Bud27 inactivation alters the phosphorylation kinetics of ribosomal protein S6, a readout of TORC1 activity. We conclude that Bud27 impacts the homeostasis of the ribosome biogenesis process by regulating the activity of the three RNA polymerases and, in this way, the synthesis of ribosomal components. This quite likely occurs through a functional connection of Bud27 with the TOR signaling pathway.

Expression of miR-135b in Psoriatic Skin and Its Association with Disease Improvement

miRNAs have been associated with psoriasis since just over a decade. However, we are far from a complete understanding of their role during the development of this disease. Our objective was to characterize the cutaneous expression of miRNAs not previously described in psoriasis, the changes induced following the treatment with biologicals and their association with disease improvement. Next generation sequencing was performed from five skin samples from psoriasis patients (lesional and non-lesional skin) and five controls, and from this cohort, 12 microRNAs were selected to be analyzed in skin samples from 44 patients with plaque psoriasis. In 15 patients, an additional sample was obtained after three months of biological treatment. MiR-9-5p, miR-133a-3p and miR-375 were downregulated in the lesional skin of psoriasis patients. After treatment, expression of miR-133a-3p, miR-375, miR-378a and miR-135b in residual lesions returned towards the levels observed in non-lesional skin. The decrease in miR-135b levels after treatment with biologics was associated with both the improvement of patients evaluated through Psoriasis Area and Severity Index score and the decrease in local inflammatory response. Moreover, basal expression of miR-135b along with age was associated with the improvement of psoriasis, suggesting its possible usefulness as a prognostic biomarker.

Indigenous Pseudomonas spp. Strains from the Olive (Olea europaea L.) Rhizosphere as Effective Biocontrol Agents against Verticillium dahliae: From the Host Roots to the Bacterial Genomes

The use of biological control agents (BCA), alone or in combination with other management measures, has gained attention over the past decades, driven by the need to seek for sustainable and eco-friendly alternatives to confront plant pathogens. The rhizosphere of olive (Olea europaea L.) plants is a source of bacteria with potential as biocontrol tools against Verticillium wilt of olive (VWO) caused by Verticillium dahliae Kleb. A collection of bacterial isolates from healthy nursery-produced olive (cultivar Picual, susceptible to VWO) plants was generated based on morphological, biochemical and metabolic characteristics, chemical sensitivities, and on their in vitro antagonistic activity against several olive pathogens. Three strains (PIC25, PIC105, and PICF141) showing high in vitro inhibition ability of pathogens' growth, particularly against V. dahliae, were eventually selected. Their effectiveness against VWO caused by the defoliating pathotype of V. dahliae was also demonstrated, strain PICF141 being the rhizobacteria showing the best performance as BCA. Genotypic and phenotypic traits traditionally associated with plant growth promotion and/or biocontrol abilities were evaluated as well (e.g., phytase, xylanase, catalase, cellulase, chitinase, glucanase activities, and siderophore and HCN production). Multi-locus sequence analyses of conserved genes enabled the identification of these strains as Pseudomonas spp. Strain PICF141 was affiliated to the “Pseudomonas mandelii subgroup,” within the “Pseudomonas fluorescens group,” Pseudomonas lini being the closest species. Strains PIC25 and PIC105 were affiliated to the “Pseudomonas aeruginosa group,” Pseudomonas indica being the closest relative. Moreover, we identified P. indica (PIC105) for the first time as a BCA. Genome sequencing and in silico analyses allowed the identification of traits commonly associated with plant-bacteria interactions. Finally, the root colonization ability of these olive rhizobacteria was assessed, providing valuable information for the future development of formulations based on these strains. A set of actions, from rhizosphere isolation to genome analysis, is proposed and discussed for selecting indigenous rhizobacteria as effective BCAs.

Transcriptional response after exposure to domoic acid-producing Pseudo-nitzschia in the digestive gland of the mussel Mytilus galloprovincialis

Bivalve molluscs are filter feeding species that can accumulate biotoxins in their body tissues during harmful algal blooms. Amnesic Shellfish Poisoning (ASP) is caused by species of the diatom genus Pseudo-nitzschia, which produces the toxin domoic acid. The Mytilus galloprovincialis digestive gland transcriptome was de novo assembled based on the sequencing of 12 cDNA libraries, six obtained from control mussels and six from mussels naturally exposed to domoic acid-producing diatom Pseudo-nitzschia australis. After de novo assembly 94,727 transcripts were obtained, with an average length of 1015 bp and a N50 length of 761 bp. The assembled transcripts were clustered (homology > 90%) into 69,294 unigenes. Differential gene expression analysis was performed (DESeq2 algorithm) in the digestive gland following exposure to the toxic algae. A total of 1158 differentially expressed unigenes (absolute fold change > 1.5 and p-value < 0.05) were detected: 686 up-regulated and 472 down-regulated. Several membrane transporters belonging to the family of the SLC (solute carriers) were over-expressed in exposed mussels. Functional enrichment was performed using Pfam annotations obtained from the genes differentially expressed, 37 Pfam families were found to be significantly (FDR adjusted p-value < 0.1) enriched. Some of these families (sulfotransferases, aldo/keto reductases, carboxylesterases, C1q domain and fibrinogen C-terminal globular domain) could be putatively involved in detoxification processes, in the response against of the oxidative stress and in immunological processes. Protein network analysis with STRING algorithm found alteration of the Notch signaling pathway under the action of domoic acid-producing Pseudo-nitzschia. In conclusion, this study provides a high quality reference transcriptome of M. galloprovincialis digestive gland and identifies potential genes involved in the response to domoic acid.

Exosomes from bulk and stem cells from human prostate cancer have a differential microRNA content that contributes cooperatively over local and pre-metastatic niche

The different prostate cancer (PCa) cell populations (bulk and cancer stem cells, CSCs) release exosomes that contain miRNAs that could modify the local or premetastatic niche. The analysis of the differential expression of miRNAs in exosomes allows evaluating the differential biological effect of both populations on the niche, and the identification of potential biomarkers and therapeutic targets. Five PCa primary cell cultures were established to originate bulk and CSCs cultures. From them, exosomes were purified by precipitation for miRNAs extraction to perform a comparative profile of miRNAs by next generation sequencing in an Illumina platform. 1839 miRNAs were identified in the exosomes. Of these 990 were known miRNAs, from which only 19 were significantly differentially expressed: 6 were overexpressed in CSCs and 13 in bulk cells exosomes. miR-100-5p and miR-21-5p were the most abundant miRNAs. Bioinformatics analysis indicated that differentially expressed miRNAs are highly related with PCa carcinogenesis, fibroblast proliferation, differentiation and migration, and angiogenesis. Besides, miRNAs from bulk cells affects osteoblast differentiation. Later, their effect was evaluated in normal prostate fibroblasts (WPMY-1) where transfection with miR-100-5p, miR-21-5p and miR-139-5p increased the expression of metalloproteinases (MMPs) -2, -9 and -13 and RANKL and fibroblast migration. The higher effect was achieved with miR21 transfection. As conclusion, miRNAs have a differential pattern between PCa bulk and CSCs exosomes that act collaboratively in PCa progression and metastasis. The most abundant miRNAs in PCa exosomes are interesting potential biomarkers and therapeutic targets.

Quantitative global studies of reactomes and metabolomes using a vectorial representation of reactions and chemical compounds

Background

Global studies of the protein repertories of organisms are providing important information on the characteristics of the protein space. Many of these studies entail classification of the protein repertory on the basis of structure and/or sequence similarities. The situation is different for metabolism. Because there is no good way of measuring similarities between chemical reactions, there is a barrier to the development of global classifications of "metabolic space" and subsequent studies comparable to those done for protein sequences and structures.

Results

In this work, we propose a vectorial representation of chemical reactions, which allows them to be compared and classified. In this representation, chemical compounds, reactions and pathways may be represented in the same vectorial space. We show that the representation of chemical compounds reflects their physicochemical properties and can be used for predictive purposes. We use the vectorial representations of reactions to perform a global classification of the reactome of the model organism E. coli.

Conclusions

We show that this unsupervised clustering results in groups of enzymes more coherent in biological terms than equivalent groupings obtained from the EC hierarchy. This hierarchical clustering produces an optimal set of 21 groups which we analyzed for their biological meaning.

Species connectivities and reaction mechanisms from neutral response experiments

We develop a new method for obtaining connectivity data for nonlinear reaction networks, based on linear response experiments. In our approach the linear response is not the result of an approximation procedure but is due to the appropriate design of the response experiments, that is (1) they are carried out with the preservation of constant values for the total (labeled plus unlabeled) input and output fluxes and (2) the labeled compounds obey a neutrality condition (i.e., they have practically the same kinetic and transport properties as the unlabeled compounds). Under these circumstances the linear response equations hold even though the kinetics of the process is highly nonlinear. On the basis of this linear response law, we develop a method for evaluating reaction connectivities in biochemical networks from stationary response experiments. Given a system in a stationary regime, a pulse of a labeled species is introduced (with conservation of the total flux) and then the response of all the species of the network is recorded. The mechanistic information is contained in a connectivity matrix, K, which can be evaluated from the response data by means of differential as well as integral methods. The approach does not require any prior knowledge of the reaction mechanism. We carried out a numerical study of the method, based on a two-step procedure. Starting from a known reaction mechanism, we generated response data sets, to which we add noise; then, we use the noisy data sets for retrieving the connectivity matrix. The calculations were done with two programs written in Mathematica: the urea cycle and the upper part of glycolysis are used as sample biochemical networks. Given enough computer power, there are no limitations concerning the number of species involved in the response experiments; on current desktop systems processing responses of teens of species would take a few hours. The method is limited by the occurrence of experimental errors: if experimental errors in the evaluation of fluxes are larger than 10%, the method may fail to reproduce the correct values of some elements of the connectivity matrix.