Genomic GC Research Articles

Convergent reductive evolution in bee-associated lactic acid bacteria.

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae [lactic acid bacteria (LAB)] family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated LAB), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts toward the bee environment in LAB. In these species, we observed significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, acetaldehyde-alcohol dehydrogenase (adhE), encodes a bifunctional aldehyde-alcohol dehydrogenase which has been associated with the evolution of fructophily, a rare phenotypic trait that is pervasive across bee-associated LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent losses of the same sets of genes.IMPORTANCESeveral LAB species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine learning-based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference for fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia.

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

Genome Survey of Male Rana dybowskii to Further Understand the Sex Determination Mechanism.

Rana dybowskii is one of the important aquaculture species in Northeast China. The fallopian tubes of female R. dybowskii are used to prepare oviductus ranae (an important traditional Chinese medicine). Therefore, R. dybowskii females have higher economical value than males. An increasing female R. dybowskii population can increase the benefits from R. dybowskii culture. However, the genome of amphibians is complex, making it difficult to investigate their sex determination mechanism. In this study, we analyzed the genome of male R. dybowskii using next-generation sequencing technology. A total of 200,046,452,400 bp of clean data were obtained, and the K-mer analysis indicated that the depth was 50×. The genome size of R. dybowskii was approximately 3585.05 M, with a heterozygosity rate, repeat sequence ratio, and genome GC content of 1.15%, 68.96%, and approximately 43.0%, respectively. In total, 270,785 contigs and 498 scaffolds were generated. The size of the contigs and scaffolds was 3,748,543,415 and 3,765,862,278 bp, respectively, with the N50 length of 31,988 and 336,385,783. The longest contig and scaffold were of the size 137,967,485 and 1,808,367,828 bp, respectively. The number of contigs and scaffolds > 10K nt was 99,620 and 451, respectively. Through annotation, 40,913 genes were obtained, including 156,609 CDS (i.e., 3.83 CDS per gene). Sequence alignment was performed with the assembled scaffolding genome in this study. Two and one fragment had high homology with two male-specific DNA molecular markers of R. dybowskii discovered previously (namely, MSM-222 and MSM-261, respectively). In addition, the Dmrt1 gene of R. dybowskii was obtained with a length of 18,893 bp by comparison and splicing. The forward primers amplifying MSM-222 and MSM-261 were located at 322-343 and 14,501-14,526 bp of Dmrt1, respectively. However, sequence alignment revealed that MSM-222 and MSM-261 were not located on Dmrt1, and only some homologous parts were observed. This indicated that in addition to Dmrt1, other important genes may play a crucial role in the sex determination mechanism of R. dybowskii. Our study provided a foundation for the subsequent high-quality genome construction and provided important genomic resources for future studies on R. dybowskii.

Compression rates of microbial genomes are associated with genome size and base composition

BackgroundTo what degree a string of symbols can be compressed reveals important details about its complexity. For instance, strings that are not compressible are random and carry a low information potential while the opposite is true for highly compressible strings. We explore to what extent microbial genomes are amenable to compression as they vary considerably both with respect to size and base composition. For instance, microbial genome sizes vary from less than 100,000 base pairs in symbionts to more than 10 million in soil-dwellers. Genomic base composition, often summarized as genomic AT or GC content due to the similar frequencies of adenine and thymine on one hand and cytosine and guanine on the other, also vary substantially; the most extreme microbes can have genomes with AT content below 25% or above 85% AT. Base composition determines the frequency of DNA words, consisting of multiple nucleotides or oligonucleotides, and may therefore also influence compressibility. Using 4,713 RefSeq genomes, we examined the association between compressibility, using both a DNA based- (MBGC) and a general purpose (ZPAQ) compression algorithm, and genome size, AT content as well as genomic oligonucleotide usage variance (OUV) using generalized additive models.ResultsWe find that genome size (p < 0.001) and OUV (p < 0.001) are both strongly associated with genome redundancy for both type of file compressors. The DNA-based MBGC compressor managed to improve compression with approximately 3% on average with respect to ZPAQ. Moreover, MBGC detected a significant (p < 0.001) compression ratio difference between AT poor and AT rich genomes which was not detected with ZPAQ.ConclusionAs lack of compressibility is equivalent to randomness, our findings suggest that smaller and AT rich genomes may have accumulated more random mutations on average than larger and AT poor genomes which, in turn, were significantly more redundant. Moreover, we find that OUV is a strong proxy for genome compressibility in microbial genomes. The ZPAQ compressor was found to agree with the MBGC compressor, albeit with a poorer performance, except for the compressibility of AT-rich and AT-poor/GC-rich genomes.

