Improved transcriptome assembly using a hybrid of long and short reads with StringTie

Genome assembly of Stewartia sinensis reveals origin and evolution of orphan genes in Theaceae

Orphan genes play crucial roles in diverse biological processes, but the evolutionary trajectories and functional divergence remain largely unexplored. The Theaceae family, including the economically and culturally important tea plant, offers a distinctive model to examine these aspects. Here, we integrated Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to decode a pseudo-chromosomal genome assembly of Stewartia sinensis, from the earliest-diverging tribe of Theaceae, spanning 2.95 Gb. Comparative genomic analysis revealed the absence of recent whole-genome duplication events in the Theaceae ancestor, highlighting tandem duplications as the predominant mechanism of gene expansion. We identified 31,331 orphan genes, some of which appear to have ancient origins, suggesting early emergence with frequent gains and losses, while others seem more specific and recent. Notably, orphan genes are distinguished by shorter lengths, fewer exons and functional domains compared to genes that originate much earlier, like transcription factors. Moreover, tandem duplication contributes significantly to the adaptive evolution and characteristic diversity of Theaceae, and it is also a major mechanism driving the origination of orphan genes. This study illuminates the evolutionary dynamics of orphan genes, providing a valuable resource for understanding the origin and evolution of tea plant flavor and enhancing genetic breeding efforts.

Long-read RNA sequencing atlas of human microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. Mapping the genetics of gene expression in human microglia has identified several loci associated with disease-associated genetic variants in microglia-specific regulatory elements. However, identifying genetic effects on splicing is challenging because of the use of short sequencing reads. Here, we present the isoform-centric microglia genomic atlas (isoMiGA), which leverages long-read RNA sequencing to identify 35,879 novel microglia isoforms. We show that these isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ancestry meta-analysis of 555 human microglia short-read RNA sequencing samples from 391 donors, and found associations with genetic risk loci in Alzheimer's and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice-site usage.

Genomic, transcriptomic, and metabolomic analyses reveal convergent evolution of oxime biosynthesis in Darwin's orchid

Angraecum sesquipedale, also known as Darwin's orchid, possesses an exceptionally long nectar spur. Charles Darwin predicted the orchid to be pollinated by a hawkmoth with a correspondingly long proboscis, later identified as Xanthopan praedicta. In this plant-pollinator interaction, the A. sesquipedale flower emits a complex blend of scent compounds dominated by diurnally regulated oximes (R1R2C = N-OH) to attract crepuscular and nocturnal pollinators. The molecular mechanism of oxime biosynthesis remains unclear in orchids. Here, we present the chromosome-level genome of A. sesquipedale. The haploid genome size is 2.10 Gb and represents 19 pseudochromosomes. Cytochrome P450 encoding genes of the CYP79 family known to be involved in oxime biosynthesis in seed plants are not present in the A. sesquipedale genome nor the genomes of other members of the orchid family. Metabolomic analysis of the A. sesquipedale flower revealed a substantial release of oximes at dusk during the blooming stage. By integrating metabolomic and transcriptomic correlation approaches, flavin-containing monooxygenases (FMOs) encoded by six tandem-repeat genes in the A. sesquipedale genome are identified as catalyzing the formation of oximes present. Further in vitro and in vivo assays confirm the function of FMOs in the oxime biosynthesis. We designate these FMOs as orchid oxime synthases 1-6. The evolutionary aspects related to the CYP79 gene losses and neofunctionalization of FMO-catalyzed biosynthesis of oximes in Darwin's orchid provide new insights into the convergent evolution of biosynthetic pathways.

Chromosome-level genome assembly of flathead asp (Pseudaspius leptocephalus)

Pseudaspius leptocephalus is an economically and ecologically important fish species endemic to the Heilongjiang River Basin. Sequencing of the complete genome of P. leptocephalus could provide new insights with implications for the breeding and conservation of this species. The aim of this study was to obtain the chromosome-level genome of P. leptocephalus using PacBio and Hi-C sequencing. PacBio sequencing yielded an assembly of 889.21 Mb, with a contig N50 size of 10.9 Mb. Hi-C sequencing data were used for chromosome mapping and ultimately yielded 25 chromosome sequences. The success rate of chromosome mapping was 95.0%, with a BUSCO evaluation integrity score of 97.3%. A total of 25,531 protein-coding genes were predicted and 24,787 genes were functionally annotated. The BUSCO evaluation integrity score for the predicted protein-coding genes was 93.7%. Repeat sequences accounted for 39.5% of the total length of the genome. This work could aid the artificial breeding of P. leptocephalus, promote the conservation of biological resources, and provide a basis for future genomic and evolutionary studies of this species.

Evidence for evolution of a new sex chromosome within the haploid-dominant Marchantiales plant lineage

Sex chromosomes have evolved independently in numerous lineages across the Tree of Life, in both diploid-dominant species, including many animals and plants, and the less studied haploid-dominant plants and algae. Strict genetic sex determination ensures that individuals reproduce by outcrossing. However, species with separate sexes (termed dioecy in diploid plants, and dioicy in haploid plants) may sometimes evolve different sex systems, and become monoicous, with the ability to self-fertilize. Here, we studied dioicy-monoicy transitions in the ancient liverwort haploid-dominant plant lineage, using three telomere-to-telomere gapless chromosome-scale reference genome assemblies from the Ricciaceae group of Marchantiales. Ancestral liverworts are believed to have been dioicous, with U and V chromosomes (chromosome 9) determining femaleness and maleness, respectively. We confirm the finding that monoicy in Ricciocarpos natans evolved from a dioicous ancestor, and most ancestrally U chromosomal genes have been retained on autosomes in this species. We also describe evidence suggesting the possible re-evolution of dioicy in the genus Riccia, with probable de novo establishment of a sex chromosome from an autosome (chromosome 5), and further translocations of genes from the new sex chromosome to autosomes. Our results also indicated that micro-chromosomes are consistent genomic features, and may have evolved independently from sex chromosomes in Ricciocarpos and Riccia lineages.

Gene novelty and gene family expansion in the early evolution of Lepidoptera

BACKGROUND

Almost 10% of all known animal species belong to Lepidoptera: moths and butterflies. To understand how this incredible diversity evolved we assess the role of gene gain in driving early lepidopteran evolution. Here, we compared the complete genomes of 115 insect species, including 99 Lepidoptera, to search for novel genes coincident with the emergence of Lepidoptera.

RESULTS

We find 217 orthogroups or gene families which emerged on the branch leading to Lepidoptera; of these 177 likely arose by gene duplication followed by extensive sequence divergence, 2 are candidates for origin by horizontal gene transfer, and 38 have no known homology outside of Lepidoptera and possibly arose via de novo gene genesis. We focus on two new gene families that are conserved across all lepidopteran species and underwent extensive duplication, suggesting important roles in lepidopteran biology. One encodes a family of sugar and ion transporter molecules, potentially involved in the evolution of diverse feeding behaviours in early Lepidoptera. The second encodes a family of unusual propeller-shaped proteins that likely originated by horizontal gene transfer from Spiroplasma bacteria; we name these the Lepidoptera propellin genes.

CONCLUSION

We provide the first insights into the role of genetic novelty in the early evolution of Lepidoptera. This gives new insight into the rate of gene gain during the evolution of the order as well as providing context on the likely mechanisms of origin. We describe examples of new genes which were retained and duplicated further in all lepidopteran species, suggesting their importance in Lepidoptera evolution.

A chromosome-level genome assembly of the Peruvian Algarrobo (Neltuma pallida) provides insights on its adaptation to its unique ecological niche

The dry forests of northern Peru are dominated by the legumous tree Neltuma pallida which is adapted to hot arid and semiarid conditions in the tropics. Despite having been successfully introduced in multiple other areas around the world, N. pallida is currently threatened in its native area, where it is invaluable for the dry forest ecosystem and human subsistence. A major tool for enhancing ecosystem conservation and understanding the adaptive properties of N. pallida to dry forest ecosystems is the construction of a reference genome sequence. Here, we report on a high-quality reference genome for N. pallida. The final genome assembly size is 403.7 Mb, consisting of 14 pseudochromosomes and 63 scaffolds with an N50 size of 26.2 Mb and a 34.3% GC content. Use of Benchmarking Universal Single Copy Orthologs revealed 99.2% complete orthologs. Long terminal repeat elements dominated the repetitive sequence content which was 51.2%. Genes were annotated using N. pallida transcripts, plant protein sequences, and ab initio predictions resulting in 22,409 protein-coding genes encoding 24,607 gene models. Comparative genomic analysis showed evidence of rapidly evolving gene families related to disease resistance, transcription factors, and signaling pathways. The chromosome-scale N. pallida reference genome will be a useful resource for understanding plant evolution in extreme and highly variable environments.

Unveiling Photoperiod-Responsive Regulatory Networks in Tropical Maize Through Transcriptome Analysis

