A novel hybrid gene prediction method employing protein multiple sequence alignments

Chromosome-level genome assembly, annotation, and population genomic resource of argali (Ovis ammon)

Argali stands as the largest species among wild sheep in Central and East Asia, with a concerning rate of decline estimated at 30%. The intraspecific taxonomy of argali remains contentious due to limited genomic data and unclear geographic separation. In this study, we constructed a chromosome-level genome assembly and annotation for the Tibetan argali (O. a. hodgsoni), together with population genomic resequencing of 32 individuals representing four subspecies. The contig-level genome was 2.64 Gb in size, with a contig N50 length of 71.69 Mb and an estimated genomic completeness of 96.01%. Using Hi-C sequencing data scaffolding, 99.90% of initially assembled sequences were mapped and oriented onto 28 pseudo-chromosomes except the Y chromosome. Annotation uncovered 21,564 protein-coding genes and 46.38% repeat sequences. The average coverage of the population resequencing data was 23.74 with mean mapping ratio up to of 97.19%. The high-quality genome assembly and annotation of the Tibetan argali, coupled with the high-depth population genomic data, will serve as a valuable genetic resource for studies on the taxonomy and conservation of argali.

Chromosome-level genome assemblies of sunflower oilseed and confectionery cultivars

The sunflower (Helianthus annuus L.), belonging to the Asteraceae family, is the world's fourth most important oil crop. Sunflower cultivars are categorized into oilseed and confectionery types. Here, we present chromosome-level genome assemblies of two Chinese sunflower cultivars-oilseed and confectionery-using PacBio HiFi and Hi-C sequencing. The oilseed cultivar, OXS, has a genome assembly spanning 3.03 Gb with 99.58% of sequences anchored to 17 chromosomes and a contig N50 length of 154.78 Mb. The first published confectionery cultivar genome, YDS, mirrors this closely with a 3.02 Gb assembly, contig N50 length of 153.87 Mb and 99.40% of sequences mapped similarly. Gene completeness reached 98.2% for OXS and 98.4% for YDS, with LTR Assembly Index scores of 24.73 and 25.85, respectively. Comparative genomics identified rapidly evolving gene families linked to synthesis, growth, and stress defense. Additionally, we found high collinearity between the YDS and OXS genomes, despite three significant inversions, and detected 15,056 large deletions and insertions. These findings lay a robust foundation for advanced genomic research and breeding innovations in sunflowers.

A novel mutation in SoIAA20 confers cross-resistance to 2,4-Dichlorophenoxyacetic acid and other auxinic herbicides in Sonchus oleraceus

BACKGROUND

2,4-Dichlorophenoxyacetic acid (2,4-D) and other auxinic herbicides are important for weed control in cropping systems globally. Weeds with resistance to 2,4-D and other auxinic herbicides have evolved, including several populations of Sonchus oleraceus from multiple sites in Australia. We report the underlying mechanism in these populations that gives rise to auxinic herbicide resistance.

RESULTS

We studied a total of three susceptible and eight resistant Sonchus oleraceus populations. All resistant populations had a deletion of three amino acids flanking the degron sequence of an Aux/IAA gene, SoIAA20, which was not found in the three susceptible populations. The eight populations with the resistant allele were also resistant to dicamba, fluroxypyr and clopyralid. The resistant plants also had reduced movement of 2,4-D out of the treated tissues compared to susceptible plants.

CONCLUSION

The paired deletion flanking the degron region of SoIAA20 likely provides resistance to 2,4-D by restricting the movement of 2,4-D from the treated tissue to the rest of the plant. We hypothesise that this deletion keeps the 2,4-D bound to the target site. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Comparative Genomics Reveals Sources of Genetic Variability in the Asexual Fungal Plant Pathogen Colletotrichum lupini

Fungal plant pathogens cause major crop losses worldwide, with many featuring compartmentalised genomes that include both core and accessory regions, which are believed to drive adaptation. The highly host-specific fungus Colletotrichum lupini greatly impacts lupin (Lupinus spp.) cultivation. This pathogen is part of clade 1 of the C. acutatum species complex and comprises four genetically uniform, presumably clonal, lineages (I-IV). Despite this, variation in virulence and morphology has been observed within these lineages. To investigate the potential sources of genetic variability in this asexual fungus, we compared the genomes of 16 C. lupini strains and 17 related Colletotrichum species. Phylogenomics confirmed the presence of four distinct lineages, but further examination based on genome size, gene content, transposable elements (TEs), and deletions revealed that lineage II could be split into two groups, II-A and II-B. TE content varied between lineages and correlated strongly with genome size variation, supporting a role for TEs in genome expansion in this species. Pangenome analysis revealed a highly variable accessory genome, including a minichromosome present in lineages II, III, and IV, but absent in lineage I. Accessory genes and effectors appeared to cluster in proximity to TEs. Presence/absence variation of putative effectors was lineage-specific, suggesting that these genes play a crucial role in determining host range. Notably, no effectors were found on the TE-rich minichromosome. Our findings shed light on the potential mechanisms generating genetic diversity in this asexual fungal pathogen that could aid future disease management.

A draft genome assembly of the reef-building coral Acropora hemprichii from the central Red Sea

Coral reef ecosystems are under threat from climate change. Thus, active interventions to spur coral conservation/restoration are critical to support reef survival, greatly informed by a molecular understanding of resilience. The genus Acropora is a species-rich and globally prevalent reef builder that has experienced dramatic declines in the Caribbean. Here we generated a draft genome of the common coral Acropora hemprichii from the central Red Sea, one of the warmest water bodies in the world. We assembled the genome using 10x Chromium sequencing with subsequent scaffolding using a reference genome and Illumina short-read sequencing contigs. The A. hemprichii genome has an assembly size of 495.6 Mb confirmed using physical size estimation, of which 247.8 Mb (50%) are repeats. The scaffold N50 is 1.38 Mb with 99.6% of BUSCO genes identified (93.7% complete, 5.9% fragmented), providing a set of 26,865 protein-coding genes. The Red Sea A. hemprichii reference genome provides a valuable resource for studies aiming to decode the genomic architecture of resilience, e.g. through comparative analyses with other Acropora genomes.

Roadmap and Considerations for Genome Editing in a Non-Model Organism: Genetic Variations and Off-Target Profiling

