Small Genome Research Articles

Molecular characterization and functional study of the NLRP3 inflammasome genes in Tetraodon nigroviridis

Nucleotide-binding oligomerization domain-like receptor protein 3 (NLRP3) inflammasome is an important inflammasome in mammals, which is of great significance to eliminate pathogens. However, the research of the NLRP3 inflammasome in teleost is limited. Tetraodon nigroviridis has the characteristics of small genome and easy feeding, which can be used as a model for the study of fish immune function. In present study, three NLRP3 inflammasome component genes (NLRP3, ASC and caspase-1) in T. nigroviridis has been cloned. Real-time fluorescence quantitative PCR showed that TnNLRP3 (T. nigroviridis NLRP3), TnASC (T. nigroviridis ASC) and Tncaspase-1 (T. nigroviridis caspase-1) mRNA in various tissues from health T. nigroviridis were highly expressed in immune-related tissues, such as spleen, gill, head kidney and intestine. After Vibrio parahemolyticus infection, the expression of TnNLRP3, TnASC and Tncaspase-1 mRNA in spleen, gill, head kidney reached a peak at 24 h, and the expression levels of these genes in intestine were the highest at 48 h. After the transfection of TnASC-pAcGFP-N1 monomer GFP plasmid into cos-7 cells, ASC specks, the activation marker of NLRP3 inflammasome, were observed. Bimolecular fluorescence complementarity and fluorescence colocation experiment showed that TnASC and Tncaspase-1 of TnNLRP3 inflammasome were co-located near the cell nucleus, and potentially interacted with each other. NLRP3 inflammasome inducer nigericin and agonist ATP could significantly induce the expression of TnNLRP3, TnASC and Tncaspase-1 mRNA, and activation of NLRP3 inflammasome could promote the generation of mature TnIL-1β (T. nigroviridis IL-1β). These results uncover that T. nigroviridis NLRP3 inflammasome could participate in the antibacterial immune response and the generation of mature TnIL-1β after activation.

Identification of SNPs and Candidate Genes Associated with Salt Stress in Two Korean Sorghum Cultivars and Understanding Selection Pressures in the Breeding Process

One of the abiotic stresses, salt stress, has an impact on the production and development of crops around the world. Sorghum is a functional genomics model crop of C4 plants due to its small genome size, and it is suitable for providing a clue to the mechanism associated with salt tolerance at the transcriptomic level. However, the mechanism of salt-related genes in sorghum has not been well described. RNA sequencing, using QuantSeq, was performed on two Korean cultivars, ‘Sodamchal’ and ‘Nampungchal’, which are known to have different intensities in response to salt stress, between a control and high-salinity treatment over a different time-course. In addition, physiological responses such as the proline, anthocyanin, chlorophyll, and reducing sugar contents were evaluated under the salt-stress treatment between these two sorghum cultivars. Moreover, differentially expressed genes (DEGs) between the Nampungchal and Sodamchal cultivars were identified in their leaves and roots, respectively. Moreover, the function of DEGs was confirmed through GO classification and KEGG pathway. We also analyzed the correlation between the selection pressure with DEGs by identifying Ka/Ks of DEGs. In the breeding process, the role of positive or negative selected genes was analyzed. Therefore, a new hypothesis on selection pressure was proposed from the breeding perspective of cultivars. A comparative analysis of the two sorghum cultivars provides candidate genes involved in the salt-stress response and may offer a better understanding of the salt-tolerance mechanism in sorghum.

Plant genome size modulates grassland community responses to multi-nutrient additions.

Grassland ecosystems cover c. 40% of global land area and contain c. 40% of soil organic carbon. Understanding the effects of adding nutrients to grasslands is essential because they provide much of our food, support diverse ecosystem services and harbor rich biodiversity. Using the meadow steppe (grassland) study site of Inner Mongolia, we manipulated seven key nutrients and a cocktail of micronutrients to examine their effects on grassland biomass productivity and diversity. The results, explained in structural equation models, link two previously disparate hypotheses in grassland ecology: (1) the light asymmetry competition hypothesis and (2) the genome size-nutrient interaction hypothesis. We show that aboveground net primary productivity increases predominantly from species with large genome sizes with the addition of nitrogen, and nitrogen plus phosphorus. This drives an asymmetric competition for light, causing a decline in species richness mainly in species with small genome sizes. This dynamic is likely to be caused by the nutrient demands of the nucleus and/or the scaling effects of nuclear size on cell size which impact water use efficiency. The model will help inform the best management approaches to reverse the rapid and unprecedented degradation of grasslands globally.