Comparative genome analysis of endophytic Bacillus amyloliquefaciens MR4: a potential biocontrol agent isolated from wild medicinal plant root tissue.

In this investigation, 396 endophytic bacterial strains from six indigenous medicinal plant species within the Xinjiang Tumor Peak National Nature Reserve were subjected to screening. The strain MR4 emerged as a noteworthy contender, demonstrating pronounced biocontrol capabilities coupled with exceptional cold tolerance. Through morphological scrutiny and comprehensive genomic sequencing, MR4 was identified as Bacillus amyloliquefaciens. Antagonistic assays revealed MR4's efficacy in suppressing the causative agents of cotton wilt and verticillium wilt, achieving inhibition rates surpassing 50%. Analyses, underpinned by PCR methodologies, indicated MR4's capacity to biosynthesize a minimum of eight distinct antimicrobial agents. The whole-genome sequencing data indicated that B. amyloliquefaciens MR4 had the genome size and GC content of 4,017,872bp and 47.14%, respectively, and 4191 coding genes were identified. The genome consists of a single chromosome and one plasmid. Moreover, it was augmented by annotations from various databases, including GO, KEGG, and COG. The pathogenicity of MR4 undergoes evaluation, while predictions concerning MR4's secondary metabolites have disclosed gene clusters for 13 varieties of these compounds, with particular emphasis on surfactins and fengycin. Comparative analyses with four paradigmatic strains shed light on MR4's genomic composition and its phylogenetic lineage within the Bacillus genus. The genomic data pertaining to MR4 have been duly submitted to the NCBI GenBank, bearing the accession numbers CP146236 (Chr1) and CP146237 (plas1). This study endeavors to furnish potent microbial resources for the biocontrol and enhancement of plant growth, thereby providing a theoretical groundwork for MR4's agronomic utilization.

Genomic characterization of an uncommon Delftia acidovorans isolate obtained from aBulgarian immunocompetent outpatient diagnosed with bronchitis.

Delftia acidovorans is an aerobic, non-fermenting Gram-negative bacterium (NFGNB), found in soil, water and hospital environments. It is rarely clinically significant, most commonly affecting hospitalized or immunocompromised patients. The present study aimed to explore the genomic characteristics of a Bulgarian clinical D. acidovorans isolate (designated Dac759) in comparison to all strains of this species with available genomes in the NCBI Genome database (n = 34). Dac759 was obtained in 2021 from the sputum of a 65-year-old female immunocompetent outpatient with bronchitis. Species identification using MALDI-TOF mass spectrometry, antimicrobial susceptibility testing, whole-genome sequencing (WGS), and phylogenomic analysis were performed. The isolate demonstrated high-level resistance to colistin (16 mg L-1); resistance to gentamicin; reduced susceptibility to piperacillin, piperacillin-tazobactam, ceftazidime, cefepime, ciprofloxacin, and levofloxacin; and susceptibility to imipenem, meropenem, amikacin, and tobramycin. The observed genome size (6.43 Mb) and GC content (66.76%) were comparable with the accessible data from sequenced D. acidovorans genomes. A limited number of resistance determinants were identified in the assembled genome as follows: blaOXA-459, emrE, oqxB, and mexCD-oprJ. The phylogenomic analysis indicated a high heterogenicity of the included D. acidovorans genomes. In conclusion, to the best of our knowledge, this is the first documented case of a clinically relevant D. acidovorans isolate in Bulgaria. Unlike the majority of reports in the literature, Dac759 affected a patient with no malignancies or other preexisting comorbidities. With this in mind, its genome sequence is a valuable resource for the fundamental study of uncommon bacterial pathogens of public health importance.