The CRISPR/Cas genome editing approach in non-model organisms poses challenges that remain to be resolved. Here, we demonstrated a generalized roadmap for a de novo genome annotation approach applied to the non-model organism Macrobrachium rosenbergii. We also addressed the typical genome editing challenges arising from genetic variations, such as a high frequency of single nucleotide polymorphisms, differences in sex chromosomes, and repetitive sequences that can lead to off-target events. For the genome editing of M. rosenbergii, our laboratory recently adapted the CRISPR/Cas genome editing approach to embryos and the embryonic primary cell culture. In this continuation study, an annotation pipeline was trained to predict the gene models by leveraging the available genomic, transcriptomic, and proteomic data, and enabling accurate gene prediction and guide design for knock-outs. A next-generation sequencing analysis demonstrated a high frequency of genetic variations in genes on both autosomal and sex chromosomes, which have been shown to affect the accuracy of editing analyses. To enable future applications based on the CRISPR/Cas tool in non-model organisms, we also verified the reliability of editing efficiency and tracked off-target frequencies. Despite the lack of comprehensive information on non-model organisms, this study provides an example of the feasibility of selecting and editing specific genes with a high degree of certainty.

Discovery of Alanomyces manoharacharyi: A Novel Fungus Identified Using Genome Sequencing and Metabolomic Analysis

In this study, a new species of Alanomyces was isolated as an endophyte from the bark of Azadirachta indica from Mulshi, Maharashtra. The identity of this isolate was confirmed based on the asexual morphological characteristics as well as multi-gene phylogeny based on the internal transcribed spacer (ITS) and large subunit (LSU) nuclear ribosomal RNA (rRNA) regions. As this was the second species to be reported in this genus, we sequenced the genome of this species to increase our knowledge about the possible applicability of this genus to various industries. Its genome length was found to be 35.01 Mb, harboring 7870 protein-coding genes as per Augustus and 8101 genes using GeMoMa. Many genes were annotated using the Clusters of Orthologous Groups (COGs) database, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), Swiss-Prot, NCBI non-redundant nucleotide sequences (NTs), and NCBI non-redundant protein sequences (NRs). The number of repeating sequences was predicted using Proteinmask and RepeatMasker; tRNA were detected using tRNAscan and snRNA were predicted using rfam_scan. The genome was also annotated using the Pathogen-Host Interactions Database (PHI-base) and AntiSMASH. To confirm the evolutionary history, average nucleotide identity (ANIb), phylogeny based on orthologous proteins, and single nucleotide polymorphisms (SNPs) were carried out. Metabolic profiling of the methanolic extract of dried biomass and ethyl acetate extract of the filtrate revealed a variety of compounds of great importance in the pharmaceutical and cosmetic industry. The characterization and genomic analysis of the newly discovered species Alanomyces manoharacharyi highlights its potential applicability across multiple industries, particularly in pharmaceuticals and cosmetics due to its diverse secondary metabolites and unique genetic features it possesses.

A pilot study to directly estimate radiation-induced mutation in large Japanese field mouse duo sample, mother and offspring, excluding unknown father, using ddRAD sequencing

DNA mutations are one of the effects of radiation exposure. A large amount of radioactive materials was released into the environment from the Fukushima Daiichi Nuclear Power Plant after a major earthquake and tsunami. Wild animals and plants living in highly radiation-contaminated areas are constantly exposed to high doses of radiation, and concerns occur about its effects on their health and the next generations. As a pilot study, double-digest restriction site-associated DNA (ddRAD) sequencing was conducted to assess the incidence of mutations in wild large Japanese field mice collected from the evacuation area. The optimal combination of restriction enzymes, encompassing the functionally important coding regions, was selected using in silico analysis. These enzymes were used for ddRAD sequence analysis of females and their fetuses to evaluate mutation rates. The results indicated that no significant differences were observed in mutation rates between mothers and fetuses in the study areas.

Fine-scale genomic analysis of the tree endophyte Annulohypoxylon sp. FPYF3050 producing monoterpene 1,8-cineole

1,8-Cineole has a potential in control crop pests and biofuels. The endophytic fungus, Annulohypoxylon sp. FPYF3050 (Neolitsea pulchella), can produce over 90% 1,8-cineole of relative area in its natural volatiles, possessed nematicide properties. The annotated genome of this strain will provide insights into potential application in biofumigation and terpene-based advanced biofuel.

Molecular evolution and genetic diversity of defective chorion 1 in Anastrepha fraterculus and Anastrepha obliqua (Diptera, Tephritidae)

The family Tephritidae comprises numerous fruit fly species, some of which are economically significant, such as several in the genus Anastrepha. Most pest species in this genus belong to the fraterculus group, characterized by closely related species that are difficult to differentiate due to recent divergence and gene flow. Identifying genetic markers for their study is paramount for understanding the group's evolution and eventual phytosanitary control. Because there is variation in eggshell morphology among species in the genus, the study of the rapidly evolving defective chorion 1 (dec-1) gene, which is crucial for chorion formation and reproduction, could provide relevant information for Anastrepha differentiation. We compared transcriptome sequences of dec-1 from two of the most important pest species in the genus, Anastrepha fraterculus and Anastrepha obliqua to dec-1 sequences from Anastrepha ludens, which was used for structure prediction. Furthermore, we amplified a conserved exon across populations of these species. These data revealed three alternative transcripts in A. fraterculus and A. obliqua, consistent with patterns found in other Tephritidae; we obtained orthologous sequences for these other tephritids from NCBI to investigate patterns of selection affecting this gene at different hierarchical levels using different methods. These analyses show a general pattern of purifying selection across the whole gene and throughout its history at different hierarchical levels, from populations to more distantly related species. That notwithstanding, we still found evidence of positive and episodic diversifying selection at different levels. Different parts of the gene have shown distinct evolutionary rates, which were associated with the diverse proproteins produced by posttranslational changes of DEC-1, with proproteins that are incorporated in the chorion earlier in egg formation being in general more conserved than others that are incorporated later. This correlation appears more evident in certain lineages, including the branch that separates Anastrepha, as well as other internal branches that differentiate species within the genus. Our data showed that this gene shows remarkable variation across its different exons, which has proven to be informative at different evolutionary levels. These changes hold promise not only for studying differentiation in Anastrepha but also for the eventual management of selected pest species.

Near complete assembly of Pyricularia penniseti infecting Cenchrus grass identified its eight core chromosomes