Rapid transcriptomic and physiological changes in the freshwater pennate diatom Mayamaea pseudoterrestris in response to copper exposure.

Diatoms function as major primary producers, accumulating large amounts of biomass in most aquatic environments. Given their rapid responses to changes in environmental conditions, diatoms are used for the biological monitoring of water quality and for performing ecotoxicological tests in aquatic ecosystems. However, the molecular basis for their toxicity to chemical compounds remains largely unknown. Here, we sequenced the genome of a freshwater diatom, Mayamaea pseudoterrestris NIES-4280, which has been proposed as an alternative strain of Navicula pelliculosa UTEX 664 for performing the Organisation for Economic Co-operation and Development ecotoxicological test. This study shows that M. pseudoterrestris has a small genome and carries the lowest number of genes among freshwater diatoms. The gene content of M. pseudoterrestris is similar to that of the model marine diatom, Phaeodactylum tricornutum. Genes related to cell motility, polysaccharide metabolism, oxidative stress alleviation, intracellular calcium signalling, and reactive compound detoxification showed rapid changes in their expression patterns in response to copper exposure. Active gliding motility was observed in response to copper addition, and copper exposure decreased intracellular calcium concentration. These findings enhance our understanding of the environmental adaptation of diatoms, and elucidate the molecular basis of toxicity of chemical compounds in algae.

Life cycle and functional genomics of the unicellular red alga Galdieria for elucidating algal and plant evolution and industrial use

Sexual reproduction is widespread in eukaryotes; however, only asexual reproduction has been observed in unicellular red algae, including Galdieria, which branched early in Archaeplastida. Galdieria possesses a small genome; it is polyextremophile, grows either photoautotrophically, mixotrophically, or heterotrophically, and is being developed as an industrial source of vitamins and pigments because of its high biomass productivity. Here, we show that Galdieria exhibits a sexual life cycle, alternating between cell-walled diploid and cell wall-less haploid, and that both phases can proliferate asexually. The haploid can move over surfaces and undergo self-diploidization or generate heterozygous diploids through mating. Further, we prepared the whole genome and a comparative transcriptome dataset between the diploid and haploid and developed genetic tools for the stable gene expression, gene disruption, and selectable marker recycling system using the cell wall-less haploid. The BELL/KNOX and MADS-box transcription factors, which function in haploid-to-diploid transition and development in plants, are specifically expressed in the haploid and diploid, respectively, and are involved in the haploid-to-diploid transition in Galdieria, providing information on the missing link of the sexual life cycle evolution in Archaeplastida. Four actin genes are differently involved in motility of the haploid and cytokinesis in the diploid, both of which are myosin independent and likely reflect ancestral roles of actin. We have also generated photosynthesis-deficient mutants, such as blue-colored cells, which were depleted in chlorophyll and carotenoids, for industrial pigment production. These features of Galdieria facilitate the understanding of the evolution of algae and plants and the industrial use of microalgae.

Achieving a more robust antiviral RNAi via subverting a viral virulence protein

As submicroscopic organisms with small genomes, viruses can only replicate inside living cells. In fact, viruses can infect all forms of life, from animals and plants to microorganisms including fungi, bacteria, and archaea. Besides well-known human diseases such as human acquired immune deficiency syndrome, common cold, influenza, hepatitis, and coronavirus disease 2019, viral pathogens also cause more than $30 billion crop yield losses annually worldwide (Chauhan et al., 2019).

Genome-wide identification of foxtail millet's TRX family and a functional analysis of SiNRX1 in response to drought and salt stresses in transgenic Arabidopsis.

Thioredoxins (TRXs) are small-molecule proteins with redox activity that play very important roles in the growth, development, and stress resistance of plants. Foxtail millet (Setaria italica) gradually became a model crop for stress resistance research because of its advantages such as its resistance to sterility and its small genome. To date, the thioredoxin (TRX) family has been identified in Arabidopsis thaliana, rice and wheat. However, studies of the TRX family in foxtail millet have not been reported, and the biological function of this family remains unclear. In this study, 35 SiTRX genes were identified in the whole genome of foxtail millet through bioinformatic analysis. According to phylogenetic analysis, 35 SiTRXs can be divided into 13 types. The chromosome distribution, gene structure, cis-elements and conserved protein motifs of 35 SiTRXs were characterized. Three nucleoredoxin (NRX) members were further identified by a structural analysis of TRX family members. The expression patterns of foxtail millet’s SiNRX members under abiotic stresses showed that they have different stress-response patterns. In addition, subcellular localization revealed that SiNRXs were localized to the nucleus, cytoplasm and membrane. Further studies demonstrated that the overexpression of SiNRX1 enhanced Arabidopsis’ tolerance to drought and salt stresses, resulting in a higher survival rate and better growth performance. Moreover, the expression levels of several known stress-related genes were generally higher in overexpressed lines than in the wild-type. Thus, this study provides a general picture of the TRX family in foxtail millet and lay a foundation for further research on the mechanism of the action of TRX proteins on abiotic stresses.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome

