In silico analysis of non-synonymous single nucleotide polymorphisms in human DAZL gene associated with male infertility

High-throughput sequencing and in-silico analysis confirm pathogenicity of novel MSH3 variants in African American colorectal cancer

The maintenance of DNA sequence integrity is critical to avoid accumulation of cancer-causing mutations. Inactivation of DNA Mismatch Repair (MMR) genes (e.g., MLH1 and MSH2) is common among many cancers, including colorectal cancer (CRC) and is the driver of classic microsatellite instability (MSI) in tumors. Somatic MSH3 alterations have been linked to a specific form of MSI called elevated microsatellite alterations at selected tetranucleotide repeats (EMAST) that is associated with patient poor prognosis and elevated among African American (AA) rectal cancer patients. Genetic variants of MSH3 and their pathogenicity vary among different populations, such as among AA, which are not well-represented in publicly available databases. Targeted exome sequencing of MSH3 among AA CRC samples followed by computational bioinformatic pipeline and molecular dynamic simulation analysis approach confirmed six identified MSH3 variants (c.G1237A, c.C2759T, c.G1397A, c.G2926A, c.C3028T, c.G3241A) that corresponded to MSH3 amino-acid changes (p.E413K; p.S466N; p.S920F; p.E976K; p.H1010Y; p.E1081K). All identified MSH3 variants were non-synonymous, novel, pathogenic, and show loss or gain of hydrogen bonding, ionic bonding, hydrophobic bonding, and disulfide bonding and have a deleterious effect on the structure of MSH3 protein. Some variants were located within the ATPase site of MSH3, affecting ATP hydrolysis that is critical for MSH3's function. Other variants were in the MSH3-MSH2 interacting domain, important for MSH3's binding to MSH2. Overall, our data suggest that these variants among AA CRC patients affect the function of MSH3 making them pathogenic and likely contributing to the development or advancement of CRC among AA. Further clarifying functional studies will be necessary to fully understand the impact of these variants on MSH3 function and CRC development in AA patients.

Cloning, expression analysis and localization of DAZL gene implicated in germ cell development of male Hezuo pig

Deleted in azoospermia-like (DAZL) is essential for mammalian testicular function and spermatogenesis. To explore the molecular characterization, expression patterns, and cellular localization of the DAZL in Hezuo pig testes, testicular tissue was isolated from Hezuo pig at five development stages including 30 days old (30 d), 90 days old (90 d), 120 days old (120 d), 180 days old (180 d), and 240 days old (240 d). DAZL cDNA was first cloned using the RT-PCR method, and its molecular characterization was analyzed using relevant bioinformatics software. Subsequently, the expression patterns and cellular localization of DAZL were evaluated using quantitative real-time PCR (qRT-PCR), Western blot, and immunohistochemistry. The cloning and sequence analysis showed that the Hezuo pig DAZL cDNA fragment contained 888 bp open reading frame (ORF) capable of encoding 295 amino acid residues and exhibited high identities with some other mammals. The qRT-PCR and Western blot results indicated that DAZL was specifically expressed in Hezuo pig testes, and DAZL levels of both mRNA and protein were expressed at all five reproductive stages of Hezuo pig testes, with extremely significant higher expression levels in 90 d, 120 d, 180 d, and 240 d than those in 30 d (p < 0.01). Additionally, immunohistochemistry results revealed that DAZL protein was mainly localized in gonocytes at 30 d testes, primary spermatocytes, and spermatozoon at other developmental stages, and Leydig cells throughout five development stages. Together, these results suggested that DAZL may play an important role by regulating the proliferation or differentiation of gonocytes, development of primary spermatocytes and spermatozoon, and functional maintenance of Leydig cells in testicular development and spermatogenesis of Hezuo pig. Nevertheless, the specific regulatory mechanisms underlying these phenomena still requires further investigated and verified.

DAZL Knockout Pigs as Recipients for Spermatogonial Stem Cell Transplantation

