TACO produces robust multisample transcriptome assemblies from RNA-seq

Near-complete assembly and comprehensive annotation of the wheat Chinese Spring genome

A complete reference genome assembly is crucial for biological research and genetic improvement. Owing to its large size and highly repetitive nature, there are numerous gaps in the globally used wheat Chinese Spring (CS) genome assembly. In this study, we generated a 14.46 Gb near-complete assembly of the CS genome, with a contig N50 of over 266 Mb and an overall base accuracy of 99.9963%. Among the 290 gaps that remained (26, 257, and 7 gaps from the A, B, and D subgenomes, respectively), 278 were extremely high-copy tandem repeats, whereas the remaining 12 were transposable-element-associated gaps. Four chromosome assemblies were completely gap-free, including chr1D, chr3D, chr4D, and chr5D. Extensive annotation of the near-complete genome revealed 151 405 high-confidence genes, of which 59 180 were newly annotated, including 7602 newly assembled genes. Except for the centromere of chr1B, which has a gap associated with superlong GAA repeat arrays, the centromeric sequences of all of the remaining 20 chromosomes were completely assembled. Our near-complete assembly revealed that the extent of tandem repeats, such as simple-sequence repeats, was highly uneven among different subgenomes. Similarly, the repeat compositions of the centromeres also varied among the three subgenomes. With the genome sequences of all six types of seed storage proteins (SSPs) fully assembled, the expression of ω-gliadin was found to be contributed entirely by the B subgenome, whereas the expression of the other five types of SSPs was most abundant from the D subgenome. The near-complete CS genome will serve as a valuable resource for genomic and functional genomic research and breeding of wheat as well as its related species.

Building better genome annotations across the tree of life

Recent technological advances in long-read DNA sequencing accompanied by reduction in costs have made the production of genome assemblies financially achievable and computationally feasible, such that genome assembly no longer represents the major hurdle to evolutionary analysis for most nonmodel organisms. Now, the more difficult challenge is to properly annotate a draft genome assembly once it has been constructed. The primary challenge to annotations is how to select from the myriad gene prediction tools that are currently available, determine what kinds of data are necessary to generate high-quality annotations, and evaluate the quality of the annotation. To determine which methods perform the best and to determine whether the inclusion of RNA-seq data is necessary to obtain a high-quality annotation, we generated annotations with 12 different methods for 21 different species spanning vertebrates, plants, and insects. We found that the annotation transfer method TOGA, BRAKER3, and the RNA-seq assembler StringTie were consistently top performers across a variety of metrics including BUSCO recovery, CDS length, and false-positive rate, with the exception that TOGA performed less well in some monocots with respect to BUSCO recovery. The choice of which of the top-performing methods will depend upon the feasibility of whole-genome alignment, availability of RNA-seq data, importance of capturing noncoding parts of the transcriptome, and, when whole-genome alignment is not feasible, the relative performance in BUSCO recovery between BRAKER3 and StringTie. When whole-genome alignment is not feasible, inclusion of RNA-seq data will lead to substantial improvements to genome annotations.

Signatures of Selection for Resistance/Tolerance to Perkinsus olseni in Grooved Carpet Shell Clam (Ruditapes decussatus) Using a Population Genomics Approach

The grooved carpet shell clam (Ruditapes decussatus) is a bivalve of high commercial value distributed throughout the European coast. Its production has suffered a decline caused by different factors, especially by the parasite Perkinsus olsenii. Improving production of R. decussatus requires genomic resources to ascertain the genetic factors underlying resistance/tolerance to P. olsenii. In this study, the first reference genome of R. decussatus was assembled through long- and short-read sequencing (1677 contigs; 1.386 Mb) and further scaffolded at chromosome level with Hi-C (19 superscaffolds; 95.4% of assembly). Repetitive elements were identified (32%) and masked for annotation of 38,276 coding- and 13,056 non-coding genes. This genome was used as a reference to develop a 2bRAD-Seq 13,438 SNP panel for a genomic screening on six shellfish beds distributed across the Atlantic Ocean and Mediterranean Sea. Beds were selected by perkinsosis prevalence and the infection level was individually evaluated in all the samples. Genetic diversity was significantly higher in the Mediterranean than in the Atlantic region. The main genetic breakage was detected between those regions (FST = 0.224), being the Mediterranean more heterogeneous than the Atlantic. Several loci under divergent selection (394 outliers; 261 genomic windows) were detected across shellfish beds. Samples were also inspected to detect signals of selection for resistance/tolerance to P. olsenii by using infection-level and population-genomics approaches, and 90 common divergent outliers for resistance/tolerance to perkinsosis were identified and used for gene mining. Candidate genes and markers identified provide invaluable information for controlling perkinsosis and for improving production of the grooved carpet shell clam.

A chromosome-level genome assembly of the polyembryonic endoparasitoid Macrocentrus cingulum

Polyembryony is a unique reproductive pattern, where multiple genetically identical offspring develop from a single egg. Macrocentrus cingulum, a polyembryonic endoparasitoid, is one of the dominant parasitoids of the globally agricultural pest, Ostrinia furnacalis, and it can serve as a valuable model for investigating the mechanisms of polyembryony. However, the previously published genome of M. cingulum remains highly fragmented, limiting comprehensive studies of its biological characteristics. Here, we present a chromosome-level genome assembly of M. cingulum using PacBio HiFi and Hi-C sequencing. This genome assembly is approximately 151.93 Mb with a contig N50 of 8.42 Mb and a scaffold N50 of 16.93 Mb, organized into 9 chromosomes. The repeat sequences constitute 25.52% of the genome assembly. A total of 14,471 protein-coding genes with 98.2% BUSCO completeness were predicted, of which 12,500 genes have been annotated in public biological databases. In conclusion, this reported genome should be a valuable genomic resource for exploring macroevolutionary mechanisms underlying polyembryony, and the molecular mechanisms of polyembryony in M. cingulum.

Chromosome-level genome assembly and annotation of Anthurium amnicola

Anthurium amnicola is in the monocot family Araceae, subfamily Pothoideae and is a contributing species in Hawaii floriculture industry hybrids. To support future genetic improvements to this commodity, we sequenced and assembled the A. amnicola genome to chromosome-scale using PacBio HiFi and short-read Hi-C sequencing. A total of 98.51% of the 4.79 Gb genome is anchored into 15 chromosomes, with 99.2% gene completeness and a high LTR assembly index (LAI) score of 21.73, indicative of a complete, high-quality assembly. Annotation reveals the presence of 20,380 protein-coding genes, with 78.52% of the genome composed of repetitive sequences, predominantly long terminal repeat retrotransposons (LTR-RT). Phylogenetic analysis identified evolutionary relationships between A. amnicola and representative species in the Araceae and other plant families. This study provides the first reference genome sequence for the neotropical genus Anthurium and insights into Araceae evolution, benefiting the floriculture industry and evolutionary studies.