Fungi from the Pyricularia genus cause blast disease in many economically important crops and grasses, such as wheat, rice, and Cenchrus grass JUJUNCAO. Structure variation associated with the gain and loss of effectors contributes largely to the adaptive evolution of this fungus towards diverse host plants. A telomere-to-telomere genome assembly would facilitate the identification of genome-wide structural variations through comparative genomics. Here, we report a telomere-to-telomere, near-complete genome assembly of a Pyricularia penniseti isolate JC-1 infecting JUJUNCAO. The assembly consists of eight core chromosomes and two supernumerary chromosomes, named mini1 and mini2, spanning 42.1 Mb. We annotated 12,156 protein-coding genes and identified 4.54% of the genome as repetitive sequences. The two supernumerary chromosomes contained fewer genes and more repetitive sequences than the core chromosomes. Our genome and results provide valuable resources for the future study in genome evolution, structure variation and host adaptation of the Pyricularia fungus.

Bank vole genomics links determinate and indeterminate growth of teeth

BACKGROUND

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Bank and prairie voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars while retaining similar size and shape, providing alternative models for studying roots.

RESULTS

We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. Bulk transcriptomics comparisons of embryonic molar development between bank voles and mice demonstrated overall conservation of gene expression levels, with species-specific differences corresponding to the accelerated and more extensive patterning of the vole molar. We leverage convergent evolution of unrooted molars across the clade to examine changes that may underlie the evolution of unrooted molars. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation.

CONCLUSIONS

Our results support ongoing evolution of dental genes across Glires and identify candidate genes for mechanistic studies of root formation. Comparative research using the bank vole as a model species can reveal the complex evolutionary background of convergent evolution for ever-growing molars.

Chromosome-level genome assembly of American sweetgum (Liquidambar styraciflua, Altingiaceae)

The deciduous American sweetgum (Liquidambar styraciflua, Altingiaceae) is a popular ornamental and economically valuable tree renowned for its sweet-smelling bark resin, abundant volatile substances, and spectacular fall leaf color. However, the absence of a reference genome hinders thorough investigations into the mechanisms underlying phenotypic variation, secondary metabolite synthesis and adaptation, both in this species and other Liquidambar members. In this study, we sequenced and constructed a chromosome-level assembly of the L. styraciflua genome, covering 662.48 Mb with a scaffold N50 of 39.54 Mb, by integrating PacBio, Illumina and chromosome conformation capture data. We identified 58.83% of the genome sequences as repetitive elements and 25,713 protein-coding genes, 97.28% of which were functionally annotated. The genome sequencing reads, assembly and annotation data have been deposited in publicly available repositories. This high-quality genome assembly provides valuable resources for further evolutionary and functional genomic studies in American sweetgum and other Liquidambar species.

Chromosome-level genome assembly and comparative genomics shed light on Helicoverpa assulta ecology and pest management

BACKGROUND

The Oriental tobacco budworm, Helicoverpa assulta, a specialist herbivorous insect that exclusively feeds on plants of the Solanaceae family, causes considerable damage to crops, such as tobacco and hot pepper. The absence of a genome sequence for this species hinders further research on its pest management and ecological adaptation.

RESULTS

Here, we present a high-quality chromosome-level genome of a Korean strain of H. assulta (Pyeongchang strain, K18). The total assembly spans 424.4 Mb with an N50 length of 14.54 Mb and 37% GC content. The assembled genome (ASM2961881v1) comprises 31 chromosomes, similar to other congeneric generalist species including H. armigera and H. zea. In terms of genomic assembly quality, the complete BUSCOs and repeat content accounted for 98.3% and 33.01% of the genome, respectively. Based on this assembly, 19 485 protein-coding genes were predicted in the genome annotation. A comparative analysis was conducted using the identified number of protein-coding genes in H. armigera (24154) and H. zea (23696). Out of the 19 485 predicted genes, 137 genes in 15 orthogroups were found to have expanded significantly in H. assulta, while 149 genes in 95 orthogroups contracted rapidly.

CONCLUSION

This study revealed specific gene expansions and contractions in H. assulta compared to those in its close relatives, indicating potential adaptations related to its specialized feeding habits. Also, the comparative genome analysis provides valuable insights for the integrated pest management of H. assulta and other globally significant pests in the Heliothinae subfamily. © 2024 Society of Chemical Industry.

Structure and Biosynthesis of Perochalasins A-C, Open-Chain Merocytochalasans Produced by the Marine-Derived Fungus Peroneutypa sp. M16

Novel open-chain merocytochalasans, perochalasins A-C (1-3), containing an unusual N-O six-membered heterocyclic moiety, were isolated from cultures of the marine-derived Peroneutypa sp. M16 fungus, along with cytochalasin Z27 (4), cytochalasin Z28 (5), [12]-cytochalasin (6), and phenochalasin B (7). The structures of compounds 1-3 were established by analysis of the spectroscopic data. Full genome sequencing of Peroneutypa sp. M16 enabled the identification of a cytochalasan biosynthetic gene cluster and a proposal for the biosynthetic assembly of perochalasins. The proposal is supported by the nonenzymatic conversion of phenochalasin B (7) into 1-3, based on isotope-labeled hydroxylamine (15NH2OH and ND2OD) feeding studies in vivo and in vitro. In contrast to other merocytochalasans, these are the first cytochalasans confirmed to arise via nucleophilic addition and at a distinct location from the reactive macrocycle olefin, potentially expanding further the range of merocytochalasans to be discovered or engineered. Cytochalasin Z27 (4) exhibited antiplasmodial activities in the low micromolar range against the chloroquine-sensitive Plasmodium falciparum 3D7 strain as well as against resistant strains of the parasite (Dd2, TM90C6B, and 3D7r_MMV848).

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

Latent infection of an active giant endogenous virus in a unicellular green alga

Latency is a common strategy in a wide range of viral lineages, but its prevalence in giant viruses remains unknown. Here we describe the activity and viral production from a 617 kbp integrated giant viral element in the model green alga Chlamydomonas reinhardtii. We resolve the integrated viral region using long-read sequencing and show that viral particles are produced and released in otherwise healthy cultures. A diverse array of viral-encoded selfish genetic elements are expressed during GEVE reactivation and produce proteins that are packaged in virions. In addition, we show that field isolates of Chlamydomonas sp. harbor latent giant viruses related to the C. reinhardtii GEVE that exhibit similar infection dynamics, demonstrating that giant virus latency is prevalent in natural host communities. Our work reports the largest temperate virus documented to date and the first active GEVE identified in a unicellular eukaryote, substantially expanding the known limits of viral latency.

