STRIDE: Species Tree Root Inference from Gene Duplication Events

Comparative transcriptomics reveals a highly polymorphic Xanthomonas HrpG virulence regulon

BACKGROUND

Bacteria of the genus Xanthomonas cause economically significant diseases in various crops. Their virulence is dependent on the translocation of type III effectors (T3Es) into plant cells by the type III secretion system (T3SS), a process regulated by the master response regulator HrpG. Although HrpG has been studied for over two decades, its regulon across diverse Xanthomonas species, particularly beyond type III secretion, remains understudied.

RESULTS

In this study, we conducted transcriptome sequencing to explore the HrpG regulons of 17 Xanthomonas strains, encompassing six species and nine pathovars, each exhibiting distinct host and tissue specificities. We employed constitutive expression of plasmid-borne hrpG*, which encodes a constitutively active form of HrpG, to induce the regulon. Our findings reveal substantial inter- and intra-specific diversity in the HrpG* regulons across the strains. Besides 21 genes directly involved in the biosynthesis of the T3SS, the core HrpG* regulon is limited to only five additional genes encoding the transcriptional activator HrpX, the two T3E proteins XopR and XopL, a major facility superfamily (MFS) transporter, and the phosphatase PhoC. Interestingly, genes involved in chemotaxis and genes encoding enzymes with carbohydrate-active and proteolytic activities are variably regulated by HrpG*.

CONCLUSIONS

The diversity in the HrpG* regulon suggests that HrpG-dependent virulence in Xanthomonas might be achieved through several distinct strain-specific strategies, potentially reflecting adaptation to diverse ecological niches. These findings enhance our understanding of the complex role of HrpG in regulating various virulence and adaptive pathways, extending beyond T3Es and the T3SS.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Genome divergence and reproductive incompatibility among populations of Ganaspis near brasiliensis

During the last decade, the spotted wing drosophila, Drosophila suzukii, has spread from eastern Asia to the Americas, Europe, and Africa. This fly attacks many species of cultivated and wild fruits with soft, thin skins, where its serrated ovipositor allows it to lay eggs in undamaged fruit. Parasitoids from the native range of D. suzukii may provide sustainable management of this polyphagous pest. Among these parasitoids, host-specificity testing has revealed a lineage of Ganaspis near brasiliensis, referred to in this paper as G1, that appears to be a cryptic species more host-specific to D. suzukii than other parasitoids. Differentiation among cryptic species is critical for introduction and subsequent evaluation of their impact on D. suzukii. Here, we present results on divergence in genomic sequences and architecture and reproductive isolation between lineages of Ganaspis near brasiliensis that appear to be cryptic species. We studied five populations, two from China, two from Japan, and one from Canada, identified as the G1 vs G3 lineages based on differences in cytochrome oxidase l sequences. We assembled and annotated the genomes of these populations and analyzed divergences in sequence and genome architecture between them. We also report results from crosses to test reproductive compatibility between the G3 lineage from China and the G1 lineage from Japan. The combined results on sequence divergence, differences in genome architectures, ortholog divergence, reproductive incompatibility, differences in host ranges and microhabitat preferences, and differences in morphology show that these lineages are different species. Thus, the decision to evaluate the lineages separately and only import and introduce the more host-specific lineage to North America and Europe was appropriate.

Oomycete Metabolism Is Highly Dynamic and Reflects Lifestyle Adaptations

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyze the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as animal and plant pathogens, with the latter causing devastating diseases with significant economic and/or ecological impacts. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defense compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolisms. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation. Numerous oomycete species are devastating plant pathogens that cause major damage in crops and natural ecosystems. Their interactions with hosts are shaped by strong selection, but how selection affects adaptation of the primary metabolism to a pathogenic lifestyle is not yet well established. By pan-genome and metabolic network analyses of distantly related oomycete pathogens and their nonpathogenic relatives, we reveal considerable lifestyle- and lineage-specific adaptations. This study contributes to a better understanding of metabolic adaptations in pathogenic oomycetes in relation to lifestyle, host, and environment, and the findings will help in pinpointing potential targets for disease control. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Giants' cooperation: a draft genome of the giant ciliate Muniziella cunhai suggests its ecological role in the capybara's digestive metabolism

Several hundred ciliate species live in animals' guts as a part of their microbiome. Among them, Muniziella cunhai (Trichostomatia, Pycnotrichidae), the largest described ciliate, is found exclusively associated with Hydrochoerus hydrochaeris (capybara), the largest known rodent reaching up to 90 kg. Here, we present the sequence, structural and functional annotation of this giant microeukaryote macronuclear genome and discuss its phylogenetic placement. The 85 Mb genome is highly AT rich (GC content 25.71 %) and encodes a total of 11 397 protein-coding genes, of which 2793 could have their functions predicted with automated functional assignments. Functional annotation showed that M. cunhai can digest recalcitrant structural carbohydrates, non-structural carbohydrates, and microbial cell walls, suggesting a role in diet metabolization and in microbial population control in the capybara's intestine. Moreover, the phylogenetic placement of M. cunhai provides insights on the origins of gigantism in the subclass Trichostomatia.

The genome assembly and annotation of the cricket Gryllus longicercus

The order Orthoptera includes insects such as grasshoppers, katydids, and crickets, among which there are important species for ecosystem stability and pollination, as well as research organisms in different fields such as neurobiology, ecology, and evolution. Crickets, with more than 2,400 described species, are emerging as novel model research organisms, for their diversity, worldwide distribution, regeneration capacity, and their characteristic acoustic communication. Here we report the assembly and annotation of the first New World cricket, that of Gryllus longicercus Weissman & Gray 2019. The genome assembly, generated by combining 44.54 Gb of long reads from PacBio and 120.44 Gb of short Illumina reads, has a length of 1.85 Gb. The genome annotation yielded 19,715 transcripts from 14,789 gene models.

Comparative genomic analysis of strain Priestia megaterium B1 reveals conserved potential for adaptation to endophytism and plant growth promotion