Chromosome-level genome assembly and annotation of the White-spotted spinefoot Siganus canaliculatus

The White-spotted spinefoot S. canaliculatus, is an economically important marine fish in South China and featured by possessing poisonous glands in its fin spines. However, the unavailability of the S. canaliculatus genome has been a serious obstacle to genetic breeding as well as basic researches such as uncovering genomic basis underlying its toxigenic glands. Here, we presented a chromosome-level genome assembly coupled with good annotation of S. canaliculatus using multiple omics technologies. The assembled genome size was 547.39 Mb, with a contig N50 and scaffold N50 length of 21.41 Mb and 21.79 Mb, respectively. Approximately 95.32% (521.76 Mb) of assembled sequences were placed into 24 pseudochromosomes with the support of Hi-C contact map. Furthermore, around 16.37% of the genome was composed of repetitive elements. The quality of the assembly assessed using BUSCO showed that 98.6% of BUSCO genes were identified as complete. 25,323 protein-coding genes were predicted after integration of three kinds of evidence, of which 96.96% were functionally annotated in at least one of nine protein databases. In sum, the chromosome-level genome assembly and annotation provide fundamental resources for genetic breeding and molecular mechanism related studies of S. canaliculatus.

Chromosome-Level Genome Assembly and Annotation of Corallium rubrum: A Mediterranean Coral Threatened by Overharvesting and Climate Change

Reference genomes are key resources in biodiversity conservation. Yet, sequencing efforts are not evenly distributed across the tree of life raising concerns over our ability to enlighten conservation with genomic data. Good-quality reference genomes remain scarce in octocorals while these species are highly relevant targets for conservation. Here, we present the first annotated reference genome in the red coral, Corallium rubrum (Linnaeus, 1758), a habitat-forming octocoral from the Mediterranean and neighboring Atlantic, impacted by overharvesting and anthropogenic warming-induced mass mortality events. Combining long reads from Oxford Nanopore Technologies (ONT), Illumina paired-end reads for improving the base accuracy of the ONT-based genome assembly, and Arima Hi-C contact data to place the sequences into chromosomes, we assembled a genome of 532 Mb (20 chromosomes, 309 scaffolds) with contig and scaffold N50 of 1.6 and 18.5 Mb, respectively. Fifty percent of the sequence (L50) was contained in seven superscaffolds. The consensus quality value of the final assembly was 42, and the single and duplicated gene completeness reported by BUSCO was 86.4% and 1%, respectively (metazoa_odb10 database). We annotated 26,348 protein-coding genes and 34,548 noncoding transcripts. This annotated chromosome-level genome assembly, one of the first in octocorals and the first in Scleralcyonacea order, is currently used in a project based on whole-genome resequencing dedicated to the conservation and management of C. rubrum.

Atlantic mackerel population structure does not support genetically distinct spawning components

Background

The Atlantic mackerel, Scomber scombrus (Linnaeus, 1758) is a commercially valuable migratory pelagic fish inhabiting the northern Atlantic Ocean and the Mediterranean Sea. Given its highly migratory behaviour for feeding and spawning, several studies have been conducted to assess differentiation among spawning components to better define management units, as well as to investigate possible adaptations to comprehend and predict recent range expansion northwards.

Methods

Here, the genome of S. scombrus was sequenced and annotated, as an increasing number of population genetic studies have proven the relevance of reference genomes to investigate genomic markers/regions potentially linked to differences at finer scale. Such reference genome was used to map Restriction-site-associated sequencing (RAD-seq) reads for SNP discovery and genotyping in more than 500 samples distributed along the species range. The resulting genotyping tables have been used to perform connectivity and adaptation analyses.

Results

The assembly of the reference genome for S. scombrus resulted in a genome of 741 Mb. Our population genetic results show that the Atlantic mackerel consist of three previously known genetically isolated units (Northwest Atlantic, Northeast Atlantic, Mediterranean), and provide no evidence for genetically distinct spawning components within the Northwest or Northeast Atlantic.

Conclusions

Therefore, our findings resolved previous uncertainties by confirming the absence of genetically isolated spawning components in each side of the northern Atlantic, thus rejecting homing behaviour and the need to redefine management boundaries in this species. In addition, no further genetic signs of ongoing adaptation were detected in this species.

Chromosome-level genome assembly and annotation of Barbel chub Squaliobarbus curriculus

The barbel chub Squaliobarbus curriculus, is an economically important freshwater fish in China. The fishery production of the wild populations has declined dramatically, making the development of aquaculture urgently needed. However, the lack of high-quality genome has impeded its artificial breeding and genetic breeding. Herein, we present a chromosome-level genome assembly for S. curriculus by combining HiFi sequencing, Hi-C sequencing, Iso-seq and short-reads RNA-seq data. This assembly was 910.27 Mb in size, with a contig N50 length of 34.70 Mb. 99.50% of the assembled sequences were placed onto 24 chromosomes supported by Hi-C contact map. Using Iso-seq and short-reads RNA-seq data, we identified 28,329 protein-coding genes based on three prediction methods. Of these genes, 27,207 genes (96.04%) were functionally annotated to at least one of the six commonly used databases. Additionally, we annotated 2,041 miRNAs, 16,426 tRNAs, 5,488 rRNAs and 1,536 snRNAs in the S. curriculus genome. Overall, the chromosome-level genome of S. curriculus will provide valuable genomic resources for genetic breeding, population genomics, sex-related marker identifications, and other future studies.

Transcriptome Assembly at Single-Cell Resolution with Beaver

Emerging single-cell RNA sequencing techniques (scRNA-seq) has enabled the study of cellular transcriptome heterogeneity, yet accurate reconstruction of full-length transcripts at single-cell resolution remains challenging due to high dropout rates and sparse coverage. While meta-assembly approaches offer promising solutions by integrating information across multiple cells, current methods struggle to balance consensus assembly with cell-specific transcriptional signatures. Here, we present Beaver, a cell-specific transcript assembler designed for short-read scRNA-seq data. Beaver implements a transcript fragment graph to organize individual assemblies and designs an efficient dynamic programming algorithm that searches for candidate full-length transcripts from the graph. Beaver in-corporates two random forest models trained on 51 meticulously engineered features that accurately estimate the likelihood of each candidate transcript being expressed in individual cells. Our experiments, performed using both real and simulated Smart-seq3 scRNA-seq data, firmly show that Beaver substantially outperforms existing meta-assemblers and single-sample assemblers. At the same level of sensitivity, Beaver achieved 32.0%-64.6%, 13.5%-36.6%, and 9.8%-36.3% higher precision in average compared to meta-assemblers Aletsch, TransMeta, and PsiCLASS, respectively, with similar improvements over single-sample assemblers Scallop2 (10.1%-43.6%) and StringTie2 (24.3%-67.0%). Beaver is freely available at https://github.com/Shao-Group/beaver . Scripts that reproduce the experimental results of this manuscript are available at https://github.com/Shao-Group/beaver-test .