Genome Sequencing of Lentinula edodes Revealed a Genomic Variant Block Associated with a Thermo-Tolerant Trait in Fruit Body Formation

The formation of multicellular fruiting bodies in basidiomycete mushrooms is a crucial developmental process for sexual reproduction and subsequent spore development. Temperature is one of the most critical factors influencing the phase transition for mushroom reproduction. During the domestication of mushrooms, traits related to fruiting bodies have significantly impacted agricultural adaptation and human preferences. Recent research has demonstrated that chromosomal variations, such as structural variants (SVs) and variant blocks (VBs), play crucial roles in agronomic traits and evolutionary processes. However, the lack of high-quality genomic information and important trait data have hindered comprehensive identification and characterization in Lentinula edodes breeding processes. In this study, the genomes of two monokaryotic L. edodes strains, characterized by thermo-tolerance and thermo-sensitivity during fruiting body formation, were reassembled at the chromosomal level. Comparative genomic studies of four thermo-tolerant and thermo-sensitive monokaryotic L. edodes strains identified a 0.56 Mbp variant block on chromosome 9. Genes associated with DNA repair or cellular response to DNA damage stimulus were enriched in this variant block. Finally, we developed eight CAPS markers from the variant block to discriminate the thermo-tolerant traits in L. edodes cultivars. Our findings show that the identified variant block is highly correlated with the thermo-tolerant trait for fruiting body formation and that alleles present in this block may have been artificially selected during L. edodes domestication.

The Chrysosplenium sinicum genome provides insights into adaptive evolution of shade plants

Chrysosplenium sinicum, a traditional Tibetan medicinal plant, can successfully thrive in low-light environments for long periods of time. To investigate the adaptive evolution of shade plants in low-light environments, we generated a chromosome-scale genome assembly (~320 Mb) for C. sinicum by combining PacBio sequencing and Hi-C technologies. Based on our results, gene families related to photosynthesis and cell respiration greatly expanded and evolved in C. sinicum genome due to intracellular DNA transfer from organelle genome to nuclear genome. Under positive selective pressure, adaptive evolution of light-harvesting complex II (LHCII) component protein CsLhcb1s resulted in the expansion of threonine residues at the phosphorylation site of STN7 kinase, potentially establishing a crucial genomic foundation for enhancing C. sinicum's adaptability in low-light environments. Through transcriptome and metabolome analysis, we identified chrysosplenol and chrysosplenoside as predominant flavonoid metabolites of C. sinicum and predicted their synthesis pathways. In addition, analysis of alternative splicing (AS) revealed that AS events help regulate state transition and flavonoid biosynthesis. The present study provides new insights into the genomes of shade plants exposed to low-light conditions and adaptive evolution of these genomes; in addition, the results improve our current knowledge on the biosynthetic and regulatory processes of chrysosplenol and chrysosplenoside.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

Evolution of the biosynthetic pathways of terpene scent compounds in roses

It is unknown why roses are terpene-rich, what the terpene biosynthetic pathways in roses are, and why only a few rose species produce the major components of rose essential oil. Here, we assembled two high-quality chromosome-level genomes for Rosa rugosa and Rosa multiflora. We also re-sequenced 132 individuals from the F1 progeny of Rosa chinensis and Rosa wichuraiana and 36 of their related species. Comparative genomics revealed that expansions of the 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) and terpene synthases (TPSs) gene families led to the enrichment of terpenes in rose scent components. We constructed a terpene biosynthesis network and discovered a TPS-independent citronellol biosynthetic pathway in roses through gene functional identification, genome-wide association studies (GWASs), and multi-omic analysis. Heterologous co-expression of rose citronellol biosynthetic genes in Nicotiana benthamiana led to citronellol production. Our genomic and metabolomic analyses suggested that the copy number of NUDX1-1a determines the citronellol content in different rose species. Our findings not only provide additional genome and gene resources and reveal the evolution of the terpene biosynthetic pathways but also present a nearly complete scenario for terpenoid metabolism that will facilitate the breeding of fragrant roses and the production of rose oil.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Hi-C metagenomics facilitate comparative genome analysis of bacteria and yeast from spontaneous beer and cider

Sequence-based analysis of fermented foods and beverages' microbiomes offers insights into their impact on taste and consumer health. High-throughput metagenomics provide detailed taxonomic and functional community profiling, but bacterial and yeast genome reconstruction and mobile genetic elements tracking are to be improved. We established a pipeline for exploring fermented foods microbiomes using metagenomics coupled with chromosome conformation capture (Hi-C metagenomics). The approach was applied to analyze a collection of spontaneously fermented beers and ciders (n = 12). The Hi-C reads were used to reconstruct the metagenome-assembled genomes (MAGs) of bacteria and yeasts facilitating subsequent comparative genomic analysis, assembly scaffolding and exploration of "plasmid-bacteria" links. For a subset of beverages, yeasts were isolated and characterized phenotypically. The reconstructed Hi-C MAGs primarily belonged to the Lactobacillaceae family in beers, along with Acetobacteraceae and Enterobacteriaceae in ciders, exhibiting improved quality compared to conventional metagenomic MAGs. Comparative genomic analysis of Lactobacillaceae Hi-C MAGs revealed clustering by niche and suggested genetic determinants of survival and probiotic potential. For Pediococcus damnosus, Hi-C-based networks of contigs enabled linking bacteria with plasmids. Analyzing phylogeny and accessory genes in the context of known reference genomes offered insights into the niche specialization of beer lactobacilli. The subspecies-level diversity of cider Tatumella spp. was disentangled using a Hi-C-based graph. We obtained highly complete yeast Hi-C MAGs primarily represented by Brettanomyces and Saccharomyces, with Hi-C-facilitated chromosome-level genome assembly for the former. Utilizing Hi-C metagenomics to unravel the genomic content of individual species can provide a deeper understanding of the ecological interactions within the food microbiome, aid in bioprospecting beneficial microorganisms, improving quality control and improving innovative fermented products.

Bioinformatics-based identification of GH12 endoxyloglucanases in citrus-pathogenic Penicillium spp

