Whole-genome Alignment Research Articles

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder (PGGB), a pipeline for constructing pangenome graphs without bias or exclusion. PGGB uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

High-quality genome of a novel Thermosynechococcaceae species from Namibia and characterization of its protein expression patterns at elevated temperatures.

Thermophilic cyanobacteria thrive in extreme environments, making their thermoresistant enzymes valuable for industrial applications. Common habitats include hot springs, which act as evolutionary accelerators for speciation due to geographical isolation. The family Thermosynechococcaceae comprises thermophilic cyanobacteria known for their ability to thrive in high-temperature environments. These bacteria are notable for their photosynthetic capabilities, significantly contributing to primary production in extreme habitats. Members of Thermosynechococcaceae exhibit unique adaptations that allow them to perform photosynthesis efficiently at elevated temperatures, making them subjects of interest for studies on microbial ecology, evolution, and potential biotechnological applications. In this study, the genome of a thermophilic cyanobacterium, isolated from a hot spring near Okahandja in Namibia, was sequenced using a PacBio Sequel IIe long-read platform. Cultivations were performed at elevated temperatures of 40, 50, and 55°C, followed by proteome analyses based on the annotated genome. Phylogenetic investigations, informed by the 16S rRNA gene and aligned nucleotide identity (ANI), suggest that the novel cyanobacterium is a member of the family Thermosynechococcaceae. Furthermore, the new species was assigned to a separate branch, potentially representing a novel genus. Whole-genome alignments supported this finding, revealing few conserved regions and multiple genetic rearrangement events. Additionally, 129 proteins were identified as differentially expressed in a temperature-dependent manner. The results of this study broaden our understanding of cyanobacterial adaptation to extreme environments, providing a novel high-quality genome of Thermosynechococcaceae cyanobacterium sp. Okahandja and several promising candidate proteins for expression and characterization studies.

Uncovering Candidate Pathogenicity Genes in Erwinia pyrifoliae YKB12327 via Tn5-insertion Mutagenesis

<i>Erwinia pyrifoliae</i> is a gram-negative bacterial pathogen that commonly causes black shoot blight in pear and apple tree. Although the pathogenicity of this bacterial species is very similar to <i>E. amylovora</i>, there is no specific explanation of its pathogenic genes and mechanisms. In this study, our investigation into <i>E. pyrifoliae</i> pathogenicity involved generating seven YKB12327 mutant strains using Tn5 transposon mutagenesis. Observations revealed weakened growth rate and loss of pathogenicity in these mutants. Whole-genome sequencing and alignment analysis identified transposon insertions within the coding sequences of five strains and in the intergenic region of two strains. Annotation analysis elucidated genes directly or indirectly associated with pathogenicity. Notably, mutant strain MT16 displayed a transposon insertion mutation in the cyclic-di-GMP phosphodiesterase (<i>pdeF</i>) gene, a key player in bacterial signaling, governing microbial behavior and adaptation to environmental changes. Our findings provide insights into the genetic regulation of <i>E. pyrifoliae</i> pathogenicity, suggesting potential avenues for further research aimed at understanding and controlling this bacterial pathogen by targeting <i>pdeF</i> to mitigate apple black shoot blight disease.

Development of Reliable Male-Specific Molecular Markers for Genetic Sex Identification in Sea Cucumber Apostichopus japonicus.

Apostichopus japonicus is an important marine aquaculture species in China, with high nutritional and economic value. In A. japonicus, there is no obvious sexual dimorphism in external appearance, and sex differentiation primarily relies on the observation of mature gonads after dissection, which leads to difficulties in sex identification. The confusion in sex identification greatly reduces breeding efficiency in the sea cucumber industry. Therefore, developing a reliable sex-specific marker is crucial. In this study, we identified 586 male-specific sequences through whole-genome assembly and sequence alignment, but did not identify any female-specific sequences, inferring an XY-type sex determination system in sea cucumbers. We developed a set of male-specific molecular markers to establish an accurate, stable, and widely adaptable genetic sex identification technique for A. japonicus. The male-specific molecular markers were validated with 100% accuracy in sea cucumber populations from six different geographical regions in China. In conclusion, this study provides further evidence for the XY-type sex determination system in A. japonicus and establishes an effective genetic sex identification method for multi-geographic populations, which benefits future study on reproductive biology and has significant implications in sea cucumber aquaculture industry.

Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies.

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting (ILS) and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and single nucleotide polymorphism-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that ILS, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Assembly and Genome Annotation of Different Strains of Apple Fruit Moth Virus (Cydia pomonella granulovirus).

