Genome Annotation and Curation Using MAKER and MAKER-P

Several secondary metabolite gene clusters in the genomes of ten Penicillium spp. raise the risk of multiple mycotoxin occurrence in chestnuts

Penicillium spp. produce a great variety of secondary metabolites, including several mycotoxins, on food substrates. Chestnuts represent a favorable substrate for Penicillium spp. development. In this study, the genomes of ten Penicillium species, virulent on chestnuts, were sequenced and annotated: P. bialowiezense. P. pancosmium, P. manginii, P. discolor, P. crustosum, P. palitans, P. viridicatum, P. glandicola, P. taurinense and P. terrarumae. Assembly size ranges from 27.5 to 36.8 Mb and the number of encoded genes ranges from 9,867 to 12,520. The total number of predicted biosynthetic gene clusters (BGCs) in the ten species is 551. The most represented families of BGCs are non ribosomal peptide synthase (191) and polyketide synthase (175), followed by terpene synthases (87). Genome-wide collections of gene phylogenies (phylomes) were reconstructed for each of the newly sequenced Penicillium species allowing for the prediction of orthologous relationships among our species, as well as other 20 annotated Penicillium species available in the public domain. We investigated in silico the presence of BGCs for 10 secondary metabolites, including 5 mycotoxins, whose production was validated in vivo through chemical analyses. Among the clusters present in this set of species we found andrastin A and its related cluster atlantinone A, mycophenolic acid, patulin, penitrem A and the cluster responsible for the synthesis of roquefortine C/glandicoline A/glandicoline B/meleagrin. We confirmed the presence of these clusters in several of the Penicillium species conforming our dataset and verified their capacity to synthesize them in a chestnut-based medium with chemical analysis. Interestingly, we identified mycotoxin clusters in some species for the first time, such as the andrastin A cluster in P. flavigenum and P. taurinense, and the roquefortine C cluster in P. nalgiovense and P. taurinense. Chestnuts proved to be an optimal substrate for species of Penicillium with different mycotoxigenic potential, opening the door to risks related to the occurrence of multiple mycotoxins in the same food matrix.

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

Draft genome sequence of Penicillium citrinum B9, a plant growth-promoting symbiont from barley rhizosphere

Penicillium citrinum strain B9 is a plant growth-promoting fungus isolated from Barley (Hordeum vulgare) rhizosphere. We report the first draft genome of P. citrinum B9 assembled using single-molecule real-time sequencing and Illumina reads. The assembled genome spans 31.3 Mb comprising nine contigs and 10,106 protein-encoding genes.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

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 haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Aureococcus anophagefferens strain CCMP1851: draft genome of a second Kratosvirus quantuckense-susceptible host strain for an emerging host-giant virus model system

Here, we report the draft genome of Aureococcus anophagefferens strain CCMP1851, which is susceptible to the virus Kratosvirus quantuckense. CCMP1851 complements an available genome for a virus-resistant strain (CCMP1850) isolated from the same bloom. Future studies can now use this genome to examine genetic hints of virus resistance and susceptibility.

Phased chromosome-scale genome assembly of an asexual, allopolyploid root-knot nematode reveals complex subgenomic structure

We present the chromosome-scale genome assembly of the allopolyploid root-knot nematode Meloidogyne javanica. We show that the M. javanica genome is predominantly allotetraploid, comprising two subgenomes, A and B, that most likely originated from hybridisation of two ancestral parental species. The assembly was annotated using full-length non-chimeric transcripts, comparison to reference databases, and ab initio prediction techniques, and the subgenomes were phased using ancestral k-mer spectral analysis. Subgenome B appears to show fission of chromosomal contigs, and while there is substantial synteny between subgenomes, we also identified regions lacking synteny that may have diverged in the ancestral genomes prior to or following hybridisation. This annotated and phased genome assembly forms a significant resource for understanding the origins and genetics of these globally important plant pathogens.

Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Genome assembly and transcriptomic analysis to elucidate the ability of Nasonovia ribisnigri to break host plant resistance

