Highly divergent satellitomes of two barley species of agronomic importance, Hordeum chilense and H. vulgare.

In this paper, we have performed an in-depth study of the complete set of the satellite DNA (satDNA) families (i.e. the satellitomes) in the genome of two barley species of agronomic value in a breeding framework, H. chilense (H1 and H7 accessions) and H. vulgare (H106 accession), which can be useful tools for studying chromosome associations during meiosis. The study has led to the analysis of a total of 18 satDNA families in H. vulgare, 25 satDNA families in H. chilense (accession H1) and 27 satDNA families in H. chilense (accession H7) that constitute 46 different satDNA families forming 36 homology groups. Our study highlights different important contributions of evolutionary and applied interests. Thus, both barley species show very divergent satDNA profiles, which could be partly explained by the differential effects of domestication versus wildlife. Divergence derives from the differential amplification of different common ancestral satellites and the emergence of new satellites in H. chilense, usually from pre-existing ones but also random sequences. There are also differences between the two H. chilense accessions, which support genetically distinct groups. The fluorescence in situ hybridization (FISH) patterns of some satDNAs yield distinctive genetic markers for the identification of specific H. chilense or H. vulgare chromosomes. Some of the satellites have peculiar structures or are related to transposable elements which provide information about their origin and expansion. Among these, we discuss the existence of different (peri)centromeric satellites that supply this region with some plasticity important for centromere evolution. These peri(centromeric) satDNAs and the set of subtelomeric satDNAs (a total of 38 different families) are analyzed in the framework of breeding as the high diversity found in the subtelomeric regions might support their putative implication in chromosome recognition and pairing during meiosis, a key point in the production of addition/substitution lines and hybrids.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria.

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

Five centuries of genome evolution and multi-host adaptation of Campylobacter jejuni in Brazil.

Consumption of raw, undercooked or contaminated animal food products is a frequent cause of Campylobacter jejuni infection. Brazil is the world's third largest producer and a major exporter of chicken meat, yet population-level genomic investigations of C. jejuni in the country remain scarce. Analysis of 221 C. jejuni genomes from Brazil shows that the overall core and accessory genomic features of C. jejuni are influenced by the identity of the human or animal source. Of the 60 sequence types detected, ST353 is the most prevalent and consists of samples from chicken and human sources. Notably, we identified the presence of diverse bla genes from the OXA-61 and OXA-184 families that confer beta-lactam resistance as well as the operon cmeABCR related to multidrug efflux pump, which contributes to resistance against tetracyclines, macrolides and quinolones. Based on limited data, we estimated the most recent common ancestor of ST353 to the late 1500s, coinciding with the time the Portuguese first arrived in Brazil and introduced domesticated chickens into the country. We identified at least two instances of ancestral chicken-to-human infections in ST353. The evolution of C. jejuni in Brazil was driven by the confluence of clinically relevant genetic elements, multi-host adaptation and clonal population growth that coincided with major socio-economic changes in poultry farming.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation.

BACKGROUND: Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family. FINDINGS: The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species. CONCLUSIONS: The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Comparative analysis of transposable elements dynamics in fish with different sex chromosome systems.

Transposable elements (TEs) are widespread genomic components with substantial roles in genome evolution and sex chromosome differentiation. In this study, we compared the TE composition of three closely related fish with different sex chromosome systems: Megaleporinus elongatus (Z1Z1Z2Z2/Z1W1Z2W2), Megaleporinus macrocephalus (ZZ/ZW) (both with highly differentiated W sex chromosomes), and Leporinus friderici (without heteromorphic sex chromosomes). We created custom TE libraries for each species using clustering methods and manual annotation and prediction, and we predicted TE temporal dynamics through divergence-based analysis. The TE abundance ranged from 16% to 21% in the three mobilomes, with L. friderici having the lowest overall. Despite the recent amplification of TEs in all three species, we observed differing expansion activities, particularly between the two genera. Both Megaleporinus recently experienced high retrotransposon activity, with a reduction in DNA TEs, which could have implications in sex chromosome composition. In contrast, L. friderici showed the opposite pattern. Therefore, despite having similar TE compositions, Megaleporinus and Leporinus exhibit distinct TE histories that likely evolved after their separation, highlighting a rapid TE expansion over short evolutionary periods.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria.

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Phylogenomic discernments into Anaerolineaceae thermal adaptations and the proposal of a candidate genus Mesolinea.

