The Ensembl gene annotation system

The genome sequence of Blair's Shoulder-knot, Lithophane leautieri (Boisduval, 1829)

We present a genome assembly from an individual male Lithophane leautieri (Blair's Shoulder-knot; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 521.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.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,254 protein coding genes.

The genome sequence of the Barred Chestnut moth, Diarsia dahlii (Hübner, 1813)

We present a genome assembly from an individual female Diarsia dahlii (the Barred Chestnut; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 683.0 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.36 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,177 protein coding genes.

Specific microRNAs in stallion spermatozoa are potential biomarkers of high functionality

Males of some species, from horses to humans, require medical help for subfertility problems. There is an urgent need for novel molecular assays that reflect spermatozoal function. In the last 25 years, studies examined RNAs in spermatozoa as a window into gene expression during their development and, more recently, for their functions in early embryo development. In clinics, more dense spermatozoa are isolated by density gradient centrifugation before use in artificial insemination to increase pregnancy rates. The objectives of the current study were to discover and quantify the microRNAs in stallion spermatozoa and identify those with differential expression levels in more dense versus less dense spermatozoa. First, spermatozoa from seven stallions were separated into more dense and less dense populations by density gradient centrifugation. Next, small RNAs were sequenced from each of the 14 RNA samples. We identified 287 different mature microRNAs within the 11,824,720 total mature miRNA reads from stallion spermatozoa. The most prevalent was miR-10a/b-5p. The less dense spermatozoa had fewer mature microRNAs and more microRNA precursor sequences than more dense spermatozoa, perhaps indicating that less dense spermatozoa are less mature. Two of the most prevalent microRNAs in more dense stallion spermatozoa were predicted to target mRNAs that encode proteins that accelerate mRNA decay. Nine microRNAs were more highly expressed in more dense spermatozoa. Three of those microRNAs were predicted to target mRNAs that encode proteins involved in protein decay. Both mRNA and protein decay are very active in late spermiogenesis but not in mature spermatozoa. The identified microRNAs may be part of the mechanism to shut down those processes. The microRNAs with greater expression in more dense spermatozoa may be useful biomarkers for spermatozoa with greater functional capabilities.

The genome sequence of the Banded Burying beetle, Nicrophorus investigator Zetterstedt, 1824

We present a genome assembly from a female Nicrophorus investigator (Banded Burying beetle; Arthropoda; Insecta; Coleoptera; Silphidae). The genome sequence is 202.3 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 23.3 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,046 protein coding genes.

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Diminished Immune Cell Adhesion in Hypoimmune ICAM-1 Knockout Pluripotent Stem Cells

Hypoimmune gene edited human pluripotent stem cells (hPSCs) are a promising platform for developing reparative cellular therapies that evade immune rejection. Existing first-generation hypoimmune strategies have used CRISPR/Cas9 editing to modulate genes associated with adaptive (e.g., T cell) immune responses, but have largely not addressed the innate immune cells (e.g., monocytes, neutrophils) that mediate inflammation and rejection processes occurring early after graft transplantation. We identified the adhesion molecule ICAM-1 as a novel hypoimmune target that plays multiple critical roles in both adaptive and innate immune responses post-transplantation. In a series of studies, we found that ICAM-1 blocking or knock-out (KO) in hPSC-derived cardiovascular therapies imparted significantly diminished binding of multiple immune cell types. ICAM-1 KO resulted in diminished T cell proliferation responses in vitro and in longer in vivo retention/protection of KO grafts following immune cell encounter in NeoThy humanized mice. The ICAM-1 KO edit was also introduced into existing first-generation hypoimmune hPSCs and prevented immune cell binding, thereby enhancing the overall hypoimmune capacity of the cells. This novel hypoimmune editing strategy has the potential to improve the long-term efficacy and safety profiles of regenerative therapies for cardiovascular pathologies and a number of other diseases.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Meta-analysis towards FSHD reveals misregulation of neuromuscular junction, nuclear envelope, and spliceosome

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most common autosomal dominant muscle disorders, yet no cure or amelioration exists. The clinical presentation is diverse, making it difficult to identify the actual driving pathomechanism among many downstream events. To unravel this complexity, we performed a meta-analysis of 13 original omics datasets (in total 171 FSHD and 129 control samples). Our approach confirmed previous findings about the disease pathology and specified them further. We confirmed increased expression of former proposed DUX4 biomarkers, and furthermore impairment of the respiratory chain. Notably, the meta-analysis provides insights about so far not reported pathways, including misregulation of neuromuscular junction protein encoding genes, downregulation of the spliceosome, and extensive alterations of nuclear envelope protein expression. Finally, we developed a publicly available shiny app to provide a platform for researchers who want to search our analysis for genes of interest in the future.

The genome sequence of the Lulworth Skipper, Thymelicus acteon (Rottemburg, 1775)

We present a genome assembly from an individual male Thymelicus acteon (the Lulworth Skipper; Arthropoda; Insecta; Lepidoptera; Hesperiidae). The genome sequence is 537.0 megabases in span. Most of the assembly is scaffolded into 28 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 17.08 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,813 protein-coding genes.

Assessing the reliability of point mutation as data augmentation for deep learning with genomic data

BACKGROUND

Deep neural networks (DNNs) have the potential to revolutionize our understanding and treatment of genetic diseases. An inherent limitation of deep neural networks, however, is their high demand for data during training. To overcome this challenge, other fields, such as computer vision, use various data augmentation techniques to artificially increase the available training data for DNNs. Unfortunately, most data augmentation techniques used in other domains do not transfer well to genomic data.

RESULTS

Most genomic data possesses peculiar properties and data augmentations may significantly alter the intrinsic properties of the data. In this work, we propose a novel data augmentation technique for genomic data inspired by biology: point mutations. By employing point mutations as substitutes for codons, we demonstrate that our newly proposed data augmentation technique enhances the performance of DNNs across various genomic tasks that involve coding regions, such as translation initiation and splice site detection.

CONCLUSION

Silent and missense mutations are found to positively influence effectiveness, while nonsense mutations and random mutations in non-coding regions generally lead to degradation. Overall, point mutation-based augmentations in genomic datasets present valuable opportunities for improving the accuracy and reliability of predictive models for DNA sequences.

The genome sequence of the marbled rockcod, Notothenia rossii Richardson, 1844

We present a genome assembly from an individual Notothenia rossii (the marbled rockcod; Chordata; Actinopterygii; Perciformes; Nototheniidae). The genome sequence is 1,042.9 megabases in span. Most of the assembly is scaffolded into 12 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 21.68 kilobases in length. Gene annotation of this assembly on Ensembl identified 24,432 protein coding genes.

The genome sequence of the Elm Groundling moth, Carpatolechia fugitivella (Zeller, 1839)

We present a genome assembly from an individual male Carpatolechia fugitivella (the Elm Groundling; Arthropoda; Insecta; Lepidoptera; Gelechiidae). The genome sequence is 493.1 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 15.26 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,721 protein coding genes.

The genome sequence of the Atlantic cod, Gadus morhua (Linnaeus, 1758)

We present a genome assembly from an individual male Gadus morhua (the Atlantic cod; Chordata; Actinopteri; Gadiformes; Gadidae). The genome sequence is 669.9 megabases in span. Most of the assembly is scaffolded into 23 chromosomal pseudomolecules. Gene annotation of this assembly on Ensembl identified 23,515 protein coding genes.

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 with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. 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,700 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.

The genome sequence of the Spruce-seed moth, Cydia strobilella (Linnaeus, 1758)

We present a genome assembly from an individual male Cydia strobilella (the spruce-seed moth; Arthropoda; Insecta; Lepidoptera; Tortricidae). The genome sequence is 542.6 megabases in span. Most of the assembly is scaffolded into 28 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 16.78 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,058 protein coding genes.

The genome sequence of Tenthredo notha Klug, 1814, a sawfly

We present a genome assembly from an individual female Tenthredo notha (Arthropoda; Insecta; Hymenoptera; Tenthredinidae). The genome sequence is 253 megabases in span. Most of the assembly (99.91%) is scaffolded into 20 chromosomal pseudomolecules. The mitochondrial genome was also assembled and is 19.8 kilobases in length. Gene annotation of this assembly on Ensembl has identified 10,235 protein coding genes.

The genome sequence of the common green furrow bee, Lasioglossum morio (Fabricius, 1793)

We present a genome assembly from an individual male Lasioglossum morio (the common green furrow bee; Arthropoda; Insecta; Hymenoptera; Halictidae). The genome sequence is 547 megabases in span. Over half of the assembly (55.79%) is scaffolded into 12 chromosomal pseudomolecules. The mitochondrial genome was also assembled, and is 16.8 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,460 protein coding genes.

The genome sequence of the Orange-tailed Clearwing, Synanthedon andrenaeformis (Laspeyres, 1801)

We present a genome assembly from an individual male Synanthedon andrenaeformis (the Orange-tailed Clearwing; Arthropoda; Insecta; Lepidoptera; Sesiidae). The genome sequence is 348.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 16.65 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,867 protein coding genes.

Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy

Allergen-specific venom immunotherapy (VIT) is a well-established therapy for Hymenoptera venom allergy (HVA). However, the precise mechanism underlying its clinical effect remains uncertain. Our study aimed to identify the molecular mechanisms associated with VIT efficiency. We prospectively included 19 patients with HVA undergoing VIT (sampled before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a sting challenge) and 9 healthy controls. RNA sequencing of whole blood was performed on an Illumina sequencing platform. Longitudinal transcriptomic profiling revealed the importance of the inhibition of the NFκB pathway and the downregulation of DUX4 transcripts for the early protection and induction of tolerance after finishing VIT. Furthermore, successful treatment was associated with inhibiting Th2, Th17, and macrophage alternative signalling pathways in synergy with the inhibition of the PPAR pathway and further silencing of the Th2 response. The immune system became activated when reaching the maintenance dose and was suppressed after finishing VIT. Finally, successful VIT restores the immune system's balance to a state similar to that of healthy individuals. Our results underline the important role of the inhibition of four pathways in the clinical effect of VIT: Th2, Th17, NFκB, and macrophage signalling. Two biomarkers specific for successful VIT, regardless of the time of sampling, were C4BPA and RPS10-NUDT3 and should be further tested as potential biomarkers.

The genome sequence of the giant clam, Tridacna gigas (Linnaeus, 1758)

We present a chromosomal-level genome assembly from an individual Tridacna gigas (the giant clam; Mollusca; Bivalvia; Veneroida; Cardiidae). The genome sequence is 1,175.9 megabases in span. Most of the assembly is scaffolded into 17 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 25.34 kilobases in length. Gene annotation of this assembly on Ensembl identified 18,177 protein coding genes.

The genome sequence of the 6-spot burnet, Zygaena filipendulae (Linnaeus, 1758)

We present a genome assembly from an individual female Zygaena filipendulae (6-spot burnet; Arthropoda; Insecta; Lepidoptera; Zygaenidae). The genome sequence is 365.9 megabases in span. The majority of the assembly (99.99%) is scaffolded into 31 chromosomal pseudomolecules, with the W and Z sex chromosomes assembled. The complete mitochondrial genome was also assembled and is 15.6 kilobases in length. Gene annotation of this assembly on Ensembl has identified 12,493 protein coding genes.

Biosurfer for systematic tracking of regulatory mechanisms leading to protein isoform diversity

Long-read RNA sequencing has shed light on transcriptomic complexity, but questions remain about the functionality of downstream protein products. We introduce Biosurfer, a computational approach for comparing protein isoforms, while systematically tracking the transcriptional, splicing, and translational variations that underlie differences in the sequences of the protein products. Using Biosurfer, we analyzed the differences in 32,799 pairs of GENCODE annotated protein isoforms, finding a majority (70%) of variable N-termini are due to the alternative transcription start sites, while only 9% arise from 5' UTR alternative splicing. Biosurfer's detailed tracking of nucleotide-to-residue relationships helped reveal an uncommonly tracked source of single amino acid residue changes arising from the codon splits at junctions. For 17% of internal sequence changes, such split codon patterns lead to single residue differences, termed "ragged codons". Of variable C-termini, 72% involve splice- or intron retention-induced reading frameshifts. We found an unusual pattern of reading frame changes, in which the first frameshift is closely followed by a distinct second frameshift that restores the original frame, which we term a "snapback" frameshift. We analyzed long read RNA-seq-predicted proteome of a human cell line and found similar trends as compared to our GENCODE analysis, with the exception of a higher proportion of isoforms predicted to undergo nonsense-mediated decay. Biosurfer's comprehensive characterization of long-read RNA-seq datasets should accelerate insights of the functional role of protein isoforms, providing mechanistic explanation of the origins of the proteomic diversity driven by the alternative splicing. Biosurfer is available as a Python package at https://github.com/sheynkman-lab/biosurfer.

Follicular Fluid-Derived Extracellular Vesicles Influence on In Vitro Maturation of Equine Oocyte: Impact on Cumulus Cell Viability, Expansion and Transcriptome

Cumulus cell (CC) expansion is pivotal for oocyte maturation, during which CCs release factors that initiate paracrine signaling within the follicular fluid (FF). The FF is abundant in extracellular vesicles (EVs) that facilitate intercellular communication. Although bovine and murine EVs can control cumulus expansion, these effects have not been observed in equines. This study aimed to assess the impact of FF-derived EVs (ffEVs) on equine CC expansion, viability, and transcriptome. Cumulus-oocyte complexes (COCs) that underwent in vitro maturation (IVM) in the presence (200 µg protein/mL) or absence (control) of ffEVs were assessed for cumulus expansion and viability. CCs were isolated after 12 h of IVM, followed by RNA extraction, cDNA library generation, and subsequent transcriptome analysis using next-generation sequencing. Confocal microscopy images illustrated the internalization of labeled ffEVs by CCs. Supplementation with ffEVs significantly enhanced cumulus expansion in both compacted (Cp, p < 0.0001) and expanded (Ex, p < 0.05) COCs, while viability increased in Cp groups (p < 0.01), but decreased in Ex groups (p < 0.05), compared to the controls. Although transcriptome analysis revealed a subtle effect on CC RNA profiles, differentially expressed genes encompassed processes (e.g., MAPK and Wnt signaling) potentially crucial for cumulus properties and, consequently, oocyte maturation.

The genome sequence of a heart cockle, Fragum fragum (Linnaeus, 1758)

We present a genome assembly from an individual specimen of Fragum fragum (a heart cockle; Mollusca; Bivalvia; Veneroida; Cardiidae). The genome sequence is 1,153.1 megabases in span. Most of the assembly is scaffolded into 19 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 22.36 kilobases in length. Gene annotation of this assembly on Ensembl identified 17,262 protein coding genes.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

The genome sequence of the Elephant Hawk-moth, Deilephila elpenor (Linnaeus, 1758)

We present a genome assembly from an individual female Deilephila elpenor (the Elephant Hawk-moth; Arthropoda; Insecta; Lepidoptera; Sphingidae). The genome sequence is 414.1 megabases in span. Most of the assembly is scaffolded into 30 chromosomal pseudomolecules, including the Z and W sex chromosomes. The mitochondrial genome has also been assembled and is 15.37 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,748 protein coding genes.

The genome sequence of the wood-carving leafcutter bee, Megachile ligniseca (Kirby, 1802)

We present a genome assembly from an individual female Megachile ligniseca (the wood-carving leafcutter bee; Arthropoda; Insecta; Hymenoptera; Megachilidae). The genome sequence is 290.0 megabases in span. Most of the assembly is scaffolded into 16 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 23.71 kilobases in length. Gene annotation of this assembly on Ensembl 11,722 protein coding genes.

The genome sequence of the big-headed mining bee, Andrena bucephala (Stephens, 1846)

We present a genome assembly from an individual female Andrena bucephala (the Big-headed Mining Bee; Arthropoda; Insecta; Hymenoptera; Andrenidae). The genome sequence is 379.8 megabases in span. Most of the assembly is scaffolded into 5 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 19.57 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,022 protein coding genes.

The genome sequence of the White-pinion Spotted, Lomographa bimaculata (Fabricius, 1775)

We present a genome assembly from an individual male Lomographa bimaculata (the White-pinion Spotted; Arthropoda; Insecta; Lepidoptera; Geometridae). The genome sequence is 554.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 16.66 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,749 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. 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 editing is a common attribute of relapsed T cell acute lymphoblastic leukemia (T-ALL) regardless of molecular subtype. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL patient-derived xenograft models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA to avoid detection by the innate immune sensor melanoma differentiation-associated protein 5 (MDA5). Moreover, we uncover that the cell-intrinsic level of MDA5 dictates the dependency on the 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 an effective therapeutic strategy for eliminating T-ALL LICs.

The genome sequence of the Mournful Wasp, Pemphredon lugubris (Fabricius, 1793)

We present a genome assembly from an individual male Pemphredon lugubris (the Mournful Wasp; Arthropoda; Insecta; Hymenoptera; Crabronidae). The genome sequence is 328.1 megabases in span. Most of the assembly is scaffolded into 5 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.88 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,335 protein coding genes.

The genome sequence of the Lobe-spurred Furrow Bee, Lasioglossum pauxillum (Schenck, 1853)

We present a genome assembly from an individual female Lasioglossum pauxillum (the Lobe-spurred Furrow Bee; Arthropoda; Insecta; Hymenoptera; Halictidae). The genome sequence is 432.0 megabases in span. Most of the assembly is scaffolded into 9 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 27.71 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,353 protein coding genes.

The genome sequence of the Hebrew Character, Orthosia gothica (Linnaeus, 1758)

We present a genome assembly from an individual female Orthosia gothica (the Hebrew character; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 1,065.1 megabases in span. Most of the assembly is scaffolded into 37 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.38 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,691 protein coding genes.

The genome sequence of the Mottled Pug, Eupithecia exiguata (Hübner, 1813)

We present a genome assembly from an individual male Eupithecia exiguata (the Mottled Pug; Arthropoda; Insecta; Lepidoptera; Geometridae). The genome sequence is 372.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 16.39 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,194 protein coding genes.

The genome sequence of the black-footed limpet, Patella depressa (Pennant, 1777)

We present a genome assembly from an individual Patella depressa (the black-footed limpet; Mollusca; Gastropoda; Patellogastropoda; Patellidae). The genome sequence is 683.7 megabases in span. Most of the assembly is scaffolded into 9 chromosomal pseudomolecules. Gene annotation of this assembly on Ensembl identified 20,502 protein coding genes.

The genome sequence of the Emperor moth, Saturnia pavonia (Linnaeus, 1758)

We present a genome assembly from an individual male Saturnia pavonia (the Emperor moth; Arthropoda; Insecta; Lepidoptera; Saturniidae). The genome sequence is 489.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 15.29 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,903 protein coding genes.

The genome sequence of the White-point, Mythimna albipuncta (Denis & Schiffermüller, 1775)

We present a genome assembly from an individual male Mythimna albipuncta (the White-point; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 698.6 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.38 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,679 protein coding genes.

The genome sequence of the spotted cranefly, Nephrotoma appendiculata (Pierre, 1919)

We present a genome assembly from an individual male Nephrotoma appendiculata (the spotted cranefly; Arthropoda; Insecta; Diptera; Tipulidae). The genome sequence is 1,138.0 megabases in span. Most of the assembly is scaffolded into 4 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 17.42 kilobases in length. Gene annotation of this assembly on Ensembl identified 17,753 protein coding genes.

The genome sequence of the Stripe-backed Dasysyrphus, Dasysyrphus albostriatus (Fallén, 1817)

We present a genome assembly from an individual female Dasysyrphus albostriatus (the Stripe-backed Dasysyrphus; Arthropoda; Insecta; Diptera; Syrphidae). The genome sequence is 662.5 megabases in span. Most of the assembly is scaffolded into 5 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 17.55 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,259 protein coding genes.

A common allele increases endometrial Wnt4 expression, with antagonistic implications for pregnancy, reproductive cancers, and endometriosis

The common human SNP rs3820282 is associated with multiple phenotypes including gestational length and likelihood of endometriosis and cancer, presenting a paradigmatic pleiotropic variant. Deleterious pleiotropic mutations cause the co-occurrence of disorders either within individuals, or across population. When adverse and advantageous effects are combined, pleiotropy can maintain high population frequencies of deleterious alleles. To reveal the causal molecular mechanisms of this pleiotropic SNP, we introduced this substitution into the mouse genome by CRISPR/Cas 9. Previous work showed that rs3820282 introduces a high-affinity estrogen receptor alpha-binding site at the Wnt4 locus. Here, we show that this mutation upregulates Wnt4 transcription in endometrial stroma, following the preovulatory estrogen peak. Effects on uterine transcription include downregulation of epithelial proliferation and induction of progesterone-regulated pro-implantation genes. We propose that these changes increase uterine permissiveness to embryo invasion, whereas they decrease resistance to invasion by cancer and endometriotic foci in other estrogen-responsive tissues.

Human whole-exome genotype data for Alzheimer's disease

The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.

SpliceProt 2.0: A Sequence Repository of Human, Mouse, and Rat Proteoforms

SpliceProt 2.0 is a public proteogenomics database that aims to list the sequence of known proteins and potential new proteoforms in human, mouse, and rat proteomes. This updated repository provides an even broader range of computationally translated proteins and serves, for example, to aid with proteomic validation of splice variants absent from the reference UniProtKB/SwissProt database. We demonstrate the value of SpliceProt 2.0 to predict orthologous proteins between humans and murines based on transcript reconstruction, sequence annotation and detection at the transcriptome and proteome levels. In this release, the annotation data used in the reconstruction of transcripts based on the methodology of ternary matrices were acquired from new databases such as Ensembl, UniProt, and APPRIS. Another innovation implemented in the pipeline is the exclusion of transcripts predicted to be susceptible to degradation through the NMD pathway. Taken together, our repository and its applications represent a valuable resource for the proteogenomics community.

Review of Computational Methods and Database Sources for Predicting the Effects of Coding Frameshift Small Insertion and Deletion Variations

Genetic variations (including substitutions, insertions, and deletions) exert a profound influence on DNA sequences. These variations are systematically classified as synonymous, nonsynonymous, and nonsense, each manifesting distinct effects on proteins. The implementation of high-throughput sequencing has significantly augmented our comprehension of the intricate interplay between gene variations and protein structure and function, as well as their ramifications in the context of diseases. Frameshift variations, particularly small insertions and deletions (indels), disrupt protein coding and are instrumental in disease pathogenesis. This review presents a succinct review of computational methods, databases, current challenges, and future directions in predicting the consequences of coding frameshift small indels variations. We analyzed the predictive efficacy, reliability, and utilization of computational methods and variant account, reliability, and utilization of database. Besides, we also compared the prediction methodologies on GOF/LOF pathogenic variation data. Addressing the challenges pertaining to prediction accuracy and cross-species generalizability, nascent technologies such as AI and deep learning harbor immense potential to enhance predictive capabilities. The importance of interdisciplinary research and collaboration cannot be overstated for devising effective diagnosis, treatment, and prevention strategies concerning diseases associated with coding frameshift indels variations.

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology

BACKGROUND

The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.

RESULTS

A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.

CONCLUSIONS

These validated results show significant improvement over current bovine genome annotations.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs

BACKGROUND

Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs.

RESULTS

The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes.

CONCLUSIONS

After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

The genome sequence of the white-footed hoverfly, Platycheirus albimanus (Fabricius, 1781)

We present a genome assembly from an individual female Platycheirus albimanus (the white-footed hoverfly; Arthropoda; Insecta; Diptera; Syrphidae). The genome sequence is 677.8 megabases in span. Most of the assembly is scaffolded into 4 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 18.17 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,568 protein coding genes.

A CNN based m5c RNA methylation predictor

Post-transcriptional modifications of RNA play a key role in performing a variety of biological processes, such as stability and immune tolerance, RNA splicing, protein translation and RNA degradation. One of these RNA modifications is m5c which participates in various cellular functions like RNA structural stability and translation efficiency, got popularity among biologists. By applying biological experiments to detect RNA m5c methylation sites would require much more efforts, time and money. Most of the researchers are using pre-processed RNA sequences of 41 nucleotides where the methylated cytosine is in the center. Therefore, it is possible that some of the information around these motif may have lost. The conventional methods are unable to process the RNA sequence directly due to high dimensionality and thus need optimized techniques for better features extraction. To handle the above challenges the goal of this study is to employ an end-to-end, 1D CNN based model to classify and interpret m5c methylated data sites. Moreover, our aim is to analyze the sequence in its full length where the methylated cytosine may not be in the center. The evaluation of the proposed architecture showed a promising results by outperforming state-of-the-art techniques in terms of sensitivity and accuracy. Our model achieve 96.70% sensitivity and 96.21% accuracy for 41 nucleotides sequences while 96.10% accuracy for full length sequences.

Hepatocyte CHRNA4 mediates the MASH-promotive effects of immune cell-produced acetylcholine and smoking exposure in mice and humans

Metabolic dysfunction-associated steatohepatitis (MASH) is a leading risk factor for liver cirrhosis and hepatocellular carcinoma. Here, we report that CHRNA4, a subunit of nicotinic acetylcholine receptors (nAChRs), is an accelerator of MASH progression. CHRNA4 also mediates the MASH-promotive effects induced by smoking. Chrna4 was expressed specifically in hepatocytes and exhibited increased levels in mice and patients with MASH. Elevated CHRNA4 levels were positively correlated with MASH severity. We further revealed that during MASH development, acetylcholine released from immune cells or nicotine derived from smoking functioned as an agonist to activate hepatocyte-intrinsic CHRNA4, inducing calcium influx and activation of inflammatory signaling. The communication between immune cells and hepatocytes via the acetylcholine-CHRNA4 axis led to the production of a variety of cytokines, eliciting inflammation in liver and promoting the pathogenesis of MASH. Genetic and pharmacological inhibition of CHRNA4 protected mice from diet-induced MASH. Targeting CHRNA4 might be a promising strategy for MASH therapeutics.

Whole-Exome Sequencing Among School-Aged Children With High Myopia

Importance

High myopia (HM) is one of the leading causes of visual impairment worldwide. Genetic factors are known to play an important role in the development of HM.

Objective

To identify risk variants in a large HM cohort and to examine the implications of genetic testing of schoolchildren with HM.

Design, Setting, and Participants

This cohort study retrospectively reviewed whole-exome sequencing (WES) results in 6215 schoolchildren with HM who underwent genetic testing between September 2019 and July 2020 in Wenzhou City, China. HM is defined as a spherical equivalent refraction (SER) of -6.00 diopters (D) or less. The study setting was a genetic testing laboratory and a multicenter school census. Data were analyzed from July 2021 to June 2022.

Main Outcomes and Measures

The frequency and distribution of positive germline variants, the percentage of individuals with HM in both eyes, and subsequent variant yield for common high myopia (CHM; -8.00 D ≤ SER ≤ -6.00 D), ultra myopia (UM; -10.00 D ≤ SER < -8.00 D), and extreme myopia (EM; SER < -10.00 D).

Results

Of the 6215 schoolchildren with HM, 3278 (52.74%) were male. Their mean (SD) age was 14.87 (2.02) years, including 355 students in primary school, 1970 in junior high school, and 3890 in senior high school. The mean (SD) SER was -7.51 (-1.36) D for the right eye and -7.46 (-1.34) D for the left eye. Among schoolchildren with HM, genetic testing yielded 271 potential pathogenic variants in 75 HM candidate genes in 964 diagnoses (15.52%). A total of 36 known variants were found in 490 HM participants (7.88%) and 235 protein-truncating variants (PTVs) in 506 participants (8.14%). Involved variant yield was significantly positively associated with SER (Cochran-Armitage test for trend Z = 2.5492; P = .01), which ranged from 7.66% in the CHM group, 8.70% in the UM group, to 11.90% in the EM group. We also found that primary school students with EM had the highest variant yield of PTVs (8 of 35 students [22.86%]), which was 1.77 and 4.78 times that of the UM and CHM, respectively.

Conclusions and Relevance

In this cohort study of WES for HM, several potential pathogenic variants were identified in a substantial number of schoolchildren with HM. The high variation frequency in younger students with EM can provide clues for genetic screening and clinical examinations of HM to promote long-term follow-up assessment.

Circadian clock crosstalks with autism

BACKGROUND

The mechanism underlying autism spectrum disorder (ASD) remains incompletely understood, but researchers have identified over a thousand genes involved in complex interactions within the brain, nervous, and immune systems, particularly during the mechanism of brain development. Various contributory environmental effects including circadian rhythm have also been studied in ASD. Thus, capturing the global picture of the ASD-clock network in combined form is critical.

METHODS

We reconstructed the protein-protein interaction network of ASD and circadian rhythm to understand the connection between autism and the circadian clock. A graph theoretical study is undertaken to evaluate whether the network attributes are biologically realistic. The gene ontology enrichment analyses provide information about the most important biological processes.

RESULTS

This study takes a fresh look at metabolic mechanisms and the identification of potential key proteins/pathways (ribosome biogenesis, oxidative stress, insulin/IGF pathway, Wnt pathway, and mTOR pathway), as well as the effects of specific conditions (such as maternal stress or disruption of circadian rhythm) on the development of ASD due to environmental factors.

CONCLUSION

Understanding the relationship between circadian rhythm and ASD provides insight into the involvement of these essential pathways in the pathogenesis/etiology of ASD, as well as potential early intervention options and chronotherapeutic strategies for treating or preventing the neurodevelopmental disorder.