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.

Genome-wide identification and evolution of the tubulin gene family in Camelina sativa.

BACKGROUND: Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, ß- and γ-subfamilies, while α- and ß-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa. RESULTS: We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 ß- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome. CONCLUSIONS: For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and ß-tubulins.

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus).

Background: Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups. Methods: A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison. Results: Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa. Discussion: This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution.

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

Exploration of the pearl millet phospholipase gene family to identify potential candidates for grain quality traits.

BACKGROUND: Phospholipases constitute a diverse category of enzymes responsible for the breakdown of phospholipids. Their involvement in signal transduction with a pivotal role in plant development and stress responses is well documented. RESULTS: In the present investigation, a thorough genome-wide analysis revealed that the pearl millet genome contains at least 44 phospholipase genes distributed across its 7 chromosomes, with chromosome one harbouring the highest number of these genes. The synteny analysis suggested a close genetic relationship of pearl millet phospholipases with that of foxtail millet and sorghum. All identified genes were examined to unravel their gene structures, protein attributes, cis-regulatory elements, and expression patterns in two pearl millet genotypes contrasting for rancidity. All the phospholipases have a high alpha-helix content and distorted regions within the predicted secondary structures. Moreover, many of these enzymes possess binding sites for both metal and non-metal ligands. Additionally, the putative promoter regions associated with these genes exhibit multiple copies of cis-elements specifically responsive to biotic and abiotic stress factors and signaling molecules. The transcriptional profiling of 44 phospholipase genes in two genotypes contrasting for rancidity across six key tissues during pearl millet growth revealed a predominant expression in grains, followed by seed coat and endosperm. Specifically, the genes PgPLD-alpha1-1, PgPLD-alpha1-5, PgPLD-delta1-7a, PgPLA1-II-1a, and PgPLD-delta1-2a exhibited notable expression in grains of both the genotypes while showing negligible expression in the other five tissues. The sequence alignment of putative promoters revealed several variations including SNPs and InDels. These variations resulted in modifications to the corresponding cis-acting elements, forming distinct transcription factor binding sites suggesting the transcriptional-level regulation for these five genes in pearl millet. CONCLUSIONS: The current study utilized a genome-wide computational analysis to characterize the phospholipase gene family in pearl millet. A comprehensive expression profile of 44 phospholipases led to the identification of five grain-specific candidates. This underscores a potential role for at least these five genes in grain quality traits including the regulation of rancidity in pearl millet. Therefore, this study marks the first exploration highlighting the possible impact of phospholipases towards enhancing agronomic traits in pearl millet.

The Expansion of Sirtuin Gene Family in Gilthead Sea Bream (Sparus aurata)-Phylogenetic, Syntenic, and Functional Insights across the Vertebrate/Fish Lineage.

The Sirtuin (SIRT1-7) family comprises seven evolutionary-conserved enzymes that couple cellular NAD availability with health, nutrition and welfare status in vertebrates. This study re-annotated the sirt3/5 branch in the gilthead sea bream, revealing three paralogues of sirt3 (sirt3.1a/sirt3.1b/sirt3.2) and two of sirt5 (sirt5a/sirt5b) in this Perciform fish. The phylogeny and synteny analyses unveiled that the Sirt3.1/Sirt3.2 dichotomy was retained in teleosts and aquatic-living Sarcopterygian after early vertebrate 2R whole genome duplication (WGD). Additionally, only certain percomorphaceae and gilthead sea bream showed a conserved tandem-duplicated synteny block involving the mammalian-clustered sirt3.1 gene (psmd13-sirt3.1a/b-drd4-cdhr5-ctsd). Conversely, the expansion of the Sirt5 branch was shaped by the teleost-specific 3R WGD. As extensively reviewed in the literature, human-orthologues (sirt3.1/sirt5a) showed a high, conserved expression in skeletal muscle that increased as development advanced. However, recent sirt3.2 and sirt5b suffered an overall muscle transcriptional silencing across life, as well as an enhanced expression on immune-relevant tissues and gills. These findings fill gaps in the ontogeny and differentiation of Sirt genes in the environmentally adaptable gilthead sea bream, becoming a good starting point to advance towards a full understanding of its neo-functionalization. The mechanisms originating from these new paralogs also open new perspectives in the study of cellular energy sensing processes in vertebrates.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis.

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety.

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Genome-wide identification of GH28 family and insight into its contributions to pod shattering resistance in Brassica napus L.