In our study, we aimed to explore the genomic and phenotypic traits of Priestia megaterium strain B1, which was isolated from root material of healthy apple plants, to adapt to the endophytic lifestyle and promote plant growth. We identified putative genes encoding proteins involved in chemotaxis, flagella biosynthesis, biofilm formation, secretory systems, detoxification, transporters, and transcription regulation. Furthermore, B1 exhibited both swarming and swimming motilities, along with biofilm formation. Both genomic and physiological analyses revealed the potential of B1 to promote plant growth through the production of indole-3-acetic acid and siderophores, as well as the solubilization of phosphate and zinc. To deduce potential genomic features associated with endophytism across members of P. megaterium strains, we conducted a comparative genomic analysis involving 27 and 31 genomes of strains recovered from plant and soil habitats, respectively, in addition to our strain B1. Our results indicated a closed pan genome and comparable genome size of strains from both habitats, suggesting a facultative host association and adaptive lifestyle to both habitats. Additionally, we performed a sparse Partial Least Squares Discriminant Analysis to infer the most discriminative functional features of the two habitats based on Pfam annotation. Despite the distinctive clustering of both groups, functional enrichment analysis revealed no significant enrichment of any Pfam domain in both habitats. Furthermore, when assessing genetic elements related to adaptation to endophytism in each individual strain, we observed their widespread presence among strains from both habitats. Moreover, all members displayed potential genetic elements for promoting plant growth.IMPORTANCEBoth genomic and phenotypic analyses yielded valuable insights into the capacity of P. megaterium B1 to adapt to the plant niche and enhance its growth. The comparative genomic analysis revealed that P. megaterium members, whether derived from soil or plant sources, possess the essential genetic machinery for interacting with plants and enhancing their growth. The conservation of these traits across various strains of this species extends its potential application as a bio-stimulant in diverse environments. This significance also applies to strain B1, particularly regarding its application to enhance the growth of plants facing apple replant disease conditions.

De Novo Assembly and Annotation of the Siganus fuscescens (Houttuyn, 1782) Genome: Marking a Pioneering Advance for the Siganidae Family

This study presents the first draft genome of Siganus fuscescens, and thereby establishes the first whole-genome sequence for a species in the Siganidae family. Leveraging both long and short read sequencing technologies, i.e., Oxford Nanopore and Illumina sequencing, we successfully assembled a mitogenome spanning 16.494 Kb and a first haploid genome encompassing 498 Mb. The assembled genome accounted for a 99.6% of the estimated genome size and was organized into 164 contigs with an N50 of 7.2 Mb. This genome assembly showed a GC content of 42.9% and a high Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness score of 99.5% using actinopterygii_odb10 lineage, thereby meeting stringent quality standards. In addition to its structural aspects, our study also examined the functional genomics of this species, including the intricate capacity to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFAs) and secrete venom. Notably, our analyses revealed various repeats elements, which collectively constituted 17.43% of the genome. Moreover, annotation of 28,351 genes uncovered both shared genetic signatures and those that are unique to S. fuscescens. Our assembled genome also displayed a moderate prevalence of gene duplication compared to other fish species, which suggests that this species has a distinctive evolutionary trajectory and potentially unique functional constraints. Taken altogether, this genomic resource establishes a robust foundation for future research on the biology, evolution, and the aquaculture potential of S. fuscescens.

De novo transcriptomes of cave and surface isopod crustaceans: insights from 11 species across three suborders

Isopods are a diverse group of crustaceans, that inhabit various environments, including terrestrial, freshwater, and marine, both on the surface and in the underground. The biological mechanisms underlying their wide range of adaptations to diverse ecological niches remain elusive. In order to unravel the molecular basis of their adaptability, we generated a comprehensive RNAseq dataset comprising 11 isopod species belonging to the three different suborders: freshwater Asellota, marine, brackish and freshwater Sphaeromatidea, and terrestrial Oniscidea, with representatives from families Asellidae, Sphaeromatidae, and Trichoniscidae, respectively. Representatives of each family were collected from both cave and surface environments, representing at least three independent cave colonization events. Three biological replicates were sequenced from each species to ensure data robustness. The 11 high-quality RNAseq datasets will serve as a valuable resource for understanding cave-specific adaptations, comparative and functional genomics, ecological annotation as well as aid in conservation efforts of these non-model organisms. Importantly, transcriptomes of eight featured species have been made publicly accessible for the first time.

Catabolism and interactions of syntrophic propionate- and acetate oxidizing microorganisms under mesophilic, high-ammonia conditions

Microbial inhibition by high ammonia concentrations is a recurring problem that significantly restricts methane formation from intermediate acids, i.e., propionate and acetate, during anaerobic digestion of protein-rich waste material. Studying the syntrophic communities that perform acid conversion is challenging, due to their relatively low abundance within the microbial communities typically found in biogas processes and disruption of their cooperative behavior in pure cultures. To overcome these limitations, this study examined growth parameters and microbial community dynamics of highly enriched mesophilic and ammonia-tolerant syntrophic propionate and acetate-oxidizing communities and analyzed their metabolic activity and cooperative behavior using metagenomic and metatranscriptomic approaches. Cultivation in batch set-up demonstrated biphasic utilization of propionate, wherein acetate accumulated and underwent oxidation before complete degradation of propionate. Three key species for syntrophic acid degradation were inferred from genomic sequence information and gene expression: a syntrophic propionate-oxidizing bacterium (SPOB) "Candidatus Syntrophopropionicum ammoniitolerans", a syntrophic acetate-oxidizing bacterium (SAOB) Syntrophaceticus schinkii and a novel hydrogenotrophic methanogen, for which we propose the provisional name "Candidatus Methanoculleus ammoniitolerans". The results revealed consistent transcriptional profiles of the SAOB and the methanogen both during propionate and acetate oxidation, regardless of the presence of an active propionate oxidizer. Gene expression indicated versatile capabilities of the two syntrophic bacteria, utilizing both molecular hydrogen and formate as an outlet for reducing equivalents formed during acid oxidation, while conserving energy through build-up of sodium/proton motive force. The methanogen used hydrogen and formate as electron sources. Furthermore, results of the present study provided a framework for future research into ammonia tolerance, mobility, aggregate formation and interspecies cooperation.

WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids

Summary

Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests.

Availability and Implementation

The tool is available at: https://fishevo.xmu.edu.cn/.

A chromosome-level haplotype-resolved genome assembly of oriental tobacco budworm (Helicoverpa assulta)

Oriental tobacco budworm (Helicoverpa assulta) and cotton bollworm (Helicoverpa armigera) are two closely related species within the genus Helicoverpa. They have similar appearances and consistent damage patterns, often leading to confusion. However, the cotton bollworm is a typical polyphagous insect, while the oriental tobacco budworm belongs to the oligophagous insects. In this study, we used Nanopore, PacBio, and Illumina platforms to sequence the genome of H. assulta and used Hifiasm to create a haplotype-resolved draft genome. The Hi-C technique helped anchor 33 primary contigs to 32 chromosomes, including two sex chromosomes, Z and W. The final primary haploid genome assembly was approximately 415.19 Mb in length. BUSCO analysis revealed a high degree of completeness, with 99.0% gene coverage in this genome assembly. The repeat sequences constituted 38.39% of the genome assembly, and we annotated 17093 protein-coding genes. The high-quality genome assembly of the oriental tobacco budworm serves as a valuable genetic resource that enhances our comprehension of how they select hosts in a complex odour environment. It will also aid in developing an effective control policy.

NmrB (AN9181) expression is activated under oxidative stress conditions acting as a metabolic repressor of Aspergillus nidulans

Aspergilli comprise a diversity of species that have been extensively studied due to their catabolic diversity, biotechnological and ecological value, and pathogenicity. An impressive level of structural and functional conservation has been shown for aspergilli, regardless of many (yet) cryptic genomic elements. We have hypothesized the existence of conserved genes responsive to stress in aspergilli. To test the hypothesis of such conserved stress regulators in aspergilli, a straightforward computational strategy integrating well-established bioinformatic tools was used as the starting point. Specifically, five transcriptome-based datasets on exposure to organic compounds were used, covering three distinct Aspergillus species. Among the identified up-regulated genes, only one gene showed the same response in all conditions, AN9181. This gene encodes a protein containing a phenylcoumaran benzylic ether reductase-like domain and a Nitrogen metabolite repressor regulator domain (NmrA). Deletion of this gene caused significant phenotypic alterations compared to that of the parental strain across diverse conditions. Specifically, the deletion of AN9181 raised the mutant's metabolic activity in different nitrogen sources. The acquired data supports that AN9181 acts by repressing (slowing down) A. nidulans growth when exposed to aromatic compounds in a concentration dependent manner. The same phenotype was observed for amphotericin B. Finally, AN9181 underwent differential upregulation under oxidative stress conditions. Collectively, the data suggest that AN9181, herein assigned as NmrB (Nitrogen Metabolite Repression Regulator B), builds up the genetic machinery of perception of oxidative stress by negatively regulating growth under such conditions.

Genome report: chromosome-scale genome assembly of the West Indian fruit fly Anastrepha obliqua (Diptera: Tephritidae)

The West Indian fruit fly, Anastrepha obliqua, is a major pest of mango in Central and South America and attacks more than 60 species of host fruits. To support current genetic and genomic research on A. obliqua, we sequenced the genome using high-fidelity long-read sequencing. This resulted in a highly contiguous contig assembly with 90% of the genome in 10 contigs. The contig assembly was placed in a chromosomal context using synteny with a closely related species, Anastrepha ludens, as both are members of the Anastrepha fraterculus group. The resulting assembly represents the five autosomes and the X chromosome which represents 95.9% of the genome, and 199 unplaced contigs representing the remaining 4.1%. Orthology analysis across the structural annotation sets of high quality tephritid genomes demonstrates the gene annotations are robust, and identified genes unique to Anastrepha species that may help define their pestiferous nature that can be used as a starting point for comparative genomics. This genome assembly represents the first of this species and will serve as a foundation for future genetic and genomic research in support of its management as an agricultural pest.

Genetic and Functional Diversity Help Explain Pathogenic, Weakly Pathogenic, and Commensal Lifestyles in the Genus Xanthomonas

The genus Xanthomonas has been primarily studied for pathogenic interactions with plants. However, besides host and tissue-specific pathogenic strains, this genus also comprises nonpathogenic strains isolated from a broad range of hosts, sometimes in association with pathogenic strains, and other environments, including rainwater. Based on their incapacity or limited capacity to cause symptoms on the host of isolation, nonpathogenic xanthomonads can be further characterized as commensal and weakly pathogenic. This study aimed to understand the diversity and evolution of nonpathogenic xanthomonads compared to their pathogenic counterparts based on their cooccurrence and phylogenetic relationship and to identify genomic traits that form the basis of a life history framework that groups xanthomonads by ecological strategies. We sequenced genomes of 83 strains spanning the genus phylogeny and identified eight novel species, indicating unexplored diversity. While some nonpathogenic species have experienced a recent loss of a type III secretion system, specifically the hrp2 cluster, we observed an apparent lack of association of the hrp2 cluster with lifestyles of diverse species. We performed association analysis on a large data set of 337 Xanthomonas strains to explain how xanthomonads may have established association with the plants across the continuum of lifestyles from commensals to weak pathogens to pathogens. Presence of distinct transcriptional regulators, distinct nutrient utilization and assimilation genes, transcriptional regulators, and chemotaxis genes may explain lifestyle-specific adaptations of xanthomonads.

Comparative genomics reveals insight into the phylogeny and habitat adaptation of novel Amycolatopsis species, an endophytic actinomycete associated with scab lesions on potato tubers