BackgroundB chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome.ResultsThrough small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome.ConclusionsDue to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

Jatropha curcas L. as a Plant Model for Studies on Vegetative Propagation of Native Forest Plants.

Highlights Jatropha curcas L. has great potential to be used as a model plant in several studies involving native forest species.The immersion in the 2,4-D solution accelerated the emission of primary roots in hardwood cuttings.Studies on vegetative propagation of native species can use Jatropha curcas L. species as a model for obtaining important information in a short time and reducing labor costs.The immersion of cuttings of native species in solutions with low concentrations of 2,4-D can favor the rooting process and vegetative propagation. Even though it is a forest native plant, there are already several studies evaluating the small genome of Jatropha curcas L., which belongs to the Euphorbiaceae family, and may be an excellent representative model for the other plants from the same family. Jatropha curcas L. plant has fast growth, precocity, and great adaptability, facilitating silvicultural studies, allowing important information to be obtained quickly, and reducing labor costs. This information justifies the use of the species as a model plant in studies involving the reproduction of native plants. This study aimed to evaluate the possibility of using Jatropha curcas L. as a model plant for studies involving native forest plants and establish possible recommendations for the vegetative propagation of the species using hardwood cuttings. The information collected can be helpful to other native forest plant species, similar to Jatropha curcas L. To this end, the effects of hardwood cutting length (10, 20, and 30 cm) and the part of the hardwood cuttings (basal, middle, and apex) were evaluated. Moreover, the influence of immersing the hardwood cuttings in solutions containing micronutrients (boron or zinc) or plant regulators (2,4-D, GA3) and a biostimulant composed of kinetin (0.09 g L−1), gibberellic acid (0.05 g L−1), and 4-indole-3-butyric acid (0.05 g L−1). The experiments were carried out in duplicates. In one duplicate, sand was used as the substrate, and rooting evaluations were made 77 days after planting. In another duplicate, a substrate composed of 50% soil, 40% poultry litter, and 10% sand was used, and the evaluations of the saplings were performed 120 days after planting. The GA3 solutions inhibited the roots’ and sprouts’ emissions, while immersion in 2,4-D solution increased the number of primary roots at 77 days after planting. The hardwood cuttings from the basal part of the branch had the best results for producing saplings.

Context dependent prediction in DNA sequence using neural networks.

One way to better understand the structure in DNA is by learning to predict the sequence. Here, we trained a model to predict the missing base at any given position, given its left and right flanking contexts. Our best-performing model was a neural network that obtained an accuracy close to 54% on the human genome, which is 2% points better than modelling the data using a Markov model. In likelihood-ratio tests, the neural network performed significantly better than any of the alternative models by a large margin. We report on where the accuracy was obtained, first observing that the performance appeared to be uniform over the chromosomes. The models performed best in repetitive sequences, as expected, although their performance far from random in the more difficult coding sections, the proportions being ~70:40%. We further explored the sources of the accuracy, Fourier transforming the predictions revealed weak but clear periodic signals. In the human genome the characteristic periods hinted at connections to nucleosome positioning. We found similar periodic signals in GC/AT content in the human genome, which to the best of our knowledge have not been reported before. On other large genomes similarly high accuracy was found, while lower predictive accuracy was observed on smaller genomes. Only in the mouse genome did we see periodic signals in the same range as in the human genome, though weaker and of a different type. This indicates that the sources of these signals are other or more than nucleosome arrangement. Interestingly, applying a model trained on the mouse genome to the human genome resulted in a performance far below that of the human model, except in the difficult coding regions. Despite the clear outcomes of the likelihood-ratio tests, there is currently a limited superiority of the neural network methods over the Markov model. We expect, however, that there is great potential for better modelling DNA using different neural network architectures.

Describe of Arabidopsis Thaliana Plant as a Model Plant in Biotechnology

