Apollo: Democratizing genome annotation

Diversified quantity, gene structure, and expression profile of OPR gene family of A. annua

Artemisia annua, the source of artemisinin production, is a traditional herb used for treating malaria for thousand years. The genetic background is of high heterozygosity and traits (plant height, biomass, artemisinin content, etc.) are diverse across different germplasms. Unraveling the key genes associated with growth and secondary metabolism is essential for the efficient production of artemisinin. The 12-oxo-phytodienoic acid reductase (OPR) genes, crucial for plant growth and development and stress resistance, remain unexplored in A. annua. In this study, nine OPR genes (named as AaOPR1 to AaOPR9) were identified in A. annua, including two pairs of genes formed from recent tandem duplications. The number of OPRs varied among different haplotype genomes, and each OPR gene exhibiting distinct expression pattern. Moreover, the OPR family displayed evolutionarily activity with significant variations in numbers and gene structures observed across different plant species. Widespread gene duplication of OPRs, observed in the majority of analyzed plant genomes, brought evolutionary potential. DBR2, a member of AaOPRs involved in artemisinin biosynthesis, had two copies (AaOPR1/DBR2.1 and AaOPR2/DBR2.2) with different expression patterns, one of which was a recently generated copy with a significant 7-amino acids truncation. Heterologous protein expression and functional characterization of the two copies of DBR2 yielded multiple isomers with identical molecular weights but different arrangements, indicating neofunctionalization of the newly generated copy. The polymorphism within the OPR gene family merely scratches the surface of the genetic diversity in A. annua, and further investigation of genetic features is needed for the screening of elite germplasm resources.

Collaborative metabolic curation of an emerging model marine bacterium, Alteromonas macleodii ATCC 27126

Inferring the metabolic capabilities of an organism from its genome is a challenging process, relying on computationally-derived or manually curated metabolic networks. Manual curation can correct mistakes in the draft network and add missing reactions based on the literature, but requires significant expertise and is often the bottleneck for high-quality metabolic reconstructions. Here, we present a synopsis of a community curation workshop for the model marine bacterium Alteromonas macleodii ATCC 27126 and its genome database in BioCyc, focusing on pathways for utilizing organic carbon and nitrogen sources. Due to the scarcity of biochemical information or gene knock-outs, the curation process relied primarily on published growth phenotypes and bioinformatic analyses, including comparisons with related Alteromonas strains. We report full pathways for the utilization of the algal polysaccharides alginate and pectin in contrast to inconclusive evidence for one-carbon metabolism and mixed acid fermentation, in accordance with the lack of growth on methanol and formate. Pathways for amino acid degradation are ubiquitous across Alteromonas macleodii strains, yet enzymes in the pathways for the degradation of threonine, tryptophan and tyrosine were not identified. Nucleotide degradation pathways are also partial in ATCC 27126. We postulate that demonstrated growth on nitrate as sole nitrogen source proceeds via a nitrate reductase pathway that is a hybrid of known pathways. Our evidence highlights the value of joint and interactive curation efforts, but also shows major knowledge gaps regarding Alteromonas metabolism. The manually-curated metabolic reconstruction is available as a "Tier-2" database on BioCyc.

A Hitchhiker's Guide to long-read genomic analysis

Over the past decade, long-read sequencing has evolved into a pivotal technology for uncovering the hidden and complex regions of the genome. Significant cost efficiency, scalability, and accuracy advancements have driven this evolution. Concurrently, novel analytical methods have emerged to harness the full potential of long reads. These advancements have enabled milestones such as the first fully completed human genome, enhanced identification and understanding of complex genomic variants, and deeper insights into the interplay between epigenetics and genomic variation. This mini-review provides a comprehensive overview of the latest developments in long-read DNA sequencing analysis, encompassing reference-based and de novo assembly approaches. We explore the entire workflow, from initial data processing to variant calling and annotation, focusing on how these methods improve our ability to interpret a wide array of genomic variants. Additionally, we discuss the current challenges, limitations, and future directions in the field, offering a detailed examination of the state-of-the-art bioinformatics methods for long-read sequencing.

Genomic Insights into Post-Domestication Expansion and Selection of Body Size in Ponies

Horse domestication revolutionizes human civilization by transforming transportation, agriculture, and warfare patterns. Despite extensive studies on modern domestic horse origins, the intricate demographic history and genetic signatures underlying pony size remain unexplored. Here, a high-quality genome assembly of the Chinese Debao pony is presented, and 452 qualified individuals from 64 horse breeds worldwide are extensively analyzed. The authors' results reveal the conservation of ancient components in East Asian horses and close relationships between Asian horses and Western pony lineages. Genetic analyses suggest an Asian paternal origin for European pony breeds. These pony-sized horses share close genetic affinities, potentially attributed to their early expansion and adaptation to local environments. In addition, promising cis-regulatory elements influencing horse withers height by regulating genes such as RFLNA and FOXO1 are identified. Overall, this study provides insightful perspectives on the dispersal history and genetic determinants underlying body size in ponies, offering broader implications for horse population management and improvement.

Complete mitochondrial genome of Paraisaria gracilioides (Hypocreales: Ophiocordycipitaceae) and phylogenetic analysis

Paraisaria gracilioides (Kobayasi) Luangsa-ard, Mongkolsamrit & Samson, originally described as Ophiocordyceps gracilioides, is a member of the Ophiocordycipitaceae, within the order Hypocreales. This study, presented the mitochondrial genome of P. gracilioides structured as a circular molecule spanning 181,140 bp, was larger than those of most Ophiocordycipitaceae species. Despite this size variation, encoding 49 genes, including 15 PCGs, 26 tRNAs, two rRNAs, and six ORFs. Phylogenetic analysis based on nucleotide sequences of 14 PCGs revealed that P. gracilioides clustered with Paraisaria gracilis, forming a subgroup within Ophiocordyceps species. The results contribute to refining the taxonomic framework of Ophiocordycipitaceae.

Genome-wide identification and expression analysis of the ADH gene family in Artemisia annua L. under UV-B stress

ADHs are key genes that catalyze the interconversion between alcohols and aldehydes, which play crucial roles in plant adaptation to a range of abiotic stresses. However, the characterization and evolutionary pathways of ADH genes in the antimalarial plant Artemisia annua are still unclear. This study identified 49 ADH genes in A. annua and conducted a detailed analysis of their structural features, conserved motifs, and duplication types, revealing that tandem and dispersed duplications are the primary mechanisms of gene expansion. Evolutionary analysis of ADH genes between A. annua (AanADH) and A. argyi (AarADH) revealed dynamic changes, with 35 genes identified deriving from their most recent common ancestor in both species. ADH1, crucial for artemisinin production, had two copies in both species, expanding via dispersed duplication in A. annua but whole-genome duplication in A. argyi. CREs and WGCNA analysis suggested that AanADH genes may be regulated by UV-B stress. Following short-term UV-B treatment, 16 DEGs were identified, including ADH1 (AanADH6 and AanADH7), and these genes were significantly downregulated after two hours treatment (UV2h) and upregulated after four hours treatment (UV4h). The expression changes of these genes were further confirmed by GO enrichment analysis and qRT-PCR experiments. Overall, this study comprehensively characterized the ADH gene family in A. annua and systematically identified AanADH genes that were responsive to UV-B stress, providing a foundation for further research on their roles in abiotic stress responses.