Background/Objectives: Maize (Zea mays L.), a crop of worldwide importance, owes its adaptability to diverse environments to its genetic variation. However, tropical maize exhibits intrinsic photoperiod sensitivity, limiting its adaptability to temperate regions. Photoperiod sensitivity significantly affects the flowering time and other agronomic traits, but the underlying molecular mechanisms remain poorly understood. In this study, the aim is to elucidate the transcriptional regulatory networks mediating photoperiod responses in tropical maize inbred line Su65, providing insights into improving photoperiod adaptability. Methods: RNA-seq analysis was carried out to investigate photoperiod-responsive genes and pathways in tropical line Su65 exposed to varying photoperiod conditions. Differential expression analysis, functional enrichment, and the construction of protein-protein interaction (PPI) networks were carried out to investigate transcriptional dynamics. Additionally, qRT-PCR was employed to confirm the expression patterns of key candidate genes and generate detailed temporal expression profiles. Results: A total of 1728 differentially expressed genes (DEGs) were identified, with significant enrichment in pathways such as stress responses, redox homeostasis, and secondary metabolite biosynthesis. A set of new key hub genes (such as Zm00001d048531, Zm00001d018821, Zm00001d034892, etc.) were identified through PPI network analysis. Temporal expression profiling of ZmPHYB1, ZmPHYC1, ZmFKF2, ZmGI2, and ZmPRR37a, the key genes involved in circadian rhythms, revealed distinct regulatory patterns of photoperiod-sensitive genes at different time points, highlighting their roles in flowering time regulation and developmental transitions. Conclusions: In this study, critical molecular networks underlying photoperiod sensitivity in tropical maize are uncovered and a foundation is provided for improving photoperiod adaptability through genetic improvement. By integrating RNA-seq and qRT-PCR, the research offers valuable insights into transcriptional dynamics and their role in maize development under photoperiodic regulation.

A pilot study of transcriptomic preimplantation genetic testing (PGT-T): towards a new step in embryo selection?

STUDY QUESTION

Is it possible to predict an euploid chromosomal constitution and identify a transcriptomic profile compatible with extended embryonic development from RNA sequencing (RNA-Seq) data?

SUMMARY ANSWER

It has been possible to obtain a karyotype comparable to preimplantation genetic testing for aneuploidy (PGT-A), in addition to a transcriptomic signature of embryos which might be suggestive of improved implantation capacity.

WHAT IS KNOWN ALREADY

Conventional assessment of embryo competence, based on morphology and morphokinetic, lacks knowledge of molecular aspects and faces controversy in predicting ploidy status. Understanding the embryonic transcriptome is crucial, as gene expression influences development and implantation. PGT has improved pregnancy rates, but problems persist when high-quality euploid embryos do not reach term. In fact, only around 50-60% implant, of which 10% result in miscarriage. Comprehensive approaches, including RNA-Seq, offer the potential to discover molecular markers of reproductive competence, and could theoretically be combined with extended-embryo culture platforms up to Day 14 that can be utilized as a proxy to study embryo development at post-implantation stages.

STUDY DESIGN, SIZE, DURATION

This prospective pilot cohort study was conducted from March 2023 to August 2023. A total of 30 vitrified human blastocysts with previous PGT-A diagnosis on Day 5 (D5) or Day 6 (D6) of development were analysed: n = 15 euploid and n = 15 aneuploid. Finally, 21 embryo samples were included in the study; the rest (n = 9) were excluded due to poor quality pre-sequencing data (n = 7) or highly discordant data (n = 2).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Following warming and re-expansion, embryos underwent a second trophectoderm (TE) biopsy. The embryos were then cultured until day 11 to assess their development. Biopsy analysis by RNA-Seq, studied the differential expressed genes (DEG) to compare embryos which did not or did attach to the plate: unattached embryos (n = 12) versus attached embryos (n = 9). Thus, we also obtained a specific transcriptomic signature of embryos with a "theoretical" capacity for sustained implantation, based on plate attachment on day 11.

MAIN RESULTS AND THE ROLE OF CHANCE

The digital karyotype obtained by RNA-Seq showed good concordance with the earlier PGT-A data, with a sensitivity of 0.81, a specificity of 0.83, a Cohen's Kappa of 0.66, and an area under the ROC of 0.9. At the gene level, 76 statistically significant DEGs were found in the comparison unattached versus attached embryos (Padj < 0.05; FC > 1). To address the functional implications of these differences, significantly deregulated pathways according to GO and KEGG categories were identified. The mural trophectoderm (TE) of the unattached blastocysts showed 63 significantly deregulated terms, displaying upregulation in autophagy, apoptosis, protein kinase and ubiquitin-like protein ligase activity, and downregulation of ribosome, spliceosome, kinetochore, segregation, and chromosome condensation processes. The overall transcriptomic signature specific to embryos still attached to the plate on day 11 (with a theoretically higher implantation capacity) consists of 501 genes, including: EMP2, AURKB, FOLR1, NOTCH3, LRP2, FZD5, MDH1, APOD, GPX8, COLEC12, HSPA1A, CMTM7, BEX3, which are related to implantation and embryonic development (raw P-value < 0.05; shrunk LFC > 1.1). These findings indicate that it might be possible to identify euploid embryos with a greater capacity for implantation and development, after excluding those embryos that present chromosomal alterations.

LIMITATIONS, REASONS FOR CAUTION

This study included a small sample size, remarkable variability between samples, and low success rate of RNA amplification. Also, structural chromosomal abnormalities were not included, and it was not possible to diagnose mosaic embryos. TE biopsy does not assure the chromosomal status of the whole embryo. The maximum day for in vitro development was Day 11, and attachment to the plate on this day does not provide a clear indication of implantation capacity and viability, which was not tested in this study.

WIDER IMPLICATIONS OF THE FINDINGS

The short-term goals following on from this pilot study is to expand the sample size with embryos of more complex abnormalities, and to perform a prospective in vitro preclinical validation. In a more distant future and with optimal results, this technique could have clinical application, thus increasing clinical outcomes by assessing both chromosomal content and transcriptomic profiling.

STUDY FUNDING/COMPETING INTEREST(S)

The Institut Valencià de Competitivitat Empresarial (IVACE) (IMIDCA/2022/39) and Generalitat Valenciana (CIACIF/2021/11) supported the present study. A.C. is an employee of JUNO Genetics. He has received honoraria for an IBSA lecture and a Merck lecture. He is also a minor shareholder of IVIRMA Global. The other authors have no conflicts of interest to declare.

TRIAL REGISTRATION NUMBER

N/A.

Skeletal Muscle mRNA Splicing Variants Association With Four Different Fitness and Energetic Measures in the GESTALT Study

BACKGROUND

Physical activity is essential for maintaining muscle mitochondrial function and aerobic capacity. The molecular mechanisms underlying such protective effects are incompletely understood, in part because it is difficult to separate the effects of disease status and physical activity. We explored the association of human skeletal muscle transcriptomic with four measures of energetics and mitochondria oxidative capacity in healthy individuals.

METHODS

Using RNA sequencing of vastus lateralis muscle biopsies from 82 GESTALT participants (52 males, aged 22-89 years), we explored gene and splicing variant expression profiles associated with self-reported physical activity, peak oxygen consumption (VO2 peak), muscle oxidative capacity (kPCr) and mitochondrial respiration (Mit-O2 flux). The effect of aging on gene expression was examined in participants with low and high VO2 peak.

RESULTS

The four measures of energetics were negative correlated with age and generally intercorrelated. We identified protein-coding genes associated with four energetic measures adjusting for age, muscle fiber-ratio, sex and batch effect. Mitochondrial pathways were overrepresented across all energetic variables, albeit with little overlap at the gene level. Alternative spliced transcript isoforms associated with energetics were primarily enriched for cytoplasmic ribonucleoprotein granules. The splicing pathway was up-regulated with aging in low but not in high fitness participants, and transcript isoforms detected in the low fitness group pertain to processes such as cell cycle regulation, RNA/protein localization, nuclear transport and catabolism.

CONCLUSIONS

A consistent mitochondrial signature emerged across all energetic measures. Alternative splicing was enhanced in older, low fitness participants supporting the energy-splicing axis hypothesis. The identified splicing variants were enriched in pathways involving the accumulation of ribonucleoproteins in cytoplasmic granules, whose function remains unclear. Further research is needed to understand the function of these proteoforms in promoting adaptation to low energy availability.

Comparative Genomics and Epigenomics of Transcriptional Regulation

Transcriptional regulation in response to diverse physiological cues involves complicated biological processes. Recent initiatives that leverage whole genome sequencing and annotation of regulatory elements significantly contribute to our understanding of transcriptional gene regulation. Advances in the data sets available for comparative genomics and epigenomics can identify evolutionarily constrained regulatory variants and shed light on noncoding elements that influence transcription in different tissues and developmental stages across species. Most epigenomic data, however, are generated from healthy subjects at specific developmental stages. To bridge the genotype-phenotype gap, future research should focus on generating multidimensional epigenomic data under diverse physiological conditions. Farm animal species offer advantages in terms of feasibility, cost, and experimental design for such integrative analyses in comparison to humans. Deep learning modeling and cutting-edge technologies in sequencing and functional screening and validation also provide great promise for better understanding transcriptional regulation in this dynamic field.

