Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves.

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the "nuclear compensation hypothesis," a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with "doubly uniparental" mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Elucidating the Genetic Relationships on the Original Old Sicilian Triticum Spp. Collection by SNP Genotyping.

Several Triticum species spread in cultivation in Sicily and neighboring regions over the centuries, which led to the establishment of a large genetic diversity. Many ancient varieties were widely cultivated until the beginning of the last century before being replaced by modern varieties. Recently, they have been reintroduced in cultivation in Sicily. Here, the genetic diversity of 115 and 11 accessions from Sicily and Calabria, respectively, belonging to Triticum species was evaluated using a high-density SNP array. Einkorn, emmer, and spelta wheat genotypes were used as outgroups for species and subspecies; five modern varieties of durum and bread wheat were used as references. A principal coordinates analysis (PCoA) and an unweighted pair group method with arithmetic mean (UPGMA) showed four distinct groups among Triticum species and T. turgidum subspecies. The population structure analysis distinguished five gene pools, among which three appeared private to the T. aestivum, T. turgidum subsp. Turgidum, and 'Timilia' group. The principal component analysis (PCA) displayed a bio-morphological trait relationship of a subset (110) of ancient wheat varieties and their wide variability within the T. turgidum subsp. durum subgroups. A discriminant analysis of principal components (DAPC) and phylogenetic analyses applied to the four durum wheat subgroups revealed that the improved varieties harbored a different gene pool compared to the most ancient varieties. The 'Russello' and 'Russello Ibleo' groups were distinguished; both displayed higher genetic variability compared to the 'Timilia' group accessions. This research represents a comprehensive approach to fingerprinting the old wheat Sicilian germplasm, which is useful in avoiding commercial fraud and sustaining the cultivation of landraces and ancient varieties.

Integrated omics approach reveals the molecular pathways activated in tomato by Kocuria rhizophila, a soil plant growth-promoting bacterium.

Plant microbial biostimulants application has become a promising and eco-friendly agricultural strategy to improve crop yields, reducing chemical inputs for more sustainable cropping systems. The soil dwelling bacterium Kocuria rhizophila was previously characterized as Plant Growth Promoting Bacteria (PGPB) for its multiple PGP traits, such as indole-3-acetic acid production, phosphate solubilization capability and salt and drought stress tolerance. Here, we evaluated by a multi-omics approach, the PGP activity of K. rhizophila on tomato, revealing the molecular pathways by which it promotes plant growth. Transcriptomic analysis showed several up-regulated genes mainly related to amino acid metabolism, cell wall organization, lipid and secondary metabolism, together with a modulation in the DNA methylation profile, after PGPB inoculation. In agreement, proteins involved in photosynthesis, cell division, and plant growth were highly accumulated by K. rhizophila. Furthermore, "amino acid and peptides", "monosaccharides", and "TCA" classes of metabolites resulted the most affected by PGPB treatment, as well as dopamine, a catecholamine neurotransmitter mediating plant growth through S-adenosylmethionine decarboxylase (SAMDC), a gene enhancing the vegetative growth, up-regulated in tomato by K. rhizophila treatment. Interestingly, eight gene modules well correlated with differentially accumulated proteins (DAPs) and metabolites (DAMs), among which two modules showed the highest correlation with nine proteins, including a nucleoside diphosphate kinase, and cytosolic ascorbate peroxidase, as well as with several amino acids and metabolites involved in TCA cycle. Overall, our findings highlighted that sugars and amino acids, energy regulators, involved in tomato plant growth, were strongly modulated by the K. rhizophila-plant interaction.

Short-term transcriptomic analysis at organ scale reveals candidate genes involved in low N responses in NUE-contrasting tomato genotypes.