Complete genome sequence of Bacillus subtilis bacteriophage Adastra

Adastra is a lytic bacteriophage that infects Bacillus subtilis. Here, we report the sequencing and annotation of the 136,306-bp genome of Adastra and its similarity to other myophages in the SPO1 family.

Long-term human influence on the demography and genetic diversity of the hyperdominant Bertholletia excelsa in the Amazon Basin

The Amazon rainforest is characterized by a limited number of hyperdominant trees that play an oversized role in its ecosystems, nutrient cycle, and rainfall production. Some of these, such as the Brazil nut, appear to have been intensively exploited and dispersed by Indigenous populations since their earliest arrival in this part of South America around 13,000 years ago. However, the genetic diversity-and geographic structure-of these species remains poorly understood, as does their exact relationship with past human land use. We use a new genome assembly for Brazil nut to analyze 270 individuals sampled at areas with varying intensities of archaeological evidence. We demonstrate that overall low genetic diversity, with a notable decrease since the Late Pleistocene, is accompanied by significant geographic structure, where evidence for improved gene flow and regeneration by long-term traditional human management is linked to increased genetic diversity. We argue that historical perspectives on the genetic diversity of key tree species, such as the Brazil nut, can support the development of more active management strategies today.

Chromosomal genome assembly resolves drug resistance loci in the parasitic nematode Teladorsagia circumcincta

The parasitic nematode Teladorsagia circumcincta is one of the most important pathogens of sheep and goats in temperate climates worldwide and can rapidly evolve resistance to drugs used to control it. To understand the genetics of drug resistance, we have generated a highly contiguous genome assembly for the UK T. circumcincta isolate, MTci2. Assembly using PacBio long-reads and Hi-C long-molecule scaffolding together with manual curation resulted in a 573 Mb assembly (N50 = 84 Mb, total scaffolds = 1,286) with five autosomal and one sex-linked chromosomal-scale scaffolds consistent with its karyotype. The genome resource was further improved via annotation of 22,948 genes, with manual curation of over 3,200 of these, resulting in a robust and near complete resource (96.3% complete protein BUSCOs) to support basic and applied research on this important veterinary pathogen. Genome-wide analyses of drug resistance, combining evidence from three distinct experiments, identified selection around known candidate genes for benzimidazole, levamisole and ivermectin resistance, as well as novel regions associated with ivermectin and moxidectin resistance. These insights into contemporary and historic genetic selection further emphasise the importance of contiguous genome assemblies in interpreting genome-wide genetic variation associated with drug resistance and identifying key loci to prioritise in developing diagnostic markers of anthelmintic resistance to support parasite control.

Genome sequences of four Ixodes species expands understanding of tick evolution

BACKGROUND

Ticks, hematophagous Acari, pose a significant threat by transmitting various pathogens to their vertebrate hosts during feeding. Despite advances in tick genomics, high-quality genomes were lacking until recently, particularly in the genus Ixodes, which includes the main vectors of Lyme disease.

RESULTS

Here, we present the genome sequences of four tick species, derived from a single female individual, with a particular focus on the European species Ixodes ricinus, achieving a chromosome-level assembly. Additionally, draft assemblies were generated for the three other Ixodes species, I. persulcatus, I. pacificus, and I. hexagonus. The quality of the four genomes and extensive annotation of several important gene families have allowed us to study the evolution of gene repertoires at the level of the genus Ixodes and of the tick group. We have determined gene families that have undergone major amplifications during the evolution of ticks, while an expression atlas obtained for I. ricinus reveals striking patterns of specialization both between and within gene families. Notably, several gene family amplifications are associated with a proliferation of single-exon genes-most strikingly for fatty acid elongases and sulfotransferases.

CONCLUSIONS

The integration of our data with existing genomes establishes a solid framework for the study of gene evolution, improving our understanding of tick biology. In addition, our work lays the foundations for applied research and innovative control targeting these organisms.

Features and evolutionary adaptations of the mitochondrial genome of Garuga forrestii W. W. Sm

Introduction

Garuga forrestii W. W. Sm. is a tree species of the Burseraceae family, endemic to China, found in hot/warm-dry valleys. This species plays a crucial role in maintaining biodiversity in these ecosystems.

Methods

We performed de novo assembly of the Garuga forrestii mitochondrial genome using PMAT (v.1.5.4), resulting in a typical circular molecule of 606,853 bp. The genome consists of 31 tRNA genes, 3 rRNA genes, 35 protein-coding genes, and 1 pseudogene. The study also investigates RNA editing sites and evolutionary patterns.

Results

The mitochondrial genome exhibits a low proportion of repetitive sequences (3.30%), suggesting a highly conserved structure. A high copy number of the trnM-CAT gene (4 copies) is noted, which may contribute to genomic rearrangement and adaptive evolution. Among the 476 RNA editing sites, hydrophilic-hydrophobic and hydrophobic-hydrophobic editing events are most common, accounting for 77.10%. Negative selection predominates among most genes (Ka/Ks < 1), while a few genes (e.g., matR, nad3, rps1, rps12, and rps4) show signs of positive selection (Ka/Ks > 1), potentially conferring evolutionary advantages. Additionally, a significant A/T bias is observed at the third codon position. Phylogenomic analysis supports the APG IV classification, with no evidence of horizontal gene transfer.

Discussion

This mitochondrial genome offers valuable insights into the adaptive mechanisms and evolutionary processes of Garuga forrestii. It enhances our understanding of the species' biogeography in tropical Southeast Asia and Southwest China, providing key information on the evolutionary history of this genus.

Comprehensive analysis of the first complete mitogenome and plastome of a traditional Chinese medicine Viola diffusa

BACKGROUND

Viola diffusa is used in the formulation of various Traditional Chinese Medicines (TCMs), including antiviral, antimicrobial, antitussive, and anti-inflammatory drugs, due to its richness in flavonoids and triterpenoids. The biosynthesis of these compounds is largely mediated by cytochrome P450 enzymes, which are primarily located in the membranes of mitochondria and the endoplasmic reticulum.

RESULTS