Aphid genomic resources enable the study of complex life history traits and provide information on vector biology, host adaption and speciation. The currant-lettuce aphid (Nasonovia ribisnigri (Hemiptera: Aphididae) (Mosley)) is a cosmopolitan pest of outdoor lettuce (Lactuca sativa (Asterales: Asteraceae) (Linnaeus)). Until recently, the use of resistant cultivars was an effective method for managing N. ribisnigri. A resistant cultivar containing a single gene (Nr-locus), introduced in the 1980s, conferred complete resistance to feeding. Overreliance of this Nr-locus in lettuce resulted in N. ribisnigri's ability to break resistance mechanism, with first reports during 2003. Our work attempts to understand which candidate gene(s) are associated with this resistance-breaking mechanism. We present two de novo draft assembles for N. ribisnigri genomes, corresponding to both avirulent (Nr-locus susceptible) and virulent (Nr-locus resistant) biotypes. Changes in gene expression of the two N. ribisnigri biotypes were investigated using transcriptomic analyses of RNA-sequencing (RNA-seq) data to understand the potential mechanisms of resistance to the Nr-locus in lettuce. The draft genome assemblies were 94.2% and 91.4% complete for the avirulent and virulent biotypes, respectively. Out of the 18,872 differentially expressed genes, a single gene/locus was identified in N. ribisnigri that was shared between two resistant-breaking biotypes. This locus was further explored and validated in Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) experiments and has predicted localisations in both the cytoplasm and nucleus. This is the first study to provide evidence that a single gene/locus is likely responsible for the ability of N. ribisnigri to overcome the Nr-locus resistance in the lettuce host.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Long-Read Sequencing Genome Assembly of Ceratocystis fimbriata Enables Development of Molecular Diagnostics for Sweetpotato Black Rot

Ceratocystis fimbriata is a destructive fungal pathogen of sweetpotato (Ipomoea batatas) that leads to losses at all stages of sweetpotato production. Accurate detection of C. fimbriata would allow for more efficient deployment of management tactics in sweetpotato production. To develop a diagnostic assay, a hybrid genome assembly of C. fimbriata isolate AS236 was generated. The resulting 31.7-MB assembly was near-chromosome level, with 18 contigs, 6,481 predicted genes, and a BUSCO completion score of 98.4% when compared with the fungus-specific lineage database. Additional Illumina DNA reads from C. manginecans, C. platani, and a second C. fimbriata isolate (C1421) were then mapped to the assembled genome using BOWTIE2 and counted using HTSeq, which identified 148 genes present only within C. fimbriata as molecular diagnostic candidates; 6 single-copy and 35 highly multi-copy (>40 BLAST hits), as determined through a self-BLAST-P alignment. Primers for PCR were designed in the 200-bp flanking region of the first exon for each candidate, and the candidates were validated against a diverse DNA panel containing Ceratocystis species, sweetpotato pathogens, and plants. After validation, two diagnostic candidates amplified only C. fimbriata DNA and were considered to be highly specific to the species. These genetic markers will serve as valuable diagnostic tools with multiple applications including the detection of C. fimbriata in seed, soil, and wash water in sweetpotato production.

Assembly and analysis of the genome of Notholithocarpus densiflorus

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

De Novo Genome Assembly for the Coppery Titi Monkey (Plecturocebus cupreus): An Emerging Nonhuman Primate Model for Behavioral Research

The coppery titi monkey (Plecturocebus cupreus) is an emerging nonhuman primate model system for behavioral and neurobiological research. At the same time, the almost entire absence of genomic resources for the species has hampered insights into the genetic underpinnings of the phenotypic traits of interest. To facilitate future genotype-to-phenotype studies, we here present a high-quality, fully annotated de novo genome assembly for the species with chromosome-length scaffolds spanning the autosomes and chromosome X (scaffold N50 = 130.8 Mb), constructed using data obtained from several orthologous short- and long-read sequencing and scaffolding techniques. With a base-level accuracy of ∼99.99% in chromosome-length scaffolds as well as benchmarking universal single-copy ortholog and k-mer completeness scores of >99.0% and 95.1% at the genome level, this assembly represents one of the most complete Pitheciidae genomes to date, making it an invaluable resource for comparative evolutionary genomics research to improve our understanding of lineage-specific changes underlying adaptive traits as well as deleterious mutations associated with disease.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Genomes of multicellular algal sisters to land plants illuminate signaling network evolution

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Nitrogen-fixing organelle in a marine alga

Symbiotic interactions were key to the evolution of chloroplast and mitochondria organelles, which mediate carbon and energy metabolism in eukaryotes. Biological nitrogen fixation, the reduction of abundant atmospheric nitrogen gas (N2) to biologically available ammonia, is a key metabolic process performed exclusively by prokaryotes. Candidatus Atelocyanobacterium thalassa, or UCYN-A, is a metabolically streamlined N2-fixing cyanobacterium previously reported to be an endosymbiont of a marine unicellular alga. Here we show that UCYN-A has been tightly integrated into algal cell architecture and organellar division and that it imports proteins encoded by the algal genome. These are characteristics of organelles and show that UCYN-A has evolved beyond endosymbiosis and functions as an early evolutionary stage N2-fixing organelle, or "nitroplast."

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