This study delves into the evolutionary history of Anaerolineaceae, a diverse bacterial family within the Chloroflexota phylum. Employing a multi-faceted approach, including phylogenetic analyses, genomic comparisons, and exploration of adaptive features, the research unveils novel insights into the family's taxonomy and evolutionary dynamics. The investigation employs metagenome-assembled genomes (MAGs), emphasizing their prevalence in anaerobic environments. Notably, a novel mesophilic lineage, tentatively named Mesolinea, emerges within Anaerolineaceae, showcasing a distinctive genomic profile and apparent adaptation to a mesophilic lifestyle. The comprehensive genomic analyses shed light on the family's complex evolutionary patterns, including the conservation of key operons in thermophiles, providing a foundation for understanding the diverse ecological roles and adaptive strategies of Anaerolineaceae members.

An evolutionary view of the Fusarium core genome.

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

A rapid and affordable amplicon-based method for next-generation genome sequencing of the infectious bursal disease virus.

The infectious bursal disease virus (IBDV) causes a severe immunosuppressive disorder in young chickens. IBDV evolution resulted in the emergence of strains with divergent genetic, antigenic, and pathogenic characteristics. Genetic classification is typically performed by sequencing the coding region of the most immunogenic region of the viral protein 2 (VP2). Sequencing both double-stranded RNA genome segments is essential to achieve a more comprehensive IBDV classification that can detect recombinants and reassortments. Here, we report the development and standardization of a tiled PCR amplicon protocol for the direct and cost-effective genome sequencing of global IBDV strains using next-generation technology. Primers for tiled PCR were designed with adapters to bypass expensive and time-consuming library preparation steps. Sequencing was performed on Illumina MiniSeq equipment, and fourteen complete genomes of field strains were assembled using reference sequences. The PCR-enrichment step was used to obtain genomes from low-titer biological samples that were difficult to amplify using traditional sequencing. Phylogenetic analyses of the obtained genomes confirmed previous strain classification. By combining the enrichment methodology with massive sequencing, it is possible to obtain IBDV genomic sequences in a fast and affordable manner. This procedure can be a valuable tool to better understand virus epidemiology.

Clonal reproduction of Moniliophthora roreri and the emergence of unique lineages with distinct genomes during range expansion.

The basidiomycete Moniliophthora roreri causes frosty pod rot of cacao (Theobroma cacao) in the western hemisphere. Moniliophthora roreri is considered asexual and haploid throughout its hemibiotrophic life cycle. To understand the processes driving genome modification, using long-read sequencing technology, we sequenced and assembled 5 high-quality M. roreri genomes out of a collection of 99 isolates collected throughout the pathogen's range. We obtained chromosome-scale assemblies composed of 11 scaffolds. We used short-read technology to sequence the genomes of 22 similarly chosen isolates. Alignments among the 5 reference assemblies revealed inversions, translocations, and duplications between and within scaffolds. Isolates at the front of the pathogens' expanding range tend to share lineage-specific structural variants, as confirmed by short-read sequencing. We identified, for the first time, 3 new mating type A locus alleles (5 in total) and 1 new potential mating type B locus allele (3 in total). Currently, only 2 mating type combinations, A1B1 and A2B2, are known to exist outside of Colombia. A systematic survey of the M. roreri transcriptome across 2 isolates identified an expanded candidate effector pool and provided evidence that effector candidate genes unique to the Moniliophthoras are preferentially expressed during the biotrophic phase of disease. Notably, M. roreri isolates in Costa Rica carry a chromosome segment duplication that has doubled the associated gene complement and includes secreted proteins and candidate effectors. Clonal reproduction of the haploid M. roreri genome has allowed lineages with unique genome structures and compositions to dominate as it expands its range, displaying a significant founder effect.

The complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. Tribe (Myrtaceae).

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral evolution was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