This study presents the complete assembly of the mitogenome and plastome of Viola diffusa. The circular mitogenome spans 474,721 bp with a GC content of 44.17% and encodes 36 unique protein-coding genes, 21 tRNA, and 3 rRNA. Except for the RSCU values of 1 observed for the start codon (AUG) and tryptophan (UGG), the mitochondrial protein-coding genes exhibited a codon usage bias, with most estimates deviating from 1, similar to patterns observed in closely related species. Analysis of repetitive sequences in the mitogenome demonstrated potential homologous recombination mediated by these repeats. Sequence transfer analysis revealed 24 homologous sequences shared between the mitogenome and plastome, including nine full-length genes. Collinearity was observed among Viola diffusa species within the other members of Malpighiales order, indicated by the presence of homologous fragments. The length and arrangement of collinear blocks varied, and the mitogenome exhibited a high frequency of gene rearrangement.

CONCLUSIONS

We present the first complete assembly of the mitogenome and plastome of Viola diffusa, highlighting its implications for pharmacological, evolutionary, and taxonomic studies. Our research underscores the multifaceted importance of comprehensive mitogenome analysis.

VPS13B is localized at the interface between Golgi cisternae and is a functional partner of FAM177A1

Mutations in VPS13B, a member of a protein family implicated in bulk lipid transport between adjacent membranes, cause Cohen syndrome. VPS13B is known to be concentrated in the Golgi complex, but its precise location within this organelle and thus the site(s) where it achieves lipid transport remains unclear. Here, we show that VPS13B is localized at the interface between proximal and distal Golgi subcompartments and that Golgi complex reformation after Brefeldin A (BFA)-induced disruption is delayed in VPS13B KO cells. This delay is phenocopied by the loss of FAM177A1, a Golgi complex protein of unknown function reported to be a VPS13B interactor and whose mutations also result in a developmental disorder. In zebrafish, the vps13b ortholog, not previously annotated in this organism, genetically interacts with fam177a1. Collectively, these findings raise the possibility that bulk lipid transport by VPS13B may play a role in the dynamics of Golgi membranes and that VPS13B may be assisted in this function by FAM177A1.

Building a FAIR data ecosystem for incorporating single-cell transcriptomics data into agricultural genome to phenome research

Introduction

The agriculture genomics community has numerous data submission standards available, but the standards for describing and storing single-cell (SC, e.g., scRNA- seq) data are comparatively underdeveloped.

Methods

To bridge this gap, we leveraged recent advancements in human genomics infrastructure, such as the integration of the Human Cell Atlas Data Portal with Terra, a secure, scalable, open-source platform for biomedical researchers to access data, run analysis tools, and collaborate. In parallel, the Single Cell Expression Atlas at EMBL-EBI offers a comprehensive data ingestion portal for high-throughput sequencing datasets, including plants, protists, and animals (including humans). Developing data tools connecting these resources would offer significant advantages to the agricultural genomics community. The FAANG data portal at EMBL-EBI emphasizes delivering rich metadata and highly accurate and reliable annotation of farmed animals but is not computationally linked to either of these resources.

Results

Herein, we describe a pilot-scale project that determines whether the current FAANG metadata standards for livestock can be used to ingest scRNA-seq datasets into Terra in a manner consistent with HCA Data Portal standards. Importantly, rich scRNA-seq metadata can now be brokered through the FAANG data portal using a semi-automated process, thereby avoiding the need for substantial expert curation. We have further extended the functionality of this tool so that validated and ingested SC files within the HCA Data Portal are transferred to Terra for further analysis. In addition, we verified data ingestion into Terra, hosted on Azure, and demonstrated the use of a workflow to analyze the first ingested porcine scRNA-seq dataset. Additionally, we have also developed prototype tools to visualize the output of scRNA-seq analyses on genome browsers to compare gene expression patterns across tissues and cell populations. This JBrowse tool now features distinct tracks, showcasing PBMC scRNA-seq alongside two bulk RNA-seq experiments.

Discussion

We intend to further build upon these existing tools to construct a scientist-friendly data resource and analytical ecosystem based on Findable, Accessible, Interoperable, and Reusable (FAIR) SC principles to facilitate SC-level genomic analysis through data ingestion, storage, retrieval, re-use, visualization, and comparative annotation across agricultural species.

Structural and evolutionary analyses of the mitochondrial genome of Spuriopimpinella brachycarpa

Introduction

Spuriopimpinella brachycarpa (Kom.) Kitag., a member of the Apiaceae family, is a perennial aromatic herb native to Northeast Asia with applications in culinary and traditional medicine. Despite its significance, most studies on S. brachycarpa have primarily focused on its phytochemical properties, with limited insights into its molecular and genomic characteristics.

Methods

This study presents the sequencing and assembly of the mitochondrial genome (mitogenome) of S. brachycarpa using second- and third-generation high-throughput sequencing technologies. Comprehensive analyses were performed on its structural organization, RNA editing sites, relative synonymous codon usage (RSCU), and repeat sequences. Comparative analyses with closely related species were also conducted.

Results

The mitogenome exhibited a multi-branched structure, with a total length of 523,512 bp and a GC content of 43.37%. Annotation revealed 30 unique protein-coding genes, 21 tRNA genes, and three rRNA genes. Comparative analysis indicated that the S. brachycarpa mitogenome contains structural variations but shares collinear features with other Apiaceae species. We identified 618 potential RNA editing sites involving C-to-U conversions and discovered 59 homologous fragments between the mitogenome and plastome, comprising 8.13% of the mitogenome.

Discussion

These results enrich the genomic database of Apiaceae, providing valuable insights into the evolutionary relationships and genetic diversity within the family.

The golden genome annotation of Ganoderma lingzhi reveals a more complex scenario of eukaryotic gene structure and transcription activity

BACKGROUND

It is generally accepted that nuclear genes in eukaryotes are located independently on chromosomes and expressed in a monocistronic manner. However, accumulating evidence suggests a more complex landscape of gene structure and transcription. Ganoderma lingzhi, a model medicinal fungus, currently lacks high-quality genome annotation, hindering genetic studies.

RESULTS

Here, we reported a golden annotation of G. lingzhi, featuring 14,147 high-confidence genes derived from extensive manual corrections. Novel characteristics of gene structure and transcription were identified accordingly. Notably, non-canonical splicing sites accounted for 1.99% of the whole genome, with the predominant types being GC-AG (1.85%), GT-AC (0.05%), and GT-GG (0.04%). 1165 pairs of genes were found to have overlapped transcribed regions, and 92.19% of which showed opposite directions of gene transcription. A total of 5,412,158 genetic variations were identified among 13 G. lingzhi strains, and the manually corrected gene sets resulted in enhanced functional annotation of these variations. More than 60% of G. lingzhi genes were alternatively spliced. In addition, we found that two or more protein-coding genes (PCGs) can be transcribed into a single RNA molecule, referred to as polycistronic genes. In total, 1272 polycistronic genes associated with 2815 PCGs were identified.

