Elucidating the origins of phycocyanobilin biosynthesis and phycobiliproteins.

Terrestrial ecosystems and human societies depend on oxygenic photosynthesis, which began to reshape our atmosphere approximately 2.5 billion years ago. The earliest known organisms carrying out oxygenic photosynthesis are the cyanobacteria, which use large complexes of phycobiliproteins as light-harvesting antennae. Phycobiliproteins rely on phycocyanobilin (PCB), a linear tetrapyrrole (bilin) chromophore, as the light-harvesting pigment that transfers absorbed light energy from phycobilisomes to the chlorophyll-based photosynthetic apparatus. Cyanobacteria synthesize PCB from heme in two steps: A heme oxygenase converts heme into biliverdin IXα (BV), and the ferredoxin-dependent bilin reductase (FDBR) PcyA then converts BV into PCB. In the current work, we examine the origins of this pathway. We demonstrate that PcyA evolved from pre-PcyA proteins found in nonphotosynthetic bacteria and that pre-PcyA enzymes are active FDBRs that do not yield PCB. Pre-PcyA genes are associated with two gene clusters. Both clusters encode bilin-binding globin proteins, phycobiliprotein paralogs that we designate as BBAGs (bilin biosynthesis-associated globins). Some cyanobacteria also contain one such gene cluster, including a BBAG, two V4R proteins, and an iron-sulfur protein. Phylogenetic analysis shows that this cluster is descended from those associated with pre-PcyA proteins and that light-harvesting phycobiliproteins are also descended from BBAGs found in other bacteria. We propose that PcyA and phycobiliproteins originated in heterotrophic, nonphotosynthetic bacteria and were subsequently acquired by cyanobacteria.

VirID: Beyond Virus Discovery - An Integrated Platform for Comprehensive RNA Virus Characterization.

RNA viruses exhibit vast phylogenetic diversity and can significantly impact public health and agriculture. However, current bioinformatics tools for viral discovery from metagenomic data frequently generate false positive virus results, overestimate viral diversity, and misclassify virus sequences. Additionally, current tools often fail to determine virus-host associations, which hampers investigation of the potential threat posed by a newly detected virus. To address these issues we developed VirID, a software tool specifically designed for the discovery and characterization of RNA viruses from metagenomic data. The basis of VirID is a comprehensive RNA-dependent RNA polymerase (RdRP) database to enhance a workflow that includes RNA virus discovery, phylogenetic analysis, and phylogeny-based virus characterization. Benchmark tests on a simulated data set demonstrated that VirID had high accuracy in profiling viruses and estimating viral richness. In evaluations with real-world samples, VirID was able to identity RNA viruses of all type, but also provided accurate estimations of viral genetic diversity and virus classification, as well as comprehensive insights into virus associations with humans, animals, and plants. VirID therefore offers a robust tool for virus discovery and serves as a valuable resource in basic virological studies, pathogen surveillance, and early warning systems for infectious disease outbreaks.

phyloBARCODER: A Web Tool for Phylogenetic Classification of Eukaryote Metabarcodes Using Custom Reference Databases.

