A non-synonymous SNP in homolog of BADH2 gene is associated with fresh pod fragrance in dolichos bean (Lablab purpureus var. lignosus (Prain) Kumari)

Fresh pods are harvestable and marketable economic product in dolichos bean. Fresh pod fragrance is one of the 'farmers' and 'consumers' preferred traits in dolichos bean varieties. The pods with high fragrance fetch a premium price in the market. In breeding programmes, pod fragrance is routinely assessed by organoleptic (sensory) means, which is highly relative and subjective. Identification of linked DNA markers not only offer an objective means but also enable selection of fragment genotypes at seedling stage itself. Betaine aldehyde dehydrogenase (BADH) is known to be the key gene responsible for fragrance in other legumes such as vegetable soybean and mung bean. Taking cues from highly conserved domains in proteins coded by BADH genes, we isolated dolichos bean homolog (LpBADH2) of soybean GmBADH2 gene using reported degenerate primers designed to conserved domains. Analysis of the translated amino acid sequence of LpBADH2 showed high degree of similarity (97.30%) with those of soybean homolog (GmBADH2). Conserved amino acid sequence of aldehyde dehydrogenase-super family were also identified in LpBADH2. Multiple sequence alignment of nucleotide sequences of LpBADH2 with those of related legumes using "ClustalW" revealed the presence of a single non-synonymous single nucleotide polymorphic (SNPs) and three synonymous SNP sites in LpBADH2. The substitution of the amino acid tyrosine in (fragrant genotypes) with phenyl alanine (non-fragrant genotypes) in protein coded by LpBADH2 appeared to be the cause for switch over from fragrance to non-fragrance in dolichos bean. These results are discussed in relation to strategies to breed dolichos bean cultivars with desired level of pod fragrance.

A single nucleotide polymorphism of the thyrotropin releasing hormone degrading ectoenzyme (TRHDE) gene is associated with post-partum anestrus in Murrah buffalo

Thyrotropin releasing hormone degrading enzyme (TRHDE) gene is implicated in Thyrotropin releasing hormone (TRH) mediated prolactin secretion. It has been shown that the prolactin secretion alters the Gonadotropin-releasinghormone(GnRH) mediated estrous cycle. Therefore, TRHDE may also regulate postpartum anestrus. Earlier studies reported the role of non-synonymous single nucleotide polymorphism (SNPs) in various pathophysiological conditions by altering the structure and function of the proteins. Hence, in the present study, we identified SNPs in the putative promoter, first exon, middle exon and 3'-UTR containing the last exon of TRHDE gene and determined their association with postpartum anestrus (PPA) in Murrah buffaloes. We found one non synonymous SNP (G > C at 118095875 bp on chromosome 4) in the first exon of TRHDE and performed its association analysis in a population sample of 50 extreme PPA (residual PPAI: 123.06 ± 12.98 days) and 50 normal (residual PPAI: -80.46 ± 3.19 days) buffaloes. The residual PPAI value was the observed PPAI adjusted for the effect of 38 non-genetic factors. The analysis showed a significant (P < 0.004167) association of this SNP with PPA in buffaloes. Molecular dynamics simulations (MDS) also supported that the C allele altering Glutamine to Histidine at the amino acid 148 of TRHDE could enhance the stability and rigidity of TRHDE protein, which may lower its activity, increase TRH and prolactin, and reduce GnRH in PPA buffaloes. The MDS analysis further strengthens the association of the SNP (G > C) in the TRHDE gene with PPA condition in Murrah buffaloes. However, further investigation is needed to prove the MDS observations.

Genome-wide analysis of 10664 SARS-CoV-2 genomes to identify virus strains in 73 countries based on single nucleotide polymorphism

Since the onslaught of SARS-CoV-2, the research community has been searching for a vaccine to fight against this virus. However, during this period, the virus has mutated to adapt to the different environmental conditions in the world and made the task of vaccine design more challenging. In this situation, the identification of virus strains is very much timely and important task. We have performed genome-wide analysis of 10664 SARS-CoV-2 genomes of 73 countries to identify and prepare a Single Nucleotide Polymorphism (SNP) dataset of SARS-CoV-2. Thereafter, with the use of this SNP data, the advantage of hierarchical clustering is taken care of in such a way so that Average Linkage and Complete Linkage with Jaccard and Hamming distance functions are applied separately in order to identify the virus strains as clusters present in the SNP data. In this regard, the consensus of both the clustering results are also considered while Silhouette index is used as a cluster validity index to measure the goodness of the clusters as well to determine the number of clusters or virus strains. As a result, we have identified five major clusters or virus strains present worldwide. Apart from quantitative measures, these clusters are also visualized using Visual Assessment of Tendency (VAT) plot. The evolution of these clusters are also shown. Furthermore, top 10 signature SNPs are identified in each cluster and the non-synonymous signature SNPs are visualised in the respective protein structures. Also, the sequence and structural homology-based prediction along with the protein structural stability of these non-synonymous signature SNPs are reported in order to judge the characteristics of the identified clusters. As a consequence, T85I, Q57H and R203M in NSP2, ORF3a and Nucleocapsid respectively are found to be responsible for Cluster 1 as they are damaging and unstable non-synonymous signature SNPs. Similarly, F506L and S507C in Exon are responsible for both Clusters 3 and 4 while Clusters 2 and 5 do not exhibit such behaviour due to the absence of any non-synonymous signature SNPs. In addition to all these, the code, SNP dataset, 10664 labelled SARS-CoV-2 strains and additional results as supplementary are provided through our website for further use.

