The Ensembl gene annotation system

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

The genome sequence of a digger wasp, Ectemnius lituratus (Panzer,1805)

We present a genome assembly from an individual female Ectemnius lituratus (a digger wasp; Arthropoda; Insecta; Hymenoptera; Crabronidae). The genome sequence is 235.1 megabases in span. Most of the assembly is scaffolded into 13 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 29.67 kilobases in length. Gene annotation of this assembly on Ensembl identified 9,724 protein coding genes.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and non-coding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared to other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared to a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

The genome sequence of the Chestnut, Conistra vaccinii (Linnaeus, 1761)

We present a genome assembly from an individual male Conistra vaccinii (the Chestnut; Arthropoda; None; Lepidoptera; Noctuidae). The genome sequence is 720.8 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.44 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,109 protein coding genes.

The genome sequence of the Beautiful Hook-tip, Laspeyria flexula (Denis & Schiffermüller, 1775)

We present a genome assembly from an individual male Laspeyria flexula (the Beautiful Hook-tip; Arthropoda; Insecta; Lepidoptera; Erebidae). The genome sequence is 450.9 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.58 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,281 protein coding genes.

The genome sequence of the Six-belted Clearwing, Bembecia ichneumoniformis (Denis & Schiffermüller, 1775)

We present a genome assembly from an individual male Bembecia ichneumoniformis (the Six-belted Clearwing; Arthropoda; Insecta; Lepidoptera; Sesiidae). The genome sequence is 511.4 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.32 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,114 protein coding genes.

The transcriptional legacy of developmental stochasticity

Genetic and environmental variation are key contributors during organism development, but the influence of minor perturbations or noise is difficult to assess. This study focuses on the stochastic variation in allele-specific expression that persists through cell divisions in the nine-banded armadillo (Dasypus novemcinctus). We investigated the blood transcriptome of five wild monozygotic quadruplets over time to explore the influence of developmental stochasticity on gene expression. We identify an enduring signal of autosomal allelic variability that distinguishes individuals within a quadruplet despite their genetic similarity. This stochastic allelic variation, akin to X-inactivation but broader, provides insight into non-genetic influences on phenotype. The presence of stochastically canalized allelic signatures represents a novel axis for characterizing organismal variability, complementing traditional approaches based on genetic and environmental factors. We also developed a model to explain the inconsistent penetrance associated with these stochastically canalized allelic expressions. By elucidating mechanisms underlying the persistence of allele-specific expression, we enhance understanding of development's role in shaping organismal diversity.

The genome sequence of a darkling beetle, Lagria hirta (Linnaeus, 1758)

We present a genome assembly from an individual male Lagria hirta (darkling beetle; Arthropoda; Insecta; Coleoptera; Tenebrionidae). The genome sequence is 336.8 megabases in span. Most of the assembly is scaffolded into 12 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 16.12 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,850protein coding genes.

The genome sequence of the Saxon Wasp, Dolichovespula saxonica (Fabricius, 1793)

We present a genome assembly from an individual male Dolichovespula saxonica (the Saxon Wasp; Arthropoda; Insecta; Hymenoptera; Vespidae). The genome sequence is 221.8 megabases in span. Most of the assembly is scaffolded into 26 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 18.97 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,856 protein coding genes.

The genome sequence of Fabricius' Nomad Bee, Nomada fabriciana (Linne, 1767)

