A testis‐specific gene within a widely expressed gene: Contrasting evolutionary patterns of two differentially expressed mammalian proteins encoded by a single gene,CAMK4

Understanding the patterns of genetic variations within fertility‐related genes and the evolutionary forces that shape such variations is crucial in predicting the fitness landscapes of subsequent generations. This study reports distinct evolutionary features of two differentially expressed mammalian proteins [CaMKIV (Ca²⁺/calmodulin‐dependent protein kinase IV) and CaS (calspermin)] that are encoded by a single gene, CAMK4. The multifunctional CaMKIV, which is expressed in multiple tissues including testis and ovary, is evolving at a relatively low rate (0.46–0.64 × 10^?9 nucleotide substitutions/site/year), whereas the testis‐specific CaS gene, which is predominantly expressed in post‐meiotic cells, evolves at least three to four times faster (1.48–1.98 × 10^?9 substitutions/site/year). Concomitantly, maximum‐likelihood‐based selection analyses revealed that the ubiquitously expressed CaMKIV is constrained by intense purifying selection and, therefore, remained functionally highly conserved throughout the mammalian evolution, whereas the testis‐specific CaS gene is under strong positive selection. The substitution rates of different mammalian lineages within both genes are positively correlated with GC content, indicating the possible influence of GC‐biased gene conversion on the estimated substitution rates. The observation of such unusually high GC content of the CaS gene (≈74%), particularly in the lineage that comprises the bovine species, suggests the possible role of GC‐biased gene conversion in the evolution of CaS that mimics positive selection.     

158 A molecular dynamics simulation-based prediction of deleterious angiogenin mutations causing amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the selective death of motor neurons leading to paralysis and death between 3–5?years of diagnosis. Through whole genome association studies, several single nucleotide polymorphisms (SNPs) encoding missense mutations in angiogenin (ANG) protein have been identified as one of the primary factors causing ALS. Structural studies of ANG show that catalytic triad comprising His13, Lys40, and His114 residues imparts ribonucleolytic activity while nuclear localization signal residues ³¹RRR³³ are responsible for nuclear translocation activity. Loss of either ribonucleolytic activity or nuclear translocation activity or both of these functions due to mutations cause ALS. However, the mechanisms of loss-of-functions of ANG mutants are not completely understood. Here, we present a cohesive and comprehensive picture of functional loss mechanisms of all known ALS-associated ANG mutants by extensive molecular dynamics (MD) simulations (Padhi, Kumar, Vasaikar, Jayaram, & Gomes, 2012 Padhi, A. K., Kumar, H., Vasaikar, S. V., Jayaram, B. and Gomes, J. 2012. Mechanisms of loss of functions of human angiogenin variants implicated in amyotrophic lateral sclerosis. PLoS One, 7(2): e32479[PubMed] , [Google Scholar]; AK, 2013 Padhi, A.K., Jayaram, B., & Gomes, J. (accepted for publication). Prediction of functional loss of human angiogenin mutants associated with ALS by molecular dynamics simulations. Scientific Reports (NPG).  [Google Scholar]). Our studies show that conformational switching of catalytic residue His114 is responsible for the loss of ribonucleolytic activity while reduction in solvent-accessible surface area (SASA) of ³¹RRR³³ as a result of local folding is responsible for the loss of nuclear translocation activity (Padhi et al., 2012 Padhi, A. K., Kumar, H., Vasaikar, S. V., Jayaram, B. and Gomes, J. 2012. Mechanisms of loss of functions of human angiogenin variants implicated in amyotrophic lateral sclerosis. PLoS One, 7(2): e32479[PubMed] , [Google Scholar]; AK, 2013 Padhi, A.K., Jayaram, B., & Gomes, J. (accepted for publication). Prediction of functional loss of human angiogenin mutants associated with ALS by molecular dynamics simulations. Scientific Reports (NPG).  [Google Scholar]). Our prediction of loss-of-functions of 17 ANG mutants correlated positively with the reported experimental results. We have subsequently developed a fast molecular dynamics method based on certain global attributes / dynamic markers that can be used to determine whether a mutation is deleterious or benign. To make our method accessible to researchers and clinicians, we created a web server-based tool, ANGDelMut, freely available at http://bioschool.iitd.ernet.in/research.htm, where a user can submit new mutations to ascertain whether they cause ALS. We hope that our method will benefit the community at large and will pave the way for the development of a successful therapy for patients suffering from ALS.     

Molecular Evolutionary and Epidemiological Dynamics of Genotypes 1G and 2B of Rubella Virus