Cydia pomonella granulovirus is a natural pathogen for Cydia pomonella that is used as a biocontrol agent of insect populations. The study of granulovirus virulence is of particular interest since the development of resistance in natural populations of C. pomonella has been observed during the long-term use of the Mexican isolate CpGV. In our study, we present the genomes of 18 CpGV strains endemic to southern Russia and from Kazakhstan, as well as a strain included in the commercial preparation "Madex Twin", which were sequenced and analyzed. We performed comparative genomic analysis using several tools. From comparisons at the level of genes and protein products that are involved in the infection process of virosis, synonymous and missense substitution variants have been identified. The average nucleotide identity has demonstrated a high similarity with other granulovirus genomes of different geographic origins. Whole-genome alignment of the 18 genomes relative to the reference revealed regions of low similarity. Analysis of gene repertoire variation has shown that BZR GV 4, BZR GV 6, and BZR GV L-7 strains have been the closest in gene content to the commercial "Madex Twin" strain. We have confirmed two deletions using read depth coverage data in regions lacking genes shown by homology analysis for granuloviruses BZR GV L-4 and BZR GV L-6; however, they are not related to the known genes causing viral pathogenicity. Thus, we have isolated novel CpGV strains and analyzed their potential as strains producing highly effective bioinsecticides against C. pomonella.

Unveiling the complete chloroplast genome of tribulus macropterus var. Arabicus (hosni) al-hemaid & j. Thomas: genome structure, comparative analysis and phylogeny

The present investigation unveiled the first complete chloroplast (cp) genome of Tribulus macropterus var. arabicus (Hosni) Al-Hemaid & J. Thomas (Zygophyllaceae), a medicinal herb, indigenous to Saudi Arabia. The cp genome, comprising a length of 158,179 bp and a GC content of 35.8%, exhibited the typical circular quadripartite arrangement of flowering plants, including two inverted repeat regions (25842 bp) separated by a large single-copy (88873 bp) and a small single-copy region (17622 bp). Genome annotation unraveled 132 genes, comprising of 87 protein-coding genes, 37 tRNAs and eight rRNAs. A comparative plastomics approach demonstrated a very similar genome structure and gene organization in closely related taxa. Whole-genome alignment indicated that the inverted repeat regions exhibited greater conservation when compared to the single copy zones. Repeat analysis of the cp genome identified 80 simple sequence repeats, with the majority (64) being mononucleotides. Among the longer repeats, forward repeats were frequent (20) followed by palindromic and reverse repeats. The nucleotide diversity endeavor identified some hypervariable sites (rpl22, cemA, trnL-UAG) in the small and large single copies which would offer opportunities to design molecular markers for potential application in taxonomic identification. Phylogenetic analysis with rbcL barcode elucidated the distinct position of T. macropterus var. arabicus compared to T. macropterus within Zygophyllaceae and further validated the assembly. The findings of this investigation provide significant insights into biological fields, particularly enhancing the current understanding of the genetic and evolutionary variations within Zygophyllaceae. Bangladesh J. Plant Taxon. 31(1): 1-14, 2024 (June)

Whole-Genome Alignment: Methods, Challenges, and Future Directions

Whole-genome alignment (WGA) is a critical process in comparative genomics, facilitating the detection of genetic variants and aiding our understanding of evolution. This paper offers a detailed overview and categorization of WGA techniques, encompassing suffix tree-based, hash-based, anchor-based, and graph-based methods. It elaborates on the algorithmic properties of these tools, focusing on performance and methodological aspects. This paper underscores the latest progress in WGA, emphasizing the increasing capacity to manage the growing intricacy and volume of genomic data. However, the field still grapples with computational and biological hurdles affecting the precision and speed of WGA. We explore these challenges and potential future solutions. This paper aims to provide a comprehensive resource for researchers, deepening our understanding of WGA tools and their applications, constraints, and prospects.

Comparative Probiogenomics Analysis of Limosilactobacillus fermentum 3872.

Limosilactobacillus fermentum is an important member of the lactic acid bacteria group and holds immense potential for probiotic properties in human health and relevant industries. In this study, a comparative probiogenomic approach was applied to analyze the genome sequence of L. fermentum 3872, which was extracted from a commercially available yogurt sample, along with 20 different publicly available strains. Results indicate that the genome size of the characterized L. fermentum 3892 strain is 2,057,839bp, with a single- and circular-type chromosome possessing a G + C content of 51.69%. The genome of L. fermentum 3892 strain comprises a total of 2120 open reading frames (ORFs), two genes encoding rRNAs, and 53 genes encoding tRNAs. Upon comparative probiogenomic analysis, two plasmid sequences were detected among the study strains, including one for the L. fermentum 3872 genome, which was found between position 1,288,203 and 1,289,237 with an identity of 80.98. The whole-genome alignment revealed 2223 identical sites and a pairwise identity of 98.9%, indicating a significant difference of 1.1% among genome strains. Comparison of amino acid encoding genes among strains included in this study suggests that the strain 3872 exhibited the highest degree of amino acids present, including glutamine, glutamate, aspartate, asparagine, lysine, threonine, methionine, and cysteine. The comparative antibiotic resistome profiling revealed that strain 3872 exhibited a high resistant capacity only to ciprofloxacin antibiotics as compared to other strains. This study provides a genomic-based evaluation approach for comparative probiotic strain analysis in commercial foods and their significance to human health.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments.

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BackgroundFungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.ResultsHere, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.ConclusionsOur results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