Background: Understanding the complex regulatory network underlying plant nitrogen (N) responses associated with high Nitrogen Use Efficiency (NUE) is one of the main challenges for sustainable cropping systems. Nitrate (NO3 -), acting as both an N source and a signal molecule, provokes very fast transcriptome reprogramming, allowing plants to adapt to its availability. These changes are genotype- and tissue-specific; thus, the comparison between contrasting genotypes is crucial to uncovering high NUE mechanisms. Methods: Here, we compared, for the first time, the spatio-temporal transcriptome changes in both root and shoot of two NUE contrasting tomato genotypes, Regina Ostuni (high-NUE) and UC82 (low-NUE), in response to short-term (within 24 h) low (LN) and high (HN) NO3 - resupply. Results: Using time-series transcriptome data (0, 8, and 24 h), we identified 395 and 482 N-responsive genes differentially expressed (DEGs) between RO and UC82 in shoot and root, respectively. Protein kinase signaling plant hormone signal transduction, and phenylpropanoid biosynthesis were the main enriched metabolic pathways in shoot and root, respectively, and were upregulated in RO compared to UC82. Interestingly, several N transporters belonging to NRT and NPF families, such as NRT2.3, NRT2.4, NPF1.2, and NPF8.3, were found differentially expressed between RO and UC82 genotypes, which might explain the contrasting NUE performances. Transcription factors (TFs) belonging to several families, such as ERF, LOB, GLK, NFYB, ARF, Zinc-finger, and MYB, were differentially expressed between genotypes in response to LN. A complementary Weighted Gene Co-expression Network Analysis (WGCNA) allowed the identification of LN-responsive co-expression modules in RO shoot and root. The regulatory network analysis revealed candidate genes that might have key functions in short-term LN regulation. In particular, an asparagine synthetase (ASNS), a CBL-interacting serine/threonine-protein kinase 1 (CIPK1), a cytokinin riboside 5'-monophosphate phosphoribohydrolase (LOG8), a glycosyltransferase (UGT73C4), and an ERF2 were identified in the shoot, while an LRR receptor-like serine/threonine-protein kinase (FEI1) and two TFs NF-YB5 and LOB37 were identified in the root. Discussion: Our results revealed potential candidate genes that independently and/or concurrently may regulate short-term low-N response, suggesting a key role played by cytokinin and ROS balancing in early LN regulation mechanisms adopted by the N-use efficient genotype RO.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families.

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

Transcriptome changes induced by Arbuscular mycorrhizal symbiosis in leaves of durum wheat (Triticum durum Desf.) promote higher salt tolerance.

The salinity of soil is a relevant environmental problem around the world, with climate change raising its relevance, particularly in arid and semiarid areas. Arbuscular Mycorrhizal Fungi (AMF) positively affect plant growth and health by mitigating biotic and abiotic stresses, including salt stress. The mechanisms through which these benefits manifest are, however, still unclear. This work aimed to identify key genes involved in the response to salt stress induced by AMF using RNA-Seq analysis on durum wheat (Triticum turgidum L. subsp. durum Desf. Husn.). Five hundred sixty-three differentially expressed genes (DEGs), many of which involved in pathways related to plant stress responses, were identified. The expression of genes involved in trehalose metabolism, RNA processing, vesicle trafficking, cell wall organization, and signal transduction was significantly enhanced by the AMF symbiosis. A downregulation of genes involved in both enzymatic and non-enzymatic oxidative stress responses as well as amino acids, lipids, and carbohydrates metabolisms was also detected, suggesting a lower oxidative stress condition in the AMF inoculated plants. Interestingly, many transcription factor families, including WRKY, NAC, and MYB, already known for their key role in plant abiotic stress response, were found differentially expressed between treatments. This study provides valuable insights on AMF-induced gene expression modulation and the beneficial effects of plant-AMF interaction in durum wheat under salt stress.

Genetic Diversity and Population Structure of Common Bean (Phaseolus vulgaris L.) Landraces in the Lazio Region of Italy.

Common bean cultivation has historically been a typical component of rural economies in Italy, particularly in mountainous and hilly zones along the Apennine ridge of the central and southern regions, where the production is focused on local landraces cultivated by small-scale farmers using low-input production systems. Such landraces are at risk of genetic erosion because of the recent socioeconomic changes in rural communities. One hundred fourteen accessions belonging to 66 landraces still being grown in the Lazio region were characterized using a multidisciplinary approach. This approach included morphological (seed traits), biochemical (phaseolin and phytohemagglutinin patterns), and molecular (microsatellite loci) analyses to investigate their genetic variation, structure, and distinctiveness, which will be essential for the implementation of adequate ex situ and in situ conservation strategies. Another objective of this study was to determine the original gene pool (Andean and Mesoamerican) of the investigated landraces and to evaluate the cross-hybridization events between the two ancestral gene pools in the P. vulgaris germplasm in the Lazio region. Molecular analyses on 456 samples (four for each of the 114 accessions) revealed that the P. vulgaris germplasm in the Lazio region exhibited a high level of genetic diversity (He = 0.622) and that the Mesoamerican and Andean gene pools were clearly differentiated, with the Andean gene pool prevailing (77%) and 12% of landraces representing putative hybrids between the two gene pools. A model-based cluster analysis based on the molecular markers highlighted three main groups in agreement with the phaseolin patterns and growth habit of landraces. The combined utilisation of morphological, biochemical, and molecular data allowed for the differentiation of all landraces and the resolution of certain instances of homonymy and synonymy. Furthermore, although a high level of homozygosity was found across all landraces, 32 of the 66 examined (49%) exhibited genetic variability, indicating that the analysis based on a single or few plants per landrace, as usually carried out, may provide incomplete information.