Dissection of the genetic basis and molecular mechanism of ovule number per ovary in oilseed rape (Brassica napus L.)

Ovule number per ovary (ONPO) determines the maximum potential of seed number per fruit that is a direct component of seed yield in crops. This study aimed to dissect the genetic basis and molecular mechanism of ONPO using a newly developed doubled haploid (DH) population in oilseed rape. In all the four investigated environments, the ONPO of 201 DH lines exhibited normal distribution with a wide variation from 22.6 to 41.8, suggesting quantitative inheritance appropriate for mapping QTL. A skeleton genetic map of 2111 markers within 19 linkage groups was developed, with a total of 1715.71 cM in length and an average of 0.82 cM between markers. Linkage mapping identified ten QTLs that were distributed on eight chromosomes and explained 7.0-15.9% of the phenotypic variance. Among these, four were identical to the reported and two were repeatedly detected with relatively large effects, highlighting their potential for marker-assisted selection. Phytohormone quantification of ovaries (at the ovule initiation stage) from two pools of high and low ONPO lines showed significant differences in the levels of nine sub-types of phytohormones, suggesting their important roles in regulating ovule number. Transcriptomic analysis identified 7689 differentially expressed genes (DEGs) between the two pools, of which nearly half were enriched into functional categories of reported genes regulating ONPO, including protein, RNA, signalling, miscellaneous, development, hormone metabolism, and tetrapyrrole synthesis. Integration of linkage QTL mapping, transcriptome sequencing and BLAST analysis identified 15 homologues of reported ovule number genes and 327 DEGs in the QTL regions, which were considered as direct and potential candidate genes. These findings propose further insights into the genetic basis and molecular mechanisms of ONPO, which will facilitate future gene cloning and genetic improvement for enhancing seed yield in oilseed rape.

Near telomere-to-telomere genome assemblies of Silkie Gallus gallus and Mallard Anas platyrhynchos restored the structure of chromosomes and "missing" genes in birds

BACKGROUND

Chickens and ducks are vital sources of animal protein for humans. Recent pangenome studies suggest that a single genome is insufficient to represent the genetic information of a species, highlighting the need for more comprehensive genomes. The bird genome has more than tens of microchromosomes, but comparative genomics, annotations, and the discovery of variations are hindered by inadequate telomere-to-telomere level assemblies. We aim to complete the chicken and duck genomes, recover missing genes, and reveal common and unique chromosomal features between birds.

RESULTS

The near telomere-to-telomere genomes of Silkie Gallus gallus and Mallard Anas platyrhynchos were successfully assembled via multiple high-coverage complementary technologies, with quality values of 36.65 and 44.17 for Silkie and Mallard, respectively; and BUSCO scores of 96.55% and 96.97% for Silkie and Mallard, respectively; the mapping rates reached over 99.52% for both assembled genomes, these evaluation results ensured high completeness and accuracy. We successfully annotated 20,253 and 19,621 protein-coding genes for Silkie and Mallard, respectively, and assembled gap-free sex chromosomes in Mallard for the first time. Comparative analysis revealed that microchromosomes differ from macrochromosomes in terms of GC content, repetitive sequence abundance, gene density, and levels of 5mC methylation. Different types of arrangements of centromeric repeat sequence centromeres exist in both Silkie and the Mallard genomes, with Mallard centromeres being invaded by CR1. The highly heterochromatic W chromosome, which serves as a refuge for ERVs, contains disproportionately long ERVs. Both Silkie and the Mallard genomes presented relatively high 5mC methylation levels on sex chromosomes and microchromosomes, and the telomeres and centromeres presented significantly higher 5mC methylation levels than the whole genome. Finally, we recovered 325 missing genes via our new genomes and annotated TNFA in Mallard for the first time, revealing conserved protein structures and tissue-specific expression.

CONCLUSIONS

The near telomere-to-telomere assemblies in Mallard and Silkie, with the first gap-free sex chromosomes in ducks, significantly enhanced our understanding of genetic structures in birds, specifically highlighting the distinctive chromosome features between the chicken and duck genomes. This foundational work also provides a series of newly identified missing genes for further investigation.

A Chromosome-Scale Genome of Trametes versicolor and Transcriptome-Based Screening for Light-Induced Genes That Promote Triterpene Biosynthesis

Trametes versicolor is an important fungus with medicinal properties and a significant role in lignocellulose degradation. In this study, we constructed a high-quality chromosome-level genome of T. versicolor using Illumina, PacBio HiFi, and Hi-C sequencing technologies. The assembled genome is 47.42 Mb in size and contains 13,307 protein-coding genes. BUSCO analysis revealed genome and gene completeness results of 95.80% and 95.90%, respectively. Phylogenetic analysis showed that T. versicolor is most closely related to T. pubescens, followed by T. cinnabarina and T. coccinea. Comparative genomic analysis identified 266 syntenic blocks between T. versicolor and Wolfiporia cocos, indicating a conserved evolutionary pattern between the two species. Gene family analysis highlighted the expansion and contraction of genes in functional categories related to the biosynthesis of secondary metabolites, including several T. versicolor-specific genes. Key genes involved in lignocellulose degradation and triterpene production were identified within the CAZyme and CYP450 gene families. Transcriptomic analysis under dark and light conditions revealed significant changes in the expression of genes related to secondary metabolism, suggesting that light signals regulate metabolic pathways. A total of 2577 transporter proteins and 2582 membrane proteins were identified and mapped in the T. versicolor genome, and 33 secondary metabolite gene clusters were identified, including two light-sensitive triterpene biosynthesis clusters. This study offers a comprehensive genomic resource for further investigation into the functional genomics, metabolic regulation, and triterpene biosynthesis of T. versicolor, providing valuable insights into fungal evolution and biotechnological applications.

Environmental community transcriptomics: strategies and struggles

Transcriptomics is the study of RNA transcripts, the portion of the genome that is transcribed, in a specific cell, tissue, or organism. Transcriptomics provides insight into gene expression patterns, regulation, and the underlying mechanisms of cellular processes. Community transcriptomics takes this a step further by studying the RNA transcripts from environmental assemblies of organisms, with the intention of better understanding the interactions between members of the community. Community transcriptomics requires successful extraction of RNA from a diverse set of organisms and subsequent analysis via mapping those reads to a reference genome or de novo assembly of the reads. Both, extraction protocols and the analysis steps can pose hurdles for community transcriptomics. This review covers advances in transcriptomic techniques and assesses the viability of applying them to community transcriptomics.

Comprehensive genome annotation of the model ciliate Tetrahymena thermophila by in-depth epigenetic and transcriptomic profiling

The ciliate Tetrahymena thermophila is a well-established unicellular model eukaryote, contributing significantly to foundational biological discoveries. Despite its acknowledged importance, current studies on Tetrahymena biology face challenges due to gene annotation inaccuracy, particularly the notable absence of untranslated regions (UTRs). To comprehensively annotate the Tetrahymena macronuclear genome, we collected extensive transcriptomic data spanning various cell stages. To ascertain transcript orientation and transcription start/end sites, we incorporated data on epigenetic marks displaying enrichment towards the 5' end of gene bodies, including H3 lysine 4 tri-methylation (H3K4me3), histone variant H2A.Z, nucleosome positioning and N6-methyldeoxyadenine (6mA). Cap-seq data was subsequently applied to validate the accuracy of identified transcription start sites. Additionally, we integrated Nanopore direct RNA sequencing (DRS), strand-specific RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) data. Using a newly developed bioinformatic pipeline, coupled with manual curation and experimental validation, our work yielded substantial improvements to the current gene models, including the addition of 2,481 new genes, updates to 23,936 existing genes, and the incorporation of 8,339 alternatively spliced isoforms. Furthermore, novel UTR information was annotated for 26,687 high-confidence genes. Intriguingly, 20% of protein-coding genes were identified to have natural antisense transcripts characterized by high diversity in alternative splicing, thus offering insights into understanding transcriptional regulation. Our work will enhance the utility of Tetrahymena as a robust genetic toolkit for advancing biological research, and provides a promising framework for genome annotation in other eukaryotes.

RNF138 contributes to cisplatin resistance in nasopharyngeal carcinoma cells

Resistance to chemotherapy is a significant concern in the treatment of nasopharyngeal carcinoma (NPC), and occurs due to various mechanisms. This study is aimed to evaluate the effects of RING finger protein 138 (RNF138) in the development of cisplatin resistance to NPC. After gene overexpression and silencing, the expression levels of RNF138 were evaluated. The impacts of RNF138 on the proliferation and apoptosis rate of NPC cells were then assessed. γ-H2AX-mediated DNA damage was determined via immunofluorescence assay. Moreover, a tumor xenograft mouse model was developed to investigate the role of RNF138 on NPC progression in vivo. Additionally, transcriptome analysis was performed in 5-8 F cells transfection with OE-RNF1138 or OE-NC.Cisplatin significantly inhibited the proliferation, and promoted apoptosis and DNA damage in NPC cells; however, overexpression of RNF138 reversed the aforementioned regulatory role of cisplatin on NPC cells. Knockdown of RNF138 resulted in contrasting phenotypic outcomes. Additionally, in nude mice, RNF138 overexpression attenuated the suppressive effects of cisplatin on the growth of xenograft tumor, while RNF138 silencing further enhanced the inhibiting role of cisplatin. We further indicated that in 5-8 F cells following RNF138 overexpression, some pathways such as PI3K-Akt signaling pathway, human papillomavirus infection and ErbB signaling pathway that have been reported to be associated with NPC progression and cisplatin resistance were significantly enriched. These findings indicate that the modulation of RNF138 could potentially address the issue of chemotherapy failure by overcoming cisplatin resistance in NPC cells, making it a promising candidate for targeted drug therapy.

