Calcisponges have a ParaHox gene and dynamic expression of dispersed NK homeobox genes

This study was funded by the Sars Centre core budget to M. Adamska. Sequencing was performed at the Norwegian High Throughput Sequencing Centre funded by the Norwegian Research Council. O.M.R. and D.E.K.F. acknowledge support from the BBSRC and the School of Biology, University of St Andrews.Sponges are simple animals with few cell types, but their genomes paradoxically contain a wide variety of developmental transcription factors1,2,3,4, including homeobox genes belonging to the Antennapedia (ANTP) class5,6, which in bilaterians encompass Hox, ParaHox and NK genes. In the genome of the demosponge Amphimedon queenslandica, no Hox or ParaHox genes are present, but NK genes are linked in a tight cluster similar to the NK clusters of bilaterians5. It has been proposed that Hox and ParaHox genes originated from NK cluster genes after divergence of sponges from the lineage leading to cnidarians and bilaterians5,7. On the other hand, synteny analysis lends support to the notion that the absence of Hox and ParaHox genes in Amphimedon is a result of secondary loss (the ghost locus hypothesis)8. Here we analysed complete suites of ANTP-class homeoboxes in two calcareous sponges, Sycon ciliatum and Leucosolenia complicata. Our phylogenetic analyses demonstrate that these calcisponges possess orthologues of bilaterian NK genes (Hex, Hmx and Msx), a varying number of additional NK genes and one ParaHox gene, Cdx. Despite the generation of scaffolds spanning multiple genes, we find no evidence of clustering of Sycon NK genes. All Sycon ANTP-class genes are developmentally expressed, with patterns suggesting their involvement in cell type specification in embryos and adults, metamorphosis and body plan patterning. These results demonstrate that ParaHox genes predate the origin of sponges, thus confirming the ghost locus hypothesis8, and highlight the need to analyse the genomes of multiple sponge lineages to obtain a complete picture of the ancestral composition of the first animal genome.PostprintPeer reviewe

Chromosomal‐level reference genome of the incense tree Aquilaria sinensis

