Human SNPs resulting in premature stop codons and protein truncation

Truncated variants of thyroid hormone receptor beta display disease-inflicting malfunctioning at cellular level

Thyroid hormone receptor β (THRβ) is a member of the nuclear receptor superfamily of ligand-modulated transcription factors. Upon ligand binding, THRβ sequentially recruits the components of transcriptional machinery to modulate target gene expression. In addition to regulating diverse physiological processes, THRβ plays a crucial role in hypothalamus-pituitary-thyroid axis feedback regulation. Anomalies in THRβ gene/protein structure are associated with onset of diverse disease states. In this study, we investigated disease-inflicting truncated variants of THRβ using in-silico analysis and cell-based assays. We examined the THRβ truncated variants on multiple test parameters, including subcellular localization, ligand-receptor interactions, transcriptional functions, interaction with heterodimeric partner RXR, and receptor-chromatin interactions. Moreover, molecular dynamic simulation approaches predicted that shortened THRβ-LBD due to point mutations contributes proportionally to the loss of structural integrity and receptor stability. Deviant subcellular localization and compromised transcriptional function were apparent with these truncated variants. Present study shows that 'mitotic bookmarking' property of some THRβ variants is also affected. The study highlights that structural and conformational attributes of THRβ are necessary for normal receptor functioning, and any deviations may contribute to the underlying cause of the inflicted diseases. We anticipate that insights derived herein may contribute to improved mechanistic understanding to assess disease predisposition.

Increased hsa-miR-100-5p Expression Improves Hepatocellular Carcinoma Prognosis in the Asian Population with PLK1 Variant rs27770A>G

Hepatocellular carcinoma (HCC) has the highest incidence and mortality in the Asian population, and race is an independent risk factor affecting survival time in liver cancer. Micro RNAs (miRNAs) are remarkably dysregulated in HCC and closely associated with HCC prognosis. Recent studies show that genetic variability between ethnic groups may result in differences in the specificity of HCC miRNA biomarkers. Here, we reveal a high expression level of hsa-miR-100-5p, an HCC prognosis-related miRNA, which improves HCC prognosis in the Asian Population with Polo-like kinase 1 (PLK1) variant rs27770A>G. In this study, we discovered that hsa-miR-100-5p was downregulated in various HCC cell lines. While mimics transient transfection and mouse liver cancer model confirmed the interaction between hsa-miR-100-5p and PLK1, a stratified analysis based on the Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) data suggest both low hsa-miR-100-5p expression level and high PLK1 expression level associated with poor HCC prognosis, especially in the Asian population. According to the 1000 Genomes Project database, the SNP rs27770 located in 3'UTR of PLK1 had a significantly higher G allele frequency in the East Asian population. Bioinformatics analysis suggested that rs27770 A>G affects PLK1 mRNA secondary structure and alters the hsa-miR-100-5p/PLK1 interaction by forming an additional seedless binding site. This racial variation caused PLK1 to be more vulnerable to hsa-miR-100-5p inhibition, resulting in hsa-miR-100-5p being more favorable for HCC prognosis in the Asian population.

Mitotic genome-bookmarking by nuclear hormone receptors: A novel dimension in epigenetic reprogramming and disease assessment

Arrival of multi-colored fluorescent proteins and advances in live cell imaging has immensely contributed to our understanding of intracellular trafficking of nuclear receptors and their roles in gene regulatory functions. These regulatory events need to be faithfully propagated from progenitor to progeny cells. This is corroborated by multiple converging mechanisms that include histone modifications and lately, the phenomenon of 'mitotic genome-bookmarking' by specific transcription factors. This phenomenon refers to the retention and feed-forward transmission of progenitor's architectural blueprint of active transcription status which is silenced and preserved during mitosis. Upon mitotic exit, this phenomenon ensures accurate reactivation of transcriptome, proteome, cellular traits and phenotypes in the progeny cells. In addition to diverse modes of genome-bookmarking by nuclear receptors, a correlation between disease-associated receptor polymorphism and disruption of this phenomenon is apparent. However, breakthrough technologies shall reveal finer details of this phenomenon to help achieve normalcy in receptor-specific diseases.

Loss-of-function variants in KCTD19 cause non-obstructive azoospermia in humans

Azoospermia is a significant cause of male infertility, with non-obstructive azoospermia (NOA) being the most severe type of spermatogenic failure. NOA is mostly caused by congenital factors, but our understanding of its genetic causes is very limited. Here, we identified a frameshift variant (c.201_202insAC, p.Tyr68Thrfs∗17) and two nonsense variants (c.1897C>T, p.Gln633∗; c.2005C>T, p.Gln669∗) in KCTD19 (potassium channel tetramerization domain containing 19) from two unrelated infertile Chinese men and a consanguineous Pakistani family with three infertile brothers. Testicular histological analyses revealed meiotic metaphase I (MMI) arrest in the affected individuals. Mice modeling KCTD19 variants recapitulated the same MMI arrest phenotype due to severe disrupted individualization of MMI chromosomes. Further analysis showed a complete loss of KCTD19 protein in both Kctd19 mutant mouse testes and affected individual testes. Collectively, our findings demonstrate the pathogenicity of the identified KCTD19 variants and highlight an essential role of KCTD19 in MMI chromosome individualization.