Identification and characterization of trait-specific SNPs using ddRAD sequencing in water buffalo

Single Nucleotide Polymorphism (SNP) is one of the important molecular markers widely used in animal breeding program for improvement of any desirable genetic traits. Considering this, the present study was carried out to identify, annotate and analyze the SNPs related to four important traits of buffalo viz. milk volume, age at first calving, post-partum cyclicity and feed conversion efficiency. We identified 246,495, 168,202, 74,136 and 194,747 genome-wide SNPs related to mentioned traits, respectively using ddRAD sequencing technique based on 85 samples of Murrah Buffaloes. Distribution of these SNPs were highest (61.69%) and lowest (1.78%) in intron and exon regions, respectively. Under coding regions, the SNPs for the four traits were further classified as synonymous (4697) and non-synonymous (3827). Moreover, Gene Ontology (GO) terms of identified genes assigned to various traits. These characterized SNPs will enhance the knowledge of cellular mechanism for enhancing productivity of water buffalo through molecular breeding.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.

Identification of candidate protective variants for common diseases and evaluation of their protective potential

BACKGROUND

Human polymorphisms with derived alleles that are protective against disease may provide powerful translational opportunities. Here we report a method to identify such candidate polymorphisms and apply it to common non-synonymous SNPs (nsSNPs) associated with common diseases. Our study also sought to establish which of the identified protective nsSNPs show evidence of positive selection, taking this as indirect evidence that the protective variant has a beneficial effect on phenotype. Further, we performed an analysis to quantify the predicted effect of each protective variant on protein function/structure.

RESULTS

An initial analysis of eight SNPs previously identified as associated with age-related macular degeneration (AMD), revealed that two of them have a derived allele that is protective against developing the disease. One is in the complement component 2 gene (C2; E318D) and the other is in the complement factor B gene (CFB; R32Q). Then, combining genomewide ancestral allele information with known common disease-associated nsSNPs from the GWAS catalog, we found 32 additional SNPs which have a derived allele that is disease protective. Out of the total 34 identified candidate protective variants (CPVs), we found that 30 show stronger evidence of positive selection than the protective variant in lipoprotein lipase (LPL; S447X), which has already been translated into gene therapy. Furthermore, 11 of these CPVs have a higher probability of affecting protein structure than the lipoprotein lipase protective variant (LPL; S447X).

CONCLUSIONS

We identify 34 CPVs from the human genome. Diseases they confer protection against include, but are not limited to, type 2 diabetes, inflammatory bowel disease, age-related macular degeneration, multiple sclerosis and rheumatoid arthritis. We propose that those 30 CPVs with evidence of stronger positive selection than the LPL protective variant, may be considered as priority candidates for therapeutic approaches. The next step towards translation will require testing the hypotheses generated by our analyses, specifically whether the CPV arose from a gain-of-function or a loss-of-function mutation.

Analysis of consequences of non-synonymous SNP in feed conversion ratio associated TGF-β receptor type 3 gene in chicken

The recent advances in high throughput sequencing technology accelerate possible ways for the study of genome wide variation in several organisms and associated consequences. In the present study, mutations in TGFBR3 showing significant association with FCR trait in chicken during exome sequencing were further analyzed. Out of four SNPs, one nsSNP p.Val451Leu was found in the coding region of TGFBR3. In silico tools such as SnpSift and PANTHER predicted it as deleterious (0.04) and to be tolerated, respectively, while I-Mutant revealed that protein stability decreased. The TGFBR3 I-TASSER model has a C-score of 0.85, which was validated using PROCHECK. Based on MD simulation, mutant protein structure deviated from native with RMSD 0.08 Å due to change in the H-bonding distances of mutant residue. The docking of TGFBR3 with interacting TGFBR2 inferred that mutant required more global energy. Therefore, the present study will provide useful information about functional SNPs that have an impact on FCR traits.

•

Investigated functional nsSNP p.Val451Leu (rs312979494) in feed conversion ratio (FCR) associated TGFBR3 of chicken

•

Computational tools (SIFT and I-Mutant 2.0) predicted that this nsSNP was deleterious.

•

Mutant structure of TGFBR3 showed high energies and RMS deviations compared to native using MD simulation.

•

Molecular docking of TGFBR3 with interacting protein TGFBR2 showed an increase in global energy of mutant compared to native.

•

We have predicted that functional SNP has an impact on TGFBR3 of chicken and thus can be treated as candidate SNP in screening.