GeMoMa: Homology-Based Gene Prediction Utilizing Intron Position Conservation and RNA-seq Data

Chromosome-level genome assembly of the two-spotted spider mite Tetranychus urticae

The two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae), is a notorious pest in agriculture that has developed resistance to almost all chemical types used for its control. Here, we assembled a chromosome-level genome for the TSSM using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled contigs had a total length of 103.94 Mb with an N50 of 3.46 Mb, with 87.7 Mb of 34 contigs anchored to three chromosomes. The chromosome-level genome assembly had a BUSCO completeness of 94.8%. We identified 15,604 protein-coding genes, with 11,435 genes that could be functionally annotated. The high-quality genome provides invaluable resources for the genetic and evolutionary study of TSSM.

Draft genome assemblies of the ponerine ant Odontoponera transversa and the carpenter ant Camponotus friedae (Hymenoptera: Formicidae)

OBJECTIVES

Ants are ecologically dominant insects in most terrestrial ecosystems, with more than 14,000 extant species in about 340 genera recorded to date. However, genomic resources are still scarce for most species, especially for species endemic in East or Southeast Asia, limiting the study of phylogeny, speciation and adaptation of this evolutionarily successful animal lineage. Here, we assemble and annotate the genomes of Odontoponera transversa and Camponotus friedae, two ant species with a natural distribution in China, to facilitate future study of ant evolution.

DATA DESCRIPTION

We obtained a total of 16 Gb and 51 Gb PacBio HiFi data for O. transversa and C. friedae, respectively, which were assembled into the draft genomes of 339 Mb for O. transversa and 233 Mb for C. friedae. Genome assessments by multiple metrics showed good completeness and high accuracy of the two assemblies. Gene annotations assisted by RNA-seq data yielded a comparable number of protein-coding genes in the two genomes (10,892 for O. transversa and 11,296 for C. friedae), while repeat annotations revealed a remarkable difference of repeat content between these two ant species (149.4 Mb for O. transversa versus 49.7 Mb for C. friedae). Besides, complete mitochondrial genomes for the two species were assembled and annotated.

A high-quality chromosome-level genome assembly of the endangered tree Kmeria septentrionalis

Kmeria septentrionalis is a critically endangered tree endemic to Guangxi, China, and is listed on the International Union for Conservation of Nature's Red List. The lack of genetic information and high-quality genome data has hindered conservation efforts and studies on this species. In this study, we present a chromosome-level genome assembly of K. septentrionalis. The genome was initially assembled to be 2.57 Gb, with a contig N50 of 11.93 Mb. Hi-C guided genome assembly allowed us to anchor 98.83% of the total length of the initial contigs onto 19 pseudochromosomes, resulting in a scaffold N50 of 135.08 Mb. The final chromosome-level genome, spaning 2.54 Gb, achieved a BUSCO completeness of 98.9% and contained 1.67 Gb repetitive elements and 35,927 coding genes. This high-quality genome assembly provides a valuable resource for understanding the genetic basis of conservation-related traits and biological properties of this endangered tree species. Furthermore, it lays a critical foundation for evolutionary studies within the Magnoliaceae family.

A near complete genome assembly of the East Friesian sheep genome

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

The crane fly glycosylated triketide δ-lactone cornicinine elicits akinete differentiation of the cyanobiont in aquatic Azolla fern symbioses

The restriction of plant-symbiont dinitrogen fixation by an insect semiochemical had not been previously described. Here we report on a glycosylated triketide δ-lactone from Nephrotoma cornicina crane flies, cornicinine, that causes chlorosis in the floating-fern symbioses from the genus Azolla. Only the glycosylated trans-A form of chemically synthesized cornicinine was active: 500 nM cornicinine in the growth medium turned all cyanobacterial filaments from Nostoc azollae inside the host leaf-cavities into akinetes typically secreting CTB-bacteriocins. Cornicinine further inhibited akinete germination in Azolla sporelings, precluding re-establishment of the symbiosis during sexual reproduction. It did not impact development of the plant Arabidopsis thaliana or several free-living cyanobacteria from the genera Anabaena or Nostoc but affected the fern host without cyanobiont. Fern-host mRNA sequencing from isolated leaf cavities confirmed high NH4-assimilation and proanthocyanidin biosynthesis in this trichome-rich tissue. After cornicinine treatment, it revealed activation of Cullin-RING ubiquitin-ligase-pathways, known to mediate metabolite signaling and plant elicitation consistent with the chlorosis phenotype, and increased JA-oxidase, sulfate transport and exosome formation. The work begins to uncover molecular mechanisms of cyanobiont differentiation in a seed-free plant symbiosis important for wetland ecology or circular crop-production today, that once caused massive CO2 draw-down during the Eocene geological past.

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Identifying genes within pathways in unannotated genomes with PaGeSearch