This is the author accepted manuscript. The final version is available on open access from Wiley via the DOI in this recordData Accessibility: The raw genome and RNA sequencing data have been deposited in the SRA under Bioproject numbers SRR10737433 and PRJNA534170. The final chromosome assembly was submitted to NCBI Assembly under accession number VZPZ00000000 in NCBI.Trees in the genus Aquilaria (Thymelaeaceae) are known as lign aloes, and are native to the forests of southeast Asia. Lign aloes produce agarwood as an antimicrobial defence. Agarwood has a long history of cultural and medicinal use, and is of considerable commercial value. However, due to habitat destruction and over collection, lign aloes are threatened in the wild. We present a chromosomal‐level assembly for Aquilaria sinensis, a lign aloe endemic to China known as the incense tree, based on Illumina short‐read, 10X Genomics linked‐read, and Hi‐C sequencing data. Our 783.8Mbp A. sinensis genome assembly is of high physical contiguity, with a scaffold N50 of 87.6Mbp, and high completeness, with a 95.8% BUSCO score for eudicotyledon genes. We include 17 transcriptomes from various plant tissues, providing a total of 35,965 gene models. We reveal the first complete set of genes involved in sesquiterpenoid production, plant defence, and agarwood production for the genus Aquilaria, including genes involved in the biosynthesis of sesquiterpenoids via the mevalonic acid (MVA), 1‐deoxy‐D‐xylulose‐5‐phosphate (DXP), and methylerythritol phosphate (MEP) pathways. We perform a detailed repeat content analysis, revealing that transposable elements account for ~61% of the genome, with major contributions from gypsy‐like and copia‐like LTR retroelements. We also provide a comparative analysis of repeat content across sequenced species in the order Malvales. Our study reveals the first chromosomal‐level genome assembly for a tree in the genus Aquilaria and provides an unprecedented opportunity to address a variety of applied, genomic and evolutionary questions in the Thymelaeaceae more widely.Chinese University of Hong KongAgriculture, Fisheries and Conservation Department, Government of the Hong Kong Special Administrative RegionBiotechnology and Biological Sciences Research Council (BBSRC

Millipede genomes reveal unique adaptations during myriapod evolution

This is the final version. Available from Public Library of Science via the DOI in this record.The final genome assemblies have been deposited to NCBI database with accession numbers JAAFCF000000000 and JAAFCE000000000. The mRNA and sRNA transcriptomic data generated in this study have been deposited to the NCBI database under the following BioProject accessions: PRJNA564202 (Helicorthomorpha holstii) and PRJNA564195 (Trigoniulus corallinus).The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs-including species-specific microRNA arm switching-providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes-and in particular the reconstruction of the myriapod ancestral situation-for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.Hong Kong Research Grants Council (RGC) General Research FundHong Kong Research Grants Council (RGC) General Research FundThe Chinese University of Hong Kong (CUHK

Chromosomal-level reference genome of the moth Heortia vitessoides (Lepidoptera: Crambidae), a major pest of agarwood-producing trees

This is the final version. Available on open access from Elsevier via the DOI in this recordAvailability of data and materials: The final chromosome assembly was submitted to NCBI Assembly under accession number JACJUM000000000 in NCBI. The raw reads generated in this study have been deposited to the NCBI database under the BioProject accessions: PRJNA654728, the genome annotation files were deposited in the Figshare (https://doi.org/10.6084/m9.figshare.19633668). The microRNA sequences of known species were obtained from both miRbase [42] and MirGeneDB [21].The moth Heortia vitessoides Moore (Lepidoptera: Crambidae) is a major pest of ecologically, commercially and culturally important agarwood-producing trees in the genus Aquilaria. In particular, H. vitessoides is one of the most destructive defoliating pests of the incense tree Aquilaria sinesis, which produces a valuable fragrant wood used as incense and in traditional Chinese medicine [33]. Nevertheless, a genomic resource for H. vitessoides is lacking. Here, we present a chromosomal-level assembly for H. vitessoides, consisting of a 517 megabase (Mb) genome assembly with high physical contiguity (scaffold N50 of 18.2 Mb) and high completeness (97.9% complete BUSCO score). To aid gene annotation, 8 messenger RNA transcriptomes from different developmental stages were generated, and a total of 16,421 gene models were predicted. Expansion of gene families involved in xenobiotic metabolism and development were detected, including duplications of cytosolic sulfotransferase (SULT) genes shared among lepidopterans. In addition, small RNA sequencing of 5 developmental stages of H. vitessoides facilitated the identification of 85 lepidopteran conserved microRNAs, 94 lineage-specific microRNAs, as well as several microRNA clusters. A large proportion of the H. vitessoides genome consists of repeats, with a 29.12% total genomic contribution from transposable elements, of which long interspersed nuclear elements (LINEs) are the dominant component (17.41%). A sharp decrease in the genome-wide percentage of LINEs with lower levels of genetic distance to family consensus sequences suggests that LINE activity has peaked in H. vitessoides. In contrast, opposing patterns suggest a substantial recent increase in DNA and LTR element activity. Together with annotations of essential sesquiterpenoid hormonal pathways, neuropeptides, microRNAs and transposable elements, the high-quality genomic and transcriptomic resources we provide for the economically important moth H. vitessoides provide a platform for the development of genomic approaches to pest management, and contribute to addressing fundamental research questions in Lepidoptera.Hong Kong Research Grant Council Collaborative Research FundGeneral Research FundChinese University of Hong Kon

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe

This is the final version. Available on open access from BMC via the DOI in this recordAvailability of data and materials: The raw reads generated in this study have been deposited to the NCBI database under the BioProject accession PRJNA664668 [110]. The final chromosome assembly was submitted to NCBI Assembly under accession number JADANM000000000 in NCBI [111]. The genome annotation files were deposited in the Figshare [112].BACKGROUND: Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored. RESULTS: Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies. CONCLUSIONS: Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.Hong Kong Research Grant Council Collaborative Research Fund CRFGeneral Research Fund GRFArea of ExcellenceChinese University of Hong KongBiotechnology and Biological Sciences Research Council (BBSRC

Myriapod genomes reveal ancestral horizontal gene transfer and hormonal gene loss in millipedes

This is the final version. Available on open access from Nature Research via the DOI in this recordData availability; The final assemblies were submitted to NCBI Assembly under accession numbers WWPM00000000 (Glomeris maerens), JAAFCF000000000 (Helicorthomorpha holstii), WWPL00000000 (Anaulaciulus tonginus), JAAIVG000000000 (Niponia nodulosa), JAAFCE000000000 (Trigoniulus corallinus), JAHWFP000000000 (Lithobius niger), JAHWFO000000000 (Rhysida immarginata) and JAFIDM000000000 (Thereuonema tuberculata) in NCBI. The raw reads generated in this study were deposited to the NCBI database under the BioProject accessions PRJNA598061 (Glomeris maerens), PRJNA564202 (Helicorthomorpha holstii), PRJNA598060 (Anaulaciulus tonginus), PRJNA606398 (Niponia nodulosa), PRJNA564195 (Trigoniulus corallinus), PRJNA738717 (Lithobius niger), PRJNA701115 (Rhysida immarginata) and PRJNA699399 (Thereuonema tuberculata). The genome annotation files were deposited in the Figshare (https://doi.org/10.6084/m9.figshare.15088722). The databases are available for download from the following websites: eggNOG http://eggnog5.embl.de/download/eggnog_5.0/, GO http://geneontology.org/, KEGG https://www.genome.jp/kegg/pathway.html, and KOG https://www.hsls.pitt.edu/obrc/index.php?page=URL1144075392. Source data are provided with this paper.Code availability: The scripts for carrying out analyses of this study were deposited in Zenodo: https://doi.org/10.5281/zenodo.571873479 and https://doi.org/10.5281/zenodo.648262594.Animals display a fascinating diversity of body plans. Correspondingly, genomic analyses have revealed dynamic evolution of gene gains and losses among animal lineages. Here we sequence six new myriapod genomes (three millipedes, three centipedes) at key phylogenetic positions within this major but understudied arthropod lineage. We combine these with existing genomic resources to conduct a comparative analysis across all available myriapod genomes. We find that millipedes generally have considerably smaller genomes than centipedes, with the repeatome being a major contributor to genome size, driven by independent large gains of transposons in three centipede species. In contrast to millipedes, centipedes gained a large number of gene families after the subphyla diverged, with gains contributing to sensory and locomotory adaptations that facilitated their ecological shift to predation. We identify distinct horizontal gene transfer (HGT) events from bacteria to millipedes and centipedes, with no identifiable HGTs shared among all myriapods. Loss of juvenile hormone O-methyltransferase, a key enzyme in catalysing sesquiterpenoid hormone production in arthropods, was also revealed in all millipede lineages. Our findings suggest that the rapid evolution of distinct genomic pathways in centipede and millipede lineages following their divergence from the myriapod ancestor, was shaped by differing ecological pressures

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasi

This is the final version. Available on open access from Oxford University Press via the DOI in this recordData Availability: The raw genome and RNA sequencing data have been deposited in the SRA under Bioproject No. PRJNA673593. The final chromosome assembly was submitted to NCBI Assembly under accession No. JADKLZ000000000. All data can also be found in the GigaScience Database [75].BACKGROUND: Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea. FINDINGS: Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage. CONCLUSION: This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely

Sex-Biased Expression of MicroRNAs in Schistosoma mansoni

Schistosomiasis is an important neglected tropical disease caused by digenean helminth parasites of the genus Schistosoma. Schistosomes are unusual in that they are dioecious and the adult worms live in the blood system. MicroRNAs play crucial roles during gene regulation and are likely to be important in sex differentiation in dioecious species. Here we characterize 112 microRNAs from adult Schistosoma mansoni individuals, including 84 novel microRNA families, and investigate the expression pattern in different sexes. By deep sequencing, we measured the relative expression levels of conserved and newly identified microRNAs between male and female samples. We observed that 13 microRNAs exhibited sex-biased expression, 10 of which are more abundant in females than in males. Sex chromosomes showed a paucity of female-biased genes, as predicted by theoretical evolutionary models. We propose that the recent emergence of separate sexes in Schistosoma had an effect on the chromosomal distribution and evolution of microRNAs, and that microRNAs are likely to participate in the sex differentiation/maintenance process