Whole genome-wide chromosome fusion and new gene birth in the Monopterus albus genome

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

Swamp eel aldehyde reductase is involved in response to nitrosative stress via regulating NO/GSH levels

Aldehyde reductase (ALR) plays key roles in the detoxification of toxic aldehyde. In this study, the authors cloned the swamp eel ALR gene using rapid amplification of cDNA ends-PCR (RACE-PCR). The recombinant protein (rALR) was expressed in Escherichia coli and purified using a Ni2+ -NTA chelating column. The rALR protein exhibited efficient reductive activity towards several aldehydes, ketones and S-nitrosoglutathione (GSNO). A spot assay suggested that the recombinant E. coli strain expressing rALR showed better resistance to formaldehyde, sodium nitrite and GSNO stress, suggesting that swamp eel ALR is crucial for redox homeostasis in vivo. Consequently, the authors investigated the effect of rALR on the oxidative parameters of the liver in swamp eels challenged with Aeromonas hydrophila. The hepatic glutathione (GSH) content significantly increased, and the hepatic NO content and levels of reactive oxygen species and reactive nitrogen species significantly decreased when rALR was administered. In addition, the mRNA expression of hepatic Alr, HO1 and Nrf2 was significantly upregulated, whereas the expression levels of NF-κB, IL-1β and NOS1 were significantly downregulated in the rALR-administered group. Collectively, these results suggest that ALR is involved in the response to nitrosative stress by regulating GSH/NO levels in the swamp eel.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

Chromosome-Scale, Haplotype-Resolved Genome Assembly of Non-Sex-Reversal Females of Swamp Eel Using High-Fidelity Long Reads and Hi-C Data

The Asian swamp eel (Monopterus albus) is an excellent model species for studying sex change and chromosome evolution. M. albus is also widely reared in East Asia and South-East Asia because of its great nutritional value. The low fecundity of this species (about 300 eggs per fish) greatly hinders fries production and breeding programs. Interestingly, about 3–5% of the eels could remain as females for 3 years and lay more than 3,000 eggs per fish, which are referred to as non-sex-reversal (NSR) females. Here, we presented a new chromosome-level genome assembly of such NSR females using Illumina, HiFi, and Hi-C sequencing technologies. The new assembly (Mal.V2_NSR) is 838.39 Mb in length, and the N50 of the contigs is 49.8 Mb. Compared with the previous assembly obtained using the continuous long-read sequencing technology (Mal.V1_CLR), we found a remarkable increase of continuity in the new assembly Mal.V2_NSR with a 20-times longer contig N50. Chromosomes 2 and 12 were assembled into a single contig, respectively. Meanwhile, two highly contiguous haplotype assemblies were also obtained, with contig N50 being 14.54 and 12.13 Mb, respectively. BUSCO and Merqury analyses indicate completeness and high accuracy of these three assemblies. A comparative genomic analysis revealed substantial structural variations (SVs) between Mal.V2_NSR and Mal.V1_CLR and two phased haplotype assemblies, as well as whole chromosome fusion events when compared with the zig-zag eel. Additionally, our newly obtained assembly provides a genomic view of sex-related genes and a complete landscape of the MHC genes. Therefore, these high-quality genome assemblies would provide great help for future breeding works of the swamp eel, and it is a valuable new reference for genetic and genomic studies of this species.

Decoding genome recombination and sex reversal

Over the past 440 years since the discovery of the medicinal value of swamp eels, much progress has been made in the study of their biology. The fish is emerging as an important model animal in sexual development, in addition to economic and pharmaceutical implications. Tracing genomic history that shapes speciation of the fish has led to discovery of the whole genome-wide chromosome fission/fusion events. Natural intersex differentiation is a compelling feature for sexual development research. Notably, identification of progenitors of germline stem cells that have bipotential to differentiate into either male or female germline stem cells provides new insight into sex reversal. Here, we review these advances that have propelled the field forwards and present unsolved issues that will guide future investigations to finally elucidate vertebrate sexual development using the new model.

Susceptibility of a cell line derived from the kidney of Chinese rice-field eel, Monopterus albus to the infection of rhabdovirus, CrERV