Rapeseed (Brassica napus L.), accounts for nearly 16% of vegetable oil, is the world's second produced oilseed. However, pod shattering has caused significant yield loses in rapeseed production, particularly during mechanical harvesting. The GH28 genes can promote pod shattering by changing the structure of the pod cell wall in Arabidopsis. However, the role of the GH28 gene family in rapeseed was largely unknown. Therefore, a genome-wide comprehensive analysis was conducted to classify the role of GH28 gene family on rapeseed pod shattering. A total of 37 BnaGH28 genes in the rapeseed genome were identified. These BnaGH28s can be divided into five groups (Group A-E), based on phylogenetic and synteny analysis. Protein property, gene structure, conserved motif, cis-acting element, and gene expression profile of BnaGH28 genes in the same group were similar. Specially, the expression level of genes in group A-D was gradually decreased, but increased in group E with the development of silique. Among eleven higher expressed genes in group E, two BnaGH28 genes (BnaA07T0199500ZS and BnaC06T0206500ZS) were significantly regulated by IAA or GA treatment. And the significant effects of BnaA07T0199500ZS variation on pod shattering resistance were also demonstrated in present study. These results could open a new window for insight into the role of BnaGH28 genes on pod shattering resistance in rapeseed.

AGO, a Framework for the Reconstruction of Ancestral Syntenies and Gene Orders.

Reconstructing ancestral gene orders from the genome data of extant species is an important problem in comparative and evolutionary genomics. In a phylogenomics setting that accounts for gene family evolution through gene duplication and gene loss, the reconstruction of ancestral gene orders involves several steps, including multiple sequence alignment, the inference of reconciled gene trees, and the inference of ancestral syntenies and gene adjacencies. For each of the steps of such a process, several methods can be used and implemented using a growing corpus of, often parameterized, tools; in practice, interfacing such tools into an ancestral gene order reconstruction pipeline is far from trivial. This chapter introduces AGO, a Python-based framework aimed at creating ancestral gene order reconstruction pipelines allowing to interface and parameterize different bioinformatics tools. The authors illustrate the features of AGO by reconstructing ancestral gene orders for the X chromosome of three ancestral Anopheles species using three different pipelines. AGO is freely available at https://github.com/cchauve/AGO-pipeline .

GhiPLATZ17 and GhiPLATZ22, zinc-dependent DNA-binding transcription factors, promote salt tolerance in upland cotton.

KEY MESSAGE: The utilization of transcriptome analysis, functional validation, VIGS, and DAB techniques have provided evidence that GhiPLATZ17 and GhiPLATZ22 play a pivotal role in improving the salt tolerance of upland cotton. PLATZ (Plant AT-rich sequences and zinc-binding proteins) are known to be key regulators in plant growth, development, and response to salt stress. In this study, we comprehensively analyzed the PLATZ family in ten cotton species in response to salinity stress. Gossypium herbaceum boasts 25 distinct PLATZ genes, paralleled by 24 in G. raimondii, 25 in G. arboreum, 46 in G. hirsutum, 48 in G. barbadense, 43 in G. tomentosum, 67 in G. mustelinum, 60 in G. darwinii, 46 in G. ekmanianum, and a total of 53 PLATZ genes attributed to G. stephensii. The PLATZ gene family shed light on the hybridization and allopolyploidy events that occurred during the evolutionary history of allotetraploid cotton. Ka/Ks analysis suggested that the PLATZ gene family underwent intense purifying selection during cotton evolution. Analysis of synteny and gene collinearity revealed a complex pattern of segmental and dispersed duplication events to expand PLATZ genes in cotton. Cis-acting elements and gene expressions revealed that GhiPLATZ exhibited salt stress resistance. Transcriptome analysis, functional validation, virus-induced gene silencing (VIGS), and diaminobenzidine staining (DAB) demonstrated that GhiPLATZ17 and GhiPLATZ22 enhance salt tolerance in upland cotton. The study can potentially advance our understanding of identifying salt-resistant genes in cotton.

SAFPred: synteny-aware gene function prediction for bacteria using protein embeddings.

MOTIVATION: Today, we know the function of only a small fraction of the protein sequences predicted from genomic data. This problem is even more salient for bacteria, which represent some of the most phylogenetically and metabolically diverse taxa on Earth. This low rate of bacterial gene annotation is compounded by the fact that most function prediction algorithms have focused on eukaryotes, and conventional annotation approaches rely on the presence of similar sequences in existing databases. However, often there are no such sequences for novel bacterial proteins. Thus, we need improved gene function prediction methods tailored for bacteria. Recently, transformer-based language models-adopted from the natural language processing field-have been used to obtain new representations of proteins, to replace amino acid sequences. These representations, referred to as protein embeddings, have shown promise for improving annotation of eukaryotes, but there have been only limited applications on bacterial genomes. RESULTS: To predict gene functions in bacteria, we developed SAFPred, a novel synteny-aware gene function prediction tool based on protein embeddings from state-of-the-art protein language models. SAFpred also leverages the unique operon structure of bacteria through conserved synteny. SAFPred outperformed both conventional sequence-based annotation methods and state-of-the-art methods on multiple bacterial species, including for distant homolog detection, where the sequence similarity to the proteins in the training set was as low as 40%. Using SAFPred to identify gene functions across diverse enterococci, of which some species are major clinical threats, we identified 11 previously unrecognized putative novel toxins, with potential significance to human and animal health. AVAILABILITY AND IMPLEMENTATION: https://github.com/AbeelLab/safpred.

