How do synonymous mutations affect fitness?

CRISPR-Cas9-assisted genome editing in E. coli elevates the frequency of unintended mutations

Cas-assisted lambda Red recombineering techniques have rapidly become a mainstay of bacterial genome editing. Such techniques have been used to construct both individual mutants and massive libraries to assess the effects of genomic changes. We have found that a commonly used Cas9-assisted editing method results in unintended mutations elsewhere in the genome in 26% of edited clones. The unintended mutations are frequently found over 200 kb from the intended edit site and even over 10 kb from potential off-target sites. We attribute the high frequency of unintended mutations to error-prone polymerases expressed in response to dsDNA breaks introduced at the edit site. Most unintended mutations occur in regulatory or coding regions and thus may have phenotypic effects. Our findings highlight the risks associated with genome editing techniques involving dsDNA breaks in E. coli and likely other bacteria and emphasize the importance of sequencing the genomes of edited cells to ensure the absence of unintended mutations.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Modulating Gene Expression within a Microbiome Based on Computational Models

Simple Summary

Development of computational biology methodologies has provided comprehensive understanding of the complexity of microbiomes, and the extensive ways in which they influence their environment. This has awakened a new research goal, aiming to not only understand the mechanisms in which microbiomes function, but to actively modulate and engineer them for various purposes. However, current microbiome engineering techniques are usually manually tailored for a specific system and neglect the different interactions between the new genetic information and the bacterial population, turning a blind eye to processes such as horizontal gene transfer, mutations, and other genetic alterations. In this work, we developed a generic computational method to automatically tune the expression of heterologous genes within a microbiome according to given preferences, to allow the functionality of the engineering process to propagate in longer periods of time. This goal was achieved by treating each part of the gene individually and considering long term fitness effects on the environment, providing computational and experimental evidence for this approach.

Abstract

Recent research in the field of bioinformatics and molecular biology has revealed the immense complexity and uniqueness of microbiomes, while also showcasing the impact of the symbiosis between a microbiome and its host or environment. A core property influencing this process is horizontal gene transfer between members of the bacterial community used to maintain genetic variation. The essential effect of this mechanism is the exposure of genetic information to a wide array of members of the community, creating an additional “layer” of information in the microbiome named the “plasmidome”. From an engineering perspective, introduction of genetic information to an environment must be facilitated into chosen species which will be able to carry out the desired effect instead of competing and inhibiting it. Moreover, this process of information transfer imposes concerns for the biosafety of genetic engineering of microbiomes as exposure of genetic information into unwanted hosts can have unprecedented ecological impacts. Current technologies are usually experimentally developed for a specific host/environment, and only deal with the transformation process itself at best, ignoring the impact of horizontal gene transfer and gene-microbiome interactions that occur over larger periods of time in uncontrolled environments. The goal of this research was to design new microbiome-specific versions of engineered genetic information, providing an additional layer of compatibility to existing engineering techniques. The engineering framework is entirely computational and is agnostic to the selected microbiome or gene by reducing the problem into the following set up: microbiome species can be defined as wanted or unwanted hosts of the modification. Then, every element related to gene expression (e.g., promoters, coding regions, etc.) and regulation is individually examined and engineered by novel algorithms to provide the defined expression preferences. Additionally, the synergistic effect of the combination of engineered gene blocks facilitates robustness to random mutations that might occur over time. This method has been validated using both computational and experimental tools, stemming from the research done in the iGEM 2021 competition, by the TAU group.

Code inside the codon: The role of synonymous mutations in regulating splicing machinery and its impact on disease

In eukaryotes, precise pre-mRNA processing, including alternative splicing, is essential to carry out the intricate protein translation process. Both point mutations (that alter the translated protein sequence) and synonymous mutations (that do not alter the translated protein sequence) are capable of affecting the splicing process. Synonymous mutations are known to affect gene expression via altering mRNA stability, mRNA secondary structure, splicing processes, and translational kinetics. In higher eukaryotes, precise splicing is regulated by three weakly conserved cis-elements, 5' and 3' splice sites and the branch site. Many other cis-acting elements (exonic/intronic splicing enhancers and silencers) and trans-acting splicing factors (serine and arginine-rich proteins and heterogeneous nuclear ribonucleoproteins) have also been found to enhance or suppress the splicing process. The appearance of synonymous mutations in cis-acting elements can alter the splicing process by changing the binding pattern of splicing factors to exonic splicing enhancers or silencer motifs. This results in exon skipping, intron retention, and various other forms of alternative splicing, eventually leading to the emergence of a wide range of diseases. The focus of this review is to elucidate the role of synonymous mutations and their impact on abnormal splicing mechanisms. Further, this study highlights the function of synonymous mutation in mediating abnormal splicing in cancer and development of X-linked, and autosomal inherited diseases.

Synonymous Variants: Necessary Nuance in our Understanding of Cancer Drivers and Treatment Outcomes

Once called "silent mutations" and assumed to have no effect on protein structure and function, synonymous variants are now recognized to be drivers for some cancers. There have been significant advances in our understanding of the numerous mechanisms by which synonymous single nucleotide variants (sSNVs) can affect protein structure and function by affecting pre-mRNA splicing, mRNA expression, stability, folding, miRNA binding, translation kinetics, and co-translational folding. This review highlights the need for considering sSNVs in cancer biology to gain a better understanding of the genetic determinants of human cancers and to improve their diagnosis and treatment. We surveyed the literature for reports of sSNVs in cancer and found numerous studies on the consequences of sSNVs on gene function with supporting in vitro evidence. We also found reports of sSNVs that have statistically significant associations with specific cancer types but for which in vitro studies are lacking to support the reported associations. Additionally, we found reports of germline and somatic sSNVs that were observed in numerous clinical studies and for which in silico analysis predicts possible effects on gene function. We provide a review of these investigations and discuss necessary future studies to elucidate the mechanisms by which sSNVs disrupt protein function and are play a role in tumorigeneses, cancer progression, and treatment efficacy. As splicing dysregulation is one of the most well recognized mechanisms by which sSNVs impact protein function, we also include our own in silico analysis for predicting which sSNVs may disrupt pre-mRNA splicing.

First Report of Distinct Bamboo mosaic virus (BaMV) Isolates Infecting Bambusa funghomii in Vietnam and the Identification of a Highly Variable Region in the BaMV Genome

New isolates of the Bamboo mosaic virus (BaMV) were identified in Bambusa funghomii bamboo in Vietnam. Sequence analyses revealed that the Vietnam isolates are distinct from all known BaMV strains, sharing the highest sequence identities (about 77%) with the Yoshi isolates reported in California, USA. Unique satellite RNAs were also found to be associated with the BaMV Vietnam isolates. A possible recombination event was detected in the genome of BaMV-VN2. A highly variable region was identified in the ORF1 gene, in between the methyl transferase domain and helicase domain. These results revealed the presence of unique BaMV isolates in an additional bamboo species in one more country, Vietnam, and provided evidence in support of the possible involvement of environmental or host factors in the diversification and evolution of BaMV.