Analyzes of pan-genome and resequencing atlas unveil the genetic basis of jujube domestication

Jujube (Ziziphus jujuba Mill.), belonging to the Rhamnaceae family, is gaining increasing prominence as a perennial fruit crop with significant economic and medicinal values. Here, we conduct de novo assembly of four reference-grade genomes, encompassing one wild and three cultivated jujube accessions. We present insights into the population structure, genetic diversity, and genomic variations within a diverse collection of 1059 jujube accessions. Analyzes of the jujube pan-genome, based on our four assemblies and four previously released genomes, reveal extensive genomic variations within domestication-associated regions, potentially leading to the discovery of a candidate gene that regulates flowering and fruit ripening. By leveraging the pan-genome and a large-scale resequencing population, we identify two candidate genes involved in domestication traits, including the seed-setting rate, the bearing-shoot length and the leaf size in jujube. These genomic resources will accelerate evolutionary and functional genomics studies of jujube.

lncRNA-gene network analysis reveals the effects of early maternal nutrition on mineral homeostasis and energy metabolism in the fetal liver transcriptome of beef heifers

Maternal nutrition during pregnancy influences fetal development; however, the regulatory markers of fetal programming across different gestational phases remain underexplored in livestock models. Herein, we investigated the regulatory role of long non-coding RNAs (lncRNAs) on fetal liver gene expression, the impacts of maternal vitamin and mineral supplementation, and the rate of maternal body weight gain during the periconceptual period. To this end, crossbred Angus heifers (n=31) were randomly assigned to a 2×2 factorial design to evaluate the main effects of the rate of weight gain (low gain [LG, avg. daily gain of 0.28 kg/day] vs. moderate gain [MG, avg. daily gain of 0.79 kg/day]) and vitamins and minerals supplementation (VTM vs. NoVTM). On day 83±0.27 of gestation, fetuses were collected for morphometric measurements, and fetal liver was collected for transcriptomic and mineral analyses. The maternal diet significantly affected fetal liver development and mineral reserves. Using an RNA-Seq approach, we identified 320 unique differentially expressed genes (DEGs) across all six comparisons (FDR <0.05). Furthermore, lncRNAs were predicted through the FEELnc pipeline, revealing 99 unique differentially expressed lncRNAs (DELs). The over-represented pathways and biological processes (BPs) were associated with energy metabolism, Wnt signaling, CoA carboxylase activity, and fatty acid metabolism. The DEL-regulated BPs were associated with metal ion transport, pyrimidine metabolism, and classical energy metabolism-related glycolytic, gluconeogenic, and TCA cycle pathways. Our findings suggest that lncRNAs regulate mineral homeostasis- and energy metabolism-related gene networks in the fetal liver in response to early maternal nutrition.

The haplotype-resolved Prymnesium parvum (type B) microalga genome reveals the genetic basis of its fish-killing toxins

The catastrophic loss of aquatic life in the Central European Oder River in 2022, caused by a toxic bloom of the haptophyte microalga Prymnesium parvum (in a wide sense, s.l.), underscores the need to improve our understanding of the genomic basis of the toxin. Previous morphological, phylogenetic, and genomic studies have revealed cryptic diversity within P. parvum s.l. and uncovered three clade-specific (types A, B, and C) prymnesin toxins. Here, we used state-of-the-art long-read sequencing and assembled the first haplotype-resolved diploid genome of a P. parvum type B from the strain responsible for the Oder disaster. Comparative analyses with type A genomes uncovered a genome-size expansion driven by repetitive elements in type B. We also found conserved synteny but divergent evolution in several polyketide synthase (PKS) genes, which are known to underlie toxin production in combination with environmental cues. We identified an approximately 20-kbp deletion in the largest PKS gene of type B that we link to differences in the chemical structure of types A and B prymnesins. Flow cytometry and electron microscopy analyses confirmed diploidy in the Oder River strain and revealed differences to closely related strains in both ploidy and morphology. Our results provide unprecedented resolution of strain diversity in P. parvum s.l. and a better understanding of the genomic basis of toxin variability in haptophytes. The reference-quality genome will enable us to better understand changes in microbial diversity in the face of increasing environmental pressures and provides a basis for strain-level monitoring of invasive Prymnesium in the future.

Heat stress reprograms herbivory-induced defense responses in potato plants

Climate change is predicted to increase the occurrence of extreme weather events such as heatwaves, which may thereby impact the outcome of plant-herbivore interactions. While elevated temperature is known to directly affect herbivore growth, it remains largely unclear if it indirectly influences herbivore performance by affecting the host plant they feed on. In this study, we investigated how transient exposure to high temperature influences plant herbivory-induced defenses at the transcript and metabolic level. To this end, we studied the interaction between potato (Solanum tuberosum) plants and the larvae of the potato tuber moth (Phthorimaea operculella) under different temperature regimes. We found that P. operculella larvae grew heavier on leaves co-stressed by high temperature and insect herbivory than on leaves pre-stressed by herbivory alone. We also observed that high temperature treatments altered phylotranscriptomic patterns upon herbivory, which changed from an evolutionary hourglass pattern, in which transcriptomic responses at early and late time points after elicitation are more variable than the ones in the middle, to a vase pattern. Specifically, transcripts of many herbivory-induced genes in the early and late defense stage were suppressed by HT treatment, whereas those in the intermediate stage peaked earlier. Additionally, we observed that high temperature impaired the induction of jasmonates and defense compounds upon herbivory. Moreover, using jasmonate-reduced (JA-reduced, irAOC) and -elevated (JA-Ile-elevated, irCYP94B3s) potato plants, we showed that high temperature suppresses JA signaling mediated plant-induced defense to herbivore attack. Thus, our study provides evidences on how temperature reprograms plant-induced defense to herbivores.

The Human Developing Cerebral Cortex Is Characterized by an Elevated De Novo Expression of Long Noncoding RNAs in Excitatory Neurons