Bioactive Metabolite Survey of Actinobacteria Showing Plant Growth Promoting Traits to Develop Novel Biofertilizers.

The use of chemical fertilizers and pesticides has caused harmful impacts on the environment with the increase in economic burden. Biofertilizers are biological products containing living microorganisms capable of improving plant growth through eco-friendly mechanisms. In this work, three actinobacterial strains Streptomyces violaceoruber, Streptomyces coelicolor, and Kocuria rhizophila were characterized for multiple plant growth promoting (PGP) traits such as indole acetic acid production, phosphate solubilization, N2-fixation, and drought and salt tolerance. Then, these strains were investigated for their secreted and cellular metabolome, revealing a rich arsenal of bioactive molecules, including antibiotics and siderophores, with S. violaceoruber being the most prolific strain. Furthermore, the in vivo assays, performed on tomato (Solanum lycopersicum L.), resulted in an improved germination index and the growth of seedlings from seeds treated with PGP actinobacteria, with a particular focus on S. violaceoruber cultures. In particular, this last strain, producing volatile organic compounds having antimicrobial activity, was able to modulate volatilome and exert control on the global DNA methylation of tomato seedlings. Thus, these results, confirming the efficacy of the selected actinobacteria strains in promoting plant growth and development by producing volatile and non-volatile bioactive molecules, can promote eco-friendly alternatives in sustainable agriculture.

WRKY Gene Family Drives Dormancy Release in Onion Bulbs.

Onion (Allium cepa L.) is an important bulb crop grown worldwide. Dormancy in bulbous plants is an important physiological state mainly regulated by a complex gene network that determines a stop of vegetative growth during unfavorable seasons. Limited knowledge on the molecular mechanisms that regulate dormancy in onion were available until now. Here, a comparison between uninfected and onion yellow dwarf virus (OYDV)-infected onion bulbs highlighted an altered dormancy in the virus-infected plants, causing several symptoms, such as leaf striping, growth reduction, early bulb sprouting and rooting, as well as a lower abscisic acid (ABA) level at the start of dormancy. Furthermore, by comparing three dormancy stages, almost five thousand four hundred (5390) differentially expressed genes (DEGs) were found in uninfected bulbs, while the number of DEGs was significantly reduced (1322) in OYDV-infected bulbs. Genes involved in cell wall modification, proteolysis, and hormone signaling, such as ABA, gibberellins (GAs), indole-3-acetic acid (IAA), and brassinosteroids (BRs), that have already been reported as key dormancy-related pathways, were the most enriched ones in the healthy plants. Interestingly, several transcription factors (TFs) were up-regulated in the uninfected bulbs, among them three genes belonging to the WRKY family, for the first time characterized in onion, were identified during dormancy release. The involvement of specific WRKY genes in breaking dormancy in onion was confirmed by GO enrichment and network analysis, highlighting a correlation between AcWRKY32 and genes driving plant development, cell wall modification, and division via gibberellin and auxin homeostasis, two key processes in dormancy release. Overall, we present, for the first time, a detailed molecular analysis of the dormancy process, a description of the WRKY-TF family in onion, providing a better understanding of the role played by AcWRKY32 in the bulb dormancy release. The TF co-expressed genes may represent targets for controlling the early sprouting in onion, laying the foundations for novel breeding programs to improve shelf life and reduce postharvest.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades.

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

A Complex Gene Network Mediated by Ethylene Signal Transduction TFs Defines the Flower Induction and Differentiation in Olea europaea L.