CONCLUSIONS

The widespread presence of polycistronic genes in G. lingzhi strongly complements the theory that polycistron is also present in eukaryotic genomes. The extraordinary gene structure and transcriptional activity uncovered through this golden annotation provide implications for the study of genes, genomes, and related studies in G. lingzhi and other eukaryotes.

Deciphering the complex organelle genomes of two Rhododendron species and insights into adaptive evolution patterns in high-altitude

BACKGROUND

The genomes within organelles are crucial for physiological functions such as respiration and photosynthesis and may also contribute to environmental adaptation. However, the limited availability of genetic resources, particularly mitochondrial genomes, poses significant challenges for in-depth investigations.

RESULTS

Here, we explored various assembly methodologies and successfully reconstructed the complex organelle genomes of two Rhododendron species: Rhododendron nivale subsp. boreale and Rhododendron vialii. The mitogenomes of these species exhibit various conformations, as evidenced by long-reads mapping. Notably, only the mitogenome of R. vialii can be depicted as a singular circular molecule. The plastomes of both species conform to the typical quadripartite structure but exhibit elongated inverted repeat (IR) regions. Compared to the high similarity between plastomes, the mitogenomes display more obvious differences in structure, repeat sequences, and codon usage. Based on the analysis of 58 organelle genomes from angiosperms inhabiting various altitudes, we inferred the genetic adaptations associated with high-altitude environments. Phylogenetic analysis revealed partial inconsistencies between plastome- and mitogenome-derived phylogenies. Additionally, evolutionary lineage was determined to exert a greater influence on codon usage than altitude. Importantly, genes such as atp4, atp9, mttB, and clpP exhibited signs of positive selection in several high-altitude species, suggesting a potential link to alpine adaptation.

CONCLUSIONS

We tested the effectiveness of different organelle assembly methods for dealing with complex genomes, while also providing and validating high-quality organelle genomes of two Rhododendron species. Additionally, we hypothesized potential strategies for high-altitude adaptation of organelles. These findings offer a reference for the assembly of complex organelle genomes, while also providing new insights and valuable resources for understanding their adaptive evolution patterns.

RNA-seq and metabolomic analyses of beneficial plant phenol biochemical pathways in red alder

Red alder (Alnus rubra) has highly desirable wood, dye pigment, and (traditional) medicinal properties which have been capitalized on for thousands of years, including by Pacific West Coast Native Americans. A rapidly growing tree species native to North American western coastal and riparian regions, it undergoes symbiosis with actinobacterium Frankia via their nitrogen-fixing root nodules. Red alder's desirable properties are, however, largely attributed to its bioactive plant phenol metabolites, including for plant defense, for its attractive wood and bark coloration, and various beneficial medicinal properties. Integrated transcriptome and metabolome data analyses were carried out using buds, leaves, stems, roots, and root nodules from greenhouse grown red alder saplings with samples collected during different time-points (Spring, Summer, and Fall) of the growing season. Pollen and catkins were collected from field grown mature trees. Overall plant phenol biochemical pathways operative in red alder were determined, with a particular emphasis on potentially identifying candidates for the long unknown gateway entry points to the proanthocyanidin (PA) and ellagitannin metabolic classes, as well as in gaining better understanding of the biochemical basis of diarylheptanoid formation, i.e. that help define red alder's varied medicinal uses, and its extensive wood and dye usage.

Chromosome-scale assembly with improved annotation provides insights into breed-wide genomic structure and diversity in domestic cats

INTRODUCTION

Comprehensive genomic resources offer insights into biological features, including traits/disease-related genetic loci. The current reference genome assembly for the domestic cat (Felis catus), Felis_Catus_9.0 (felCat9), derived from sequences of the Abyssinian cat, may inadequately represent the general cat population, limiting the extent of deducible genetic variations.

OBJECTIVES

The goal was to develop Anicom American Shorthair 1.0 (AnAms1.0), a reference-grade chromosome-scale cat genome assembly.

METHODS

In contrast to prior assemblies relying on Abyssinian cat sequences, AnAms1.0 was constructed from the sequences of more popular American Shorthair breed, which is related to more breeds than the Abyssinian cat. By combining advanced genomics technologies, including PacBio long-read sequencing and Hi-C- and optical mapping data-based sequence scaffolding, we compared AnAms1.0 to existing Felidae genome assemblies (20 scaffolds, scaffolds N50 > 150 Mbp). Homology-based and ab initio gene annotation through Iso-Seq and RNA-Seq was used to identify new coding genes and splice variants.

RESULTS

AnAms1.0 demonstrated superior contiguity and accuracy than existing Felidae genome assemblies. Using AnAms1.0, we identified over 1.5 thousand structural variants and 29 million repetitions compared to felCat9. Additionally, we identified > 1,600 novel protein-coding genes. Notably, olfactory receptor structural variants and cardiomyopathy-related variants were identified.

CONCLUSION