A novel endophytic actinomycete, strain MEP2-6T, was isolated from scab tissues of potato tubers collected from Mae Fag Mai Sub-district, San Sai District, Chiang Mai Province, Thailand. Strain MEP2-6T is a gram-positive filamentous bacteria characterized by meso-diaminopimelic acid in cell wall peptidoglycan and arabinose, galactose, glucose, and ribose in whole-cell hydrolysates. Diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, and hydroxy-phosphatidylethanolamine were the major phospholipids, of which MK-9(H6) was the predominant menaquinone, whereas iso-C16:0 and iso-C15:0 were the major cellular fatty acids. The genome of the strain was 10,277,369 bp in size with a G + C content of 71.7%. The 16S rRNA gene phylogenetic and core phylogenomic analyses revealed that strain MEP2-6T was closely related to Amycolatopsis lexingtonensis NRRL B-24131T (99.4%), A. pretoriensis DSM 44654T (99.3%), and A. eburnea GLM-1T (98.9%). Notably, strain MEP2-6T displayed 91.7%, 91.8%, and 87% ANIb and 49%, 48.8%, and 35.4% dDDH to A. lexingtonensis DSM 44653T (=NRRL B-24131T), A. eburnea GLM-1T, and A. pretoriensis DSM 44654T, respectively. Based on phenotypic, chemotaxonomic, and genomic data, strain MEP2-6T could be officially assigned to a novel species within the genus Amycolatopsis, for which the name Amycolatopsis solani sp. nov. has been proposed. The type of strain is MEP2-6T (=JCM 36309T = TBRC 17632T = NBRC 116395T). Amycolatopsis solani MEP2-6T was strongly proven to be a non-phytopathogen of potato scab disease because stunting of seedlings and necrotic lesions on potato tuber slices were not observed, and there were no core biosynthetic genes associated with the BGCs of phytotoxin-inducing scab lesions. Furthermore, comparative genomics can provide a better understanding of the genetic mechanisms that enable A. solani MEP2-6T to adapt to the plant endosphere. Importantly, the strain smBGCs accommodated 33 smBGCs encoded for several bioactive compounds, which could be beneficially applied in the fields of agriculture and medicine. Consequently, strain MEP2-6T is a promising candidate as a novel biocontrol agent and antibiotic producer.

Unveiling the evolutionary history of lingonberry (Vaccinium vitis-idaea L.) through genome sequencing and assembly of European and North American subspecies

Lingonberry (Vaccinium vitis-idaea L.) produces tiny red berries that are tart and nutty in flavor. It grows widely in the circumpolar region, including Scandinavia, northern parts of Eurasia, Alaska, and Canada. Although cultivation is currently limited, the plant has a long history of cultural use among indigenous communities. Given its potential as a food source, genomic resources for lingonberry are significantly lacking. To advance genomic knowledge, the genomes for 2 subspecies of lingonberry (V. vitis-idaea ssp. minus and ssp. vitis-idaea var. 'Red Candy') were sequenced and de novo assembled into contig-level assemblies. The assemblies were scaffolded using the bilberry genome (Vaccinium myrtillus) to generate a chromosome-anchored reference genome consisting of 12 chromosomes each with a total length of 548.07 Mb [contig N50 = 1.17 Mb, BUSCO (C%) = 96.5%] for ssp. vitis-idaea and 518.70 Mb [contig N50 = 1.40 Mb, BUSCO (C%) = 96.9%] for ssp. minus. RNA-seq-based gene annotation identified 27,243 and 25,718 genes on the respective assembly, and transposable element detection methods found that 45.82 and 44.58% of the genome were repeats. Phylogenetic analysis confirmed that lingonberry was most closely related to bilberry and was more closely related to blueberries than cranberries. Estimates of past effective population size suggested a continuous decline over the past 1-3 MYA, possibly due to the impacts of repeated glacial cycles during the Pleistocene leading to frequent population fragmentation. The genomic resource created in this study can be used to identify industry-relevant genes (e.g. anthocyanin production), infer phylogeny, and call sequence-level variants (e.g. SNPs) in future research.

High-quality genome assembly and genetic transformation system of Lasiodiplodia theobromae strain LTTK16-3, a fungal pathogen of Chinese hickory

Lasiodiplodia theobromae, as one of the causative agents associated with Chinese hickory trunk cankers, has caused huge economic losses to the Chinese hickory industry. Although the biological characteristics of this pathogen and the occurrence pattern of this disease have been well studied, few studies have addressed the related mechanisms due to the poor molecular and genetic study basis of this fungus. In this study, we sequenced and assembled L. theobromae strain LTTK16-3, isolated from a Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province, China. Phylogenetic analysis and comparative genomics analysis presented crucial cues in the prediction of LTTK16-3, which shared similar regulatory mechanisms of transcription, DNA replication, and DNA damage response with the other four Chinese hickory trunk canker-associated Botryosphaeria strains including, Botryosphaeria dothidea, Botryosphaeria fabicerciana, Botryosphaeria qingyuanensis, and Botryosphaeria corticis. Moreover, it contained 18 strain-specific protein clusters (not conserved in the other L. theobromae strains, AM2As and CITRA15), with potential roles in specific host-pathogen interactions during the Chinese hickory infection. Additionally, an efficient system for L. theobromae protoplast preparation and polyethylene glycol (PEG) -mediated genetic transformation was firstly established as the foundation for its future mechanisms study. Collectively, the high-quality genome data and the efficient transformation system of L. theobromae here set up the possibility of targeted molecular improvements for Chinese hickory canker control.IMPORTANCEFungi with disparate genomic features are physiologically diverse, possessing species-specific survival strategies and environmental adaptation mechanisms. The high-quality genome data and related molecular genetic studies are the basis for revealing the mechanisms behind the physiological traits that are responsible for their environmental fitness. In this study, we sequenced and assembled the LTTK16-3 strain, the genome of Lasiodiplodia theobromae first obtained from a diseased Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province, China. Further phylogenetic analysis and comparative genomics analysis provide crucial cues in the prediction of the proteins with potential roles in specific host-pathogen interactions during the Chinese hickory infection. An efficient PEG-mediated genetic transformation system of L. theobromae was established as the foundation for the future mechanisms exploration. The above genetic information and tools set up valuable clues to study L. theobromae pathogenesis and assist in Chinese hickory canker control.

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.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Chromosome-level genome assembly and demographic history of Euryodendron excelsum in monotypic genus endemic to China

