Complete Genome Assembly Research Articles

The complete mitogenome of an unidentified Oikopleura species

Appendicularians are planktonic tunicates abundant all over the world. Currently, only two complete annotated mitochondrial genome assemblies are available for appendicularians, both for cryptic species of Oikopleura dioica. This underrepresentation of available appendicularian mitochondrial genomes limits environmental DNA sequencing (eDNA) studies that rely on mitochondrial markers as a taxonomic barcode. We report the complete mitochondrial genome assembly and annotation of an unknown appendicularian species isolated from the Amami Oshima island, Kagoshima prefecture, Japan, that has significant sequence difference with other currently available assemblies and will serve as a useful resource for ecological studies and further mitochondrial studies of appendicularians.

Genomic and Untargeted Metabolomic Analysis of Secondary Metabolites in the Streptomyces griseoaurantiacus Strain MH191 Shows Media-Based Dependency for the Production of Bioactive Compounds with Potential Antifungal Activity.

Streptomyces species can form beneficial relationships with hosts as endophytes, including the phytopathogen-inhibiting strain, Streptomyces griseoaurantiacusMH191, isolated from wheat plants. Using genomic characterization and untargeted metabolomics, we explored the capacity of strain MH191 to inhibit a range of fungal phytopathogens through the production of secondary metabolites. Complete genome assembly of strain MH191 predicted 24 biosynthetic gene clusters. Secondary metabolite production was assessed following culture on six different media, with the detection of 205 putative compounds. Members of the manumycin family, undecylprodigiosin, and desferrioxamine were identified as the predominant metabolites. Antifungal activity was validated for undecylprodigiosin and manumycin. These compounds were produced from different BGCs, which showed similarity to asukamycin, undecylprodigiosin, and FW0622 gene clusters, respectively. The growth of strain MH191 on different media illustrated the metabolic regulation of these gene clusters and the strain's extended chemical potential, with the asukamycin gene cluster alone, producing a variety of antifungal metabolites. The study highlights the extended chemical capability of strain MH191, which could be exploited as a biological control agent for designing future crop protection solutions.

GCI: a continuity inspector for complete genome assembly.