We developed phyloBARCODER (https://github.com/jun-inoue/phyloBARCODER), a new web tool that can identify short DNA sequences to the species level using metabarcoding. phyloBARCODER estimates phylogenetic trees based on the uploaded anonymous DNA sequences and reference sequences from databases. Without such phylogenetic contexts, alternative, similarity-based methods independently identify species names and anonymous sequences of the same group by pairwise comparisons between queries and database sequences, with the caveat that they must match exactly or very closely. By putting metabarcoding sequences into a phylogenetic context, phyloBARCODER accurately identifies (i) species or classification of query sequences and (ii) anonymous sequences associated with the same species or even with populations of query sequences, with clear and accurate explanations. Version 1 of phyloBARCODER stores a database comprising all eukaryotic mitochondrial gene sequences. Moreover, by uploading their own databases, phyloBARCODER users can conduct species identification specialized for sequences obtained from a local geographic region or those of nonmitochondrial genes, e.g. ITS or rbcL.

HGTphyloDetect: facilitating the identification and phylogenetic analysis of horizontal gene transfer.

BACKGROUND: Horizontal gene transfer (HGT) is an important driver in genome evolution, gain-of-function, and metabolic adaptation to environmental niches. Genome-wide identification of putative HGT events has become increasingly practical, given the rapid growth of genomic data. However, existing HGT analysis toolboxes are not widely used, limited by their inability to perform phylogenetic reconstruction to explore potential donors, and the detection of HGT from both evolutionarily distant and closely related species. RESULTS: In this study, we have developed HGTphyloDetect, which is a versatile computational toolbox that combines high-throughput analysis with phylogenetic inference, to facilitate comprehensive investigation of HGT events. Two case studies with Saccharomyces cerevisiae and Candida versatilis demonstrate the ability of HGTphyloDetect to identify horizontally acquired genes with high accuracy. In addition, HGTphyloDetect enables phylogenetic analysis to illustrate a likely path of gene transmission among the evolutionarily distant or closely related species. CONCLUSIONS: The HGTphyloDetect computational toolbox is designed for ease of use and can accurately find HGT events with a very low false discovery rate in a high-throughput manner. The HGTphyloDetect toolbox and its related user tutorial are freely available at https://github.com/SysBioChalmers/HGTphyloDetect.

Assessing the Adequacy of Morphological Models using Posterior Predictive Simulations.

Reconstructing the evolutionary history of different groups of organisms provides insight into how life originated and diversified on Earth. Phylogenetic trees are commonly used to estimate this evolutionary history. Within Bayesian phylogenetics a major step in estimating a tree is in choosing an appropriate model of character evolution. While the most common character data used is molecular sequence data, morphological data remains a vital source of information. The use of morphological characters allows for the incorporation fossil taxa, and despite advances in molecular sequencing, continues to play a significant role in neontology. Moreover, it is the main data source that allows us to unite extinct and extant taxa directly under the same generating process. We therefore require suitable models of morphological character evolution, the most common being the Mk Lewis model. While it is frequently used in both palaeobiology and neontology, it is not known whether the simple Mk substitution model, or any extensions to it, provide a sufficiently good description of the process of morphological evolution. In this study we investigate the impact of different morphological models on empirical tetrapod data sets. Specifically, we compare unpartitioned Mk models with those where characters are partitioned by the number of observed states, both with and without allowing for rate variation across sites and accounting for ascertainment bias. We show that the choice of substitution model has an impact on both topology and branch lengths, highlighting the importance of model choice. Through simulations, we validate the use of the model adequacy approach, posterior predictive simulations, for choosing an appropriate model. Additionally, we compare the performance of model adequacy with Bayesian model selection. We demonstrate how model selection approaches based on marginal likelihoods are not appropriate for choosing between models with partition schemes that vary in character state space (i.e., that vary in Q-matrix state size). Using posterior predictive simulations, we found that current variations of the Mk model are often performing adequately in capturing the evolutionary dynamics that generated our data. We do not find any preference for a particular model extension across multiple data sets, indicating that there is no 'one size fits all' when it comes to morphological data and that careful consideration should be given to choosing models of discrete character evolution. By using suitable models of character evolution, we can increase our confidence in our phylogenetic estimates, which should in turn allow us to gain more accurate insights into the evolutionary history of both extinct and extant taxa.

Deciphering the genetic basis of resistome and virulome diversity among multidrug-resistant Mycoplasma hominis.

Mycoplasma hominis, a commensal bacterium that commonly inhabits the genital tract, leading to infections in both the genitourinary and extragenital regions. However, the antimicrobial resistance and pathogenic mechanisms of M. hominis isolated from extra-urogenital cystic abscess is largely unknown. This study reports the genomic epidemiological characteristics of a M. hominis isolate recovered from a pelvic abscess sample in China. Genomic DNA was extracted and sequenced using Illumina HiSeq X Ten platform. De novo assembly was performed and in silico analysis was accomplished by multiple bioinformatics tools. For phylogenomic analysis, publicly available M. hominis genomes were retrieved from NCBI GenBank database. Whole genome sequencing data showed that the genome size of M. hominis MH4246 was calculated as 679,746 bp, with 558 protein-coding sequences and a G + C content of 26.9%. M. hominis MH4246 is resistant to fluoroquinolones and macrolides, harboring mutations in the quinolone resistance-determining regions (QRDRs) (GyrA S153L, ParC S91I and ParE V417I) and 23S rRNA gene (G280A, C1500T, T1548C and T2218C). Multiple virulence determinants, such as tuf, hlyA, vaa, oppA, MHO_0730 and alr genes, were identified. Phylogenetic analysis showed that the closest relative of M. hominis MH4246 was the strain MH-1 recovered from China, which differed by 3490 SNPs. Overall, this study contributes to the comprehension of genomic characteristics, antimicrobial resistance patterns, and the mechanisms underlying the pathogenicity of this pathogen.

Assembly and comparative analysis of the complete mitochondrial genome of white towel gourd (Luffa cylindrica).

Mitochondria play an important role in the energy production of plant cells through independent genetic systems. This study has aimed to assemble and annotate the functions of the mitochondrial (mt) genome of Luffa cylindrica. The mt genome of L. cylindrica contained two chromosomes with lengths of 380,879 bp and 67,982 bp, respectively. Seventy-seven genes including 39 protein-coding genes, 34 tRNA genes, 3 rRNA genes, and 1 pseudogene, were identified. About 90.63% of the codons ended with A or U bases, and 98.63% of monomers contained A/T, which contributed to the high A/T content (55.91%) of the complete mt genome. Six genes (ATP8, CCMFC, NAD4, RPL10, RPL5 and RPS4) showed positive selection. Phylogenetic analysis indicates that L. cylindrica is closely related to L. acutangula. The present results provide the mt genome of L. cylindrica, which may facilitate possible genetic variation, evolutionary, and molecular breeding studies of L. cylindrica.

Characterizing the etiology of recurrent tuberculosis using whole genome sequencing-Alaska, USA, 2008-2020.

BACKGROUND: Understanding the etiology of recurrent tuberculosis (rTB) is important for effective TB control. Prior to the advent of whole genome sequencing (WGS), attributing rTB to relapse or reinfection using genetic information was complicated by the limited resolution of conventional genotyping methods. METHODS: We applied a systematic method of evaluating whole genome single nucleotide polymorphism (wgSNP) distances and results of phylogenetic analyses to characterize the etiology of rTB in American Indian and Alaska Native (AIAN) persons in Alaska during 2008-2020. We contextualized our findings through descriptive analyses of surveillance data and results of a literature search for investigations that characterized rTB etiology using WGS. RESULTS: The percentage of TB cases in AIAN persons in Alaska classified as recurrent episodes (11.8%) was three times the national percentage (3.9%). Of 38 recurrent episodes included in genetic analyses, we attributed 25 (65.8%) to reinfection based on wgSNP distances and phylogenetic analyses; this proportion was the highest among 16 published point estimates identified through the literature search. By comparison, we attributed 11 of 38 (28.9%) and 6 of 38 (15.8%) recurrent episodes to reinfection based on wgSNP distances alone and on conventional genotyping methods, respectively. CONCLUSIONS: WGS and attribution criteria involving genetic distances and patterns of relatedness can provide an effective means of elucidating rTB etiology. Our findings indicate that rTB occurs at high proportions among AIAN persons in Alaska and is frequently attributable to reinfection, reinforcing the importance of active surveillance and control measures to limit the spread of TB disease in Alaskan AIAN communities.

Fern cell walls and the evolution of arabinogalactan proteins in streptophytes.

Significant changes have occurred in plant cell wall composition during evolution and diversification of tracheophytes. As the sister lineage to seed plants, knowledge on the cell wall of ferns is key to track evolutionary changes across tracheophytes and to understand seed plant-specific evolutionary innovations. Fern cell wall composition is not fully understood, including limited knowledge of glycoproteins such as the fern arabinogalactan proteins (AGPs). Here, we characterize the AGPs from the leptosporangiate fern genera Azolla, Salvinia, and Ceratopteris. The carbohydrate moiety of seed plant AGPs consists of a galactan backbone including mainly 1,3- and 1,3,6-linked pyranosidic galactose, which is conserved across the investigated fern AGPs. Yet, unlike AGPs of angiosperms, those of ferns contained the unusual sugar 3-O-methylrhamnose. Besides terminal furanosidic arabinose, Ara (Araf), the main linkage type of Araf in the ferns was 1,2-linked Araf, whereas in seed plants 1,5-linked Araf is often dominating. Antibodies directed against carbohydrate epitopes of AGPs supported the structural differences between AGPs of ferns and seed plants. Comparison of AGP linkage types across the streptophyte lineage showed that angiosperms have rather conserved monosaccharide linkage types; by contrast bryophytes, ferns, and gymnosperms showed more variability. Phylogenetic analyses of glycosyltransferases involved in AGP biosynthesis and bioinformatic search for AGP protein backbones revealed a versatile genetic toolkit for AGP complexity in ferns. Our data reveal important differences across AGP diversity of which the functional significance is unknown. This diversity sheds light on the evolution of the hallmark feature of tracheophytes: their elaborate cell walls.

GERDH: an interactive multi-omics database for cross-species data mining in horticultural crops.

Horticultural plants contribute immensely to the quality of human's life. The rapid development of omics studies on horticultural plants has resulted in large volumes of valuable growth- and development-related data. Genes that are essential for growth and development are highly conserved in evolution. Cross-species data mining reduces the impact of species heterogeneity and has been extensively used for conserved gene identification. Owing to the lack of a comprehensive database for cross-species data mining using multi-omics data from all horticultural plant species, the current resources in this field are far from satisfactory. Here, we introduce GERDH (https://dphdatabase.com), a database platform for cross-species data mining among horticultural plants, based on 12 961 uniformly processed publicly available omics libraries from more than 150 horticultural plant accessions, including fruits, vegetables and ornamental plants. Important and conserved genes that are essential for a specific biological process can be obtained by cross-species analysis module with interactive web-based data analysis and visualization. Moreover, GERDH is equipped with seven online analysis tools, including gene expression, in-species analysis, epigenetic regulation, gene co-expression, enrichment/pathway and phylogenetic analysis. By interactive cross-species analysis, we identified key genes contributing to postharvest storage. By gene expression analysis, we explored new functions of CmEIN3 in flower development, which was validated by transgenic chrysanthemum analysis. We believe that GERDH will be a useful resource for key gene identification and will allow for omics big data to be more available and accessible to horticultural plant community members.

Structural Classification Insights Into the Plant Defensive Peptides.

This study presents a structural phylogenetic analysis of plant defensive peptides, revealing their evolutionary relationships, structural diversification, and functional adaptations. Utilizing a robust dataset comprising both experimental and predicted structures sourced from the RCSB Protein Data Bank and AlphaFold DB, we constructed a detailed phylogenetic tree to elucidate the distinct evolutionary paths of plant defensive peptide families. Our findings showcase the evolutionary intricacies of defensive peptides, highlighting their diversity and the conservation of key structural motifs critical to their antimicrobial or defensive functions. The results also underscore the adaptive significance of defensive peptides in plant evolution, highlighting their roles in responding to ecological pressures and pathogen interactions.

Substrate selectivity and unique sequence signatures in SWEET/semiSWEET homologs of four taxonomic groups: Sequence analysis and phylogenetic studies.

The recently discovered SWEET (Sugar Will Eventually be Exported Transporter) proteins are involved in the selective transport of monosaccharides and disaccharides. The prokaryotic counterparts, semiSWEETs, form dimers with each monomer forming a triple-helix transmembrane bundle (THB). The longer eukaryotic SWEETs have seven transmembrane helices with two THBs and a linker helix. Structures of semiSWEETs/SWEETs have been determined experimentally. Experimental studies revealed the role of plant SWEETs in vital physiological processes and identified residues responsible for substrate selectivity. However, SWEETs/semiSWEETs from metazoans and bacteria are not characterized. In this study, we used structure-based sequence alignment and compared more than 2000 SWEET/semiSWEETs from four different taxonomic groups. Conservation of residue/chemical property was examined at all positions. Properties of clades/subclades of phylogenetic trees from each taxonomic group were analyzed. Conservation pattern of known residues in the selectivity-filter was used to predict the substrate preference of plant SWEETs and some clusters of metazoans and bacteria. Some residues at the gating and substrate-binding regions, pore-facing positions and at the helix-helix interface are conserved across all taxonomic groups. Conservation of polar/charged residues at specific pore-facing positions, helix-helix interface and in loops seems to be unique for plant SWEETs. Overall, the number of conserved residues is less in metazoan SWEETs. Plant and metazoan SWEETs exhibit high conservation of four and three proline residues respectively in "proline tetrad." Further experimental studies can validate the predicted substrate selectivity and significance of conserved polar/charged/aromatic residues at structurally and functionally important positions of SWEETs/semiSWEETs in plants, metazoans and bacteria.

Vertical Transmission of Hepatitis C Virus Among Women With a History of Injection Opioid Use.

We evaluated vertical transmission and linkage to care in women with hepatitis C virus (HCV) and history of injection drug use employing co-localized testing and treatment. Transmission occurred in 1 of 23 infants, with mother-infant genetic distance of 1.26%. Rates for infant testing, maternal linkage, and cure were 77%, 52%, and 100%, respectively.

Complete mitochondrial genome and phylogenetic analysis of Dollfustrema vaneyi (Trematoda: Bucephalidae).

BACKGROUND: The Bucephalidae is a large family of digenean trematodes but most previous analyses of its phylogenetic position have relied on a single mitochondrial gene or morphological features. Mitochondrial genomes (mitogenomes) remain unavailable for the entire family. To address this, we sequenced the complete mitogenome of Dollfustrema vaneyi and analyzed the phylogenetic relationships with other trematodes. RESULTS: The circular genome of Dollfustrema vaneyi spanned 14,959 bp and contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a major non-coding region. We used concatenated amino acid and nucleotide sequences of all 36 genes for phylogenetic analyses, conducted using MrBayes, IQ-TREE and PhyloBayes. We identified pronounced topological instability across different analyses. The addition of recently sequenced two mitogenomes for the Aspidogastrea subclass along with the use of a site-heterogeneous model stabilized the topology, particularly the positions of Azygiidae and Bucephalidae. The stabilized results indicated that Azygiidae was the closest lineage to Bucephalidae in the available dataset, and together, they clustered at the base of the Plagiorchiida. CONCLUSIONS: Our study provides the first comprehensive description and annotation of the mitochondrial genome for the Bucephalidae family. The results indicate a close phylogenetic relationship between Azygiidae and Bucephalidae, and reveal their basal placement within the order Plagiorchiida. Furthermore, the inclusion of Aspidogastrea mitogenomes and the site-heterogeneous model significantly improved the topological stability. These data will provide key molecular resources for future taxonomic and phylogenetic studies of the family Bucephalidae.

Comparative analysis of complete chloroplast genomes of Synotis species (Asteraceae, Senecioneae) for identification and phylogenetic analysis.

BACKGROUND: The Synotis (C. B. Clarke) C. Jeffrey & Y. L. Chen is an ecologically important genus of the tribe Senecioneae, family Asteraceae. Because most species of the genus bear similar morphology, traditional morphological identification methods are very difficult to discriminate them. Therefore, it is essential to develop a reliable and effective identification method for Synotis species. In this study, the complete chloroplast (cp.) genomes of four Synotis species, S. cavaleriei (H.Lév.) C. Jeffrey & Y.L. Chen, S. duclouxii (Dunn) C. Jeffrey & Y.L. Chen, S. nagensium (C.B. Clarke) C. Jeffrey & Y.L. Chen and S. erythropappa (Bureau & Franch.) C. Jeffrey & Y. L. Chen had been sequenced using next-generation sequencing technology and reported here. RESULTS: These four cp. genomes exhibited a typical quadripartite structure and contained the large single-copy regions (LSC, 83,288 to 83,399 bp), the small single-copy regions (SSC, 18,262 to 18,287 bp), and the inverted repeat regions (IR, 24,837 to 24,842 bp). Each of the four cp. genomes encoded 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (ycf1 and rps19). The highly variable regions (trnC-GCA-petN, ccsA-psaC, trnE-UUC-rpoB, ycf1, ccsA and petN) may be used as potential molecular barcodes. The complete cp. genomes sequence of Synotis could be used as the potentially effective super-barcode to accurately identify Synotis species. Phylogenetic analysis demonstrated that the four Synotis species were clustered into a monophyletic group, and they were closed to the Senecio, Crassocephalum and Dendrosenecio in tribe Senecioneae. CONCLUSIONS: This study will be useful for further species identification, evolution, genetic diversity and phylogenetic studies within this genus Synotis and the tribe Senecioneae.

Analysis of the recombination and evolution of the new type mutant pseudorabies virus XJ5 in China.

Pseudorabies have caused enormous economic losses in China's pig industry and have recurred on many large pig farms since late 2011. The disease is caused by highly pathogenic, antigenic variant pseudorabies virus (vPRV) strains. Our laboratory isolated a pseudorabies virus in 2015 and named it XJ5. The pathogenic ability of this mutant strain was much stronger than that of the original isolate. After we sequenced its whole genome (GenBank accession number: OP512542), we found that its overall structure was not greatly changed compared with that of the previous strain Ea (KX423960.1). The whole genome alignment showed that XJ5 had a strong genetic relationship with the strains isolated in China after 2012 reported in GenBank. Based on the isolation time of XJ5 and the mutation and recombination analysis of programs, we found that the whole genome homology of XJ5 and other strains with Chinese isolates was greater than 95%, while the homology with strains outside Asia was less than 94%, which indicated that there may be some recombination and mutation patterns. We found that virulent PRV isolates emerged successively in China in 2011 and formed two different evolutionary clades from foreign isolates. At the same time, this may be due to improper immunization and the presence of wild strains in the field, and recent reports have confirmed that Bartha vaccine strains recombine with wild strains to obtain new pathogenic strains. We performed genetic evolution analysis of XJ5 isolated and sequenced in our laboratory to trace its possible mutations and recombination. We found that XJ5 may be the result of natural mutation of a virus in a branch of mutant strains widely existing in China.

Solanum aculeatissimum and Solanum torvum chloroplast genome sequences: a comparative analysis with other Solanum chloroplast genomes.

BACKGROUND: Solanum aculeatissimum and Solanum torvum belong to the Solanum species, and they are essential plants known for their high resistance to diseases and adverse conditions. They are frequently used as rootstocks for grafting and are often crossbred with other Solanum species to leverage their resistance traits. However, the phylogenetic relationship between S. aculeatissimum and S. torvum within the Solanum genus remains unclear. Therefore, this paper aims to sequence the complete chloroplast genomes of S. aculeatissimum and S. torvum and analyze them in comparison with 29 other previously published chloroplast genomes of Solanum species. RESULTS: We observed that the chloroplast genomes of S. aculeatissimum and S. torvum possess typical tetrameric structures, consisting of one Large Single Copy (LSC) region, two reverse-symmetric Inverted Repeats (IRs), and one Small Single Copy (SSC) region. The total length of these chloroplast genomes ranged from 154,942 to 156,004 bp, with minimal variation. The highest GC content was found in the IR region, while the lowest was in the SSC region. Regarding gene content, the total number of chloroplast genes and CDS genes remained relatively consistent, ranging from 128 to 134 and 83 to 91, respectively. Nevertheless, there was notable variability in the number of tRNA genes and rRNAs. Relative synonymous codon usage (RSCU) analysis revealed that both S. aculeatissimum and S. torvum preferred codons that utilized A and U bases. Analysis of the IR boundary regions indicated that contraction and expansion primarily occurred at the junction between SSC and IR regions. Nucleotide polymorphism analysis and structural variation analysis demonstrated that chloroplast variation in Solanum species mainly occurred in the LSC and SSC regions. Repeat sequence analysis revealed that A/T was the most frequent base pair in simple repeat sequences (SSR), while Palindromic and Forward repeats were more common in long sequence repeats (LSR), with Reverse and Complement repeats being less frequent. Phylogenetic analysis indicated that S. aculeatissimum and S. torvum belonged to the same meristem and were more closely related to Cultivated Eggplant. CONCLUSION: These findings enhance our comprehension of chloroplast genomes within the Solanum genus, offering valuable insights for plant classification, evolutionary studies, and potential molecular markers for species identification.

Investigating pigeon circovirus infection in a pigeon farm: molecular detection, phylogenetic analysis and complete genome analysis.

BACKGROUND: Pigeon circovirus infections in pigeons (Columba livia domestica) have been reported worldwide. Pigeons should be PiCV-free when utilized as qualified experimental animals. However, pigeons can be freely purchased as experimental animals without any clear guidelines to follow. Herein, we investigated the status quo of PiCV infections on a pigeon farm in Beijing, China, which provides pigeons for experimental use. RESULTS: PiCV infection was verified in at least three types of tissues in all forty pigeons tested. A total of 29 full-length genomes were obtained and deposited in GenBank. The whole genome sequence comparison among the 29 identified PiCV strains revealed nucleotide homologies of 85.8-100%, and these sequences exhibited nucleotide homologies of 82.7-98.9% as compared with those of the reference sequences. The cap gene displayed genetic diversity, with a wide range of amino acid homologies ranging from 64.5% to 100%. Phylogenetic analysis of the 29 full-genome sequences revealed that the PiCV strains in this study could be further divided into four clades: A (17.2%), B (10.4%), C (37.9%) and D (34.5%). Thirteen recombination events were also detected in 18 out of the 29 PiCV genomes obtained in this study. Phylogenetic research using the rep and cap genes verified the recombination events, which occurred between clades A/F, A/B, C/D, and B/D among the 18 PiCV strains studied. CONCLUSIONS: In conclusion, PiCV infection, which is highly genetically varied, is extremely widespread on pigeon farms in Beijing. These findings indicate that if pigeons are to be used as experimental animals, it is necessary to evaluate the impact of PiCV infection on the results.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant.

BACKGROUND: Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported. RESULTS: In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted. CONCLUSIONS: This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome.

BACKGROUND: Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear. RESULTS: We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species. CONCLUSIONS: These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.