Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies.

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus).

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

High-quality Japanese flounder genome aids in identifying stress-related genes using gene coexpression network.

The Japanese flounder is one of the most economically important marine flatfish. However, due to the increased frequency of extreme weather events and high-density industrial farming, an increasing number of environmental stresses have become severe threats to the healthy development of the Japanese flounder culture industry. Herein, we produced a high-quality chromosome-scale Japanese flounder genome using PacBio Circular Consensus Sequencing technologies. The assembled Japanese flounder genome spanned 588.22 Mb with a contig N50 size of 24.35 Mb. In total, 105.89 Mb of repetitive sequences and 22,565 protein-coding genes were identified by genome annotation. In addition, 67 candidate genes responding to distinct stresses were identified by gene coexpression network analysis based on 16 published stress-related RNA-seq datasets encompassing 198 samples. A high-quality chromosome-scale Japanese flounder genome and candidate stress-related gene set will not only serve as key resources for genomics studies and further research on the underlying stress responsive molecular mechanisms in Japanese flounder but will also advance the progress of genetic improvement and comprehensive stress-resistant molecular breeding of Japanese flounder.

Identification of stress-related genes by co-expression network analysis based on the improved turbot genome.

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties.

Draft genomes of female and male turbot Scophthalmus maximus.

Turbot (Scophthalmus maximus) is a commercially important flatfish species in aquaculture. It has a drastic sexual dimorphism, with females growing faster than males. In the present study, we sequenced and de novo assembled female and male turbot genomes. The assembled female genome was 568 Mb (scaffold N50, 6.2 Mb, BUSCO 97.4%), and the male genome was 584 Mb (scaffold N50, 5.9 Mb, BUSCO 96.6%). Using two genetic maps, we anchored female scaffolds representing 535 Mb onto 22 chromosomes. Annotation of the female anchored genome identified 87.8 Mb transposon elements and 20,134 genes. We identified 17,936 gene families, of which 369 gene families were flatfish specific. Phylogenetic analysis showed that the turbot, Japanese flounder and Chinese tongue sole form a clade that diverged from other teleosts approximately 78 Mya. This report of female and male turbot draft genomes and annotated genes provides a new resource for identifying sex determination genes, elucidating the evolution of adaptive traits in flatfish and developing genetic techniques to increase the sustainability of turbot aquaculture.

Genome-wide identification, characterization, and expression analysis of the TRAF gene family in Chinese tongue sole (Cynoglossus semilaevis).

Tumor necrosis factor (TNF) receptor-associated factors (TRAFs) play crucial roles as signaling mediators for the TNF receptor (TNFR) superfamily and the interleukin-1 receptor/Toll-like receptor (IL-1R/TLR) superfamily. TRAFs collectively play important roles in multiple biological processes and organismal immunity. However, systematic identification of the TRAF gene family in teleost fish has not yet been reported, and there is little available information about its roles in innate immunity in Chinese tongue sole (Cynoglossus semilaevis), an aquaculture fish of high economic value. In the present study, we identified and characterized seven TRAF genes, namely, CsTRAF2a, CsTRAF2b, CsTRAF3, CsTRAF4, CsTRAF5, CsTRAF6 and CsTRAF7, in Chinese tongue sole, and the complete ORFs of the CsTRAFs were cloned. Sequence analysis revealed various genomic structures of the CsTRAFs and showed that they contain typical conserved domains compared with mammalian TRAFs. Phylogenetic analysis indicated the evolutionary relationships of TRAF family members in teleost fish and revealed an absence of TRAF1 in most species and TRAF5 in some species of teleosts. Analysis of the gene structures and motifs showed the diversity and distribution of exon-intron structures and conserved motifs in Chinese tongue sole and several other teleost species. Real-time quantitative PCR was used to investigate the expression patterns of CsTRAF genes in tissues of healthy fish and in the gills, livers and spleens of fish after bacterial infection with Vibrio harveyi. The results indicate that only CsTRAF2a is relatively highly expressed in the brain and that the other CsTRAFs are highly expressed in immune-related tissues and may participate in the immune response after infection with pathogenic bacteria. Functional analysis of CsTRAF3, CsTRAF4 and CsTRAF6 revealed that only CsTRAF6 could strongly activate the NF-кB pathway after overexpression of CsTRAF3, CsTRAF4 and CsTRAF6 in HEK-293T cells. This systematic analysis provided valuable information about the diverse roles of TRAFs in the innate immune response to pathogenic bacterial infection in teleost fish and will contribute to the functional characterization of CsTRAF genes in further research.

Evolutionary dynamics of lincRNA transcription in nine citrus species.