Earl Grey: A Fully Automated User-Friendly Transposable Element Annotation and Analysis Pipeline

Transposable elements (TEs) are major components of eukaryotic genomes and are implicated in a range of evolutionary processes. Yet, TE annotation and characterization remain challenging, particularly for nonspecialists, since existing pipelines are typically complicated to install, run, and extract data from. Current methods of automated TE annotation are also subject to issues that reduce overall quality, particularly (i) fragmented and overlapping TE annotations, leading to erroneous estimates of TE count and coverage, and (ii) repeat models represented by short sections of total TE length, with poor capture of 5' and 3' ends. To address these issues, we present Earl Grey, a fully automated TE annotation pipeline designed for user-friendly curation and annotation of TEs in eukaryotic genome assemblies. Using nine simulated genomes and an annotation of Drosophila melanogaster, we show that Earl Grey outperforms current widely used TE annotation methodologies in ameliorating the issues mentioned above while scoring highly in benchmarking for TE annotation and classification and being robust across genomic contexts. Earl Grey provides a comprehensive and fully automated TE annotation toolkit that provides researchers with paper-ready summary figures and outputs in standard formats compatible with other bioinformatics tools. Earl Grey has a modular format, with great scope for the inclusion of additional modules focused on further quality control and tailored analyses in future releases.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

Chromosomal-Level Reference Genome for the Chinese Endemic Pygmy Grasshopper, Zhengitettix transpicula, Sheds Light on Tetrigidae Evolution and Advancing Conservation Efforts

The pygmy grasshopper, Zhengitettix transpicula, is a Chinese endemic species with an exceedingly limited distribution and fragile population structure, rendering it vulnerable to extinction. We present a high-continuity, chromosome-scale reference genome assembly to elucidate this species' distinctive biology and inform conservation. Employing an integrated sequencing approach, we achieved a 970.40 Mb assembly with 96.32% coverage across seven pseudo-chromosomes and impressive continuity (N50 > 220 Mb). Genome annotation achieves identification with 99.2% BUSCO completeness, supporting quality. Comparative analyses with 14 genomes from Orthoptera-facilitated phylogenomics and revealed 549 significantly expanded gene families in Z. transpicula associated with metabolism, stress response, and development. However, genomic analysis exposed remarkably low heterozygosity (0.02%), implying a severe genetic bottleneck from small, fragmented populations, characteristic of species vulnerable to extinction from environmental disruptions. Elucidating the genetic basis of population dynamics and specialization provides an imperative guideline for habitat conservation and restoration of this rare organism. Moreover, divergent evolution analysis of the CYP305m2 gene regulating locust aggregation highlighted potential structural and hence functional variations between Acrididae and Tetrigidae. Our chromosomal genomic characterization of Z. transpicula advances Orthopteran resources, establishing a framework for evolutionary developmental explorations and applied conservation genomics, reversing the trajectory of this unique grasshopper lineage towards oblivion.

Human pangenome analysis of sequences missing from the reference genome reveals their widespread evolutionary, phenotypic, and functional roles

Nonreference sequences (NRSs) are DNA sequences present in global populations but absent in the current human reference genome. However, the extent and functional significance of NRSs in the human genomes and populations remains unclear. Here, we de novo assembled 539 genomes from five genetically divergent human populations using long-read sequencing technology, resulting in the identification of 5.1 million NRSs. These were merged into 45284 unique NRSs, with 29.7% being novel discoveries. Among these NRSs, 38.7% were common across the five populations, and 35.6% were population specific. The use of a graph-based pangenome approach allowed for the detection of 565 transcript expression quantitative trait loci on NRSs, with 426 of these being novel findings. Moreover, 26 NRS candidates displayed evidence of adaptive selection within human populations. Genes situated in close proximity to or intersecting with these candidates may be associated with metabolism and type 2 diabetes. Genome-wide association studies revealed 14 NRSs to be significantly associated with eight phenotypes. Additionally, 154 NRSs were found to be in strong linkage disequilibrium with 258 phenotype-associated SNPs in the GWAS catalogue. Our work expands the understanding of human NRSs and provides novel insights into their functions, facilitating evolutionary and biomedical researches.

First whole-genome sequence and assembly of the Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps), a critically endangered species, using Oxford Nanopore Technologies

The Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps) is currently considered one of the most endangered primates in the world and is classified as critically endangered [International union for conservation of nature (IUCN)]. It faces multiple threats, the most significant one being habitat loss due to deforestation in western Ecuador. Genomic tools are keys for the management of endangered species, but this requires a reference genome, which until now was unavailable for A. f. fusciceps. The present study reports the first whole-genome sequence and assembly of A. f. fusciceps generated using Oxford Nanopore long reads. DNA was extracted from a subadult male, and libraries were prepared for sequencing following the Ligation Sequencing Kit SQK-LSK112 workflow. Sequencing was performed using a MinION Mk1C sequencer. The sequencing reads were processed to generate a genome assembly. Two different assemblers were used to obtain draft genomes using raw reads, of which the Flye assembly was found to be superior. The final assembly has a total length of 2.63 Gb and contains 3,861 contigs, with an N50 of 7,560,531 bp. The assembly was analyzed for annotation completeness based on primate ortholog prediction using a high-resolution database, and was found to be 84.3% complete, with a low number of duplicated genes indicating a precise assembly. The annotation of the assembly predicted 31,417 protein-coding genes, comparable with other mammal assemblies. A reference genome for this critically endangered species will allow researchers to gain insight into the genetics of its populations and thus aid conservation and management efforts of this vulnerable species.

Whole-genome sequence and annotation of Penstemon davidsonii

Penstemon is the most speciose flowering plant genus endemic to North America. Penstemon species' diverse morphology and adaptation to various environments have made them a valuable model system for studying evolution. Here, we report the first full reference genome assembly and annotation for Penstemon davidsonii. Using PacBio long-read sequencing and Hi-C scaffolding technology, we constructed a de novo reference genome of 437,568,744 bases, with a contig N50 of 40 Mb and L50 of 5. The annotation includes 18,199 gene models, and both the genome and transcriptome assembly contain over 95% complete eudicot BUSCOs. This genome assembly will serve as a valuable reference for studying the evolutionary history and genetic diversity of the Penstemon genus.

High-Quality Genome Assemblies of 4 Members of the Podospora anserina Species Complex

The filamentous fungus Podospora anserina is a model organism used extensively in the study of molecular biology, senescence, prion biology, meiotic drive, mating-type chromosome evolution, and plant biomass degradation. It has recently been established that P. anserina is a member of a complex of 7 closely related species. In addition to P. anserina, high-quality genomic resources are available for 2 of these taxa. Here, we provide chromosome-level annotated assemblies of the 4 remaining species of the complex, as well as a comprehensive data set of annotated assemblies from a total of 28 Podospora genomes. We find that all 7 species have genomes of around 35 Mb arranged in 7 chromosomes that are mostly collinear and less than 2% divergent from each other at genic regions. We further attempt to resolve their phylogenetic relationships, finding significant levels of phylogenetic conflict as expected from a rapid and recent diversification.

Genetic modification of Candida maltosa, a non-pathogenic CTG species, reveals EFG1 function

Candida maltosa is closely related to important pathogenic Candida species, especially C. tropicalis and C. albicans, but it has been rarely isolated from humans. For this reason, through comparative studies, it could be a powerful model to understand the genetic underpinnings of the pathogenicity of Candida species. Here, we generated a cohesive assembly of the C. maltosa genome and developed genetic engineering tools that will facilitate studying this species at a molecular level. We used a combination of short and long-read sequencing to build a polished genomic draft composed of 14 Mbp, 45 contigs and close to 5700 genes. This assembly represents a substantial improvement from the currently available sequences that are composed of thousands of contigs. Genomic comparison with C. albicans and C. tropicalis revealed a substantial reduction in the total number of genes in C. maltosa. However, gene loss seems not to be associated to the avirulence of this species given that most genes that have been previously associated with pathogenicity were also present in C. maltosa. To be able to edit the genome of C. maltosa we generated a set of triple auxotrophic strains so that gene deletions can be performed similarly to what has been routinely done in pathogenic Candida species. As a proof of concept, we generated gene knockouts of EFG1, a gene that encodes a transcription factor that is essential for filamentation and biofilm formation in C. albicans and C. tropicalis. Characterization of these mutants showed that Efg1 also plays a role in biofilm formation and filamentous growth in C. maltosa, but it seems to be a repressor of filamentation in this species. The genome assembly and auxotrophic mutants developed here are a key step forward to start using C. maltosa for comparative and evolutionary studies at a molecular level.

Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Assembly and annotation of 2 high-quality columbid reference genomes from sequencing of a Columba livia × Columba guinea F1 hybrid

