Genome Biology and Evolution最新文献

Recent advances in high-throughput sequencing have exponentially increased the number of genomic data available for animals (Metazoa) in the last decades, with high-quality chromosome-level genomes being published almost daily. Nevertheless, generating a new genome is not an easy task due to the high cost of genome sequencing, the high complexity of assembly, and the lack of standardized protocols for genome annotation. The lack of consensus in the annotation and publication of genome files hinders research by making researchers lose time in reformatting the files for their purposes but can also reduce the quality of the genetic repertoire for an evolutionary study. Thus, the use of transcriptomes obtained using the same pipeline as a proxy for the genetic content of species remains a valuable resource that is easier to obtain, cheaper, and more comparable than genomes. In a previous study, we presented the Metazoan Assemblies from Transcriptomic Ensembles database (MATEdb), a repository of high-quality transcriptomic and genomic data for the two most diverse animal phyla, Arthropoda and Mollusca. Here, we present the newest version of MATEdb (MATEdb2) that overcomes some of the previous limitations of our database: (i) we include data from all animal phyla where public data are available, and (ii) we provide gene annotations extracted from the original GFF genome files using the same pipeline. In total, we provide proteomes inferred from high-quality transcriptomic or genomic data for almost 1,000 animal species, including the longest isoforms, all isoforms, and functional annotation based on sequence homology and protein language models, as well as the embedding representations of the sequences. We believe this new version of MATEdb will accelerate research on animal phylogenomics while saving thousands of hours of computational work in a plea for open, greener, and collaborative science.

Fungi are well-known producers of bioactive secondary metabolites (SMs), which have been exploited for decades by humankind for various medical applications like therapeutics and antibiotics. SMs are synthesized by biosynthetic gene clusters (BGCs)-physically co-localized and co-regulated genes. Because BGCs are often regulated by histone post-translational modifications (PTMs), it was suggested that their chromosomal location is important for their expression. Studies in a few fungal species indicated an enrichment of BGCs in sub-telomeric regions; however, there is no evidence that BGCs with distinct genomic localization are regulated by different histone PTMs. Here, we used 174 Aspergillus species covering 22 sections to determine the correlation between BGC genomic localization, gene expression, and histone PTMs. We found a high abundance and diversity of SM backbone genes across the Aspergillus genus, with notable unique genes within sections. Being unique or conserved in many species, BGCs showed a strong bias for being localized in low-synteny regions, regardless of their position in chromosomes. Using chromosome-level assemblies, we also confirmed a significantly biased localization in sub-telomeric regions. Notably, SM backbone genes in sub-telomeric regions and about half of those in low-synteny regions exhibit higher gene expression variability, likely due to the similar higher variability in H3K4me3 and H3K36me3 histone PTMs; while variations in histone H3 acetylation and H3K9me3 are not correlated to genomic localization and expression variation, as analyzed in two Aspergillus species. Expression variability across four Aspergillus species further supports that BGCs tend to be located in low-synteny regions and that regulation of expression in those regions likely involves different histone PTMs than the most commonly studied modifications.