Arabidopsis thaliana plant has appeared as a promising model because it is self-fertilizing, which makes it preserve its characteristics and can produce thousands of seeds, which makes it ideal for conducting mutagenesis experiments. The identification of its selected traits is more clear because it has a relatively small genome. Also, the organization of the genome in it is simple and is especially suitable in genetic and biology experiments, as these plants represent small-sized and fast-rotating factories. The current study attempted to highlight the most important aspects and applications related to the use of this plant in various fields of biotechnology and to clarify its most important advantages.

Genomic characterization of two metagenome-assembled genomes of Tropheryma whipplei from China.

Whipple’s disease is a rare chronic systemic disease that affects almost any organ system of the body caused by the intracellular bacterium Tropheryma whipplei, which is found ubiquitously in the environment. Sequencing of the T. whipplei genome has revealed that it has a reduced genome (0.93 Mbp), a characteristic shared with other intracellular bacteria. Until our research started, 19 T. whipplei strains had been sequenced from cultures originated in France, Canada, and Germany. The genome of T. whipplei bacterium has not been studied in Asia yet. Here, two metagenome-assembled genomes (MAGs) of T. whipplei from China were reconstructed through metagenomic next-generation sequencing (mNGS) and genome binning. We also provided genomic insights into the geographical role and genomic features by analyzing the whole genome. The whole-genome phylogenetic tree was constructed based on single-nucleotide polymorphism (SNP) distance calculations and then grouped by distance similarity. The phylogenetic tree shows inconsistencies with geographic origins, thus suggesting that the variations in geographical origins cannot explain the phylogenetic relationships among the 21 T. whipplei strains. The two Chinese strains were closely related to each other, and also found to be related to strains from Germany (T. whipplei TW08/27) and France (T. whipplei Bcu26 and T. whipplei Neuro1). Furthermore, the Average Nucleotide Identity (ANI) matrix also showed no association between geographic origins and genomic similarities. The pan-genome analysis revealed that T. whipplei has a closed pan-genome composed of big core-genomes and small accessory genomes, like other intracellular bacteria. By examining the genotypes of the sequenced strains, all 21 T. whipplei strains were found to be resistant to fluoroquinolones, due to the genetic mutations in genes gyrA, gyrB, parC, and parE. The 21 T. Whipplei strains shared the same virulence factors, except for the alpC gene, which existed in 7 out of the 21 T. whipplei strains. When comparing 21 entire T. whipplei pan-genomes from various nations, it was discovered that the bacterium also possessed a closed genome, which was a trait shared by intracellular pathogens.

Reference genome of the Woolly Sculpin, Clinocottus analis.

Sculpins (Family Cottidae) are generally cold-temperate intertidal reef fishes most commonly found in the North Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Woolly Sculpin, Clinocottus analis, to establish a genomic model for understanding phylogeographic structure of inshore marine taxa along the California coast. These patterns, in turn, should further inform the design of marine protected areas using dispersal models based on genomic data. The small genome of C. analis is typical of marine fishes at less than 1 Gb (genome size = 538 Mb), and our assembly is near-chromosome level (contig N50 = 9.1 Mb, scaffold N50 = 21 Mb, BUSCO completeness = 97.9%). Within the context of the CCGP, the Woolly Sculpin genome will be used as a reference for future whole-genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of marine protected areas across the state.

Rice grain yield and quality improvement via CRISPR/Cas9 system: an updated review

Rice (Oryza sativa L.) is an important staple food crop worldwide. To meet the growing nutritional requirements of the increasing population in the face of climate change, qualitative and quantitative traits of rice need to be improved. During recent years, genome editing has played a great role in the development of superior varieties of grain crops. Genome editing and speed breeding have improved the accuracy and pace of rice breeding. New breeding technologies including genome editing have been established in rice, expanding the potential for crop improvement. Over a decade, site-directed mutagenesis tools like Zinc Finger Nucleases (ZFN), Transcriptional activator-like Effector Nucleases (TALENs), and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) System were used and have played a great role in rice yield and quality enhancement. In addition, most recently other genome editing techniques like prime editing and base editors have also been used for efficient genome editing in rice. Since rice is an excellent model system for functional studies due to its small genome and close synthetic relationships with other cereal crops, new genome-editing technologies continue to be developed for use in rice. Genomic alteration employing genome editing technologies (GETs) like CRISPR/Cas9 for reverse genetics has opened new avenues in agricultural sciences such as rice yield and grain quality improvement. Currently, CRISPR/Cas9 technology is widely used by researchers for genome editing to achieve the desired biological objectives, because of its simple targeting, easy-to-design, cost-effective, and versatile tool for precise and efficient plant genome editing. Over the past few years many genes related to rice grain quality and yield enhancement have been successfully edited via CRISPR/Cas9 technology method to cater to the growing demand for food worldwide. The effectiveness of these methods is being verified by the researchers and crop scientists worldwide. In this review we focus on genome-editing tools for rice improvement to address the progress made and provide examples of genome editing in rice. We also discuss safety concerns and methods for obtaining transgene-free crops.