Predicting heterosis via genetic distance and the number of SNPs in selected segments of chromosomes in maize

A reliable method is needed for predicting heterosis to help maize (Zea mays L.) breeders develop new hybrids more efficiently. The objectives of this study were to 1) investigate if the numbers of selected PEUS SNPs (the SNP in the Promoters (1 kb upstream of the start codon), Exons, Untranslated region (UTR), and Stop codons) could be used for predicting MPH or BPH of GY; 2) if the number of PEUS SNPs is a better predictor of MPH and/or BPH of GY than genetic distance (GD). A line × tester experiment was conducted with 19 elite maize inbreds from three heterotic groups, which were crossed with five testers. The multi-location trial data on GY were recorded. Whole-genome resequencing of the 24 inbreds was carried out. After filtration, a total of 58,986,791 SNPs were called with high confidence. Selected SNPs in the promoters, exons, untranslated region (UTRs), and stop codons (PEUS SNPs) were counted, and the GD was calculated. The correlation between heterozygous PEUS SNPs/GD and mean MPH, BPH of GY revealed that 1) both the number of heterozygous PEUS SNP and the GD were highly correlated to both MPH_GY and BPH_GY at p<0.01 with correlation coefficients for the number of heterozygous PEUS SNP being higher than that for GD; 2) the mean number of heterozygous PEUS SNPs was also highly correlated with mean BPH_GY or mean MPH_GY (p<0.05) in the 95 crosses grouped by either male or female parents, implying that inbreds can be selected before making the actual crosses in the field. We concluded that the number of heterozygous PEUS SNPs would be a better predictor of MPH_GY and BPH_GY than GD. Hence, maize breeders could use heterozygous PEUS SNPs to select inbreds with high heterosis potential before actually making the crosses, thus improving the breeding efficiency.

Beta-defensins as marker for male fertility: a comprehensive review†

Bovine male fertility in animals has a direct impact on the productivity of dairy herds. The epididymal sperm maturations involve extensive sperm surface modifications to gain the fertilizing ability, especially by absorptions of the plethora of biomolecules, including glycoprotein beta-defensins (BDs), enzymes, organic ions, protein, and phospholipids. Defensins are broad-range nonspecific antimicrobial peptides that exhibit strong relations with innate and adaptive immunity, but their roles in male fertility are relatively recently identified. In the course of evolution, BD genes give rise to different clusters with specific functions, especially reproductive functions, by undergoing duplications and nonsynonymous mutations. BD polymorphisms have been reported with milk compositions, disease resistance, and antimicrobial activities. However, in recent decades, the link of BD polymorphisms with fertility has emerged as an appealing improvement of reproductive performance such as sperm motility, membrane integrity, cervical mucus penetration, evading of uterus immunosurveillance, oviduct cell attachment, and egg recognition. The reproductive-specific glycosylated BD class-A BDs (CA-BDs) have shown age- and sex-specific expressions in male reproductive organs, signifying their physiological pleiotropism, especially in the sperm maturation and sperm transport in the female reproductive tract. By considering adult male reproductive organ-specific BD expressions, importance in sperm functionalities, and bioinformatic analysis, we have selected two bovine BBD126 and BBD129 genes as novel potential biomarkers of bovine male fertility. Despite the importance of BDs, however, genomic characterization of most BD genes across most livestock and nonmodel organisms remains predictive/incomplete. The current review discusses our understanding of BD pleiotropic functions, polymorphism, and genomic structural attributes concerning the fertilizability of the male gamete in dairy animals.

CRISPR nuclease off-target activity and mitigation strategies

The discovery of CRISPR has allowed site-specific genomic modification to become a reality and this technology is now being applied in a number of human clinical trials. While this technology has demonstrated impressive efficacy in the clinic to date, there remains the potential for unintended on- and off-target effects of CRISPR nuclease activity. A variety of in silico-based prediction tools and empirically derived experimental methods have been developed to identify the most common unintended effect-small insertions and deletions at genomic sites with homology to the guide RNA. However, large-scale aberrations have recently been reported such as translocations, inversions, deletions, and even chromothripsis. These are more difficult to detect using current workflows indicating a major unmet need in the field. In this review we summarize potential sequencing-based solutions that may be able to detect these large-scale effects even at low frequencies of occurrence. In addition, many of the current clinical trials using CRISPR involve ex vivo isolation of a patient's own stem cells, modification, and re-transplantation. However, there is growing interest in direct, in vivo delivery of genome editing tools. While this strategy has the potential to address disease in cell types that are not amenable to ex vivo manipulation, in vivo editing has only one desired outcome-on-target editing in the cell type of interest. CRISPR activity in unintended cell types (both on- and off-target) is therefore a major safety as well as ethical concern in tissues that could enable germline transmission. In this review, we have summarized the strengths and weaknesses of current editing and delivery tools and potential improvements to off-target and off-tissue CRISPR activity detection. We have also outlined potential mitigation strategies that will ensure that the safety of CRISPR keeps pace with efficacy, a necessary requirement if this technology is to realize its full translational potential.