The outstanding human cognitive capacities are computed in the cerebral cortex, a mammalian-specific brain region and the place of massive biological innovation. Long noncoding RNAs have emerged as gene regulatory elements with higher evolutionary turnover than mRNAs. The many long noncoding RNAs identified in neural tissues make them candidates for molecular sources of cerebral cortex evolution and disease. Here, we characterized the genomic and cellular shifts that occurred during the evolution of the long noncoding RNA repertoire expressed in the developing cerebral cortex and explored putative roles for these long noncoding RNAs in the evolution of the human brain. Using transcriptomics and comparative genomics, we comprehensively annotated the cortical transcriptomes of humans, rhesus macaques, mice, and chickens and classified human cortical long noncoding RNAs into evolutionary groups as a function of their predicted minimal ages. Long noncoding RNA evolutionary groups showed differences in expression levels, splicing efficiencies, transposable element contents, genomic distributions, and transcription factor binding to their promoters. Furthermore, older long noncoding RNAs showed preferential expression in germinative zones, outer radial glial cells, and cortical inhibitory (GABAergic) neurons. In comparison, younger long noncoding RNAs showed preferential expression in cortical excitatory (glutamatergic) neurons, were enriched in primate and human-specific gene co-expression modules, and were dysregulated in neurodevelopmental disorders. These results suggest different evolutionary routes for older and younger cortical long noncoding RNAs, highlighting old long noncoding RNAs as a possible source of molecular evolution of conserved developmental programs; conversely, we propose that the de novo expression of primate- and human-specific young long noncoding RNAs is a putative source of molecular evolution and dysfunction of cortical excitatory neurons, warranting further investigation.

Accurate assembly of multiple RNA-seq samples with Aletsch

MOTIVATION

High-throughput RNA sequencing has become indispensable for decoding gene activities, yet the challenge of reconstructing full-length transcripts persists. Traditional single-sample assemblers frequently produce fragmented transcripts, especially in single-cell RNA-seq data. While algorithms designed for assembling multiple samples exist, they encounter various limitations.

RESULTS

We present Aletsch, a new assembler for multiple bulk or single-cell RNA-seq samples. Aletsch incorporates several algorithmic innovations, including a "bridging" system that can effectively integrate multiple samples to restore missed junctions in individual samples, and a new graph-decomposition algorithm that leverages "supporting" information across multiple samples to guide the decomposition of complex vertices. A standout feature of Aletsch is its application of a random forest model with 50 well-designed features for scoring transcripts. We demonstrate its robust adaptability across different chromosomes, datasets, and species. Our experiments, conducted on RNA-seq data from several protocols, firmly demonstrate Aletsch's significant outperformance over existing meta-assemblers. As an example, when measured with the partial area under the precision-recall curve (pAUC, constrained by precision), Aletsch surpasses the leading assemblers TransMeta by 22.9%-62.1% and PsiCLASS by 23.0%-175.5% on human datasets.

AVAILABILITY AND IMPLEMENTATION

Aletsch is freely available at https://github.com/Shao-Group/aletsch. Scripts that reproduce the experimental results of this manuscript is available at https://github.com/Shao-Group/aletsch-test.

Identification of antennal alternative splicing by combining genome and full-length transcriptome analysis in Bactrocera dorsalis

Alternative splicing is an essential post-transcriptional regulatory mechanism that diversifies gene function by generating multiple protein isoforms from a single gene and act as a crucial role in insect environmental adaptation. Olfaction, a key sense for insect adaptation, relies heavily on the antennae, which are the primary olfactory organs expressing most of the olfactory genes. Despite the extensive annotation of olfactory genes within insect antennal tissues facilitated by high-throughput sequencing technology advancements, systematic analyses of alternative splicing are still relatively less. In this study, we focused on the oriental fruit fly (Bactrocera dorsalis), a significant pest of fruit crops. We performed a detailed analysis of alternative splicing in its antennae by utilizing the full-length transcriptome of its antennal tissue and the insect's genome. The results revealed 8600 non-redundant full-length transcripts identified in the oriental fruit fly antennal full-length transcriptome, spanning 4,145 gene loci. Over 40% of these loci exhibited multiple isoforms. Among these, 161 genes showed sex-biased isoform switching, involving seven different types of alternative splicing. Notably, events involving alternative transcription start sites (ATSS) and alternative transcription termination sites (ATTS) were the most common. Of all the genes undergoing ATSS and ATTS alternative splicing between male and female, 32 genes were alternatively spliced in protein coding regions, potentially affecting protein function. These genes were categorized based on the length of the sex-biased isoforms, with the highest difference in isoform fraction (dIF) associated with the ATSS type, including genes such as BdorABCA13, BdorCAT2, and BdorTSN3. Additionally, transcription factor binding sites for doublesex were identified upstream of both BdorABCA13 and BdorCAT2. Besides being expressed in the antennal tissues, BdorABCA13 and BdorCAT2 are also expressed in the mouthparts, legs, and genitalia of both female and male adults, suggesting their functional diversity. This study reveals alternative splicing events in the antennae of Bactrophora dorsalis from two aspects: odorant receptor genes and other types of genes expressed in the antennae. This study not only provides a research foundation for understanding the regulation of gene function by alternative splicing in the oriental fruit fly but also offers new insights for utilizing olfaction-based behavioral manipulation techniques to manage this pest.

PvGeneExpDB: An integrative gene expression database for in-depth understanding on the Pacific white shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei) is a high-valued economic farming species. With the development of high-throughput sequencing technology, cumulative large-scale transcriptomic studies have been revealing molecular landscape of various biological conditions including genetic selection, breeding, evolution, disease landscape, etc. However, no single experiment or databases allow thorough investigations of transcriptomic dynamics for these progressions. Meanwhile, the available datasets are often scattered and lack management. Here, we have established PvGeneExpDB, the first gene expression database for L. vannamei (www.bio-marine-scau.com/pv_ex/), which encompasses gene expression profiles, differential expression, and co-expression analyses under various biological conditions. Based on the analyses of 7 datasets, which include 53 samples with accurate and detailed records, PvGeneExpDB identifies 20,599 novel transcripts, shows expression profiles of a total of 20,817 genes, and implements Gene Ontology (GO) reconstruction of 76.7 % of these genes. Besides, 26 co-expressed groups were first identified by large-scale, cross-sample Weighted Gene Co-expression Network Analysis (WGCNA). By integrating the gene expression data in the database, our goal is to deepen the biological understanding of L. vannamei.

Diversity and host interaction of the gut microbiota in specific pathogen-free pigs

Pigs are widely used as animal models in various studies related to humans. The interaction between the gut microbiota and the host has significant effects on the host's health and disease status. However, although there have been many studies investigating the pig gut microbiota, the findings have been inconsistent due to variations in rearing conditions. Interactions between the gut microbiota and host have not been fully explored in pigs. Specific pathogen-free (SPF) pigs are ideal non-primate large animals to study the interactions between the gut microbiota and the host. In this study, we performed high-throughput sequencing analysis of the gut microbiota and the gut tissue transcriptome of six SPF pigs to provide a systematic understanding of the composition, function, and spatial distribution of gut microbiota in SPF pigs. We identified significant differences in microbial diversity and functionality among different gastrointestinal tract sites. Metagenomics data analysis revealed significant differences in alpha diversity and beta diversity of microbiota in different gastrointestinal sites of SPF pigs. Additionally, transcriptomic data indicated significant differences in gene expression as well as KEGG and GO functional enrichment between the small intestine and large intestine. Furthermore, by combining microbial metagenomics and host transcriptomics analyses, specific correlations were found between gut microbiota and host genes. These included a negative correlation between the TCN1 gene and Prevotella dentalis, possibly related to bacterial metabolic pathways involving vitamin B12, and a positive correlation between the BDH1 gene and Roseburia hominis, possibly because both are involved in fatty acid metabolism. These findings lay the groundwork for further exploration of the co-evolution between the microbiota and the host, specifically in relation to nutrition, metabolism, and immunity. In conclusion, we have elucidated the diversity of the gut microbiota in SPF pigs and conducted a detailed investigation into the interactions between the gut microbiota and host gene expression. These results contribute to our understanding of the intricate dynamics between the gut microbiota and the host, offering important references for advancements in life science research, bioproduct production, and sustainable development in animal husbandry.