In biological research, the identification and comparison of genes within specific pathways across the genomes of various species are invaluable. However, annotating the entire genome is resource intensive, and sequence similarity searches often yield results that are not actually genes. To address these limitations, we introduce Pathway Gene Search (PaGeSearch), a tool designed to identify genes from predefined lists, especially those in specific pathways, within genomes. The tool uses an initial sequence similarity search to identify relevant genomic regions, followed by targeted gene prediction and neural network-based result filtering. PaGeSearch suggests the regions that are most likely the orthologs of the genes in the query and is designed to be applicable for species within five classes: mammals, fish, birds, eudicotyledons, and Liliopsida. Compared with GeMoMa and miniprot, PaGeSearch generally outperforms in terms of sensitivity and positive predictive value, as well as negative predictive value. Also, the exon coverage of gene models from PaGeSearch is higher compared with those in GeMoMa and miniprot. Although its performance shows increased variability when applied to actual biological pathways, it nonetheless maintains an acceptable level of accuracy. Evaluating PaGeSearch across different assembly levels, chromosome, scaffold, and contig shows minimal variation in outcomes, indicating that PaGeSearch is resilient to variations in assembly quality.

Chromosome-level genome assembly of Euphorbia tirucalli (Euphorbiaceae), a highly stress-tolerant oil plant

Euphorbia, one of the largest genera of flowering plants, is well-known for containing many biofuel crops. Euphorbia tirucalli, an evergreen succulent mainly native to the Africa continent but cultivated worldwide, is a promising petroleum plant with high tolerance to drought and salt stress. However, the exploration of such an important plant resource is severely hampered by the lack of a reference genome. Here, we present the chromosome-level genome assembly of E. tirucalli using PacBio HiFi sequencing and Hi-C technology. Its genome size was approximately 745.62 Mb, with a contig N50 of 74.16 Mb. A total of 743.63 Mb (99.73%) of the assembled sequences were anchored to 10 chromosomes with a complete BUSCO score of 97.80%. Genome annotation revealed 26,304 protein-coding genes, and 76.37% of the genome was identified as repeat elements. The high-quality genome provides valuable genetic resources that would be useful for unraveling the genetic mechanisms of biofuel synthesis and evolutionary adaptation of E. tirucalli.

Chromosome-level genome assembly of ridgetail white shrimp Exopalaemon carinicauda

Exopalaemon carinicauda, a eurythermal and euryhaline shrimp, contributes one third of the total biomass production of polyculture ponds in eastern China and is considered as a potential ideal experimental animal for research on crustaceans. We conducted a high-quality chromosome-level genome assembly of E. carinicauda combining PacBio HiFi and Hi-C sequencing data. The total assembly size was 5.86 Gb, with a contig N50 of 235.52 kb and a scaffold N50 of 138.24 Mb. Approximately 95.29% of the assembled sequences were anchored onto 45 pseudochromosomes. BUSCO analysis revealed that 92.89% of 1,013 single-copy genes were highly conserved orthologs. A total of 44, 288 protein-coding genes were predicted, of which 70.53% were functionally annotated. Given its high heterozygosity (2.62%) and large proportion of repeat sequences (71.49%), it is one of the most complex genome assemblies. This chromosome-scale genome will be a valuable resource for future molecular breeding and functional genomics research on E. carinicauda.

Genome-wide sequencing reveals geographical variations in the thermal adaptation of an aquaculture species with frequent seedling introductions

Climate change and human activity are essential factors affecting marine biodiversity and aquaculture, and understanding the impacts of human activities on the genetic structure to increasing high temperatures is crucial for sustainable aquaculture and marine biodiversity conservation. As a commercially important bivalve, the Manila clam Ruditapes philippinarum is widely distributed along the coast of China, and it has been frequently introduced from Fujian Province, China, to other regions for aquaculture. In this study, we collected four populations of Manila clams from different areas to evaluate their thermal tolerance by measuring cardiac performance and genetic variations using whole-genome resequencing. The upper thermal limits of the clams showed high variations within and among populations. Different populations displayed divergent genetic compositions, and the admixed population was partly derived from the Zhangzhou population in Fujian Province, implying a complex genomic landscape under the influence of local genetic sources and human introductions. Multiple single nucleotide polymorphisms (SNPs) were associated with the cardiac functional traits, and some of these SNPs can affect the codon usage and the structural stability of the resulting protein. This study shed light on the importance of establishing long-term ecological and genetic monitoring programs at the local level to enhance resilience to future climate change.