Chinese rice-field eels rhabdovirus (CrERV), belonging to the genus Perhabdovirus in the family Rhabdoviridae, is the causative agent of the haemorrhagic disease of Chinese rice-field eels, Monopterus albus. The present study aims to establish a cell line derived from the kidney of Chinese rice-field eel (CrEK) for the further study of the pathogenic virus. CrEK cells were epithelioid-like and grew well in M199 medium supplemented with 10% foetal bovine serum at 28°C, and the cell line has been subcultured for more than 80 times. Karyotyping analysis of CrEK cells at 25th passage indicated a modal chromosome number of 24. Significant cytopathic effect (CPE) was observed in CrEK cells after infection with CrERV, and the virus titre reached 10^{7.8 ± 0.45} TCID₅₀ /mL. The transmission electron microscopy revealed that there were a large number of virus particles in the cytoplasm of cells. The virus infection in cells was also assayed by using indirect immunofluorescence assay (IFA), fluorescence in situ hybridization (FISH), reverse transcription PCR (RT-PCR) and quantitative real-time reverse transcription-PCR (qRT-PCR). In experimental infection, CrERV cultured by cells could cause over 90% mortality in fish. CrEK represents the first kidney cell line originated from Chinese rice-field eels and be a potential material for investigating the mechanism of virus infection in this fish and the control methods for the disease.

Chromosome restructuring and number change during the evolution of Morus notabilis and Morus alba

Mulberry (Morus spp.) is an economically important plant as the main food plant used for rearing domesticated silkworm and it has multiple uses in traditional Chinese medicine. Two basic chromosome numbers (Morus notabilis, n = 7, and Morus alba, n = 14) have been reported in the genus Morus, but the evolutionary history and relationship between them remain unclear. In the present study, a 335-Mb high-quality chromosome-scale genome was assembled for the wild mulberry species M. notabilis. Comparative genomic analyses indicated high chromosomal synteny between the 14 chromosomes of cultivated M. alba and the six chromosomes of wild M. notabilis. These results were successfully verified by fluorescence in situ hybridization. Chromosomal fission/fusion events played crucial roles in the chromosome restructuring process between M. notabilis and M. alba. The activity of the centromere was another key factor that ensured the stable inheritance of chromosomes. Our results also revealed that long terminal repeat retrotransposons were a major driver of the genome divergence and evolution of the mulberry genomes after they diverged from each other. This study provides important insights and a solid foundation for studying the evolution of mulberry, allowing the accelerated genetic improvement of cultivated mulberry species.

Histological alterations, oxidative stress, and inflammatory response in the liver of swamp eel (Monopterus albus) acutely exposed to copper

Copper (Cu) is widely used as an essential trace element in diets as well as a therapeutic chemical. However, excessive Cu has deleterious effects on organisms, including teleosts. Although numerous toxic effects of Cu have been reported, the effects of Cu exposure on the swamp eel (Monopterus albus) as well as the underlying mechanisms have not yet been elucidated. In this study, swamp eels were acutely exposed to 100, 200, and 400 μg/L of Cu for 96 h to evaluate liver histopathology, oxidative stress, and inflammation. Dissolution of hepatocyte membrane, vacuolar degeneration, and inflammatory cell infiltration were detected in the livers of the Cu-treated swamp eels, especially in the 400 μg Cu/L group. Cu-induced hepatic dysfunction was further verified by the elevated activities of glutamate oxaloacetate transaminase (GOT) and glutamate pyruvate transaminase (GPT) and transcript levels of GOT and GPT genes. In addition, Cu exposure decreased the activities of total superoxide dismutase T-SOD and catalase (CAT) and the contents of glutathione (GSH) and total antioxidant capacity (T-AOC) and increased the levels of malondialdehyde (MDA). Cu exposure also significantly decreased the transcript levels of glutathione synthetase (GSS) and increased the transcript levels of SOD1, SOD2, CAT, and heme oxygenase-1 (HO-1) genes. Furthermore, pro-inflammatory genes such as interleukin (IL)-1β, tumor necrosis factor-α (TNF-α), and IL-8 were significantly upregulated. These results indicate that Cu induces oxidative stress and inflammatory response and causes pathological changes in the liver of the swamp eel.

Pattern of New Gene Origination in a Special Fish Lineage, the Flatfishes