Effects of Selection at Linked Sites on Patterns of Genetic Variability

Patterns of variation and evolution at a given site in a genome can be strongly influenced by the effects of selection at genetically linked sites. In particular, the recombination rates of genomic regions correlate with their amount of within-population genetic variability, the degree to which the frequency distributions of DNA sequence variants differ from their neutral expectations, and the levels of adaptation of their functional components. We review the major population genetic processes that are thought to lead to these patterns, focusing on their effects on patterns of variability: selective sweeps, background selection, associative overdominance, and Hill–Robertson interference among deleterious mutations. We emphasize the difficulties in distinguishing among the footprints of these processes and disentangling them from the effects of purely demographic factors such as population size changes. We also discuss how interactions between selective and demographic processes can significantly affect patterns of variability within genomes.

Evolutionary Patterns of Codon Usage in Major Lineages of Porcine Reproductive and Respiratory Syndrome Virus in China

Porcine reproductive and respiratory syndrome virus (PRRSV) is economically important and characterized by its extensive variation. The codon usage patterns and their influence on viral evolution and host adaptation among different PRRSV strains remain largely unknown. Here, the codon usage of ORF5 genes from lineages 1, 3, 5, and 8, and MLV strains of type 2 PRRSV in China was analyzed. A compositional property analysis of ORF5 genes revealed that nucleotide C is most frequently used at the third position of codons, accompanied by rich GC3s. The effective number of codon (ENC) and codon pair bias (CPB) values indicate that all ORF5 genes have low codon bias and the differences in CPB scores among four lineages are almost not significant. When compared with host codon usage patterns, lineage 1 strains show higher CAI and SiD values, with a high similarity to pig, which might relate to its predominant epidemic propensity in the field. The CAI, RCDI, and SiD values of ORF5 genes from different passages of MLV JXA1R indicate no relation between attenuation and CPB or codon adaptation decrease during serial passage on non-host cells. These findings provide a novel way of understanding the PRRSV's evolution, related to viral survival, host adaptation, and virulence.

Life-stage-dependent supergene haplotype frequencies and metapopulation neutral genetic patterns of Atlantic cod, Gadus morhua, from Canada's Northern cod stock region and adjacent areas

Among highly migratory fish species, nursery areas occupied by juveniles often differ from adult habitats. To better understand the spatial dynamics of Canada's Northern cod stock, juveniles caught off the east coast of Newfoundland and Labrador were compared to adults from the same region as well as individuals from other areas in Atlantic Canada using double-digest restriction site-associated DNA sequencing-derived single nucleotide polymorphisms. A reduced proportion of homozygotes with a chromosomal inversion located in linkage group 1 (LG1) was detected between juvenile and adult samples in the Northern cod stock region, potentially indicating age-dependent habitat use or ontogenetic selection for attributes associated with the many genes located in LG1. No selectively neutral genetic differences were found between samples from the Northern cod stock; nevertheless, significant differences were found between some of these samples and cod collected from St. Pierre Bank, Bay of Fundy, Browns Bank and the southern Scotian Shelf. Clustering analysis of variants at neutral loci provided evidence for three major genetic units: (a) the Newfoundland Atlantic Coast, (b) eastern and southern Gulf of St. Lawrence and Burgeo Bank and (c) the Bay of Fundy, Browns Bank and southern Scotian Shelf. Both adaptive and neutral population structure within the Northern cod stock should be considered by managers to promote demographic rebuilding of the stock.

Defect in Mitochondrial NADH-Dehydrogenase Genes in Canine Mast Cell Tumours

Recent studies have demonstrated a significant role of mitochondrial DNA (mtDNA) defects in the pathogenesis of many human and some canine tumours. The aim of this study was to identify mutations in the ND2 and ND4 mitochondrial genes in canine mast cell tumours and determine their association with the process of neoplastic transformation and the phenotypic traits of dogs. In total, 136 gene sequences from 68 biological samples, including blood and neoplastic tissue samples from 34 dogs with diagnosed MCTs, were analysed. The study consisted in DNA sequencing of the ND2 and ND4 genes as well as bioinformatics and statistical analyses. For the first time, mutations in NADH-dehydrogenase genes were detected in dogs with MCTs. In total, 22 polymorphic loci and 19 mutations in the ND2 and ND4 genes were identified. The majority of the identified mutations were homoplasmic, and tumour heteroplasmy was detected in eight nucleotide positions in three dogs. Seven of the ND2 mutations and two of the ND4 mutations caused an amino acid change. The changes in non-synonymous protein-coding SNPs did not exert an adverse effect on proteins. A statistically significant correlation of the presence of mutations/polymorphisms with the sex, age, and size of the dogs and the tumour location was demonstrated. Polymorphisms and mutations in NADH-dehydrogenase genes, including mastocyte-specific changes, in canine mast cell tumours that had not been reported earlier in the literature were identified. Some of these changes may imply that these are the hotspot mutations in canine mast cell tumours. It cannot be excluded that the molecular changes are directly associated with the development of mast cell tumours, and further investigations are needed to verify whether they can become molecular markers of MCTs in the future.

Contribution of mRNA Splicing to Mismatch Repair Gene Sequence Variant Interpretation

Functional assays that assess mRNA splicing can be used in interpretation of the clinical significance of sequence variants, including the Lynch syndrome-associated mismatch repair (MMR) genes. The purpose of this study was to investigate the contribution of splicing assay data to the classification of MMR gene sequence variants. We assayed mRNA splicing for 24 sequence variants in MLH1, MSH2, and MSH6, including 12 missense variants that were also assessed using a cell-free in vitro MMR activity (CIMRA) assay. Multifactorial likelihood analysis was conducted for each variant, combining CIMRA outputs and clinical data where available. We collated these results with existing public data to provide a dataset of splicing assay results for a total of 671 MMR gene sequence variants (328 missense/in-frame indel), and published and unpublished repair activity measurements for 154 of these variants. There were 241 variants for which a splicing aberration was detected: 92 complete impact, 33 incomplete impact, and 116 where it was not possible to determine complete versus incomplete splicing impact. Splicing results mostly aided in the interpretation of intronic (72%) and silent (92%) variants and were the least useful for missense substitutions/in-frame indels (10%). MMR protein functional activity assays were more useful in the analysis of these exonic variants but by design they were not able to detect clinically important splicing aberrations identified by parallel mRNA assays. The development of high throughput assays that can quantitatively assess impact on mRNA transcript expression and protein function in parallel will streamline classification of MMR gene sequence variants.

Polymorphisms in JAK2 Gene are Associated with Production Traits and Mastitis Resistance in Dairy Cattle

