Analysis of phylogeny and codon usage bias and relationship of GC content, amino acid composition with expression of the structural nif genes

Using codon usage analysis to speculate potential animal hosts of hepatitis E virus: An exploratory study

There has been an increase in the reported number of animals worldwide that carry the hepatitis E virus (HEV). This study aimed to explore potential animal hosts for HEV through codon usage analysis. Full-length HEV sequences of six genotypes as well as codon usage of potential animal hosts were collected. Moreover, nucleotide composition and codon usage bias were compared across HEV genotypes and animal hosts. Based on the analysis for human HEV-1 and humans, the results were basically consistent with epidemiology evidence. Among 17 potential animal hosts, all HEV genotypes exhibited a preference for guanine/cytosine in the third position of synonymous codons. Furthermore, non-human primates and humans have large high-frequency codons identical to HEV in addition to a high correlation of codon fraction with HEV. Some animals in close contact with humans showed high preference for HEV, including cattle, dogs, and rats with HEV-A, cats, dogs, and swine with HEV-C1. Codon usage bias has limited efficiency in determining the hosts for HEV, but it may provide indicative clues for potential animal hosts when combined with experimental and epidemiological evidence.

Virtual screening and docking analysis of novel ligands for selective enhancement of tea (Camellia sinensis) flavonoids

•

Tea-specific flavonoid biosynthetic pathway (FBP) was retrieved from KEGG.

•

Putative ligands were predicted to enhance enzymes-substrate binding affinity.

•

FBP genes showed moderately higher expression & relatively strong codon adaptation.

•

Most of the genes were AT-rich and biased to A/U-ending synonymous codons.

•

Mutational selection was determining the selective constraints on codon bias.

Flavour of tea is mainly contributed by a group of polyphenols – flavonoids. However, the content of flavonoid fluctuates seasonally and is found to be higher in the first flush of tea, when compared to the second flush. This disparity in the flavonoid content, and hence taste, incurs heavy economic losses to the tea plantation industry each harvest season. For our present study, four key product-specific enzymes (PAL, FNS, FLS and ANS) of the tea-specific flavonoid pathway were selected to perform molecular docking studies with specific virtually screened allosteric modulators. Results of docking analyses showed Naringenin, 2-Morpholin-4-ium-4-ylethanesulfonate, 6-C-Glucosylquercetin, 2-Oxoglutaric acid, 3,5,7,3′,4′-pentahydroxyflavone to be capable of improving the spontaneity of the enzyme-substrate reactions in terms of docking score, RMSD values, and non-covalent interactions (H-bond,hydrophobic interaction, Π-stacking, salt bridge, etc.). Further, the evolutionary relationship of tea flavonoid pathway enzymes was constructed and compared with related taxa. The codon usage-based of tea flavonoid biosynthetic genes indicated the non-biasness of their nucleotide composition. Overall this study will provide a direction towards putative ligand-dependent enhancement of flavonoid content, irrespective of seasonal variation.

Microbiome Analysis via OTU and ASV-Based Pipelines-A Comparative Interpretation of Ecological Data in WWTP Systems

Linking community composition and ecosystem function via the cultivation-independent analysis of marker genes, e.g., the 16S rRNA gene, is a staple of microbial ecology and dependent disciplines. The certainty of results, independent of the bioinformatic handling, is imperative for any advances made within the field. In this work, thermophilic anaerobic co-digestion experimental data, together with primary and waste-activated sludge prokaryotic community data, were analyzed with two pipelines that apply different principles when dealing with technical, sequencing, and PCR biases. One pipeline (VSEARCH) employs clustering methods, generating individual operational taxonomic units (OTUs), while the other (DADA2) is based on sequencing error correction algorithms and generates exact amplicon sequence variants (ASVs). The outcomes of both pipelines were compared within the framework of ecological-driven data analysis. Both pipelines provided comparable results that would generally allow for the same interpretations. Yet, the two approaches also delivered community compositions that differed between 6.75% and 10.81% between pipelines. Inconsistencies were also observed linked to biologically driven variability in the samples, which affected the two pipelines differently. These pipeline-dependent differences in taxonomic assignment could lead to different conclusions and interfere with any downstream analysis made for such mis- or not-identified species, e.g., network analysis or predictions of their respective ecosystem service.