Euryodendron excelsum is in a monotypic genus Euryodendron, endemic to China. It has intermediate morphisms in the Pentaphylacaceae or Theaceae families, which make it distinct. Due to anthropogenic disturbance, E. excelsum is currently found in very restricted and fragmented areas with extremely small populations. Although much research and effort has been applied towards its conservation, its long-term survival mechanisms and evolutionary history remain elusive, especially from a genomic aspect. Therefore, using a combination of long/short whole genome sequencing, RNA sequencing reads, and Hi-C data, we assembled and annotated a high-quality genome for E. excelsum. The genome assembly of E. excelsum comprised 1,059,895,887 bp with 99.66% anchored into 23 pseudo-chromosomes and a 99.0% BUSCO completeness. Comparative genomic analysis revealed the expansion of terpenoid and flavonoid secondary metabolite genes, and displayed a tandem and/or proximal duplication framework of these genes. E. excelsum also displayed genes associated with growth, development, and defence adaptation from whole genome duplication. Demographic analysis indicated that its fluctuations in population size and its recent population decline were related to cold climate changes. The E. excelsum genome assembly provides a highly valuable resource for evolutionary and ecological research in the future, aiding its conservation, management, and restoration.

Highly resolved genome assembly and comparative transcriptome profiling reveal genes related to developmental stages of tapeworm Ligula intestinalis

Ligula intestinalis (Cestoda: Diphyllobothriidae) is an emerging model organism for studies on parasite population biology and host-parasite interactions. However, a well-resolved genome and catalogue of its gene content has not been previously developed. Here, we present the first genome assembly of L. intestinalis, based on Oxford Nanopore Technologies, Illumina and Omni-C sequencing methodologies. We use transcriptome profiling to compare plerocercoid larvae and adult worms and identify differentially expressed genes (DEGs) associated with these life stages. The genome assembly is 775.3 mega (M)bp in size, with scaffold N50 value of 118 Mbp and encodes 27 256 predicted protein-coding sequences. Over 60% of the genome consists of repetitive sequences. Synteny analyses showed that the 10 largest scaffolds representing 75% of the genome display high correspondence to full chromosomes of cyclophyllidean tapeworms. Mapping RNA-seq data to the new reference genome, we identified 3922 differentially expressed genes in adults compared with plerocercoids. Gene ontology analyses revealed over-represented genes involved in reproductive development of the adult stage (e.g. sperm production) and significantly enriched DEGs associated with immune evasion of plerocercoids in their fish host. This study provides the first insights into the molecular biology of L. intestinalis and provides the most highly contiguous assembly to date of a diphyllobothriid tapeworm useful for population and comparative genomic investigations of parasitic flatworms.

Uncovering Evolutionary Adaptations in Common Warthogs through Genomic Analyses

In the Suidae family, warthogs show significant survival adaptability and trait specificity. This study offers a comparative genomic analysis between the warthog and other Suidae species, including the Luchuan pig, Duroc pig, and Red River hog. By integrating the four genomes with sequences from the other four species, we identified 8868 single-copy orthologous genes. Based on 8868 orthologous protein sequences, phylogenetic assessments highlighted divergence timelines and unique evolutionary branches within suid species. Warthogs exist on different evolutionary branches compared to DRCs and LCs, with a divergence time preceding that of DRC and LC. Contraction and expansion analyses of warthog gene families have been conducted to elucidate the mechanisms of their evolutionary adaptations. Using GO, KEGG, and MGI databases, warthogs showed a preference for expansion in sensory genes and contraction in metabolic genes, underscoring phenotypic diversity and adaptive evolution direction. Associating genes with the QTLdb-pigSS11 database revealed links between gene families and immunity traits. The overlap of olfactory genes in immune-related QTL regions highlighted their importance in evolutionary adaptations. This work highlights the unique evolutionary strategies and adaptive mechanisms of warthogs, guiding future research into the distinct adaptability and disease resistance in pigs, particularly focusing on traits such as resistance to African Swine Fever Virus.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

The genomes of Dahlia pinnata, Cosmos bipinnatus, and Bidens alba in tribe Coreopsideae provide insights into polyploid evolution and inulin biosynthesis

BACKGROUND

The Coreopsideae tribe, a subset of the Asteraceae family, encompasses economically vital genera like Dahlia, Cosmos, and Bidens, which are widely employed in medicine, horticulture, ecology, and food applications. Nevertheless, the lack of reference genomes hinders evolutionary and biological investigations in this tribe.

RESULTS

Here, we present 3 haplotype-resolved chromosome-level reference genomes of the tribe Coreopsideae, including 2 popular flowering plants (Dahlia pinnata and Cosmos bipinnatus) and 1 invasive weed plant (Bidens alba), with assembled genome sizes 3.93 G, 1.02 G, and 1.87 G, respectively. We found that Gypsy transposable elements contribute mostly to the larger genome size of D. pinnata, and multiple chromosome rearrangements have occurred in tribe Coreopsideae. Besides the shared whole-genome duplication (WGD-2) in the Heliantheae alliance, our analyses showed that D. pinnata and B. alba each underwent an independent recent WGD-3 event: in D. pinnata, it is more likely to be a self-WGD, while in B. alba, it is from the hybridization of 2 ancestor species. Further, we identified key genes in the inulin metabolic pathway and found that the pseudogenization of 1-FEH1 and 1-FEH2 genes in D. pinnata and the deletion of 3 key residues of 1-FFT proteins in C. bipinnatus and B. alba may probably explain why D. pinnata produces much more inulin than the other 2 plants.

CONCLUSIONS

Collectively, the genomic resources for the Coreopsideae tribe will promote phylogenomics in Asteraceae plants, facilitate ornamental molecular breeding improvements and inulin production, and help prevent invasive weeds.

Enterococcus montenegrensis sp. nov., isolated from artisanal Montenegrin dry sausage

