Chromosome-scale Genome Assembly Research Articles

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully functional new organs, including new brains, de novo. The regeneration of a new brain requires the formation of diverse neural cell types and their assembly into an organized structure with correctly wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neural subpopulations, however, how these transcriptional programs are initiated in response to injury remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia, to study wound-induced transcriptional regulatory events that lead to the production of neurons and subsequently a functional brain. Footprinting analysis using chromatin accessibility data on a chromosome-scale genome assembly revealed that binding sites for the Nuclear Factor Y (NFY) transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal function. Single-cell transcriptome analysis combined with functional studies identified soxC+ stem cells as a putative progenitor population for multiple neural subtypes. Further, we found that wound-induced soxC expression is likely under direct transcriptional control by NFY, uncovering a mechanism for the initiation of a neural differentiation pathway by early wound-induced binding of a transcriptional regulator.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

A chromosome-scale genome assembly of the nipa palm hispid beetle Octodonta nipae

Nipa palm hispid beetle (Octodonta nipae) is an insect species that is native to Malaysia but has spread to southern China and beyond, seriously threatening palm production. A lack of high-quality genome resources has hindered understanding of the insect’s invasive characteristics and ecological adaptations. Here, we combined Illumina short read, PacBio long-read, and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to generate a high-quality, chromosome-scale genome assembly of nipa palm hispid beetle. The genome assembly was 1.31 Gb in size, consisting of nine chromosomes. The contig and scaffold N50 values were 1.022 Mb and 148.6 Mb, respectively. The genome assembly completeness was estimated at 99.1%. Annotation revealed 16,305 protein-coding genes and 62.16% repeat sequences. This high-quality genome assembly is a valuable resource that will contribute to understanding of the genetic factors underlying the invasive characteristics of nipa palm hispid beetle, ultimately promoting development of efficient control policies.

Chromosome-scale Genome Assembly and Characterization of Top-quality Japanese Green Tea Cultivar 'Seimei'.

Japanese green tea, an essential beverage in Japanese culture, is characterized by the initial steaming of freshly harvested leaves during production. This process efficiently inactivates endogenous enzymes such as polyphenol oxidases, resulting in the production of sencha, gyokuro, and matcha that preserves the vibrant green color of young leaves. Although genome sequences of several tea cultivars and germplasms have been published, no reference genome sequences are available for Japanese green tea cultivars. Here, we constructed a reference genome sequence of the cultivar 'Seimei', which is used to produce high-quality Japanese green tea. Using the PacBio HiFi and Hi-C technologies for chromosome-scale genome assembly, we obtained 15 chromosome sequences with a total genome size of 3.1 Gb and an N50 of 214.9 Mb. By analyzing the genomic diversity of 23 Japanese tea cultivars and lines, including the leading green tea cultivars 'Yabukita' and 'Saemidori', revealed several candidate genes that could be related to the characteristics of Japanese green tea. The reference genome of 'Seimei' and information on genomic diversity of Japanese green tea cultivars should provide crucial information for effective breeding of such cultivars in the future.

A fully haplotype-resolved and nearly gap-free genome assembly of wheat stripe rust fungus

Stripe rust fungus Puccinia striiformis f. sp. tritici (Pst) is a destructive pathogen of wheat worldwide. Pst has a macrocyclic-heteroecious lifecycle, in which one-celled urediniospores are dikaryotic, each nucleus containing one haploid genome. We successfully generated the first fully haplotype-resolved and nearly gap-free chromosome-scale genome assembly of Pst by combining PacBio HiFi sequencing and trio-binning strategy. The genome size of the two haploid assemblies was 75.59 Mb and 75.91 Mb with contig N50 of 4.17 Mb and 4.60 Mb, and both had 18 pseudochromosomes. The high consensus quality values of 55.57 and 59.02 for both haplotypes confirmed the correctness of the assembly. Of the total 18 chromosomes, 15 and 16 were gapless while there were only five and two gaps for the remaining chromosomes of the two haplotypes, respectively. In total, 15,046 and 15,050 protein-coding genes were predicted for the two haplotypes, and the complete BUSCO scores achieved 97.7% and 97.9%, respectively. The genome will lay the foundation for further research on genetic variations and the evolution of rust fungi.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named “Guangchenpi”, has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3’,4’-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B.juncea (AABB).

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B.carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

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.

Chromosome-scale genome assembly clarifies the mechanism of flooding tolerance and evolutionary history of Myricaria laxiflora