The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in the JAK2 gene on the production and mastitis related traits in dairy cattle. Blood and milk samples were collected from 201 lactating dairy cattle of three breeds, i.e. Holstein Friesian (HF), Jersey (J) and Achai (A) and their crosses maintained at well-established dairy farms in Khyber Pakhtunkhwa, Pakistan. Generalized linear model was used to evaluate the association between genotypes and the studied traits. A DNA pool was made from randomly selected 30 samples which revealed three SNPs, i.e. SNP 1 in 5’ upstream region (G>A, rs379754157), SNP 2 in intron 15 (A>G, rs134192265), and SNP 3 in exon 20 (A>G, rs110298451) that were further validated in the population under study using SNaPshot technique. Of the three SNPs, SNP 1 did not obey Hardy-Weinberg equilibrium (P<0.05). SNP 2 and SNP 3 were found to be in strong linkage disequilibrium and allele G was highly prevalent compared to allele A in these SNPs. in SNP 1, the GG genotype was associated with significantly (P<0.01) higher SCC, whereas SNP 2 and SNP 3 were significantly (P<0.01) associated with higher lactose percentage compared to the other geno-types. The haplogroups association analysis revealed that H1H2 (GG GG AG) has significantly lower SCC than H2H2 (GG GG GG). The results infer that JAK2 could be an important candidate gene and the studied SNPs might be useful genetic markers for production and mastitis related traits.

Comprehensive codon usage analysis of porcine deltacoronavirus

Porcine deltacoronavirus (PDCoV) is a newly identified coronavirus of pigs that was first reported in Hong Kong in 2012. Since then, many PDCoV isolates have been identified worldwide. In this study, we analyzed the codon usage pattern of the S gene using complete coding sequences and complete PDCoV genomes to gain a deeper understanding of their genetic relationships and evolutionary history. We found that during evolution three groups evolved with a relatively low codon usage bias (effective number of codons (ENC) of 52). The factors driving bias were complex. However, the primary element influencing the codon bias of PDCoVs was natural selection. Our results revealed that different natural environments may have a significant impact on the genetic characteristics of the strains. In the future, more epidemiological surveys are required to examine the factors that resulted in the emergence and outbreak of this virus.

Compositional properties and codon usage of TP73 gene family

The TP73 gene is considered as one of the members of TP53 gene family and shows much homology to p53 gene. TP73 gene plays a pivotal role in cancer studies in addition to other biological functions. Codon usage bias (CUB) is the phenomenon of unequal usage of synonymous codons for an amino acid wherein some codons are more frequently used than others and it reveals the evolutionary relationship of a gene. Here, we report the pattern of codon usage in TP73 gene using various bioinformatic tools as no work was reported yet. Nucleotide composition analysis suggested that the mean nucleobase C was the highest, followed by G and the gene was GC rich. Correlation analysis between codon usage and GC3 suggested that most of the GC-ending codons showed positive correlation while most of the AT-ending codons showed negative correlation with GC3 in the coding sequences of TP73 gene variants in human. The CUB is moderate in human TP73 gene as evident from intrinsic codon deviation index (ICDI) analysis. Nature selected against two codons namely ATA (isoleucine) and AGA (arginine) in the coding sequences of TP73 gene during the course of evolution. A significant correlation (p < 0.05) was found between overall nucleotide composition and its composition at the 3rd codon position, indicating that both mutation pressure and natural selection might influence the CUB. The correlation analysis between ICDI and biochemical properties of protein suggested that variation of CUB was associated with degree of hydrophobicity and length of protein.

Accounting for Programmed Ribosomal Frameshifting in the Computation of Codon Usage Bias Indices

Experimental evidence shows that synonymous mutations can have important consequences on genetic fitness. Many organisms display codon usage bias (CUB), where synonymous codons that are translated into the same amino acid appear with distinct frequency. Within genomes, CUB is thought to arise from selection for translational efficiency and accuracy, termed the translational efficiency hypothesis (TEH). Indeed, CUB indices correlate with protein expression levels, which is widely interpreted as evidence for translational selection. However, these tests neglect -1 programmed ribosomal frameshifting (-1 PRF), an important translational disruption effect found across all organisms of the tree of life. Genes that contain -1 PRF signals should cost more to express than genes without. Thus, CUB indices that do not consider -1 PRF may overestimate genes’ true adaptation to translational efficiency and accuracy constraints. Here, we first investigate whether -1 PRF signals do indeed carry such translational cost. We then propose two corrections for CUB indices for genes containing -1 PRF signals. We retest the TEH in Saccharomyces cerevisiae under these corrections. We find that the correlation between corrected CUB index and protein expression remains intact for most levels of uniform -1 PRF efficiencies, and tends to increase when these efficiencies decline with protein expression. We conclude that the TEH is strengthened and that -1 PRF events constitute a promising and useful tool to examine the relationships between CUB and selection for translation efficiency and accuracy.

The fitness landscape of the codon space across environments

Fitness landscapes map the relationship between genotypes and fitness. However, most fitness landscape studies ignore the genetic architecture imposed by the codon table and thereby neglect the potential role of synonymous mutations. To quantify the fitness effects of synonymous mutations and their potential impact on adaptation on a fitness landscape, we use a new software based on Bayesian Monte Carlo Markov Chain methods and re-estimate selection coefficients of all possible codon mutations across 9 amino acid positions in Saccharomyces cerevisiae Hsp90 across 6 environments. We quantify the distribution of fitness effects of synonymous mutations and show that it is dominated by many mutations of small or no effect and few mutations of larger effect. We then compare the shape of the codon fitness landscape across amino acid positions and environments, and quantify how the consideration of synonymous fitness effects changes the evolutionary dynamics on these fitness landscapes. Together these results highlight a possible role of synonymous mutations in adaptation and indicate the potential mis-inference when they are neglected in fitness landscape studies.

Genetic and evolutionary analysis of emerging H3N2 canine influenza virus

The H3N2 canine influenza virus (CIV) originated from an avian species. Since its emergence, it has circulated in multiple states and has caused pandemics among dog populations; however, no comprehensive studies have explored the causes driving these ongoing cases. The study of the codon usage patterns of viruses can reveal the genetic changes required for the viruses to adapt to new hosts and the external environment. Here we performed a thorough genetic, evolutionary, and codon usage analysis. We identified three evolutionary H3N2 CIV clades from a timescaled phylogenetic tree, namely, Origin, China, and Korea/USA, by principal component analysis (PCA). Additionally, we found a low codon usage bias and that mutation pressure, natural selection, and dinucleotide abundance shape the codon usage bias of H3N2 CIVs, with natural selection being more crucial than the others. Moreover, the human codon adaptation index was similar to that of dogs (the natural host) and cats. In addition, the H3N2 CIV similarity index values were higher than those of the avian influenza virus (AIV), suggesting viral adaptation to the host. Therefore, H3N2 CIVs may pose a potential risk to public health in the future, and further epidemiologic, evolutionary, and pathogenetic studies are required.

The novel T755C mutation in <i>BMP15</i> is associated with the litter size of Iranian Afshari, Ghezel, and Shal breeds