Millions of tons of citrus peel waste are produced every year as a byproduct of the juice industry. Citrus peel is rich in pectin and xyloglucan, but while the pectin is extracted for use in the food industry, the xyloglucan is currently not valorized. To target hydrolytic degradation of citrus peel xyloglucan into oligosaccharides, we have used bioinformatics to identify three glycoside hydrolase 12 (GH12) endoxyloglucanases (EC 3.2.1.151) from the citrus fruit pathogens Penicillium italicum GL-Gan1 and Penicillium digitatum Pd1 and characterized them on xyloglucan obtained by alkaline extraction from citrus peel. The enzymes displayed pH-temperature optima of pH 4.6-5.3 and 35-37°C. PdGH12 from P. digitatum and PiGH12A from P. italicum share 84% sequence identity and displayed similar kinetics, although kcat was highest for PdGH12. In contrast, PiGH12B from P. italicum, which has the otherwise conserved Trp in subsite -4 replaced with a Tyr, displayed a 3 times higher KM and a 4 times lower kcat/KM than PiGH12A, but was the most thermostable enzyme of the three Penicillium-derived endoxyloglucanases. The benchmark enzyme AnGH12 from Aspergillus nidulans was more thermally stable and had a higher pH-temperature optimum than the enzymes from Penicillum spp. The difference in structure of the xyloglucan oligosaccharides extracted from citrus peel xyloglucan and tamarind xyloglucan by the new endoxyloglucanases was determined by LC-MS. The inclusion of citrus peel xyloglucan demonstrated that the endoxyloglucanases liberated fucosylated xyloglucan oligomers, implying that these enzymes have the potential to upgrade citrus peel residues to produce oligomers useful as intermediates or bioactive compounds.

Identification of homologs of the Chlamydia trachomatis effector CteG reveals a family of Chlamydiaceae type III secreted proteins that can be delivered into host cells

Chlamydiae are a large group of obligate endosymbionts of eukaryotes that includes the Chlamydiaceae family, comprising several animal pathogens. Among Chlamydiaceae, Chlamydia trachomatis causes widespread ocular and urogenital infections in humans. Like many bacterial pathogens, all Chlamydiae manipulate host cells by injecting them with type III secretion effector proteins. We previously characterized the C. trachomatis effector CteG, which localizes at the host cell Golgi and plasma membrane during distinct phases of the chlamydial infectious cycle. Here, we show that CteG is a Chlamydiaceae-specific effector with over 60 homologs phylogenetically categorized into two distinct clades (CteG I and CteG II) and exhibiting several inparalogs and outparalogs. Notably, cteG I homologs are syntenic to C. trachomatis cteG, whereas cteG II homologs are syntenic among themselves but not with C. trachomatis cteG. This indicates a complex evolution of cteG homologs, which is unique among C. trachomatis effectors, marked by numerous events of gene duplication and loss. Despite relatively modest sequence conservation, nearly all tested CteG I and CteG II proteins were identified as type III secretion substrates using Yersinia as a heterologous bacterial host. Moreover, most of the type III secreted CteG I and CteG II homologs were delivered by C. trachomatis into host cells, where they localized at the Golgi region and cell periphery. Overall, this provided insights into the evolution of bacterial effectors and revealed a Chlamydiaceae family of type III secreted proteins that underwent substantial divergence during evolution while conserving the capacity to localize at specific host cell compartments.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

The de novo, chromosome-level genome assembly of the sweet chestnut (Castanea sativa Mill.) Cv. Marrone Di Chiusa Pesio

OBJECTIVES

The sweet chestnut Castanea sativa Mill. is the only native Castanea species in Europe, and it is a tree of high economic value that provides appreciated fruits and valuable wood. In this study, we assembled a high-quality nuclear genome of the ancient Italian chestnut variety 'Marrone di Chiusa Pesio' using a combination of Oxford Nanopore Technologies long reads, whole-genome and Omni-C Illumina short reads.

DATA DESCRIPTION

The genome was assembled into 238 scaffolds with an N50 size of 21.8 Mb and an N80 size of 7.1 Mb for a total assembled sequence of 750 Mb. The BUSCO assessment revealed that 98.6% of the genome matched the embryophyte dataset, highlighting good completeness of the genetic space. After chromosome-level scaffolding, 12 chromosomes with a total length of 715.8 and 713.0 Mb were constructed for haplotype 1 and haplotype 2, respectively. The repetitive elements represented 37.3% and 37.4% of the total assembled genome in haplotype 1 and haplotype 2, respectively. A total of 57,653 and 58,146 genes were predicted in the two haplotypes, and approximately 73% of the genes were functionally annotated using the EggNOG-mapper. The assembled genome will be a valuable resource and reference for future chestnut breeding and genetic improvement.

A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

The ancestral karyotype of the Heliantheae Alliance, herbicide resistance, and human allergens: Insights from the genomes of common and giant ragweed

Ambrosia artemisiifolia and Ambrosia trifida (Asteraceae) are important pest species and the two greatest sources of aeroallergens globally. Here, we took advantage of a hybrid to simplify genome assembly and present chromosome-level assemblies for both species. These assemblies show high levels of completeness with Benchmarking Universal Single-Copy Ortholog (BUSCO) scores of 94.5% for A. artemisiifolia and 96.1% for A. trifida and long terminal repeat (LTR) Assembly Index values of 26.6 and 23.6, respectively. The genomes were annotated using RNA data identifying 41,642 genes in A. artemisiifolia and 50,203 in A. trifida. More than half of the genome is composed of repetitive elements, with 62% in A. artemisiifolia and 69% in A. trifida. Single copies of herbicide resistance-associated genes PPX2L, HPPD, and ALS were found, while two copies of the EPSPS gene were identified; this latter observation may reveal a possible mechanism of resistance to the herbicide glyphosate. Ten of the 12 main allergenicity genes were also localized, some forming clusters with several copies, especially in A. artemisiifolia. The evolution of genome structure has differed among these two species. The genome of A. trifida has undergone greater rearrangement, possibly the result of chromoplexy. In contrast, the genome of A. artemisiifolia retains a structure that makes the allotetraploidization of the most recent common ancestor of the Heliantheae Alliance the clearest feature of its genome. When compared to other Heliantheae Alliance species, this allowed us to reconstruct the common ancestor's karyotype-a key step for furthering of our understanding of the evolution and diversification of this economically and allergenically important group.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Pullenvalenes A-D: Nitric Oxide-Mediated Transcriptional Activation (NOMETA) Enables Discovery of Triterpene Aminoglycosides from Australian Termite Nest-Derived Fungi