AnAms1.0 facilitates the discovery of novel genes related to normal and disease phenotypes in domestic cats. The analyzed data are publicly accessible on Cats-I (https://cat.annotation.jp/), which we established as a platform for accumulating and sharing genomic resources to discover novel genetic traits and advance veterinary medicine.

Telomere-to-telomere assemblies of cattle and sheep Y-chromosomes uncover divergent structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. Here, we assemble complete and gapless telomere to telomere (T2T) Y chromosomes for these species. We find that the pseudo-autosomal regions are similar in length, but the total chromosome size is substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity is accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 19MYA. The centromeres also differ dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosomes have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Chamaeleo calyptratus (veiled chameleon) chromosome-scale genome assembly and annotation provides insights into the evolution of reptiles and developmental mechanisms

The family Chamaeleonidae comprises 228 species, boasting an extensive geographic spread and an array of evolutionary novelties and adaptations, but a paucity of genetic and molecular analyses. Veiled chameleon (Chamaeleo calyptratus) has emerged as a tractable research organism for the study of squamate early development and evolution. Here we report a chromosomal-level assembly and annotation of the veiled chameleon genome. We note a remarkable chromosomal conservation across squamates, but comparisons to more distant genomes reveal GC peaks correlating with ancestral chromosome fusion events. We subsequently identified the XX/XY region on chromosome 5, confirming environmental-independent sex determination in veiled chameleons. Furthermore, our analysis of the Hox gene family indicates that veiled chameleons possess the most complete set of 41 Hox genes, retained from an amniote ancestor. Lastly, the veiled chameleon genome has retained both ancestral paralogs of the Nodal gene, but is missing Dand5 and several other genes, recently associated with the loss of motile cilia during the establishment of left-right patterning. Thus, a complete veiled chameleon genome provides opportunities for novel insights into the evolution of reptilian genomes and the molecular mechanisms driving phenotypic variation and ecological adaptation.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Differentiation in detoxification gene complements, including neofunctionalization of duplicated cytochrome P450 genes, between lineages of cotton bollworm, Helicoverpa armigera

Here we investigate the evolutionary dynamics of five enzyme superfamilies (CYPs, GSTs, UGTs, CCEs and ABCs) involved in detoxification in Helicoverpa armigera. The reference assembly for an African isolate of the major lineages, H. a. armigera, has 373 genes in the five superfamilies. Most of its CYPs, GSTs, UGTs and CCEs and a few of its ABCs occur in blocks and most of the clustered genes are in subfamilies specifically implicated in detoxification. Most of the genes have orthologues in the reference genome for the Oceania lineage, H. a. conferta. However, clustered orthologues and subfamilies specifically implicated in detoxification show greater sequence divergence and less constraint on non-synonymous differences between the two assemblies than do other members of the five superfamilies. Two duplicated CYPs, which were found in the H. a. armigera but not H. a. conferta reference genome, were also missing in 16 Chinese populations spanning two different lineages of H. a. armigera. The enzyme produced by one of these duplicates has higher activity against esfenvalerate than a previously described chimeric CYP mutant conferring pyrethroid resistance. Various transposable elements were found in the introns of most detoxification genes, generating diverse gene structures. Extensive resequencing data for the Chinese H. a. armigera and H. a. conferta lineages also revealed complex copy number polymorphisms in 17 CCE001s in a cluster also implicated in pyrethroid metabolism, with substantial haplotype differences between all three lineages. Our results suggest that cotton bollworm has a versatile complement of detoxification genes which are evolving in diverse ways across its range.

Genetic factors acting prior to dormancy in sour cherry influence bloom time the following spring

Understanding the process of Prunus species floral development is crucial for developing strategies to manipulate bloom time and prevent crop loss due to climate change. Here, we present a detailed examination of flower development from initiation until bloom for early- and late-blooming sour cherries (Prunus cerasus) from a population segregating for a major bloom time QTL on chromosome 4. Using a new staging system, we show floral buds from early-blooming trees were persistently more advanced than those from late-blooming siblings. A genomic DNA coverage analysis revealed the late-blooming haplotype of this QTL, k, is located on a subgenome originating from the late-blooming P. fruticosa progenitor. Transcriptome analyses identified many genes within this QTL as differentially expressed between early- and late-blooming trees during the vegetative-to-floral transition. From these, we identified candidate genes for the late bloom phenotype, including multiple transcription factors homologous to Reproductive Meristem B3 domain-containing proteins. Additionally, we determined that the basis of k in sour cherry is likely separate from candidate genes found in sweet cherry-suggesting several major regulators of bloom time are located on Prunus chromosome 4.

Genome-wide identification and expression analysis revealed key transcription factors as potential regulators of high-temperature adaptation of Coriolopsis trogii

Transcription factors (TFs) play a crucial role in gene expression, and studying them can lay the foundation for future research on the functional characterization of TFs involved in various biological processes. In this study, we conducted a genome-wide identification and analysis of TFs in the thermotolerant basidiomycete fungus, Coriolopsis trogii. The TF repertoire of C. trogii consisted of 350 TFs, with C2H2 and Zn2C6 being the largest TF families. When the mycelia of C. trogii were cultured on PDA and transferred from 25 to 35 °C, 14 TFs were up-regulated and 14 TFs were down-regulated. By analyzing RNA-seq data from mycelia cultured at different temperatures and under different carbon sources, we identified 22 TFs that were differentially expressed in more than three comparisons. Co-expression analysis revealed that seven differentially expressed TFs, including four Zn2C6s, one Hap4_Hap_bind, one HMG_box, and one Zinc_knuckle, showed significant correlation with 729 targeted genes. Overall, this study provides a comprehensive characterization of the TF family and systematically screens TFs involved in the high-temperature adaptation of C. trogii, laying the groundwork for further research into the specific roles of TFs in the heat tolerance mechanisms of filamentous fungi.

Reciprocal conversion between annual and polycarpic perennial flowering behavior in the Brassicaceae

The development of perennial crops holds great promise for sustainable agriculture and food security. However, the evolution of the transition between perenniality and annuality is poorly understood. Here, using two Brassicaceae species, Crucihimalaya himalaica and Erysimum nevadense, as polycarpic perennial models, we reveal that the transition from polycarpic perennial to biennial and annual flowering behavior is a continuum determined by the dosage of three closely related MADS-box genes. Diversification of the expression patterns, functional strengths, and combinations of these genes endows species with the potential to adopt various life-history strategies. Remarkably, we find that a single gene among these three is sufficient to convert winter-annual or annual Brassicaceae plants into polycarpic perennial flowering plants. Our work delineates a genetic basis for the evolution of diverse life-history strategies in plants and lays the groundwork for the generation of diverse perennial Brassicaceae crops in the future.

mRNA cap-binding protein eIF4E1 is a novel regulator of Toxoplasma gondii latency

The protozoan parasite Toxoplasma gondii causes serious opportunistic disease due to its ability to persist in patients as latent tissue cysts. The molecular mechanisms coordinating conversion between proliferative parasites (tachyzoites) and latent cysts (bradyzoites) are not fully understood. We previously showed that phosphorylation of eIF2α accompanies bradyzoite formation, suggesting that this clinically relevant process involves regulation of mRNA translation. In this study, we investigated the composition and role of eIF4F multi-subunit complexes in translational control. Using CLIPseq, we find that the cap-binding subunit, eIF4E1, localizes to the 5'-end of all tachyzoite mRNAs, many of which show evidence of stemming from heterogeneous transcriptional start sites. We further show that eIF4E1 operates as the predominant cap-binding protein in two distinct eIF4F complexes. Using genetic and pharmacological approaches, we found that eIF4E1 deficiency triggers efficient spontaneous formation of bradyzoites without stress induction. Consistent with this result, we also show that stress-induced bradyzoites exhibit reduced eIF4E1 expression. Overall, our findings establish a novel role for eIF4F in translational control required for parasite latency and microbial persistence.

IMPORTANCE

Toxoplasma gondii is an opportunistic pathogen important to global human and animal health. There are currently no chemotherapies targeting the encysted form of the parasite. Consequently, a better understanding of the mechanisms controlling encystation is required. Here we show that the mRNA cap-binding protein, eIF4E1, regulates the encystation process. Encysted parasites reduce eIF4E1 levels, and depletion of eIF4E1 decreases the translation of ribosome-associated machinery and drives Toxoplasma encystation. Together, these data reveal a new layer of mRNA translational control that regulates parasite encystation and latency.

Evidence-based unification of potato gene models with the UniTato collaborative genome browser

Potato (Solanum tuberosum) is the most popular tuber crop and a model organism. A variety of gene models for potato exist, and despite frequent updates, they are not unified. This hinders the comparison of gene models across versions, limits the ability to reuse experimental data without significant re-analysis, and leads to missing or wrongly annotated genes. Here, we unify the recent potato double monoploid v4 and v6 gene models by developing an automated merging protocol, resulting in a Unified poTato genome model (UniTato). We subsequently established an Apollo genome browser (unitato.nib.si) that enables public access to UniTato and further community-based curation. We demonstrate how the UniTato resource can help resolve problems with missing or misplaced genes and can be used to update or consolidate a wider set of gene models or genome information. The automated protocol, genome annotation files, and a comprehensive translation table are provided at github.com/NIB-SI/unitato.

The desert woodrat (Neotoma lepida) induces a diversity of biotransformation genes in response to creosote bush resin

Liver biotransformation enzymes have long been thought to enable animals to feed on diets rich in xenobiotic compounds. However, despite decades of pharmacological research in humans and rodents, little is known about hepatic gene expression in specialized mammalian herbivores feeding on toxic diets. Leveraging a recently identified population of the desert woodrat (Neotoma lepida) found to be highly tolerant to toxic creosote bush (Larrea tridentata), we explored the expression changes of suites of biotransformation genes in response to diets enriched with varying amounts of creosote resin. Analysis of hepatic RNA-seq data indicated a dose-dependent response to these compounds, including the upregulation of several genes encoding transcription factors and numerous phase I, II, and III biotransformation families. Notably, elevated expression of five biotransformation families - carboxylesterases, cytochromes P450, aldo-keto reductases, epoxide hydrolases, and UDP-glucuronosyltransferases - corresponded to species-specific duplication events in the genome, suggesting that these genes play a prominent role in N. lepida's adaptation to creosote bush. Building on pharmaceutical studies in model rodents, we propose a hypothesis for how the differentially expressed genes are involved in the biotransformation of creosote xenobiotics. Our results provide some of the first details about how these processes likely operate in the liver of a specialized mammalian herbivore.

What is new in FungiDB: a web-based bioinformatics platform for omics-scale data analysis for fungal and oomycete species

FungiDB (https://fungidb.org) serves as a valuable online resource that seamlessly integrates genomic and related large-scale data for a wide range of fungal and oomycete species. As an integral part of the VEuPathDB Bioinformatics Resource Center (https://veupathdb.org), FungiDB continually integrates both published and unpublished data addressing various aspects of fungal biology. Established in early 2011, the database has evolved to support 674 datasets. The datasets include over 300 genomes spanning various taxa (e.g. Ascomycota, Basidiomycota, Blastocladiomycota, Chytridiomycota, Mucoromycota, as well as Albuginales, Peronosporales, Pythiales, and Saprolegniales). In addition to genomic assemblies and annotation, over 300 extra datasets encompassing diverse information, such as expression and variation data, are also available. The resource also provides an intuitive web-based interface, facilitating comprehensive approaches to data mining and visualization. Users can test their hypotheses and navigate through omics-scale datasets using a built-in search strategy system. Moreover, FungiDB offers capabilities for private data analysis via the integrated VEuPathDB Galaxy platform. FungiDB also permits genome improvements by capturing expert knowledge through the User Comments system and the Apollo genome annotation editor for structural and functional gene curation. FungiDB facilitates data exploration and analysis and contributes to advancing research efforts by capturing expert knowledge for fungal and oomycete species.

WormBase 2024: status and transitioning to Alliance infrastructure

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

The Arabidopsis Information Resource in 2024

Since 1999, The Arabidopsis Information Resource (www.arabidopsis.org) has been curating data about the Arabidopsis thaliana genome. Its primary focus is integrating experimental gene function information from the peer-reviewed literature and codifying it as controlled vocabulary annotations. Our goal is to produce a "gold standard" functional annotation set that reflects the current state of knowledge about the Arabidopsis genome. At the same time, the resource serves as a nexus for community-based collaborations aimed at improving data quality, access, and reuse. For the past decade, our work has been made possible by subscriptions from our global user base. This update covers our ongoing biocuration work, some of our modernization efforts that contribute to the first major infrastructure overhaul since 2011, the introduction of JBrowse2, and the resource's role in community activities such as organizing the structural reannotation of the genome. For gene function assessment, we used gene ontology annotations as a metric to evaluate: (1) what is currently known about Arabidopsis gene function and (2) the set of "unknown" genes. Currently, 74% of the proteome has been annotated to at least one gene ontology term. Of those loci, half have experimental support for at least one of the following aspects: molecular function, biological process, or cellular component. Our work sheds light on the genes for which we have not yet identified any published experimental data and have no functional annotation. Drawing attention to these unknown genes highlights knowledge gaps and potential sources of novel discoveries.

Updates to the Alliance of Genome Resources central infrastructure

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

From tradition to innovation: conventional and deep learning frameworks in genome annotation

Following the milestone success of the Human Genome Project, the 'Encyclopedia of DNA Elements (ENCODE)' initiative was launched in 2003 to unearth information about the numerous functional elements within the genome. This endeavor coincided with the emergence of numerous novel technologies, accompanied by the provision of vast amounts of whole-genome sequences, high-throughput data such as ChIP-Seq and RNA-Seq. Extracting biologically meaningful information from this massive dataset has become a critical aspect of many recent studies, particularly in annotating and predicting the functions of unknown genes. The core idea behind genome annotation is to identify genes and various functional elements within the genome sequence and infer their biological functions. Traditional wet-lab experimental methods still rely on extensive efforts for functional verification. However, early bioinformatics algorithms and software primarily employed shallow learning techniques; thus, the ability to characterize data and features learning was limited. With the widespread adoption of RNA-Seq technology, scientists from the biological community began to harness the potential of machine learning and deep learning approaches for gene structure prediction and functional annotation. In this context, we reviewed both conventional methods and contemporary deep learning frameworks, and highlighted novel perspectives on the challenges arising during annotation underscoring the dynamic nature of this evolving scientific landscape.

Genome annotations for the ascomycete fungi Trichoderma harzianum, Trichoderma aggressivum, and Purpureocillium lilacinum

We sequenced and annotated the genomes of the ascomycete fungi Trichoderma harzianum, Trichoderma aggressivum f. europaeum, and Purpureocillium lilacinum. Moreover, we developed a website to allow users to interactively analyze the assemblies, gene predictions, and functional annotations of these species and 70+ previously sequenced fungi.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa

KEY MESSAGE

We reported the mitochondrial genome of Ventilago leiocarpa for the first time. Two and one sites lead to the generation of stop and stat codon through editing were verified. Ventilago leiocarpa, a member of the Rhamnaceae family, is frequently utilized in traditional medicine due to the medicinal properties of its roots. In this study, we successfully assembled the mitogenome of V. leiocarpa using both BGI short reads and Nanopore long reads. This mitogenome has a total length of 331,839 bp. The annotated results showed 36 unique protein-coding, 16 tRNA and 3 rRNA genes in this mitogenome. Furthermore, we confirmed the presence of a branched structure through the utilization of long reads mapping, PCR amplification, and Sanger sequencing. Specifically, the ctg1 can form a single circular molecule or combine with ctg4 to form a linear molecule. Likewise, ctg2 can form a single circular molecule or can be connected to ctg4 to form a linear molecule. Subsequently, through a comparative analysis of the mitogenome and cpgenome sequences, we identified ten mitochondrial plastid sequences (MTPTs), including two complete protein-coding genes and five complete tRNA genes. The existence of MTPTs was verified by long reads. Colinear analysis showed that the mitogenomes of Rosales were highly divergent in structure. Finally, we identified 545 RNA editing sites involving 36 protein-coding genes by Deepred-mt. To validate our findings, we conducted PCR amplification and Sanger sequencing, which confirmed the generation of stop codons in atp9-223 and rps10-391, as well as the generation of a start codon in nad4L-2. This project reported the complex structure and RNA editing event of the V. Leiocarpa mitogenome, which will provide valuable information for the study of mitochondrial gene expression.

A haplotype-like, chromosome-level assembled and annotated genome of Biomphalaria glabrata, an important intermediate host of schistosomiasis and the best studied model of schistosomiasis vector snails

Schistosomiasis is one of the world's most devastating parasitic diseases, afflicting 251 million people globally. The Neotropical snail Biomphalaria glabrata is an important intermediate host of the human blood fluke Schistosoma mansoni and a predominant model for schistosomiasis research. To fully exploit this model snail for biomedical research, here we report a haplotype-like, chromosome-level assembled and annotated genome of the homozygous iM line of B. glabrata that we developed at the University of New Mexico. Using multiple sequencing platforms, including Illumina, PacBio, and Omni-C sequencing, 18 sequence contact matrices representing 18 haploid chromosomes (2n = 36) were generated (337x genome coverage), and 96.5% of the scaffold sequences were anchored to the 18 chromosomes. Protein-coding genes (n = 34,559), non-coding RNAs (n = 2,406), and repetitive elements (42.52% of the genome) were predicted for the whole genome, and detailed annotations for individual chromosomes were also provided. Using this genomic resource, we have investigated the genomic structure and organization of the Toll-like receptor (TLR) and fibrinogen-domain containing protein (FReD) genes, the two important immune-related gene families. Notably, TLR-like genes are scattered on 13 chromosomes. In contrast, almost all (39 of 40) fibrinogen-related genes (FREPs) (immunoglobulin superfamily (IgSF) + fibrinogen (FBG)) are clustered within a 5-million nucleotide region on chromosome 13, yielding insight into mechanisms involved in the diversification of FREPs. This is the first genome of schistosomiasis vector snails that has been assembled at the chromosome level, annotated, and analyzed. It serves as a valuable resource for a deeper understanding of the biology of vector snails, especially Biomphalaria snails.

Insights into the multi-chromosomal mitochondrial genome structure of the xero-halophytic plant Haloxylon Ammodendron (C.A.Mey.) Bunge ex Fenzl

BACKGROUND

Haloxylon ammodendron holds significance as an ecological plant, showcasing remarkable adaptability to desert conditions, halophytic environments, and sand fixation. With its potential for carbon sequestration, it emerges as a promising candidate for environmental sustainability. Furthermore, it serves as a valuable C4 plant model, offering insights into the genetic foundations of extreme drought tolerance. Despite the availability of plastid and nuclear genomes, the absence of a mitochondrial genome (mitogenome or mtDNA) hinders a comprehensive understanding of its its mtDNA structure, organization, and phylogenetic implications.

RESULTS

In the present study, the mitochondrial genome of H. ammodendron was assembled and annotated, resulting in a multi-chromosomal configuration with two circular chromosomes. The mtDNA measured 210,149 bp in length and contained 31 protein-coding genes, 18 tRNA and three rRNA. Our analysis identified a total of 66 simple sequence repeats along with 27 tandem repeats, 312 forward repeats, and 303 palindromic repeats were found. Notably, 17 sequence fragments displayed homology between the mtDNA and chloroplast genome (cpDNA), spanning 5233 bp, accounting for 2.49% of the total mitogenome size. Additionally, we predicted 337 RNA editing sites, all of the C-to-U conversion type. Phylogenetic inference confidently placed H. ammodendron in the Amaranthacea family and its close relative, Suaeda glacum.

CONCLUSIONS

H. ammodendron mtDNA showed a multi-chromosomal structure with two fully circularized molecules. This newly characterized mtDNA represents a valuable resource for gaining insights into the basis of mtDNA structure variation within Caryophyllales and the evolution of land plants, contributing to their identification, and classification.

Complete mitochondrial genome of an oleaginous microalga Vischeria punctata (Eustigmatophyceae: Chlorobotryaceae) and phylogenetic analysis

Vischeria punctata, as first described by Vischer in 1945, is a member of the family Chlorobotryaceae, within the order Eustigmatales. This species is recognized for its potential as a source of biofuels and other high-value products. In the present investigation, the whole genome of V. punctata was sequenced utilizing the Illumina HiSeq 4000 platform, enabling the assembly and annotation of its complete mitochondrial genome. The resulting circular genome spans 41,528 base pairs (bp) with a guanine-cytosine (GC) content of 27.3%. This genome encompasses 36 protein-coding genes, alongside 28 transfer RNA (tRNA), and three ribosomal RNA (rRNA) genes. The evolutionary trajectory of V. punctata was further explored by constructing a phylogenetic tree derived from the mitochondrial 33 gene dataset of 16 Ochrophyta species. Comparative analysis reveals that V. punctata bears closer ties to Vischeria sp. CAUP Q202 than to Vischeria stellata strain SAG 33.83, suggesting shared evolutionary pathways and phenotypic traits. This investigation constitutes the inaugural study into the mitochondrial evolution and phylogenetic patterning of the mitogenome in V. punctata. The outcomes from this research bolster our understanding of the genetic diversity and evolutionary processes within the class Eustigmatophyceae. In particular, the mitochondrial genome of V. punctata serves as a valuable resource in elucidating these aspects.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

Diaci v3.0: chromosome-level assembly, de novo transcriptome, and manual annotation of Diaphorina citri, insect vector of Huanglongbing

BACKGROUND

Diaphorina citri is an insect vector of "Candidatus Liberibacter asiaticus" (CLas), the gram-negative bacterial pathogen associated with citrus greening disease. Control measures rely on pesticides with negative impacts on the environment, natural ecosystems, and human and animal health. In contrast, gene-targeting methods have the potential to specifically target the vector species and/or reduce pathogen transmission.

RESULTS

To improve the genomic resources needed for targeted pest control, we assembled a D. citri genome based on PacBio long reads followed by proximity ligation-based scaffolding. The 474-Mb genome has 13 chromosomal-length scaffolds. In total, 1,036 genes were manually curated as part of a community annotation project, composed primarily of undergraduate students. We also computationally identified a total of 1,015 putative transcription factors (TFs) and were able to infer motifs for 337 TFs (33%). In addition, we produced a genome-independent transcriptome and genomes for D. citri endosymbionts.

CONCLUSIONS

Manual annotation provided more accurate gene models for use by researchers and provided an excellent training opportunity for students from multiple institutions. All resources are available on CitrusGreening.org and NCBI. The chromosomal-length D. citri genome assembly serves as a blueprint for the development of collaborative genomics projects for other medically and agriculturally significant insect vectors.

Comparative complete chloroplast genome of Geum japonicum: evolution and phylogenetic analysis

Geum japonicum (Rosaceae) has been widely used in China as a traditional herbal medicine due to its high economic and medicinal value. However, the appearance of Geum species is relatively similar, making identification difficult by conventional phenotypic methods, and the studies of genomics and species evolution are lacking. To better distinguish the medicinal varieties and fill this gap, we carried out relevant research on the chloroplast genome of G. japonicum. Results show a typical quadripartite structure of the chloroplast genome of G. japonicum with a length of 156,042 bp. There are totally 131 unique genes in the genome, including 87 protein-coding genes, 36 tRNA genes, and 8 rRNA genes, and there were also 87 SSRs identified and mostly mononucleotide Adenine-Thymine. We next compared the plastid genomes among four Geum species and obtained 14 hypervariable regions, including ndhF, psbE, trnG-UCC, ccsA, trnQ-UUG, rps16, psbK, trnL-UAA, ycf1, ndhD, atpA, petN, rps14, and trnK-UUU. Phylogenetic analysis revealed that G. japonicum is most closely related to Geum aleppicum, and possibly has some evolutionary relatedness with an ancient relic plant Taihangia rupestris. This research enriched the genome resources and provided fundamental insights for evolutionary studies and the phylogeny of Geum.

VPS13B is localized at the cis-trans Golgi complex interface and is a functional partner of FAM177A1

Mutations in VPS13B, a member of a protein family implicated in bulk lipid transport between adjacent membranes, cause Cohen syndrome. VPS13B is known to be concentrated in the Golgi complex, but its precise location within this organelle and thus the site(s) where it achieves lipid transport remains unclear. Here we show that VPS13B is localized at the interface between cis and trans Golgi sub-compartments and that Golgi complex re-formation after Brefeldin A (BFA) induced disruption is delayed in VPS13B KO cells. This delay is phenocopied by loss of FAM177A1, a Golgi complex protein of unknown function reported to be a VPS13B interactor and whose mutations also result in a developmental disorder. In zebrafish, the vps13b orthologue, not previously annotated in this organism, genetically interacts with fam177a1. Collectively, these findings raise the possibility that bulk lipid transport by VPS13B may play a role in expanding Golgi membranes and that VPS13B may be assisted in this function by FAM177A1.

The alliance of genome resources: transforming comparative genomics

Comparing genomic and biological characteristics across multiple species is essential to using model systems to investigate the molecular and cellular mechanisms underlying human biology and disease and to translate mechanistic insights from studies in model organisms for clinical applications. Building a scalable knowledge commons platform that supports cross-species comparison of rich, expertly curated knowledge regarding gene function, phenotype, and disease associations available for model organisms and humans is the primary mission of the Alliance of Genome Resources (the Alliance). The Alliance is a consortium of seven model organism knowledgebases (mouse, rat, yeast, nematode, zebrafish, frog, fruit fly) and the Gene Ontology resource. The Alliance uses a common set of gene ortholog assertions as the basis for comparing biological annotations across the organisms represented in the Alliance. The major types of knowledge associated with genes that are represented in the Alliance database currently include gene function, phenotypic alleles and variants, human disease associations, pathways, gene expression, and both protein-protein and genetic interactions. The Alliance has enhanced the ability of researchers to easily compare biological annotations for common data types across model organisms and human through the implementation of shared programmatic access mechanisms, data-specific web pages with a unified "look and feel", and interactive user interfaces specifically designed to support comparative biology. The modular infrastructure developed by the Alliance allows the resource to serve as an extensible "knowledge commons" capable of expanding to accommodate additional model organisms.

Updates to the Alliance of Genome Resources Central Infrastructure Alliance of Genome Resources Consortium

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively-studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, C. elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and APIs. Here we focus on developments over the last two years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse. We describe our progress towards a central persistent database to support curation, the data modeling that underpins harmonization, and progress towards a state-of-the art literature curation system with integrated Artificial Intelligence and Machine Learning (AI/ML).

Genome-wide characterization of regulator of chromosome condensation 1 (RCC1) gene family in Artemisia annua L. revealed a conservation evolutionary pattern

BACKGROUND

Artemisia annua is the major source for artemisinin production. The artemisinin content in A. annua is affected by different types of light especially the UV light. UVR8, a member of RCC1 gene family was found to be the UV-B receptor in plants. The gene structures, evolutionary history and expression profile of UVR8 or RCC1 genes remain undiscovered in A. annua.

RESULTS

Twenty-two RCC1 genes (AaRCC1) were identified in each haplotype genome of two diploid strains of A. annua, LQ-9 and HAN1. Varied gene structures and sequences among paralogs were observed. The divergence of most RCC1 genes occurred at 46.7 - 51 MYA which overlapped with species divergence of core Asteraceae during the Eocene, while no recent novel RCC1 members were found in A. annua genome. The number of RCC1 genes remained stable among eudicots and RCC1 genes underwent purifying selection. The expression profile of AaRCC1 is analogous to that of Arabidopsis thaliana (AtRCC1) when responding to environmental stress.

CONCLUSIONS

This study provided a comprehensive characterization of the AaRCC1 gene family and suggested that RCC1 genes were conserved in gene number, structures, constitution of amino acids and expression profiles among eudicots.