Diversity patterns and speciation processes in a two-island system with continuous migration.

Geographic isolation is a central mechanism of speciation, but perfect isolation of populations is rare. Although speciation can be hindered if gene flow is large, intermediate levels of migration can enhance speciation by introducing genetic novelty in the semi-isolated populations or founding small communities of migrants. Here, we consider a two-island neutral model of speciation with continuous migration and study diversity patterns as a function of the migration probability, population size, and number of genes involved in reproductive isolation (dubbed as genome size). For small genomes, low levels of migration induce speciation on the islands that otherwise would not occur. Diversity, however, drops sharply to a single species inhabiting both islands as the migration probability increases. For large genomes, sympatric speciation occurs even when the islands are strictly isolated. Then species richness per island increases with the probability of migration, but the total number of species decreases as they become cosmopolitan. For each genome size, there is an optimal migration intensity for each population size that maximizes the number of species. We discuss the observed modes of speciation induced by migration and how they increase species richness in the insular system while promoting asymmetry between the islands and hinderingendemism.

Intron size minimisation in teleosts

BackgroundSpliceosomal introns are parts of primary transcripts that are removed by RNA splicing. Although introns apparently do not contribute to the function of the mature transcript, in vertebrates they comprise the majority of the transcribed region increasing the metabolic cost of transcription. The persistence of long introns across evolutionary time suggests functional roles that can offset this metabolic cost. The teleosts comprise one of the largest vertebrate clades. They have unusually compact and variable genome sizes and provide a suitable system for analysing intron evolution.ResultsWe have analysed intron lengths in 172 vertebrate genomes and show that teleost intron lengths are relatively short, highly variable and bimodally distributed. Introns that were long in teleosts were also found to be long in mammals and were more likely to be found in regulatory genes and to contain conserved sequences. Our results argue that intron length has decreased in parallel in a non-random manner throughout teleost evolution and represent a deviation from the ancestral state.ConclusionOur observations indicate an accelerated rate of intron size evolution in the teleosts and that teleost introns can be divided into two classes by their length. Teleost intron sizes have evolved primarily as a side-effect of genome size evolution and small genomes are dominated by short introns (<256 base pairs). However, a non-random subset of introns has resisted this process across the teleosts and these are more likely have functional roles in all vertebrate clades.

Applications of In Vitro Tissue Culture Technologies in Breeding and Genetic Improvement of Wheat.

Sources of new genetic variability have been limited to existing germplasm in the past. Wheat has been studied extensively for various agronomic traits located throughout the genome. The large size of the chromosomes and the ability of its polyploid genome to tolerate the addition or loss of chromosomes facilitated rapid progress in the early study of wheat genetics using cytogenetic techniques. At the same time, its large genome size has limited the progress in genetic characterization studies focused on diploid species, with a small genome and genetic engineering procedures already developed. Today, the genetic transformation and gene editing procedures offer attractive alternatives to conventional techniques for breeding wheat because they allow one or more of the genes to be introduced or altered into an elite cultivar without affecting its genetic background. Recently, significant advances have been made in regenerating various plant tissues, providing the essential basis for regenerating transgenic plants. In addition, Agrobacterium-mediated, biolistic, and in planta particle bombardment (iPB) gene delivery procedures have been developed for wheat transformation and advanced transgenic wheat development. As a result, several useful genes are now available that have been transferred or would be helpful to be transferred to wheat in addition to the current traditional effort to improve trait values, such as resistance to abiotic and biotic factors, grain quality, and plant architecture. Furthermore, the in planta genome editing method will significantly contribute to the social implementation of genome-edited crops to innovate the breeding pipeline and leverage unique climate adaptations.

EfGD: the Erianthus fulvus genome database.