Metabolite-based genome-wide association studies enable the dissection of the genetic bases of flavonoids, betaine and spermidine in wolfberry (Lycium)

Wolfberry is a plant with medicinal and food values. However, its bioactive ingredients and the corresponding genetic bases have not been determined. Here, we de novo generated a chromosome-level genome assembly for wolfberry, yielding a genome sequence of ~1.77 Gb with contig N50 of 50.55 Mb and 39 224 predicted gene models. A variation map, using 307 re-sequenced accessions, was called based on this genome assembly. Furthermore, the fruit metabolome of these accessions was profiled using 563 annotated metabolites, which separated Lycium barbarum L. and non-L. barbarum L. The flavonoids, coumarins, alkaloids and nicotinic acid contents were higher in the former than in the latter. A metabolite-based genome-wide association study mapped 156 164 significant single nucleotide polymorphisms corresponding to 340 metabolites. This included 19 219 unique lead single nucleotide polymorphisms in 1517 significant association loci, of which three metabolites, flavonoids, betaine and spermidine, were highlighted. Two candidate genes, LbUGT (evm.TU.chr07.2692) and LbCHS (evm.TU.chr07.2738), with non-synonymous mutations, were associated with the flavonoids content. LbCHS is a structural gene that interacts with a nearby MYB transcription factor (evm.TU.chr07.2726) both in L. barbarum and L. ruthenicum. Thus, these three genes might be involved in the biosynthesis/metabolism of flavonoids. LbSSADH (evm.TU.chr09.627) was identified as possibly participating in betaine biosynthesis/metabolism. Four lycibarbarspermidines (E-G and O) were identified, and only the lycibarbarspermidines O content was higher in L. barbarum varieties than in non-L. barbarum varieties. The evm.TU.chr07.2680 gene associated with lycibarbarspermidines O was annotated as an acetyl-CoA-benzylalcohol acetyltransferase, suggesting that it is a candidate gene for spermidine biosynthesis. These results provide novel insights into the specific metabolite profile of non-L. barbarum L. and the genetic bases of flavonoids, betaine and spermidine biosynthesis/metabolism.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

What one genus of showy moths can say about migration, adaptation, and wing pattern

The Ornate Moth, Utetheisa ornatrix, has served as a model species in chemical ecology studies for decades. Like in the widely publicized stories of the Monarch and other milkweed butterflies, the Ornate Moth and its relatives are tropical insects colonizing whole continents assisted by their chemical defenses. With the recent advances in genomic techniques and evo-devo research, it is becoming a model for studies in other areas, from wing pattern development to phylogeography, from toxicology to epigenetics. We used a genomic approach to learn about Utetheisa's evolution, detoxification, dispersal abilities, and wing pattern diversity. We present an evolutionary genomic analysis of the worldwide genus Utetheisa, then focusing on U. ornatrix. Our reference genome of U. ornatrix reveals gene duplications in the regions possibly associated with detoxification abilities, which allows them to feed on toxic food plants. Finally, comparative genomic analysis of over 100 U. ornatrix specimens from the museum with apparent differences in wing patterns suggest the potential roles of cortex and lim3 genes in wing pattern formation of Lepidoptera and the utility of museum-preserved collection specimens for wing pattern research.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899)

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Chromosome-scale genome assemblies of Himalopsyche anomala and Eubasilissa splendida (Insecta: Trichoptera)

Trichoptera is one of the most evolutionarily successful aquatic insect lineages and is highly valued value in adaptive evolution research. This study presents the chromosome-level genome assemblies of Himalopsyche anomala and Eubasilissa splendida achieved using PacBio, Illumina, and Hi-C sequencing. For H. anomala and E. splendida, assembly sizes were 663.43 and 859.28 Mb, with scaffold N50 lengths of 28.44 and 31.17 Mb, respectively. In H. anomala and E. splendida, we anchored 24 and 29 pseudochromosomes, and identified 11,469 and 10,554 protein-coding genes, respectively. The high-quality genomes of H. anomala and E. splendida provide critical genomic resources for understanding the evolution and ecology of Trichoptera and performing comparative genomics analyses.

The genome of Przewalski's horse (Equus ferus przewalskii)