Abstract. The present study was conducted to determine the molecular genetic variation in bone morphogenetic protein 15 (BMP15) of some selected Iranian sheep breeds and aims to provide relevant genetic information for twin-birth lambing. Therefore, a total of 44 rams and 213 ewes with single-, twin-, and triplet-birth lambing were used from Shal, Ghezel, Afshari and Lori-Bakhtiari sheep breeds. Three primer pairs were used for amplification of exons 1 and 2 of the BMP15 gene. The sequencing result of distinctive single-strand conformation polymorphism (SSCP) patterns showed segregation of the T755C mutation and the ensuing non-conservative substitution of L252P in the BMP15 propeptide. All three ewes with triplet-birth lambing and two sterile ewes were heterozygotic for this mutation. Association analysis revealed that this mutation was significantly related to the ewe's litter size, so that ewes with the CT genotype had 0.24 and 0.30 more lambs than those with the CC (p < 0.05) and TT (p < 0.01) genotypes, respectively. The mutation effect on the litter size was incomplete dominance in Afshari breed and over-dominance in Shal and Ghezel breeds. In conclusion, the findings in this experiment suggested synergistic effects of the T755C mutation and other unknown mutations in other effective genes which cause infertility and higher litter size in two different ways.

Comprehensive analysis of the codon usage patterns in the envelope glycoprotein E2 gene of the classical swine fever virus

The classical swine fever virus (CSFV), circulating worldwide, is a highly contagious virus. Since the emergence of CSFV, it has caused great economic loss in swine industry. The envelope glycoprotein E2 gene of the CSFV is an immunoprotective antigen that induces the immune system to produce neutralizing antibodies. Therefore, it is essential to study the codon usage of the E2 gene of the CSFV. In this study, 140 coding sequences of the E2 gene were analyzed. The value of effective number of codons (ENC) showed low codon usage bias in the E2 gene. Our study showed that codon usage could be described mainly by mutation pressure ENC plot analysis combined with principal component analysis (PCA) and translational selection-correlation analysis between the general average hydropathicity (Gravy) and aromaticity (Aroma), and nucleotides at the third position of codons (A3s, T3s, G3s, C3s and GC3s). Furthermore, the neutrality analysis, which explained the relationship between GC12s and GC3s, revealed that natural selection had a key role compared with mutational bias during the evolution of the E2 gene. These results lay a foundation for further research on the molecular evolution of CSFV.

A new and updated resource for codon usage tables

Background

Due to the degeneracy of the genetic code, most amino acids can be encoded by multiple synonymous codons. Synonymous codons naturally occur with different frequencies in different organisms. The choice of codons may affect protein expression, structure, and function. Recombinant gene technologies commonly take advantage of the former effect by implementing a technique termed codon optimization, in which codons are replaced with synonymous ones in order to increase protein expression. This technique relies on the accurate knowledge of codon usage frequencies. Accurately quantifying codon usage bias for different organisms is useful not only for codon optimization, but also for evolutionary and translation studies: phylogenetic relations of organisms, and host-pathogen co-evolution relationships, may be explored through their codon usage similarities. Furthermore, codon usage has been shown to affect protein structure and function through interfering with translation kinetics, and cotranslational protein folding.

Results

Despite the obvious need for accurate codon usage tables, currently available resources are either limited in scope, encompassing only organisms from specific domains of life, or greatly outdated. Taking advantage of the exponential growth of GenBank and the creation of NCBI’s RefSeq database, we have developed a new database, the High-performance Integrated Virtual Environment-Codon Usage Tables (HIVE-CUTs), to present and analyse codon usage tables for every organism with publicly available sequencing data. Compared to existing databases, this new database is more comprehensive, addresses concerns that limited the accuracy of earlier databases, and provides several new functionalities, such as the ability to view and compare codon usage between individual organisms and across taxonomical clades, through graphical representation or through commonly used indices. In addition, it is being routinely updated to keep up with the continuous flow of new data in GenBank and RefSeq.

Conclusion

Given the impact of codon usage bias on recombinant gene technologies, this database will facilitate effective development and review of recombinant drug products and will be instrumental in a wide area of biological research. The database is available at hive.biochemistry.gwu.edu/review/codon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1793-7) contains supplementary material, which is available to authorized users.

Genetic and codon usage bias analyses of polymerase genes of equine influenza virus and its relation to evolution

BACKGROUND

Equine influenza is a major health problem of equines worldwide. The polymerase genes of influenza virus have key roles in virus replication, transcription, transmission between hosts and pathogenesis. Hence, the comprehensive genetic and codon usage bias of polymerase genes of equine influenza virus (EIV) were analyzed to elucidate the genetic and evolutionary relationships in a novel perspective.

RESULTS

The group - specific consensus amino acid substitutions were identified in all polymerase genes of EIVs that led to divergence of EIVs into various clades. The consistent amino acid changes were also detected in the Florida clade 2 EIVs circulating in Europe and Asia since 2007. To study the codon usage patterns, a total of 281,324 codons of polymerase genes of EIV H3N8 isolates from 1963 to 2015 were systemically analyzed. The polymerase genes of EIVs exhibit a weak codon usage bias. The ENc-GC3s and Neutrality plots indicated that natural selection is the major influencing factor of codon usage bias, and that the impact of mutation pressure is comparatively minor. The methods for estimating host imposed translation pressure suggested that the polymerase acidic (PA) gene seems to be under less translational pressure compared to polymerase basic 1 (PB1) and polymerase basic 2 (PB2) genes. The multivariate statistical analysis of polymerase genes divided EIVs into four evolutionary diverged clusters - Pre-divergent, Eurasian, Florida sub-lineage 1 and 2.

CONCLUSIONS

Various lineage specific amino acid substitutions observed in all polymerase genes of EIVs and especially, clade 2 EIVs underwent major variations which led to the emergence of a phylogenetically distinct group of EIVs originating from Richmond/1/07. The codon usage bias was low in all the polymerase genes of EIVs that was influenced by the multiple factors such as the nucleotide compositions, mutation pressure, aromaticity and hydropathicity. However, natural selection was the major influencing factor in defining the codon usage patterns and evolution of polymerase genes of EIVs.

Synonymous Codons: Choose Wisely for Expression

The genetic code, which defines the amino acid sequence of a protein, also contains information that influences the rate and efficiency of translation. Neither the mechanisms nor functions of codon-mediated regulation were well understood. The prevailing model was that the slow translation of codons decoded by rare tRNAs reduces efficiency. Recent genome-wide analyses have clarified several issues. Specific codons and codon combinations modulate ribosome speed and facilitate protein folding. However, tRNA availability is not the sole determinant of rate; rather, interactions between adjacent codons and wobble base pairing are key. One mechanism linking translation efficiency and codon use is that slower decoding is coupled to reduced mRNA stability. Changes in tRNA supply mediate biological regulationfor instance,, changes in tRNA amounts facilitate cancer metastasis.

The importance of mRNA structure in determining the pathogenicity of synonymous and non-synonymous mutations in haemophilia