We report on the use of nitric oxide-mediated transcriptional activation (NOMETA) as an innovative means to detect and access new classes of microbial natural products encoded within silent biosynthetic gene clusters. A small library of termite nest- and mangrove-derived fungi and actinomyces was subjected to cultivation profiling using a miniaturized 24-well format approach (MATRIX) in the presence and absence of nitric oxide, with the resulting metabolomes subjected to comparative chemical analysis using UPLC-DAD and GNPS molecular networking. This strategy prompted study of Talaromyces sp. CMB-TN6F and Coccidiodes sp. CMB-TN39F, leading to discovery of the triterpene glycoside pullenvalenes A-D (1-4), featuring an unprecedented triterpene carbon skeleton and rare 6-O-methyl-N-acetyl-d-glucosaminyl glycoside residues. Structure elucidation of 1-4 was achieved by a combination of detailed spectroscopic analysis, chemical degradation, derivatization and synthesis, and biosynthetic considerations.

A new gene finding tool GeneMark-ETP significantly improves the accuracy of automatic annotation of large eukaryotic genomes

Large-scale genomic initiatives, such as the Earth BioGenome Project, require efficient methods for eukaryotic genome annotation. Here we present an automatic gene finder, GeneMark-ETP, integrating genomic-, transcriptomic- and protein-derived evidence that has been developed with a focus on large plant and animal genomes. GeneMark-ETP first identifies genomic loci where extrinsic data is sufficient for making gene predictions with 'high confidence'. The genes situated in the genomic space between the high confidence genes are predicted in the next stage. The set of high confidence genes serves as an initial training set for the statistical model. Further on, the model parameters are iteratively updated in the rounds of gene prediction and parameter re-estimation. Upon reaching convergence, GeneMark-ETP makes the final predictions and delivers the whole complement of predicted genes. GeneMark-ETP outperformed gene finders using a single type of extrinsic evidence. Comparisons with gene finders utilizing both transcript- and protein-derived extrinsic evidence, MAKER2, and TSEBRA, demonstrated that GeneMark-ETP delivered state-of-the-art gene prediction accuracy with the margin of outperforming existing approaches increasing in its applications to larger and more complex eukaryotic genomes.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

Genomic and Transcriptomic Comparisons of the Twig Blight Pathogen, Passalora sequoiae, with Mycosphaerellaceae Foliar and Conifer Pathogens

Passalora sequoiae is a foliar pathogen to conifer tree species. In this study, we conducted whole-genome and transcriptome analyses on isolates of P. sequoiae collected from symptomatic Leyland cypress leaves from a Christmas tree farm in Mississippi. The objectives for this research were to elucidate the pathogenicity mechanisms of P. sequoiae by characterizing the genome and transcriptome and possibly identify unique and shared predicted genes in comparison with non-conifer/canker and foliar pathogens in the family Mycosphaerellaceae. P. sequoiae was found to be similar to other foliar Mycosphaerellaceae pathogens and likely represents a hemibiotrophic lifestyle based on comparisons across pathogens. The genome and in planta transcriptome highlighted some unique features of P. sequoiae: the significant presence of chitin synthases and fructose-degrading carbohydrate-degrading enzymes, trans-AT PKS genes, and antibiotic gene clusters that were unique to P. sequoiae compared with the other Mycosphaerellaceae species genomes. Several transcripts that were highly expressed in planta were identified as effectors, yet the functions were not characterized. These targets provide ample resources to continue to characterize pathogen-conifer host interactions in conifer foliar pathogens. Furthermore, this research helps build genomic resources for an important plant pathogen on Leyland cypress that will further our ability to develop novel management practices that could begin with breeding for resistance.

Phylosymbiosis and metabolomics resolve phenotypically plastic and cryptic sponge species in the genus Agelas across the Caribbean basin

Fundamental to holobiont biology is recognising how variation in microbial composition and function relates to host phenotypic variation. Sponges often exhibit considerable phenotypic plasticity and also harbour dense microbial communities that function to protect and nourish hosts. One of the most prominent sponge genera on Caribbean coral reefs is Agelas. Using a comprehensive set of morphological (growth form, spicule), chemical and molecular data on 13 recognised species of Agelas in the Caribbean basin, we were able to define only five species (=clades) and found that many morphospecies designations were incongruent with phylogenomic and population genetic analyses. Microbial communities were also strongly differentiated between phylogenetic species, showing little evidence of cryptic divergence and relatively low correlation with morphospecies assignment. Metagenomic analyses also showed strong correspondence to phylogenetic species, and to a lesser extent, geographical and morphological characters. Surprisingly, the variation in secondary metabolites produced by sponge holobionts was explained by geography and morphospecies assignment, in addition to phylogenetic species, and covaried significantly with a subset of microbial symbionts. Spicule characteristics were highly plastic, under greater impact from geographical location than phylogeny. Our results suggest that while phenotypic plasticity is rampant in Agelas, morphological differences within phylogenetic species affect functionally important ecological traits, including the composition of the symbiotic microbial communities and metabolomic profiles.

Uncovering microbiomes of the rice phyllosphere using long-read metagenomic sequencing

The plant microbiome is crucial for plant growth, yet many important questions remain, such as the identification of specific bacterial species in plants, their genetic content, and location of these genes on chromosomes or plasmids. To gain insights into the genetic makeup of the rice-phyllosphere, we perform a metagenomic analysis using long-read sequences. Here, 1.8 Gb reads are assembled into 26,067 contigs including 142 circular sequences. Within these contigs, 669 complete 16S rRNA genes are clustered into 166 bacterial species, 121 of which show low identity (<97%) to defined sequences, suggesting novel species. The circular contigs contain novel chromosomes and a megaplasmid, and most of the smaller circular contigs are defined as novel plasmids or bacteriophages. One circular contig represents the complete chromosome of a difficult-to-culture bacterium Candidatus Saccharibacteria. Our findings demonstrate the efficacy of long-read-based metagenomics for profiling microbial communities and discovering novel sequences in plant-microbiome studies.

Genomics of the expanding pine pathogen Lecanosticta acicola reveals patterns of ongoing genetic admixture

Lecanosticta acicola is the causal agent for brown spot needle blight that affects pine trees across the northern hemisphere. Based on marker genes and microsatellite data, two distinct lineages have been identified that were introduced into Europe on two separate occasions. Despite their overall distinct geographic distribution, they have been found to coexist in regions of northern Spain and France. Here, we present the first genome-wide study of Lecanosticta acicola, including assembly of the reference genome and a population genomics analysis of 70 natural isolates from northern Spain. We show that most of the isolates belong to the southern lineage but show signs of introgression with northern lineage isolates, indicating mating between the two lineages. We also identify phenotypic differences between the two lineages based on the activity profiles of 20 enzymes, with introgressed strains being more phenotypically similar to members of the southern lineage. In conclusion, we show undergoing genetic admixture between the two main lineages of L. acicola in a region of recent expansion.

