Resource For Comparative Genomics Research Articles

Telomere-to-telomere Genome Assembly of two representative Asian and European pear cultivars

As the third most important temperate fruit, Pear (Pyrus spp.) exhibits a remarkable genetic diversity and is classified into two mainly categories known as Asian pear and European pear. Although several pear genomes are available, most of the released versions are fragmented and not chromosome-level high-quality. In this study, we report two high-quality genomes for Pyrus bretschneideri Rhed. cv. ‘Danshansuli’ (DS) and Pyrus communis L. cv. ‘Conference’ (KFL), which represent the predominant Asian and European cultivars, respectively, with nearly telomere-to-telomere (T2T) gap-free level. The finally assembled genome sizes for DS and KFL were 510.98 Mb and 510.71 Mb, respectively, with Contig N50 of 29.47 Mb and 30.47 Mb, where each chromosome was represented by a single contig. The DS and KFL genomes yielded a total of 46,394 and 44,702 protein-coding genes, respectively. Among these genes, the functional annotation accounted for 96.47% and 96.46% in the DS and KFL genomes. The two novels nearly T2T genomic information offers an invaluable resource for comparative genomics, genetic diversity analysis, molecular breeding strategies, and functional exploration.

Genome resequencing and genome-wide polymorphisms in mosquito vectors Aedes aegypti and Aedes albopictus from south India

Aedes aegypti and Aedes albopictus mosquitoes spread major vector-borne viral diseases in tropical and sub-tropical regions of the globe. In this study, we sequenced the genome of Indian Ae. aegypti and Ae. albopictus and mapped to their reference genomes. Comparative genomics were performed between our strain and the reference strains. A total of 14,416,484 single nucleotide polymorphisms (SNPs) and 156,487 insertions and deletions (InDels) were found in Ae. aegypti, and 28,940,433 SNPs and 188,987 InDels in Ae. albopictus. Particular emphasis was given to gene families involved in mosquito digestion, development, and innate immunity, which could be putative candidates for vector control. Serine protease cascades and their inhibitors called serpins, play a central role in these processes. We extracted high-impact variants in genes associated with serine proteases and serpins. This study reports for the first time a high coverage genome sequence data of an Indian Ae. albopictus mosquito. The results from this study will provide insights into Indian Aedes specific polymorphisms and the evolution of immune related genes in mosquitoes, which can serve as a resource for future comparative genomics and those pursuing the development of targeted biopesticides for effective mosquito control strategies.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life.

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Population Genomics analysis of Leptosphaeria biglobosa Associated with Brassica napus in China Reveals That Geographical Distribution Influences Its Genetic Polymorphism.

Blackleg disease, a major threat to Brassica crops worldwide, is primarily caused by the pathogen Leptosphaeria biglobosa. Investigating the genetic variation of L. biglobosa is crucial for managing and preventing the disease in Brassica napus. To date, there is scarce genomic variation information available for populations of L. biglobosa in China. In this study, 73 L. biglobosa strains of canola stalks were collected from 12 provinces in China and subjected to re-sequencing. The 73 assemblies averaged 1340 contigs, 72,123 bp N50, and 30.17 Mb in size. In total, 9409 core orthogroups and 867 accessory orthogroups were identified. A total of 727,724 mutation loci were identified, including 695,230 SNPs and 32,494 indels. Principal component analysis (PCA) and population structure analysis showed that these strains could be divided into seven subgroups. The strains in most provinces were clustered into a single subgroup, suggesting a strong influence of the geographic environment on strain variation. The average nucleotide diversity (θπ) of all strains was 1.03 × 10-3, indicating important genetic diversity among strains from different regions of China. This study provides valuable resources for future comparative genomics, gives new insights into the population evolution of L. biglobosa, and supports the development of strategies for managing blackleg disease in canola.

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.

A chromosome-level genome assembly of the disco clam, Ctenoides ales.

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

A chromosome-level genome assembly of East Asia endemic minnow Zacco platypus

Zacco platypus is an endemic colorful freshwater minnow that is intensively distributed in East Asia. In this study, two adult female individuals collected from Haihe River basin were used for karyotypic study and genome sequencing, respectively. The karyotype formula of Z. platypus is 2N = 48 = 18 M + 24SM/ST + 6 T. We used PacBio long-read sequencing and Hi-C technology to assemble a chromosome-level genome of Z. platypus. As a result, an 814.87 Mb genome was assembled with the PacBio long reads. Subsequently, 98.64% assembled sequences were anchored into 24 chromosomes based on the Hi-C data. The chromosome-level assembly contained 54 scaffolds with a N50 length of 32.32 Mb. Repeat elements accounted for 52.35% in genome, and 24,779 protein-coding genes were predicted, with 92.11% were functionally annotated with the public databases. BUSCO analysis yielded a completeness score of 96.5%. This high-quality genome assembly provides valuable resources for future functional genomic research, comparative genomics, and evolutionary studies of genus Zacco.

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

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level. Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina. Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers. The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754. Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