The complete chloroplast genome of Stauntonia chinensis and compared analysis revealed adaptive evolution of subfamily Lardizabaloideae species in China

BACKGROUND

Stauntonia chinensis DC. belongs to subfamily Lardizabaloideae, which is widely grown throughout southern China. It has been used as a traditional herbal medicinal plant, which could synthesize a number of triterpenoid saponins with anticancer and anti-inflammatory activities. However, the wild resources of this species and its relatives were threatened by over-exploitation before the genetic diversity and evolutionary analysis were uncovered. Thus, the complete chloroplast genome sequences of Stauntonia chinensis and comparative analysis of chloroplast genomes of Lardizabaloideae species are necessary and crucial to understand the plastome evolution of this subfamily.

RESULTS

A series of analyses including genome structure, GC content, repeat structure, SSR component, nucleotide diversity and codon usage were performed by comparing chloroplast genomes of Stauntonia chinensis and its relatives. Although the chloroplast genomes of eight Lardizabaloideae plants were evolutionary conserved, the comparative analysis also showed several variation hotspots, which were considered as highly variable regions. Additionally, pairwise Ka/Ks analysis showed that most of the chloroplast genes of Lardizabaloideae species underwent purifying selection, whereas 25 chloroplast protein coding genes were identified with positive selection in this subfamily species by using branch-site model. Bayesian and ML phylogeny on CCG (complete chloroplast genome) and CDs (coding DNA sequences) produced a well-resolved phylogeny of Lardizabaloideae plastid lineages.

CONCLUSIONS

This study enhanced the understanding of the evolution of Lardizabaloideae and its relatives. All the obtained genetic resources will facilitate future studies in DNA barcode, species discrimination, the intraspecific and interspecific variability and the phylogenetic relationships of subfamily Lardizabaloideae.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

The Complete Mitochondrial Genome of endemic giant tarantula, Lyrognathus crotalus (Araneae: Theraphosidae) and comparative analysis

The complete mitochondrial genome of Lyrognathus crotalus is sequenced, annotated and compared with other spider mitogenomes. It is 13,865 bp long and featured by 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and a control region (CR). Most of the PCGs used ATN start codon except cox3, and nad4 with TTG. Comparative studies indicated the use of TTG, TTA, TTT, GTG, CTG, CTA as start codons by few PCGs. Most of the tRNAs were truncated and do not fold into the typical cloverleaf structure. Further, the motif (CATATA) was detected in CR of nine species including L. crotalus. The gene arrangement of L. crotalus compared with ancestral arthropod showed the transposition of five tRNAs and one tandem duplication random loss (TDRL) event. Five plesiomophic gene blocks (A-E) were identified, of which, four (A, B, D, E) retained in all taxa except family Salticidae. However, block C was retained in Mygalomorphae and two families of Araneomorphae (Hypochilidae and Pholcidae). Out of 146 derived gene boundaries in all taxa, 15 synapomorphic gene boundaries were identified. TreeREx analysis also revealed the transposition of trnI, which makes three derived boundaries and congruent with the result of the gene boundary mapping. Maximum likelihood and Bayesian inference showed similar topologies and congruent with morphology, and previously reported multi-gene phylogeny. However, the Gene-Order based phylogeny showed sister relationship of L. crotalus with two Araneomorphae family members (Hypochilidae and Pholcidae) and other Mygalomorphae species.

Loss of phosphatase activity in PTEN (phosphatase and tensin homolog deleted on chromosome ten) results in endometrial carcinoma in humans: An in-silico study