Spermatogonial stem cell (SSC) transplantation into the testis of a germ cell (GC)-depleted surrogate allows transmission of donor genotype via donor-derived sperm produced by the recipient. Transplantation of gene-edited SSCs provides an approach to propagate gene-edited large animal models. DAZL is a conserved RNA-binding protein important for GC development, and DAZL knockout (KO) causes defects in GC commitment and differentiation. We characterized DAZL-KO pigs as SSC transplantation recipients. While there were GCs in 1-week-old (wko) KO, complete GC depletion was observed by 10 wko. Donor GCs were transplanted into 18 DAZL-KO recipients at 10-13 wko. At sexual maturity, semen and testes were evaluated for transplantation efficiency and spermatogenesis. Approximately 22% of recipient seminiferous tubules contained GCs, including elongated spermatids and proliferating spermatogonia. The ejaculate of 89% of recipients contained sperm, exclusively from donor origin. However, sperm concentration was lower than the wild-type range. Testicular protein expression and serum hormonal levels were comparable between DAZL-KO and wild-type. Intratesticular testosterone and Leydig cell volume were increased, and Leydig cell number decreased in transplanted DAZL-KO testis compared to wild-type. In summary, DAZL-KO pigs support donor-derived spermatogenesis following SSC transplantation, but low spermatogenic efficiency currently limits their use for the production of offspring.

Al-Safar M, Mahmood SA. Men's ND1 gene genetic makeup Toxoplasmosis and Oligospermia affecting couples' infertility

Introduction: Toxoplasma infection was higher in infertile couples than fertile couples, probably due to anti-sperm antibodies that were higher in couples with Toxoplasmosis. Investigations of T. gondii infections in men with infertility showed that among 100 cases of men's infertility, 36% were serologically positive for Toxoplasma-IgG and IgM. It has been concluded that T. gondii can affect men's fertility and result in infertility. Materials and Methods: Selective infertile males were asked about days of sexual abstinence. Seminal fluid samples were collected following a minimum of 2 days and a maximum of 7 days from abstinence. Every patient was given a clean, wide mouth, sterile, dry, graduated plastic and warm disposable container. The samples were obtained by masturbation in a private room near the semen analysis lab to reduce seminal exposure to temperature fluctuations and control the time from collection to analysis. Results: For the ND1 gene, samples of 8 different fertility groups have been sequenced. These sequences have been compared to reference sequences taken from the NCBI database. Several mutations in various nucleotide positions of the ND1 regions have been detected in samples from multiple groups. The base substitution has been positioned on the nucleotides (nts) 3480, 3567, 3591, 3693, and 4216. The T to C evolution was notorious at nt 3480 in ND1 genes. The SNP was detected in an asthenospermia human (Sample code: 010480). Keywords: Sequence, ND1 gen, Oligospermia Toxoplasmosis, Couples infertility

An In silico Approach towards Finding the Cancer-Causing Mutations in Human MET Gene

Mesenchymal-epithelial transition (MET) factor is a proto-oncogene encoding tyrosine kinase receptor with hepatocyte growth factor (HGF) or scatter factor (SF). It is found on the human chromosome number 7 and regulates the diverse cellular mechanisms of the human body. The impact of mutations occurring in the MET gene is demonstrated by their detrimental effects on normal cellular functions. These mutations can change the structure and function of MET leading to different diseases such as lung cancer, neck cancer, colorectal cancer, and many other complex syndromes. Hence, the current study focused on finding deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) and their subsequent impact on the protein's structure and functions, which may contribute to the emergence of cancers. These nsSNPs were first identified utilizing computational tools like SIFT, PROVEAN, PANTHER-PSEP, PolyPhen-2, I-Mutant 2.0, and MUpro. A total of 45359 SNPs of MET gene were accumulated from the database of dbSNP, and among them, 1306 SNPs were identified as non-synonymous or missense variants. Out of all 1306 nsSNPs, 18 were found to be the most deleterious. Moreover, these nsSNPs exhibited substantial effects on structure, binding affinity with ligand, phylogenetic conservation, secondary structure, and post-translational modification sites of MET, which were evaluated using MutPred2, RaptorX, ConSurf, PSIPRED, and MusiteDeep, respectively. Also, these deleterious nsSNPs were accompanied by changes in properties of MET like residue charge, size, and hydrophobicity. These findings along with the docking results are indicating the potency of the identified SNPs to alter the structure and function of the protein, which may lead to the development of cancers. Nonetheless, Genome-wide association study (GWAS) studies and experimental research are required to confirm the analysis of these nsSNPs.