A Telomere-to-Telomere Genome Assembly Resource of Bipolaris sorokiniana, the Fungal Pathogen Causing Spot Blotch and Common Root Rot Diseases in Barley and Wheat

Bipolaris sorokiniana ( Cochliobolus sativus) is a fungal pathogen that causes spot blotch, common root rot, and kernel blight in barley and wheat. Four pathotypes (0, 1, 2, and 7) of the fungus were previously identified based on their virulence on three barley differential lines (Bowman, ND 5883, and NDB 112). Although several genome assemblies have been reported for B. sorokiniana, a telomere-to-telomere genome assembly is still lacking for this fungus. In this study, we assembled the genome of a pathotype 0 isolate (ND93-1) of B. sorokiniana using both PacBio HiFi reads and ultralong Oxford Nanopore Technology reads corrected with Illumina paired-end 100-bp short reads. The combined genome assembly of ND93-1 has an estimated size of 35.7 Mb and consists of 16 scaffolds, each having two telomeres. A total of 11,564 protein-coding genes were predicted, including 1,633 genes encoding for secretory proteins and 473 genes for effectors. This telomere-to-telomere genome assembly provides an important resource for comparative genomics and understanding molecular biology of B. sorokiniana and related Bipolaris species. [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 .

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

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

Genome Assembly and Microsatellite Marker Development Using Illumina and PacBio Sequencing in the Carex pumila (Cyperaceae) from Korea.

This study is the first to report the characterization of Carex pumila genomic information. Assembly of the genome generated a draft of C. pumila based on PacBio Sequel II and Illumina paired-end sequencing, which was assembled from 2941 contigs with an estimated genome size of 0.346 Gb. The estimate of repeats in the genome was 31.0%, and heterozygosity ranged from 0.426 to 0.441%. The integrity evaluation of the assembly revealed 1481 complete benchmarked universal single-copy orthologs (BUSCO) (91.76%), indicating the high quality of the draft assembly. A total of 23,402 protein-coding genes were successfully predicted and annotated in the protein database. UpsetR plots showed that 7481 orthogroups were shared by all species. The phylogenetic tree showed that C. pumila is a close but distant relative of Ananas comosus. C. pumila had greater contraction (3154) than expansion (392). Among the extended gene families, aquaporins have been found to be enriched. Primers for microsatellite markers determined 30 polymorphic markers out of 100. The average number of alleles amplified by these 30 polymorphic markers was 4 to 12, with an average polymorphism information content (PIC) value of 0.660. In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of C. pumila.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

A chromosome-level genome assembly of the darkbarbel catfish Pelteobagrus vachelli

The darkbarbel catfish (Pelteobagrus vachelli), an economically important aquaculture species in China, is extensively employed in hybrid yellow catfish production due to its superior growth rate. However, information on its genome has been limited, constraining further genetic studies and breeding programs. Leveraging the power of PacBio long-read sequencing and Hi-C technologies, we present a high-quality, chromosome-level genome assembly for the darkbarbel catfish. The resulting assembly spans 692.10 Mb, with an impressive 99.9% distribution over 26 chromosomes. The contig N50 and scaffold N50 are 13.30 Mb and 27.55 Mb, respectively. The genome is predicted to contain 22,109 protein-coding genes, with 96.1% having functional annotations. Repeat elements account for approximately 35.79% of the genomic landscape. The completeness of darkbarbel catfish genome assembly is highlighted by a BUSCO score of 99.07%. This high-quality genome assembly provides a critical resource for future hybrid catfish breeding, comparative genomics, and evolutionary studies in catfish and other related species.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis.

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

A chromosome-level genome assembly of Stenchaetothrips biformis and comparative genomic analysis highlights distinct host adaptations among thrips

Insects have a limited host range due to genomic adaptation. Thysanoptera, commonly known as thrips, occupies distinct feeding habitats, but there is a lack of comparative genomic analyses and limited genomic resources available. In this study, the chromosome-level genome of Stenchaetothrips biformis, an oligophagous pest of rice, is assembled using multiple sequencing technologies, including PacBio, Illumina short-reads, and Hi-C technology. A 338.86 Mb genome is obtained, consisting of 1269 contigs with a contig N50 size of 381 kb and a scaffold N50 size of 18.21 Mb. Thereafter, 17,167 protein-coding genes and 36.25% repetitive elements are annotated. Comparative genomic analyses with two other polyphagous thrips, revealing contracted chemosensory-related and expanded stress response and detoxification gene families in S. biformis, potentially facilitating rice adaptation. In the polyphagous thrips species Frankliniella occidentalis and Thrips palmi, expanded gene families are enriched in metabolism of aromatic and anthocyanin-containing compounds, immunity against viruses, and detoxification enzymes. These expansion gene families play crucial roles not only in adapting to hosts but also in development of pesticide resistance, as evidenced by transcriptome results after insecticides treatment. This study provides a chromosome-level genome assembly and lays the foundation for further studies on thrips evolution and pest management.

Efficient and highly continuous chromosome-level genome assembly of the first chameleon genome.

Most amniote genomes are diploid, moderate in size (approximately 1-6 Gbp), and contain a large proportion of repetitive sequences. The development of next-generation sequencing technology, especially the emergence of high-fidelity (HiFi) long-read data, has made it feasible to resolve high-quality genome assembly for non-model species efficiently. However, reference genomes for squamate reptiles has lagged behind other amniote lineages. Here we de novo assembled the first genome from the Chameleonidae family, the panther chameleon (Furcifer pardalis). We obtained telomere-to-telomere contigs using only HiFi data, reaching a contig N50 of 158.72 Mbp. The final chromosome-level assembly is 1.61 Gbp in size and 100% of primary contigs were placed to pseudochromosomes using Hi-C interaction data. We also found that sequencing depth > 30 folds can ensure both the integrity and accuracy of the genome, while insufficient depth led to false increase in genome size and proportion of duplicated genes. We provide a high-quality reference genome valuable for evolutionary and ecological studies in chameleons as well as providing comparative genomic resources for squamate reptiles.

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Evolutionary constraint and innovation across hundreds of placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

A near complete genome assembly of chia assists in identification of key fatty acid desaturases in developing seeds.

Chia is an annual crop whose seeds have the highest content of α-linolenic acid (ALA) of any plant known to date. We generated a high-quality assembly of the chia genome using circular consensus sequencing (CCS) of PacBio. The assembled six chromosomes are composed of 21 contigs and have a total length of 361.7 Mb. Genome annotation revealed a 53.5% repeat content and 35,850 protein-coding genes. Chia shared a common ancestor with Salvia splendens ~6.1 million years ago. Utilizing the reference genome and two transcriptome datasets, we identified candidate fatty acid desaturases responsible for ALA biosynthesis during chia seed development. Because the seed of S. splendens contains significantly lower proportion of ALA but similar total contents of unsaturated fatty acids, we suggest that strong expression of two ShFAD3 genes are critical for the high ALA content of chia seeds. This genome assembly will serve as a valuable resource for breeding, comparative genomics, and functional genomics studies of chia.

Genome Assembly and Comparative Analysis of the Egg Parasitoid Wasp Trichogramma dendrolimi Shed Light on the Composition and Evolution of Olfactory Receptors and Venoms

Trichogramma dendrolimi is one of the most successfully industrialized Trichogramma species used to control agricultural and forestry pests in China. However, the molecular mechanisms underlying its host recognition and parasitism remain largely unknown, partially due to the limited genome information of this parasitoid wasp. Here, we present a high-quality de novo assembly of T. dendrolimi through a combination of Illumina and PacBio sequencing technologies. The final assembly had a length of 215.2 Mb and contains 316 scaffolds with a scaffold N50 size of 1.41 Mb. Repetitive sequences with a length of 63.4 Mb and 12,785 protein-coding genes were identified. Significantly expanded gene families were identified to be involved in the development and regulatory processes, while remarkably contracted gene families were involved in the transport processes in T. dendrolimi. The olfactory and venom-associated genes were identified in T. dendrolimi and 24 other hymenopteran species, using uniform methods combining BLAST and HMM profiling. The identified venom genes of T. dendrolimi were enriched in antioxidant activity, tricarboxylic acid cycle, response to oxidative stress and cell redox homeostasis. Our study provides an important resource for comparative genomics and functional studies to interpret the molecular mechanisms underlying host recognition and parasitism of Trichogramma species.