Uneven distribution of prokaryote-derived horizontal gene transfer in fungi: a lifestyle-dependent phenomenon

Horizontal gene transfer (HGT) in fungi is less understood despite its significance in prokaryotes. In this study, we systematically searched for HGT events in 829 representative fungal genomes. Using a combination of sequence similarity and phylogeny-based approaches, we detected 20,093 prokaryotic-derived transferred genes across 750 fungal genomes, via 8,815 distinct HGT events. Notably, our analysis revealed that eight lifestyle-related traits significantly influence HGT diversity in fungi. For instance, parasites exhibited the highest estimated number of HGT-acquired genes, followed by saprotrophs, with symbionts showing the lowest. HGT-acquired genes were predominantly associated with metabolism and cellular functions and underwent purifying selection. Moreover, horizontally transferred genes with introns have significantly higher expression levels compared to intron-lacking genes, suggesting a probable role of intron gains in the adaptation of HGT-acquired genes. Overall, our findings highlight the influence of lifestyle on HGT diversity in fungi and underscore the substantial contribution of HGT to fungal adaptation.

IMPORTANCE

This study sheds new light on the role of horizontal gene transfer (HGT) in fungi, an area that has remained relatively unexplored compared to its well-established prevalence in bacteria. By analyzing 829 fungal genomes, we identified over 20,000 genes that fungi acquired from prokaryotes, revealing the significant impact of HGT on fungal evolution. Our findings highlight that fungal lifestyle traits, such as being parasitic or saprotrophic, play a key role in determining the extent of HGT, with parasites showing the highest gene acquisition rates. We also uncovered unique patterns of HGT occurrence based on fungal morphology and reproduction. Importantly, genes with introns, which are more highly expressed, appear to play a crucial role in fungal adaptation. This research deepens our understanding of how HGT contributes to the metabolic diversity and ecological success of fungi, and it underscores the broader significance of gene transfer in shaping fungal evolution.

Genomic evidence for hybridization and introgression between blue peafowl and endangered green peafowl and molecular foundation of leucistic plumage of blue peafowl

INTRODUCTION

The blue peafowl (Pavo cristatus) and the green peafowl (Pavo muticus) have garnered significant public affection due to their stunning appearance, although the green peafowl is currently endangered. The causative mutation that causes the leucistic plumage of the blue peafowl (also called white peafowl) remains unknown.

RESULTS

In this study, we generated a chromosome-level reference genome of the blue peafowl with a contig N50 of 30.6 Mb, including the autosomes, Z and W sex chromosomes, and a complete mitochondria DNA sequence. Data from 77 peafowl whole genomes, 76 peafowl mitochondrial genomes, and 33 peafowl W chromosomes genomes provided the first substantial genetic evidence for recent hybridization between green peafowls and blue peafowls. We found 3 hybrid green peafowls in zoo samples rather than in the wild samples, with a blue peafowl genomic content of 16-34%. Maternal genetic analysis showed that 2 of the hybrid female green peafowls contained complete blue peafowl mitochondrial genomes and W chromosomes. Some animal protection agencies release captive green peafowls in order to maintain the wild population of green peafowls. Therefore, to better protect the endangered green peafowl, we suggest that purebred identification must be carried out before releasing green peafowls from zoos into the wild in order to prevent the hybrid green peafowl from contaminating the wild green peafowl. In addition, we also found that there were historical introgression events of green peafowl to blue peafowl in 4 zoo blue peafowl individuals. The introgressed genomic regions contain IGFBP1 and IGFBP3 genes that could affect blue peafowl body size. Finally, we identified that the nonsense mutation (g.4:12583552G>A) in the EDNRB2 gene is the genetic causative mutation for leucistic plumage of blue peafowl, preventing melanocytes from being transported into plumage, thereby inhibiting melanin deposition.

CONCLUSION

Our research provides both theoretical and empirical support for the conservation of the endangered green peafowl. The high-quality genome and genomic data also provide a valuable resource for blue peafowl genomics-assisted breeding.

Genomic insights into ecological adaptation of oaks revealed by phylogenomic analysis of multiple species

Understanding the ecological adaptation of tree species can not only reveal the evolutionary potential but also benefit biodiversity conservation under global climate change. Quercus is a keystone genus in Northern Hemisphere forests, and its wide distribution in diverse ecosystems and long evolutionary history make it an ideal model for studying the genomic basis of ecological adaptations. Here we used a newly sequenced genome of Quercus gilva, an evergreen oak species from East Asia, with 18 published Fagales genomes to determine how Fagaceae genomes have evolved, identify genomic footprints of ecological adaptability in oaks in general, as well as between evergreen and deciduous oaks. We found that oak species exhibited a higher degree of genomic conservation and stability, as indicated by the absence of large-scale chromosomal structural variations or additional whole-genome duplication events. In addition, we identified expansion and tandem repetitions within gene families that contribute to plant physical and chemical defense (e.g., cuticle biosynthesis and oxidosqualene cyclase genes), which may represent the foundation for the ecological adaptation of oak species. Circadian rhythm and hormone-related genes may regulate the habits of evergreen and deciduous oaks. This study provides a comprehensive perspective on the ecological adaptations of tree species based on phylogenetic, genome evolutionary, and functional genomic analyses.

Purification of micrococcal nuclease for use in ribosomal profiling of high-salinity extremophiles

Nucleases, that is, enzymes that catalyze the hydrolysis of phosphodiester bonds in nucleic acids, are essential tools in molecular biology and biotechnology. Staphylococcus aureus nuclease is particularly interesting due to its thermostability and Ca2+ dependence, making it the prime choice for applications where nuclease modulation is critical, such as ribosome profiling in bacteria and halophilic archaea. The latter poses a technical and economical challenge: high salt reaction conditions are essential for maintaining ribosome integrity but negatively impact the micrococcal nuclease (MNase) activity, necessitating using large amounts of nuclease to achieve efficient cleavage. Here, we set out to generate an optimized production protocol for two forms of MNase-fully processed MNaseA and the 19 amino acid propeptide-containing MNaseB-and to biochemically benchmark them against a commercial nuclease. Our results show that both MNases are highly active in normal reaction conditions, but MNaseA maintains higher enzymatic activity in high salt concentrations than MNaseB. MNaseA also retains >90% of its activity after multiple freeze-thaw cycles when stored at -80 °C in a buffer containing 5% glycerol. Importantly, ribosome profiling experiments in the haloarchaeon Haloferax volcanii demonstrated that MNaseA produces ribosome footprints and hallmarks of active translation highly comparable to those obtained with the commercial nuclease, making it a suitable alternative for high-salt ribosome profiling applications. In conclusion, our method can be easily implemented for efficient MNaseA production, thereby providing access to an effective, robust, and cost-efficient alternative to commercial nucleases, as well as facilitating future translation studies into halophilic organisms.

The Isodon serra genome sheds light on tanshinone biosynthesis and reveals the recursive karyotype evolutionary histories within Lamiales

Lamiales is one of the largest orders of angiosperms with a complex evolutionary history and plays a significant role in human life. However, the polyploidization and chromosome evolution histories within this group remain in mystery. Among Lamiales, Isodon serra (Maxim.) Kudô shines for its abundance of diterpenes, notably tanshinones, long used in East Asia to combat toxicity and inflammation. Yet, the genes driving its biosynthesis and the factors governing its regulation linger in obscurity. Here, we present the telomere-to-telomere genome assembly of I. serra and, through gene-to-metabolite network analyses, pinpoint the pivotal tanshinone biosynthesis genes and their co-expressed transcription factors. Particularly, through luciferase (LUC) assays, we speculate that IsMYB-13 and IsbHLH-8 may upregulate IsCYP76AH101, which is the key step in the biosynthesis of the tanshinone precursor. Among Lamiales, Oleaceae, Gesneriaceae and Plantaginaceae successively sister to a clade of seven Lamiales families, all sharing a recent whole-genome duplication (designated as α event). By reconstructing the ancestral Lamiales karyotypes (ALK) and post-α event (ALKα), we trace chromosomal evolution trajectories across Lamiales species. Notably, one chromosomal fusion is detected from ALK to ALKα, and three shared chromosomal fusion events are detected sequentially from ALKα to I. serra, which fully supports the phylogeny constructed using single-copy genes. This comprehensive study illuminates the genome evolution and chromosomal dynamics of Lamiales, further enhancing our understanding of the biosynthetic mechanisms underlying the medicinal properties of I. serra.

Telomere-to-telomere sheep genome assembly identifies variants associated with wool fineness

