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.

Gdnf, a germ cell-derived factor, regulates zebrafish germ cell stemness through the creation of new spermatogonial niches (germ and Sertoli cells) and inhibition of spermatogonial differentiation in an autocrine and paracrine manners

Glial cell line-derived neurotrophic factor (GDNF) and its receptor (GDNF Family Receptor α1 - GFRα1) are well known to mediate spermatogonial stem cell (SSC) proliferation and survival in the mammalian testes. In nonmammalian species, Gdnf and Gfrα1 orthologs have been found but their functions remain poorly investigated in the testis. Considering this background, this study aimed to understand the roles of Gdnf-Gfrα1 signaling pathway in the zebrafish testis by combining in vivo, in silico and ex vivo approaches. Our analysis showed that zebrafish exhibited two paralogs of Gndf (gdnfa and gdnfb) and its receptor, Gfrα1 (gfrα1a and gfrα1b), in agreement with the teleost-specific third round (3R) of whole genome duplication. Expression analysis further revealed that gdnfa and gfrα1a were the most expressed copies in the zebrafish adult testes. Subsequently, we demonstrated that gdnfa is expressed in the germ cells, while Gfrα1a was detected in early spermatogonia (mainly in types Aund and Adiff) and Sertoli cells. Functional ex vivo analysis showed that Gdnf promoted the creation of new available niches by stimulating proliferation of both type Aund spermatogonia and their surrounding Sertoli cells, but without changing pou5f3 mRNA levels. Strikingly, Gdnf also inhibited late spermatogonial differentiation as shown by the decrease of type B spermatogonia and down-regulation of dazl in the co-treatment with Fsh. Altogether, our data revealed for the first time that a germ cell-derived factor is associated with maintaining germ cell stemness through the creation of new available niches, supporting development of differentiating spermatogonial cysts and inhibiting late spermatogonial differentiation in autocrine and paracrine manners.

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.

B Chromosomes in Free-Living Flatworms of the Genus Macrostomum (Platyhelminthes, Macrostomorpha)

B chromosomes (Bs) or supernumerary chromosomes are extra chromosomes in the species karyotype that can vary in its copy number. Bs are widespread in eukaryotes. Usually, the Bs of specimens collected from natural populations are the object of the B chromosome studies. We applied another approach analyzing the Bs in animals maintained under the laboratory conditions as lines and cultures. In this study, three species of the Macrostomum genus that underwent a recent whole-genome duplication (WGD) were involved. In laboratory lines of M. lignano and M. janickei, the frequency of Bs was less than 1%, while in the laboratory culture of M. mirumnovem, it was nearer 30%. Their number in specimens of the culture varied from 1 to 14. Mosaicism on Bs was discovered in parts of these animals. We analyzed the distribution of Bs among the worms of the laboratory cultures during long-term cultivation, the transmission rates of Bs in the progeny obtained from crosses of worms with different numbers of Bs, and from self-fertilized isolated worms. The DNA content of the Bs in M. mirumnovem was analyzed with the chromosomal in situ suppression (CISS) hybridization of microdissected DNA probes derived from A chromosomes (As). Bs mainly consisted of repetitive DNA. The cytogenetic analysis also revealed the divergence and high variation in large metacentric chromosomes (LMs) containing numerous regions enriched for repeats. The possible mechanisms of the appearance and evolution of Bs and LMs in species of the Macrostomum genus were also discussed.

DNA methylation in transposable elements disrupts the connection between three-dimensional chromatin organization and gene expression upon rice genome duplication.

Polyploidy serves as a major force in plant evolution and domestication of cultivated crops. However, the relationship and underlying mechanism between three-dimensional (3D) chromatin organization and gene expression upon rice genome duplication is largely unknown. Here we compared the 3D chromatin structures between diploid (2C) and autotetraploid (4C) rice by high-throughput chromosome conformation capture analysis, and found that 4C rice presents weakened intra-chromosomal interactions compared to its 2C progenitor. Moreover, we found that changes of 3D chromatin organizations including chromatin compartments, topologically associating domain (TAD) and loops uncouple from gene expression. Moreover, DNA methylations in the regulatory sequences of genes in compartment A/B switched regions and TAD boundaries are not related to their expressions. Importantly, in contrast to that there was no significant difference of methylation levels in TEs in promoters of differentially expressed genes (DEGs) and non-DEGs between 2C and 4C rice, we found that the hypermethylated transposable elements across genes in compartment A/B switched regions and TAD boundaries suppress the expression of these genes. We propose that the rice genome doubling might modulate TE methylation which results in the disconnection between the alteration of 3D chromatin structure and gene expression.

Polyploidy in Industrial Crops: Applications and Perspectives in Plant Breeding

Polyploidisation is an important process in the evolution of many plant species. An additional set of chromosomes can be derived from intraspecific genome duplication (autopolyploidy) or hybridising divergent genomes and chromosome doubling (allopolyploidy). Special forms of polyploidy are autoallopolyploidy and segmental allopolyploidy. Polyploidy arises from two basic processes: spontaneously occurring disturbances of meiotic division and induced by antimitotic agents’ disruption of mitosis. The first involves the induction and fusion of unreduced gametes, resulting in the formation of triploids and tetraploids. The second process uses antimitotics that disrupt cellular microtubules and prevent chromosome’s sister chromatids motion during anaphase. Colchicine, oryzalin, and trifluralin are the most commonly used antimitotics for inducing polyploids in plants. The exposure time and concentration of the antimitotics and the species, cultivar, genotype, and tissue type affect the efficiency of genome duplication. Polyploids are distinguished from diploids by increased cell size and vegetative parts of plants and increased content of secondary metabolites. Genome duplication generates several changes at the epigenetic level resulting in altered gene expression. Polyploidisation is used in plant breeding to overcome the non-viability and infertility of interspecific hybrids, obtain seedless polyploid cultivars and increase resistance/tolerance to biotic and abiotic factors.