In Silico Identification and Characterization of Satellite DNAs in 23 Drosophila Species from the Montium Group.

Satellite DNA (satDNA) is a class of tandemly repeated non-protein coding DNA sequences which can be found in abundance in eukaryotic genomes. They can be functional, impact the genomic architecture in many ways, and their rapid evolution has consequences for species diversification. We took advantage of the recent availability of sequenced genomes from 23 Drosophila species from the montium group to study their satDNA landscape. For this purpose, we used publicly available whole-genome sequencing Illumina reads and the TAREAN (tandem repeat analyzer) pipeline. We provide the characterization of 101 non-homologous satDNA families in this group, 93 of which are described here for the first time. Their repeat units vary in size from 4 bp to 1897 bp, but most satDNAs show repeat units < 100 bp long and, among them, repeats ≤ 10 bp are the most frequent ones. The genomic contribution of the satDNAs ranges from ~1.4% to 21.6%. There is no significant correlation between satDNA content and genome sizes in the 23 species. We also found that at least one satDNA originated from an expansion of the central tandem repeats (CTRs) present inside a Helitron transposon. Finally, some satDNAs may be useful as taxonomic markers for the identification of species or subgroups within the group.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome.

BACKGROUND: B chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome. RESULTS: Through small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome. CONCLUSIONS: Due to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

Translation Comes First: Ancient and Convergent Selection of Codon Usage Bias Across Prokaryotic Genomes.

Codon usage is the outcome of different evolutionary processes and can inform us about the conditions in which organisms live and evolve. Here, we present R_ENC', which is an improvement to the original S index developed by dos Reis et al. (2004). Our index is less sensitive to G+C content, which greatly affects synonymous codon usage in prokaryotes, making it better suited to detect selection acting on codon usage. We used R_ENC' to estimate the extent of selected codon usage bias in 1800 genomes representing 26 prokaryotic phyla. We found that Gammaproteobacteria, Betaproteobacteria, Actinobacteria, and Firmicutes are the phyla/subphyla showing more genomes with selected codon usage bias. In particular, we found that several lineages within Gammaproteobacteria and Firmicutes show a similar set of functional terms enriched in genes under selected codon usage bias, indicating convergent evolution. We also show that selected codon usage bias tends to evolve in genes coding for the translation machinery before other functional GO terms. Finally, we discuss the possibility to use R_ENC' to predict whether lineages evolved in copiotrophic or oligotrophic environments.

Genomic analysis reveals the role of integrative and conjugative elements in plant pathogenic bacteria.

BACKGROUND: ICEs are mobile genetic elements found integrated into bacterial chromosomes that can excise and be transferred to a new cell. They play an important role in horizontal gene transmission and carry accessory genes that may provide interesting phenotypes for the bacteria. Here, we seek to research the presence and the role of ICEs in 300 genomes of phytopathogenic bacteria with the greatest scientific and economic impact. RESULTS: Seventy-eight ICEs (45 distinct elements) were identified and characterized in chromosomes of Agrobacterium tumefaciens, Dickeya dadantii, and D. solani, Pectobacterium carotovorum and P. atrosepticum, Pseudomonas syringae, Ralstonia solanacearum Species Complex, and Xanthomonas campestris. Intriguingly, the co-occurrence of four ICEs was observed in some P. syringae strains. Moreover, we identified 31 novel elements, carrying 396 accessory genes with potential influence on virulence and fitness, such as genes coding for functions related to T3SS, cell wall degradation and resistance to heavy metals. We also present the analysis of previously reported data on the expression of cargo genes related to the virulence of P. atrosepticum ICEs, which evidences the role of these genes in the infection process of tobacco plants. CONCLUSIONS: Altogether, this paper has highlighted the potential of ICEs to affect the pathogenicity and lifestyle of these phytopathogens and direct the spread of significant putative virulence genes in phytopathogenic bacteria.

What We Talk About When We Talk About "Junk DNA".