Seminal plasma extracellular vesicles tRF-Val-AAC-010 can serve as a predictive factor of successful microdissection testicular sperm extraction in patients with non-obstructive azoospermia

BACKGROUND

There is a lack of biomarkers for distinguishing non-obstructive azoospermia (NOA) patients with successful sperm retrieval (Sp+) from those with failed sperm retrieval (Sp-). This study aimed to determine the potential of extracellular vesicles tRNA-derived small RNA (tsRNA) as a novel non-invasive biomarker for successful sperm retrieval by microdissection testicular sperm extraction (mTESE).

METHODS

The study included 18 patients with NOA with successful sperm retrieval (Sp+) and 23 patients with NOA with failed sperm retrieval (Sp-), 15 obstructive azoospermia (OA) patients, 5 idiopathic oligospermia (IO) patients, and 12 healthy people. Seminal plasma extracellular vesicles tsRNA levels were used in a two-stage case-control study (screened by tsRNA sequencing on Illumina NextSeq instrument and validated by qRT-PCR). The bioinformatic analysis was performed to determine the role of tsRNA in the pathogenesis of non-obstructive azoospermia.

RESULTS

Two tsRNAs (tRF-Val-AAC-010: AUC = 0.96, specificity = 80%, sensitivity = 95%; tRF-Pro-AGG-003: AUC = 0.96, specificity = 87%, sensitivity = 95%) were found to have high predictive accuracy for distinguishing the origin of azoospermia. In addition, the extracellular vesicles tRF-Val-AAC-010 resulted in high predictive ability (AUC = 0.89, sensitivity = 72%, specificity = 91%, P < 0.0001) in predicting the presence of sperm in non-obstructive azoospermia undergoing mTESE. Finally, bioinformatic analysis revealed that tRF-Val-AAC-010 were involved in spermatogenesis.

CONCLUSIONS

This study identified that the extracellular vesicles tRF-Val-AAC-010 and tRF-Pro-AGG-003 are biomarkers for the diagnosis of non-obstructive azoospermia, and that tRF-Val-AAC-010 as a potential non-invasive biomarker for predicting the presence of sperm in non-obstructive azoospermia testicular tissue.

Computational Insights into the Structural and Functional Impacts of nsSNPs of Bone Morphogenetic Proteins

BMPs (bone morphogenetic proteins) are multipurpose (transforming growth factor)TGF-superfamily released cytokines. These glycoproteins, acting as disulfide-linked homo- or heterodimers, are highly potent regulators of bone and cartilage production and repair, cell proliferation throughout embryonic development, and bone homeostasis in the adults. Due to the fact that genetic variation might influence structural functions, this study is aimed to determine the pathogenic effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in BMP genes. The implications of these variations, investigated using computational analysis and molecular models of the mature TGF-β domain, revealed the impact of modifications on the function of BMP protein. The three-dimensional (3D) structure analysis was performed on the nsSNP Y316S, V386G, E387G, C389G, and C391G nsSNP in the TGF-β domain of chicken BMP2 and H344P, S347P, V357A nsSNP in the TGF-β domain of chicken BMP4 protein that was anticipated to be harmful and of high risk. The ability of the proteins to perform variety of tasks interact with other molecules depends on their tertiary structural composition. The current analysis revealed the four most damaging variants (Y316S, V386G, E387G, C389G, and C391G), highly conserved and functional and are located in the TGF-beta domain of BMP2 and BMP4. The amino acid substitutions E387G, C389G, and C391G are discovered in the binding region. It was observed that the mutations in the TGF-beta domain caused significant changes in its structural organization including the substrate binding sites. Current findings will assist future research focused on the role of these variants in BMP function loss and their role in skeletal disorders, and this will possibly help to develop practical strategies for treating bone-related conditions.

Autosomal dominant polycystic kidney disease (ADPKD) in Tunisia: From molecular genetics to the development of prognostic tools