Rubella Virus (RV), which causes measles-like rashes in children, puts millions of infants at risk of congenital defects across the globe. Employing phylogenetic approaches to the whole genome sequence data and E1 glycoprotein sequence data, the present study reports the substitution rates and dates of emergence of all thirteen previously described rubella genotypes, and gains important insights into the epidemiological dynamics of two geographically widely distributed genotypes 1G and 2B. The overall nucleotide substitution rate of this non-vector-borne RV is in the order of 10⁻³ substitutions/site/year, which is considerably higher than the substitution rates previously reported for the vector-borne alphaviruses within the same family. Currently circulating strains of RV share a common ancestor that existed within the last 150 years, with 95% Highest Posterior Density values ranging from 1868 to 1926 AD. Viral strains within the respective genotypes began diverging between the year 1930 s and 1980 s. Both genotype 1G and 2B have shown a decline in effective number of infections since 1990 s, a period during which mass immunization programs against RV were adapted across the globe. Although both genotypes showed some extent of spatial genetic structuring, the analyses also depicted an inter-continental viral dispersal. Such a viral dispersal pattern could be related to the migration of infected individuals across the regions coupled with a low coverage of MMR vaccination.     

Computational and Functional Characterization of Angiogenin Mutations,and Correlation with Amyotrophic Lateral Sclerosis

The Angiogenin (ANG) gene is frequently mutated in patients suffering from the neurodegenerative disease - amyotrophic lateral sclerosis (ALS). Most of the ALS-causing mutations in Angiogenin affect either its ribonucleolytic or nuclear translocation activity. Here we report the functional characterization of two previously uncharacterized missense mutations in Angiogenin - D22G and L35P. We predict the nature of loss-of-function(s) in these mutants through our previously established Molecular Dynamics (MD) simulation extended to 100 ns, and show that the predictions are entirely validated through biochemical studies with wild-type and mutated proteins. Based on our studies, we provide a biological explanation for the loss-of-function of D22G-Angiogenin leading to ALS, and suggest that the L35P-Angiogenin mutation would probably cause ALS symptoms in individuals harboring this mutation. Our study thus highlights the strength of MD simulation-based predictions, and suggests that this method can be used for correlating mutations in Angiogenin or other effector proteins with ALS symptoms.     

A Link between Arabinose Utilization and Oxalotrophy in Bradyrhizobium japonicum

Rhizobia have a versatile catabolism that allows them to compete successfully with other microorganisms for nutrients in the soil and in the rhizosphere of their respective host plants. In this study, Bradyrhizobium japonicum USDA 110 was found to be able to utilize oxalate as the sole carbon source. A proteome analysis of cells grown in minimal medium containing arabinose suggested that oxalate oxidation extends the arabinose degradation branch via glycolaldehyde. A mutant of the key pathway genes oxc (for oxalyl-coenzyme A decarboxylase) and frc (for formyl-coenzyme A transferase) was constructed and shown to be (i) impaired in growth on arabinose and (ii) unable to grow on oxalate. Oxalate was detected in roots and, at elevated levels, in root nodules of four different B. japonicum host plants. Mixed-inoculation experiments with wild-type and oxc-frc mutant cells revealed that oxalotrophy might be a beneficial trait of B. japonicum at some stage during legume root nodule colonization.     

Pattern of nucleotide substitution and divergence of prophenoloxidase in decapods

Despite the unprecedented development in identification and characterization of prophenoloxidase (proPO) in commercially important decapods, little is known about the evolutionary relationship, rate of amino acid replacement and differential selection pressures operating on proPO of different species of decapods. Here we report the evolutionary relationship among these nine decapod species based on proPO gene and types of selective pressures operating on proPO codon sites. Our analyses revealed that all the nine decapod species shared a common ancestor. The mean percentage sequence divergence at proPO gene was 34.4+/-0.6%. Pairwise estimates of nonsynonymous to synonymous ratio (omega) for Homarus americanus-H. gammarus is greater than one, therefore indicating adaptive evolution (functional diversification) of proPO in these two species. In contrast, strong purifying selection (omega<1) was observed in all other species pairs. However, phylogenetically closely related decapods revealed relatively higher omega value (omega=0.15+/-0.3) than the distantly related species pairs (omega=0.0075+/-0.005). These discrepancies could be due to higher fixation probability of beneficial mutation in closely related species. Maximum likelihood-based codon substitution analyses revealed a strong purifying selection operating on most of the codon sites, therefore suggesting proPO is functionally constrained (purifying selection). Codon substitution analyses have also revealed the evidence of strong purifying selection in haemocyanin subunits of decapods.