Genomic factors shaping codon usage across the Saccharomycotina subphylum.

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy greater than 90% using a machine-learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Deciphering the complete chloroplast genome sequence of Meconopsis torquata Prain: Insights into genome structure, comparative analysis and phylogenetic relationship

In the present study, we have characterized the complete chloroplast (Cp) genome of Meconopsis torquata Prain (family Papaveraceae), revealing the plastome size of 153,290 bp, and a GC content of 38.72 %. The cp genome features the typical circular quadripartite structure found in flowering plants, including a pair of inverted repeat regions (25,816 bp), isolated by a small single-copy region (17,740 bp) and a large single-copy (83,918 bp). Genome annotation revealed 132 genes: 87 protein-coding genes, 37 tRNAs and eight rRNAs. This comparative study demonstrated that the genome structure, gene number and GC ratio are consistent with several other cp genomes of Meconopsis and Papaver genera. A total of 120 SSRs were detected in the plastome, the majority (111) of which were mononucleotide repeats. Among the longer repeats, palindromic sequences were most common, followed by forward, reverse, and complement repeats. The whole genome alignment revealed the conserved nature of the inverted repeat region over single-copy zones. Nucleotide diversity unveiled hypervariable sites (ycf1, rps16, accD, atpB and psbD) in both the small and large single-copy regions, which could be useful for designing molecular markers for taxonomic identification. Phylogenetic analysis revealed a close alliance of M. torquata with other Meconopsis species, such as M. pinnatifolia and M. paniculata, with strong bootstrap support. Molecular dating suggests that M. torquata originated during the Tortonian age of the Miocene epoch of the Cenozoic era. These findings provide valuable insights for biological research, especially in understanding the genetic and evolutionary divergence within the Papaveraceae family.

Description and Whole-Genome Sequencing of Mariniflexile litorale sp. nov., Isolated from the Shallow Sediments of the Sea of Japan.

A Gram-negative, aerobic, rod-shaped, non-motile, yellow-pigmented bacterium, KMM 9835T, was isolated from the sediment sample obtained from the Amur Bay of the Sea of Japan seashore, Russia. Phylogenetic analyses based on the 16S rRNA gene and whole genome sequences positioned the novel strain KMM 9835T in the genus Mariniflexile as a separate line sharing the highest 16S rRNA gene sequence similarities of 96.6% and 96.2% with Mariniflexile soesokkakense RSSK-9T and Mariniflexile fucanivorans SW5T, respectively, and similarity values of <96% to other recognized Mariniflexile species. The average nucleotide identity and digital DNA-DNA hybridization values between strain KMM 9835T and M. soesokkakense KCTC 32427T, Mariniflexile gromovii KCTC 12570T, M. fucanivorans DSM 18792T, and M. maritimum M5A1MT were 83.0%, 82.5%, 83.4%, and 78.3% and 30.7%, 29.6%, 29.5%, and 24.4%, respectively. The genomic DNA GC content of strain KMM 9835T was 32.5 mol%. The dominant menaquinone was MK-6, and the major fatty acids were iso-C15:0, iso-C15:1ω10c, and C15:0. The polar lipids of strain KMM 9835T consisted of phosphatidylethanolamine, two unidentified aminolipids, an unidentified phospholipid, and six unidentified lipids. A pan-genome analysis showed that the KMM 9835T genome encoded 753 singletons. The annotated singletons were more often related to transport protein systems (SusC), transcriptional regulators (AraC, LytTR, LacI), and enzymes (glycosylases). The KMM 9835T genome was highly enriched in CAZyme-encoding genes, the proportion of which reached 7.3%. Moreover, the KMM 9835T genome was characterized by a high abundance of CAZyme gene families (GH43, GH28, PL1, PL10, CE8, and CE12), indicating its potential to catabolize pectin. This may represent part of an adaptation strategy facilitating microbial consumption of plant polymeric substrates in aquatic environments near shorelines and freshwater sources. Based on the combination of phylogenetic and phenotypic characterization, the marine sediment strain KMM 9835T (=KCTC 92792T) represents a novel species of the genus Mariniflexile, for which the name Mariniflexile litorale sp. nov. is proposed.