Pigeons and doves (family Columbidae) are one of the most diverse extant avian lineages, and many species have served as key models for evolutionary genomics, developmental biology, physiology, and behavioral studies. Building genomic resources for columbids is essential to further many of these studies. Here, we present high-quality genome assemblies and annotations for 2 columbid species, Columba livia and Columba guinea. We simultaneously assembled C. livia and C. guinea genomes from long-read sequencing of a single F1 hybrid individual. The new C. livia genome assembly (Cliv_3) shows improved completeness and contiguity relative to Cliv_2.1, with an annotation incorporating long-read IsoSeq data for more accurate gene models. Intensive selective breeding of C. livia has given rise to hundreds of breeds with diverse morphological and behavioral characteristics, and Cliv_3 offers improved tools for mapping the genomic architecture of interesting traits. The C. guinea genome assembly is the first for this species and is a new resource for avian comparative genomics. Together, these assemblies and annotations provide improved resources for functional studies of columbids and avian comparative genomics in general.

The chromosome-scale genome and the genetic resistance machinery against insect herbivores of the Mexican toloache, Datura stramonium

Plant resistance refers to the heritable ability of plants to reduce damage caused by natural enemies, such as herbivores and pathogens, either through constitutive or induced traits like chemical compounds or trichomes. However, the genetic architecture-the number and genome location of genes that affect plant defense and the magnitude of their effects-of plant resistance to arthropod herbivores in natural populations remains poorly understood. In this study, we aimed to unveil the genetic architecture of plant resistance to insect herbivores in the annual herb Datura stramonium (Solanaceae) through quantitative trait loci mapping. We achieved this by assembling the species' genome and constructing a linkage map using an F2 progeny transplanted into natural habitats. Furthermore, we conducted differential gene expression analysis between undamaged and damaged plants caused by the primary folivore, Lema daturaphila larvae. Our genome assembly resulted in 6,109 scaffolds distributed across 12 haploid chromosomes. A single quantitative trait loci region on chromosome 3 was associated with plant resistance, spanning 0 to 5.17 cM. The explained variance by the quantitative trait loci was 8.44%. Our findings imply that the resistance mechanisms of D. stramonium are shaped by the complex interplay of multiple genes with minor effects. Protein-protein interaction networks involving genes within the quantitative trait loci region and overexpressed genes uncovered the key role of receptor-like cytoplasmic kinases in signaling and regulating tropane alkaloids and terpenoids, which serve as powerful chemical defenses against D. stramonium herbivores. The data generated in our study constitute important resources for delving into the evolution and ecology of secondary compounds mediating plant-insect interactions.

Genome sequencing and functional analysis of a multipurpose medicinal herb Tinospora cordifolia (Giloy)

Tinospora cordifolia (Willd.) Hook.f. & Thomson, also known as Giloy, is among the most important medicinal plants that have numerous therapeutic applications in human health due to the production of a diverse array of secondary metabolites. To gain genomic insights into the medicinal properties of T. cordifolia, the genome sequencing was carried out using 10× Genomics linked read and Nanopore long-read technologies. The draft genome assembly of T. cordifolia was comprised of 1.01 Gbp, which is the genome sequenced from the plant family Menispermaceae. We also performed the genome size estimation for T. cordifolia, which was found to be 1.13 Gbp. The deep sequencing of transcriptome from the leaf tissue was also performed. The genome and transcriptome assemblies were used to construct the gene set, resulting in 17,245 coding gene sequences. Further, the phylogenetic position of T. cordifolia was also positioned as basal eudicot by constructing a genome-wide phylogenetic tree using multiple species. Further, a comprehensive comparative evolutionary analysis of gene families contraction/expansion and multiple signatures of adaptive evolution was performed. The genes involved in benzyl iso-quinoline alkaloid, terpenoid, lignin and flavonoid biosynthesis pathways were found with signatures of adaptive evolution. These evolutionary adaptations in genes provide genomic insights into the presence of diverse medicinal properties of this plant. The genes involved in the common symbiosis signalling pathway associated with endosymbiosis (Arbuscular Mycorrhiza) were found to be adaptively evolved. The genes involved in adventitious root formation, peroxisome biogenesis, biosynthesis of phytohormones, and tolerance against abiotic and biotic stresses were also found to be adaptively evolved in T. cordifolia.

Genome assembly and microsatellite marker development using Illumina and PacBio sequencing in Persicaria maackiana (Polygonaceae) from Korea