IMPORTANCE

Lecanosticta acicola is a fungal pathogen causing severe defoliation, growth reduction, and even death in more than 70 conifer species. Despite the increasing incidence of this species, little is known about its population dynamics. Two divergent lineages have been described that have now been found together in regions of France and Spain, but it is unknown how these mixed populations evolve. Here we present the first reference genome for this important plant pathogenic fungi and use it to study the population genomics of 70 isolates from an affected forest in the north of Spain. We find signs of introgression between the two main lineages, indicating that active mating is occurring in this region which could propitiate the appearance of novel traits in this species. We also study the phenotypic differences across this population based on enzymatic activities on 20 compounds.

Decryption of the survival "black box": gene family expansion promotes the encystment in ciliated protists

BACKGROUND

Encystment is an important survival strategy extensively employed by microbial organisms to survive unfavorable conditions. Single-celled ciliated protists (ciliates) are popular model eukaryotes for studying encystment, whereby these cells degenerate their ciliary structures and develop cyst walls, then reverse the process under more favorable conditions. However, to date, the evolutionary basis and mechanism for encystment in ciliates is largely unknown. With the rapid development of high-throughput sequencing technologies, genome sequencing and comparative genomics of ciliates have become effective methods to provide insights into above questions.

RESULTS

Here, we profiled the MAC genome of Pseudourostyla cristata, a model hypotrich ciliate for encystment studies. Like other hypotrich MAC genomes, the P. cristata MAC genome is extremely fragmented with a single gene on most chromosomes, and encodes introns that are generally small and lack a conserved branch point for pre-mRNA splicing. Gene family expansion analyses indicate that multiple gene families involved in the encystment are expanded during the evolution of P. cristata. Furthermore, genomic comparisons with other five representative hypotrichs indicate that gene families of phosphorelay sensor kinase, which play a role in the two-component signal transduction system that is related to encystment, show significant expansion among all six hypotrichs. Additionally, cyst wall-related chitin synthase genes have experienced structural changes that increase them from single-exon to multi-exon genes during evolution. These genomic features potentially promote the encystment in hypotrichs and enhance their ability to survive in adverse environments during evolution.

CONCLUSIONS

We systematically investigated the genomic structure of hypotrichs and key evolutionary phenomenon, gene family expansion, for encystment promotion in ciliates. In summary, our results provided insights into the evolutionary mechanism of encystment in ciliates.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

The identification of small molecule inhibitors with anthelmintic activities that target conserved proteins among ruminant gastrointestinal nematodes

Gastrointestinal nematode (GIN) infections are a major concern for the ruminant industry worldwide and result in significant production losses. Naturally occurring polyparasitism and increasing drug resistance that potentiate disease outcomes are observed among the most prevalent GINs of veterinary importance. Within the five major taxonomic clades, clade Va represents a group of GINs that predominantly affect the abomasum or small intestine of ruminants. However, the development of effective broad-spectrum anthelmintics against ruminant clade Va GINs has been challenged by a lack of comprehensive druggable genome resources. Here, we first assembled draft genomes for three clade Va species (Cooperia oncophora, Trichostrongylus colubriformis, and Ostertagia ostertagi) and compared them with closely related ruminant GINs. Genome-wide phylogenetic reconstruction showed a relationship among ruminant GINs structured by taxonomic classification. Orthogroup (OG) inference and functional enrichment analyses identified 220 clade Va-specific and Va-conserved OGs, enriched for functions related to cell cycle and cellular senescence. Further transcriptomic analysis identified 61 taxonomically and functionally conserved clade Va OGs that may function as drug targets for new broad-spectrum anthelmintics. Chemogenomic screening identified 11 compounds targeting homologs of these OGs, thus having potential anthelmintic activity. In in vitro phenotypic assays, three kinase inhibitors (digitoxigenin, K-252a, and staurosporine) exhibited broad-spectrum anthelmintic activities against clade Va GINs by obstructing the motility of exsheathed L3 (xL3) or molting of xL3 to L4. These results demonstrate valuable applications of the new ruminant GIN genomes in gaining better insights into their life cycles and offer a contemporary approach to discovering the next generation of anthelmintics.IMPORTANCEGastrointestinal nematode (GIN) infections in ruminants are caused by parasites that inhibit normal function in the digestive tract of cattle, sheep, and goats, thereby causing morbidity and mortality. Coinfection and increasing drug resistance to current therapeutic agents will continue to worsen disease outcomes and impose significant production losses on domestic livestock producers worldwide. In combination with ongoing therapeutic efforts, advancing the discovery of new drugs with novel modes of action is critical for better controlling GIN infections. The significance of this study is in assembling and characterizing new GIN genomes of Cooperia oncophora, Ostertagia ostertagi, and Trichostrongylus colubriformis for facilitating a multi-omics approach to identify novel, biologically conserved drug targets for five major GINs of veterinary importance. With this information, we were then able to demonstrate the potential of commercially available compounds as new anthelmintics.

Short macrocyclic peptides in sponge genomes

Sponges (Porifera) contain many peptide-specialized metabolites with potent biological activities and significant roles in shaping marine ecology. It is well established that symbiotic bacteria produce bioactive "sponge" peptides, both on the ribosome (RiPPs) and nonribosomally. Here, we demonstrate that sponges themselves also produce many bioactive macrocyclic peptides, such as phakellistatins and related proline-rich macrocyclic peptides (PRMPs). Using the Stylissa carteri sponge transcriptome, methods were developed to find sequences encoding 46 distinct RiPP-type core peptides, of which ten encoded previously identified PRMP sequences. With this basis set, the genome and transcriptome of the sponge Axinella corrugata was interrogated to find 35 PRMP precursor peptides encoding 31 unique core peptide sequences. At least 11 of these produced cyclic peptides that were present in the sponge and could be characterized by mass spectrometry, including stylissamides A-D and seven previously undescribed compounds. Precursor peptides were encoded in the A. corrugata genome, confirming their animal origin. The peptides contained signal peptide sequences and highly repetitive recognition sequence-core peptide elements with up to 25 PRMP copies in a single precursor. In comparison to sponges without PRMPs, PRMP sponges are incredibly enriched in potentially secreted polypeptides, with >23,000 individual signal peptide encoding genes found in a single transcriptome. The similarities between PRMP biosynthetic genes and neuropeptides in terms of their biosynthetic logic suggest a fundamental biology linked to circular peptides, possibly indicating a widespread and underappreciated diversity of signaling peptide post-translational modifications across the animal kingdom.