Assembly and Comparative Analyses of the Chloroplast Genomes of the Threatened Plant Rosa anemoniflora

Due to insufficient molecular biology and genetic research on Rosa anemoniflora, this endangered plant has not yet received effective protection. Therefore, the complete chloroplast genome sequence of R. anemoniflora, along with comparative analysis of the chloroplast genomes of related species, is necessary and crucial for reconstructing phylogenetic relationships and developing genetic markers to conserve these species. A series of analyses, including genome structure, GC content, gene number, selection pressure, and nucleotide diversity, were performed by comparing the chloroplast genomes of R. anemoniflora and its relatives. The results indicate that the chloroplast genomes of R. anemoniflora and its close relatives are highly conserved in all genome characteristics, and all protein-coding genes in R. anemoniflora have not experienced significant positive or negative selection pressures. Comparative analysis revealed several variation hotspots, such as the atpH-atpI region, which can serve as a DNA barcode for distinguishing R. anemoniflora from its close relatives. Finally, the results confirmed that R. anemoniflora belongs to Rosa section Synstylae and that R. anemoniflora and its close relatives likely originated from the Fujian and Taiwan regions of China and diverged approximately 3.24 million years ago. This study provides crucial information for future biodiversity conservation and genetic resource management.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Identification and classification of the genomes of novel microviruses in poultry slaughterhouse.

Microviridae is a family of phages with circular ssDNA genomes and they are widely found in various environments and organisms. In this study, virome techniques were employed to explore potential members of Microviridae in a poultry slaughterhouse, leading to the identification of 98 novel and complete microvirus genomes. Using a similarity clustering network classification approach, these viruses were found to belong to at least 6 new subfamilies within Microviridae and 3 higher-level taxonomic units. Genome size, GC content and genome structure of these new taxa showed evident regularities, validating the rationality of our classification method. Our method can divide microviruses into about 45 additional detailed clusters, which may serve as a new standard for classifying Microviridae members. Furthermore, by addressing the scarcity of host information for microviruses, the current study significantly broadened their host range and discovered over 20 possible new hosts, including important pathogenic bacteria such as Helicobacter pylori and Vibrio cholerae, as well as different taxa demonstrated different host specificities. The findings of this study effectively expand the diversity of the Microviridae family, providing new insights for their classification and identification. Additionally, it offers a novel perspective for monitoring and controlling pathogenic microorganisms in poultry slaughterhouse environments.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BackgroundGenomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.ResultsHere, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.ConclusionsOur results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

Comparative Genomics Reveals Genetic Diversity and Variation in Metabolic Traits in Fructilactobacillus sanfranciscensis Strains.