A chromosome-level haplotype-resolved genome assembly of oriental tobacco budworm (Helicoverpa assulta)

Oriental tobacco budworm (Helicoverpa assulta) and cotton bollworm (Helicoverpa armigera) are two closely related species within the genus Helicoverpa. They have similar appearances and consistent damage patterns, often leading to confusion. However, the cotton bollworm is a typical polyphagous insect, while the oriental tobacco budworm belongs to the oligophagous insects. In this study, we used Nanopore, PacBio, and Illumina platforms to sequence the genome of H. assulta and used Hifiasm to create a haplotype-resolved draft genome. The Hi-C technique helped anchor 33 primary contigs to 32 chromosomes, including two sex chromosomes, Z and W. The final primary haploid genome assembly was approximately 415.19 Mb in length. BUSCO analysis revealed a high degree of completeness, with 99.0% gene coverage in this genome assembly. The repeat sequences constituted 38.39% of the genome assembly, and we annotated 17093 protein-coding genes. The high-quality genome assembly of the oriental tobacco budworm serves as a valuable genetic resource that enhances our comprehension of how they select hosts in a complex odour environment. It will also aid in developing an effective control policy.

SAFB restricts contact domain boundaries associated with L1 chimeric transcription

Long interspersed element-1 (LINE-1 or L1) comprises 17% of the human genome, continuously generates genetic variations, and causes disease in certain cases. However, the regulation and function of L1 remain poorly understood. Here, we uncover that L1 can enrich RNA polymerase IIs (RNA Pol IIs), express L1 chimeric transcripts, and create contact domain boundaries in human cells. This impact of L1 is restricted by a nuclear matrix protein scaffold attachment factor B (SAFB) that recognizes transcriptionally active L1s by binding L1 transcripts to inhibit RNA Pol II enrichment. Acute inhibition of RNA Pol II transcription abolishes the domain boundaries associated with L1 chimeric transcripts, indicating a transcription-dependent mechanism. Deleting L1 impairs domain boundary formation, and L1 insertions during evolution have introduced species-specific domain boundaries. Our data show that L1 can create RNA Pol II-enriched regions that alter genome organization and that SAFB regulates L1 and RNA Pol II activity to preserve gene regulation.

Leaf morphology related genes revealed by integrating Pan-transcriptome, GWAS and eQTL analyses in a Liriodendron population

Leaf plays an indispensable role in plant development and growth. Although many known genes related to leaf morphology development have been identified, elucidating the complex genetic basis of leaf morphological traits remains a challenge. Liriodendron plants are common ornamental trees due to their unique leaf shapes, while the molecular mechanism underlying Liriodendron leaf morphogenesis has remained unknown. Herein, we firstly constructed a population-level pan-transcriptome of Liriodendron from 81 accessions to explore the expression presence or absence variations (ePAVs), global expression differences at the population level, as well as differentially expressed genes (DEGs) between the Liriodendron chinense and Liriodendron tulipifera accessions. Subsequently, we integrated a genome-wide association study (GWAS), expression quantitative trait loci (eQTL), and transcriptome-wide association study (TWAS) to identify candidate genes related to leaf morphology. Through GWAS analysis, we identified 18 and 17 significant allelic loci in the leaf size and leaf shape modules, respectively. In addition, we discerned 16 candidate genes in relation to leaf morphological traits via TWAS. Further, integrating the co-localization results of GWAS and eQTL, we determined two regulatory hotspot regions, hot88 and hot758, related to leaf size and leaf shape, respectively. Finally, co-expression analysis, eQTL, and linkage mapping together demonstrated that Lchi_4g10795 regulate their own expression levels through cis-eQTL to affect the expression of downstream genes and cooperatively participate in the development of Liriodendron leaf morphology. These findings will improve our understanding of the molecular regulatory mechanism of Liriodendron leaf morphogenesis and will also accelerate molecular breeding of Liriodendron.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

ADNP modulates SINE B2-derived CTCF-binding sites during blastocyst formation in mice

CTCF is crucial for chromatin structure and transcription regulation in early embryonic development. However, the kinetics of CTCF chromatin occupation in preimplantation embryos have remained unclear. In this study, we used CUT&RUN technology to investigate CTCF occupancy in mouse preimplantation development. Our findings revealed that CTCF begins binding to the genome prior to zygotic genome activation (ZGA), with a preference for CTCF-anchored chromatin loops. Although the majority of CTCF occupancy is consistently maintained, we identified a specific set of binding sites enriched in the mouse-specific short interspersed element (SINE) family B2 that are restricted to the cleavage stages. Notably, we discovered that the neuroprotective protein ADNP counteracts the stable association of CTCF at SINE B2-derived CTCF-binding sites. Knockout of Adnp in the zygote led to impaired CTCF binding signal recovery, failed deposition of H3K9me3, and transcriptional derepression of SINE B2 during the morula-to-blastocyst transition, which further led to unfaithful cell differentiation in embryos around implantation. Our analysis highlights an ADNP-dependent restriction of CTCF binding during cell differentiation in preimplantation embryos. Furthermore, our findings shed light on the functional importance of transposable elements (TEs) in promoting genetic innovation and actively shaping the early embryo developmental process specific to mammals.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Pig H3K4me3, H3K27ac, and gene expression profiles reveal reproductive tissue-specific activity of transposable elements