The origin and structural evolution of de novo genes in Drosophila

Recent studies reveal that de novo gene origination from previously non-genic sequences is a common mechanism for gene innovation. These young genes provide an opportunity to study the structural and functional origins of proteins. Here, we combine high-quality base-level whole-genome alignments and computational structural modeling to study the origination, evolution, and protein structures of lineage-specific de novo genes. We identify 555 de novo gene candidates in D. melanogaster that originated within the Drosophilinae lineage. Sequence composition, evolutionary rates, and expression patterns indicate possible gradual functional or adaptive shifts with their gene ages. Surprisingly, we find little overall protein structural changes in candidates from the Drosophilinae lineage. We identify several candidates with potentially well-folded protein structures. Ancestral sequence reconstruction analysis reveals that most potentially well-folded candidates are often born well-folded. Single-cell RNA-seq analysis in testis shows that although most de novo gene candidates are enriched in spermatocytes, several young candidates are biased towards the early spermatogenesis stage, indicating potentially important but less emphasized roles of early germline cells in the de novo gene origination in testis. This study provides a systematic overview of the origin, evolution, and protein structural changes of Drosophilinae-specific de novo genes.

Graphylo: A deep learning approach for predicting regulatory DNA and RNA sites from whole-genome multiple alignments

This study focuses on enhancing the prediction of regulatory functional sites in DNA and RNA sequences, a crucial aspect of gene regulation. Current methods, such as motif overrepresentation and machine learning, often lack specificity. To address this issue, the study leverages evolutionary information and introduces Graphylo, a deep-learning approach for predicting transcription factor binding sites in the human genome. Graphylo combines Convolutional Neural Networks for DNA sequences with Graph Convolutional Networks on phylogenetic trees, using information from placental mammals' genomes and evolutionary history. The research demonstrates that Graphylo consistently outperforms both single-species deep learning techniques and methods that incorporate inter-species conservation scores on a wide range of datasets. It achieves this by utilizing a species-based attention model for evolutionary insights and an integrated gradient approach for nucleotide-level model interpretability. This innovative approach offers a promising avenue for improving the accuracy of regulatory site prediction in genomics.

A super-pangenome of the North American wild grape species

BackgroundCapturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.ResultsHere we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce’s disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.ConclusionsThis study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

Genomic investigation of the emergence of vanD vancomycin-resistant Enterococcus faecium.

Vancomycin-resistant Enterococcus (VRE) is an increasingly identified cause of human disease, with most infections resulting from the vanA and vanB genotypes; less is known about other clinically relevant genotypes. Here we report a genomic exploration of a vanD VRE faecium (VREfm), which arose de novo during a single infectious episode. The genomes of the vancomycin-susceptible E. faecium (VSEfm) recipient and resulting VREfm were subjected to long-read sequencing and closed, with whole-genome alignments, cross-mapping and orthologue clustering used to identify genomic variation. Three key differences were identified. (i) The VREfm chromosome gained a 142.6 kb integrative conjugative element (ICE) harbouring the vanD locus. (ii) The native ligase (ddl) was disrupted by an ISEfm1 insertion. (iii) A large 1.74 Mb chromosomal inversion of unknown consequence occurred. Alignment and phylogenetic-based comparisons of the VREfm with a global collection of vanD-harbouring genomes identified strong similarities in the 120-160 kb genomic region surrounding vanD, suggestive of a common mobile element and integration site, irrespective of the diverse taxonomic, geographical and host origins of the isolates. This isolate diversity revealed that this putative ICE (and its source) is globally disseminated and is capable of being acquired by different genera. Although the incidence of vanD VREfm is low, understanding its emergence and potential for spread is crucial for the ongoing efforts to reduce antimicrobial resistance.

Genetically Related Avian Influenza H7N9 Viruses Exhibit Different Pathogenicity in Mice

Avian influenza viruses can cross species barriers and adapt to mammals. The H7N9 subtype AIV that emerged in China in 2013 caused 1568 human infections, with a mortality rate of nearly 40%. We conducted a retrospective analysis of H7N9 viruses that were isolated in live poultry markets in 2013. We found that two avian-origin H7N9 isolates, A/chicken/Eastern China/JTC4/2013 and A/chicken/Eastern China/JTC11/2013, have a similar genetic background but exhibit different pathogenicity in mice. Whole-genome alignment of the two H7N9 viruses was carried out, and only six amino acid differences mapped in five genes, including the well-known virulence molecular marker PB2-E627K. Our retrospective analysis highlighted the importance of monitoring the adaptive mutations in avian influenza viruses with zoonotic potential.