Fructilactobacillus sanfranciscensis is a significant and dominant bacterial species of sourdough microbiota from ecological and functional perspectives. Despite the remarkable prevalence of different strains of this species in sourdoughs worldwide, the drivers behind the genetic diversity of this species needed to be clarified. In this research, 14 F. sanfranciscensis strains were isolated from sourdough samples to evaluate the genetic diversity and variation in metabolic traits. These 14 and 31 other strains (obtained from the NCBI database) genomes were compared. The values for genome size and GC content, on average, turned out to 1.31 Mbp and 34.25%, respectively. In 45 F. sanfranciscensis strains, there were 162 core genes and 0 to 51 unique genes present in each strain. The primary functions of core genes were related to nucleotide, lipid transport, and amino acid, as well as carbohydrate metabolism. The size of core genes accounted for 41.18% of the pan-genome size in 14 F. sanfranciscensis strains, i.e., 0.70 Mbp of 1.70 Mbp. There were genetic variations among the 14 strains involved in carbohydrate utilization and antibiotic resistance. Moreover, exopolysaccharides biosynthesis-related genes were annotated, including epsABD, wxz, wzy. The Type IIA & IE CRISPR-Cas systems, pediocin PA-1 and Lacticin_3147_A1 bacteriocins operons were also discovered in F. sanfranciscensis. These findings can help to select desirable F. sanfranciscensis strains to develop standardized starter culture for sourdough fermentation, and expect to provide traditional fermented pasta with a higher quality and nutritional value for the consumers.

Genomic comparative analysis of Ophiocordyceps unilateralis sensu lato.

Ophiocordyceps unilateralis sensu lato is a common pathogenic fungus of ants. A new species, O. fusiformispora, was described based on morphology and phylogenetic evidence from five genes (SSU, LSU, TEF1α, RPB1, and RPB2). The whole genomes of O. fusiformispora, O. contiispora, O. subtiliphialida, O. satoi, O. flabellata, O. acroasca, and O. camponoti-leonardi were sequenced and annotated and compared with whole genome sequences of other species in O. unilateralis sensu lato. The basic genome-wide characteristics of the 12 species showed that the related species had similar GC content and genome size. AntiSMASH and local BLAST analyses revealed that the number and types of putative SM BGCs, NPPS, PKS, and hybrid PKS-NRPS domains for the 12 species differed significantly among different species in the same genus. The putative BGC of five compounds, namely, NG-391, lucilactaene, higginsianin B, pyripyropene A, and pyranonigrin E were excavated. NG-391 and lucilactaene were 7-desmethyl analogs of fusarin C. Furthermore, the 12 genomes had common domains, such as KS-AT-DH-MT-ER-KR-ACP and SAT-KS-AT-PT-ACP-ACP-Te. The ML and BI trees of SAT-KS-AT-PT-ACP-ACP-Te were highly consistent with the multigene phylogenetic tree in the 12 species. This study provided a method to obtain the living culture of O. unilateralis sensu lato species and its asexual formed on the basis of living culture, which was of great value for further study of O. unilateralis sensu lato species in the future, and also laid a foundation for further analysis of secondary metabolites of O. unilateralis sensu lato.

Phylogeny Influences Genome Size and GC but not Sequence and Organismal Complexity in Staphylococci

Phylogeny Influences Genome Size and GC but not Sequence and Organismal Complexity in Staphylococci

Abstract 318: Demystifying tumor heterogeneity with a fully automated, high-throughput single-cell DNA-Seq (scDNA-seq) workflow