A novel, Gram-positive, facultative anaerobe, coccoid and non-motile bacterium, designated as CoE-012-22T was isolated from dried beef sausage (the original name in Montenegro is Govedji Kulen) manufactured in the municipality of Rozaje (Montenegro) in 2021. Cells of this strain were oxidase- and catalase-negative. Growth occurred at 4-50 °C, at pH 5.0-8.0 and with 0-6.5 % (w/v) NaCl in diverse growth media. MALDI-TOF analysis identified the strain as Enterococcus canintestini (log score 2). Phylogenetic analysis of the 16S rRNA gene and whole genome sequences assigned the strain to the genus Enterococcus. The closest relatives were E. canintestini DSM 21207T and E. dispar ATCC 51266T with 16S rRNA gene sequence pairwise similarities of 99.34 and 98.59 %, respectively. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between isolate CoE-012-22T and other enterococci species were below the thresholds for species delineation thresholds (95.0 % ANI; 70.0 % dDDH) with maximum identities of 84.13 % (ANIb), 86.43 % (ANIm) and 28.4 % (dDDH) to E. saigonensis JCM 31193T and 70.97 % (ANIb), 88.99 % (ANIm) and 32.4 % (dDDH) to E. malodoratus ATCC 43197T. Two unknown Enterococcus isolates, Enterococcus sp. MJM12 and Enterococcus SMC-9, showed identities of 99.87 and 99.94 % (16S rRNA), 98.57 and 98.65 % (ANIb), 98.93 and 99.02 % (ANIm), and 89.8 and 90.0 % (dDDH) to strain CoE-012-22T and can therefore be regarded as the same species. Based on the characterization results, strain CoE-012-22T was considered to represent a novel species, for which the name Enterococcus montenegrensis sp. nov. is proposed. The type strain is CoE-012-22T (=DSM 115843T=NCIMB 15468T).

Spirochaete genome identified in red abalone sample represents a novel genus Candidatus Haliotispira gen. nov. within the order Spirochaetales

A fully assembled spirochaete genome was identified as a contaminating scaffold in our red abalone (Haliotis rufescens) genome assembly. In this paper, we describe the analysis of this bacterial genome. The assembled spirochaete genome is 3.25 Mb in size with 48.5 mol% G+C content. The proteomes of 38 species were compared with the spirochaete genome and it was discovered to form an independent branch within the family Spirochaetaceae on the phylogenetic tree. The comparison of 16S rRNA sequences and average nucleotide identity scores between the spirochaete genome with known species of different families in Spirochaetia indicate that it is an unknown species. Further, the percentage of conserved proteins compared to neighbouring taxa confirm that it does not belong to a known genus within Spirochaetaceae. We propose the name Candidatus Haliotispira prima gen. nov., sp. nov. based on its taxonomic placement and origin. We also tested for the presence of this species in different species of abalone and found that it is also present in white abalone (Haliotis sorenseni). In addition, we highlight the need for better classification of taxa within the class Spirochaetia.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae)

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

Identification of Incomplete Annotations of Biosynthesis Pathways in Rhodophytes Using a Multi-Omics Approach

Rhodophytes (red algae) are an important source of natural products and are, therefore, a current research focus in terms of metabolite production. The recent increase in publicly available Rhodophyte whole genome and transcriptome assemblies provides the resources needed for in silico metabolic pathway analysis. Thus, this study aimed to create a Rhodophyte multi-omics resource, utilising both genomes and transcriptome assemblies with functional annotations to explore Rhodophyte metabolism. The genomes and transcriptomes of 72 Rhodophytes were functionally annotated and integrated with metabolic reconstruction and phylogenetic inference, orthology prediction, and gene duplication analysis to analyse their metabolic pathways. This resource was utilised via two main investigations: the identification of bioactive sterol biosynthesis pathways and the evolutionary analysis of gene duplications for known enzymes. We report that sterol pathways, including campesterol, β-sitosterol, ergocalciferol and cholesterol biosynthesis pathways, all showed incomplete annotated pathways across all Rhodophytes despite prior in vivo studies showing otherwise. Gene duplication analysis revealed high rates of duplication of halide-associated haem peroxidases in Florideophyte algae, which are involved in the biosynthesis of drug-related halogenated secondary metabolites. In summary, this research revealed trends in Rhodophyte metabolic pathways that have been under-researched and require further functional analysis. Furthermore, the high duplication of haem peroxidases and other peroxidase enzymes offers insight into the potential drug development of Rhodophyte halogenated secondary metabolites.

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.

Comparative transcriptomics reveals divergence in pathogen response gene families amongst 20 forest tree species

Forest trees provide critical ecosystem services for humanity that are under threat due to ongoing global change. Measuring and characterizing genetic diversity are key to understanding adaptive potential and developing strategies to mitigate negative consequences arising from climate change. In the area of forest genetic diversity, genetic divergence caused by large-scale changes at the chromosomal level has been largely understudied. In this study, we used the RNA-seq data of 20 co-occurring forest trees species from genera including Acer, Alnus, Amelanchier, Betula, Cornus, Corylus, Dirca, Fraxinus, Ostrya, Populus, Prunus, Quercus, Ribes, Tilia, and Ulmus sampled from Upper Peninsula of Michigan. These data were used to infer the origin and maintenance of gene family variation, species divergence time, as well as gene family expansion and contraction. We identified a signal of common whole genome duplication events shared by core eudicots. We also found rapid evolution, namely fast expansion or fast contraction of gene families, in plant-pathogen interaction genes amongst the studied diploid species. Finally, the results lay the foundation for further research on the genetic diversity and adaptive capacity of forest trees, which will inform forest management and conservation policies.

Whole genome assembly and annotation of the endangered Caribbean coral Acropora cervicornis