High-quality genome resource of a basidiomycetous yeast, Moesziomyces antarcticus isolate RS1, isolated from rice seed

Moesziomyces antarcticus (anamorph: Pseudozyma antarctica) is a basidiomycetous yeast in the Ustilaginaceae family and is a core member of the rice seed microbiome. M. antarcticus RS1 was isolated from surface-sterilized rice seeds. This 18.287 Mb draft genome of M. antarcticus RS1 is comprised of a 60.8% GC content and 6,817 protein-coding genes.

A Case of Okur-Chung Neurodevelopmental Syndrome with a Novel, de novo Variant on the CSNK2A1 Gene in a Turkish Patient

Introduction

Okur-Chung neurodevelopmental syndrome (OCNDS; #617062) has been associated with heterozygous mutations in the CSNK2A1 gene (*115440) mapped on the chromosome's 20p13 region.

Case Presentation

The analysis was performed on a 2-year-old patient who was admitted to our genetic diseases evaluation center by his family with a complaint of hypotonia. We detected a heterozygous NM_177559.3 (CSNK2A1):c.1139_1140dupGG (p.Met381GlyfsTer32) variant in the CSNK2A1 gene from a whole-exome sequence analysis.

Conclusion

The variant that we detected has not been reported in open-access databases to date, so it was evaluated as a novel likely pathogenic variant according to the ACMG-2015 criteria. No variant was detected upon segregation analysis of the patient's parents; therefore, the related variant was evaluated as de novo. In this study, we offer the first report of a pathogenic frameshift variant in the CSNK2A1 gene that has a relationship with OCNDS.

Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

The polyketide to fatty acid transition in the evolution of animal lipid metabolism

Animals synthesize simple lipids using a distinct fatty acid synthase (FAS) related to the type I polyketide synthase (PKS) enzymes that produce complex specialized metabolites. The evolutionary origin of the animal FAS and its relationship to the diversity of PKSs remain unclear despite the critical role of lipid synthesis in cellular metabolism. Recently, an animal FAS-like PKS (AFPK) was identified in sacoglossan molluscs. Here, we explore the phylogenetic distribution of AFPKs and other PKS and FAS enzymes across the tree of life. We found AFPKs widely distributed in arthropods and molluscs (>6300 newly described AFPK sequences). The AFPKs form a clade with the animal FAS, providing an evolutionary link bridging the type I PKSs and the animal FAS. We found molluscan AFPK diversification correlated with shell loss, suggesting AFPKs provide a chemical defense. Arthropods have few or no PKSs, but our results indicate AFPKs contributed to their ecological and evolutionary success by facilitating branched hydrocarbon and pheromone biosynthesis. Although animal metabolism is well studied, surprising new metabolic enzyme classes such as AFPKs await discovery.

Genomic exploration of the endangered oriental stork, Ciconia boyciana, sheds light on migration adaptation and future conservation

BACKGROUND

The oriental stork, Ciconia boyciana, is an endangered migratory bird listed on the International Union for Conservation of Nature's Red List. The bird population has experienced a rapid decline in the past decades, with nest locations and stop-over sites largely degraded due to human-bird conflicts. Multipronged conservation efforts are required to secure the future of oriental storks. We propose that a thorough understanding of the genome-wide genetic background of this threatened bird species is critical to make future conservation strategies.

FINDINGS

In this study, the first chromosome-scale reference genome was presented for the oriental stork with high quality, contiguity, and accuracy. The assembled genome size was 1.24 Gb with a scaffold N50 of 103 Mb, and 1.23 Gb contigs (99.32%) were anchored to 35 chromosomes. Population genomic analysis did not show a genetic structure in the wild population. Genome-wide genetic diversity (π = 0.0012) of the oriental stork was at a moderate to high level among threatened bird species, and the inbreeding risk was also not significant (FROH = 5.56% ± 5.30%). Reconstruction of demographic history indicated a rapid recent population decline likely driven by human activities. Genes that were under positive selection associated with the migratory trait were identified in relation to the long-term potentiation, photoreceptor cell organization, circadian rhythm, muscle development, and energy metabolism, indicating the essential interplay between genetic and ecological adaptation.

CONCLUSIONS

Our study presents the first chromosome-scale genome assembly of the oriental stork and provides a genomic basis for understanding a genetic background of the oriental stork, the population's extinction risks, and the migratory characteristics, which will facilitate the decision of future conservation plans for this species.

Comparative Genome Annotation

Newly sequenced genomes are being added to the tree of life at an unprecedented fast pace. A large proportion of such new genomes are phylogenetically close to previously sequenced and annotated genomes. In other cases, whole clades of closely related species or strains ought to be annotated simultaneously. Often, in subsequent studies, differences between the closely related species or strains are in the focus of research when the shared gene structures prevail. We here review methods for comparative structural genome annotation. The reviewed methods include classical approaches such as the alignment of protein sequences or protein profiles against the genome and comparative gene prediction methods that exploit a genome alignment to annotate either a single target genome or all input genomes simultaneously. We discuss how the methods depend on the phylogenetic placement of genomes, give advice on the choice of methods, and examine the consistency between gene structure annotations in an example. Furthermore, we provide practical advice on genome annotation in general.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes.

Chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder)

Lithocarpus, with >320 species, is the second largest genus of Fagaceae. However, the lack of a reference genome limits the molecular biology and functional study of Lithocarpus species. Here, we report the chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder), the first Lithocarpus species to be sequenced to date. Sweet tea has a 952-Mb genome, with a 21.4-Mb contig N50 value and 98.6% complete BUSCO score. In addition, the per-base consensus accuracy and completeness of the genome were estimated at 60.6 and 81.4, respectively. Genome annotation predicted 37,396 protein-coding genes, with repetitive sequences accounting for 64.2% of the genome. The genome did not undergo whole-genome duplication after the gamma (γ) hexaploidy event. Phylogenetic analysis showed that sweet tea diverged from the genus Quercus approximately at 59 million years ago. The high-quality genome assembly and gene annotation resources enrich the genomics of sweet tea, and will facilitate functional genomic studies in sweet tea and other Fagaceae species.