Chromosome-level genome assembly of the Chinese three-keeled pond turtle (Mauremys reevesii) provides insights into freshwater adaptation

Genome sequencing and transcriptome analysis provide an insight into the hypoxia resistance of Channa asiatica

Channa asiatica is an economically valuable fish species and excellent model for studying hypoxic tolerance. However, the underlying genetic and molecular mechanisms are poorly understood. In this study, we assembled a high-quality C. asiatica genome (23 chromosomes, totaling 722 Mb) using a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing. Repetitive elements accounted for 28.39%of the C. asiatica genome, and 23,949 protein-coding genes were predicted, with 96.63 % of these functionally annotated. Moreover, a comparative genomic analysis of 12 fish genomes showed that gene families associated with oxygen binding and transport were expanded in C. asiatica. In addition, transcriptome analysis revealed that multiple oxidative stress pathways were activated when C. asiatica was exposed to air. In conclusion, this study provided high-quality genome assembly and transcriptome data, both serving as critical resources for researching the genetic basis of hypoxic tolerance in C. asiatica.

De novo genome assemblies of two cryptodiran turtles with ZZ/ZW and XX/XY sex chromosomes provide insights into patterns of genome reshuffling and uncover novel 3D genome folding in amniotes

Understanding the evolution of chromatin conformation among species is fundamental to elucidate the architecture and plasticity of genomes. Nonrandom interactions of linearly distant loci regulate gene function in species-specific patterns, affecting genome function, evolution, and, ultimately, speciation. Yet, data from nonmodel organisms are scarce. To capture the macroevolutionary diversity of vertebrate chromatin conformation, here we generate de novo genome assemblies for two cryptodiran (hidden-neck) turtles via Illumina sequencing, chromosome conformation capture, and RNA-seq: Apalone spinifera (ZZ/ZW, 2n = 66) and Staurotypus triporcatus (XX/XY, 2n = 54). We detected differences in the three-dimensional (3D) chromatin structure in turtles compared to other amniotes beyond the fusion/fission events detected in the linear genomes. Namely, whole-genome comparisons revealed distinct trends of chromosome rearrangements in turtles: (1) a low rate of genome reshuffling in Apalone (Trionychidae) whose karyotype is highly conserved when compared to chicken (likely ancestral for turtles), and (2) a moderate rate of fusions/fissions in Staurotypus (Kinosternidae) and Trachemys scripta (Emydidae). Furthermore, we identified a chromosome folding pattern that enables "centromere-telomere interactions" previously undetected in turtles. The combined turtle pattern of "centromere-telomere interactions" (discovered here) plus "centromere clustering" (previously reported in sauropsids) is novel for amniotes and it counters previous hypotheses about amniote 3D chromatin structure. We hypothesize that the divergent pattern found in turtles originated from an amniote ancestral state defined by a nuclear configuration with extensive associations among microchromosomes that were preserved upon the reshuffling of the linear genome.

The Evolution of Transglutaminases Underlies the Origin and Loss of Cornified Skin Appendages in Vertebrates

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

Zona pellucida family genes in Chinese pond turtle: identification, expression profiles, and role in the spermatozoa acrosome reaction†

The zona pellucida (ZP) is an extracellular matrix that surrounds all vertebrate eggs, and it is involved in fertilization and species-specific recognition. Numerous in-depth studies of the ZP proteins of mammals, birds, amphibians, and fishes have been conducted, but systematic investigation of the ZP family genes and their role during fertilization in reptiles has not been reported to date. In this study, we identified six turtle ZP (Tu-ZP) gene subfamilies (Tu-ZP1, Tu-ZP2, Tu-ZP3, Tu-ZP4, Tu-ZPD, and Tu-ZPAX) based on whole genome sequence data from Mauremys reevesii. We found that Tu-ZP4 had large segmental duplication and was distributed on three chromosomes, and we also detected gene duplication in the other Tu-ZP genes. To evaluate the role of Tu-ZP proteins in sperm-egg binding, we assessed the expression pattern of these Tu-ZP proteins and their ability to induce the spermatozoa acrosome reaction in M. reevesii. In conclusion, this is the first report of the existence of gene duplication of Tu-ZP genes and that Tu-ZP2, Tu-ZP3, and Tu-ZPD can induce acrosome exocytosis of spermatogenesis in the reptile.

Probiotics and Immunostimulant modulate intestinal flora diversity in Reeves pond tortoise (Mauremys reevesii) and effects of Clostridium butyricum on its spleen transcriptome