Ongoing efforts to improve sheep reference genome assemblies still leave many gaps and incomplete regions, resulting in a few common failures and errors in genomic studies. Here, we report a 2.85-Gb gap-free telomere-to-telomere genome of a ram (T2T-sheep1.0), including all autosomes and the X and Y chromosomes. This genome adds 220.05 Mb of previously unresolved regions and 754 new genes to the most updated reference assembly ARS-UI_Ramb_v3.0; it contains four types of repeat units (SatI, SatII, SatIII and CenY) in centromeric regions. T2T-sheep1.0 has a base accuracy of more than 99.999%, corrects several structural errors in previous reference assemblies and improves structural variant detection in repetitive sequences. Alignment of whole-genome short-read sequences of global domestic and wild sheep against T2T-sheep1.0 identifies 2,664,979 new single-nucleotide polymorphisms in previously unresolved regions, which improves the population genetic analyses and detection of selective signals for domestication (for example, ABCC4) and wool fineness (for example, FOXQ1).

Near telomere-to-telomere assembly of the Tarim pigeon (Columba livia) genome

Pigeons serve as important model animals and commercial poultry. The Tarim pigeon, as a breed of Columba livia, is a locally indigenous breed unique to China. While the genome of C. livia was published in 2013, its assembly was fragmented and incomplete. In this study, we generated a near telomere-to-telomere assembly of the pigeon genome using the sequencing platform of PacBio HiFi, Nanopore long reads and Hi-C. The assembled genome spans 1295.8 Mb, with a contig N50 size of 49 Mb and a scaffold N50 size of 85 Mb. Approximately 98.4% of the assembly is anchored onto 41 chromosomes, with a BUSCO assessment indicating a completeness of 97.2%. And telomeres were identified at both ends of the four chromosomes. A total of 21,450 genes were annotated. The genome assembly of C. livia lays the foundation for understanding their genetic composition and evolutionary history and contributes to the pigeon breeding industry. Additionally, it will provide a basis for further management and conservation of pigeon breed diversity.

Transposable element exonization generates a reservoir of evolving and functional protein isoforms

Alternative splicing enhances protein diversity in different ways, including through exonization of transposable elements (TEs). Recent transcriptomic analyses identified thousands of unannotated spliced transcripts with exonizing TEs, but their contribution to the proteome and biological relevance remains unclear. Here, we use transcriptome assembly, ribosome profiling, and proteomics to describe a population of 1,227 unannotated TE exonizing isoforms generated by mRNA splicing and recurrent in human populations. Despite being shorter and lowly expressed, these isoforms are shared between individuals and efficiently translated. Functional analyses show stable expression, specific cellular localization, and, in some cases, modified functions. Exonized TEs are rich in ancient genes, whereas the involved splice sites are recent and can be evolutionarily conserved. In addition, exonized TEs contribute to the secondary structure of the emerging isoforms, supporting their functional relevance. We conclude that TE-spliced isoforms represent a diversity reservoir of functional proteins on which natural selection can act.

The CCL5/CCR5/SHP2 axis sustains Stat1 phosphorylation and activates NF-κB signaling promoting M1 macrophage polarization and exacerbating chronic prostatic inflammation

BACKGROUND AND OBJECTIVE

Chronic prostatitis (CP) is a condition markered by persistent prostate inflammation, yet the specific cytokines driving its progression remain largely undefined. This study aims to identify key cytokines involved in CP and investigate their role in driving inflammatory responses through mechanistic and therapeutic exploration.

METHODS

A 48-cytokine panel test was conducted to compare the plasma cytokine profiles between participants with CP-like symptoms (CP-LS) and healthy controls. Experimental autoimmune prostatitis (EAP) models were used for functional validation, with further mechanistic studies performed through in vivo and in vitro assays. Pharmacological inhibition was applied using maraviroc, and pathway inhibitors to assess therapeutic potential.

RESULTS

Our analysis identified CCL5 as one of the most prominently elevated cytokines in CP-LS patients. Further validation in the EAP model mice confirmed elevated CCL5 levels, highlighting its role in driving prostatic inflammation. Mechanistic studies revealed that CCL5 interacts with the CCR5 receptor, promoting M1 macrophage polarization and activating key inflammatory signaling pathways, including Stat1 and NF-κB, as indicated by increased phosphorylation of Stat1 and p65. In vitro, CCL5 combined with LPS stimulation amplified these effects, further promoting M1 polarization. CCL5 also sustained Stat1 activation by inhibiting its dephosphorylation through reduced interaction with SHP2, leading to prolonged inflammatory signaling. Single-cell transcriptomics confirmed high CCR5 expression in macrophages, correlating with inflammatory pathways. Pharmacological inhibition of CCR5, or its downstream signaling, significantly reduced macrophage-driven inflammation both in vivo and in vitro.

CONCLUSION

These findings establish the CCL5/CCR5 axis as a critical driver of persistant prostatic inflammation and present it as a potential therapeutic target for CP.

Human introns contain conserved tissue-specific cryptic poison exons

Eukaryotic cells express a large number of transcripts from a single gene due to alternative splicing. Despite hundreds of thousands of splice isoforms being annotated in databases, it has been reported that the current exon catalogs remain incomplete. At the same time, introns of human protein-coding (PC) genes contain a large number of evolutionarily conserved elements with unknown function. Here, we explore the possibility that some of them represent cryptic exons that are expressed in rare conditions. We identified a group of cryptic exons that are similar to the annotated exons in terms of evolutionary conservation and RNA-seq read coverage in the Genotype-Tissue Expression dataset. Most of them were poison, i.e. generated an nonsense-mediated decay (NMD) isoform upon inclusion, and many showed signs of tissue-specific and cancer-specific expression and regulation. We performed RNA-seq in A549 cell line treated with cycloheximide to inactivate NMD and confirmed using quantitative polymerase chain reaction that seven of eight exons tested are, indeed, expressed. This study shows that introns of human PC genes contain cryptic poison exons, which reside in conserved intronic regions and remain not fully annotated due to insufficient representation in RNA-seq libraries.

A novel small RNA PhaS contributes to polymyxin B-heteroresistance in carbapenem-resistant Klebsiella pneumoniae

In recent years, polymyxin has been used as a last-resort therapy for carbapenem-resistant bacterial infections. The emergence of heteroresistance (HR) to polymyxin hampers the efficacy of polymyxin treatment by amplifying resistant subpopulation. However, the mechanisms behind polymyxin HR remain unclear. Small noncoding RNAs (sRNAs) play an important role in regulating drug resistance. The purpose of this study was to investigate the effects and mechanisms of sRNA on polymyxin B (PB)-HR in carbapenem-resistant Klebsiella pneumoniae. In this study, a novel sRNA PhaS was identified by transcriptome sequencing. PhaS expression was elevated in the PB heteroresistant subpopulation. Overexpression and deletion of PhaS were constructed in three carbapenem-resistant K. pneumoniae strains. Population analysis profiling, growth curve, and time-killing curve analysis showed that PhaS enhanced PB-HR. In addition, we verified that PhaS directly targeted phoP through the green fluorescent protein reporter system. PhaS promoted the expression of phoP, thereby encouraging the expression of downstream genes pmrD and arnT. This upregulation of arnT promoted the 4-amino-4-deoxyL-arabinosaccharide (L-Ara4N) modification of lipid A in PhaS overexpressing strains, thus enhancing PB-HR. Further, within the promoter region of PhaS, specific PhoP recognition sites were identified. ONPG assays and RT-qPCR analysis confirmed that PhaS expression was positively modulated by PhoP and thus up-regulated by PB stimulation. To sum up, a novel sRNA enhancing PB-HR was identified and a positive feedback regulatory pathway of sRNA-PhoP/Q was demonstrated in the study. This helps to provide a more comprehensive and clear understanding of the underlying mechanisms behind polymyxin HR in carbapenem-resistant K. pneumoniae.

Diurnal control of H3K27me1 deposition shapes expression of a subset of cell cycle and DNA damage response genes

Rhythmic oscillation of biological processes helps organisms adapt their physiological responses to the most appropriate time of the day. Chromatin remodeling has been described as one of the molecular mechanisms controlling these oscillations. The importance of these changes in transcriptional activation as well as in the maintenance of heterochromatic regions has been widely demonstrated. However, little is still known on how diurnal changes can impact the global status of chromatin modifications and, hence, control gene expression. In plants, the repressive mark H3K27me1, deposited by ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 and 6 (ATXR5 and 6) methyltransferases, is largely associated with transposable elements but also covers lowly expressed genes. Here we show that this histone modification is preferentially deposited during the night. In euchromatic regions, it is found along the bodies of DNA damage response genes (DDR), where it is needed for their proper expression. The absence of H3K27me1 translates into an enhanced expression of DDR genes that follows a rhythmic oscillation pattern. This evidences a link between chromatin modifications and their synchronization with the diurnal cycle in order to accurately modulate the activation of biological processes to the most appropriate time of the day.

Analysis of genetic structure and identification of important genes associated with muscle growth in Fujian Muscovy duck