Coral species in the genus Acropora are key ecological components of coral reefs worldwide and represent the most diverse genus of scleractinian corals. While key species of Indo-Pacific Acropora have annotated genomes, no annotated genome has been published for either of the two species of Caribbean Acropora. Here we present the first fully annotated genome of the endangered Caribbean staghorn coral, Acropora cervicornis. We assembled and annotated this genome using high-fidelity nanopore long-read sequencing with gene annotations validated with mRNA sequencing. The assembled genome size is 318 Mb, with 28,059 validated genes. Comparative genomic analyses with other Acropora revealed unique features in A. cervicornis, including contractions in immune pathways and expansions in signaling pathways. Phylogenetic analysis confirms previous findings showing that A. cervicornis diverged from Indo-Pacific relatives around 41 million years ago, with the closure of the western Tethys Sea, prior to the primary radiation of Indo-Pacific Acropora. This new A. cervicornis genome enriches our understanding of the speciose Acropora and addresses evolutionary inquiries concerning speciation and hybridization in this diverse clade.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Predation-resistant Pseudomonas bacteria engage in symbiont-like behavior with the social amoeba Dictyostelium discoideum

The soil amoeba Dictyostelium discoideum acts as both a predator and potential host for diverse bacteria. We tested fifteen Pseudomonas strains that were isolated from transiently infected wild D. discoideum for ability to escape predation and infect D. discoideum fruiting bodies. Three predation-resistant strains frequently caused extracellular infections of fruiting bodies but were not found within spores. Furthermore, infection by one of these species induces secondary infections and suppresses predation of otherwise edible bacteria. Another strain can persist inside of amoebae after being phagocytosed but is rarely taken up. We sequenced isolate genomes and discovered that predation-resistant isolates are not monophyletic. Many Pseudomonas isolates encode secretion systems and toxins known to improve resistance to phagocytosis in other species, as well as diverse secondary metabolite biosynthetic gene clusters that may contribute to predation resistance. However, the distribution of these genes alone cannot explain why some strains are edible and others are not. Each lineage may employ a unique mechanism for resistance.

What can be lost? Genomic perspective on the lipid metabolism of Mucoromycota

Mucoromycota is a phylum of early diverging fungal (EDF) lineages, of mostly plant-associated terrestrial fungi. Some strains have been selected as promising biotechnological organisms due to their ability to produce polyunsaturated fatty acids and efficient conversion of nutrients into lipids. Others get their lipids from the host plant and are unable to produce even the essential ones on their own. Following the advancement in EDF genome sequencing, we carried out a systematic survey of lipid metabolism protein families across different EDF lineages. This enabled us to explore the genomic basis of the previously documented ability to produce several types of lipids within the fungal tree of life. The core lipid metabolism genes showed no significant diversity in distribution, however specialized lipid metabolic pathways differed in this regard among different fungal lineages. In total 165 out of 202 genes involved in lipid metabolism were present in all tested fungal lineages, while remaining 37 genes were found to be absent in some of fungal lineages. Duplications were observed for 69 genes. For the first time we demonstrate that ergosterol is not being produced by several independent groups of plant-associated fungi due to the losses of different ERG genes. Instead, they possess an ancestral pathway leading to the synthesis of cholesterol, which is absent in other fungal lineages. The lack of diacylglycerol kinase in both Mortierellomycotina and Blastocladiomycota opens the question on sterol equilibrium regulation in these organisms. Early diverging fungi retained most of beta oxidation components common with animals including Nudt7, Nudt12 and Nudt19 pointing at peroxisome divergence in Dikarya. Finally, Glomeromycotina and Mortierellomycotina representatives have a similar set of desaturases and elongases related to the synthesis of complex, polyunsaturated fatty acids pointing at an ancient expansion of fatty acid metabolism currently being explored by biotechnological studies.

QR-STAR: A Polynomial-Time Statistically Consistent Method for Rooting Species Trees Under the Coalescent

We address the problem of rooting an unrooted species tree given a set of unrooted gene trees, under the assumption that gene trees evolve within the model species tree under the multispecies coalescent (MSC) model. Quintet Rooting (QR) is a polynomial time algorithm that was recently proposed for this problem, which is based on the theory developed by Allman, Degnan, and Rhodes that proves the identifiability of rooted 5-taxon trees from unrooted gene trees under the MSC. However, although QR had good accuracy in simulations, its statistical consistency was left as an open problem. We present QR-STAR, a variant of QR with an additional step and a different cost function, and prove that it is statistically consistent under the MSC. Moreover, we derive sample complexity bounds for QR-STAR and show that a particular variant of it based on "short quintets" has polynomial sample complexity. Finally, our simulation study under a variety of model conditions shows that QR-STAR matches or improves on the accuracy of QR. QR-STAR is available in open-source form on github.

Chromosome-level genome assembly and population genomic analysis provide novel insights into the immunity and evolution of Sogatella furcifera

Sogatella furcifera is a destructive agricultural pest causing large threats to rice production in China and Southeast Asian countries. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in S. furcifera. In present study, a chromosome-level assembly of the S. furcifera genome was completed (0.64 GB), comprising 15 chromosomes covered 95.04% of the estimated genome size, along with other 624 small scaffolds making up the remaining 4.96% of the genome of S. furcifera. A total of 24,669 protein-coding genes, 1211 long noncoding RNA and 7595 circular RNA transcripts were predicted in this study. Comparative genomic analysis revealed rapidly evolved genes were associated with multiple immune-related pathways in S. furcifera. Genome resequencing of 44 individuals from 12 geographic populations revealed frequent gene flow among populations. The systemic genomic analysis will provide more insights into the understanding of the immunity and evolutionary adaptation of S. furcifera.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

New reference genomes of honey bee-associated bacteria Paenibacillus melissococcoides, Paenibacillus dendritiformis, and Paenibacillus thiaminolyticus

We sequenced the genomes of recently discovered Paenibacillus melissococcoides (CCOS 2000) and of the type strains of closely related P. thiaminolyticus (DSM 7262) and P. dendritiformis (LMG 21716). The three genomes set the basis to unambiguous diagnostic of these honey bee associated Paenibacillus bacteria.

CACO: A Core-Attachment Method With Cross-Species Functional Ortholog Information to Detect Human Protein Complexes