In this study, we investigated the effects of Clostridium butyricum (group A), Bacillus subtilis (group B), and the immune enhancer algal β-1,3 glucan (group C) on the intestinal flora of Mauremys reevesii and the effects of C. butyricum on the transcriptome of M. reevesii splenic immune tissues. M. reevesii were assigened to four groups, each containing three replicates from 18 samples. Juvenile turtles with an initial weight of 106.35 ± 0.03 g were fed a basic diet containing no probiotics (control group D), or a basic diet containing the A, B, or C supplement. After the turtles had been fed for 60, 90, and 120 d of the experimental period, high-throughput sequencing of the 16S rRNA gene revealed no significant difference in alpha diversity among the four groups at 60 days of feeding (P > 0.05), and at 90 days, the alpha diversity in group A was significantly different (P < 0.05), with an increase of 26.62% in the Shannon index and a decrease of 83.33% in the Simpson index; at 120 d, the alpha diversity (Shannon index)showed a decreasing trend in order for groups A, B, and C, At the phylum level, the abundance of Bacteroidetes, Proteobacteria, and Fusobacteria in group A increased significantly with increasing feeding time (P < 0.05)，At the genus level, the abundance of Ruminococcaceae and Anaerotruncus in group A increased significantly compared with that in the other three groups (P < 0.05). Transcriptome analysis showed that 384 genes were differentially expressed in the spleen of M. reevesii, 195 genes were upregulated and 189 genes were downregulated, and C. butyricum TF201120 regulated the hematopoietic cell lineage signaling pathway in the spleen of M. reevesii (P < 0.05). The regulation of several identified immune-related genes was confirmed by qPCR, These results showed that C. butyricum, B. subtilis and the immune enhancer algal β-1,3 glucan can improve the intestinal flora of M. reevesii, with C. butyricum TF20 being the most effective and significantly enhancing the immunity of M. reevesii.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation

BACKGROUND

The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation.

RESULTS

The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude.

CONCLUSIONS

Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

Chromosome-Level Genome Assembly of a Fragrant Japonica Rice Cultivar 'Changxianggeng 1813' Provides Insights into Genomic Variations between Fragrant and Non-Fragrant Japonica Rice

East Asia has an abundant resource of fragrant japonica rice that is gaining increasing interest among both consumers and producers. However, genomic resources and in particular complete genome sequences currently available for the breeding of fragrant japonica rice are still scarce. Here, integrating Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods, we presented a high-quality chromosome-level genome assembly (~378.78 Mb) for a new fragrant japonica cultivar ‘Changxianggeng 1813’, with 31,671 predicated protein-coding genes. Based on the annotated genome sequence, we demonstrated that it was the badh2-E2 type of deletion (a 7-bp deletion in the second exon) that caused fragrance in ‘Changxianggeng 1813’. Comparative genomic analyses revealed that multiple gene families involved in the abiotic stress response were expanded in the ‘Changxianggeng 1813’ genome, which further supported the previous finding that no generalized loss of abiotic stress tolerance associated with the fragrance phenotype. Although the ‘Changxianggeng 1813’ genome showed high genomic synteny with the genome of the non-fragrant japonica rice cultivar Nipponbare, a total of 289,970 single nucleotide polymorphisms (SNPs), 96,093 small insertion-deletion polymorphisms (InDels), and 8690 large structure variants (SVs, >1000 bp) were identified between them. Together, these genomic resources will be valuable for elucidating the mechanisms underlying economically important traits and have wide-ranging implications for genomics-assisted breeding in fragrant japonica rice.

Chronology of embryonic and gonadal development in the Reeves' turtle, Mauremys reevesii

Temperature-dependent sex determination (TSD) is a mechanism in which environmental temperature, rather than innate zygotic genotype, determines the fate of sexual differentiation during embryonic development. Reeves' turtle (also known as the Chinese three-keeled pond turtle, Mauremys reevesii) exhibits TSD and is the only species whose genome has been determined in Geoemydidae to date. Thus, M. reevesii occupy phylogenetically important position for the study of TSD and can be compared to other TSD species to elucidate the underlying molecular mechanism of this process. Nevertheless, neither embryogenesis nor gonadogenesis has been described in this species. Therefore, herein, we investigated the chronology of normal embryonic development and gonadal structures in M. reevesii under both female- and male-producing incubation temperatures (FPT 31 °C or MPT 26 °C, respectively). External morphology remains indistinct between the two temperature regimes throughout the studied embryonic stages. However, the gonadal ridges present on the mesonephros at stage 16 develop and sexually differentiate at FPT and MPT. Ovarian and testicular structures begin to develop at stages 18-19 at FPT and stages 20-21 at MPT, respectively, and thus, the sexual differentiation of gonadal structures began earlier in the embryos at FPT than at MPT. Our results suggest that temperature sensitive period, at which the gonadal structures remain sexually undifferentiated, spans from stage 16 (or earlier) to stages 18-19 at FPT and to stages 20-21 at MPT. Understanding the temperature-dependent differentiation in gonadal structures during embryonic development is a prerequisite for investigating molecular basis underlying TSD. Thus, the result of the present study will facilitate further developmental studies on TSD in M. reevesii.