The Przewalski's horse (Equus ferus przewalskii) is an endangered equid native to the steppes of central Asia. After becoming extinct in the wild, multiple conservation efforts convened to preserve the species including captive breeding programs, reintroduction and monitoring systems, protected lands, and cloning. Availability of a highly contiguous reference genome is essential to support these continued efforts. We used Oxford Nanopore sequencing to produce a scaffold-level 2.5 Gb nuclear assembly and 16,002 bp mitogenome from a captive Przewalski's mare. All assembly drafts were generated from 111 Gb of sequence from a single PromethION R10.4.1 flow cell. The mitogenome contained 37 genes in the standard mammalian configuration and was 99.63% identical to the domestic horse (Equus caballus). The nuclear assembly, EquPr2, contained 2,146 scaffolds with an N50 of 85.1 Mb, 43X mean depth, and BUSCO quality score of 98.92%. EquPr2 successfully improves upon the existing Przewalski's horse reference genome (Burgud), with 25-fold fewer scaffolds, a 166-fold larger N50, and phased pseudohaplotypes. Modified basecalls revealed 79.5% DNA methylation and 2.1% hydroxymethylation globally. Allele-specific methylation analysis between pseudohaplotypes revealed 226 differentially methylated regions (DMRs) in known imprinted genes and loci not previously reported as imprinted. The heterozygosity rate of 0.165% matches previous estimates for the species and compares favorably to other endangered animals. This improved Przewalski's horse assembly will serve as a valuable resource for conservation efforts and comparative genomics investigations.

Genomic insights into the evolution of flavonoid biosynthesis and O-methyltransferase and glucosyltransferase in Chrysanthemum indicum