Biliary atresia combined Wilson disease identified by whole exome sequencing in Vietnamese patient with severe liver failure

RATIONALE

Hepatobiliary diseases such as biliary atresia (BA), Wilson disease, and progressive familial intrahepatic cholestasis are common causes of morbidity and mortality in young children. Affected patients progress rapidly to end-stage cirrhosis and require liver transplantation or die. Mutations in many genes have been identified to play an important role in the pathogenesis of hepatobiliary diseases.

PATIENT CONCERNS AND DIAGNOSIS

In this study, we identified mutations in an 8-year-old girl who had severe liver failure. The patient was first diagnosed with BA at 2.5 months of age and has undergone Kasai surgery to connect the umbilical cord and jejunum. After that, the patient suddenly had unusual developments with symptoms of jaundice, acute liver failure with hemolysis. She was tested and diagnosed with Wilson disease.

INTERVENTIONS AND OUTCOMES

She was treated according to the regimen for a patient with Wilson disease but had abnormal progress leading to severe liver failure. Genetic analysis was performed by whole exome sequencing and Sanger sequencing methods. The genetic analysis revealed that the patient had a homozygous mutation (p.Gly17Glyfs77∗) in the KRT18 gene, a double heterozygous mutation (p.Ser105∗ and p.Pro992Leu) in the ATP7B gene, and a homozygous variant (p.Val444Ala) in the ABCB11 gene. In silico prediction of mutations indicated that these mutations are the cause of the severe liver failure in the patient.

LESSON

This is a rare clinical case of a BA patient combined with Wilson disease. Our results suggested that whole exome sequencing is an effective diagnostic tool and emphasizes the importance of early diagnosis and appropriate management to save lives and prevent serious complications in the patient.

Hypothalamic kisspeptin and kisspeptin receptors: Species variation in reproduction and reproductive behaviours

Kisspeptin, encoded by the KISS1 gene, was first discovered as a potential metastasis suppressor gene. The prepro-kisspeptin precursor is cleaved into shorter mature bioactive peptides of varying sizes that bind to the G protein-coupled receptor GPR54 (=KISS1R). Over the last two decades, multiple types of Kiss and KissR genes have been discovered in mammalian and non-mammalian vertebrate species, but they are remarkably absent in birds. Kiss neuronal populations are distributed mainly in the hypothalamus. The KissRs are widely distributed in the brain, including the hypothalamic and non-hypothalamic regions, such as the hippocampus, amygdala, and habenula. The role of KISS1-KISS1R in humans and Kiss1-Kiss1R in rodents is associated with puberty, gonadal maturation, and the reproductive axis. However, recent gene deletion studies in zebrafish and medaka have provided controversial results, suggesting that the reproductive role of kiss is dispensable. This review highlights the evolutionary history, localisation, and significance of Kiss-KissR in reproduction and reproductive behaviours in mammalian and non-mammalian vertebrates.

The link of ERCC2 rs13181 and ERCC4 rs2276466 polymorphisms with breast cancer in the Bangladeshi population

BACKGROUND

Breast cancer (BC) is the most common disease in women and the leading cause of death from cancer globally. Epidemiological studies examined that nucleotide excision repair genes ERCC2 (rs13181) and ERCC4 (rs2276466) SNPs might increase cancer risk. Based on the previous investigation, this study was conducted to explore the correlation between these polymorphisms and BC susceptibility in Bangladeshi women.

METHODS AND RESULTS

Between January 2019 and January 2020, 140 blood samples were collected from female patients histologically diagnosed with BC, and 111 female controls were recruited from non-cancer subjects. Genotyping was performed applying the PCR-RFLP method, and all statistical analyzes were conducted using SPSS, version 25.0. Comparison of characteristics and clinicopathological features between ERCC2 rs13181 and ERCC4 rs2276466 carriers and non-carriers showed no significant link with BC. Analysis of ERCC2 rs13181 with the risk of BC showed that the GG genotype and G allele carriers showed a fourfold and 1.78-fold higher risk (OR 4.00, P = 0.001 and OR 1.78, P = 0.002) that are statistically significant. In addition, the patients with combined TG+GG genotype revealed a 2.09-fold increased chance (OR 2.09, P = 0.020) BC development. Analysis of recessive model (GG vs. TT+TG) also depicted 2.74-times significantly higher risk (OR 2.74, P = 0.002). On the other hand, ERCC4 rs2276466 polymorphism did not show any significant association with BC (P > 0.05).

CONCLUSIONS

Our findings show that ERCC2 rs13181 is linked to an elevated risk of BC. Our study also shows that ERCC4 rs2276466 polymorphism has no substantial risk of BC in the Bangladeshi population.

Whole-Genome Profiles of Malay Colorectal Cancer Patients with Intact MMR Proteins