BACKGROUND

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level.

OBJECTIVE

Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina.

METHODS

Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers.

RESULTS

The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754.

CONCLUSION

Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

The genome assembly and annotation of the white-lipped tree pit viper Trimeresurus albolabris

Trimeresurus albolabris, also known as the white-lipped pit viper or white-lipped tree viper, is a highly venomous snake distributed across Southeast Asia and the cause of many snakebite cases. In this study, we report the first whole genome assembly of T. albolabris obtained with next-generation sequencing from a specimen collected in Mengzi, Yunnan, China. After genome sequencing and assembly, the genome of this male T. albolabris individual was 1.51 Gb in length and included 38.42% repeat-element content. Using this genome, 21,695 genes were identified, and 99.17% of genes could be annotated using gene functional databases. Our genome assembly and annotation process was validated using a phylogenetic tree, which included six species and focused on single-copy genes of nuclear genomes. This research will contribute to future studies on Trimeresurus biology and the genetic basis of snake venom.

Integration of Yeast Episomal/Integrative Plasmid Causes Genotypic and Phenotypic Diversity and Improved Sesquiterpene Production in Metabolically Engineered Saccharomyces cerevisiae

The variability in phenotypic outcomes among biological replicates in engineered microbial factories presents a captivating mystery. Establishing the association between phenotypic variability and genetic drivers is important to solve this intricate puzzle. We applied a previously developed auxin-inducible depletion of hexokinase 2 as a metabolic engineering strategy for improved nerolidol production in Saccharomyces cerevisiae, and biological replicates exhibit a dichotomy in nerolidol production of either 3.5 or 2.5 g L-1 nerolidol. Harnessing Oxford Nanopore's long-read genomic sequencing, we reveal a potential genetic cause─the chromosome integration of a 2μ sequence-based yeast episomal plasmid, encoding the expression cassettes for nerolidol synthetic enzymes. This finding was reinforced through chromosome integration revalidation, engineering nerolidol and valencene production strains, and generating a diverse pool of yeast clones, each uniquely fingerprinted by gene copy numbers, plasmid integrations, other genomic rearrangements, protein expression levels, growth rate, and target product productivities. Τhe best clone in two strains produced 3.5 g L-1 nerolidol and ∼0.96 g L-1 valencene. Comparable genotypic and phenotypic variations were also generated through the integration of a yeast integrative plasmid lacking 2μ sequences. Our work shows that multiple factors, including plasmid integration status, subchromosomal location, gene copy number, sesquiterpene synthase expression level, and genome rearrangement, together play a complicated determinant role on the productivities of sesquiterpene product. Integration of yeast episomal/integrative plasmids may be used as a versatile method for increasing the diversity and optimizing the efficiency of yeast cell factories, thereby uncovering metabolic control mechanisms.

A standard workflow for community-driven manual curation of Strongyloides genome annotations

Advances in the functional genomics and bioinformatics toolkits for Strongyloides species have positioned these species as genetically tractable model systems for gastrointestinal parasitic nematodes. As community interest in mechanistic studies of Strongyloides species continues to grow, publicly accessible reference genomes and associated genome annotations are critical resources for researchers. Genome annotations for multiple Strongyloides species are broadly available via the WormBase and WormBase ParaSite online repositories. However, a recent phylogenetic analysis of the receptor-type guanylate cyclase (rGC) gene family in two Strongyloides species highlights the potential for errors in a large percentage of current Strongyloides gene models. Here, we present three examples of gene annotation updates within the Strongyloides rGC gene family; each example illustrates a type of error that may occur frequently within the annotation data for Strongyloides genomes. We also extend our analysis to 405 previously curated Strongyloides genes to confirm that gene model errors are found at high rates across gene families. Finally, we introduce a standard manual curation workflow for assessing gene annotation quality and generating corrections, and we discuss how it may be used to facilitate community-driven curation of parasitic nematode biodata. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Haplotype-resolved chromosome-scale genomes of the Asian and African Savannah Elephants