Flavonoids are a class of secondary metabolites widely distributed in plants. Regiospecific modification by methylation and glycosylation determines flavonoid diversity. A rare flavone glycoside, diosmin (luteolin-4'-methoxyl-7-O-glucosyl-rhamnoside), occurs in Chrysanthemum indicum. How Chrysanthemum plants evolve new biosynthetic capacities remains elusive. Here, we assemble a 3.11-Gb high-quality C. indicum genome with a contig N50 value of 4.39 Mb and annotate 50,606 protein-coding genes. One (CiCOMT10) of the tandemly repeated O-methyltransferase genes undergoes neofunctionalization, preferentially transferring the methyl group to the 4'-hydroxyl group of luteolin with ortho-substituents to form diosmetin. In addition, CiUGT11 (UGT88B3) specifically glucosylates 7-OH group of diosmetin. Next, we construct a one-pot cascade biocatalyst system by combining CiCOMT10, CiUGT11, and our previously identified rhamnosyltransferase, effectively producing diosmin with over 80% conversion from luteolin. This study clarifies the role of transferases in flavonoid diversity and provides important gene elements essential for producing rare flavone.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

Pedigree genome data of an early-matured Geng/japonica glutinous rice mega variety Longgeng 57

By using PacBio HiFi technology, we produced over 700 Gb of long-read sequencing (LRS) raw data; and by using Illumina paired-end whole-genome shotgun (WGS) sequencing technology, we generated more than 70 Gb of short-read sequencing (SRS) data. With LRS data, we assembled one genome and then generate a set of annotation data for an early-matured Geng/japonica glutinous rice mega variety genome, Longgeng 57 (LG57), which carries multiple elite traits including good grain quality and wide adaptability. Together with the SRS data from three parents of LG57, pedigree genome variations were called for three representative types of genes. These data sets can be used for deep variation mining, aid in the discovery of new insights into genome structure, function, and evolution, and help to provide essential support to biological research in general.

Genomic and Comparative Transcriptomic Analyses Reveal Key Genes Associated with the Biosynthesis Regulation of Okaramine B in Penicillium daleae NBP-49626

Restricted production of fungal secondary metabolites hinders the ability to conduct comprehensive research and development of novel biopesticides. Okaramine B from Penicillium demonstrates remarkable insecticidal efficacy; however, its biosynthetic yield is low, and its regulatory mechanism remains unknown. The present study found that the yield difference was influenced by fermentation modes in okaramine-producing strains and performed genomic and comparative transcriptome analysis of P. daleae strain NBP-49626, which exhibits significant features. The NBP-49626 genome is 37.4 Mb, and it encodes 10,131 protein-encoding genes. Up to 5097 differentially expressed genes (DEGs) were identified during the submerged and semi-solid fermentation processes. The oka gene cluster, lacking regulatory and transport genes, displayed distinct transcriptional patterns in response to the fermentation modes and yield of Okaramine B. Although transcription trends of most known global regulatory genes are inconsistent with those of oka, this study identified five potential regulatory genes, including two novel Zn(II)2Cys6 transcription factors, Reg2 and Reg19. A significant correlation was also observed between tryptophan metabolism and Okaramine B yields. In addition, several transporter genes were identified as DEGs. These results were confirmed using real-time quantitative PCR. This study provides comprehensive information regarding the regulatory mechanism of Okaramine B biosynthesis in Penicillium and is critical to the further yield improvement for the development of insecticides.

Two chromosome-level genomes of Smittia aterrima and Smittia pratorum (Diptera, Chironomidae)

Chironomids are one of the most abundant aquatic insects and are widely distributed in various biological communities. However, the lack of high-quality genomes has hindered our ability to study the evolution and ecology of this group. Here, we used Nanopore long reads and Hi-C data to produce two chromosome-level genomes from mixed genomic data. The genomes of Smittia aterrima (SateA) and Smittia pratorum (SateB) were assembled into three chromosomes, with sizes of 78.45 Mb and 71.56 Mb, scaffold N50 lengths of 25.73 and 23.53 Mb, and BUSCO completeness of 98.5% and 97.8% (n = 1,367), 5.68 Mb (7.24%) and 1.94 Mb (2.72%) of repetitive elements, and predicted 12,330 (97.70% BUSCO completeness) and 11,250 (97.40%) protein-coding genes, respectively. These high-quality genomes will serve as valuable resources for comprehending the evolution and environmental adaptation of chironomids.

De novo Phased Genome Assembly, Annotation and Population Genotyping of Alectoris Chukar

The Alectoris Chukar (chukar) is the most geographically widespread partridge species in the world, demonstrating exceptional adaptability to diverse ecological environments. However, the scarcity of genetic resources for chukar has hindered research into its adaptive evolution and molecular breeding. In this study, we have sequenced and assembled a high-quality, phased chukar genome that consists of 31 pairs of relatively complete diploid chromosomes. Our BUSCO analysis reported a high completeness score of 96.8% and 96.5%, with respect to universal single-copy orthologs and a low duplication rate (0.3% and 0.5%) for two assemblies. Through resequencing and population genomic analyses of six subspecies, we have curated invaluable genotype data that underscores the adaptive evolution of chukar in response to both arid and high-altitude environments. These data will significantly contribute to research on how chukars adaptively evolve to cope with desertification and alpine climates.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Chromosome-scale Genome assembly of the critically endangered White-eared Night-Heron (Gorsachius magnificus)

The White-eared Night-Heron (Gorsachius magnificus, G. magnificus) is a critically endangered heron that is very poorly known and only found in southern China and northern Vietnam, with an estimated population of 250 to 999 mature individuals. However, the lack of a reference genome has hindered the implementation of conservation management efforts. In this study, we present the first high-quality chromosome-scale reference genome, which was assembled by integrating PacBio long-reads sequencing, Illumina paired-end sequencing, and Hi-C technology. The genome has a total length of 1.176 Gb, with a scaffold N50 of 84.77 Mb and a contig N50 of 18.46 Mb. Utilizing Hi-C data, we anchored 99.89% of the scaffold sequences onto 29 pairs of chromosomes. Additionally, we identified 18,062 protein-coding genes in the genome, with 95.00% of which were functionally annotated. Notably, BUSCO assessment confirmed the presence of 97.2% of highly conserved Aves genes within the genome. This chromosome-level genome assembly and annotation will be valuable for future investigating the G. magnificus's evolutionary adaptation and conservation.

Gene modelling and annotation for the Hawaiian bobtail squid, Euprymna scolopes

Coleoid cephalopods possess numerous complex, species-specific morphological and behavioural adaptations, e.g., a uniquely structured nervous system that is the largest among the invertebrates. The Hawaiian bobtail squid (Euprymna scolopes) is one of the most established cephalopod species. With its recent publication of the chromosomal-scale genome assembly and regulatory genomic data, it also emerges as a key model for cephalopod gene regulation and evolution. However, the latest genome assembly has been lacking a native gene model set. Our manuscript describes the generation of new long-read transcriptomic data and, made using this combined with a plethora of publicly available transcriptomic and protein sequence data, a new reference annotation for E. scolopes.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

Telomere-to-Telomere Genome Assembly of Tibetan Medicinal Mushroom Ganoderma leucocontextum and the First Copia Centromeric Retrotransposon in Macro-Fungi Genome

A complete telomere-to-telomere (T2T) genome has been a longstanding goal in the field of genomic research. By integrating high-coverage and precise long-read sequencing data using multiple assembly strategies, we present here the first T2T gap-free genome assembly of Ganoderma leucocontextum strain GL72, a Tibetan medicinal mushroom. The T2T genome, with a size of 46.69 Mb, consists 13 complete nuclear chromosomes and typical telomeric repeats (CCCTAA)n were detected at both ends of 13 chromosomes. The high mapping rate, uniform genome coverage, a complete BUSCOs of 99.7%, and base accuracy exceeding 99.999% indicate that this assembly represents the highest level of completeness and quality. Regions characterized by distinct structural attributes, including highest Hi-C interaction intensity, high repeat content, decreased gene density, low GC content, and minimal or no transcription levels across all chromosomes may represent potential centromeres. Sequence analysis revealed the first Copia centromeric retrotransposon in macro-fungi genome. Phylogenomic analysis identified that G. leucocontextum and G. tsugae diverged from the other Ganoderma species approximately 9.8-17.9 MYA. The prediction of secondary metabolic clusters confirmed the capability of this fungus to produce a substantial quantity of metabolites. This T2T gap-free genome will contribute to the genomic 'dark matter' elucidation and server as a great reference for genetics, genomics, and evolutionary studies of G. leucocontextum.

Polyploid Genome Assembly Provides Insights into Morphological Development and Ascorbic Acid Accumulation of Sauropus androgynus

Sauropus androgynus (S. androgynus) (2n = 4x = 52) is one of the most popular functional leafy vegetables in South and Southeast Asia. With its rich nutritional and pharmaceutical values, it has traditionally had widespread use for dietary and herbal purposes. Here, the genome of S. androgynus was sequenced and assembled, revealing a genome size of 1.55 Gb with 26 pseudo-chromosomes. Phylogenetic analysis traced back the divergence of Sauropus from Phyllanthus to approximately 29.67 million years ago (Mya). Genome analysis revealed that S. androgynus polyploidized around 20.51 Mya and shared a γ event about 132.95 Mya. Gene function analysis suggested that the expansion of pathways related to phloem development, lignin biosynthesis, and photosynthesis tended to result in the morphological differences among species within the Phyllanthaceae family, characterized by varying ploidy levels. The high accumulation of ascorbic acid in S. androgynus was attributed to the high expression of genes associated with the L-galactose pathway and recycling pathway. Moreover, the expanded gene families of S. androgynus exhibited multiple biochemical pathways associated with its comprehensive pharmacological activity, geographic adaptation and distinctive pleasurable flavor. Altogether, our findings represent a crucial genomic asset for S. androgynus, casting light on the intricate ploidy within the Phyllanthaceae family.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

A genome assembly for Orinus kokonorica provides insights into the origin, adaptive evolution and further diversification of two closely related grass genera

Ancient whole-genome duplication (WGD) or polyploidization is prevalent in plants and has played a crucial role in plant adaptation. However, the underlying genomic basis of ecological adaptation and subsequent diversification after WGD are still poorly understood in most plants. Here, we report a chromosome-scale genome assembly for the genus Orinus (Orinus kokonorica as representative) and preform comparative genomics with its closely related genus Cleistogenes (Cleistogenes songorica as representative), both belonging to a newly named subtribe Orininae of the grass subfamily Chloridoideae. The two genera may share one paleo-allotetraploidy event before 10 million years ago, and the two subgenomes of O. kokonorica display neither fractionation bias nor global homoeolog expression dominance. We find substantial genome rearrangements and extensive structural variations (SVs) between the two species. With comparative transcriptomics, we demonstrate that functional innovations of orthologous genes may have played an important role in promoting adaptive evolution and diversification of the two genera after polyploidization. In addition, copy number variations and extensive SVs between orthologs of flower and rhizome related genes may contribute to the morphological differences between the two genera. Our results provide new insights into the adaptive evolution and subsequent diversification of the two genera after polyploidization.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Gene duplications facilitate C4-CAM compatibility in common purslane

Common purslane (Portulaca oleracea) integrates both C4 and crassulacean acid metabolism (CAM) photosynthesis pathways and is a promising model plant to explore C4-CAM plasticity. Here, we report a high-quality chromosome-level genome of nicotinamide adenine dinucleotide (NAD)-malic enzyme (ME) subtype common purslane that provides evidence for 2 rounds of whole-genome duplication (WGD) with an ancient WGD (P-β) in the common ancestor to Portulacaceae and Cactaceae around 66.30 million years ago (Mya) and another (Po-α) specific to common purslane lineage around 7.74 Mya. A larger number of gene copies encoding key enzymes/transporters involved in C4 and CAM pathways were detected in common purslane than in related species. Phylogeny, conserved functional site, and collinearity analyses revealed that the Po-α WGD produced the phosphoenolpyruvate carboxylase-encoded gene copies used for photosynthesis in common purslane, while the P-β WGD event produced 2 ancestral genes of functionally differentiated (C4- and CAM-specific) beta carbonic anhydrases involved in the C4 + CAM pathways. Additionally, cis-element enrichment analysis in the promoters showed that CAM-specific genes have recruited both evening and midnight circadian elements as well as the Abscisic acid (ABA)-independent regulatory module mediated by ethylene-response factor cis-elements. Overall, this study provides insights into the origin and evolutionary process of C4 and CAM pathways in common purslane, as well as potential targets for engineering crops by integrating C4 or CAM metabolism.

Genome Assembly of Cordia subcordata, a Coastal Protection Species in Tropical Coral Islands

Cordia subcordata trees or shrubs, belonging to the Boraginaceae family, have strong resistance and have adapted to their habitat on a tropical coral island in China, but the lack of genome information regarding its genetic background is unclear. In this study, the genome was assembled using both short/long whole genome sequencing reads and Hi-C reads. The assembled genome was 475.3 Mb, with 468.7 Mb (99.22%) of the sequences assembled into 16 chromosomes. Repeat sequences accounted for 54.41% of the assembled genome. A total of 26,615 genes were predicted, and 25,730 genes were functionally annotated using different annotation databases. Based on its genome and the other 17 species, phylogenetic analysis using 336 single-copy genes obtained from ortholog analysis showed that C. subcordata was a sister to Coffea eugenioides, and the divergence time was estimated to be 77 MYA between the two species. Gene family evolution analysis indicated that the significantly expanded gene families were functionally related to chemical defenses against diseases. These results can provide a reference to a deeper understanding of the genetic background of C. subcordata and can be helpful in exploring its adaptation mechanism on tropical coral islands in the future.

Gap-free genome assembly of Salangid icefish Neosalanx taihuensis

Neosalanx taihuensis is widely distributed in freshwater and brackish water areas in China. Due to its high commercial value, it has been artificially introduced into many lakes and reservoirs, showing strong ecological adaptability. Here, a gap-free chromosome-level reference genome was constructed by combining short reads, PacBio HiFi long reads, Nanopore ultralong reads and Hi-C data. The reference genome of N. taihuensis was 397.29 Mb with a contig N50 of 15.61 Mb. The assembled sequences were anchored to 28 chromosomes. Furthermore, 20,024 protein-coding genes and 98.16% of the predicted genes were annotated in publicly available biological databases. This high-quality gap-free assembled genome will provide an essential reference for studying the evolution and ecological adaptability of N. taihuensis.

A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa)