Background: This study aimed to identify new genes associated with CRC in patients with normal mismatch repair (MMR) protein expression. Method: Whole-genome sequencing (WGS) was performed in seven early-age-onset Malay CRC patients. Potential germline genetic variants, including single-nucleotide variations and insertions and deletions (indels), were prioritized using functional and predictive algorithms. Results: An average of 3.2 million single-nucleotide variations (SNVs) and over 800 indels were identified. Three potential candidate variants in three genes—IFNE, PTCH2 and SEMA3D—which were predicted to affect protein function, were identified in three Malay CRC patients. In addition, 19 candidate genes—ANKDD1B, CENPM, CLDN5, MAGEB16, MAP3K14, MOB3C, MS4A12, MUC19, OR2L8, OR51Q1, OR51AR1, PDE4DIP, PKD1L3, PRIM2, PRM3, SEC22B, TPTE, USP29 and ZNF117—harbouring nonsense variants were prioritised. These genes are suggested to play a role in cancer predisposition and to be associated with cancer risk. Pathway enrichment analysis indicated significant enrichment in the olfactory signalling pathway. Conclusion: This study provides a new spectrum of insights into the potential genes, variants and pathways associated with CRC in Malay patients.

Nonsense suppression therapies in human genetic diseases

About 11% of all human disease-associated gene lesions are nonsense mutations, resulting in the introduction of an in-frame premature translation-termination codon (PTC) into the protein-coding gene sequence. When translated, PTC-containing mRNAs originate truncated and often dysfunctional proteins that might be non-functional or have gain-of-function or dominant-negative effects. Therapeutic strategies aimed at suppressing PTCs to restore deficient protein function-the so-called nonsense suppression (or PTC readthrough) therapies-have the potential to provide a therapeutic benefit for many patients and in a broad range of genetic disorders, including cancer. These therapeutic approaches comprise the use of translational readthrough-inducing compounds that make the translational machinery recode an in-frame PTC into a sense codon. However, most of the mRNAs carrying a PTC can be rapidly degraded by the surveillance mechanism of nonsense-mediated decay (NMD), thus decreasing the levels of PTC-containing mRNAs in the cell and their availability for PTC readthrough. Accordingly, the use of NMD inhibitors, or readthrough-compound potentiators, may enhance the efficiency of PTC suppression. Here, we review the mechanisms of PTC readthrough and their regulation, as well as the recent advances in the development of novel approaches for PTC suppression, and their role in personalized medicine.

CRISPR-Cas9 cytidine and adenosine base editing of splice-sites mediates highly-efficient disruption of proteins in primary and immortalized cells

CRISPR-Cas9 cytidine and adenosine base editors (CBEs and ABEs) can disrupt genes without introducing double-stranded breaks by inactivating splice sites (BE-splice) or by introducing premature stop (pmSTOP) codons. However, no in-depth comparison of these methods or a modular tool for designing BE-splice sgRNAs exists. To address these needs, we develop SpliceR ( http://z.umn.edu/spliceR ) to design and rank BE-splice sgRNAs for any Ensembl annotated genome, and compared disruption approaches in T cells using a screen against the TCR-CD3 MHC Class I immune synapse. Among the targeted genes, we find that targeting splice-donors is the most reliable disruption method, followed by targeting splice-acceptors, and introducing pmSTOPs. Further, the CBE BE4 is more effective for disruption than the ABE ABE7.10, however this disparity is eliminated by employing ABE8e. Collectively, we demonstrate a robust method for gene disruption, accompanied by a modular design tool that is of use to basic and translational researchers alike.

Novel SCN5A variants identified in a group of Iranian Brugada syndrome patients

Brugada syndrome (BrS) is a rare hereditary arrhythmia syndrome that increases an individual's risk for sudden cardiac death (SCD) due to ventricular fibrillation. This disorder is regarded as a notable cause of death in individuals aged less than 40 years, responsible for up to 40% of sudden deaths in cases without structural heart disease, and is reported to be an endemic in Asian countries. Mutations in SCN5A are found in approximately 30% of patients with Brugada syndrome. This study aimed to investigate mutations in the SCN5A gene in a group of Iranian Brugada syndrome patients. Nine probands (n = 9, male, mean age = 39) diagnosed with Brugada syndrome were enrolled in this study. Exon 2 to 29 were amplified by PCR and subjected to direct sequencing. Eight in silico prediction tools were used to anticipate the effects of non-synonymous variants. Seven known polymorphisms and 2 previously reported disease-causing mutations, including H558R and G1406R, were found in the studied cases. Twenty novel variants were identified: 15 missense, 2 frameshift, 2 synonymous, and one nonsense variants. In silico tools predicted 11 non-synonymous variants to have damaging effects, whereas frameshift and nonsense variants were considered inherently pathogenic. The novel variants identified in this study, alongside previously reported mutations, are highly likely to be the cause of the Brugada syndrome phenotype observed in the patient group. Further analysis is required to understand the physiological effects caused by these variants.

Genetic variants of BIRC3 and NRG1 in the NLRP3 inflammasome pathway are associated with non-small cell lung cancer survival

The nod-like receptor protein 3 (NLRP3) is one of the most characterized inflammasomes, and its genetic variation and functional dysregulation are involved in pathogenesis of several cancers. To systematically evaluate the role of NLRP3 in predicting outcomes of patients with non-small cell lung cancer (NSCLC), we performed a two-phase analysis for associations between genetic variants in NLRP3 inflammasome pathway genes and NSCLC survival by using a published genome-wide association study (GWAS) dataset from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial. We used multivariate Cox proportional hazards regression analysis with Bayesian false discovery probability (≤0.80) for multiple testing correction to evaluate associations between 20,730 single-nucleotide polymorphisms (SNPs) in 176 genes and overall survival of 1,185 NSCLC patients from the PLCO trial. We further validated the identified significant SNPs in another GWAS dataset with survival data from 984 NSCLC patients of the Harvard Lung Cancer Susceptibility (HLCS) study. The results showed that two independent SNPs in two different genes (i.e., BIRC3 rs11225211 and NRG1 rs4733124) were significantly associated with the NSCLC overall survival, with a combined hazards ratio (HR) of 0.83 [95% confidence interval (CI) = 0.74-0.93 and P = 0.0009] and 1.18 (95% CI = 1.06-1.31) and P = 0.002], respectively. However, further expression quantitative trait loci (eQTL) analysis showed no evidence for correlations between the two SNPs and mRNA expression levels of corresponding genes. These results indicated that genetic variants in the NLRP3 imflammasome pathway gene-sets might be predictors of NSCLC survival, but the molecular mechanisms underlying the observed associations warrant further investigations.