Myricaria laxiflora, belonging to the Tamaricaceae family, is a valuable plant resource naturally occurring along the riverbanks of the Yangtze River Basin. In this study, we constructed a chromosome-level M. laxiflora genome which was 1.34 Gb in size, with 12 chromosomes and a contig N50 of 11.69 Mb. In total, 23,719 protein-coding genes were discovered across the M. laxiflora genome. Phylogenetic investigation indicated that Tamaricaceae and Polygonaceae diverged between 65 and 81 million years ago (Mya). M. laxiflora underwent one whole genome duplication (WGD) between 19.74 and 24.60 Mya, and experienced a long terminal repeat (LTR) insertion around 0.15 Mya. These genomic changes may be related to environmental upheavals (such as the uplift of the Qinghai-Tibet Plateau) experienced by M. laxiflora during its evolution, and also may have provided a foundation for its evolution and adaptation. To analyze the flooding adaption mechanism of M. laxiflora, we compared morphological, physiological, and transcriptomic changes under control and submerged conditions. Overall, M. laxiflora was found to adopt a "quiescence" strategy when completely submerged. Specifically, in response to submersion, photosynthesis weakens, growth stagnates, glycolysis is activated, and the antioxidant enzyme system is enhanced. The “ethylene response pathway” was ascertained as a likely regulator of this process, and a MlERF-VII gene (Mla11G026380) may be the key gene regulating this process. Our findings shed light on the underlying mechanisms behind the flooding resistance displayed by M. laxiflora, while also establishing a basis for the conservation of its genetic resources.

Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps.

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

Chromosome-level genome assembly of the predatory stink bug Arma custos

The stink bug Arma custos (Hemiptera: Pentatomidae) is a predatory enemy successfully used for biocontrol of lepidopteran and coleopteran pests in notorious invasive species. In this study, a high-quality chromosome-scale genome assembly of A. custos was achieved through a combination of Illumina sequencing, PacBio HiFi sequencing, and Hi-C scaffolding techniques. The final assembled genome was 969.02 Mb in size, with 935.94 Mb anchored to seven chromosomes, and a scaffold N50 length of 135.75 Mb. This genome comprised 52.78% repetitive elements. The detected complete BUSCO score was 99.34%, indicating its completeness. A total of 13,708 protein-coding genes were predicted in the genome, and 13219 of them were annotated. This genome provides an invaluable resource for further research on various aspects of predatory bugs, such as biology, genetics, and functional genomics.

The genomic and enzymatic basis for iridoid biosynthesis in cat thyme (Teucrium marum).

Iridoids are non-canonical monoterpenoids produced by both insects and plants. An example is the cat-attracting and insect-repelling volatile iridoid nepetalactone, produced by Nepeta sp. (catmint) and aphids. Recently, both nepetalactone biosynthetic pathways were elucidated, showing a remarkable convergent evolution. The iridoid, dolichodial, produced by Teucrium marum (cat thyme) and multiple insect species, has highly similar properties to nepetalactone but its biosynthetic origin remains unknown. We set out to determine the genomic, enzymatic, and evolutionary basis of iridoid biosynthesis in T. marum. First, we generated a de novo chromosome-scale genome assembly for T. marum using Oxford Nanopore Technologies long reads and proximity-by-ligation Hi-C reads. The 610.3 Mb assembly spans 15 pseudomolecules with a 32.9 Mb N50 scaffold size. This enabled identification of iridoid biosynthetic genes, whose roles were verified via activity assays. Phylogenomic analysis revealed that the evolutionary history of T. marum iridoid synthase, the iridoid scaffold-forming enzyme, is not orthologous to typical iridoid synthases but is derived from its conserved paralog. We discovered an enzymatic route from nepetalactol to diverse iridoids through the coupled activity of an iridoid oxidase cytochrome P450 and acetyltransferases, via an inferred acylated intermediate. This work provides a genomic resource for specialized metabolite research in mints and demonstration of the role of acetylation in T. marum iridoid diversity. This work will enable future biocatalytic or biosynthetic production of potent insect repellents, as well as comparative studies into iridoid biosynthesis in insects.

Genome assembly and multi-omic analyses reveal the mechanisms underlying flower color formation in Torenia fournieri.

Torenia fournieri Lind. is an ornamental plant that is popular for its numerous flowers and variety of colors. However, its genomic evolutionary history and the genetic and metabolic bases of flower color formation remain poorly understood. Here, we report the first T. fournieri reference genome, which was resolved to the chromosome scale and was 164.4 Mb in size. Phylogenetic analyses clarified relationships with other plant species, and a comparative genomic analysis indicated that the shared ancestor of T. fournieri and Antirrhinum majus underwent a whole genome duplication event. Joint transcriptomic and metabolomic analyses identified many metabolites related to pelargonidin, peonidin, and naringenin production in rose (TfR)-colored flowers. Samples with blue (TfB) and deep blue (TfD) colors contained numerous derivatives of petunidin, cyanidin, quercetin, and malvidin; differences in the abundances of these metabolites and expression levels of the associated genes were hypothesized to be responsible for variety-specific differences in flower color. Furthermore, the genes encoding flavonoid 3-hydroxylase, anthocyanin synthase, and anthocyanin reductase were differentially expressed between flowers of different colors. Overall, we successfully identified key genes and metabolites involved in T. fournieri flower color formation. The data provided by the chromosome-scale genome assembly establish a basis for understanding the differentiation of this species and will facilitate future genetic studies and genomic-assisted breeding.