As plentiful high-quality genome assemblies have been accumulated, reference-guided genome assembly can be a good approach to reconstruct a high-quality assembly. Here, we present a chromosome-level genome assembly of the Korean crossbred pig called Nanchukmacdon (the NCMD assembly) using the reference-guided assembly approach with short and long reads. The NCMD assembly contains 20 chromosome-level scaffolds with a total size of 2.38 Gbp (N50: 138.77 Mbp). Its BUSCO score is 93.1%, which is comparable to the pig reference assembly, and a total of 20,588 protein-coding genes, 8,651 non-coding genes, and 996.14 Mbp of repetitive elements are annotated. The NCMD assembly was also used to close many gaps in the pig reference assembly. This NCMD assembly and annotation provide foundational resources for the genomic analyses of pig and related species.

Whole Genome Assembly and Annotation of Blackstripe Livebearer Poeciliopsis prolifica

The blackstripe livebearer Poeciliopsis prolifica is a live-bearing fish belonging to the family Poeciliidae with high level of postfertilization maternal investment (matrotrophy). This viviparous matrotrophic species has evolved a structure similarly to the mammalian placenta. Placentas have independently evolved multiple times in Poeciliidae from nonplacental ancestors, which provide an opportunity to study the placental evolution. However, there is a lack of high-quality reference genomes for the placental species in Poeciliidae. In this study, we present a 674 Mb assembly of P. prolifica in 504 contigs with excellent continuity (contig N50 7.7 Mb) and completeness (97.2% Benchmarking Universal Single-Copy Orthologs [BUSCO] completeness score, including 92.6% single-copy and 4.6% duplicated BUSCO score). A total of 27,227 protein-coding genes were annotated from the merged datasets based on bioinformatic prediction, RNA sequencing and homology evidence. Phylogenomic analyses revealed that P. prolifica diverged from the guppy (Poecilia reticulata) ∼19 Ma. Our research provides the necessary resources and the genomic toolkit for investigating the genetic underpinning of placentation.