The Proboscidea, which includes modern elephants, were once the largest terrestrial animals among extant species. They suffered mass extinction during the Ice Age. As a unique branch on the evolutionary tree, the Proboscidea are of great significance for the study of living animals. In this study, we generate chromosome-scale and haplotype-resolved genome assemblies for two extant Proboscidea species (Asian Elephant, Elephas maximus and African Savannah Elephant, Loxodonta africana) using Pacbio, Hi-C, and DNBSEQ technologies. The assembled genome sizes of the Asian and African Savannah Elephant are 3.38 Gb and 3.31 Gb, with scaffold N50 values of 130 Mb and 122 Mb, respectively. Using Hi-C technology ~97% of the scaffolds are anchored to 29 pseudochromosomes. Additionally, we identify ~9 Mb Y-linked sequences for each species. The high-quality genome assemblies in this study provide a valuable resource for future research on ecology, evolution, biology and conservation of Proboscidea species.

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.

Developmental transcriptomes predict adult social behaviours in the socially flexible sweat bee, Lasioglossum baleicum

Natural variation can provide important insights into the genetic and environmental factors that shape social behaviour and its evolution. The sweat bee, Lasioglossum baleicum, is a socially flexible bee capable of producing both solitary and eusocial nests. We demonstrate that within a single nesting aggregation, soil temperatures are a strong predictor of the social structure of nests. Sites with warmer temperatures in the spring have a higher frequency of social nests than cooler sites, perhaps because warmer temperatures provide a longer reproductive window for those nests. To identify the molecular correlates of this behavioural variation, we generated a de novo genome assembly for L. baleicum, and we used transcriptomic profiling to compare adults and developing offspring from eusocial and solitary nests. We find that adult, reproductive females have similar expression profiles regardless of social structure in the nest, but that there are strong differences between reproductive females and workers from social nests. We also find substantial differences in the transcriptomic profiles of stage-matched pupae from warmer, social-biased sites compared to cooler, solitary-biased sites. These transcriptional differences are strongly predictive of adult reproductive state, suggesting that the developmental environment may set the stage for adult behaviours in L. baleicum. Together, our results help to characterize the molecular mechanisms shaping variation in social behaviour and highlight a potential role of environmental tuning during development as a factor shaping adult behaviour and physiology in this socially flexible bee.

Genome assembly and annotation of the king ratsnake, Elaphe carinata

The king ratsnake (Elaphe carinata) of the genus Elaphe is a common large, non-venomous snake widely distributed in Southeast and East Asia. It is an economically important farmed species. As a non-venomous snake, the king ratsnake predates venomous snakes, such as cobras and pit vipers. However, the immune and digestive mechanisms of the king ratsnake remain unclear. Despite their economic and research importance, we lack genomic resources that would benefit toxicology, phylogeography, and immunogenetics studies. Here, we used single-tube long fragment read sequencing to generate the first contiguous genome of a king ratsnake from Huangshan City, Anhui province, China. The genome size is 1.56 GB with a scaffold N50 of 6.53M. The total length of the genome is approximately 621 Mb, while the repeat content is 42.26%. Additionally, we predicted 22,339 protein-coding genes, including 22,065 with functional annotations. Our genome is a potentially useful addition to those available for snakes.

Evolution of ACCase-inhibitor resistance in Chloris virgata is conferred by a Trp2027Cys mutation in the herbicide target site

BACKGROUND

Chloris virgata is a troublesome weed in tropical regions. With the evolution of glyphosate resistance in key grass species, acetyl CoA carboxylase (ACCase) inhibitors have become a commonly used tool in soybean production areas in Brazil. We assessed if suspected resistant populations exhibited cross resistance to the different classes of ACCase inhibitors and investigated the resistance mechanisms in C. virgata.

RESULTS

Dose-response experiments revealed resistance to haloxyfop-methyl and pinoxaden, with 432- and 3-fold resistance, respectively, compared to susceptible populations. Due to the lack of genetic resources for C. virgata, we sequenced, assembled, and annotated the genome using short-read Illumina technology. The k-mer analysis estimated a genome size of approximately 336 Mbp, with BUSCO completeness of 97%, and over 36 000 gene models were annotated. We examined if ACCase copy number variation and increased gene expression were involved in the resistance phenotype and found no difference when compared to a susceptible population. A mutation was detected in ACCase that encodes for amino acid position 2027, resulting in a tryptophan-to-cysteine (Trp2027Cys) substitution. We found the resistant population absorbed 11.4% less herbicide and retained 21% more herbicide on the treated leaf compared to the susceptible population. We developed a genotyping assay targeting the resistance-endowing Trp2027Cys substitution for quick resistance diagnosis.

CONCLUSION

A Trp2027Cys amino acid substitution in ACCase confers resistance to haloxyfop and pinoxaden in C. virgata. We provide important insights into the evolutionary history of C. virgata and a draft genome as a useful resource to further our understanding of the biology in the genus Chloris. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Chromosome-Level Genome Assembly of the Blacktail Brush Lizard, Urosaurus nigricaudus, Reveals Dosage Compensation in an Endemic Lizard