"Junk DNA" is a popular yet controversial concept that states that organisms carry in their genomes DNA that has no positive impact on their fitness. Nonetheless, biochemical functions have been identified for an increasing fraction of DNA elements traditionally seen as "Junk DNA". These findings have been interpreted as fundamentally undermining the "Junk DNA" concept. Here, we reinforce previous arguments that this interpretation relies on an inadequate concept of biological function that does not consider the selected effect of a given genomic structure, which is central to the "Junk DNA" concept. Next, we suggest that another (though ignored) confounding factor is that the discussion about biological functions includes two different dimensions: a horizontal, ecological dimension that reflects how a given genomic element affects fitness in a specific time, and a vertical, temporal dimension that reflects how a given genomic element persisted along time. We suggest that "Junk DNA" should be used exclusively relative to the horizontal dimension, while for the vertical dimension, we propose a new term, "Spam DNA", that reflects the fact that a given genomic element may persist in the genome even if not selected for on their origin. Importantly, these concepts are complementary. An element can be both "Spam DNA" and "Junk DNA", and "Spam DNA" can also be recruited to perform evolved biological functions, as illustrated in processes of exaptation or constructive neutral evolution.

New insights into plant natriuretic peptide evolution: From the lysogenic conversion in Xanthomonas to the lateral transfer to the whitefly Bemisia tabaci.

Plant natriuretic peptide-like (PNP) are signaling molecules related to adaptive responses to stress. The Arabidopsis thaliana PNP (AtPNP-A) is capable of modulating catalase 2 (CAT2) and rubisco activase (RCA) activity in some circumstances. Interestingly, many plant-pathogens co-opted PNP-like molecules to their benefit. For instance, the citrus pathogen Xanthomonas citri carries a PNP-like (XacPNP) that can mimic and regulate plant homeostasis, and many phytopathogenic fungi carry effectors (e.g., Ave1 and AvrLm6) that are indeed PNP-like homologs. This work investigates the PNP-like evolution across the tree of life, revealing many parallel gains and duplications in plant and fungi kingdoms. All PNP-like proteins in the final dataset are structurally similar, containing the AtPNP-A active domains modulating CAT2 activity and RCA interaction. Comparative genomics evinced that XacPNP is a lysogenic conversion factor associated with a Myoviridae-like prophage identified in many Xanthomonas species. Surprisingly, a PNP-like homolog was identified in Bemisia tabaci, an important agricultural pest, being to date the second example of lateral gene transfer (LGT) from plant to the whitefly. Moreover, the Bemisia PNP-like homolog can also be considered a potential new effector of this phloem-feeding insect. Noteworthy, the whiteflies infest many plants carrying PNP-like copies and interact with some of their bacterial and fungal pathogens, strongly suggesting complex recipient/donor traits of PNP by LGT and bringing new insights into the evolution of host-pathogen arms race across the tree of life.

Dynamic development of AT-rich heterochromatin has followed diversification and genome expansion of psammophilous Mycetophylax (Formicidae: Attini: Attina).

Heterochromatin is an important genome constituent comprised by a high density of repetitive DNA sequences that mediate chromosome structure and function. The species Mycetophylax morschi currently harbours three cytotypes: 2n = 26, 2n = 28 and 2n = 30 chromosomes. However, Mycetophylax conformis and Mycetophylax simplex harbour 2n = 30 and 2n = 36 chromosomes, respectively. None of the cytotypes of M. morschi showed any AT-positive blocks, whereas the karyotypes of M. conformis and M. simplex revealed AT-rich blocks around the pericentromeric region and on the short arm of several chromosomes. This AT-rich pattern is coincident with the known heterochromatin distribution of psammophilous Mycetophylax, confirming that heterochromatin is AT-rich, in line with the genome size and AT%. Our results demonstrated that genome size among psammophilous Mycetophylax is correlated with the proportion of base pairs, biased to adenine and thymine. Thus, genome size and the proportion of adenine and thymine in the species studied here suggest that the genome changes in psammophilous Mycetophylax are related to the expansion of repetitive DNA in AT-rich heterochromatin. Considering the phylogenetic relationship of psammophilous Mycetophylax, the dynamic development of AT-rich heterochromatin and karyotype repatterning encompasses the diversification of such ants.