Regulatory sequences and transposable elements (TEs) account for a large proportion of the genomic sequences of species; however, their roles in gene transcription, especially tissue-specific expression, remain largely unknown. Pigs serve as an excellent animal model for studying genomic sequence biology due to the extensive diversity among their wild and domesticated populations. Here, we conducted an integrated analysis using H3K27ac ChIP-seq, H3K4me3 ChIP-seq, and RNA-seq data from 10 different tissues of seven fetuses and eight closely related adult pigs. We aimed to annotate the regulatory elements and TEs to elucidate their associations with histone modifications and mRNA expression across different tissues and developmental stages. Based on correlation analysis between mRNA expression and H3K27ac and H3K4me3 peak activity, results indicated that H3K27ac exhibited stronger associations with gene expression than H3K4me3. Furthermore, 1.45% of TEs overlapped with either the H3K27ac or H3K4me3 peaks, with the majority displaying tissue-specific activity. Notably, a TE subfamily (LTR4C_SS), containing binding motifs for SIX1 and SIX4, showed specific enrichment in the H3K27ac peaks of the adult and fetal ovaries. RNA-seq analysis also revealed widespread expression of TEs in the exons or promoters of genes, including 4 688 TE-containing transcripts with distinct development stage-specific and tissue-specific expression. Of note, 1 967 TE-containing transcripts were enriched in the testes. We identified a long terminal repeat (LTR), MLT1F1, acting as a testis-specific alternative promoter in SRPK2 (a cell cycle-related protein kinase) in our pig dataset. This element was also conserved in humans and mice, suggesting either an ancient integration of TEs in genes specifically expressed in the testes or parallel evolutionary patterns. Collectively, our findings demonstrate that TEs are deeply embedded in the genome and exhibit important tissue-specific biological functions, particularly in the reproductive organs.

A chromosome-level reference genome for the common octopus, Octopus vulgaris (Cuvier, 1797)

Cephalopods are emerging animal models and include iconic species for studying the link between genomic innovations and physiological and behavioral complexities. Coleoid cephalopods possess the largest nervous system among invertebrates, both for cell counts and brain-to-body ratio. Octopus vulgaris has been at the center of a long-standing tradition of research into diverse aspects of cephalopod biology, including behavioral and neural plasticity, learning and memory recall, regeneration, and sophisticated cognition. However, no chromosome-scale genome assembly was available for O. vulgaris to aid in functional studies. To fill this gap, we sequenced and assembled a chromosome-scale genome of the common octopus, O. vulgaris. The final assembly spans 2.8 billion basepairs, 99.34% of which are in 30 chromosome-scale scaffolds. Hi-C heatmaps support a karyotype of 1n = 30 chromosomes. Comparisons with other octopus species' genomes show a conserved octopus karyotype and a pattern of local genome rearrangements between species. This new chromosome-scale genome of O. vulgaris will further facilitate research in all aspects of cephalopod biology, including various forms of plasticity and the neural machinery underlying sophisticated cognition, as well as an understanding of cephalopod evolution.

Characterization of H3K9me3 and DNA methylation co-marked CpG-rich regions during mouse development

BACKGROUND

H3K9me3 and DNA methylation co-marked CpG-rich regions (CHMs) are functionally important in mouse pre-implantation embryos, but their characteristics in other biological processes are still largely unknown.

RESULTS

In this study, we performed a comprehensive analysis to characterize CHMs during 6 mouse developmental processes, identifying over 2,600 CHMs exhibiting stable co-mark of H3K9me3 and DNA methylation patterns at CpG-rich regions. We revealed the distinctive features of CHMs, including elevated H3K9me3 signals and a significant presence in euchromatin and the potential role in silencing younger long terminal repeats (LTRs), especially in some ERVK subfamilies. The results highlight the distinct nature of universal CHMs compared to CpG-rich nonCHMs in terms of location, LTR enrichment, and DNA sequence features, enhancing our understanding of CpG-rich regions' regulatory roles.

CONCLUSIONS

This study characterizes the features of CHMs in multiple developmental processes and broadens our understanding of the regulatory roles of CpG-rich regions.

Protocol for transcriptome assembly by the TransBorrow algorithm

High-throughput RNA-seq enables comprehensive analysis of the transcriptome for various purposes. However, this technology generally generates massive amounts of sequencing reads with a shorter read length. Consequently, fast, accurate, and flexible tools are needed for assembling raw RNA-seq data into full-length transcripts and quantifying their expression levels. In this protocol, we report TransBorrow, a novel transcriptome assembly software specifically designed for short RNA-seq reads. TransBorrow is employed in conjunction with a splice-aware alignment tool (e.g. Hisat2 and Star) and some other transcriptome assembly tools (e.g. StringTie, Cufflinks, and Scallop). The protocol encompasses all necessary steps, starting from downloading and processing raw sequencing data to assembling the full-length transcripts and quantifying their expressed abundances. The execution time of the protocol may vary depending on the sizes of processed datasets and computational platforms.

Chromosome-scale assemblies of S. malaccense, S. aqueum, S. jambos, and S. syzygioides provide insights into the evolution of Syzygium genomes

Syzygium is a large and diverse tree genus in the Myrtaceae family. Genome assemblies for clove (Syzygium aromaticum, 370 Mb) and sea apple (Syzygium grande, 405 Mb) provided the first insights into the genomic features and evolution of the Syzygium genus. Here, we present additional de novo chromosome-scale genome assemblies for Syzygium malaccense, Syzygium aqueum, Syzygium jambos, and Syzygium syzygioides. Genome profiling analyses show that S. malaccense, like S. aromaticum and S. grande, is diploid (2n = 2x = 22), while the S. aqueum, S. jambos, and S. syzygioides specimens are autotetraploid (2n = 4x = 44). The genome assemblies of S. malaccense (430 Mb), S. aqueum (392 Mb), S. jambos (426 Mb), and S. syzygioides (431 Mb) are highly complete (BUSCO scores of 98%). Comparative genomics analyses showed conserved organization of the 11 chromosomes with S. aromaticum and S. grande, and revealed species-specific evolutionary dynamics of the long terminal repeat retrotransposon elements belonging to the Gypsy and Copia lineages. This set of Syzygium genomes is a valuable resource for future structural and functional comparative genomic studies on Myrtaceae species.

Chromosome-level assembly and annotation of the Xyrichtys novacula (Linnaeus, 1758) genome

The pearly razorfish (Xyrichtys novacula), commonly known as raor in the Balearic Islands, is a wrasse within the family Labridae. This fish species has particular biological and socio-cultural characteristics making it an ideal model organism in the fields of behavioural ecology, molecular ecology and conservation biology. In this study, we present the first annotated chromosome-level assembly for this species. Sequencing involved a combination of long reads with Oxford Nanopore Technologies, Illumina paired-end short reads (2 × 151 bp), Hi-C and RNA-seq from different tissues. The nuclear genome assembly has a scaffold N50 of 34.33 Mb, a total assembly span of 775.53 Mb and 99.63% of the sequence assembled into 24 superscaffolds, consistent with its known karyotype. Quality metrics revealed a consensus accuracy (QV) of 42.92 and gene completeness > 98%. The genome annotation resulted in 26,690 protein-coding genes and 12,737 non-coding transcripts. The coding regions encoded 39,613 unique protein products, 93% of them with assigned function. Overall, the publication of the X. novacula's reference genome will broaden the scope and impact of genomic research conducted on this iconic and colourful species.

Dysregulation of murine long noncoding single-cell transcriptome in nonalcoholic steatohepatitis and liver fibrosis