INTRODUCTION

Mutational analysis is commonly used to support the diagnosis and management of haemophilia. This has allowed for the generation of large mutation databases which provide unparalleled insight into genotype-phenotype relationships. Haemophilia is associated with inversions, deletions, insertions, nonsense and missense mutations. Both synonymous and non-synonymous mutations influence the base pairing of messenger RNA (mRNA), which can alter mRNA structure, cellular half-life and ribosome processivity/elongation. However, the role of mRNA structure in determining the pathogenicity of point mutations in haemophilia has not been evaluated.

AIM

To evaluate mRNA thermodynamic stability and associated RNA prediction software as a means to distinguish between neutral and disease-associated mutations in haemophilia.

METHODS

Five mRNA structure prediction software programs were used to assess the thermodynamic stability of mRNA fragments carrying neutral vs. disease-associated and synonymous vs. non-synonymous point mutations in F8, F9 and a third X-linked gene, DMD (dystrophin).

RESULTS

In F8 and DMD, disease-associated mutations tend to occur in more structurally stable mRNA regions, represented by lower MFE (minimum free energy) levels. In comparing multiple software packages for mRNA structure prediction, a 101-151 nucleotide fragment length appears to be a feasible range for structuring future studies.

CONCLUSION

mRNA thermodynamic stability is one predictive characteristic, which when combined with other RNA and protein features, may offer significant insight when screening sequencing data for novel disease-associated mutations. Our results also suggest potential utility in evaluating the mRNA thermodynamic stability profile of a gene when determining the viability of interchanging codons for biological and therapeutic applications.

Evolutionary interpretations of mycobacteriophage biodiversity and host-range through the analysis of codon usage bias

In an genomics course sponsored by the Howard Hughes Medical Institute (HHMI), undergraduate students have isolated and sequenced the genomes of more than 1,150 mycobacteriophages, creating the largest database of sequenced bacteriophages able to infect a single host, Mycobacterium smegmatis, a soil bacterium. Genomic analysis indicates that these mycobacteriophages can be grouped into 26 clusters based on genetic similarity. These clusters span a continuum of genetic diversity, with extensive genomic mosaicism among phages in different clusters. However, little is known regarding the primary hosts of these mycobacteriophages in their natural habitats, nor of their broader host ranges. As such, it is possible that the primary host of many newly isolated mycobacteriophages is not M. smegmatis, but instead a range of closely related bacterial species. However, determining mycobacteriophage host range presents difficulties associated with mycobacterial cultivability, pathogenicity and growth. Another way to gain insight into mycobacteriophage host range and ecology is through bioinformatic analysis of their genomic sequences. To this end, we examined the correlations between the codon usage biases of 199 different mycobacteriophages and those of several fully sequenced mycobacterial species in order to gain insight into the natural host range of these mycobacteriophages. We find that UPGMA clustering tends to match, but not consistently, clustering by shared nucleotide sequence identify. In addition, analysis of GC content, tRNA usage and correlations between mycobacteriophage and mycobacterial codon usage bias suggests that the preferred host of many clustered mycobacteriophages is not M. smegmatis but other, as yet unknown, members of the mycobacteria complex or closely allied bacterial species.

Relevance of Molecular Changes in the ND4 Gene in German Shepherd Dog Tumours

The aim of the study was to identify polymorphisms and mutations in the mitochondrial ND4 gene and to analyse the associations between the occurrence of molecular changes in mtDNA and phenotypic traits in tumours in German Shepherd dogs. Fifty samples obtained from blood and tumour tissues of German Shepherd dogs with diagnosed tumours were analysed. DNA extraction, amplification, and sequencing of the mtDNA ND4 gene, and bioinformatics, statistical, and in silico protein coding SNP analyses were performed. ND4 mutations and/or polymorphisms were noted in eleven nucleotide positions in nearly half of the examined dogs. All the changes were substitution mutations. A majority of the changes identified were homoplasmic. In one dog with osteosarcoma, blood heteroplasmy was detected. In two positions of the ND4 gene, presence of non-synonymous mutations leading to amino acid changes in the ND4 protein was reported. Analyses carried out to determine the deleterious effect of mutations indicated an almost 97 and 62% probability that a single amino acid substitution (p.G239V and p.I401T, respectively) in the protein has a negative impact on its function. The results of statistical analyses indicate a significant association between the occurrence of mutations in three loci of the ND4 gene and the location of tumours. The mutations identified may be a result of cell adaptation to the changes in the environment occurring during carcinogenesis. The high frequency of mutations in the tumours may indicate genetic instability of mtDNA, which may also play a role in carcinogenesis.

The functional relevance of somatic synonymous mutations in melanoma and other cancers

Recent technological advances in sequencing have flooded the field of cancer research with knowledge about somatic mutations for many different cancer types. Most cancer genomics studies focus on mutations that alter the amino acid sequence, ignoring the potential impact of synonymous mutations. However, accumulating experimental evidence has demonstrated clear consequences for gene function, leading to a widespread recognition of the functional role of synonymous mutations and their causal connection to various diseases. Here, we review the evidence supporting the direct impact of synonymous mutations on gene function via gene splicing; mRNA stability, folding, and translation; protein folding; and miRNA-based regulation of expression. These results highlight the functional contribution of synonymous mutations to oncogenesis and the need to further investigate their detection and prioritization for experimental assessment.

A case-control study of innate immunity pathway gene polymorphisms in Puerto Ricans reveals association of toll-like receptor 2 +596 variant with asthma

BACKGROUND

For many years, African Dust Storms (ADE) has been thought to be associated with high prevalence of asthma in Puerto Rico (PR). Endotoxins (ENX) have been associated with ADE particulate matter (PM) and are known to promote pro-inflammatory responses in lung cells of susceptible individuals through the Toll-like receptor (TLR2/4) signaling pathways. Genetic variants are plausible contributors to such susceptibility. Therefore, we have evaluated a series of nine single nucleotide polymorphisms (SNPs) in TLR genes, which have been correlated positive and negatively to asthma prevalence and/or risk, in the Puerto Rican asthmatic population.

METHODS

The following SNPs were evaluated in 62 asthmatics and 61 controls through Taqman® Real Time PCR Assay: TLR4 (+896A/G, +1196C/T, -6687A/G); TLR2 (+596C/T, -16934 T/A, +399A/G, +1349C/T) and CD14 (-159C/T, +1188C/G). Genotypes were assessed for asthma association employing an odds ratio (OR) analysis.

RESULTS

Minor allele frequencies (n = 123) were determined for those variants as 0.07, 0.06, 0.35, 0.35, 0.37, 0.29, 0.04, 0.35 and 0.11, respectively. Two (+596C/T, +399A/G) TLR2 SNPs showed to be more represented in the asthmatic group by 89 % and 65 %, respectively. TLR4 SNP +896A/G analysis revealed only 1 G/G genotype (2 %) on the asthmatic group. The CD14 SNPs were similarly represented in the Puerto Rican population. Only the TLR2 +596 SNP was found to be significantly associated to asthma (OR = 3.24 for CT, 2.71 for TT) and particularly to females.