Abstract Objective: Genetic heterogeneity is a key factor underlying tumor drug resistance and metastatic potential. Understanding this heterogeneity is therefore vital to improving both prognosis and treatment. Next-generation sequencing is a valuable tool for analyzing the genetic makeup of tumors. However, bulk sequencing methods lack the sensitivity to fully resolve tumor heterogeneity. While single cell methods provide a powerful approach to dissect such heterogeneity, to date, these methods have been limited in their throughput. To address this bottleneck, we have developed a fully automated workflow for generating scDNA-seq libraries based on PicoPLEX® whole genome amplification (WGA) technology. This high-throughput method, which has been optimized on our ICELL8® cx Single-Cell System, enables the generation of WGA libraries for &amp;gt;1,000 single cells within one day. Methods: We first sought to demonstrate the ability of this new high-throughput method to generate WGA libraries of comparable quality to the standard PicoPLEX workflow in terms of genome coverage, GC bias, and other typical quality metrics. We additionally assessed copy number variant (CNV) detection sensitivity using two cell lines with well-characterized small segmental CNVs: GM22601 (~25Mb deletion on chromosome 4) and GM05067 (~45Mb gain on chromosome 9). We also analyzed a lymphoblastoid line (K562) that carries a range of chromosomal aneuploidies. A total of 1,288 single-cell WGA libraries were generated and sequenced to a depth of 250K paired-end reads per cell. Data analysis was done using the Ginkgo CNV pipeline with an average bin size of 500 kb. As a final proof of principal, we also generated single cell data using tumor and adjacent normal tissue from two clear cell renal cell carcinoma (ccRCC) samples. Results: Libraries generated using our high-throughput workflow had a high mapping rate, with 94.1% of the reads being uniquely mapped. Additionally, the libraries were comparable to those generated with the standard PicoPLEX workflow in terms of coverage uniformity, GC bias and other metrics. Moreover, the segmental aneuploidies in both GM22601 and GM05067 were reliably detected in &amp;gt;90% of cells at a read depth as low as 250,000 reads per cell. Analysis of the ccRCC samples revealed subclonal heterogeneity with various CNVs common to ccRCC, including deletion of chr. 3p, amplification of chr. 5q, and duplication of chr. 2. Conclusion: By adapting PicoPLEX technology to high throughput using the ICELL8 cx Single-Cell System, we have obtained single-cell WGA libraries from up to 1,200 cells at once and enabled the reliable detection of CNVs and tumor subclones at a shallow sequencing depth. In addition, by leveraging the automated nanoliter-dispensing capabilities of ICELL8 cx system this method provides a significant reduction in reagent use and labor compared to plate-based methods. Citation Format: Xuan Li, Raymond Mendoza, Samantha Leong, Hima Anbunathan, Mike Covington, Mohammad Fallahi, Bryan Bell, Shuwen Chen, Yue Yun, Andrew Farmer. Demystifying tumor heterogeneity with a fully automated, high-throughput single-cell DNA-Seq (scDNA-seq) workflow [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 318.

An evolutionary view of the Fusarium core genome

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

Comparative genomics and phylogenomics of the genus Glycyrrhiza (Fabaceae) based on chloroplast genomes.

Glycyrrhiza (Fabaceae) species are rich in metabolites and widely used in medicine. Research on the chloroplast genome of Glycyrrhiza is important for understanding its phylogenetics, biogeography, genetic diversity, species identification, and medicinal properties. In this study, comparative genomics and phylogenomics of Glycyrrhiza were analyzed based on the chloroplast genome. The chloroplast genomes of six Glycyrrhiza species were obtained using various assembly and annotation tools. The final assembled chloroplast genome sizes for the six Glycyrrhiza species ranged from 126,380 bp to 129,115 bp, with a total of 109-110 genes annotated. Comparative genomics results showed that the chloroplast genomes of Glycyrrhiza showed typically lacking inverted repeat regions, and the genome length, structure, GC content, codon usage, and gene distribution were highly similar. Bioinformatics analysis revealed the presence of 69-96 simple sequence repeats and 61-138 long repeats in the chloroplast genomes. Combining the results of mVISTA and nucleotide diversity, four highly variable regions were screened for species identification and relationship studies. Selection pressure analysis indicated overall purifying selection in the chloroplast genomes of Glycyrrhiza, with a few positively selected genes potentially linked to environmental adaptation. Phylogenetic analyses involving all tribes of Fabaceae with published chloroplast genomes elucidated the evolutionary relationships, and divergence time estimation estimated the chronological order of species differentiations within the Fabaceae family. The results of phylogenetic analysis indicated that species from the six subfamilies formed distinct clusters, consistent with the classification scheme of the six subfamilies. In addition, the inverted repeat-lacking clade in the subfamily Papilionoideae clustered together, and it was the last to differentiate. Co-linear analysis confirmed the conserved nature of Glycyrrhiza chloroplast genomes, and instances of gene rearrangements and inversions were observed in the subfamily Papilionoideae.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BackgroundFirmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.ResultsThe F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.ConclusionsThe characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.