LncRNAs comprise a heterogeneous class of RNA-encoding genes typified by low expression, nuclear enrichment, high tissue-specificity, and functional diversity, but the vast majority remain uncharacterized. Here, we assembled the mouse liver noncoding transcriptome from >2000 bulk RNA-seq samples and discovered 48,261 liver-expressed lncRNAs, a majority novel. Using these lncRNAs as a single-cell transcriptomic reference set, we elucidated lncRNA dysregulation in mouse models of high fat diet-induced nonalcoholic steatohepatitis and carbon tetrachloride-induced liver fibrosis. Trajectory inference analysis revealed lncRNA zonation patterns across the liver lobule in each major liver cell population. Perturbations in lncRNA expression and zonation were common in several disease-associated liver cell types, including nonalcoholic steatohepatitis-associated macrophages, a hallmark of fatty liver disease progression, and collagen-producing myofibroblasts, a central feature of liver fibrosis. Single-cell-based gene regulatory network analysis using bigSCale2 linked individual lncRNAs to specific biological pathways, and network-essential regulatory lncRNAs with disease-associated functions were identified by their high network centrality metrics. For a subset of these lncRNAs, promoter sequences of the network-defined lncRNA target genes were significantly enriched for lncRNA triplex formation, providing independent mechanistic support for the lncRNA-target gene linkages predicted by the gene regulatory networks. These findings elucidate liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. A subset of the liver disease-associated regulatory lncRNAs identified have human orthologs and are promising candidates for biomarkers and therapeutic targets.

Chromosome-level genome assembly of Lilford's wall lizard, Podarcis lilfordi (Günther, 1874) from the Balearic Islands (Spain)

The Mediterranean lizard Podarcis lilfordi is an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high-quality chromosome-level assembly and annotation of the P. lilfordi genome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes translating into 38,615 proteins. Comparison to the genome of the related species Podarcis muralis revealed substantial similarity in genome size, annotation metrics, repeat content, and a strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics.

The landscape of the long non-coding RNAs in developing mouse retinas

BACKGROUND

The long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes. Nevertheless, a global view of its expression and function in the mouse retina, a crucial model for neurogenesis study, still needs to be made available.

RESULTS

Herein, by integrating the established gene models and the result from ab initio prediction using short- and long-read sequencing, we characterized 4,523 lncRNA genes (MRLGs) in developing mouse retinas (from the embryonic day of 12.5 to the neonatal day of P28), which was so far the most comprehensive collection of retinal lncRNAs. Next, derived from transcriptomics analyses of different tissues and developing retinas, we found that the MRLGs were highly spatiotemporal specific in expression and played essential roles in regulating the genesis and function of mouse retinas. In addition, we investigated the expression of MRLGs in some mouse mutants and revealed that 97 intergenic MRLGs might be involved in regulating differentiation and development of retinal neurons through Math5, Isl1, Brn3b, NRL, Onecut1, or Onecut2 mediated pathways.

CONCLUSIONS

In summary, this work significantly enhanced our knowledge of lncRNA genes in mouse retina development and provided valuable clues for future exploration of their biological roles.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

Comprehensive RNA-Seq Analysis Pipeline for Non-Model Organisms and Its Application in Schmidtea mediterranea

RNA sequencing (RNA-seq) is a high-throughput technology that provides in-depth information on transcriptome. The advancement and dropping costs of RNA sequencing, accompanied by more available reference genomes for different species, make transcriptome analysis in non-model organisms possible. Current obstacles in analyzing RNA-seq data include a lack of functional annotation, which may complicate the process of linking genes to corresponding functions. Here, we provide a one-stop RNA-seq analysis pipeline, PipeOne-NM, for transcriptome functional annotation, non-coding RNA identification, and transcripts alternative splicing analysis of non-model organisms, intended for use with Illumina platform-based RNA-seq data. We performed PipeOne-NM on 237 Schmidtea mediterranea RNA-seq runs and assembled a transcriptome with 84,827 sequences from 49,320 genes, identifying 64,582 mRNA from 35,485 genes, 20,217 lncRNA from 17,084 genes, and 3481 circRNAs from 1103 genes. In addition, we performed a co-expression analysis of lncRNA and mRNA and identified that 1319 lncRNA co-express with at least one mRNA. Further analysis of samples from S. mediterranea sexual and asexual strains revealed the role of sexual reproduction in gene expression profiles. Samples from different parts of asexual S. mediterranea revealed that differential expression profiles of different body parts correlated with the function of conduction of nerve impulses. In conclusion, PipeOne-NM has the potential to provide comprehensive transcriptome information for non-model organisms on a single platform.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Complexity and regulation of age-dependent alternative splicing in Brachypodium distachyon

Alternative splicing (AS) is a gene regulatory mechanism that generates multiple transcripts of the same gene precursor by the spliceosome complex, promoting messenger RNA complexity and proteome diversity. Although AS is extensively studied in response to environmental stresses, whether it mediates age-dependent development and how it is adjusted by growth transitions are largely unknown. Here, we comprehensively explored the AS landscape at different development stages in the grass model plant Brachypodium (Brachypodium distachyon). We identified abundant coding genes and non-coding transcripts subject to dynamic AS regulation during juvenile, adult, and reproductive transitions. Moreover, we revealed that SC35-LIKE SPLICING FACTOR 33 (SCL33), a serine/arginine-rich splicing factor in spliceosomes, plays a redundant and antagonistic role with its putative paralog, SCL33L, in regulating intron assembly across distinct developmental stages. In addition, we determined global AS variations in microRNA156 (miR156)-overproducing plants, in which growth transitions are delayed, and found that SPLs were regulated by miR156 in intron retention alteration in addition to mRNA clearance and translation inhibition manners. Finally, we demonstrated a complex regulatory process of age-dependent AS events in B. distachyon that were coincidently or separately regulated by miR156 and SCL33/SCL33L. These results illustrate a substantial -machinery of AS that mediates phase transitions in plants.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Genome Resources of Verticillium dahliae VdGL16: The Causal Agent of Vascular Wilt on the Invasive Species Ailanthus altissima

Verticillium species are known as plant pathogens responsible for wilt diseases in a large variety of dicotyledon plants and crops in many parts of the world. Here we present the draft genome sequence of Verticillium dahliae Kleb. (strain VdGL16) isolated in Italy from the invasive alien species Ailanthus altissima (Mill.; commonly known as tree-of-heaven) showing Verticillium wilt symptoms. The comparison between the newly sequenced genome with those publicly available revealed candidate genes putatively involved in pathogenicity. The genome represents a new useful source for future research on Verticillium genetics and biology as well as research on novel approaches in the control of A. altissima.

Cell-specific clock-controlled gene expression program regulates rhythmic fiber cell growth in cotton

BACKGROUND

The epidermis of cotton ovule produces fibers, the most important natural cellulose source for the global textile industry. However, the molecular mechanism of fiber cell growth is still poorly understood.

RESULTS