Chromosome-scale genome assembly of Poa trivialis and population genomics reveal widespread gene flow in a cool-season grass seed production system.

Poa trivialis (L.) is a cool-season grass species found in various environments worldwide. In addition to being a desired turfgrass species, it is a common weed of agricultural systems and natural areas. As a weed, it is an important contaminant of commercial cool-season grass seed lots, resulting in widespread gene flow facilitated by human activities and causing significant economic losses to farmers. To better understand and manage infestations, we assembled and annotated a haploid genome of P. trivialis and studied troublesome field populations from Oregon, the largest cool-season grass seed producing region in the United States. The genome assembly resulted in 1.35 Gb of DNA sequence distributed among seven chromosome-scale scaffolds, revealing a high content of transposable elements, conserved synteny with Poa annua, and a close relationship with other C3 grasses. A reduced-representation sequencing analysis of field populations revealed limited genetic diversity and suggested potential gene flow and human-assisted dispersal in the region. The genetic resources and insights into P. trivialis provided by this study will improve weed management strategies and enable the development of molecular detection tests for contaminated seed lots to limit seed-mediated gene flow. These resources should also be beneficial for turfgrass breeders seeking to improve desirable traits of commercial P. trivialis varieties and help to guide breeding efforts in other crops to enhance the resiliency of agricultural ecosystems under climate change. Significance Statement: The chromosome-scale assembly of Poa trivialis and population genomic analyses provide crucial insights into the gene flow of weedy populations in agricultural systems and contribute a valuable genomic resource for the plant science community.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species’ adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BackgroundThe genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite’s lifestyle.ResultsUsing the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument—the tissue located at the host-parasite physical interface—of intramolluscan life stages.ConclusionsWe propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

De Novo Genome Assembly at Chromosome-Scale of Hermetia illucens (Diptera Stratiomyidae) via PacBio and Omni-C Proximity Ligation Technology.

Hermetia illucens is a species of great interest for numerous industrial applications. A high-quality reference genome is already available for H. illucens. However, the worldwide maintenance of numerous captive populations of H. illucens, each with its own genotypic and phenotypic characteristics, made it of interest to perform a de novo genome assembly on one population of H. illucens to define a chromosome-scale genome assembly. By combining the PacBio and the Omni-C proximity ligation technologies, a new H. illucens chromosome-scale genome of 888.59 Mb, with a scaffold N50 value of 162.19 Mb, was assembled. The final chromosome-scale assembly obtained a BUSCO completeness of 89.1%. By exploiting the Omni-C proximity ligation technology, topologically associated domains and other topological features that play a key role in the regulation of gene expression were identified. Further, 65.62% of genomic sequences were masked as repeated sequences, and 32,516 genes were annotated using the MAKER pipeline. The H. illucens Lsp-2 genes that were annotated were further characterized, and the three-dimensional organization of the encoded proteins was predicted. A new chromosome-scale genome assembly of good quality for H. illucens was assembled, and the genomic annotation phase was initiated. The availability of this new chromosome-scale genome assembly enables the further characterization, both genotypically and phenotypically, of a species of interest for several biotechnological applications.

Chromosome-scale genome sequence of Nothopassalora personata (syn. Cercosporidium personatum), a devastating fungal pathogen of peanut

Nothopassalora personata, the causal agent of late leaf spot, is an economically important fungus that poses a significant threat to peanut production worldwide. Currently, studies on N. personata are focused on epidemiology, host resistance, early environmental detection, and molecular studies based on a limited number of genes for classification. A high-quality whole genome sequence can be a valuable tool for studying complex species and identifying molecular mechanisms in pathogenicity and adaptation. This study reports a high-quality, chromosome-scale genome assembly for N. personata with a total length of 54.6 Mb, an N50 of 2.54 Mb, a longest contig length of 4.47 Mb and a BUSCO genome completeness score of 99.5% with mostly single copy genes. Genome annotation revealed 11,295 protein-encoding genes and 12,121 transcripts in which 411 of the 995 predicted proteins were putative effectors. Genome comparison with the closely related fungi Dothistroma septosporum NZE10 and Fulvia fulva (syn. Passalora fulva, Cladosporium fulvum) isolate Race 5 revealed 13 chromosomes. This genome sequence resource can be used as a reference for aligning sequences of multiple isolates to provide relevant information on the diversity and population genetic structure of N. personata to further inform peanut resistance breeding programs and devise the best disease management strategies.

The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes

BackgroundUnderstanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood.MethodsIn this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus.ResultsWe then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes.ConclusionThe improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

A chromosome-scale assembly for 'd'Anjou' pear.

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.