Non-synonymous polymorphisms in candidate gene associated with growth traits in Channel catfish (Ictalurus punctatus, Rafinesque, 1818)

Resequencing of Myostatin, Growth Hormone, Follistatin-A-like, Insulin-like Growth Factor I (IGF-I) and Myogenin (MYOG) genes was completed to discover novel genetic variations and assess non synonymous (ns) polymorphisms (SNPs) effect on growth related traits of channel catfish. Wild and farmed animals were used as a discovering population. Resequencing lead to the identification of 59 new variants in the five analyzed genes; 66% found in introns and 34% in coding regions. From coding regions, 14 variants were synonyms and six were ns variations. A mutation rate of one in 129 bp was estimated. Four ns variations were selected for validation and association analysis. In IGF-I two ns polymorphisms, at IGF-I19 the G wild type allele was fixed in population and for IGF-I63 the C allele had a frequency of 0.972 and for mutate allele G of 0.027. In MYOG two ns SNPs were assessed. MYOG131 presented a frequency of alleles T and A, of 0.754 and 0.246, respectively and MYOG233, with a frequency of G and C of 0.775 and 0.225, respectively. Only MYOG131 (g.529T>A) was significantly associated (P < 0.04) to some growth traits. Results suggest MYOG131 g.529T>A as candidate locus for genetic enhancement of growth traits in channel catfish.

The association of two common polymorphisms in miRNAs with diabetes mellitus: A meta-analysis

BACKGROUND

MicroRNAs (miRNAs) are small noncoding single-stranded RNAs with a length of ∼21 nucleotides. Single nucleotide polymorphisms (SNPs) may affect the function of miRNAs, resulting in a variety of disorders in vivo. Recently, diabetes mellitus (DM) has become a global healthcare problem, and several studies have reported that 2 common polymorphisms (miRNA 146a rs2910164 and miRNA 27a rs895819) are related to susceptibility to diabetes. Given that no consensus had been reached regarding the association of the 2 polymorphisms with diabetes, we conducted this meta-analysis.

METHODS

Four databases (PubMed, EMBASE, Cochrane, and Web of Science) were searched up to January 9, 2019. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to evaluate the association strength. Subgroup and sensitivity analyses were also performed.

RESULTS

Six studies involving 2585 cases and 2435 controls for miR146a rs2910164 and 5 studies involving 2922 cases and 2781 controls for miR27a rs895819 were ultimately analyzed in our meta-analysis. Based on pooled results, no statistical significance in association between rs2910164 and diabetes in Caucasians, Asians, or type 2 diabetes was observed in any genetic models. Nevertheless, we found a significant correlation between miRNA27a rs895819 and diabetes in the homozygote model (CC vs TT: OR = 0.58, 95%CI [0.35,0.98]) and recessive model (CC vs CT + TT: OR = 0.59, 95%CI [0.36,0.97]). By performing subgroup analysis, we also observed that C allele conveyed a significant protective effect against diabetes development in Caucasians (C vs T: OR = 0.67, 95%CI [0.52,0.85]).

CONCLUSION

In conclusion, this meta-analysis indicated that miRNA27a rs895819 might play a protective role in diabetes, and miRNA146a rs2910164 likely had no association with diabetes.

Genome-wide discovery of DNA polymorphisms among chickpea cultivars with contrasting seed size/weight and their functional relevance

