A high-resolution comparative atlas across 74 fish genomes illuminates teleost evolution after whole-genome duplication

Teleost fish are one of the most species-rich and diverse clades amongst vertebrates, which makes them an outstanding model group for evolutionary, ecological and functional genomics. Yet, despite a growing number of sequence reference genomes, large-scale comparative analysis remains challenging in teleosts due to the specifics of their genomic organization. As legacy of a whole genome duplication dated 320 million years ago, a large fraction of teleost genomes remain in duplicate paralogous copies. This ancestral polyploidy confounds the detailed identification of orthologous genomic regions across teleost species. Here, we combine tailored gene phylogeny methodology together with the state-of-the art ancestral karyotype reconstruction to establish the first high resolution comparative atlas of paleopolyploid regions across 74 teleost fish genomes. We show that this atlas represents a unique, robust and reliable resource for fish genomics. We then use the comparative atlas to study the tetraploidization and rediploidization mechanisms that affected the ancestor of teleosts. Although the polyploid history of teleost genomes appears complex, we uncover that meiotic recombination persisted between duplicated chromosomes for over 60 million years after polyploidization, suggesting that the teleost ancestor was an autotetraploid.

Whole-Genome Duplication and Host Genotype Affect Rhizosphere Microbial Communities

Plants influence the composition of their associated microbial communities, yet the underlying host-associated genetic determinants are typically unknown. Genome duplication events are common in the evolutionary history of plants and affect many plant traits.

The Preservation of PPARγ Genome Duplicates in Some Teleost Lineages: Insights into Lipid Metabolism and Xenobiotic Exploitation

Three peroxisome proliferator-activated receptor paralogues (PPARα, -β and -γ) are currently recognized in vertebrate genomes. PPARγ is known to modulate nutrition, adipogenesis and immunity in vertebrates. Natural ligands of PPARγ have been proposed; however, the receptor also binds synthetic ligands such as endocrine disruptors. Two paralogues of PPARα and PPARβ have been documented in teleost species, a consequence of the 3R WGD. Recently, two PPARγ paralogue genes were also identified in Astyanax mexicanus. We aimed to determine whether the presence of two PPARγ paralogues is prevalent in other teleost genomes, through genomic and phylogenetic analysis. Our results showed that besides Characiformes, two PPARγ paralogous genes were also identified in other teleost taxa, coinciding with the teleost-specific, whole-genome duplication and with the retention of both genes prior to the separation of the Clupeocephala. To functionally characterize these genes, we used the European sardine (Sardina pilchardus) as a model. PPARγA and PPARγB display a different tissue distribution, despite the similarity of their functional profiles: they are unresponsive to tested fatty acids and other human PPARγ ligands yet yield a transcriptional response in the presence of tributyltin (TBT). This observation puts forward the relevance of comparative analysis to decipher alternative binding architectures and broadens the disruptive potential of man-made chemicals for aquatic species.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

AnchorWave: Sensitive alignment of genomes with high sequence diversity, extensive structural polymorphism, and whole-genome duplication

Millions of species are currently being sequenced, and their genomes are being compared. Many of them have more complex genomes than model systems and raise novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences. Moreover, alignment using many-to-many or reciprocal best hit approaches conflicts with well-studied patterns between species with different rounds of whole-genome duplication. Here, we introduce Anchored Wavefront alignment (AnchorWave), which performs whole-genome duplication–informed collinear anchor identification between genomes and performs base pair–resolved global alignment for collinear blocks using a two-piece affine gap cost strategy. This strategy enables AnchorWave to precisely identify multikilobase indels generated by transposable element (TE) presence/absence variants (PAVs). When aligning two maize genomes, AnchorWave successfully recalled 87% of previously reported TE PAVs. By contrast, other genome alignment tools showed low power for TE PAV recall. AnchorWave precisely aligns up to three times more of the genome as position matches or indels than the closest competitive approach when comparing diverse genomes. Moreover, AnchorWave recalls transcription factor–binding sites at a rate of 1.05- to 74.85-fold higher than other tools with significantly lower false-positive alignments. AnchorWave complements available genome alignment tools by showing obvious improvement when applied to genomes with dispersed repeats, active TEs, high sequence diversity, and whole-genome duplication variation.

Evolution of the Parvalbumin Genes in Teleost Fishes after the Whole-Genome Duplication

Parvalbumin is considered a major fish allergen. Here, we report the molecular evolution of the parvalbumin genes in bony fishes based on 19 whole genomes and 70 transcriptomes. We found unexpectedly high parvalbumin diversity in teleosts; three main gene types (pvalb-α, pvalb-β1, and pvalb-β2, including oncomodulins) originated at the onset of vertebrates. Teleosts have further multiplied the parvalbumin gene repertoire up to nine ancestral copies—two copies of pvalb-α, two copies of pvalb-β1, and five copies of pvalb-β2. This gene diversity is a result of teleost-specific whole-genome duplication. Two conserved parvalbumin genomic clusters carry pvalb-β1 and β2 copies, whereas pvalb-α genes are located separately in different linkage groups. Further, we investigated parvalbumin gene expression in 17 tissues of the common carp (Cyprinus carpio), a species with 21 parvalbumin genes in its genome. Two pvalb-α and eight pvalb-β2 copies are highly expressed in the muscle, while two alternative pvalb-α copies show expression in the brain and the testes, and pvalb-β1 is dominant in the retina and the kidney. The recent pairs of muscular pvalb-β2 genes show differential expression in this species. We provide robust genomic evidence of the complex evolution of the parvalbumin genes in fishes.

Integrating transcriptomes to investigate genes associated with adaptation to aquatic environments, and assess phylogenetic conflict and whole-genome duplications in Alismatales

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times beyond. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptation to aquatic environments by including 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptation, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as ((Tofieldiaceae, Araceae), core Alismatids). There is phylogenetic conflict among the backbone of the three main clades that could be due to incomplete lineage sorting and introgression. We identified 18 putative WGD events. One of them had occurred at the most recent common ancestor of core Alismatids, and four occurred at seagrass lineages. Other events are distributed in terrestrial, emergent, and submersed life-forms and seagrasses across Alismatales. We also found that lineage and life-form were each important for different evolutionary patterns for the genes related to freshwater/marine adaptation. For example, some light or ethylene-related genes were lost in the seagrass Zosteraceae, but present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments, phylogeny, and whole-genome duplication of Alismatales.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L.

Backgrounds Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae. Results In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists. Conclusions The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.