Fujian Muscovy duck is a well-known meat waterfowl in Fujian Province due to its high meat production, superior breeding potential, and strong resistance. To fully explore the genetic characteristics of these advantages, Fujian black Muscovy duck and white Muscovy duck were used for whole-genome re-sequencing and transcriptome analyses. Population structure analysis showed significant differentiation between the two feather strains. Runs of homozygosity analysis indicated a stronger artificial influence on the black-feathered strain, with ROH island genes notably enriched in muscle tissue-related terms and pathways. Selective sweep and transcriptome analysis revealed a significant enrichment of genes linked to muscle tissue and muscle fiber-related terms and pathways. Key candidate genes identified, such as MEF2C, MYOZ2, and METTL21C, are believed to play crucial roles in meat production in Fujian Muscovy duck. This study offers a new perspective on improving meat production in Fujian Muscovy duck, which can benefit breeding strategies and production management.

Bach2 repression of CD36 regulates lipid-metabolism-linked effector functions in follicular B cells

The transcription repressor Bach2 plays a crucial role in shaping humoral immunity, but its cell-autonomous function remains elusive. Here, we reveal the mechanism by which Bach2 regulates effector cell maturation in peripheral B cells. In response to Toll-like receptor (TLR) agonists, Bach2 deficiency promotes the differentiation of follicular, but not marginal zone, B cells into effector cells, producing interleukin (IL)-6 and antibodies. This phenomenon is associated with changes in lipid metabolism, such as increases in CD36 expression, lipid influx, and fatty acid oxidation. Consistent with this, Bach2-deficient B cells exhibit elevated levels of mitochondrial oxidative stress, lipid peroxidation, and p38 activation. Mechanistically, Bach2 acts as a repressor of Cd36, and inhibition of CD36 or fatty acid oxidation reduces the differentiation of naive B cells into IL-6- and antibody-secreting cells. These results indicate Bach2 as a key metabolic checkpoint regulator crucial for maintaining a functionally quiescent state of follicular B cells.

Construction of a Dataset for All Expressed Transcripts for Alzheimer's Disease Research

Accurate identification and functional annotation of splicing isoforms and non-coding RNAs (lncRNAs), alongside full-length protein-encoding transcripts, are critical for understanding gene (mis)regulation and metabolic reprogramming in Alzheimer's disease (AD). This study aims to provide a comprehensive and accurate transcriptome resource to improve existing AD transcript databases. Background/Objectives: Gene mis-regulation and metabolic reprogramming play a key role in AD, yet existing transcript databases lack accurate and comprehensive identification of splicing isoforms and lncRNAs. This study aims to generate a refined transcriptome dataset, expanding the understanding of AD onset and progression. Methods: Publicly available RNA-seq data from pre-AD and AD tissues were utilized. Advanced bioinformatics tools were applied to assemble and annotate full-length transcripts, including splicing isoforms and lncRNAs, with an emphasis on correcting errors and enhancing annotation accuracy. Results: A significantly improved transcriptome dataset was generated, which includes detailed annotations of splicing isoforms and lncRNAs. This dataset expands the scope of existing AD transcript databases and provides new insights into the molecular mechanisms underlying AD. The findings demonstrate that the refined dataset captures more relevant details about AD progression compared to publicly available data. Conclusions: The newly developed transcriptome resource and the associated analysis tools offer a valuable contribution to AD research, providing deeper insights into the disease's molecular mechanisms. This work supports future research into gene regulation and metabolic reprogramming in AD and serves as a foundation for exploring novel therapeutic targets.

Identification of susceptibility modules and genes for peri-implantitis compared to periodontitis within the same host environment using weighted gene co-expression network analysis

PURPOSE

This study aimed to identify new susceptibility modules and genes by analyzing the transcriptional profiles of peri-implantitis and periodontitis within the same host environment, using weighted gene co-expression network analysis (WGCNA).

METHODS

Gingival tissue samples were collected from 10 patients, each presenting with both periodontitis and peri-implantitis sites, and were used for RNA sequencing. We conducted WGCNA to identify key modules that showed distinct transcriptional expression profiles between periodontitis and peri-implantitis. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses were carried out using R software. Genes with an adjusted P value greater than 0.05 were excluded from gene selection using the Pearson correlation method.

RESULTS

A total of 2,226 regulated genes were identified, and those with similar expression patterns were grouped into 5 color-coded functional modules using WGCNA. Among these, 3 modules showed distinct differences in expression profiles between peri-implantitis and periodontitis. The turquoise and yellow modules were associated with upregulation in peri-implantitis, while the blue module was linked to periodontitis. This finding suggests that peri-implantitis and periodontitis have significantly different transcriptional signatures. Over-representation analysis was conducted to explore the component genes of the established modules. The top-ranked genes, selected based on their network connectivity within the modules, were identified using DESeq2 and were considered hub genes.

CONCLUSIONS

WGCNA revealed distinct modular gene patterns in peri-implantitis and periodontitis, highlighting transcriptional differences between the 2 conditions. Notably, we identified 10 key genes from each of the 3 modules-the blue module associated with periodontitis-dominant pathways, and the turquoise and yellow modules associated with peri-implantitis-dominant pathways. The hub genes and pathways unveiled in this research are likely key contributors to the progression of peri-implantitis and warrant further exploration as promising candidates.

Haplotype-resolved genome and population genomics of the threatened garden dormouse in Europe

Genomic resources are important for evaluating genetic diversity and supporting conservation efforts. The garden dormouse (Eliomys quercinus) is a small rodent that has experienced one of the most severe modern population declines in Europe. We present a high-quality haplotype-resolved reference genome for the garden dormouse, and combine comprehensive short and long-read transcriptomics data sets with homology-based methods to generate a highly complete gene annotation. Demographic history analysis of the genome reveal a sharp population decline since the last interglacial, indicating an association between colder climates and population declines before anthropogenic influence. Using our genome and genetic data from 100 individuals, largely sampled in a citizen-science project across the contemporary range, we conduct the first population genomic analysis for this species. We find clear evidence for population structure across the species' core Central European range. Notably, our data show that the Alpine population, characterized by strong differentiation and reduced genetic diversity, is reproductively isolated from other regions and likely represents a differentiated evolutionary significant unit (ESU). The predominantly declining Eastern European populations also show signs of recent isolation, a pattern consistent with a range expansion from Western to Eastern Europe during the Holocene, leaving relict populations now facing local extinction. Overall, our findings suggest that garden dormouse conservation may be enhanced in Europe through the designation of ESUs.

Telomere-to-telomere genome assembly of a male goat reveals variants associated with cashmere traits

A complete goat (Capra hircus) reference genome enhances analyses of genetic variation, thus providing insights into domestication and selection in goats and related species. Here, we assemble a telomere-to-telomere (T2T) gap-free genome (2.86 Gb) from a cashmere goat (T2T-goat1.0), including a Y chromosome of 20.96 Mb. With a base accuracy of >99.999%, T2T-goat1.0 corrects numerous genome-wide structural and base errors in previous assemblies and adds 288.5 Mb of previously unresolved regions and 446 newly assembled genes to the reference genome. We sequence the genomes of five representative goat breeds for PacBio reads, and use T2T-goat1.0 as a reference to identify a total of 63,417 structural variations (SVs) with up to 4711 (7.42%) in the previously unresolved regions. T2T-goat1.0 was applied in population analyses of global wild and domestic goats, which revealed 32,419 SVs and 25,397,794 SNPs, including 870 SVs and 545,026 SNPs in the previously unresolved regions. Also, our analyses reveal a set of selective variants and genes associated with domestication (e.g., NKG2D and ABCC4) and cashmere traits (e.g., ABCC4 and ASIP).

Nascent transcript O-MAP reveals the molecular architecture of a single-locus subnuclear compartment built by RBM20 and the TTN RNA

Eukaryotic nuclei adopt a highly compartmentalized architecture that influences nearly all genomic processes. Understanding how this architecture impacts gene expression has been hindered by a lack of tools for elucidating the molecular interactions at individual genomic loci. Here, we adapt oligonucleotide-mediated proximity-interactome mapping (O-MAP) to biochemically characterize discrete, micron-scale nuclear neighborhoods. By targeting O-MAP to introns within the TTN pre-mRNA, we systematically map the chromatin loci, RNAs, and proteins within a muscle-specific RNA factory organized around the TTN locus. This reveals an unanticipated compartmental architecture that organizes cis - and trans -interacting chromosomal domains, including a hub of transcriptionally silenced chromatin. The factory also recruits dozens of unique RNA-binding and chromatin-scaffolding factors, including QKI and SAFB, along with their target transcripts. Loss of the cardiac-specific splicing factor RBM20-a master regulator of TTN splicing that is mutated in dilated cardiomyopathy-remodels nearly every facet of this architecture. This establishes O-MAP as a pioneering method for probing single-locus, microcompartment-level interactions that are opaque to conventional tools. Our findings suggest new mechanisms by which coding genes can "moonlight" in nuclear-architectural roles.

Nanopore-Based Sequencing of the Full-Length Transcriptome of Male and Female Cleavage-Stage Embryos of the Chinese Mitten Crab (Eriocheir sinensis)