Population genomics identifies genetic signatures of carrot domestication and improvement and uncovers the origin of high-carotenoid orange carrots

Here an improved carrot reference genome and resequencing of 630 carrot accessions were used to investigate carrot domestication and improvement. The study demonstrated that carrot was domesticated during the Early Middle Ages in the region spanning western Asia to central Asia, and orange carrot was selected during the Renaissance period, probably in western Europe. A progressive reduction of genetic diversity accompanied this process. Genes controlling circadian clock/flowering and carotenoid accumulation were under selection during domestication and improvement. Three recessive genes, at the REC, Or and Y2 quantitative trait loci, were essential to select for the high α- and β-carotene orange phenotype. All three genes control high α- and β-carotene accumulation through molecular mechanisms that regulate the interactions between the carotenoid biosynthetic pathway, the photosynthetic system and chloroplast biogenesis. Overall, this study elucidated carrot domestication and breeding history and carotenoid genetics at a molecular level.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

The haplotype-resolved T2T reference genome highlights structural variation underlying agronomic traits of melon

Melon (Cucumis melo L.) is an important vegetable crop that has an extensive history of cultivation. However, the genome of wild and semi-wild melon types that can be used for the analysis of agronomic traits is not yet available. Here we report a chromosome-level T2T genome assembly for 821 (C. melo ssp. agrestis var. acidulus), a semi-wild melon with two haplotypes of ~373 Mb and ~364 Mb, respectively. Comparative genome analysis discovered a significant number of structural variants (SVs) between melo (C. melo ssp. melo) and agrestis (C. melo ssp. agrestis) genomes, including a copy number variation located in the ToLCNDV resistance locus on chromosome 11. Genome-wide association studies detected a significant signal associated with climacteric ripening and identified one candidate gene CM_ac12g14720.1 (CmABA2), encoding a cytoplasmic short chain dehydrogenase/reductase, which controls the biosynthesis of abscisic acid. This study provides valuable genetic resources for future research on melon breeding.