CONCLUSIONS

The identification of TLR SNPs will reveal potential candidates for gene-environment interactions in Puerto Ricans. As far as we know this is the first study to evaluate this type of TLR gene polymorphisms in Puerto Rican asthmatics, contributing to the current knowledge in the Hispanic population.

Impact of ABCB1 1236C > T-2677G > T-3435C > T polymorphisms on the anti-proliferative activity of imatinib, nilotinib, dasatinib and ponatinib

Overexpression of ABCB1 (also called P-glycoprotein) confers resistance to multiple anticancer drugs, including tyrosine kinase inhibitors (TKIs). Several ABCB1 single nucleotide polymorphisms affect the transporter activity. The most common ABCB1 variants are 1236C > T, 2677G > T, 3435C > T and have been associated with clinical response to imatinib in chronic myelogenous leukaemia (CML) in some studies. We evaluated the impact of these polymorphisms on the anti-proliferative effect and the intracellular accumulation of TKIs (imatinib, nilotinib, dasatinib and ponatinib) in transfected HEK293 and K562 cells. ABCB1 overexpression increased the resistance of cells to doxorubicin, vinblastine and TKIs. Imatinib anti-proliferative effect and accumulation were decreased to a larger extent in cells expressing the ABCB1 wild-type protein compared with the 1236T-2677T-3435T variant relatively to control cells. By contrast, ABCB1 polymorphisms influenced the activity of nilotinib, dasatinib and ponatinib to a much lesser extent. In conclusion, our data suggest that wild-type ABCB1 exports imatinib more efficiently than the 1236T-2677T-3435T variant protein, providing a molecular basis for the reported association between ABCB1 polymorphisms and the response to imatinib in CML. Our results also point to a weaker impact of ABCB1 polymorphisms on the activity of nilotinib, dasatinib and ponatinib.

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.

A Comparative Study of Ethylene Emanation upon Nitrogen Deficiency in Natural Accessions of Arabidopsis thaliana

An original approach to develop sustainable agriculture with less nitrogen fertilizer inputs is to tackle the cross-talk between nitrogen nutrition and plant growth regulators. In particular the gaseous hormone, ethylene, is a prime target for that purpose. The variation of ethylene production in natural accessions of the model species Arabidopsis thaliana was explored in response to the nitrate supply. Ethylene was measured with a laser-based photoacoustic detector. First, experimental conditions were established with Columbia-0 (Col-0) accession, which was grown in vitro on horizontal plates across a range of five nitrate concentrations (0.5, 1, 2.5, 5, or 10 mM). The concentrations of 1 and 10 mM nitrate were retained for further characterization. Along with a decrease of total dry biomass and higher biomass allocation to the roots, the ethylene production was 50% more important at 1 mM than at 10 mM nitrate. The total transcript levels of 1-AMINOCYCLOPROPANE-1-CARBOXYLIC ACID SYNTHASES (ACS) in roots and those of ACC OXIDASES (ACO) in shoots increased by 100% between the same treatments. This was mainly due to higher transcript levels of ACS6 and of ACO2 and ACO4 respectively. The assumption was that during nitrogen deficiency, the greater biomass allocation in favor of the roots was controlled by ethylene being released in the shoots after conversion of ACC originating from the roots. Second, biomass and ethylene productions were measured in 20 additional accessions. Across all accessions, the total dry biomass and ethylene production were correlated negatively at 1 mM but positively at 10 mM nitrate. Furthermore, polymorphism was surveyed in ACC and ethylene biosynthesis genes and gene products among accessions. Very few substitutions modifying the amino acids properties in conserved motifs of the enzymes were found in the accessions. Natural variation of ethylene production could be further explored to improve Nitrogen Use Efficiency (NUE), in particular by manipulating features like the biomass production and the timing of senescence upon nitrogen limitation.

A synonymous EGFR polymorphism predicting responsiveness to anti-EGFR therapy in metastatic colorectal cancer patients

Genetic factors are known to affect the efficiency of therapy with monoclonal antibodies (mAbs) targeting the epidermal growth factor receptor (EGFR) in patients with metastatic colorectal cancer (mCRC). At present, the only accepted molecular marker predictive of the response to anti-EGFR mAbs is the somatic mutation of KRAS and NRAS as a marker of resistance to anti-EGFR. However, only a fraction of KRAS wild-type patients benefit from that treatment. In this study, we show that the EGFR gene polymorphism rs1050171 defines, independently of RAS mutational status, a sub-population of 11 % of patients with a better clinical outcome after anti-EGFR treatment. Median PFS for patients with the GG genotype was 10.17 months compared to 5.37 of those with AG + AA genotypes. Taken together, our findings could be used to better define CRC populations responding to anti-EGFR therapy. Further studies in larger independent cohorts are necessary to validate the present observation that a synonymous polymorphism in EGFR gene impacts on clinical responsiveness.

Immune mediators in the brain and peripheral tissues in autism spectrum disorder

Increasing evidence points to a central role for immune dysregulation in autism spectrum disorder (ASD). Several ASD risk genes encode components of the immune system and many maternal immune system-related risk factors--including autoimmunity, infection and fetal reactive antibodies--are associated with ASD. In addition, there is evidence of ongoing immune dysregulation in individuals with ASD and in animal models of this disorder. Recently, several molecular signalling pathways--including pathways downstream of cytokines, the receptor MET, major histocompatibility complex class I molecules, microglia and complement factors--have been identified that link immune activation to ASD phenotypes. Together, these findings indicate that the immune system is a point of convergence for multiple ASD-related genetic and environmental risk factors.

Signs of Selection in Synonymous Sites of the Mitochondrial Cytochrome b Gene of Baikal Oilfish (Comephoridae) by mRNA Secondary Structure Alterations

Studies over the past decade have shown a significant role of synonymous mutations in posttranscriptional regulation of gene expression, which is particularly associated with messenger RNA (mRNA) secondary structure alterations. Most studies focused on prokaryote genomes and the nuclear genomes of eukaryotes while little is known about the regulation of mitochondrial DNA (mtDNA) gene expression. This paper reveals signs of selection in synonymous sites of the mitochondrial cytochrome b gene (Cytb) of Baikal oilfish or golomyankas (Comephoridae) directed towards altering the secondary structure of the mRNA and probably altering the character of mtDNA gene expression. Our findings are based on comparisons of intraspecific genetic variation patterns of small golomyanka (Comephorus dybowski) and two genetic groups of big golomyanka (Comephorus dybowskii). Two approaches were used: (i) analysis of the distribution of synonymous mutations between weak-AT (W) and strong-GC (S) nucleotides within species and groups in accordance with mutation directions from central to peripheral haplotypes and (ii) approaches based on the predicted mRNA secondary structure.