The promise and pitfalls of synteny in phylogenomics.

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Genome assembly of Melilotus officinalis provides a new reference genome for functional genomics.

BACKGROUND: Sweet yellow clover (Melilotus officinalis) is a diploid plant (2n = 16) that is native to Europe. It is an excellent legume forage. It can both fix nitrogen and serve as a medicine. A genome assembly of Melilotus officinalis that was collected from Best corporation in Beijing is available based on Nanopore sequencing. The genome of Melilotus officinalis was sequenced, assembled, and annotated. RESULTS: The latest PacBio third generation HiFi assembly and sequencing strategies were used to produce a Melilotus officinalis genome assembly size of 1,066 Mbp, contig N50 = 5 Mbp, scaffold N50 = 130 Mbp, and complete benchmarking universal single-copy orthologs (BUSCOs) = 96.4%. This annotation produced 47,873 high-confidence gene models, which will substantially aid in our research on molecular breeding. A collinear analysis showed that Melilotus officinalis and Medicago truncatula shared conserved synteny. The expansion and contraction of gene families showed that Melilotus officinalis expanded by 565 gene families and shrank by 56 gene families. The contacted gene families were associated with response to stimulus, nucleotide binding, and small molecule binding. Thus, it is related to a family of genes associated with peptidase activity, which could lead to better stress tolerance in plants. CONCLUSIONS: In this study, the latest PacBio technology was used to assemble and sequence the genome of the Melilotus officinalis and annotate its protein-coding genes. These results will expand the genomic resources available for Melilotus officinalis and should assist in subsequent research on sweet yellow clover plants.

OrthoRefine: automated enhancement of prior ortholog identification via synteny.

BACKGROUND: Identifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification. RESULTS: We developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC. CONCLUSION: OrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

Extensive Loss and Gain of Conserved Noncoding Elements During Early Teleost Evolution.

Conserved noncoding elements in vertebrates are enriched around transcription factor loci associated with development. However, loss and rapid divergence of conserved noncoding elements has been reported in teleost fish, albeit taking only few genomes into consideration. Taking advantage of the recent increase in high-quality teleost genomes, we focus on studying the evolution of teleost conserved noncoding elements, carrying out targeted genomic alignments and comparisons within the teleost phylogeny to detect conserved noncoding elements and reconstruct the ancestral teleost conserved noncoding elements repertoire. This teleost-centric approach confirms previous observations of extensive vertebrate conserved noncoding elements loss early in teleost evolution, but also reveals massive conserved noncoding elements gain in the teleost stem-group over 300 million years ago. Using synteny-based association to link conserved noncoding elements to their putatively regulated target genes, we show the most teleost gained conserved noncoding elements are found in the vicinity of orthologous loci involved in transcriptional regulation and embryonic development that are also associated with conserved noncoding elements in other vertebrates. Moreover, teleost and vertebrate conserved noncoding elements share a highly similar motif and transcription factor binding site vocabulary. We suggest that early teleost conserved noncoding element gains reflect a restructuring of the ancestral conserved noncoding element repertoire through both extreme divergence and de novo emergence. Finally, we support newly identified pan-teleost conserved noncoding elements have potential for accurate resolution of teleost phylogenetic placements in par with coding sequences, unlike ancestral only elements shared with spotted gar. This work provides new insight into conserved noncoding element evolution with great value for follow-up work on phylogenomics, comparative genomics, and the study of gene regulation evolution in teleosts.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species.

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Whole-genome assembly of a novel invertebrate herpesvirus from the gastropod Babylonia areolata.

Molluscan herpesviruses cause disease in species of major importance to aquaculture and are the only known herpesviruses to infect invertebrates, which lack an adaptive immune system. Understanding the evolution of malacoherpesviruses in relation to their hosts will likely require comparative genomic studies on multiple phylogenetic scales. Currently, only two malacoherpesvirus species have genomes that have been fully assembled, which limits the ability to perform comparative genomic studies on this family of viruses. In the present study, we fully assemble a herpesvirus from Illumina and Nanopore sequence data that were previously used to assemble the genome of the gastropod Babylonia areolata. We tentatively assign this novel herpesvirus to the genus Aurivirus within the family Malacoherpesviridae based on a phylogenetic analysis of DNA polymerase. While structurally similar to other malacoherpesvirus genomes, a synteny analysis of the novel herpesvirus with another Aurivirus species indicates that genomic rearrangements might be an important process in the evolution of this genus. We anticipate that future complete assemblies of malacoherpesviruses will be a valuable resource in comparative herpesvirus research.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa.

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence.

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.