An integrated transcriptomic and metabolic phenotype analysis to uncover the metabolic characteristics of a genetically engineered Candida utilis strain expressing δ-zein gene

Introduction

Candida utilis (C. utilis) has been extensively utilized as human food or animal feed additives. With its ability to support heterologous gene expression, C. utilis proves to be a valuable platform for the synthesis of proteins and metabolites that possess both high nutritional and economic value. However, there remains a dearth of research focused on the characteristics of C. utilis through genomic, transcriptomic and metabolic approaches.

Methods

With the aim of unraveling the molecular mechanism and genetic basis governing the biological process of C. utilis, we embarked on a de novo sequencing endeavor to acquire comprehensive sequence data. In addition, an integrated transcriptomic and metabolic phenotype analysis was performed to compare the wild-type C. utilis (WT) with a genetically engineered strain of C. utilis that harbors the heterologous δ-zein gene (RCT).

Results

δ-zein is a protein rich in methionine found in the endosperm of maize. The integrated analysis of transcriptomic and metabolic phenotypes uncovered significant metabolic diversity between the WT and RCT C. utilis. A total of 252 differentially expressed genes were identified, primarily associated with ribosome function, peroxisome activity, arginine and proline metabolism, carbon metabolism, and fatty acid degradation. In the experimental setup using PM1, PM2, and PM4 plates, a total of 284 growth conditions were tested. A comparison between the WT and RCT C. utilis demonstrated significant increases in the utilization of certain carbon source substrates by RCT. Gelatin and glycogen were found to be significantly utilized to a greater extent by RCT compared to WT. Additionally, in terms of sulfur source substrates, RCT exhibited significantly increased utilization of O-Phospho-L-Tyrosine and L-Methionine Sulfone when compared to WT.

Discussion

The introduction of δ-zein gene into C. utilis may lead to significant changes in the metabolic substrates and metabolic pathways, but does not weaken the activity of the strain. Our study provides new insights into the transcriptomic and metabolic characteristics of the genetically engineered C. utilis strain harboring δ-zein gene, which has the potential to advance the utilization of C. utilis as an efficient protein feed in agricultural applications.

Multi-omics data provide insight into the adaptation of the glasshouse plant Rheum nobile to the alpine subnival zone

Subnival glasshouse plants provide a text-book example of high-altitude adaptation with reproductive organs enclosed in specialized semi-translucent bracts, monocarpic reproduction and continuous survival under stress. Here, we present genomic, transcriptomic and metabolomic analyses for one such plant, the Noble rhubarb (Rheum nobile). Comparative genomic analyses show that an expanded number of genes and retained genes from two recent whole-genome duplication events are both relevant to subnival adaptation of this species. Most photosynthesis genes are downregulated within bracts compared to within leaves, and indeed bracts exhibit a sharp reduction in photosynthetic pigments, indicating that the bracts no longer perform photosynthesis. Contrastingly, genes related to flavonol synthesis are upregulated, providing enhanced defense against UV irradiation damage. Additionally, anatomically abnormal mesophyll combined with the downregulation of genes related to mesophyll differentiation in bracts illustrates the innovation and specification of the glass-like bracts. We further detect substantial accumulation of antifreeze proteins (e.g. AFPs, LEAs) and various metabolites (e.g. Proline, Protective sugars, procyanidins) in over-wintering roots. These findings provide new insights into subnival adaptation and the evolution of glasshouse alpine plants.