A high prevalence of genetic kidney disease in Tunisia has been detected, and their study provides very important clinical and genetic information. Autosomal dominant polycystic kidney disease (ADPKD) is one of the main causes of morbidity and mortality associated with the kidneys in Tunisia. We present here clinical and genetic characteristics of a cohort of Tunisian patients with ADPKD. Nineteen Tunisian patients with ADPKD, among 4 familial cases and 11 sporadic cases, and 50 Healthy individuals were included in this cohort. Genetic studies of PKD1/2 were carried on using Sanger sequencing and MLPA. In our study, the mean age at diagnosis was 47 ± 18 years. In addition, 84.21% of cases present a family history of ADPKD. Overall, 57.89% of the affected individuals had HTA and 26.31% patients had hematuria. 15.78 % of the patient has extra-renal cysts i.e. one patient with splenic cysts and two patients had liver cysts. 57.89 % of patients were diagnosed with various extra-renal clinical presentations i.e. myopia, hernia, deafness, intracranial aneurysm, respiratory distress, hyperthyroidism, urinary tract infection and lower back pains. The PKD1 genotype showed earlier onset of ESRD compared to PKD2 genotype (43 vs. 55 years old). Six mutations have been detected in PKD1 gene. Among them, three were novels e.g. c.688 T>G, p.C230G and c.690C>G, p.C230W among exon 5 and c.8522A>G, p.N2841S among exon 23. In addition, thirteen single nucleotides polymorphisms have been reported in PKD1 gene. Among them, eleven previously reported in heterozygous state and two novel single nucleotides polymorphisms in heterozygous and homozygous state and predicted to be probable polymorphisms by computational tools: c.496C>T, p.L166= among the exon 4, and c.10165G>C and p.E3389Gln among the exon 31. Only three single nucleotides polymorphisms previously reported in ADPKD database have been identified in PKD2 gene. The description and analysis of our cohort can help in rapid and reliable diagnosis for early management of patients in Tunisia. Indeed, predictive genetic testing can facilitate donor evaluation and increase living related kidney transplantation.

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

Functional Analysis of Single Nucleotide Polymorphism in ZUFSP Protein and Implication in Pathogenesis

Researches have revealed that functional non-synonymous Single Nucleotide Polymorphism (nsSNPs) present in the Zinc-finger with UFM1-Specific Peptidase domain protein (ZUFSP) may be involved in genetic instability and carcinogenesis. For the first time, we employed in-silico approach using predictive tools to identify and validate potential nsSNPs that could be pathogenic. Our result revealed that 8 nsSNPs (rs 112738382, rs 140094037, rs 201652589, rs 201847265, rs 202076827, rs 373634906, rs 375114528, rs 772591104) are pathogenic after being subjected to rigorous filtering process. The structural impact of the nsSNPs on ZUFSP structure indicated that the nsSNPs affect the stability of the protein by lowering ZUFSP protein stability. Furthermore, conservation analysis showed that rs 201652589, rs 140094037, rs 201847265, and rs 772591104 were highly conserved. Interestingly, the protein-protein affinity between ZUFSP and Ubiquitin was altered rs 201652589, rs 140094037, rs 201847265, and rs 772591104 had a binding affinity of - 0.46, - 0.83, - 1.62, and - 1.12 kcal/mol respectively. Our study has been able to identify potential nsSNPs that could be used as genetic biomarkers for some diseases arising as a result of aberration in the ZUFSP structure, however, being a predictive study, the identified nsSNPs need to be experimentally investigated.

In silico analysis of a novel causative mutation in Cadherin23 gene identified in an Omani family with hearing loss

Background

Hereditary hearing loss is a heterogeneous group of complex disorders with an overall incidence of one in every 500 newborns presented as syndromic and non-syndromic forms. Cadherin-related 23 (CDH23) is one of the listed deafness causative genes. It is found to be expressed in the stereocilia of hair cells and in the retina photoreceptor cells. Defective CDH23 have been associated mostly with prelingual severe-to-profound sensorineural hearing loss (SNHL) in either syndromic (USH1D) or non-syndromic SNHL (DFNB12) deafness. The purpose of this study was to identify causative mutations in an Omani family diagnosed with severe-profound sensorineural hearing loss by whole exome sequencing technique and analyzing the detected variant in silico for pathogenicity using several in silico mutation prediction software.

Results