Origination of new genes are of inherent interest of evolutionary geneticists for decades, but few studies have addressed the general pattern in a fish lineage. Using our recent released whole genome data of flatfishes, which evolved one of the most specialized body plans in vertebrates, we identified 1541 (6.9% of the starry flounder genes) flatfish-lineage-specific genes. The origination pattern of these flatfish new genes is largely similar to those observed in other vertebrates, as shown by the proportion of DNA-mediated duplication (1317; 85.5%), RNA-mediated duplication (retrogenes; 96; 6.2%), and de novo-origination (128; 8.3%). The emergence rate of species-specific genes is 32.1 per Mya and the whole average level rate for the flatfish-lineage-specific genes is 20.9 per Mya. A large proportion (31.4%) of these new genes have been subjected to selection, in contrast to the 4.0% in primates, while the old genes remain quite similar (66.4% vs. 65.0%). In addition, most of these new genes (70.8%) are found to be expressed, indicating their functionality. This study not only presents one example of systematic new gene identification in a teleost taxon based on comprehensive phylogenomic data, but also shows that new genes may play roles in body planning.

Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides)

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1 594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that, similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex-determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Cellular fate of intersex differentiation

Infertile ovotestis (mixture of ovary and testis) often occurs in intersex individuals under certain pathological and physiological conditions. However, how ovotestis is formed remains largely unknown. Here, we report the first comprehensive single-cell developmental atlas of the model ovotestis. We provide an overview of cell identities and a roadmap of germline, niche, and stem cell development in ovotestis by cell lineage reconstruction and a uniform manifold approximation and projection. We identify common progenitors of germline stem cells with two states, which reveal their bipotential nature to differentiate into both spermatogonial stem cells and female germline stem cells. Moreover, we found that ovotestis infertility was caused by degradation of female germline cells via liquid-liquid phase separation of the proteasomes in the nucleus, and impaired histone-to-protamine replacement in spermatid differentiation. Notably, signaling pathways in gonadal niche cells and their interaction with germlines synergistically determined distinct cell fate of both male and female germlines. Overall, we reveal a cellular fate map of germline and niche cell development that shapes cell differentiation direction of ovotestis, and provide novel insights into ovotestis development.

DNA methylation modification is associated with gonadal differentiation in Monopterus albus

Background

Both testis and ovary can be produced sequentially in an individual with the same genome when sex reversal occurs in the teleost Monopterus albus, and epigenetic modification is supposed to be involved in gonadal differentiation. However, DNA methylation regulation mechanism underlying the gonadal differentiation remains unclear.

Results

Here, we used liquid chromatography-electrospray ionization tandem mass spectrometry (LC–ESI–MS/MS) to simultaneously determine endogenous levels of both 5-methyl-2′-deoxycytidine (m⁵dC) and 5-hydroxymethyl-2′-deoxycytidine (hm⁵dC) during gonadal differentiation. Overall DNA methylation level was upregulated from ovary to testis via ovotestis. As a de novo methylase, dnmt3aa expression was also upregulated in the process. Notably, we determined transcription factor Foxa1 for dnmt3aa gene expression. Site-specific mutations and chromatin immunoprecipitation showed that Foxa1 can bind to and activate the dnmt3aa promoter. Furthermore, DNA methylation levels of key genes foxl2 (forkhead box L2) and cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a) in regulation of female hormone synthesis were consistently upregulated during gonadal differentiation.

Conclusions

These data suggested that dynamic change of DNA methylation modification is associated with gonadal differentiation.

Regulation of arachidonic acid oxidation and metabolism by lipid electrophiles

Arachidonic acid (AA) is a precursor of enzymatic and non-enzymatic oxidized products such as prostaglandins, thromboxanes, leukotrienes, lipoxins, and isoprostanes. These products may exert signaling or damaging roles during physiological and pathological conditions, some of them being markers of oxidative stress linked to inflammation. Recent data support the concept that cyclooxygenases (COX), lipoxygenases (LOX), and cytochrome P450 (CYP450) followed by cytosolic and microsomal dehydrogenases can convert AA to lipid-derived electrophiles (LDE). Lipid-derived electrophiles are fatty acid derivatives bearing an electron-withdrawing group that can react with nucleophiles at proteins, DNA, and small antioxidant molecules exerting potent signaling properties. This review aims to describe the formation, sources, and electrophilic anti-inflammatory actions of key mammalian LDE.