We present a genome assembly from an individual female Nomada fabriciana (Fabricius' Nomad Bee; Arthropoda; Insecta; Hymenoptera; Apidae). The genome sequence is 233.6 megabases in span. Most of the assembly is scaffolded into 12 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 19.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 9,700 protein coding genes.

The genome sequence of a digger wasp, Ectemnius continuus (Fabricius, 1804)

We present a genome assembly from an individual female Ectemnius continuus (digger wasp; Arthropoda; Insecta; Hymenoptera; Crabronidae). The genome sequence is 260.3 megabases in span. Most of the assembly is scaffolded into 14 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 27.05 kilobases in length. Gene annotation of this assembly on Ensembl identified 9,835 protein coding genes.

The genome sequence of the small wasp-sawfly, Tenthredo distinguenda (R. Stein, 1885)

We present a genome assembly from an individual male Tenthredo distinguenda (the small wasp-sawfly; Arthropoda; Insecta; Hymenoptera; Tenthredinidae). The genome sequence is 229.4 megabases in span. Most of the assembly is scaffolded into 9 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 31.6 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,332 protein coding genes.

The genome sequence of the Ingrailed Clay, Diarsia mendica (Fabricius, 1775)

We present a genome assembly from an individual male Diarsia mendica (the Ingrailed Clay; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 727.9 megabases in span. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.37 kilobases in length. Gene annotation of this assembly on Ensembl identified 14,077 protein coding genes.

Transcriptomics Meta-Analysis Reveals Phagosome and Innate Immune System Dysfunction as Potential Mechanisms in the Cortex of Alzheimer's Disease Mouse Strains

Immune-related pathways can affect the immune system directly, such as the chemokine signaling pathway, or indirectly, such as the phagosome pathway. Alzheimer's disease (AD) is reportedly associated with several immune-related pathways. However, exploring its underlying mechanism is challenging in animal studies because AD mouse strains differentially express immune-related pathway characteristics. To overcome this problem, we performed a meta-analysis to identify significant and consistent immune-related AD pathways that are expressed in different AD mouse strains. Next-generation RNA sequencing (RNA-seq) and microarray datasets for the cortex of AD mice from different strains such as APP/PSEN1, APP/PS2, 3xTg, TREM, and 5xFAD were collected from the NCBI GEO database. Each dataset's quality control and normalization were already processed from each original study source using various methods depending on the high-throughput analysis platform (FastQC, median of ratios, RMA, between array normalization). Datasets were analyzed using DESeq2 for RNA-seq and GEO2R for microarray to identify differentially expressed (DE) genes. Significantly DE genes were meta-analyzed using Stouffer's method, with significant genes further analyzed for functional enrichment. Ten datasets representing 20 conditions were obtained from the NCBI GEO database, comprising 116 control and 120 AD samples. The DE analysis identified 284 significant DE genes. The meta-analysis identified three significantly enriched immune-related AD pathways: phagosome, the complement and coagulation cascade, and chemokine signaling. Phagosomes-related genes correlated with complement and immune system. Meanwhile, phagosomes and chemokine signaling genes overlapped with B cells receptors pathway genes indicating potential correlation between phagosome, chemokines, and adaptive immune system as well. The transcriptomic meta-analysis showed that AD is associated with immune-related pathways in the brain's cortex through the phagosome, complement and coagulation cascade, and chemokine signaling pathways. Interestingly, phagosome and chemokine signaling pathways had potential correlation with B cells receptors pathway.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Lmx1a is a master regulator of the cortical hem

Development of the nervous system depends on signaling centers - specialized cellular populations that produce secreted molecules to regulate neurogenesis in the neighboring neuroepithelium. In some cases, signaling center cells also differentiate to produce key types of neurons. The formation of a signaling center involves its induction, the maintenance of expression of its secreted molecules, and cell differentiation and migration events. How these distinct processes are coordinated during signaling center development remains unknown. By performing studies in mice, we show that Lmx1a acts as a master regulator to orchestrate the formation and function of the cortical hem (CH), a critical signaling center that controls hippocampus development. Lmx1a co-regulates CH induction, its Wnt signaling, and the differentiation and migration of CH-derived Cajal-Retzius neurons. Combining RNAseq, genetic, and rescue experiments, we identified major downstream genes that mediate distinct Lmx1a-dependent processes. Our work revealed that signaling centers in the mammalian brain employ master regulatory genes and established a framework for analyzing signaling center development.

The genome sequence of the Diamondback Moth, Plutella xylostella (Linnaeus, 1758)

We present a genome assembly from an individual male Plutella xylostella (the Diamondback Moth; Arthropoda; Insecta; Lepidoptera; Plutellidae). The genome sequence is 323.3 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 35.12 kilobases in length. Gene annotation of this assembly on Ensembl identified 17,190 protein coding genes.

Jian Yun Qing Hua Decoction inhibits malignant behaviors of gastric carcinoma cells via COL12A1 mediated ferroptosis signal pathway

BACKGROUND

Jian Yun Qing Hua Decoction (JYQHD), a traditional Chinese medicine decoction, which has been applied in the treatment of gastric cancer (GC). We attempt to confirm the anti-gastric cancer effect of JYQHD and explore the mechanism of JYQHD.

METHODS

Acute toxicity test was used to understand the toxicity of JYQHD. We studied the expression and prognostic outcome of COL12A1 within GC tissues through the network databases. Using several web-based databases, we analyzed the major components and targets of JYQHD, as well as known therapeutic targets in gastric cancer. The Venn diagram was utilized to obtain the overlapped genes. Lentiviral vector, shRNAs and plasmids, were used to transfect GC cells. Cell counting kit-8 (CCK8), sphere formation, malondialdehyde (MDA), glutathione (GSH), reactive oxygen species (ROS), Fe2+, transmission electron microscopy (TEM), quantitative Real-Time Polymerase Chain Reaction (qRT-PCR), Western-Blot (WB), and immunohistochemical (IHC) assays were employed to investigate the role and mechanism of COL12A1 and JYQHD in GC.

RESULTS

The results showed that JYQHD was non-toxic and safe. JYQHD inhibited growth and sphere formation ability through inducing the ferroptosis of GC cells, and suppressed the GC cells induced subcutaneous xenograft tumor growth. COL12A1 was highly expressed in gastric cancer tissues, indicating poor prognosis. COL12A1 specifically enhanced GC cell progression and stemness via suppressing ferroptosis. JYQHD down-regulated COL12A1 in order to suppress the stemness of GC cells via inducing ferroptosis.

CONCLUSION

COL12A1 inhibited ferroptosis and enhanced stemness in GC cells. JYQHD inhibited the development of GC cells by inhibiting cancer cell stemness via the ferroptosis pathway mediated by COL12A1.

The genome sequence of the Six-striped Rustic, Xestia sexstrigata (Haworth, 1809)

We present a genome assembly from an individual female Xestia sexstrigata (the Six-striped Rustic; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 638.3 megabases in span. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 15.36 kilobases in length. Gene annotation of this assembly on Ensembl identified 15,104 protein coding genes.

The genome sequence of the Orange-tailed Mining Bee, Andrena haemorrhoa (Fabricius, 1781)

We present a genome assembly from an individual female Andrena haemorrhoa (the Orange-tailed Mining Bee; Arthropoda; Insecta; Hymenoptera; Andrenidae). The genome sequence is 330.7 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 16.46 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,908 protein coding genes.

The genome sequence of the painted nomad bee, Nomada fucata (Panzer, 1798)

We present a genome assembly from an individual female Nomada fucata (the painted nomad bee; Arthropoda; Insecta; Hymenoptera; Apidae). The genome sequence is 264.1 megabases in span. Most of the assembly is scaffolded into 8 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 22.06 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,173 protein coding genes.

Evaluating cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico

BACKGROUND

Preclinical models like cancer cell lines and patient-derived xenografts (PDXs) are vital for studying disease mechanisms and evaluating treatment options. It is essential that they accurately recapitulate the disease state of interest to generate results that will translate in the clinic. Prior studies have demonstrated that preclinical models do not recapitulate all biological aspects of human tissues, particularly with respect to the tissue of origin gene expression signatures. Therefore, it is critical to assess how well preclinical model gene expression profiles correlate with human cancer tissues to inform preclinical model selection and data analysis decisions.

AIMS

Here we evaluated how well preclinical models recapitulate human cancer and non-diseased tissue gene expression patterns in silico with respect to the full gene expression profile as well as subsetting by the most variable genes, genes significantly correlated with tumor purity, and tissue-specific genes.

METHODS

By using publicly available gene expression profiles across multiple sources, we evaluated cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico.

RESULTS

We found that using the full gene set improves correlations between preclinical model and tissue global gene expression profiles, confirmed that glioblastoma (GBM) PDX global gene expression correlation to GBM tumor global gene expression outperforms GBM cell line to GBM tumor global gene expression correlations, and demonstrated that preclinical models in our study often failed to reproduce tissue-specific expression. While including additional genes for global gene expression comparison between cell lines and tissues decreases the overall correlation, it improves the relative rank between a cell line and its tissue of origin compared to other tissues. Our findings underscore the importance of using the full gene expression set measured when comparing preclinical models and tissues and confirm that tissue-specific patterns are better preserved in GBM PDX models than in GBM cell lines.

CONCLUSION

Future studies can build on these findings to determine the specific pathways and gene sets recapitulated by particular preclinical models to facilitate model selection for a given study design or goal.

The genome sequence of the dark-edged bee fly, Bombylius major (Linnaeus, 1758)

We present a genome assembly from an individual male Bombylius major (the dark-edged bee fly; Arthropoda; Insecta; Diptera; Bombyliidae). The genome sequence is 304.3 megabases in span. The whole assembly is scaffolded into 7 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.8 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,852 protein coding genes.

The genome sequence of the short-fringed mining bee, Andrena dorsata (Kirby, 1802)

We present a genome assembly from an individual female Andrena dorsata (the short-fringed mining bee; Arthropoda; Insecta; Hymenoptera; Andrenidae). The genome sequence is 277.3 megabases in span. Most of the assembly is scaffolded into 3 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 16.11 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,916 protein coding genes.

The genome sequence of the Beautiful Golden Y, Autographa pulchrina (Haworth, 1809)

We present a genome assembly from an individual female Autographa pulchrina (the Beautiful Golden Y; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 426.2 megabases in span. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the Z and W sex chromosomes. The mitochondrial genome has also been assembled and is 15.25 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,916 protein coding genes.

The genome sequence of the red compost earthworm, Lumbricus rubellus (Hoffmeister, 1843)

We present a genome assembly from an individual Lumbricus rubellus (the red compost earthworm; Annelida; Clitellata; Haplotaxida; Lumbricidae). The genome sequence is 787.5 megabases in span. Most of the assembly is scaffolded into 18 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.81 kilobases in length. Gene annotation of this assembly on Ensembl identified 33,426 protein coding genes.

The genome sequence of the Thick-legged Hoverfly, Syritta pipiens (Linnaeus, 1758)

We present a genome assembly from an individual female Syritta pipiens (the Thick-legged Hoverfly; Arthropoda; Insecta; Diptera; Syrphidae). The genome sequence is 318.5 megabases in span. Most of the assembly is scaffolded into 5 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.76 kilobases in length. Gene annotation of this assembly on Ensembl identified 18,405 protein coding genes.

The genome sequence of the Brown Argus, Aricia agestis (Denis & Schiffermüller, 1775)

We present genome assemblies from two male Aricia agestis specimens (the Brown Argus; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequences are 435.3 and 437.4 megabases in span. Each assembly is scaffolded into 23 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genomes were assembled and are 15.47 and 15.45 kilobases in length. Gene annotation of these assemblies on Ensembl identified 12,688 and 12,654 protein coding genes.

Characterization of a strain-specific CD-1 reference genome reveals potential inter- and intra-strain functional variability

BACKGROUND

CD-1 is an outbred mouse stock that is frequently used in toxicology, pharmacology, and fundamental biomedical research. Although inbred strains are typically better suited for such studies due to minimal genetic variability, outbred stocks confer practical advantages over inbred strains, such as improved breeding performance and low cost. Knowledge of the full genetic variability of CD-1 would make it more useful in toxicology, pharmacology, and fundamental biomedical research.

RESULTS

We performed deep genomic DNA sequencing of CD-1 mice and used the data to identify genome-wide SNPs, indels, and germline transposable elements relative to the mm10 reference genome. We used multiple genome-wide sequencing data types and previously published CD-1 SNPs to validate our called variants. We used the called variants to construct a strain-specific CD-1 reference genome, which we show can improve mappability and reduce experimental biases from genome-wide sequencing data derived from CD-1 mice. Based on previously published ChIP-seq and ATAC-seq data, we find evidence that genetic variation between CD-1 mice can lead to alterations in transcription factor binding. We also identified a number of variants in the coding region of genes which could have effects on translation of genes.

CONCLUSIONS

We have identified millions of previously unidentified CD-1 variants with the potential to confound studies involving CD-1. We used the identified variants to construct a CD-1-specific reference genome, which can improve accuracy and reduce bias when aligning genomics data derived from CD-1 mice.

The genome sequence of the Locust Fly, Stomorhina lunata (Fabricius, 1805)

We present a genome assembly from an individual female Stomorhina lunata (the Locust Fly; Arthropoda; Insecta; Diptera; Rhiniidae). The genome sequence is 728.1 megabases in span. Most of the assembly is scaffolded into 6 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 16.49 kilobases in length. Gene annotation of this assembly on Ensembl identified 18,358 protein coding genes.

Context-aware transcript quantification from long-read RNA-seq data with Bambu

Most approaches to transcript quantification rely on fixed reference annotations; however, the transcriptome is dynamic and depending on the context, such static annotations contain inactive isoforms for some genes, whereas they are incomplete for others. Here we present Bambu, a method that performs machine-learning-based transcript discovery to enable quantification specific to the context of interest using long-read RNA-sequencing. To identify novel transcripts, Bambu estimates the novel discovery rate, which replaces arbitrary per-sample thresholds with a single, interpretable, precision-calibrated parameter. Bambu retains the full-length and unique read counts, enabling accurate quantification in presence of inactive isoforms. Compared to existing methods for transcript discovery, Bambu achieves greater precision without sacrificing sensitivity. We show that context-aware annotations improve quantification for both novel and known transcripts. We apply Bambu to quantify isoforms from repetitive HERVH-LTR7 retrotransposons in human embryonic stem cells, demonstrating the ability for context-specific transcript expression analysis.

Impact of the SARS-CoV-2 nucleocapsid 203K/204R mutations on the inflammatory immune response in COVID-19 severity

BACKGROUND

The excessive inflammatory responses provoked by SARS-CoV-2 infection are critical factors affecting the severity and mortality of COVID-19. Previous work found that two adjacent co-occurring mutations R203K and G204R (KR) on the nucleocapsid (N) protein correlate with increased disease severity in COVID-19 patients. However, links with the host immune response remain unclear.

METHODS

Here, we grouped nasopharyngeal swab samples of COVID-19 patients into two cohorts based on the presence and absence of SARS-CoV-2 nucleocapsid KR mutations. We performed nasopharyngeal transcriptome analysis of age, gender, and ethnicity-matched COVID-19 patients infected with either SARS-CoV-2 with KR mutations in the N protein (KR patients n = 39) or with the wild-type N protein (RG patients n = 39) and compared to healthy controls (n = 34). The impact of KR mutation on immune response was further characterized experimentally by transcriptomic and proteomic profiling of virus-like-particle (VLP) incubated cells.

RESULTS

We observed markedly elevated expression of proinflammatory cytokines, chemokines, and interferon-stimulated (ISGs) genes in the KR patients compared to RG patients. Using nasopharyngeal transcriptome data, we found significantly higher levels of neutrophils and neutrophil-to-lymphocyte (NLR) ratio in KR patients than in the RG patients. Furthermore, transcriptomic and proteomic profiling of VLP incubated cells confirmed a similar hyper-inflammatory response mediated by the KR variant.

CONCLUSIONS

Our data demonstrate an unforeseen connection between nucleocapsid KR mutations and augmented inflammatory immune response in severe COVID-19 patients. These findings provide insights into how mutations in SARS-CoV-2 modulate host immune output and pathogenesis and may contribute to more efficient therapeutics and vaccine development.

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

The genome sequence of a tachinid fly, Tachina lurida (Fabricius, 1781)

We present a genome assembly from an individual female Tachina lurida (a tachinid fly; Arthropoda; Insecta; Diptera; Tachinidae). The genome sequence is 899.2 megabases in span. Most of the assembly is scaffolded into 6 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 17.3 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,127 protein coding genes.

Rapid genome functional annotation pipeline anchored to the house sparrow (Passer domesticus, Linnaeus 1758) genome reannotation

The house sparrow (Passer domesticus) is a valuable avian model for studying evolutionary genetics, development, neurobiology, physiology, behavior, and ecology, both in laboratory and field-based settings. The current annotation of the P. domesticus genome available at the Ensembl Rapid Release site is primarily focused on gene set building and lacks functional information. In this study, we present the first comprehensive functional reannotation of the P. domesticus genome using intestinal Illumina RNA sequencing (RNA-Seq) libraries. Our revised annotation provides an expanded view of the genome, encompassing 38592 transcripts compared to the current 23574 transcripts in Ensembl. We also predicted 14717 protein-coding genes, achieving 96.4% completeness for Passeriformes lineage BUSCOs. A substantial improvement in this reannotation is the accurate delineation of untranslated region (UTR) sequences. We identified 82.7% and 93.8% of the transcripts containing 5'- and 3'-UTRs, respectively. These UTR annotations are crucial for understanding post-transcriptional regulatory processes. Our findings underscore the advantages of incorporating additional specific RNA-Seq data into genome annotation, particularly when leveraging fast and efficient data processing capabilities. This functional reannotation enhances our understanding of the P. domesticus genome, providing valuable resources for future investigations in various research fields.

The genome sequence of the Dark Spectacle, Abrostola triplasia (Linnaeus, 1758)

We present a genome assembly from an individual male Abrostola triplasia (the Dark Spectacle; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 362.7 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.34 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,532 protein coding genes.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

The genome sequence of the narrow-cheeked clusterfly, Pollenia angustigena (Wainwright, 1940)

We present a genome assembly from an individual female Pollenia angustigena (the narrow-cheeked clusterfly; Arthropoda; Insecta; Diptera; Polleniidae). The genome sequence is 1370.5 megabases in span. Most of the assembly is scaffolded into 6 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 21.01 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,930 protein coding genes.

The genome sequence of the Tawny Mining Bee, Andrena fulva (Müller, 1766)

We present a genome assembly from an individual female Andrena fulva (the Tawny Mining Bee; Arthropoda; Insecta; Hymenoptera; Andrenidae). The genome sequence is 461.7 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 14.7 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,011 protein coding genes.

Malignant A-to-I RNA editing by ADAR1 drives T-cell acute lymphoblastic leukemia relapse via attenuating dsRNA sensing

Leukemia initiating cells (LICs) are regarded as the origin of leukemia relapse and therapeutic resistance. Identifying direct stemness determinants that fuel LIC self-renewal is critical for developing targeted approaches to eliminate LICs and prevent relapse. Here, we show that the RNA editing enzyme ADAR1 is a crucial stemness factor that promotes LIC self-renewal by attenuating aberrant double-stranded RNA (dsRNA) sensing. Elevated adenosine-to-inosine (A-to-I) editing is a common attribute of relapsed T-ALL regardless of molecular subtypes. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL PDX models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA and retains unedited nuclear dsRNA to avoid detection by the innate immune sensor MDA5. Moreover, we uncovered that the cell intrinsic level of MDA5 dictates the dependency on ADAR1-MDA5 axis in T-ALL. Collectively, our results show that ADAR1 functions as a self-renewal factor that limits the sensing of endogenous dsRNA. Thus, targeting ADAR1 presents a safe and effective therapeutic strategy for eliminating T-ALL LICs.

Evaluating cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico

Preclinical models like cancer cell lines and patient-derived xenografts (PDXs) are vital for studying disease mechanisms and evaluating treatment options. It is essential that they accurately recapitulate the disease state of interest to generate results that will translate in the clinic. Prior studies have demonstrated that preclinical models do not recapitulate all biological aspects of human tissues, particularly with respect to the tissue of origin gene expression signatures. Therefore, it is critical to assess how well preclinical model gene expression profiles correlate with human cancer tissues to inform preclinical model selection and data analysis decisions. Here we evaluated how well preclinical models recapitulate human cancer and non-diseased tissue gene expression patterns in silico with respect to the full gene expression profile as well as subsetting by the most variable genes, genes significantly correlated with tumor purity, and tissue-specific genes by using publicly available gene expression profiles across multiple sources. We found that using the full gene set improves correlations between preclinical model and tissue global gene expression profiles, confirmed that GBM PDX global gene expression correlation to GBM tumor global gene expression outperforms GBM cell line to GBM tumor global gene expression correlations, and demonstrated that preclinical models in our study often failed to reproduce tissue-specific expression. While including additional genes for global gene expression comparison between cell lines and tissues decreases the overall correlation, it improves the relative rank between a cell line and its tissue of origin compared to other tissues. Our findings underscore the importance of using the full gene expression set measured when comparing preclinical models and tissues and confirm that tissue-specific patterns are better preserved in GBM PDX models than in GBM cell lines. Future studies can build on these findings to determine the specific pathways and gene sets recapitulated by particular preclinical models to facilitate model selection for a given study design or goal.

The genome sequence of a carabid beetle, Nebria salina (Fairmaire & Laboulbène, 1854)

We present a genome assembly from an individual female Nebria salina (a carabid beetle; Arthropoda; Insecta; Coleoptera; Carabidae). The genome sequence is 256.7 megabases in span. Most of the assembly is scaffolded into 21 chromosomal pseudomolecules, including the assembled X sex chromosome. The mitochondrial genome has also been assembled and is 24.7 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,671 protein coding genes.

Specific DMPK-promoter targeting by CRISPRi reverses myotonic dystrophy type 1-associated defects in patient muscle cells

Myotonic dystrophy type 1 (DM1) is a neuromuscular disease that originates from an expansion of CTG microsatellites in the 3' untranslated region of the DMPK gene, thus leading to the expression of transcripts containing expanded CUG repeats (CUGexp). The pathophysiology is explained by a toxic RNA gain of function where CUGexp RNAs form nuclear aggregates that sequester and alter the function of MBNL splicing factors, triggering splicing misregulation linked to the DM1 symptoms. There is currently no cure for DM1, and most therapeutic strategies aim at eliminating CUGexp-DMPK transcripts. Here, we investigate a DMPK-promoter silencing strategy using CRISPR interference as a new alternative approach. Different sgRNAs targeting the DMPK promoter are evaluated in DM1 patient muscle cells. The most effective guides allowed us to reduce the level of DMPK transcripts and CUGexp-RNA aggregates up to 80%. The CUGexp-DMPK repression corrects the overall transcriptome, including spliceopathy, and reverses a physiological parameter in DM1 muscle cells. Its action is specific and restricted to the DMPK gene, as confirmed by genome-wide expression analysis. Altogether, our findings highlight DMPK-promoter silencing by CRISPRi as a promising therapeutic approach for DM1.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Spheroids as a model for endometriotic lesions

The development and progression of endometriotic lesions are poorly understood, but immune cell dysfunction and inflammation are closely associated with the pathophysiology of endometriosis. There is a need for 3D in vitro models to permit the study of interactions between cell types and the microenvironment. To address this, we developed endometriotic spheroids (ES) to explore the role of epithelial-stromal interactions and model peritoneal invasion associated with lesion development. Using a nonadherent microwell culture system, spheroids were generated with immortalized endometriotic epithelial cells (12Z) combined with endometriotic stromal (iEc-ESC) or uterine stromal (iHUF) cell lines. Transcriptomic analysis found 4,522 differentially expressed genes in ES compared with spheroids containing uterine stromal cells. The top increased gene sets were inflammation-related pathways, and an overlap with baboon endometriotic lesions was highly significant. Finally, to mimic invasion of endometrial tissue into the peritoneum, a model was developed with human peritoneal mesothelial cells in an extracellular matrix. Invasion was increased in the presence of estradiol or pro-inflammatory macrophages and suppressed by a progestin. Taken together, our results strongly support the concept that ES are an appropriate model for dissecting mechanisms that contribute to endometriotic lesion development.

The genome sequence of the variegated flesh fly, Sarcophaga variegata (Scopoli, 1763)

We present a genome assembly from an individual male Sarcophaga variegata (the variegated flesh fly; Arthropoda; Insecta; Diptera; Sarcophagidae). The genome sequence is 718.5 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 18.7 kilobases in length. Gene annotation of this assembly on Ensembl identified 16,660 protein coding genes.

The genome sequence of the Lichen Button, Acleris literana (Linnaeus, 1758)

We present a genome assembly from an individual male Acleris literana (the Lichen Button; Arthropoda; Insecta; Lepidoptera; Tortricidae). The genome sequence is 674.9 megabases in span. Most of the assembly is scaffolded into 30 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 16.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,577 protein coding genes.

Microarray analysis of hub genes, non-coding RNAs and pathways in lung after whole body irradiation in a mouse model

PURPOSE

Previous research has highlighted the impact of radiation damage, with cancer patients developing acute disorders including radiation induced pneumonitis or chronic disorders including pulmonary fibrosis months after radiation therapy ends. We sought to discover biomarkers that predict these injuries and develop treatments that mitigate this damage and improve quality of life.

MATERIALS AND METHODS

Six- to eight-week-old female C57BL/6 mice received 1, 2, 4, 8, 12 Gy or sham whole body irradiation. Animals were euthanized 48 h post exposure and lungs removed, snap frozen and underwent RNA isolation. Microarray analysis was performed to determine dysregulation of messenger RNA (mRNA), microRNA (miRNA), and long non-coding RNA (lncRNA) after radiation injury.

RESULTS

We observed sustained dysregulation of specific RNA markers including: mRNAs, lncRNAs, and miRNAs across all doses. We also identified significantly upregulated genes that can indicate high dose exposure, including Cpt1c, Pdk4, Gdf15, and Eda2r, which are markers of senescence and fibrosis. Only three miRNAs were significantly dysregulated across all radiation doses: miRNA-142-3p and miRNA-142-5p were downregulated and miRNA-34a-5p was upregulated. IPA analysis predicted inhibition of several molecular pathways with increasing doses of radiation, including: T cell development, Quantity of leukocytes, Quantity of lymphocytes, and Cell viability.

CONCLUSIONS

These RNA biomarkers might be highly relevant in the development of treatments and in predicting normal tissue injury in patients undergoing radiation treatment. We are conducting further experiments in our laboratory, which includes a human lung-on-a-chip model, to develop a decision tree model using RNA biomarkers.

The genome sequence of the surf clam, Spisula solida (Linnaeus, 1758)

We present a genome assembly from an individual Spisula solida (the surf clam; Mollusca; Bivalvia; Venerida; Mactridae). The genome sequence is 932.1 megabases in span. Most of the assembly is scaffolded into 19 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 19.3 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,833 protein coding genes.

LCAT: an isoform-sensitive error correction for transcriptome sequencing long reads

As the carrier of genetic information, RNA carries the information from genes to proteins. Transcriptome sequencing technology is an important way to obtain transcriptome sequences, and it is also the basis for transcriptome research. With the development of third-generation sequencing, long reads can cover full-length transcripts and reflect the composition of different isoforms. However, the high error rate of third-generation sequencing affects the accuracy of long reads and downstream analysis. The current error correction methods seldom consider the existence of different isoforms in RNA, which makes the diversity of isoforms a serious loss. Here, we introduce LCAT (long-read error correction algorithm for transcriptome sequencing data), a wrapper algorithm of MECAT, to reduce the loss of isoform diversity while keeping MECAT's error correction performance. The experimental results show that LCAT can not only improve the quality of transcriptome sequencing long reads but also retain the diversity of isoforms.