The origins and early evolution of animals are subjects with many outstanding questions. One problem faced by researchers trying to answer them is the absence of a comprehensive database with sequences from nonbilaterians. Publicly available data are plentiful but scattered and often not associated with proper metadata. A new database presented in this paper, LukProt, is an attempt at solving this issue. The database contains protein sequences obtained mostly from genomic, transcriptomic, and metagenomic studies and is an extension of EukProt (Richter DJ, Berney C, Strassert JFH, Poh Y-P, Herman EK, Muñoz-Gómez SA, Wideman JG, Burki F, de Vargas C. EukProt: a database of genome-scale predicted proteins across the diversity of eukaryotes. Peer Community J. 2022:2:e56. https://doi.org/10.24072/pcjournal.173). LukProt adopts the EukProt naming conventions and includes data from 216 additional animals. The database is associated with a taxonomic grouping (taxogroup) scheme suitable for studying early animal evolution. Minor updates to the database will contain species additions or metadata corrections, whereas major updates will synchronize LukProt to each new version of EukProt, and releases are permanently stored on Zenodo (https://doi.org/10.5281/zenodo.7089120). A BLAST server to search the database is available at: https://lukprot.hirszfeld.pl/. Users are invited to participate in maintaining and correcting LukProt. As it can be searched without downloading locally, the database aims to be a convenient resource not only for evolutionary biologists, but for the broader scientific community as well.

When environmental change is rapid or unpredictable, phenotypic plasticity can facilitate adaptation to new or stressful environments to promote population persistence long enough for adaptive evolution to occur. However, the underlying genetic mechanisms that contribute to plasticity and its role in adaptive evolution are generally unknown. Two main opposing hypotheses dominate-genetic compensation and genetic assimilation. Here, we predominantly find evidence for genetic compensation over assimilation in adapting the freshwater algae Chlamydomonas reinhardtii to 36 g/L salt environments over 500 generations. More canalized genes in the high-salt (HS) lines displayed a pattern of genetic compensation (63%) fixing near or at the ancestral native expression level, rather than genetic assimilation of the salt-induced level, suggesting that compensation was more common during adaptation to salt. Network analysis revealed an enrichment of genes involved in energy production and salt-resistance processes in HS lines, while an increase in DNA repair mechanisms was seen in ancestral strains. In addition, whole-transcriptome similarity among ancestral and HS lines displayed the evolution of a similar plastic response to salt conditions in independently reared HS lines. We also found more cis-acting regions in the HS lines; however, the expression patterns of most genes did not mimic that of their inherited sequence. Thus, the expression changes induced via plasticity offer temporary relief, but downstream changes are required for a sustainable solution during the evolutionary process.

The Antarctic sea urchin Sterechinus neumayeri (Echinoida; Echinidae) is routinely used as a model organism for Antarctic biology. Here, we present a high-quality genome of S. neumayeri. This chromosomal-level assembly was generated using PacBio long-read sequencing and Hi-C chromatin conformation capture sequencing. This 885.3-Mb assembly exhibits high contiguity with a scaffold length N50 of 36.7 Mb assembled into 20 chromosomal length scaffolds. These putative chromosomes exhibit a high degree of synteny compared to other sea urchin models. We used transcript evidence gene modeling combined with sequence homology to identify 21,638 gene models that capture 97.4% of BUSCO orthologs. Among these, we were able to identify and annotate conserved developmental gene regulatory network orthologs, positioning S. neumayeri as a tractable model for comparative studies on evolution and development.

Urbanization is a leading factor effecting global biodiversity, driving rapid evolutionary processes in the local biota. Species that adapt and proliferate in city environments can become pests, with human activities facilitating their dispersal and excessive outbreaks. Here we present the first genome data of Plecia longiforceps, a lovebug pest in Eastern Asia with intensive aggregations recently occurring in the Seoul Metropolitan Area of Korea. PacBio HiFi and ONT Pore-C sequencing data were used to construct a highly continuous assembly with a total size of 707 Mb and 8 major pseudochromosomes, its integrity supported by the N50 length of 98.1 Mb and 96.8% BUSCO completeness. Structural and functional annotation using transcriptome data and ab initio predictions revealed a high proportion (69.3%) of repeat sequences, and synteny analysis with Bibio marci showed high levels of genomic collinearity. The genome will serve as an essential resource for both population genomics and molecular research on lovebug dispersal and outbreaks, and also implement studies on the eco-evolutionary processes of insects in urbanizing habitats.

Bacteria lose and gain repair genes as they evolve. Here, we investigate the consequences of gain and loss of 11 DNA repair genes across a broad range of bacteria. Using synonymous polymorphisms from bacteria and a set of 50 phylogenetically independent contrasts, we find no evidence that the presence or absence of these 11 genes affects either the overall level of diversity or the pattern of mutation. Using phylogenetic generalized linear squares yields a similar conclusion. It seems likely that the lack of an effect is due to variation in the genetic background and the environment which obscures any effects that the presence or absence of individual genes might have.

Genomic regions that play a role in parasite defense are often found to be highly variable, with the major histocompatibility complex serving as an iconic example. Single nucleotide polymorphisms may represent only a small portion of this variability, with Indel polymorphisms and copy number variation further contributing. In extreme cases, haplotypes may no longer be recognized as orthologous. Understanding the evolution of such highly divergent regions is challenging because the most extreme variation is not visible using reference-assisted genomic approaches. Here we analyze the case of the Pasteuria Resistance Complex in the crustacean Daphnia magna, a defense complex in the host against the common and virulent bacterium Pasteuria ramosa. Two haplotypes of this region have been previously described, with parts of it being nonhomologous, and the region has been shown to be under balancing selection. Using pan-genome analysis and tree reconciliation methods to explore the evolution of the Pasteuria Resistance Complex and its characteristics within and between species of Daphnia and other Cladoceran species, our analysis revealed a remarkable diversity in this region even among host species, with many nonhomologous hyper-divergent haplotypes. The Pasteuria Resistance Complex is characterized by extensive duplication and losses of Fucosyltransferase (FuT) and Galactosyltransferase (GalT) genes that are believed to play a role in parasite defense. The Pasteuria Resistance Complex region can be traced back to common ancestors over 250 million years. The unique combination of an ancient resistance complex and a dynamic, hyper-divergent genomic environment presents a fascinating opportunity to investigate the role of such regions in the evolution and long-term maintenance of resistance polymorphisms. Our findings offer valuable insights into the evolutionary forces shaping disease resistance and adaptation, not only in the genus Daphnia, but potentially across the entire Cladocera class.

Molecular evolution analysis typically involves identifying selection pressure and reconstructing evolutionary trends. This process usually requires access to specific data related to a target gene or gene family within a particular group of organisms. While recent advancements in high-throughput sequencing techniques have resulted in the rapid accumulation of extensive genomics and transcriptomics data and the creation of new databases in public repositories, extracting valuable insights from such vast data sets remains a significant challenge for researchers. Here, we elucidated the evolutionary history of THI1, a gene responsible for encoding thiamine thiazole synthase. The thiazole ring is a precursor for vitamin B1 and a crucial cofactor in primary metabolic pathways. A thorough search of complete genomes available within public repositories reveals 702 THI1 homologs of Archaea and Eukarya. Throughout its diversification, the plant lineage has preserved the THI1 gene by incorporating the N-terminus and targeting the chloroplasts. Likewise, evolutionary pressures and lifestyle appear to be associated with retention of TPP riboswitch sites and consequent dual posttranscriptional regulation of the de novo biosynthesis pathway in basal groups. Multicopy retention of THI1 is not a typical plant pattern, even after successive genome duplications. Examining cis-regulatory sites in plants uncovers two shared motifs across all plant lineages. A data mining of 484 transcriptome data sets supports the THI1 homolog expression under a light/dark cycle response and a tissue-specific pattern. Finally, the work presented brings a new look at public repositories as an opportunity to explore evolutionary trends to THI1.

The multigene family of the major histocompatibility complex (MHC) codes for the key antigen-presenting molecules of the vertebrate immune system. In birds, duplicated MHC class II (MHC-II) genes are highly homogenized by concerted evolution, and thus, identification of their orthologous relationships across long evolutionary timescales remains challenging. Relatively low evolutionary rate of avian MHC class IIA genes has been expected to provide a promising avenue to allow such inferences, but availability of MHC-IIA sequences in nonmodel bird species has been limited until recently. Here, taking advantage from accumulating genomic resources, we identified and analyzed MHC-IIA sequences from the most basal lineage of extant birds (Palaeognathae). Conserved region of the MHC-IIA membrane-proximal domain was used to search for orthologous relationships between palaeognath birds and nonavian reptiles. First, analyses of palaeognath sequences revealed the presence of a separate MHC-IIA gene lineage (DAA3) in kiwis, which did not cluster with previously described avian MHC-IIA lineages (DAA1 and DAA2). Next, phylogenetic reconstruction showed that kiwi DAA3 sequences form a single well-supported cluster with turtle MHC-IIA. High similarity of these sequences most likely reflects their remarkable evolutionary conservation and retention of ancient orthologous relationships, which can be traced back to basal archosauromorphs ca. 250 million years ago. Our analyses offer novel insights into macroevolutionary history of the MHC and reinforce the view that rapid accumulation of high-quality genome assemblies across divergent nonmodel species can substantially advance our understanding of gene evolution.