Seed size/weight is a major agronomic trait which determine crop productivity in legumes. To understand the genetic basis of seed size determination, we sought to identify DNA polymorphisms between two small (Himchana 1 and Pusa 362) and two large-seeded (JGK 3 and PG 0515) chickpea cultivars via whole genome resequencing. We identified a total of 75535 single nucleotide polymorphisms (SNPs), 6486 insertions and deletions (InDels), 1938 multi-nucleotide polymorphisms (MNPs) and 5025 complex variants between the two small and two large-seeded chickpea cultivars. Our analysis revealed 814, 244 and 72 seed-specific genes harboring DNA polymorphisms in promoter or non-synonymous and large-effect DNA polymorphisms, respectively. Gene ontology analysis revealed enrichment of cell growth and division related terms in these genes. Among them, at least 22 genes associated with quantitative trait loci, and those involved in cell growth and division and encoding transcription factors harbored promoter and/or large-effect/non-synonymous DNA polymorphisms. These also showed higher expression at late-embryogenesis and/or mid-maturation stages of seed development in the large-seeded cultivar, suggesting their role in seed size/weight determination in chickpea. Altogether, this study provided a valuable resource for large-scale genotyping applications and a few putative candidate genes that might play crucial role in governing seed size/weight in chickpea.

Single-Nucleotide Variations of the Human Nuclear Hormone Receptor Genes in 60,000 Individuals

Nuclear hormone receptors (NRs) mediate biologic actions of lipophilic molecules to gene transcription and are phylogenetically and functionally categorized into seven subfamilies and three groups, respectively. Single-nucleotide variations (SNVs) or polymorphisms are genetic changes influencing individual response to environmental factors and susceptibility to various disorders, and are part of the genetic diversification and basis for evolution. We sorted out SNVs of the human NR genes from 60,706 individuals, calculated three parameters (percentage of all variants, percentage of loss-of-function variants, and ratio of nonsynonymous/synonymous variants in their full protein-coding or major domain-coding sequences), and compared them with several valuables. Comparison of these parameters between NRs and control groups identified that NRs form a highly conserved gene family. The three parameters for the full coding sequence are positively correlated with each other, whereas four NR genes are distinct from the others with much higher tolerance to protein sequence-changing variants. DNA-binding domain and N-terminal domain are respectively those bearing the least and the most variation. NR subfamilies based on their phylogenetic proximity or functionality as well as diversity of tissue distribution and numbers of partner molecules are all not correlated with the variation parameters, whereas their gene age demonstrates an association. Our results suggest that the natural selection driving the NR family evolution still operates in humans. Gene age and probably the potential to adapt to various new ligands, but not current functional diversity, are major determinants for SNVs of the human NR genes.

NMD Classifier: A reliable and systematic classification tool for nonsense-mediated decay events

Nonsense-mediated decay (NMD) degrades mRNAs that include premature termination codons to avoid the translation and accumulation of truncated proteins. This mechanism has been found to participate in gene regulation and a wide spectrum of biological processes. However, the evolutionary and regulatory origins of NMD-targeted transcripts (NMDTs) have been less studied, partly because of the complexity in analyzing NMD events. Here we report NMD Classifier, a tool for systematic classification of NMD events for either annotated or de novo assembled transcripts. This tool is based on the assumption of minimal evolution/regulation–an event that leads to the least change is the most likely to occur. Our simulation results indicate that NMD Classifier can correctly identify an average of 99.3% of the NMD-causing transcript structural changes, particularly exon inclusions/exclusions and exon boundary alterations. Researchers can apply NMD Classifier to evolutionary and regulatory studies by comparing NMD events of different biological conditions or in different organisms.

SNPs in microRNA target sites and their potential role in human disease