The olive tree (Olea europaea L.) is a typical Mediterranean crop, important for olive and oil production. The high tendency to bear fruits in an uneven manner, defined as irregular or alternate bearing, results in a significant economic impact for the high losses in olives and oil production. Buds from heavy loaded ('ON') and unloaded ('OFF') branches of a unique olive tree were collected in July and the next March to compare the transcriptomic profiles and get deep insight into the molecular mechanisms regulating floral induction and differentiation. A wide set of DEGs related to ethylene TFs and to hormonal, sugar, and phenylpropanoid pathways was identified in buds collected from 'OFF' branches. These genes could directly and indirectly modulate different pathways, suggesting their key role during the lateral bud transition to flowering stage. Interestingly, several genes related to the flowering process appeared as over-expressed in buds from March 'OFF' branches and they could address the buds towards flower differentiation. By this approach, interesting candidate genes related to the switch from vegetative to reproductive stages were detected and analyzed. The functional analysis of these genes will provide tools for developing breeding programs to obtain olive trees characterized by more constant productivity over the years.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes.

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

Comparative Transcriptomics in Two Bivalve Species Offers Different Perspectives on the Evolution of Sex-Biased Genes.

Comparative genomics has become a central tool for evolutionary biology, and a better knowledge of understudied taxa represents the foundation for future work. In this study, we characterized the transcriptome of male and female mature gonads in the European clam Ruditapes decussatus, compared with that in the Manila clam Ruditapes philippinarum providing, for the first time in bivalves, information about transcription dynamics and sequence evolution of sex-biased genes. In both the species, we found a relatively low number of sex-biased genes (1,284, corresponding to 41.3% of the orthologous genes between the two species), probably due to the absence of sexual dimorphism, and the transcriptional bias is maintained in only 33% of the orthologs. The dN/dS is generally low, indicating purifying selection, with genes where the female-biased transcription is maintained between the two species showing a significantly higher dN/dS. Genes involved in embryo development, cell proliferation, and maintenance of genome stability show a faster sequence evolution. Finally, we report a lack of clear correlation between transcription level and evolutionary rate in these species, in contrast with studies that reported a negative correlation. We discuss such discrepancy and call into question some methodological approaches and rationales generally used in this type of comparative studies.

Burrowers from the Past: Mitochondrial Signatures of Ordovician Bivalve Infaunalization.

Bivalves and gastropods are the two largest classes of extant molluscs. Despite sharing a huge number of features, they do not share a key ecological one: gastropods are essentially epibenthic, although most bivalves are infaunal. However, this is not the ancestral bivalve condition; Cambrian forms were surface crawlers and only during the Ordovician a fundamental infaunalization process took place, leading to bivalves as we currently know them. This major ecological shift is linked to the exposure to a different redox environoments (hypoxic or anoxic) and with the Lower Devonian oxygenation event. We investigated selective signatures on bivalve and gastropod mitochondrial genomes with respect to a time calibrated mitochondrial phylogeny by means of dN/dS ratios. We were able to detect 1) a major signal of directional selection between the Ordovician and the Lower Devonian for bivalve mitochondrial Complex I, and 2) an overall higher directional selective pressure on bivalve Complex V with respect to gastropods. These and other minor dN/dS patterns and timings are discussed, showing that the Ordovician infaunalization event left heavy traces in bivalve mitochondrial genomes.

Comparative Large-Scale Mitogenomics Evidences Clade-Specific Evolutionary Trends in Mitochondrial DNAs of Bivalvia.

Despite the figure of complete bivalve mitochondrial genomes keeps growing, an assessment of the general features of these genomes in a phylogenetic framework is still lacking, despite the fact that bivalve mitochondrial genomes are unusual under different aspects. In this work, we constructed a dataset of one hundred mitochondrial genomes of bivalves to perform the first systematic comparative mitogenomic analysis, developing a phylogenetic background to scaffold the evolutionary history of the class' mitochondrial genomes. Highly conserved domains were identified in all protein coding genes; however, four genes (namely, atp6, nad2, nad4L, and nad6) were found to be very divergent for many respects, notwithstanding the overall purifying selection working on those genomes. Moreover, the atp8 gene was newly annotated in 20 mitochondrial genomes, where it was previously declared as lacking or only signaled. Supernumerary mitochondrial proteins were compared, but it was possible to find homologies only among strictly related species. The rearrangement rate on the molecule is too high to be used as a phylogenetic marker, but here we demonstrate for the first time in mollusks that there is correlation between rearrangement rates and evolutionary rates. We also developed a new index (HERMES) to estimate the amount of mitochondrial evolution. Many genomic features are phylogenetically congruent and this allowed us to highlight three main phases in bivalve history: the origin, the branching of palaeoheterodonts, and the second radiation leading to the present-day biodiversity.