Long intergenic non-coding RNAs (lincRNAs) play important roles in various biological processes in plants. However, little information is known about the evolutionary characteristics of lincRNAs among closely related plant species. Here, we present a large-scale comparative study of lincRNA transcription patterns in nine citrus species. By strand-specific RNA-sequencing, we identified 18 075 lincRNAs (14 575 lincRNA loci) from 34 tissue samples. The results indicated that the evolution of lincRNA transcription is more rapid than that of mRNAs. In total, 82.8-97.6% of sweet orange (Citrus sinensis) lincRNA genes were shown to have homologous sequences in other citrus genomes. However, only 15.5-28.8% of these genes had transcribed homologous lincRNAs in these citrus species, presenting a strong contrast to the high conservation of mRNA transcription (81.6-84.7%). Moreover, primitive and modern citrus lincRNAs were preferentially expressed in reproductive and vegetative organs, respectively. Evolutionarily conserved lincRNAs showed higher expression levels and lower tissue specificity than species-specific lincRNAs. Notably, we observed a similar tissue expression pattern of homologous lincRNAs in sweet orange and pummelo (Citrus grandis), suggesting that these lincRNAs may be functionally conserved and selectively maintained. We also identified and validated a lincRNA with the highest expression in fruit that acts as an endogenous target mimic (eTM) of csi-miR166c, and two lincRNAs that act as a precursor and target of csi-miR166c, respectively. These lincRNAs together with csi-miR166c could form an eTM166-miR166c-targeted lincRNA regulatory network that possibly affects citrus fruit development.

Contributions of Zea mays subspecies mexicana haplotypes to modern maize.

Maize was domesticated from lowland teosinte (Zea mays ssp. parviglumis), but the contribution of highland teosinte (Zea mays ssp. mexicana, hereafter mexicana) to modern maize is not clear. Here, two genomes for Mo17 (a modern maize inbred) and mexicana are assembled using a meta-assembly strategy after sequencing of 10 lines derived from a maize-teosinte cross. Comparative analyses reveal a high level of diversity between Mo17, B73, and mexicana, including three Mb-size structural rearrangements. The maize spontaneous mutation rate is estimated to be 2.17 × 10-8 ~3.87 × 10-8 per site per generation with a nonrandom distribution across the genome. A higher deleterious mutation rate is observed in the pericentromeric regions, and might be caused by differences in recombination frequency. Over 10% of the maize genome shows evidence of introgression from the mexicana genome, suggesting that mexicana contributed to maize adaptation and improvement. Our data offer a rich resource for constructing the pan-genome of Zea mays and genetic improvement of modern maize varieties.

Functional analysis of long intergenic non-coding RNAs in phosphate-starved rice using competing endogenous RNA network.

Long intergenic non-coding RNAs (lincRNAs) may play widespread roles in gene regulation and other biological processes, however, a systematic examination of the functions of lincRNAs in the biological responses of rice to phosphate (Pi) starvation has not been performed. Here, we used a computational method to predict the functions of lincRNAs in Pi-starved rice. Overall, 3,170 lincRNA loci were identified using RNA sequencing data from the roots and shoots of control and Pi-starved rice. A competing endogenous RNA (ceRNA) network was constructed for each tissue by considering the competing relationships between lincRNAs and genes, and the correlations between the expression levels of RNAs in ceRNA pairs. Enrichment analyses showed that most of the communities in the networks were related to the biological processes of Pi starvation. The lincRNAs in the two tissues were individually functionally annotated based on the ceRNA networks, and the differentially expressed lincRNAs were biologically meaningful. For example, XLOC_026030 was upregulated from 3 days after Pi starvation, and its functional annotation was 'cellular response to Pi starvation'. In conclusion, we systematically annotated lincRNAs in rice and identified those involved in the biological response to Pi starvation.

The draft genome of sweet orange (Citrus sinensis).

Oranges are an important nutritional source for human health and have immense economic value. Here we present a comprehensive analysis of the draft genome of sweet orange (Citrus sinensis). The assembled sequence covers 87.3% of the estimated orange genome, which is relatively compact, as 20% is composed of repetitive elements. We predicted 29,445 protein-coding genes, half of which are in the heterozygous state. With additional sequencing of two more citrus species and comparative analyses of seven citrus genomes, we present evidence to suggest that sweet orange originated from a backcross hybrid between pummelo and mandarin. Focused analysis on genes involved in vitamin C metabolism showed that GalUR, encoding the rate-limiting enzyme of the galacturonate pathway, is significantly upregulated in orange fruit, and the recent expansion of this gene family may provide a genomic basis. This draft genome represents a valuable resource for understanding and improving many important citrus traits in the future.