In the post-genomic era, the goal of personalized medicine is to determine the correlation between genotype and phenotype. Developing high-throughput genotyping technologies such as genome-wide association studies (GWAS) and the 1000 Genomes Project (http://www.internationalgenome.org/about/#1000G_PROJECT) has dramatically enhanced our ability to map where changes in the genome occur on a population level by identifying millions of single nucleotide polymorphisms (SNPs). Polymorphisms, particularly those within the coding regions of proteins and at splice junctions, have received the most attention, but it is also now clear that polymorphisms in the non-coding regions are important. In these non-coding regions, the enhancer and promoter regions have received the most attention, whereas the 3'-UTR regions have until recently been overlooked. In this review, we examine how SNPs affect microRNA-binding sites in these regions, and how mRNA stability changes can lead to disease pathogenesis.

Revisiting the morbid genome of Mendelian disorders

BACKGROUND

The pathogenicity of many Mendelian variants has been challenged by large-scale sequencing efforts. However, many rare and benign "disease mutations" are difficult to analyze due to their rarity. The Saudi Arabian variome is enriched for homozygosity due to inbreeding, a key advantage that can be exploited for the critical examination of previously published variants.

RESULTS

We collated all "disease-related mutations" listed in the Human Gene Mutation Database (HGMD) and ClinVar, including "variants of uncertain significance" (VOUS). We find that the use of public databases including 1000 Genomes, ExAC, and Kaviar can reclassify many of these variants as likely benign. Our Saudi Human Genome Program (SHGP) can reclassify many variants that are rare in public databases. Furthermore, SGPD allows us to observe many previously reported variants in the homozygous state and our extensive phenotyping of participants makes it possible to demonstrate the lack of phenotype for these variants, thus challenging their pathogenicity despite their rarity. We also find that 18 VOUS BRCA1 and BRCA2 variants that are listed in BRCA Exchange are present at least once in the homozygous state in patients who lack features of Fanconi anemia. Reassuringly, we could reciprocally demonstrate that none of those labeled as "pathogenic" were observed in the homozygous statue in individuals who lack Fanconi phenotype in our database.

CONCLUSION

Our study shows the importance of revisiting disease-related databases using public resources as well as of population-specific resources to improve the specificity of the morbid genome of Mendelian diseases in humans.

The distribution pattern of genetic variation in the transcript isoforms of the alternatively spliced protein-coding genes in the human genome

By enabling the transcription of multiple isoforms from the same gene locus, alternative-splicing mechanisms greatly expand the diversity of the human transcriptome and proteome. Currently, the alternatively spliced transcripts from each protein-coding gene locus in the human genome can be classified as either principal or non-principal isoforms, providing that they differ with respect to cross-species conservation or biological features. By mapping the variants from the 1000 Genomes Project onto the coding region of each isoform, an interesting pattern of the genetic variation distributions of the coding regions for these two types of transcript isoforms was revealed on a whole-genome scale: compared with the principal isoform-specific coding regions, the non-principal isoform-specific coding regions are significantly enriched in amino acid-changing variants, particularly those that have a strong impact on protein function and have higher derived allele frequencies, suggesting that non-principal isoform-specific substitutions are less likely to be related to phenotype changes or disease. The results herein can help us better understand the potential consequences of alternatively spliced products from a population perspective.

Functional analysis of deep intronic SNP rs13438494 in intron 24 of PCLO gene

The single nucleotide polymorphism (SNP) rs13438494 in intron 24 of PCLO was significantly associated with bipolar disorder in a meta-analysis of genome-wide association studies. In this study, we performed functional minigene analysis and bioinformatics prediction of splicing regulatory sequences to characterize the deep intronic SNP rs13438494. We constructed minigenes with A and C alleles containing exon 24, intron 24, and exon 25 of PCLO to assess the genetic effect of rs13438494 on splicing. We found that the C allele of rs13438494 reduces the splicing efficiency of the PCLO minigene. In addition, prediction analysis of enhancer/silencer motifs using the Human Splice Finder web tool indicated that rs13438494 induces the abrogation or creation of such binding sites. Our results indicate that rs13438494 alters splicing efficiency by creating or disrupting a splicing motif, which functions by binding of splicing regulatory proteins, and may ultimately result in bipolar disorder in affected people.

RNA-Seq approach for genetic improvement of meat quality in pig and evolutionary insight into the substrate specificity of animal carbonyl reductases

Changes in meat quality traits are strongly associated with alterations in postmortem metabolism which depend on genetic variations, especially nonsynonymous single nucleotide variations (nsSNVs) having critical effects on protein structure and function. To selectively identify metabolism-related nsSNVs, next-generation transcriptome sequencing (RNA-Seq) was carried out using RNAs from porcine liver, which contains a diverse range of metabolic enzymes. The multiplex SNV genotyping analysis showed that various metabolism-related genes had different nsSNV alleles. Moreover, many nsSNVs were significantly associated with multiple meat quality traits. Particularly, ch7:g.22112616A>G SNV was identified to create a single amino acid change (Thr/Ala) at the 145th residue of H1.3-like protein, very close to the putative 147th threonine phosphorylation site, suggesting that the nsSNV may affect multiple meat quality traits by affecting the epigenetic regulation of postmortem metabolism-related gene expression. Besides, one nonsynonymous variation, probably generated by gene duplication, led to a stop signal in porcine testicular carbonyl reductase (PTCR), resulting in a C-terminal (E281-A288) deletion. Molecular docking and energy minimization calculations indicated that the binding affinity of wild-type PTCR to 5α-DHT, a C₂₁-steroid, was superior to that of C-terminal-deleted PTCR or human carbonyl reductase, which was very consistent with experimental data, reported previously. Furthermore, P284 was identified as an important residue mediating the specific interaction between PTCR and 5α-DHT, and phylogenetic analysis showed that P284 is an evolutionarily conserved residue among animal carbonyl reductases, which suggests that the C-terminal tails of these reductases may have evolved under evolutionary pressure to increase the substrate specificity for C₂₁-steroids and facilitate metabolic adaptation. Altogether, our RNA-Seq revealed that selective nsSNVs were associated with meat quality traits that could be useful for successful marker-assisted selection in pigs and also represents a useful resource to enhance understanding of protein folding, substrate specificity, and the evolution of enzymes such as carbonyl reductase.

A comprehensive catalogue of functional genetic variations in the EGFR pathway: protein-protein interaction analysis reveals novel genes and polymorphisms important for cancer research

The EGFR pathway is a critical signaling pathway deregulated in many solid tumors. In addition to the initiation and progression of cancer, the EGFR pathway is also implicated in variable treatment responses and prognoses. Genetic variation in the form of Single Nucleotide Polymorphisms (SNPs) can affect the function/expression of the EGFR pathway genes. Here, we applied a systematic and comprehensive approach utilizing diverse public databases and in silico analysis tools to select putative functional genetic variations from 244 genes involved in the EGFR pathway. Our data comprises 649 SNPs. Three hundred sixty SNPs are predicted to have biological consequences (functional SNPs). These SNPs can be directly used in further studies to test their association with risk, treatment response and prognosis in cancer. To systematically cover the EGFR pathway, we also performed a network-based analysis to further select putative functional SNPs from the genes whose protein products physically interact with the EGFR pathway proteins. We utilized protein-protein interaction information and focused on 14 proteins that have a high degree of connectivity (interacting with > or = 10 proteins) with the EGFR pathway genes identified to have functional SNPs (f-EGFR genes). Two of these proteins (FYN and LCK) had interactions with 17 of the f-EGFR genes, yet both lacked any putative functional SNP. However, our analysis indicated the presence of potentially functional SNPs in 9 other highly interactive proteins. The genes and their SNPs identified in the network-based analysis represent potential candidates for gene-gene and SNP-SNP interaction studies in cancer research.

Studying genetic variations in cancer prognosis (and risk): a primer for clinicians

Rare, high-penetrance genetic variations account for a small portion of genetic cancer syndromes. In contrast, most cancers develop from a combination of minor genetic influences and environmental factors. There are numerous publications on cancer susceptibility. In contrast, genetic studies in treatment response and outcome analyses are a rapidly emerging field. Approaches used in disease susceptibility can be adapted for genetic outcome studies. In this review, we summarize the current knowledge on how candidate genes and genetic variations are selected to evaluate gene-outcome, gene-prognosis, and gene-treatment response relationships as applicable to the practicing oncologist.

Human olfaction: from genomic variation to phenotypic diversity

The sense of smell is a complex molecular device, encompassing several hundred olfactory receptor proteins (ORs). These receptors, encoded by the largest human gene superfamily, integrate odorant signals into an accurate 'odor image' in the brain. Widespread phenotypic diversity in human olfaction is, in part, attributable to prevalent genetic variation in OR genes, owing to copy number variation, deletion alleles and deleterious single nucleotide polymorphisms. The development of new genomic tools, including next generation sequencing and CNV assays, provides opportunities to characterize the genetic variations of this system. The advent of large-scale functional screens of expressed ORs, combined with genetic association studies, has the potential to link variations in ORs to human chemosensory phenotypes. This promises to provide a genome-wide view of human olfaction, resulting in a deeper understanding of personalized odor coding, with the potential to decipher flavor and fragrance preferences.

A genome-wide survey of the prevalence and evolutionary forces acting on human nonsense SNPs

Nonsense SNPs introduce premature termination codons into genes and can result in the absence of a gene product or in a truncated and potentially harmful protein, so they are often considered disadvantageous and are associated with disease susceptibility. As such, we might expect the disrupted allele to be rare and, in healthy people, observed only in a heterozygous state. However, some, like those in the CASP12 and ACTN3 genes, are known to be present at high frequencies and to occur often in a homozygous state and seem to have been advantageous in recent human evolution. To evaluate the selective forces acting on nonsense SNPs as a class, we have carried out a large-scale experimental survey of nonsense SNPs in the human genome by genotyping 805 of them (plus control synonymous SNPs) in 1,151 individuals from 56 worldwide populations. We identified 169 genes containing nonsense SNPs that were variable in our samples, of which 99 were found with both copies inactivated in at least one individual. We found that the sampled humans differ on average by 24 genes (out of about 20,000) because of these nonsense SNPs alone. As might be expected, nonsense SNPs as a class were found to be slightly disadvantageous over evolutionary timescales, but a few nevertheless showed signs of being possibly advantageous, as indicated by unusually high levels of population differentiation, long haplotypes, and/or high frequencies of derived alleles. This study underlines the extent of variation in gene content within humans and emphasizes the importance of understanding this type of variation.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Protein interactions in human genetic diseases

A method is presented to identify residues that form part of an interaction interface, leading to the prediction that 1,428 OMIM mutations are related to an interaction defect.

We present a novel method that combines protein structure information with protein interaction data to identify residues that form part of an interaction interface. Our prediction method can retrieve interaction hotspots with an accuracy of 60% (at a 20% false positive rate). The method was applied to all mutations in the Online Mendelian Inheritance in Man (OMIM) database, predicting 1,428 mutations to be related to an interaction defect. Combining predicted and hand-curated sets, we discuss how mutations affect protein interactions in general.

SNP2NMD: A database of human single nucleotide polymorphisms causing nonsense-mediated mRNA decay

Elucidating the effects of genetic polymorphisms on genes and gene networks is an important step in disease association studies. We developed the SNP2NMD database for human SNPs (single nucleotide polymorphisms) that result in PTCs (premature termination codons) and trigger nonsense-mediated mRNA decay (NMD). The SNP2NMD Web interfaces provide extensive genetic information on and graphical views of the queried SNP, gene, and disease terms.

AVAILABILITY

SNP2NMD is available from http://variome.net, or directly from http://bioportal.kobic.re.kr/SNP2NMD.

SUPPLEMENTARY INFORMATION

http://bioportal.kobic.re.kr/SNP2NMD/Wiki.jsp?page=Statistics.