The cleavage stage plays a crucial role in embryo development, characterized by a swift surge in cell proliferation alongside the accurate genetic material transmission to offspring. To delve into the characteristics of sex development during the cleavage stage of embryos, we generated the full-length transcriptome of Eriocheir sinensis male and female cleavage-stage embryos using Oxford Nanopore Technologies (ONT). Notably, this investigation represents the first sequencing effort distinguishing between genders in E. sinensis embryos. In the transcriptome structure analysis, male and female cleavage-stage embryos, while not clustered, exhibited a comparable frequency of alternative splicing (AS) occurrences. We also successfully identified 2875 transcription factors (TFs). The quantitative analysis showed the top 150 genes, in which the highly expressed genes in male embryos predominantly related to protein synthesis and metabolism. Further investigation unveiled 500 differentially expressed genes (DEGs), of which 7 male-biased ribosomal protein genes (RPGs) were particularly noteworthy and further confirmed. These analyses suggest that there may be a more active protein synthesis process in male E. sinensis cleavage-stage embryos. Furthermore, among the 2875 identified TFs, we predicted that 18 TFs could regulate the differentially expressed RPGs, with most TFs belonging to the zf-C2H2 and Homeobox families, which are crucial for embryonic development. During the cleavage stage of E. sinensis, the differential RPGs between genders were intricately linked to energy metabolism. We proposed that these RPGs exert regulatory effects on gene expression in E. sinensis, thereby regulating the difference of development between male and females. Our research sheds light on the developmental mechanisms of E. sinensis during the embryo stage and establishes a groundwork for a deeper understanding of sex development in E. sinensis. The results also provide comprehensive full-length transcriptome data for future gene expression and genetic studies in E. sinensis.

Transcriptomic responses of Solanum tuberosum cv. Pirol to arbuscular mycorrhiza and potato virus Y (PVY) infection

Arbuscular mycorrhizal fungi (AMF) serve as both plant symbionts and allies in resisting pathogens and environmental stresses. Mycorrhizal colonization of plant roots can influence the outcomes of plant-pathogen interactions by enhancing specific host defense mechanisms. The transcriptional responses induced by AMF in virus-infected plants remain largely unexplored. In the presented study, we employed a comprehensive transcriptomic approach and qPCR to investigate the molecular determinants underlying the interaction between AMF and potato virus Y (PVY) in Solanum tuberosum L. Our primary goal was to identify the symbiosis- and defense-related determinants activated in mycorrhizal potatoes facing PVY. Through a comparative analysis of mRNA transcriptomes in experimental treatments comprising healthy and PVY-infected potatoes colonized by two AMF species, Rhizophagus regularis or Funneliformis mosseae, we unveiled the overexpression of genes associated with mycorrhiza, including nutrient exchange, lipid transfer, and cell wall remodeling. Furthermore, we identified several differentially expressed genes upregulated in all mycorrhizal treatments that encoded pathogenesis-related proteins involved in plant immune responses, thus verifying the bioprotective role of AMF. We investigated the relationship between mycorrhiza levels and PVY levels in potato leaves and roots. We found accumulation of the virus in the leaves of mycorrhizal plants, but our studies additionally showed a reduced PVY content in potato roots colonized by AMF, which has not been previously demonstrated. Furthermore, we observed that a virus-dependent reduction in nutrient exchange could occur in mycorrhizal roots in the presence of PVY. These findings provide an insights into the interplay between virus and AMF.

Integration of mRNA and miRNA Analysis Sheds New Light on the Muscle Response to Heat Stress in Spotted Sea Bass (Lateolabrax maculatus)

Temperature is a crucial environmental factor for fish. Elevated temperatures trigger various physiological and molecular responses designed to maintain internal environmental homeostasis and ensure the proper functioning of the organism. In this study, we measured biochemical parameters and performed mRNA-miRNA integrated transcriptomic analysis to characterize changes in gene expression profiles in the muscle tissue of spotted sea bass (Lateolabrax maculatus) under heat stress. The measurement of biochemical parameters revealed that the activities of nine biochemical enzymes (ALP, γ-GT, AST, GLU, CK, ALT, TG, LDH and TC) were significantly affected to varying degrees by elevated temperatures. A total of 1940 overlapping differentially expressed genes (DEGs) were identified among the five comparisons in the muscle tissue after heat stress. Protein-protein interaction (PPI) analysis of DEGs indicated that heat shock protein genes (HSPs) were deeply involved in the response to heat stress. In addition, we detected 462 differential alternative splicing (DAS) events and 618 DAS genes, which are closely associated with sarcomere assembly in muscle, highlighting the role of alternative splicing in thermal response regulation. Moreover, 32 differentially expressed miRNAs (DEMs) were identified in response to heat stress, and 599 DEGs were predicted as potential target genes of those DEMs, generating 846 DEG-DEM negative regulatory pairs potentially associated with thermal response. Function enrichment analysis of the target genes suggested that lipid metabolism-related pathways and genes were regulated by miRNAs. By analyzing PPIs of target genes, we identified 28 key negative regulatory pairs, including 13 miRNAs (such as lma-miR-122, lma-miR-200b-5p and novel-miR-444) and 15 target genes (such as hspa13, dnaja1, and dnajb1a). This study elucidates the molecular mechanisms of response to high-temperature stress and offers valuable information for the selection and breeding of heat-tolerant strains of spotted sea bass.

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy >90% using a machine learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Exploring genetic diversity, population structure, and subgenome differences in the allopolyploid Camelina sativa: implications for future breeding and research studies

Camelina (Camelina sativa), an allohexaploid species, is an emerging aviation biofuel crop that has been the focus of resurgent interest in recent decades. To guide future breeding and crop improvement efforts, the community requires a deeper comprehension of subgenome dominance, often noted in allopolyploid species, "alongside an understanding of the genetic diversity" and population structure of material present within breeding programs. We conducted population genetic analyses of a C. sativa diversity panel, leveraging a new genome, to estimate nucleotide diversity and population structure, and analyzed for patterns of subgenome expression dominance among different organs. Our analyses confirm that C. sativa has relatively low genetic diversity and show that the SG3 subgenome has substantially lower genetic diversity compared to the other two subgenomes. Despite the low genetic diversity, our analyses identified 13 distinct subpopulations including two distinct wild populations and others putatively representing founders in existing breeding populations. When analyzing for subgenome composition of long non-coding RNAs, which are known to play important roles in (a)biotic stress tolerance, we found that the SG3 subgenome contained significantly more lincRNAs compared to other subgenomes. Similarly, transcriptome analyses revealed that expression dominance of SG3 is not as strong as previously reported and may not be universal across all organ types. From a global analysis, SG3 "was only significant higher expressed" in flower, flower bud, and fruit organs, which is an important discovery given that the crop yield is associated with these organs. Collectively, these results will be valuable for guiding future breeding efforts in camelina.

Multiomic characterization of RNA microenvironments by oligonucleotide-mediated proximity-interactome mapping

RNA molecules form complex networks of molecular interactions that are central to their function and to cellular architecture. But these interaction networks are difficult to probe in situ. Here, we introduce Oligonucleotide-mediated proximity-interactome MAPping (O-MAP), a method for elucidating the biomolecules near an RNA of interest, within its native context. O-MAP uses RNA-fluorescence in situ hybridization-like oligonucleotide probes to deliver proximity-biotinylating enzymes to a target RNA in situ, enabling nearby molecules to be enriched by streptavidin pulldown. This induces exceptionally precise biotinylation that can be easily optimized and ported to new targets or sample types. Using the noncoding RNAs 47S, 7SK and Xist as models, we develop O-MAP workflows for discovering RNA-proximal proteins, transcripts and genomic loci, yielding a multiomic characterization of these RNAs' subcellular compartments and new regulatory interactions. O-MAP requires no genetic manipulation, uses exclusively off-the-shelf parts and requires orders of magnitude fewer cells than established methods, making it accessible to most laboratories.

Subtranscriptome analysis of phospholipases D in Loxosceles venom glands: Confirmation of predominance, intra-species diversity, and description of novel isoforms

Spiders of Loxosceles genus, or Brown spiders produce a potent venom with minimal volume and protein content. Among its toxins, phospholipases D (PLDs) are notable for causing primary local and systemic manifestations observed following envenomation. They degrade cellular phospholipids, mainly sphingomyelin and lysophosphatidylcholine. We present a robust and detailed analysis of PLD transcripts from venom glands of three major clinically relevant South American species-L. intermedia, L. laeta, and L. gaucho-using next-generation sequencing. Results confirmed that PLDs are the most highly expressed toxins, accounting for 65.4 % of expression in L. intermedia, 71.8 % in L. gaucho, and 50.4 % in L. laeta. These findings further support the idea that these enzymes form a protein family both within and across species. Eighteen contigs for PLDs were found for L. gaucho, 24 for L. intermedia, and 21 for L. laeta. A detailed analysis revealed that, although all contigs display conserved amino acid residues directly involved in catalysis, magnesium coordination, and substrate affinity, they also possess distinct primary sequences with important substitutions. Such data reinforces the hypothesis that these toxins may act synergistically. Furthermore, new PLD sequences were identified within the contigs. For L. intermedia, 14 potential new isoforms were identified; 16 for L gaucho; and 16 novel sequences for L. laeta. This indicates that there is still a wealth of undisclosed information about these toxins. These data will help identify structural and functional differences among these proteins, support future functional studies, and to the comprehensive understanding of the mechanism of action of PLDs.