Erianthus fulvus (TaxID: 154759) is a valuable germplasm resource in sugarcane breeding and research and has excellent agronomic traits, such as drought resistance, cold resistance, barren tolerance and high brix. With a stable chromosome number (2n = 20) and a small genome (0.9 Gb), it is an ideal candidate for research on sugarcane. Next-generation sequencing technology has enabled a growing number of studies to focus on genomics. Due to the large amount of omics data available, a centralized platform is necessary for ensuring the consistency, independence and maintainability of these large-scale datasets through storage, analysis and integration. Here, we present a comprehensive database for the E. fulvus genome, EfGD. By using the new high-quality reference genome and its annotations, the EfGD provides the largest whole-genome sequencing reference dataset for E. fulvus, which archives 27 165 protein-coding genes and 55 564 488 SNPs from 202 newly resequenced genomes. Furthermore, we created a user-friendly graphical interface for visualizing genomic diversity, population structure and evolution and provided other tools on an open platform. Database URL: https://efgenome.ynau.edu.cn

Genome dynamics, codon usage patterns and influencing factors in Aeromonas hydrophila phages

In the present study, genome characteristics and codon usage patterns of 44 Aeromonas hydrophila phages were studied. Phage genomes varied from 30.8 to 262.0 kb with mean±SD and median values of 111.3 ± 81.4 kb and 79.4 kb, respectively. Though the great variation in phage GC contents (35.1–62.2%) was observed, GC contents of all phages (except two phages) were significantly less than the GC content (62.4 ± 5.6%) of the host. The effective number of codons (ENC) values of phage genes ranged from 27.7 to 61 with a mean±SD value of 47.4 ± 6.8. Out of a total 5773 phage genes, 207 (3.6%), 3,528 (61.1%) and 2,012 (34.9%) genes had strong (ENC < 35), moderate (35 < ENC < 50) and low (ENC ≥ 50) codon usage bias, respectively. During relative synonymous codon usage (RSCU) analysis, shared usage of preferred codons was also observed between the phages and host. During codon adaptation index (CAI) analysis, 1028 (17.8%) phage genes showed significant adaptation towards the host. Among these genes, 797 (78.0%) genes encoded hypothetical proteins or proteins of unknown function; whereas 118 (12%) genes encoded the phage structural and packaging proteins. Segregation of ENC, RSCU and CAI analysis results based on genome size also indicated that codon usage bias was more prominent in phages with small genomes. Correlation, neutrality and GC3 versus ENC analyzes indicated a more dominant role of natural selection in shaping the codon usage patterns of A. hydrophila phages. The findings of the current study could be useful from evolutionary and host-pathogen interaction perspectives leading to efficient utilization of phages for therapeutic and other applications.

Multiple Viral Protein Genome-Linked Proteins Compensate for Viral Translation in a Positive-Sense Single-Stranded RNA Virus Infection.

Viral protein genome-linked (VPg) protein plays an essential role in protein-primed replication of plus-stranded RNA viruses. VPg is covalently linked to the 5' end of the viral RNA genome via a phosphodiester bond typically at a conserved amino acid. Whereas most viruses have a single VPg, some viruses have multiple VPgs that are proposed to have redundant yet undefined roles in viral replication. Here, we use cricket paralysis virus (CrPV), a dicistrovirus that has four nonidentical copies of VPg, as a model to characterize the role of VPg copies in infection. Dicistroviruses contain two main open reading frames (ORFs) that are driven by distinct internal ribosome entry sites (IRESs). We systematically generated single and combinatorial deletions and mutations of VPg1 to VPg4 within the CrPV infectious clone and monitored viral yield in Drosophila S2 cells. Deletion of one to three VPg copies progressively decreased viral yield and delayed viral replication, suggesting a threshold number of VPgs for productive infection. Mass spectrometry analysis of CrPV VPg-linked RNAs revealed viral RNA linkage to either a serine or threonine in VPg, mutations of which in all VPgs attenuated infection. Mutating serine 4 in a single VPg abolished viral infection, indicating a dominant negative effect. Using viral minigenome reporters that monitor dicistrovirus 5' untranslated (UTR) and IRES translation revealed a relationship between VPg copy number and the ratio of distinct IRES translation activities. We uncovered a novel viral strategy whereby VPg copies in dicistrovirus genomes compensate for the relative IRES translation efficiencies to promote infection. IMPORTANCE Genetic duplication is exceedingly rare in small RNA viral genomes, as there is selective pressure to prevent RNA genomes from expanding. However, some small RNA viruses encode multiple copies of a viral protein, most notably an unusual viral protein that is linked to the viral RNA genome. Here, we investigate a family of viruses that contains multiple viral protein genome-linked proteins and reveal a novel viral strategy whereby viral protein copy number counterbalances differences in viral protein synthesis mechanisms.