A novel homozygous missense variant, c.A7436C (p. D2479A), in exon 53 of CDH23 was detected in the family while the control samples were all negative for the detected variant. In silico mutation prediction analysis showed the novel substituted D2479A to be deleterious and protein destabilizing mutation at a conserved site on CDH23 protein.

Conclusion

In silico mutation prediction analysis might be used as a useful molecular diagnostic tool benefiting both genetic counseling and mutation verification. The aspartic acid 2479 alanine missense substitution might be the main disease-causing mutation that damages CDH23 function and could be used as a genetic hearing loss marker for this particular Omani family.

Gene expression patterns and protein cellular localization suggest a novel role for DAZL in developing Tibetan sheep testes

Deleted in azoospermia-like (DAZL) is essential for mammalian spermatogenesis as it regulates proliferation, development, maturation and functional maintenance of male germ cells. Its expression and regulation vary with different species or in the same animal at different developmental stages, and despite its importance, very little is known about its roles in sheep, especially Tibetan sheep. To investigate the expression patterns and regulatory roles of DZAL in Tibetan sheep testis, testicular tissue was isolated from sheep at three crucial development stages: 3 months old, 1 year old and 3 years old. Using quantitative real-time PCR and Western blot, we found that DAZL mRNA first decreased and then increased with advancing age, while DAZL protein exhibited an opposite expression pattern, with first increased and subsequently decreased levels. Immunohistochemistry and immunofluorescence revealed that DAZL protein was located predominantly in the cytoplasm of Leydig cells and in both the cytoplasm and nucleus of spermatids. ELISA indicated that testosterone content within developing testes was first enhanced and then declined. Our results, taken together, demonstrate, for the first time, that DAZL gene is involved in Tibetan sheep spermatogenesis by regulating the development of spermatids in post-pubertal rams, along with a novel role in functional maintenance of Leydig cells in postnatal rams.

Computational Analysis of High-Risk SNPs in Human DBY Gene Responsible for Male Infertility: A Functional and Structural Impact

BACKGROUND

DEAD-box helicase 3, Y-linked (DBY) is a candidate gene of the AZF region which is involved in spermatogenesis process. Mutations in the DBY gene may disrupt the spermatogenesis and lead to infertility in men. Identification of functionally neutral mutation from the disease-causing mutation is the biggest challenge in human genetic variation analysis. Owing to the importance of DBY in male infertility, functional analysis was carried out to reveal the association between genetic mutation and phenotypic variation through various in silico approaches.

METHODS

The present study analyzed the functional consequences of the nsSNPs in human DBY gene using SIFT, PolyPhen 2, PROVEAN, SNAP2, PMut, nsSNPAnalyzer, PhD-SNP and SNPs&GO along with stability analysis through I-Mutant2.0, MuPro and iPTREE-STAB. The conservational analysis of amino acid residues, biophysical properties and conserved domains of the DBY protein was analyzed using various computational tools. The 3D structure of the protein was generated using SPARKS-X and validated using RAMPAGE.

RESULTS

Out of 1130 SNPs reported in dbSNP, only one nsSNP (G300D) was found to have a functional effect on stability as well as the function of the DBY protein. The results showed the presence of G300 in the putative structure of DBY domain.

CONCLUSION

To the best of our knowledge, this is the first study to detect pathologically significant nsSNPs (G300D) through a computational approach in the DBY which can be useful for development in potent drug discovery studies.

Azoospermia Factor C Subregion of the Y Chromosome

The azoospermia factor (AZF) region on the Y chromosome consists of genes required for spermatogenesis. Among the three subregions, the AZFc subregion located at the distal portion of AZF is the driver for genetic variation in Y chromosome. The candidate gene of AZFc is known as deleted in azoospermia gene, which is studied with interest because it is involved in germ cell development and most frequently deleted genes leading to oligozoospermia and azoospermia. Recently, two partial deletions in AZFc gr/gr and b2/b3 are characterized at the molecular level which showed homologous recombination between amplicons, affecting spermatogenesis process. There are novel methods and commercially available kits for accurate screening and characterization of microdeletions. It is important to detect the AZFc microdeletions through genetic screening and counseling those infertile men who planned to avail assisted reproduction techniques such as undergoing intracytoplasmic sperm injection or in vitro fertilization.