Quantitative Trait Loci and Candidate Genes for Neutrophil Recruitment in Sterile Inflammation Mapped in AXB-BXA Recombinant Inbred Mice

Neutrophil recruitment (NR) to sites of sterile inflammation plays a key role in tissue damage and healing potential of lesions characteristic to non-infectious inflammatory diseases. Previous studies suggested significant genetic control of neutrophil survival, function, and migration in inflammatory responses to endogenous and exogenous stimuli. We have mapped the murine genome for quantitative trait loci (QTLs) harbouring genetic determinants that regulate NR in SI using a murine model of chemically-induced peritonitis. NR was quantified in 16 AXB-BXA recombinant inbred strains and their progenitors, A/J (A) and C57BL/6J (B). A continuous distribution of NR was found among the strains, with parent B showing higher NR and parent A showing lower NR (3.0-fold difference, p=0.05). Within the progeny strains, a 5.5-fold difference in NR was observed between the lowest, BXA1, and the highest responders AXB19 (p<0.001). This data was analyzed using GeneNetwork, which linked NR to one significant QTL on chromosome 12 (Peritoneal Neutrophil Recruitment 1, PNR1) and two suggestive QTLs (PNR2, PNR3) on chromosomes 12 and 16 respectively. Sixty-four candidate genes within PNR1 were cross-referenced with currently published data, mRNA expression from two NR microarrays, and single nucleotide polymorphism analysis. The present study brings new light into the genetics of NR in response to cell injury and highlights potential candidate genes Hif1α, Fntb, and Prkch and their products for further studies on neutrophil infiltration and inflammation resolution in sterile inflammation.

Análisis computacional del efecto de polimorfismos de genes del sistema μ-Calpaína/Calpastatina sobre la calidad de la carne bovina

<p>Los genes del sistema de enzimas μ-Calpaína/Calpastatina han sido ampliamente evaluados en estudios de asociación respecto de parámetros de calidad cárnica como la terneza; previamente se han identificado varios polimorfismos asociados con la variación fenotípica en poblaciones no relacionadas de bovinos. Usando herramientas computacionales se logró postular la asociación de cuatro polimorfismos encontrados en μ-Calpaína y 11 en Calpastatina que producen una alteración de los parámetros físico-químicos, tanto del ARNm (estabilidad y polimorfismo conformacional), como de la proteína (punto isoeléctrico, potencial electroestático y superficie molecular). Es importante poder establecer el soporte biológico de polimorfismos genéticos asociados con parámetros fenotípicos que mejoren la productividad animal, lo que hace que la aproximación in silico se convierta en una herramienta útil para tal fin.</p>

Local correlations in codon preferences do not support a model of tRNA recycling

It has been proposed that patterns in the usage of synonymous codons provide evidence that individual tRNA molecules are recycled through the ribosome, translating several occurrences of the same amino acid before diffusing away. The claimed evidence is based on counting the frequency with which pairs of synonymous codons are used at nearby occurrences of the same amino acid, as compared to the frequency expected if each codon were chosen independently from a single genome-wide distribution. We show that such statistics simply measure variation in codon preferences across a genome. As a negative control on the potential contribution of pressure to exploit tRNA recycling on these signals, we examine correlations in the usage of codons that encode different amino acids. We find that these controls are statistically as strong as the claimed evidence and conclude that there is no informatic evidence that tRNA recycling is a force shaping codon usage.

Mitochondrial DNA variation analysis in cervical cancer

This study was undertaken to investigate the mitochondrial DNA (mtDNA) variation in non-malignant and malignant cervical tissue samples. We have identified 229 and 739 variations non-malignant and malignant tissues respectively distributed over 321 locations in the D-loop (50 in non-malignant and 166 in malignant; 216 variations), coding region (139 in non-malignant and 455 in malignant; 594 variations) tRNA and rRNA genes (39 in non-malignant and 119 in malignant; 158 variations). Besides, 77 novel and 34 various other disease associated variations were identified in non-malignant and malignant samples. A total of 236 tumor specific variations in 201 locations representing 30.1% in D-loop, 59.3% in coding regions and 10.6% in RNA genes were also identified. Our study shows that D loop (in 67 locations) is highly altered followed by ND5 (35 locations) region. Moreover, mtDNA alterations were significantly higher in malignant samples by two tailed Fisher's exact test (P≤0.05) with decreased mtDNA copy numbers. Bioinformatic analysis of 59 non-synonymous changes predicted several variations as damaging leading to decreased stability of the proteins. Taken together, mtDNA is highly altered in cervical cancer and functional studies are needed to be investigated to understand the consequence of these variations in cervical carcinogenesis and their potential application as biomarkers.

Mutational pressure dictates synonymous codon usage in freshwater unicellular α - cyanobacterial descendant Paulinella chromatophora and β - cyanobacterium Synechococcus elongatus PCC6301

Background

Comparative study of synonymous codon usage variations and factors influencing its diversification in α - cyanobacterial descendant Paulinella chromatophora and β - cyanobacterium Synechococcus elongatus PCC6301 has not been reported so far. In the present study, we investigated various factors associated with synonymous codon usage in the genomes of P. chromatophora and S. elongatus PCC6301 and findings were discussed.

Results

Mutational pressure was identified as the major force behind codon usage variation in both genomes. However, correspondence analysis revealed that intensity of mutational pressure was higher in S. elongatus than in P. chromatophora. Living habitats were also found to determine synonymous codon usage variations across the genomes of P. chromatophora and S. elongatus.

Conclusions

Whole genome sequencing of α-cyanobacteria in the cyanobium clade would certainly facilitate the understanding of synonymous codon usage patterns and factors contributing its diversification in presumed ancestors of photosynthetic endosymbionts of P. chromatophora.

Molecular signature of high yield (growth) influenza a virus reassortants prepared as candidate vaccine seeds

Background

Human influenza virus isolates generally grow poorly in embryonated chicken eggs. Hence, gene reassortment of influenza A wild type (wt) viruses is performed with a highly egg adapted donor virus, A/Puerto Rico/8/1934 (PR8), to provide the high yield reassortant (HYR) viral ‘seeds’ for vaccine production. HYR must contain the hemagglutinin (HA) and neuraminidase (NA) genes of wt virus and one to six ‘internal’ genes from PR8. Most studies of influenza wt and HYRs have focused on the HA gene. The main objective of this study is the identification of the molecular signature in all eight gene segments of influenza A HYR candidate vaccine seeds associated with high growth in ovo.

Methodology

The genomes of 14 wt parental viruses, 23 HYRs (5 H1N1; 2, 1976 H1N1-SOIV; 2, 2009 H1N1pdm; 2 H2N2 and 12 H3N2) and PR8 were sequenced using the high-throughput sequencing pipeline with big dye terminator chemistry.

Results