FISH painting for chromosome identification of aneuploid cauliflower (Brassica oleracea L. var. botrytis)

A common problem in the cultivation and breeding of cauliflower (Brassica oleracea L. var. botrytis) is the occurrence of aneuploids in offspring families. To reveal the chromosomal cause of such numerical variants, it was necessary to develop karyotype tools with which chromosomes can be easily identified. Since mitotic chromosomes in this crop are morphologically similar and lack differentiating banding patterns, we tested two Fluorescent in situ Hybridization (FISH) procedures for chromosome identification: (1) FISH painting with diagnostic repetitive DNA patterns and (2) cross-species chromosome painting. The first method consists of a five-colour FISH with 5s rDNA, 45S rDNA, and two Brassica rapa centromere-specific repeats, and a B. rapa BAC (KBrH092N02) containing a dispersed repeat of an unknown class. The second method is an advanced FISH technology based on hybridising DNA probes of a related species under adapted stringency conditions to identify their homoeologous loci. To this end, we applied four pools of BACs from Arabidopsis thaliana in a multicolour FISH for a banding pattern on the chromosomes of cauliflower (Brassica oleracea L. var. botrytis). Due to the genome triplication and various chromosome rearrangements of Brassica oleracea compared to Arabidopsis, we used MUMmer whole-genome alignment plot information to select Arabidopsis BAC pools with which all cauliflower chromosomes could be identified. In a sample of 21 plants with aberrant phenotypes, we demonstrated primary trisomy for chromosomes 1–6 and 8, and telo-trisomy for chromosomes 7 and 9. Finally, we discuss the advantages and drawbacks of the two painting methods and eventual alternatives for demonstrating numerical aberrations in the cauliflower populations.Graphical

GET_PANGENES: calling pangenes from plant genome alignments confirms presence-absence variation

Crop pangenomes made from individual cultivar assemblies promise easy access to conserved genes, but genome content variability and inconsistent identifiers hamper their exploration. To address this, we define pangenes, which summarize a species coding potential and link back to original annotations. The protocol get_pangenes performs whole genome alignments (WGA) to call syntenic gene models based on coordinate overlaps. A benchmark with small and large plant genomes shows that pangenes recapitulate phylogeny-based orthologies and produce complete soft-core gene sets. Moreover, WGAs support lift-over and help confirm gene presence-absence variation. Source code and documentation: https://github.com/Ensembl/plant-scripts.

Maast: genotyping thousands of microbial strains efficiently

Existing single nucleotide polymorphism (SNP) genotyping algorithms do not scale for species with thousands of sequenced strains, nor do they account for conspecific redundancy. Here we present a bioinformatics tool, Maast, which empowers population genetic meta-analysis of microbes at an unrivaled scale. Maast implements a novel algorithm to heuristically identify a minimal set of diverse conspecific genomes, then constructs a reliable SNP panel for each species, and enables rapid and accurate genotyping using a hybrid of whole-genome alignment and k-mer exact matching. We demonstrate Maast’s utility by genotyping thousands of Helicobacter pylori strains and tracking SARS-CoV-2 diversification.

Isolation, biochemical characterization, and genome sequencing of two high-quality genomes of a novel chitinolytic Jeongeupia species.

Chitin is the second most abundant polysaccharide worldwide as part of arthropods' exoskeletons and fungal cell walls. Low concentrations in soils and sediments indicate rapid decomposition through chitinolytic organisms in terrestrial and aquatic ecosystems. The enacting enzymes, so-called chitinases, and their products, chitooligosaccharides, exhibit promising characteristics with applications ranging from crop protection to cosmetics, medical, textile, and wastewater industries. Exploring novel chitinolytic organisms is crucial to expand the enzymatical toolkit for biotechnological chitin utilization and to deepen our understanding of diverse catalytic mechanisms. In this study, we present two long-read sequencing-based genomes of highly similar Jeongeupia species, which have been screened, isolated, and biochemically characterized from chitin-amended soil samples. Through metabolic characterization, whole-genome alignments, and phylogenetic analysis, we could demonstrate how the investigated strains differ from the taxonomically closest strain Jeongeupia naejangsanensis BIO-TAS4-2T (DSM 24253). In silico analysis and sequence alignment revealed a multitude of highly conserved chitinolytic enzymes in the investigated Jeongeupia genomes. Based on these results, we suggest that the two strains represent a novel species within the genus of Jeongeupia, which may be useful for environmentally friendly N-acetylglucosamine production from crustacean shell or fungal biomass waste or as a crop protection agent.