Protein complexes play an essential role in living cells. Detecting protein complexes is crucial to understand protein functions and treat complex diseases. Due to high time and resource consumption of experiment approaches, many computational approaches have been proposed to detect protein complexes. However, most of them are only based on protein-protein interaction (PPI) networks, which heavily suffer from the noise in PPI networks. Therefore, we propose a novel core-attachment method, named CACO, to detect human protein complexes, by integrating the functional information from other species via protein ortholog relations. First, CACO constructs a cross-species ortholog relation matrix and transfers GO terms from other species as a reference to evaluate the confidence of PPIs. Then, a PPI filter strategy is adopted to clean the PPI network and thus a weighted clean PPI network is constructed. Finally, a new effective core-attachment algorithm is proposed to detect protein complexes from the weighted PPI network. Compared to other thirteen state-of-the-art methods, CACO outperforms all of them in terms of F-measure and Composite Score, showing that integrating ortholog information and the proposed core-attachment algorithm are effective in detecting protein complexes.

Phylogenetic placement of Ceratophyllum submersum based on a complete plastome sequence derived from nanopore long read sequencing data

OBJECTIVE

Eutrophication poses a mounting concern in today's world. Ceratophyllum submersum L. is one of many plants capable of living in eutrophic conditions, therefore it could play a critical role in addressing the problem of eutrophication. This study aimed to take a first genomic look at C. submersum.

RESULTS

Sequencing of gDNA from C. submersum yielded enough reads to assemble a plastome. Subsequent annotation and phylogenetic analysis validated existing information regarding angiosperm relationships and the positioning of Ceratophylalles in a wider phylogenetic context.

Sphingopyxis Species Isolated from Sand Filter Biofilm at an Australian Drinking Water Treatment Works

Three strains isolated by geosmin enrichment from a sand filter in an Australian drinking water treatment works were genome sequenced to identify their taxonomic placement, and a bench-scale batch experiment confirmed their geosmin-degrading capability. Using the average nucleotide identity based on the MUMmer algorithm (ANIm), pairwise digital DNA-DNA hybridization (dDDH), and phylogenomic analyses, the strains were identified as Sphingopyxis species.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

First Identification of bla NDM-1 Producing Escherichia coli ST 9499 Isolated from Musca domestica in the Urban Center of Rio de Janeiro, Brazil

The present study presents phenotypic and molecular characterization of a multidrug-resistant strain of Escherichia coli (Lemef26), belonging to sequence type ST9499 carrying a blaNDM-1 carbapenem resistance gene. The bacterium was isolated from a specimen of Musca domestica, collected in proximity to a hospital in Rio de Janeiro City, Brazil. The strain was identified as E. coli by matrix-assisted laser desorption-ionization time of flight mass spectrometry (Maldi-TOF-MS) and via genotypic analysis (Whole-Genome Sequencing-WGS), followed by phylogenetic analysis, antibiotic resistance profiling (using phenotypic and genotypic methods) and virulence genotyping. Interestingly, the blaNDM-1 was the only resistance determinant detected using a panel of common resistance genes, as evaluated by PCR. In contrast, WGS detected genes conferring resistance to aminoglycosides, fluoroquinolones, quinolones, trimethoprim, beta-lactams, chloramphenicol, macrolides, sulfonamide, tetracycline, lincosamide and streptogramin B. Conjugation experiments demonstrated the transfer of carbapenem resistance, via acquisition of the blaNDM-1 sequence, to a sensitive receptor strain of E. coli, indicating that blaNDM-1 is located on a conjugative plasmid (most likely of the IncA/C incompatibility group, in association with the transposon Tn3000). Phylogenetic analyses placed Lemef26 within a clade of strains exhibiting allelic and environment diversity, with the greatest level of relatedness recorded with a strain isolated from a human source suggesting a possible anthropogenic origin. Analysis of the virulome revealed the presence of fimbrial and pilus genes, including a CFA/I fimbriae (cfaABCDE), common pilus (ecpABCDER), laminin-bind fimbrae (elfADG), hemorrhagic pilus (hcpABC) and fimbrial adherence determinants (stjC) indicates the ability of strain Lemef26 to colonize animal hosts. To the best of our knowledge, this study represents the first report of blaNDM-1 carbapenemase gene in an E. coli strain isolated from M. domestica. In concordance with the findings of previous studies on the carriage of MDR bacteria by flies, the data presented herein provide support to the idea that flies may represent a convenient means (as sentinel animals) for the monitoring of environmental contamination with multidrug-resistant bacteria.

Genomic and Metabolomic Analysis of the Endophytic Fungus Fusarium sp. VM-40 Isolated from the Medicinal Plant Vinca minor

The genus Fusarium is well-known to comprise many pathogenic fungi that affect cereal crops worldwide, causing severe damage to agriculture and the economy. In this study, an endophytic fungus designated Fusarium sp. VM-40 was isolated from a healthy specimen of the traditional European medicinal plant Vinca minor. Our morphological characterization and phylogenetic analysis reveal that Fusarium sp. VM-40 is closely related to Fusarium paeoniae, belonging to the F. tricinctum species complex (FTSC), the genomic architecture and secondary metabolite profile of which have not been investigated. Thus, we sequenced the whole genome of Fusarium sp. VM-40 with the new Oxford Nanopore R10.4 flowcells. The assembled genome is 40 Mb in size with a GC content of 47.72%, 15 contigs (≥50,000 bp; N 50~4.3 Mb), and 13,546 protein-coding genes, 691 of which are carbohydrate-active enzyme (CAZyme)-encoding genes. We furthermore predicted a total of 56 biosynthetic gene clusters (BGCs) with antiSMASH, 25 of which showed similarity with known BGCs. In addition, we explored the potential of this fungus to produce secondary metabolites through untargeted metabolomics. Our analyses reveal that this fungus produces structurally diverse secondary metabolites of potential pharmacological relevance (alkaloids, peptides, amides, terpenoids, and quinones). We also employed an epigenetic manipulation method to activate cryptic BGCs, which led to an increased abundance of several known compounds and the identification of several putative new compounds. Taken together, this study provides systematic research on the whole genome sequence, biosynthetic potential, and metabolome of the endophytic fungus Fusarium sp. VM-40.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.