Contrasting and combining transcriptome complexity captured by short and long RNA sequencing reads

Mapping transcriptomic variations using either short- or long-read RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, whereas short reads still provide improved coverage and error rates. Yet, open questions remain, such as how to quantitatively compare the technologies, can we combine them, and what is the benefit of such a combined view? We tackle these questions by first creating a pipeline to assess matched long- and short-read data using a variety of transcriptome statistics. We find that across data sets, algorithms, and technologies, matched short-read data detects ∼30% more splice junctions, such that ∼10%-30% of the splice junctions included at ≥20% by short reads are missed by long reads. In contrast, long reads detect many more intron-retention events and can detect full isoforms, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software, to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long-read technology or algorithm and can be combined with short-read data for improved transcriptome analysis.

Deciphering transcript architectural complexity in bacteria and archaea

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Safety Assessment of Camelid-Derived Single-Domain Antibody as Feed Additive for Juvenile Whiteleg Shrimp (Litopenaeus vannamei) Against White Spot Syndrome Virus

A six-week feeding trial was conducted to assess the safety of single-domain antibodies (sdAbs) derived from camelids against the white spot syndrome virus (WSSV) (WSSVvp28 was used as the antigen), focusing on the whole-organism responses and molecular-level changes in juvenile whiteleg shrimp (Litopenaeus vannamei). Five experimental diets with varying levels of sdAbs were formulated: CON (no sdAb supplementation); SDA8.2 (8.20% of sdAbs); SDA16.4 (16.40% of sdAbs); SDA24.6 (24.60% of sdAbs); and SDA32.8 (32.80% of sdAbs). In the CON diet, 450 mL of water per kg of diet (45%) was used to form a feed dough, while sdAbs were used to replace the water in the treatment diets. A total of 450 shrimp, with an initial body weight of 3.27 ± 0.02 g (mean ± SEM), were randomly distributed in 15 tanks (30 shrimp per tank; three tanks per treatment). Each tank was filled with 30 L of seawater (77 L capacity) in an indoor semi-recirculating system with a constant water flow rate of 1.2 L min-1. The photoperiod was maintained at 12 h of light and 12 h of dark. The water temperature, pH, salinity, and dissolved oxygen were 27.3 ± 0.1 °C, 7.61 ± 0.01, 34 ± 1 ppt, and 5.94 ± 0.04 mg L-1, respectively. During the feeding trial, the shrimp were fed the experimental diet (40% protein and 11% lipid) three times a day for six weeks. Following the feeding trial, an acute cold-water-temperature stress test was conducted by abruptly exposing the shrimp from each treatment to 15 °C for 4 h, down from 27 °C. The results showed no significant differences in the growth performance (weight gain, feed utilization efficiency, survival, etc.), plasma metabolites (aspartate aminotransferase activity, alanine aminotransferase activity, total protein, and glucose), or antioxidant enzymes (superoxide dismutase and glutathione peroxidase) among all the experimental diets (p > 0.05). In the acute cold-temperature stress test, there was no significant interaction between sdAb supplementation and temperature stress, nor any main effect from either factor, except for the main effect of temperature stress on the glucose levels, which was significantly higher in shrimp exposed to cold-temperature stress (p < 0.05). The next-generation sequencing of differentially expressed genes (DEGs) in the hepatopancreases of shrimp fed the CON, SDA16.4, and SDA32.8 diets, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, indicated that DEGs were significantly enriched in signaling pathways associated with growth, cold stress, and antioxidant systems. Overall, the results from conventional measurements suggest that the use of sdAbs against the WSSV may be safe for juvenile whiteleg shrimp. However, findings from the sophisticated analysis indicate that further research is needed to understand the molecular mechanisms underlying the observed changes, and to evaluate the long-term effects of sdAb supplementation in shrimp diets.

Non-coding cause of congenital heart defects: Abnormal RNA splicing with multiple isoforms as a mechanism for heterotaxy

Heterotaxy is a disorder characterized by severe congenital heart defects (CHDs) and abnormal left-right patterning in other thoracic or abdominal organs. Clinical and research-based genetic testing has previously focused on evaluation of coding variants to identify causes of CHDs, leaving non-coding causes of CHDs largely unknown. Variants in the transcription factor zinc finger of the cerebellum 3 (ZIC3) cause X-linked heterotaxy. We identified an X-linked heterotaxy pedigree without a coding variant in ZIC3. Whole-genome sequencing revealed a deep intronic variant (ZIC3 c.1224+3286A>G) predicted to alter RNA splicing. An in vitro minigene splicing assay confirmed the variant acts as a cryptic splice acceptor. CRISPR-Cas9 served to introduce the ZIC3 c.1224+3286A>G variant into human embryonic stem cells demonstrating pseudoexon inclusion caused by the variant. Surprisingly, Sanger sequencing of the resulting ZIC3 c.1224+3286A>G amplicons revealed several isoforms, many of which bypass the normal coding sequence of the third exon of ZIC3, causing a disruption of a DNA-binding domain and a nuclear localization signal. Short- and long-read mRNA sequencing confirmed these initial results and identified additional splicing patterns. Assessment of four isoforms determined abnormal functions in vitro and in vivo while treatment with a splice-blocking morpholino partially rescued ZIC3. These results demonstrate that pseudoexon inclusion in ZIC3 can cause heterotaxy and provide functional validation of non-coding disease causation. Our results suggest the importance of non-coding variants in heterotaxy and the need for improved methods to identify and classify non-coding variation that may contribute to CHDs.

The loach haplotype-resolved genome and the identification of Mex3a involved in fish air breathing

Fish air breathing is crucial for the transition of vertebrates from water to land. So far, the genes involved in fish air breathing have not been well identified. Here, we performed gene enrichment analysis of positively selected genes (PSGs) in loach (Misgurnus anguillicaudatus, an air-breathing fish) in comparison to Triplophysa tibetana (a non-air-breathing fish), haplotype-resolved genome assembly of the loach, and gene evolutionary analysis of air-breathing and non-air-breathing fishes and found that the PSG mex3a originated from ancient air-breathing fish species. Deletion of Mex3a impaired loach air-breathing capacity by inhibiting angiogenesis through its interaction with T-box transcription factor 20. Mex3a overexpression significantly promoted angiogenesis. Structural analysis and point mutation revealed the critical role of the 201st amino acid in loach Mex3a for angiogenesis. Our findings innovatively indicate that the ancient mex3a is a fish air-breathing gene, which holds significance for understanding fish air breathing and provides a valuable resource for cultivating hypoxia-tolerant fish varieties.

Haplotype-resolved and near-T2T genome assembly of the African catfish (Clarias gariepinus)

Airbreathing catfish are stenohaline freshwater fish capable of withstanding various environmental conditions and farming practices, including breathing atmospheric oxygen. This unique ability has enabled them to thrive in semi-terrestrial habitats. However, the genomic mechanisms underlying their adaptation to adverse ecological environments remain largely unexplored, primarily due to the limited availability of high-quality genomic resources. Here, we present a haplotype-resolved and near telomere-to-telomere (T2T) genome assembly of the African catfish (Clarias gariepinus), utilizing Oxford Nanopore, PacBio HiFi, Illumina and Hi-C sequencing technologies. The primary assembly spans 969.62 Mb with only 47 contigs, achieving a contig N50 of 33.71 Mb. Terminal telomeric signals were detected in 22 of 47 contigs, suggesting T2T assembled chromosomes. BUSCO analysis confirmed gene space completeness of 99% against the Actinopterygii dataset, highlighting the high quality of the assembly. Genome annotation identified 25,655 protein-coding genes and estimated 43.94% genome-wide repetitive elements. This data provides valuable genomic resources to advance aquaculture practices and to explore the genomic underpinnings of the ecological resilience of airbreathing catfish and related teleosts.

Evolution of translational control and the emergence of genes and open reading frames in human and non-human primate hearts

Evolutionary innovations can be driven by changes in the rates of RNA translation and the emergence of new genes and small open reading frames (sORFs). In this study, we characterized the transcriptional and translational landscape of the hearts of four primate and two rodent species through integrative ribosome and transcriptomic profiling, including adult left ventricle tissues and induced pluripotent stem cell-derived cardiomyocyte cell cultures. We show here that the translational efficiencies of subunits of the mitochondrial oxidative phosphorylation chain complexes IV and V evolved rapidly across mammalian evolution. Moreover, we discovered hundreds of species-specific and lineage-specific genomic innovations that emerged during primate evolution in the heart, including 551 genes, 504 sORFs and 76 evolutionarily conserved genes displaying human-specific cardiac-enriched expression. Overall, our work describes the evolutionary processes and mechanisms that have shaped cardiac transcription and translation in recent primate evolution and sheds light on how these can contribute to cardiac development and disease.