Recent advances in long-read sequencing technologies have significantly facilitated the production of high-quality genome assembly. The telomere-to-telomere (T2T) gapless assembly has become the new golden standard of genome assembly efforts. Several recent efforts have claimed to produce T2T-level reference genomes. However, a universal standard is still missing to qualify a genome assembly to be at T2T standard. Traditional genome assembly assessment metrics (N50 and its derivatives) have no capacity in differentiating between nearly T2T assembly and the truly T2T assembly in continuity either globally or locally. Additionally, these metrics are independent of raw reads, making them inflated easily by artificial operations. Therefore, a gaplessness evaluation tool at single-nucleotide resolution to reflect true completeness is urgently needed in the era of complete genomes. Here, we present a tool called Genome Continuity Inspector (GCI), designed to assess genome assembly continuity at single-base resolution, and evaluate how close an assembly is to the T2T level. GCI utilizes multiple aligners to map long reads from various sequencing platforms back to the assembly. By incorporating curated mapping coverage of high-confidence read alignments, GCI identifies potential assembly issues. Meanwhile, it provides GCI scores that quantify overall assembly continuity on the whole genome or chromosome scales. The open-source GCI code is freely available on Github (https://github.com/yeeus/GCI) under the MIT license.

Genome-Wide Identification of MKK Gene Family and Response to Hormone and Abiotic Stress in Rice.

Mitogen-activated protein kinase (MAPK/MPK) cascades are pivotal and highly conserved signaling modules widely distributed in eukaryotes; they play essential roles in plant growth and development, as well as biotic and abiotic stress responses. With the development of sequencing technology, the complete genome assembly of rice without gaps, T2T (Telomere-to-Telomere)-NIP (version AGIS-1.0), has recently been released. In this study, we used bioinformatic approaches to identify and analyze the rice MPK kinases (MKKs) based on the complete genome. A total of seven OsMKKs were identified, and their physical and chemical properties, chromosome localization, gene structure, subcellular localization, phylogeny, family evolution, and cis-acting elements were evaluated. OsMKKs can be divided into four subgroups based on phylogenetic relationships, and the family members located in the same evolutionary branch have relatively similar gene structures and conserved domains. Quantitative real-time PCR (qRT-PCR) revealed that all OsMKKs were highly expressed in rice seedling leaves. The expression levels of all OsMKKs were more or less altered under exogenous hormone and abiotic stress treatments, with OsMKK1, OsMKK6, and OsMKK3 being induced under almost all treatments, while the expression of OsMKK4 and OsMKK10-2 was repressed under salt and drought treatments and IAA treatment, respectively. In this study, we also summarized the recent progress in rice MPK cascades, highlighted their diverse functions, and outlined the potential MPK signaling network, facilitating further studies on OsMKK genes and rice MPK cascades.

Assembly of Mitochondrial Genomes Using Nanopore Long-Read Technology in Three Sea Chubs (Teleostei: Kyphosidae).

Complete mitochondrial genomes have become markers of choice to explore phylogenetic relationships at multiple taxonomic levels and they are often assembled using whole genome short-read sequencing. Herein, using three species of sea chubs as an example, we explored the accuracy of mitochondrial chromosomes assembled using Oxford Nanopore Technology (ONT) Kit 14 R10.4.1 long reads at different sequencing depths (high, low and very low or genome skimming) by comparing them to 'gold' standard reference mitochondrial genomes assembled using Illumina NovaSeq short reads. In two species of sea chubs, Girella nigricans and Kyphosus azureus, ONT long-read assembled mitochondrial genomes at high sequencing depths (> 25× whole [nuclear] genome) were identical to their respective short-read assembled mitochondrial genomes. Not a single 'homopolymer insertion', 'homopolymer deletion', 'simple substitution', 'single insertion', 'short insertion', 'single deletion' or 'short deletion' were detected in the long-read assembled mitochondrial genomes after aligning each one of them to their short-read counterparts. In turn, in a third species, Medialuna californiensis, a 25× sequencing depth long-read assembled mitochondrial genome was 14 nucleotides longer than its short-read counterpart. The difference in total length between the latter two assemblies was due to the presence of a short motif 14 bp long that was repeated (twice) in the long read but not in the short-read assembly. Read subsampling at a sequencing depth of 1× resulted in the assembly of partial or complete mitochondrial genomes with numerous errors, including, among others, simple indels, and indels at homopolymer regions. At 3× and 5× subsampling, genomes were identical (perfect) or almost identical (quasiperfect, 99.5% over 16,500 bp) to their respective Illumina assemblies. The newly assembled mitochondrial genomes exhibit identical gene composition and organisation compared with cofamilial species and a phylomitogenomic analysis based on translated protein-coding genes suggested that the family Kyphosidae is not monophyletic. The same analysis detected possible cases of misidentification of mitochondrial genomes deposited in GenBank. This study demonstrates that perfect (complete and fully accurate) or quasiperfect (complete but with a single or a very few errors) mitochondrial genomes can be assembled at high (> 25×) and low (3-5×) but not very low (1×, genome skimming) sequencing depths using ONT long reads and the latest ONT chemistries (Kit 14 and R10.4.1 flowcells with SUP basecalling). The newly assembled and annotated mitochondrial genomes can be used as a reference in environmental DNA studies focusing on bioprospecting and biomonitoring of these and other coastal species experiencing environmental insult. Given the small size of the sequencing device and low cost, we argue that ONT technology has the potential to improve access to high-throughput sequencing technologies in low- and moderate-income countries.

Complete genome assembly of Candida auris representative strains of three geographical clades.

The complete genome assembly of Candida auris strains B11103, B11221, and B11244 is reported in this manuscript. These strains represent the three geographical clades, namely, South Asian (Clade I), South African (Clade III), and South American (Clade IV).

Complete analysis of G-quadruplex forming sequences in the gapless assembly of human chromosome Y

Recent advancements have finally delivered a complete human genome assembly, including the elusive Y chromosome. This accomplishment closes a significant knowledge gap. Prior efforts were hampered by challenges in sequencing repetitive DNA structures such as direct and inverted repeats. We used the G4Hunter algorithm to analyze the presence of G-quadruplex forming sequences (G4s) within the current human reference genome (GRCh38) and the new telomere-to-telomere (T2T) Y chromosome assemblies. This analysis served a dual purpose: identifying the location of potential G4s within the genomes and exploring their association with functionally annotated sequences. Compared to GRCh38, the T2T assembly exhibited a significantly higher prevalence of G-quadruplex forming sequences. Notably, these repeats were abundantly located around precursor RNA, exons, genes, and within protein binding sites. This remarkable co-occurrence of G4-forming sequences with these critical regulatory regions suggests their role in fundamental DNA regulation processes. Our findings indicate that the current human reference genome significantly underestimated the number of G4s, potentially overlooking their functional importance.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BackgroundThe Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.ResultsThe 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.ConclusionsOur results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Genome sequencing and mining expand the naturalproduct repertoire of Lysobacter.

Compounds produced by living organisms serve as an important source of inspiration for the development of pharmaceuticals. A potential source of new natural products are bacteria from a genus with species that are known to produce bioactive natural products, but are relatively understudied. Lysobacter is a genus of bacteria that have attracted attention as possible biocontrol agents and are known to produce antibiotic natural products. To further explore the biosynthetic potential of Lysobacter, we sequenced the genomes of two species and performed genome mining studies on those and publicly available genomes. In this study we produced draft genome sequences for Lysobacter firmicutimachus and Lysobacter yananisis. We additionally examined 113 publicly available Lysobacter genomes and found that biosynthetic potential of individual species ranges broadly, with species having between 1 and nearly 20 biosynthetic gene clusters. Filtering for more complete genome assemblies and 9 or more biosynthetic gene clusters, we performed genome mining on 24 Lysobacter genomes. Within these genomes we identified 21 unique nonribosomal peptide, 11 unique hybrid polyketide/nonribosomal peptide, 4 unique polyketide, and 27 unique lanthipeptide biosynthetic gene clusters that produce uncharacterized compounds. Additionally, we tentatively identified the biosynthetic gene cluster in L. rmicutimachus responsible for producing plusbacins, which has not been previously identified. This study demonstrated that Lysobacter have a large repertoire of natural products that remain to be characterized. Additionally, we found that some Lysobacter species are substantially more biosynthetically gifted than others and that strains of the same species of Lysobacter have similar biosynthetic capacities.

A genome assembly for the Chryxus Arctic (Oeneis chryxus), the highest butterfly in North America.

We describe a highly contiguous and complete diploid genome assembly for the Chryxus Arctic, Oeneis chryxus (E. Doubleday, [1849]), a butterfly species complex spanning much of northern and western North America. One subspecies, the Ivallda Arctic (O. c. ivallda), is endemic to California's Sierra Nevada and of particular biogeographic interest and conservation concern. Extreme alpine habitats occupied by this subspecies include the summit of Mt. Whitney, California, representing the highest elevation butterfly population in North America. The assembly presented here consists of two haplotypes, 738.92 and 770.85Mb in length, with contig N50 values of 10.49 and 10.13Mb, scaffold N50 values of 25.35 and 25.69Mb, scaffold L50 values of 13 and 14, and BUSCO completeness scores of 96.5 and 98.3%, respectively. More than 97% of the assembly is organized into 29 scaffolds, which likely represent whole chromosomes. This assembly is the first major genomic resource for Oeneis, providing a foundational reference for future genomic studies on the taxonomy, evolutionary history, and conservation of the genus. As part of the California Conservation Genomics Project, we will use this assembly in conjunction with short-read resequencing to resolve patterns of evolutionary differentiation, adaptive genomic variation, and gene flow among remaining O. c. ivallda populations. These data can and will be used to inform the subspecies' conservation as warming climatic conditions continue to lead to the loss and fragmentation of alpine habitats. We also provide genome assemblies for the O. chryxus mitochondrion and a Wolbachia endosymbiont.

Complete Annotated Genome Assembly of Flax Pathogen Colletotrichum lini.

Colletotrichum lini is a fungal pathogen of flax that can cause significant yield and quality losses. In this work, we obtained the first complete annotated genome assembly of the highly virulent C. lini strain #394-2. The nuclear genome consisted of ten core and two accessory chromosomes and had a length of 53.7 Mb. The mitochondrial genome was 39.1 kb. The assembly was obtained by the Canu-Racon ×2-Medaka-Polca algorithm using Oxford Nanopore Technologies and Illumina data. As a result of the annotation with the Illumina RNA-Seq data, 12,449 genes were identified. Potential signaling proteins were tested for effector functions and 550 effector proteins were predicted using EffectorP. The visualization of the effector protein localization revealed that the presence of effector proteins was associated with repeat-rich regions. A comparison of the genomic structure of C. lini with chromosome-level and complete assemblies of the genus Colletotrichum representatives revealed that the genomes of Colletotrichum species differed by the presence of chromosomal rearrangements. The obtained assembly expands the knowledge of the genomic structure of Colletotrichum species and provides the basis for further studies of C. lini, which will help to understand the virulence mechanisms and protect flax from anthracnose.

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development.

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

Utilizing Deep Neural Networks to Fill Gaps in Small Genomes.

With the widespread adoption of next-generation sequencing technologies, the speed and convenience of genome sequencing have significantly improved, and many biological genomes have been sequenced. However, during the assembly of small genomes, we still face a series of challenges, including repetitive fragments, inverted repeats, low sequencing coverage, and the limitations of sequencing technologies. These challenges lead to unknown gaps in small genomes, hindering complete genome assembly. Although there are many existing assembly software options, they do not fully utilize the potential of artificial intelligence technologies, resulting in limited improvement in gap filling. Here, we propose a novel method, DLGapCloser, based on deep learning, aimed at assisting traditional tools in further filling gaps in small genomes. Firstly, we created four datasets based on the original genomes of Saccharomyces cerevisiae, Schizosaccharomyces pombe, Neurospora crassa, and Micromonas pusilla. To further extract effective information from the gene sequences, we also added homologous genomes to enrich the datasets. Secondly, we proposed the DGCNet model, which effectively extracts features and learns context from sequences flanking gaps. Addressing issues with early pruning and high memory usage in the Beam Search algorithm, we developed a new prediction algorithm, Wave-Beam Search. This algorithm alternates between expansion and contraction phases, enhancing efficiency and accuracy. Experimental results showed that the Wave-Beam Search algorithm improved the gap-filling performance of assembly tools by 7.35%, 28.57%, 42.85%, and 8.33% on the original results. Finally, we established new gap-filling standards and created and implemented a novel evaluation method. Validation on the genomes of Saccharomyces cerevisiae, Schizosaccharomyces pombe, Neurospora crassa, and Micromonas pusilla showed that DLGapCloser increased the number of filled gaps by 8.05%, 15.3%, 1.4%, and 7% compared to traditional assembly tools.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P.setigerum (2n= 44), P.somniferum (2n= 22), P.rhoeas (2n= 14), and P.bracteatum (2n= 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P.somniferum and P.setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

North American crayfish harbour diverse members of the Nudiviridae

Three novel crayfish-infecting nudiviruses from crayfish in North America represent the first genomic confirmation of nudiviruses in crayfish: Faxonius propinquus nudivirus (FpNV), Faxonius rusticus nudivirus (FrNV), and Faxonius virilis nudivirus (FvNV). Histopathology and electron microscopy revealed nuclear infections, including nuclear hypertrophy in hepatopancreatic epithelial cells and the presence of membrane-bound bacilliform virions. Metagenomic sequencing resulted in complete circular genome assembly, and phylogenetic analyses (based on nudivirus core genes) placed these viruses within the unofficial Epsilonnudivirus genus. One of the nudiviruses was detected in the antennal gland of its host, and another is correlated with invasive crayfish decline in one infected lake ecosystem - suggesting a potential route for viral transmission through water, and possible population level impact. This study highlights the importance of genomic and ecological data in elucidating the diversity and evolutionary relationships of the Nudiviridae, while expanding their known diversity and range of host species.

Highly contiguous genome assembly of Drosophila prolongata-a model for evolution of sexual dimorphism and male-specific innovations.

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of Southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2 Mb) genome assembly of 220 Mb. The repetitive content of the genome is 24.6%, the plurality of which are long terminal repeats retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, genome-wide association studies (GWAS), and other genomic analyses.

A near complete genome assembly of the East Friesian sheep genome

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Insights into the evolution and fruit color change‐related genes of chromosome doubled sweet cherry from an updated complete T2T genome assembly

Insights into the evolution and fruit color change‐related genes of chromosome doubled sweet cherry from an updated complete T2T genome assembly

The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain.

Dermacoccus barathri is the first reported pathogen within the Dermacoccus genus to cause a catheter-related bloodstream infection, which occurred in 2015. In this study, the complete genome assembly of Dermacoccus barathri was constructed, and the complete genome of Dermacoccus barathri FBCC-B549 consists of a single chromosome (3,137,745 bp) without plasmids. The constructed genome of D. barathri was compared with those of two closely related species within the Dermacoccus genus. D. barathri exhibited a pattern similar to Dermacoccus abyssi in terms of gene clusters and synteny analysis. Contrary to previous studies, biosynthetic gene cluster (BGC) analysis for predicting secondary metabolites revealed the presence of the LAP biosynthesis pathway in the complete genome of D. barathri, predicting the potential synthesis of the secondary metabolite plantazolicin. Furthermore, an analysis to investigate the potential pathogenicity of D. barathri did not reveal any antibiotic resistance genes; however, nine virulence factors were identified in the Virulence Factor Database (VFDB). According to these matching results in the VFDB, despite identifying a few factors involved in biofilm formation, further research is required to determine the actual impact of D. barathri on pathogenicity. The complete genome of D. barathri is expected to serve as a valuable resource for future studies on D. barathri, which currently lack sufficient genomic sequence information.

Genome Sequence Comparisons between Small and Large Colony Phenotypes of Equine Clinical Isolates of Arcanobacterium hippocoleae.

Arcanobacterium hippocoleae is a Gram-positive fastidious bacterium and is occasionally isolated from the reproductive tract of apparently healthy mares (Equus caballus) or from mares with reproductive tract abnormalities. Apart from a few 16S rRNA gene-based GenBank sequences and one recent report on complete genome assembly, detailed genomic sequence and clinical experimental data are not available on the bacterium. Recently, we observed an unusual increase in the detection of the organism from samples associated with mare reproductive failures in Atlantic Canada. Two colony morphotypes (i.e., small, and large) were detected in culture media, which were identified as A. hippocoleae by MALDI-TOF mass spectrometry and 16S rRNA gene sequencing. Here, we report the whole genome sequencing and characterization of the morphotype variants. The genome length of the large phenotypes was between 2.42 and 2.43, and the small phenotype was 1.99 Mbs. The orthologous nucleotide identity between the large colony phenotypes was ~99%, and the large and small colony phenotypes was between 77.86 and 78.52%, which may warrant the classification of the two morphotypes into different species. Phylogenetic analysis based on 16S rRNA genes or concatenated housekeeping genes grouped the small and large colony variants into two different genotypic clusters. The UvrA protein, which is part of the nucleotide excision repair (NER) system, and 3-isopropoylmalate dehydratase small subunit protein expressed by the leuD gene were identified as potential virulence factors in the large and small colony morphotypes, respectively. However, detailed functional studies will be required to determine the exact roles of these and other identified hypothetical proteins in the cellular metabolism and potential pathogenicity of A. hippocoleae in mares.