Silent and coding mutations were found in all internal genes derived from PR8 with the exception of the M gene. The M gene derived from PR8 was invariant in all 23 HYRs underlining the critical role of PR8 M in high yield phenotype. None of the wt virus derived internal genes had any silent change(s) except the PB1 gene in X-157. The highest number of recurrent silent and coding mutations was found in NS. With respect to the surface antigens, the majority of HYRs had coding mutations in HA; only 2 HYRs had coding mutations in NA.

Significance

In the era of application of reverse genetics to alter influenza A virus genomes, the mutations identified in the HYR gene segments associated with high growth in ovo may be of great practical benefit to modify PR8 and/or wt virus gene sequences for improved growth of vaccine ‘seed’ viruses.

Rapid evolution of a few members of nasuta-albomicans complex of Drosophila: study on two candidate genes, Sod1 and Rpd3

Drosophila nasuta nasuta (2n = 8) and D. n. albomicans (2n = 6) are morphologically identical, cross fertile and karyotypically dissimilar pair of chromosomal races belonging to nasuta subgroup of immigrans group of Drosophila. Interracial hybridization between these two races yielded karyotypically stabilized newly evolved Cytoraces with new combinations of chromosomes and DNA content, and are called nasuta-albomicans complex of Drosophila. Along with many other features, striking plasticity in the lifespan has been observed in the karyotypically stabilized members of nasuta-albomicans complex of Drosophila. These findings provide a strong background to understand any changes at the molecular levels. In view of this, we cloned and characterized Sod1 and Rpd3 in the members of nasuta-albomicans complex of Drosophila. The evolution of Sod1 and Rpd3 in D. n. nasuta and D. n. albomicans is contrasting with the other species of Drosophila, at the level of synonymous mutations, intron variation, InDels and secondary structure changes in protein. In the members of NAC of Drosophila there were synonymous changes, variations in intron sequences of Sod1, whereas, in Rpd3, synonymous, nonsynonymous, intron variation, and secondary structure changes in protein were observed. The contrasting differences in the levels of Rpd3 (and Sir2) proteins were also noticed among short-lived and long-lived Cytoraces. The Cytoraces have exhibited not only specific changes in Sod1 and Rpd3, but also show pronounced changes in the levels of synthesis of these proteins, which indicates rapid evolution of these Cytoraces in laboratory. Further these Cytoraces have become a model system to understand the process of anagenesis.

Chapter 15: disease gene prioritization

Disease-causing aberrations in the normal function of a gene define that gene as a disease gene. Proving a causal link between a gene and a disease experimentally is expensive and time-consuming. Comprehensive prioritization of candidate genes prior to experimental testing drastically reduces the associated costs. Computational gene prioritization is based on various pieces of correlative evidence that associate each gene with the given disease and suggest possible causal links. A fair amount of this evidence comes from high-throughput experimentation. Thus, well-developed methods are necessary to reliably deal with the quantity of information at hand. Existing gene prioritization techniques already significantly improve the outcomes of targeted experimental studies. Faster and more reliable techniques that account for novel data types are necessary for the development of new diagnostics, treatments, and cure for many diseases.

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Good codons, bad transcript: large reductions in gene expression and fitness arising from synonymous mutations in a key enzyme

Biased codon usage in protein-coding genes is pervasive, whereby amino acids are largely encoded by a specific subset of possible codons. Within individual genes, codon bias is stronger at evolutionarily conserved residues, favoring codons recognized by abundant tRNAs. Although this observation suggests an overall pattern of selection for translation speed and/or accuracy, other work indicates that transcript structure or binding motifs drive codon usage. However, our understanding of codon bias evolution is constrained by limited experimental data on the fitness effects of altering codons in functional genes. To bridge this gap, we generated synonymous variants of a key enzyme-coding gene in Methylobacterium extorquens. We found that mutant gene expression, enzyme production, enzyme activity, and fitness were all significantly lower than wild-type. Surprisingly, encoding the gene using only rare codons decreased fitness by 40%, whereas an allele coded entirely by frequent codons decreased fitness by more than 90%. Increasing gene expression restored mutant fitness to varying degrees, demonstrating that the fitness disadvantage of synonymous mutants arose from a lack of beneficial protein rather than costs of protein production. Protein production was negatively correlated with the frequency of motifs with high affinity for the anti-Shine-Dalgarno sequence, suggesting ribosome pausing as the dominant cause of low mutant fitness. Together, our data support the idea that, although a particular set of codons are favored on average across a genome, in an individual gene selection can either act for or against codons depending on their local context.

Engineered mutations change the structure and stability of a virus-like particle

The single-coat protein (CP) of bacteriophage Qβ self-assembles into T = 3 icosahedral virus-like particles (VLPs), of interest for a wide range of applications. These VLPs are very stable, but identification of the specific molecular determinants of this stability is lacking. To investigate these determinants along with manipulations that confer more capabilities to our VLP material, we manipulated the CP primary structure to test the importance of various putative stabilizing interactions. Optimization of a procedure to incorporate fused CP subunits allowed for good control over the average number of covalent dimers in each VLP. We confirmed that the disulfide linkages are the most important stabilizing elements for the capsid and that acidic conditions significantly enhance the resistance of VLPs to thermal degradation. Interdimer interactions were found to be less important for VLP assembly than intradimer interactions. Finally, a single point mutation in the CP resulted in a population of smaller VLPs in three distinct structural forms.

A simple method for analyzing exome sequencing data shows distinct levels of nonsynonymous variation for human immune and nervous system genes

To measure the strength of natural selection that acts upon single nucleotide variants (SNVs) in a set of human genes, we calculate the ratio between nonsynonymous SNVs (nsSNVs) per nonsynonymous site and synonymous SNVs (sSNVs) per synonymous site. We transform this ratio with a respective factor f that corrects for the bias of synonymous sites towards transitions in the genetic code and different mutation rates for transitions and transversions. This method approximates the relative density of nsSNVs (rdnsv) in comparison with the neutral expectation as inferred from the density of sSNVs. Using SNVs from a diploid genome and 200 exomes, we apply our method to immune system genes (ISGs), nervous system genes (NSGs), randomly sampled genes (RSGs), and gene ontology annotated genes. The estimate of rdnsv in an individual exome is around 20% for NSGs and 30-40% for ISGs and RSGs. This smaller rdnsv of NSGs indicates overall stronger purifying selection. To quantify the relative shift of nsSNVs towards rare variants, we next fit a linear regression model to the estimates of rdnsv over different SNV allele frequency bins. The obtained regression models show a negative slope for NSGs, ISGs and RSGs, supporting an influence of purifying selection on the frequency spectrum of segregating nsSNVs. The y-intercept of the model predicts rdnsv for an allele frequency close to 0. This parameter can be interpreted as the proportion of nonsynonymous sites where mutations are tolerated to segregate with an allele frequency notably greater than 0 in the population, given the performed normalization of the observed nsSNV to sSNV ratio. A smaller y-intercept is displayed by NSGs, indicating more nonsynonymous sites under strong negative selection. This predicts more monogenically inherited or de-novo mutation diseases that affect the nervous system.