The tumour suppressor gene, PTEN (Phosphatase and Tensin homolog deleted on chromosome Ten), can act as both protein phosphatase and lipid phosphatase, is known to play a vital role in Pi3k signalling pathway. In humans, it is located at 10q23. Loss of its phosphatase and catalytic activity is associated with various types of cancers. This study focuses on evolution, understanding the somatic missense mutation in a particular residue of PTEN and understanding the molecular mechanism that leads to endometrial carcinoma through molecular docking. Mutational analysis of H123 position indicates that the missense mutation at first position of the codon CAC by G or T, result in aspartic acid or tyrosine instead of histidine and can have negative effect on the function of PTEN. Alongside, structural analysis showed mutated PTEN has lower stability than the normal. Additionally, SNPs dataset for endometrial carcinoma suggests H123 as strongly mutated residue. The mutation in phosphatase domain of PTEN along with its effect and interaction with substrate TLA1352 were systematically studied through molecular docking. Molecular interaction study reveals that the optimal substrate binding site in PTEN is unable to interact with the substrate in the mutated condition. This observation drew attention on the impact of mutation on disease biology and enabled us to conduct follow-up studies to retrieve novel molecular targets, such as mutated protein domain and modified Asp and Tyr sites, to design effective therapies to either prevent endometrial carcinoma or impede its progression.

Base composition is the primary factor responsible for the variation of amino acid usage in zebra finch (Taeniopygia guttata)

In the present study, we carried out an examination of the amino acid usage in the zebra finch (Taeniopygia guttata) proteome. We found that tRNA abundance, base composition, hydrophobicity and aromaticity, protein second structure, cysteine residue (Cys) content and protein molecular weight had significant impact on the amino acid usage of the zebra finch. The above factors explained the total variability of 22.85%, 25.37%, 10.91%, 5.06%, 4.21%, and 3.14%, respectively. Altogether, approximately 70% of the total variability in zebra finch could be explained by such factors. Comparison of the amino acid usage between zebra finch, chicken (Gallus gallus) and human (Homo sapiens) suggested that the average frequency of various amino acid usage is generally consistent among them. Correspondence analysis indicated that base composition was the primary factor affecting the amino acid usage in zebra finch. This trend was different from chicken, but similar to human. Other factors affecting the amino acid usage in zebra finch, such as isochore structure, protein second structure, Cys frequency and protein molecular weight also showed the similar trends with human. We do not know whether the similar amino acid usage trend between human and zebra finch is related to the distinctive neural and behavioral traits, but it is worth studying in depth.

Comprehensive Analysis of Codon Usage on Rabies Virus and Other Lyssaviruses

Rabies virus (RABV) and other lyssaviruses can cause rabies and rabies-like diseases, which are a persistent public health threat to humans and other mammals. Lyssaviruses exhibit distinct characteristics in terms of geographical distribution and host specificity, indicative of a long-standing diversification to adapt to the environment. However, the evolutionary diversity of lyssaviruses, in terms of codon usage, is still unclear. We found that RABV has the lowest codon usage bias among lyssaviruses strains, evidenced by its high mean effective number of codons (ENC) (53.84 ± 0.35). Moreover, natural selection is the driving force in shaping the codon usage pattern of these strains. In summary, our study sheds light on the codon usage patterns of lyssaviruses, which can aid in the development of control strategies and experimental research.

Part by Part: Synthetic Biology Parts Used in Solventogenic Clostridia

The solventogenic Clostridia are of interest to the chemical industry because of their natural ability to produce chemicals such as butanol, acetone and ethanol from diverse feedstocks. Their use as whole cell factories presents multiple metabolic engineering targets that could lead to improved sustainability and profitability of Clostridium industrial processes. However, engineering efforts have been held back by the scarcity of genetic and synthetic biology tools. Over the past decade, genetic tools to enable transformation and chromosomal modifications have been developed, but the lack of a broad palette of synthetic biology parts remains one of the last obstacles to the rapid engineered improvement of these species for bioproduction. We have systematically reviewed existing parts that have been used in the modification of solventogenic Clostridia, revealing a narrow range of empirically chosen and nonengineered parts that are in current use. The analysis uncovers elements, such as promoters, transcriptional terminators and ribosome binding sites where increased fundamental knowledge is needed for their reliable use in different applications. Together, the review provides the most comprehensive list of parts used and also presents areas where an improved toolbox is needed for full exploitation of these industrially important bacteria.

Codon usage variation of Zika virus: The potential roles of NS2B and NS4A in its global pandemic