Here, we develop an optimized protoplasting method, and integrate single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) to systematically characterize the cells of the outer integument of ovules from wild type and fuzzless/lintless (fl) cotton (Gossypium hirsutum). By jointly analyzing the scRNA-seq data from wildtype and fl, we identify five cell populations including the fiber cell type and construct the development trajectory for fiber lineage cells. Interestingly, by time-course diurnal transcriptomic analysis, we demonstrate that the primary growth of fiber cells is a highly regulated circadian rhythmic process. Moreover, we identify a small peptide GhRALF1 that circadian rhythmically controls fiber growth possibly through oscillating auxin signaling and proton pump activity in the plasma membrane. Combining with scATAC-seq, we further identify two cardinal cis-regulatory elements (CREs, TCP motif, and TCP-like motif) which are bound by the trans factors GhTCP14s to modulate the circadian rhythmic metabolism of mitochondria and protein translation through regulating approximately one third of genes that are highly expressed in fiber cells.

CONCLUSIONS

We uncover a fiber-specific circadian clock-controlled gene expression program in regulating fiber growth. This study unprecedentedly reveals a new route to improve fiber traits by engineering the circadian clock of fiber cells.

Molecular signatures of diapause in the Asian longhorned beetle: Gene expression

Most previous studies on gene expression during insect diapause do not address among-tissue variation in physiological processes. We measured transcriptomic changes during larval diapause in the Asian longhorned beetle, Anoplophora glabripennis (Coleoptera: Cerambycidae). We conducted RNA-seq on fat body, the supraesophageal ganglion, midgut, hindgut, and Malpighian tubules during pre-diapause, diapause maintenance, post-diapause quiescence, and post-diapause development. We observed a small, but consistent, proportion of genes within each gene expression profile that were shared among tissues, lending support for a core set of diapause-associated genes whose expression is tissue-independent. We evaluated the overarching hypotheses that diapause would be associated with cell cycle arrest, developmental arrest, and increased stress tolerance and found evidence of repressed TOR and insulin signaling, reduced cell cycle activity and increased capacity of stress response via heat shock protein expression and remodeling of the cytoskeleton. However, these processes varied among tissues, with the brain and fat body appearing to maintain higher levels of cellular activity during diapause than the midgut or Malpighian tubules. We also observed temperature-dependent changes in gene expression during diapause maintenance, particularly in genes related to the heat shock response and MAPK, insulin, and TOR signaling pathways. Additionally, we provide evidence for epigenetic reorganization during the diapause/post-diapause quiescence transition and expression of genes involved in post-translational modification, highlighting the need for investigations of the protein activity of these candidate genes and processes. We conclude that diapause development is coordinated via diverse tissue-specific gene expression profiles and that canonical diapause phenotypes vary among tissues.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

Chromosome genome assembly for the meagre, Argyrosomus regius, reveals species adaptations and sciaenid sex-related locus evolution

The meagre, Argyrosomus regius, has recently become a species of increasing economic interest for the Mediterranean aquaculture and there is ongoing work to boost production efficiency through selective breeding. Access to the complete genomic sequence will provide an essential resource for studying quantitative trait-associated loci and exploring the genetic diversity of different wild populations and aquaculture stocks in more detail. Here, we present the first complete genome for A. regius, produced through a combination of long and short read technologies and an efficient in-house developed pipeline for assembly and polishing. Scaffolding using previous linkage map data allowed us to reconstruct a chromosome level assembly with high completeness, complemented with gene annotation and repeat masking. The 696 Mb long assembly has an N50 = 27.87 Mb and an L50 = 12, with 92.85% of its length placed in 24 chromosomes. We use this new resource to study the evolution of the meagre genome and other Sciaenids, via a comparative analysis of 25 high-quality teleost genomes. Combining a rigorous investigation of gene duplications with base-wise conservation analysis, we identify candidate loci related to immune, fat metabolism and growth adaptations in the meagre. Following phylogenomic reconstruction, we show highly conserved synteny within Sciaenidae. In contrast, we report rapidly evolving syntenic rearrangements and gene copy changes in the sex-related dmrt1 neighbourhood in meagre and other members of the family. These novel genomic datasets and findings will add important new tools for aquaculture studies and greatly facilitate husbandry and breeding work in the species.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

Analysis of mRNA and lncRNA Expression Profiles of Breast Muscle during Pigeon (Columbalivia) Development

The breast muscle is essential for flight and determines the meat yield and quality of the meat type in pigeons. At present, studies about long non-coding RNA (lncRNA) expression profiles in skeletal muscles across the postnatal development of pigeons have not been reported. Here, we used transcriptome sequencing to examine the White-King pigeon breast muscle at four different ages (1 day, 14 days, 28 days, and 2 years old). We identified 12,918 mRNAs and 9158 lncRNAs (5492 known lncRNAs and 3666 novel lncRNAs) in the breast muscle, and 7352 mRNAs and 4494 lncRNAs were differentially expressed in the process of development. We found that highly expressed mRNAs were mainly related to cell-basic and muscle-specific functions. Differential expression and time-series analysis showed that differentially expressed genes were primarily associated with muscle development and functions, blood vessel development, cell cycle, and energy metabolism. To further predict the possible role of lncRNAs, we also conducted the WGCNA and trans/cis analyses. We found that differentially expressed lncRNAs such as lncRNA-LOC102093252, lncRNA-G12653, lncRNA-LOC110357465, lncRNA-G14790, and lncRNA-LOC110360188 might respectively target UBE2B, Pax7, AGTR2, HDAC1, Sox8 and participate in the development of the muscle. Our study provides a valuable resource for studying the lncRNAs and mRNAs of pigeon muscles and for improving the understanding of molecular mechanisms in muscle development.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

The landscape of hervRNAs transcribed from human endogenous retroviruses across human body sites

BACKGROUND

Human endogenous retroviruses (HERVs), the remnants of ancient retroviruses, account for 8% of the human genome, but most have lost their transcriptional abilities under physiological conditions. However, mounting evidence shows that several expressed HERVs do exert biological functions. Here, we systematically characterize physiologically expressed HERVs and examine whether they may give insight into the molecular fundamentals of human development and disease.

RESULTS

We systematically identify 13,889 expressed HERVs across normal body sites and demonstrate that they are expressed in body site-specific patterns and also by sex, ethnicity, and age. Analyzing cis-ERV-related quantitative trait loci, we find that 5435 hervRNAs are regulated by genetic variants. Combining this with a genome-wide association study, we elucidate that the dysregulation of expressed HERVs might be associated with various complex diseases, particularly neurodegenerative and psychiatric diseases. We further find that physiologically activated hervRNAs are associated with histone modifications rather than DNA demethylation.

CONCLUSIONS

Our results present a locus-specific landscape of physiologically expressed hervRNAs, which represent a hidden layer of genetic architecture in development and disease.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.