Urosaurus nigricaudus is a phrynosomatid lizard endemic to the Baja California Peninsula in Mexico. This work presents a chromosome-level genome assembly and annotation from a male individual. We used PacBio long reads and HiRise scaffolding to generate a high-quality genomic assembly of 1.87 Gb distributed in 327 scaffolds, with an N50 of 279 Mb and an L50 of 3. Approximately 98.4% of the genome is contained in 14 scaffolds, with 6 large scaffolds (334-127 Mb) representing macrochromosomes and 8 small scaffolds (63-22 Mb) representing microchromosomes. Using standard gene modeling and transcriptomic data, we predicted 17,902 protein-coding genes on the genome. The repeat content is characterized by a large proportion of long interspersed nuclear elements that are relatively old. Synteny analysis revealed some microchromosomes with high repeat content are more prone to rearrangements but that both macro- and microchromosomes are well conserved across reptiles. We identified scaffold 14 as the X chromosome. This microchromosome presents perfect dosage compensation where the single X of males has the same expression levels as two X chromosomes in females. Finally, we estimated the effective population size for U. nigricaudus was extremely low, which may reflect a reduction in polymorphism related to it becoming a peninsular endemic.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Multi-omics analysis reveals the molecular response to heat stress in a "red tide" dinoflagellate

BACKGROUND

"Red tides" are harmful algal blooms caused by dinoflagellate microalgae that accumulate toxins lethal to other organisms, including humans via consumption of contaminated seafood. These algal blooms are driven by a combination of environmental factors including nutrient enrichment, particularly in warm waters, and are increasingly frequent. The molecular, regulatory, and evolutionary mechanisms that underlie the heat stress response in these harmful bloom-forming algal species remain little understood, due in part to the limited genomic resources from dinoflagellates, complicated by the large sizes of genomes, exhibiting features atypical of eukaryotes.

RESULTS

We present the de novo assembled genome (~ 4.75 Gbp with 85,849 protein-coding genes), transcriptome, proteome, and metabolome from Prorocentrum cordatum, a globally abundant, bloom-forming dinoflagellate. Using axenic algal cultures, we study the molecular mechanisms that underpin the algal response to heat stress, which is relevant to current ocean warming trends. We present the first evidence of a complementary interplay between RNA editing and exon usage that regulates the expression and functional diversity of biomolecules, reflected by reduction in photosynthesis, central metabolism, and protein synthesis. These results reveal genomic signatures and post-transcriptional regulation for the first time in a pelagic dinoflagellate.

CONCLUSIONS

Our multi-omics analyses uncover the molecular response to heat stress in an important bloom-forming algal species, which is driven by complex gene structures in a large, high-G+C genome, combined with multi-level transcriptional regulation. The dynamics and interplay of molecular regulatory mechanisms may explain in part how dinoflagellates diversified to become some of the most ecologically successful organisms on Earth.

Genome sequencing of Syzygium cumini (jamun) reveals adaptive evolution in secondary metabolism pathways associated with its medicinal properties

Syzygium cumini, also known as jambolan or jamun, is an evergreen tree widely known for its medicinal properties, fruits, and ornamental value. To understand the genomic and evolutionary basis of its medicinal properties, we sequenced S. cumini genome for the first time from the world's largest tree genus Syzygium using Oxford Nanopore and 10x Genomics sequencing technologies. We also sequenced and assembled the transcriptome of S. cumini in this study. The tetraploid and highly heterozygous draft genome of S. cumini had a total size of 709.9 Mbp with 61,195 coding genes. The phylogenetic position of S. cumini was established using a comprehensive genome-wide analysis including species from 18 Eudicot plant orders. The existence of neopolyploidy in S. cumini was evident from the higher number of coding genes and expanded gene families resulting from gene duplication events compared to the other two sequenced species from this genus. Comparative evolutionary analyses showed the adaptive evolution of genes involved in the phenylpropanoid-flavonoid (PF) biosynthesis pathway and other secondary metabolites biosynthesis such as terpenoid and alkaloid in S. cumini, along with genes involved in stress tolerance mechanisms, which was also supported by leaf transcriptome data generated in this study. The adaptive evolution of secondary metabolism pathways is associated with the wide range of pharmacological properties, specifically the anti-diabetic property, of this species conferred by the bioactive compounds that act as nutraceutical agents in modern medicine.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.