A comprehensive demonstration of Zika virus (ZIKV) molecular evolution is essential for understanding its adaptation and expansion in its recent pandemics. Despite several studies on mutations and codon usage in ZIKVs, the variations in codon usage patterns across individual genes and their biological implication remains unclear. Here, we performed a gene-by-gene comparison of the codon usage variation in ZIKVs of the African and Asian lineages. We found that besides the evidence of positive selection (Ka/Ks >1) in the Asian lineage of the ZIKV genome, codon usage patterns were gene-specific and codon usage variation of ZIKV genes, was possibly constrained by their individual functional features, such as transmembrane domains, or antigenicity. In particular, the NS2B and NS4A genes showed distinct codon usage patterns, clearly separating them from the clusters of other genes in the correspondence analysis (CA). In the Asian lineage, the NS2B and NS4A genes showed the highest codon usage bias (ENC values: 51.01 ± 0.72 and 48.89 ± 0.99 respectively), and were subjected to the highest translation selection (ENC_obs/ENC_exp ratio: 0.847 ± 0.0297 and 0.828 ± 0.0233 respectively) in comparison to the African lineages of ZIKV. The CpG frequency of the NS2B showed a gradual ascending trend in the Asian ZIKV lineages, while in NS4A it was constrained along with the expansion of the Asian lineage. Furthermore, between the African and Asian lineages, differentiated and specific over-represented codons were more prominent in the NS2B and NS4A. Together, our study implies that ZIKVs are in the process of evolutionary fine tuning their codon as seen in the recent pandemics, and NS2B and NS4A could have played a potential role in the molecular evolution of the Asian lineage and their establishment.

New insights into the factors affecting synonymous codon usage in human infecting Plasmodium species

Codon usage bias is due to the non-random usage of synonymous codons for coding amino acids. The synonymous sites are under weak selection, and codon usage bias is maintained by the equilibrium in mutational bias, genetic drift and selection pressure. The differential codon usage choices are also relevant to human infecting Plasmodium species. Recently, P. knowlesi switches its natural host, long-tailed macaques, and starts infecting humans. This review focuses on the comparative analysis of codon usage choices among human infecting P. falciparum and P. vivax along with P. knowlesi species taking their coding sequence data. The variation in GC content, amino acid frequencies, effective number of codons and other factors plays a crucial role in determining synonymous codon choices. Within species codon choices are more similar for P. vivax and P. knowlesi in comparison with P. falciparum species. This study suggests that synonymous codon choice modulates the gene expression level, mRNA stability, ribosome speed, protein folding, translation efficiency and its accuracy in Plasmodium species, and provides a valuable information regarding the codon usage pattern to facilitate gene cloning as well as expression and transfection studies for malaria causing species.

Genome-wide sequential, evolutionary, organizational and expression analyses of phenylpropanoid biosynthesis associated MYB domain transcription factors in Arabidopsis

The MYB gene family represents one of the largest groups of transcription factors in plants. Recent evidences have also demonstrated key role of MYB transcription factors in regulating the expression of major genes involved in the biosynthesis of phenylpropanoid compounds which confer biotic and abiotic stress tolerance in plant species. However, no comprehensive genome-wide analysis of the phenylpropanoid pathway-associated MYB transcription factors has been reported thus far. In this study, 11 Arabidopsis MYB proteins, such as MYB3, MYB4, MYB7, MYB11, MYB12, MYB32, MYB75, MYB90, MYB111, MYB113, and MYB114 were initially identified considering their reported regulatory function in phenylpropanoid biosynthesis pathway. Subsequent genome-wide analysis have identified the corresponding homologues from Glycine max, Vigna radiata, Oryza sativa, and Zea mays, while homologous of Arabidopsis MYB75, MYB90, MYB113, and MYB114 were not detected in rice and maize genomes. The identified MYB proteins were classified into three groups (I-III) based on phylogeny. Sequence and domain analysis revealed presence of two conserved DNA binding MYB domains in the selected MYB proteins. Promoter analysis indicated presence of cis-regulatory elements related to light signaling, development, and stress response. Expression analysis of selected Arabidopsis MYB genes revealed their function in plant development and abiotic stress response, consistent with gene ontology annotations. Together, these results provide a useful framework for further experimental studies for the functional characterization of the target MYB genes in the context of regulation of phenylpropanoid biosynthesis and plant stress response.